Balancing management and preservation of mediterranean scattered oak woodlands (Dehesas) in human-dominated landscapes

Los sistemas agroforestales mediterráneos tienen una gran importancia ecológica y socioeconómica, y mantienen altos valores medioambientales y de diversidad biológica a la vez que producen importantes servicios ecosistémicos. Estos sistemas han sido testigos de diversos cambios rápidos y drásticos en su gestión y aprovechamiento durante el último siglo. La mayor parte de la investigación desarrollada en esta tesis doctoral ha sido llevada a cabo en las dehesas españolas. Esta tesis nos muestra: i) la evidencia de la existencia de un cambio global del estrato arbóreo y del manejo del pastoreo en el todo el área de distribución de la dehesa durante los últimos 60 años; ii) la importancia del papel que juega el arbolado disperso y el adecuado manejo del ganado en la mejora de la producción, calidad y diversidad de las comunidades herbáceas, que a su vez, un pasto herbáceo bien desarrollado es importante para la rentabilidad del sistema, evaluando estos efectos bajo distintos escenarios de clima y calidad de estación; y iii) la evidencia de la falta de regeneración en sistemas agroforestales mediterráneos bajo distintos tipos de manejo del pastoreo, y además se evalúa el crecimiento y desarrollo de las pocas plántulas existentes que serán las que aseguren la viabilidad y persistencia y de estos sistemas. El arbolado disperso de estos sistemas ha experimentado una reducción importante en su densidad arbórea y fracción de cabida cubierta durante el periodo entre 1950-1980 donde tuvieron lugar importantes transformaciones en la actividad agropecuaria. La cabaña ganadera de ovino disminuyó drásticamente en los años 70 en comparación a la de bovino que desde entonces ha aumentado progresivamente hasta la actualidad. Por otro lado, el mismo tipo de manejo del ganado doméstico (especialmente bovino) durante bastante tiempo (mínimo 30 años) provocó una reducción significativa de la densidad de las plántulas. Además la probabilidad de ocurrencia y la intensidad de daños por herbivoría fue mayor bajo pastoreo bovino (con daños más intensos y consistentes) que bajo pastoreo ovino o sin pastoreo doméstico (presencia de ciervos). También el patrón de crecimiento de las plantas jóvenes estuvo afectado por el tipo de manejo, generando plántulas achaparradas en el caso del bovino y plántulas esbeltas favoreciendo el crecimiento en altura en el caso del ovino. La presencia de un arbolado disperso generó una mayor diversidad y variación en la producción de las comunidades herbáceas según las condiciones de disponibilidad de agua. Especialmente, el ecotono como microhábitat sostuvo altos valores de diversidad herbácea. La presencia del ganado bajo pastoreo continuo de intensidad moderada a alta, especialmente el bovino, incrementó los rendimientos de producción y diversidad del estrato herbáceo. Los resultados de esta tesis nos muestran la importancia que tiene la existencia de un equilibrio entre la producción y la conservación de los sistemas agroforestales mediterráneos para obtener una producción sostenible de servicios ecosistémicos mientras se asegura la perpetuación del sistema a largo plazo. Es crucial diseñar planes de gestión incorporando objetivos de conservación que integren técnicas silvopastorales apropiadas para poder aplicar en los sistemas agroforestales mediterráneos. ABSTRACT Mediterranean scattered oak woodlands have great ecological and socio-economic importance, supporting high environmental and amenity values, and relatively rich biological diversity while producing important ecosystem services. They have been witnesses of different and fast changes developed in the last century. Most of the research developed in this dissertation has conducted within dehesas. This thesis provides: i) the global change evidence of the tree layer and grazing management experienced in the land-use range of a Mediterranean scattered oak woodland (dehesa) over the last 60 years; ii) the important role of scattered trees and adequate management grazing in the improvement of grassland yield, quality and diversity - which it is important, in turn, for the system profitability - under different climate scenarios and site quality; and iii) the lack of oak regeneration evidence under some given representative management regimes and how is the growth development of these plants to assure the viability and persistence of Mediterranean scattered oak woodlands. Tree layer experienced a significant reduction in dehesas during 1950-1980 period where the highest human impacts took place. Sheep herd decreased drastically during the 1970s and, in contrast, cattle have been increasing gradually since then. On the other hand, same livestock grazing management (especially cattle) during long time (minimum 30 years) within Mediterranean scattered oak woodlands reduced strongly the density of young oak plants and showed high probability of herbivory occurrence and intensity. Young plant growth pattern was greatly modified by livestock. Cattle grazing generated stunted plants and sheep grazing generated slender plants favoring the height growth. Microsites created by large trees modified the herbaceous yield according the water availability of the year and generated high plant diversity within herbaceous communities. Especially, ecotone microsite supported high values of herbaceous diversity. The presence of livestock species increased the herbaceous yield and maintained a more diverse community under continuous grazing at both moderate and high intensities; especially cattle. Thus, around the influence of scattered trees there is a high amount of different interactions among livestock, trees and grasslands maintaining and enhancing the quality of whole dehesa system. The results of this thesis highlight how important is balancing management and preservation of Mediterranean scattered oak woodlands to obtain the optimum ecosystem services while the system conservation is assured for a long-term. It is crucial to design management plans with conservation goals that include appropriate silvopastoral practices in Mediterranean scattered oak woodlands.

Balancing management and preservation of mediterranean scattered oak woodlands (Dehesas) in human-dominated landscapes

Los sistemas agroforestales mediterráneos tienen una gran importancia ecológica y socioeconómica, y mantienen altos valores medioambientales y de diversidad biológica a la vez que producen importantes servicios ecosistémicos. Estos sistemas han sido testigos de diversos cambios rápidos y drásticos en su gestión y aprovechamiento durante el último siglo. La mayor parte de la investigación desarrollada en esta tesis doctoral ha sido llevada a cabo en las dehesas españolas. Esta tesis nos muestra: i) la evidencia de la existencia de un cambio global del estrato arbóreo y del manejo del pastoreo en el todo el área de distribución de la dehesa durante los últimos 60 años; ii) la importancia del papel que juega el arbolado disperso y el adecuado manejo del ganado en la mejora de la producción, calidad y diversidad de las comunidades herbáceas, que a su vez, un pasto herbáceo bien desarrollado es importante para la rentabilidad del sistema, evaluando estos efectos bajo distintos escenarios de clima y calidad de estación; y iii) la evidencia de la falta de regeneración en sistemas agroforestales mediterráneos bajo distintos tipos de manejo del pastoreo, y además se evalúa el crecimiento y desarrollo de las pocas plántulas existentes que serán las que aseguren la viabilidad y persistencia y de estos sistemas. El arbolado disperso de estos sistemas ha experimentado una reducción importante en su densidad arbórea y fracción de cabida cubierta durante el periodo entre 1950-1980 donde tuvieron lugar importantes transformaciones en la actividad agropecuaria. La cabaña ganadera de ovino disminuyó drásticamente en los años 70 en comparación a la de bovino que desde entonces ha aumentado progresivamente hasta la actualidad. Por otro lado, el mismo tipo de manejo del ganado doméstico (especialmente bovino) durante bastante tiempo (mínimo 30 años) provocó una reducción significativa de la densidad de las plántulas. Además la probabilidad de ocurrencia y la intensidad de daños por herbivoría fue mayor bajo pastoreo bovino (con daños más intensos y consistentes) que bajo pastoreo ovino o sin pastoreo doméstico (presencia de ciervos). También el patrón de crecimiento de las plantas jóvenes estuvo afectado por el tipo de manejo, generando plántulas achaparradas en el caso del bovino y plántulas esbeltas favoreciendo el crecimiento en altura en el caso del ovino. La presencia de un arbolado disperso generó una mayor diversidad y variación en la producción de las comunidades herbáceas según las condiciones de disponibilidad de agua. Especialmente, el ecotono como microhábitat sostuvo altos valores de diversidad herbácea. La presencia del ganado bajo pastoreo continuo de intensidad moderada a alta, especialmente el bovino, incrementó los rendimientos de producción y diversidad del estrato herbáceo. Los resultados de esta tesis nos muestran la importancia que tiene la existencia de un equilibrio entre la producción y la conservación de los sistemas agroforestales mediterráneos para obtener una producción sostenible de servicios ecosistémicos mientras se asegura la perpetuación del sistema a largo plazo. Es crucial diseñar planes de gestión incorporando objetivos de conservación que integren técnicas silvopastorales apropiadas para poder aplicar en los sistemas agroforestales mediterráneos. ABSTRACT Mediterranean scattered oak woodlands have great ecological and socio-economic importance, supporting high environmental and amenity values, and relatively rich biological diversity while producing important ecosystem services. They have been witnesses of different and fast changes developed in the last century. Most of the research developed in this dissertation has conducted within dehesas. This thesis provides: i) the global change evidence of the tree layer and grazing management experienced in the land-use range of a Mediterranean scattered oak woodland (dehesa) over the last 60 years; ii) the important role of scattered trees and adequate management grazing in the improvement of grassland yield, quality and diversity - which it is important, in turn, for the system profitability - under different climate scenarios and site quality; and iii) the lack of oak regeneration evidence under some given representative management regimes and how is the growth development of these plants to assure the viability and persistence of Mediterranean scattered oak woodlands. Tree layer experienced a significant reduction in dehesas during 1950-1980 period where the highest human impacts took place. Sheep herd decreased drastically during the 1970s and, in contrast, cattle have been increasing gradually since then. On the other hand, same livestock grazing management (especially cattle) during long time (minimum 30 years) within Mediterranean scattered oak woodlands reduced strongly the density of young oak plants and showed high probability of herbivory occurrence and intensity. Young plant growth pattern was greatly modified by livestock. Cattle grazing generated stunted plants and sheep grazing generated slender plants favoring the height growth. Microsites created by large trees modified the herbaceous yield according the water availability of the year and generated high plant diversity within herbaceous communities. Especially, ecotone microsite supported high values of herbaceous diversity. The presence of livestock species increased the herbaceous yield and maintained a more diverse community under continuous grazing at both moderate and high intensities; especially cattle. Thus, around the influence of scattered trees there is a high amount of different interactions among livestock, trees and grasslands maintaining and enhancing the quality of whole dehesa system. The results of this thesis highlight how important is balancing management and preservation of Mediterranean scattered oak woodlands to obtain the optimum ecosystem services while the system conservation is assured for a long-term. It is crucial to design management plans with conservation goals that include appropriate silvopastoral practices in Mediterranean scattered oak woodlands.

Modelo para la integración del trazado de infraestructuras lineales en el paisaje y el medio ambiente basado en SIG

El análisis de las diferentes alternativas en la planificación y diseño de corredores y trazados de carreteras debe basarse en la correcta definición de variables territoriales que sirvan como criterios para la toma de decisión y esto requiere un análisis ambiental preliminar de esas variables de calidad. En España, los estudios de viabilidad de nuevas carreteras y autovías están asociados a una fase del proceso de decisión que se corresponde con el denominado Estudio Informativo, el cual establece condicionantes físicos, ambientales, de uso del suelo y culturales que deben ser considerados en las primeras fases de la definición del trazado de un corredor de carretera. Así, la metodología más frecuente es establecer diferentes niveles de capacidad de acogida del territorio en el área de estudio con el fin de resumir las variables territoriales en mapas temáticos y facilitar el proceso de trazado de las alternativas de corredores de carretera. El paisaje es un factor limitante a tener en cuenta en la planificación y diseño de carreteras y, por tanto, deben buscarse trazados más sostenibles en relación con criterios estéticos y ecológicos del mismo. Pero este factor no es frecuentemente analizado en los Estudios Informativos e incluso, si es considerado, los estudios específicos de la calidad del paisaje (estético y ecológico) y de las formas del terreno no incorporan las recomendaciones de las guías de trazado para evitar o reducir los impactos en el paisaje. Además, los mapas de paisaje que se generan en este tipo de estudios no se corresponden con la escala de desarrollo del Estudio Informativo (1:5.000). Otro déficit común en planificación de corredores y trazados de carreteras es que no se tiene en cuenta la conectividad del paisaje durante el proceso de diseño de la carretera para prevenir la afección a los corredores de fauna existentes en el paisaje. Este déficit puede originar un posterior efecto barrera en los movimientos dispersivos de la fauna y la fragmentación de sus hábitats debido a la ocupación parcial o total de las teselas de hábitats con importancia biológica para la fauna (o hábitats focales) y a la interrupción de los corredores de fauna que concentran esos movimientos dispersivos de la fauna entre teselas. El objetivo principal de esta tesis es mejorar el estudio del paisaje para prevenir su afección durante el proceso de trazado de carreteras, facilitar la conservación de los corredores de fauna (o pasillos verdes) y la localización de medidas preventivas y correctoras en términos de selección y cuantificación de factores de idoneidad a fin de reducir los impactos visuales y ecológicos en el paisaje a escala local. Concretamente, la incorporación de valores cuantitativos y bien justificados en el proceso de decisión permite incrementar la transparencia en el proceso de diseño de corredores y trazados de carreteras. Con este fin, se han planteado cuatro preguntas específicas en esta investigación (1) ¿Cómo se seleccionan y evalúan los factores territoriales limitantes para localizar una nueva carretera por los profesionales españoles de planificación del territorio en relación con el paisaje? (2) ¿Cómo pueden ser definidos los corredores de fauna a partir de factores del paisaje que influyen en los movimientos dispersivos de la fauna? (3) ¿Cómo pueden delimitarse y evaluarse los corredores de fauna incluyendo el comportamiento parcialmente errático en los movimientos dispersivos de la fauna y el efecto barrera de los elementos antrópicos a una escala local? (4) ¿Qué y cómo las recomendaciones de diseño de carreteras relacionadas con el paisaje y las formas del terreno pueden ser incluidas en un modelo de Sistemas de Información Geográfica (SIG) para ayudar a los ingenieros civiles durante el proceso de diseño de un trazado de carreteras bajo el punto de vista de la sostenibilidad?. Esta tesis doctoral propone nuevas metodologías que mejoran el análisis visual y ecológico del paisaje utilizando indicadores y modelos SIG para obtener alternativas de trazado que produzcan un menor impacto en el paisaje. Estas metodologías fueron probadas en un paisaje heterogéneo con una alta tasa de densidad de corzo (Capreolus capreolus L.), uno de los grandes mamíferos más atropellados en la red de carreteras españolas, y donde está planificada la construcción de una nueva autovía que atravesará la mitad del área de distribución del corzo. Inicialmente, se han analizado las variables utilizadas en 22 estudios de proyectos de planificación de corredores de carreteras promovidos por el Ministerio de Fomento entre 2006 y 2008. Estas variables se agruparon según condicionantes físicos, ambientales, de usos del suelo y culturales con el fin de comparar los valores asignados de capacidad de acogida del territorio a cada variable en los diferentes estudios revisados. Posteriormente, y como etapa previa de un análisis de conectividad, se construyó un mapa de resistencia de los movimientos dispersivos del corzo en base a la literatura y al juicio de expertos. Usando esta investigación como base, se le asignó un valor de resistencia a cada factor seleccionado para construir la matriz de resistencia, ponderándolo y combinándolo con el resto de factores usando el proceso analítico jerárquico y los operadores de lógica difusa como métodos de análisis multicriterio. Posteriormente, se diseñó una metodología SIG para delimitar claramente la extensión física de los corredores de fauna de acuerdo a un valor umbral de ancho geométrico mínimo, así como la existencia de múltiples potenciales conexiones entre cada par de teselas de hábitats presentes en el paisaje estudiado. Finalmente, se realizó un procesado de datos Light Detection and Ranging (LiDAR) y un modelo SIG para calcular la calidad del paisaje (estético y ecológico), las formas del terreno que presentan características similares para trazar una carretera y la acumulación de vistas de potenciales conductores y observadores de los alrededores de la nueva vía. Las principales contribuciones de esta investigación al conocimiento científico existente en el campo de la evaluación del impacto ambiental en relación al diseño de corredores y trazados de carreteras son cuatro. Primero, el análisis realizado de 22 Estudios Informativos de planificación de carreteras reveló que los métodos aplicados por los profesionales para la evaluación de la capacidad de acogida del territorio no fue suficientemente estandarizada, ya que había una falta de uniformidad en el uso de fuentes cartográficas y en las metodologías de evaluación de la capacidad de acogida del territorio, especialmente en el análisis de la calidad del paisaje estético y ecológico. Segundo, el análisis realizado en esta tesis destaca la importancia de los métodos multicriterio para estructurar, combinar y validar factores que limitan los movimientos dispersivos de la fauna en el análisis de conectividad. Tercero, los modelos SIG desarrollados Generador de alternativas de corredores o Generator of Alternative Corridors (GAC) y Eliminador de Corredores Estrechos o Narrow Corridor Eraser (NCE) pueden ser aplicados sistemáticamente y sobre una base científica en análisis de conectividad como una mejora de las herramientas existentes para la comprensión el paisaje como una red compuesta por nodos y enlaces interconectados. Así, ejecutando los modelos GAC y NCE de forma iterativa, pueden obtenerse corredores alternativos con similar probabilidad de ser utilizados por la fauna y sin que éstos presenten cuellos de botella. Cuarto, el caso de estudio llevado a cabo de prediseño de corredores y trazado de una nueva autovía ha sido novedoso incluyendo una clasificación semisupervisada de las formas del terreno, filtrando una nube de puntos LiDAR e incluyendo la nueva geometría 3D de la carretera en el Modelo Digital de Superficie (MDS). El uso combinado del procesamiento de datos LiDAR y de índices y clasificaciones geomorfológicas puede ayudar a los responsables encargados en la toma de decisiones a evaluar qué alternativas de trazado causan el menor impacto en el paisaje, proporciona una visión global de los juicios de valor más aplicados y, en conclusión, define qué medidas de integración paisajística correctoras deben aplicarse y dónde. ABSTRACT The assessment of different alternatives in road-corridor planning and layout design must be based on a number of well-defined territorial variables that serve as decision-making criteria, and this requires a high-quality preliminary environmental analysis of those quality variables. In Spain, feasibility studies for new roads and motorways are associated to a phase of the decision procedure which corresponds with the one known as the Informative Study, which establishes the physical, environmental, land-use and cultural constraints to be considered in the early stages of defining road corridor layouts. The most common methodology is to establish different levels of Territorial Carrying Capacity (TCC) in the study area in order to summarize the territorial variables on thematic maps and facilitate the tracing process of road-corridor layout alternatives. Landscape is a constraint factor that must be considered in road planning and design, and the most sustainable layouts should be sought based on aesthetic and ecological criteria. However this factor is not often analyzed in Informative Studies and even if it is, baseline studies on landscape quality (aesthetic and ecological) and landforms do not usually include the recommendations of road tracing guides designed to avoid or reduce impacts on the landscape. The resolution of the landscape maps produced in this type of studies does not comply with the recommended road design scale (1:5,000) in the regulations for the Informative Study procedure. Another common shortcoming in road planning is that landscape ecological connectivity is not considered during road design in order to avoid affecting wildlife corridors in the landscape. In the prior road planning stage, this issue could lead to a major barrier effect for fauna dispersal movements and to the fragmentation of their habitat due to the partial or total occupation of habitat patches of biological importance for the fauna (or focal habitats), and the interruption of wildlife corridors that concentrate fauna dispersal movements between patches. The main goal of this dissertation is to improve the study of the landscape and prevent negative effects during the road tracing process, and facilitate the preservation of wildlife corridors (or green ways) and the location of preventive and corrective measures by selecting and quantifying suitability factors to reduce visual and ecological landscape impacts at a local scale. Specifically the incorporation of quantitative and well-supported values in the decision-making process provides increased transparency in the road corridors and layouts design process. Four specific questions were raised in this research: (1) How are territorial constraints selected and evaluated in terms of landscape by Spanish land-planning practitioners before locating a new road? (2) How can wildlife corridors be defined based on the landscape factors influencing the dispersal movements of fauna? (3) How can wildlife corridors be delimited and assessed to include the partially erratic movements of fauna and the barrier effect of the anthropic elements at a local scale? (4) How recommendations of road design related to landscape and landforms can be included in a Geographic Information System (GIS) model to aid civil engineers during the road layout design process and support sustainable development? This doctoral thesis proposes new methodologies that improve the assessment of the visual and ecological landscape character using indicators and GIS models to obtain road layout alternatives with a lower impact on the landscape. These methodologies were tested on a case study of a heterogeneous landscape with a high density of roe deer (Capreolus capreolus L.) –one of the large mammals most commonly hit by vehicles on the Spanish road network– and where a new motorway is planned to pass through the middle of their distribution area. We explored the variables used in 22 road-corridor planning projects sponsored by the Ministry of Public Works between 2006 and 2008. These variables were grouped into physical, environmental, land-use and cultural constraints for the purpose of comparing the TCC values assigned to each variable in the various studies reviewed. As a prior stage in a connectivity analysis, a map of resistance to roe deer dispersal movements was created based on the literature and experts judgment. Using this research as a base, each factor selected to build the matrix was assigned a resistance value and weighted and combined with the rest of the factors using the analytic hierarchy process (AHP) and fuzzy logic operators as multicriteria assessment (MCA) methods. A GIS methodology was designed to clearly delimit the physical area of wildlife corridors according to a geometric threshold width value, and the multiple potential connections between each pair of habitat patches in the landscape. A Digital Surface Model Light Detection and Ranging (LiDAR) dataset processing and a GIS model was performed to determine landscape quality (aesthetic and ecological) and landforms with similar characteristics for the road layout, and the cumulative viewshed of potential drivers and observers in the area surrounding the new motorway. The main contributions of this research to current scientific knowledge in the field of environmental impact assessment for road corridors and layouts design are four. First, the analysis of 22 Informative Studies on road planning revealed that the methods applied by practitioners for assessing the TCC were not sufficiently standardized due to the lack of uniformity in the cartographic information sources and the TCC valuation methodologies, especially in the analysis of the aesthetic and ecological quality of the landscape. Second, the analysis in this dissertation highlights the importance of multicriteria methods to structure, combine and validate factors that constrain wildlife dispersal movements in the connectivity analysis. Third, the “Generator of Alternative Corridors (GAC)” and “Narrow Corridor Eraser (NCE)” GIS models developed can be applied systematically and on a scientific basis in connectivity analyses to improve existing tools and understand landscape as a network composed of interconnected nodes and links. Thus, alternative corridors with similar probability of use by fauna and without bottlenecks can be obtained by iteratively running GAC and NCE models. Fourth, our case study of new motorway corridors and layouts design innovatively included semi-supervised classification of landforms, filtering of LiDAR point clouds and new 3D road geometry on the Digital Surface Model (DSM). The combined used of LiDAR data processing and geomorphological indices and classifications can help decision-makers assess which road layouts produce lower impacts on the landscape, provide an overall insight into the most commonly applied value judgments, and in conclusion, define which corrective measures should be applied in terms of landscaping, and where.

Caracterización y simulación de arborizaciones dentríticas con redes bayesianas incluyendo variables angulares

El funcionamiento interno del cerebro es todavía hoy en día un misterio, siendo su comprensión uno de los principales desafíos a los que se enfrenta la ciencia moderna. El córtex cerebral es el área del cerebro donde tienen lugar los procesos cerebrales de más alto nivel, cómo la imaginación, el juicio o el pensamiento abstracto. Las neuronas piramidales, un tipo específico de neurona, suponen cerca del 80% de los cerca de los 10.000 millones de que componen el córtex cerebral, haciendo de ellas un objetivo principal en el estudio del funcionamiento del cerebro. La morfología neuronal, y más específicamente la morfología dendrítica, determina cómo estas procesan la información y los patrones de conexión entre neuronas, siendo los modelos computacionales herramientas imprescindibles para el estudio de su rol en el funcionamiento del cerebro. En este trabajo hemos creado un modelo computacional, con más de 50 variables relativas a la morfología dendrítica, capaz de simular el crecimiento de arborizaciones dendríticas basales completas a partir de reconstrucciones de neuronas piramidales reales, abarcando desde el número de dendritas hasta el crecimiento los los árboles dendríticos. A diferencia de los trabajos anteriores, nuestro modelo basado en redes Bayesianas contempla la arborización dendrítica en su conjunto, teniendo en cuenta las interacciones entre dendritas y detectando de forma automática las relaciones entre las variables morfológicas que caracterizan la arborización. Además, el análisis de las redes Bayesianas puede ayudar a identificar relaciones hasta ahora desconocidas entre variables morfológicas. Motivado por el estudio de la orientación de las dendritas basales, en este trabajo se introduce una regularización L1 generalizada, aplicada al aprendizaje de la distribución von Mises multivariante, una de las principales distribuciones de probabilidad direccional multivariante. También se propone una distancia circular multivariante que puede utilizarse para estimar la divergencia de Kullback-Leibler entre dos muestras de datos circulares. Comparamos los modelos con y sin regularizaci ón en el estudio de la orientación de la dendritas basales en neuronas humanas, comprobando que, en general, el modelo regularizado obtiene mejores resultados. El muestreo, ajuste y representación de la distribución von Mises multivariante se implementa en un nuevo paquete de R denominado mvCircular.---ABSTRACT---The inner workings of the brain are, as of today, a mystery. To understand the brain is one of the main challenges faced by current science. The cerebral cortex is the region of the brain where all superior brain processes, like imagination, judge and abstract reasoning take place. Pyramidal neurons, a specific type of neurons, constitute approximately the 80% of the more than 10.000 million neurons that compound the cerebral cortex. It makes the study of the pyramidal neurons crucial in order to understand how the brain works. Neuron morphology, and specifically the dendritic morphology, determines how the information is processed in the neurons, as well as the connection patterns among neurons. Computational models are one of the main tools for studying dendritic morphology and its role in the brain function. We have built a computational model that contains more than 50 morphological variables of the dendritic arborizations. This model is able to simulate the growth of complete dendritic arborizations from real neuron reconstructions, starting with the number of basal dendrites, and ending modeling the growth of dendritic trees. One of the main diferences between our approach, mainly based on the use of Bayesian networks, and other models in the state of the art is that we model the whole dendritic arborization instead of focusing on individual trees, which makes us able to take into account the interactions between dendrites and to automatically detect relationships between the morphologic variables that characterize the arborization. Moreover, the posterior analysis of the relationships in the model can help to identify new relations between morphological variables. Motivated by the study of the basal dendrites orientation, a generalized L1 regularization applied to the multivariate von Mises distribution, one of the most used distributions in multivariate directional statistics, is also introduced in this work. We also propose a circular multivariate distance that can be used to estimate the Kullback-Leibler divergence between two circular data samples. We compare the regularized and unregularized models on basal dendrites orientation of human neurons and prove that regularized model achieves better results than non regularized von Mises model. Sampling, fitting and plotting functions for the multivariate von Mises are implemented in a new R packaged called mvCircular.

Riesgos de taludes y terraplenes en obras lineales de la República Dominicana

Partiendo de la idea de la concurrencia de un fenómeno y de las acciones que se deben llevar a cabo, se debe entender la Gestión de Riesgos como una determinación de relaciones entre lo que se asume como vulnerable, y la forma en la cual se estimará la probabilidad de un determinado hecho. A medida en que pasa el tiempo, se torna más importante el tema de vulnerabilidad y de riesgo a nivel mundial. La vulnerabilidad se vuelve más notoria ante la presencia de determinados peligros naturales como lo son los siguientes: deslizamientos de tierras, movimientos sísmicos, desbordes de ríos e inundaciones. (Pérez Soriano, 2014) La vulnerabilidad se ve en su mayoría más afectada por acciones antrópicas como lo son la construcción de viviendas en lugares de alto riesgo, como por ejemplo, viviendas en los taludes de obras lineales. Estas viviendas están condicionadas por la localización, el uso del suelo, la infraestructura, la distribución y la densidad de la población, su capacidad de organización, etc. La Gestión de Riesgos se vuelve cada vez más exigente en cuanto a la calidad de los servicios ofertados, además de también, cumplir con la responsabilidad jurídica que implica la concepción, diseño y construcción de taludes y terraplenes en obras lineales. (Pérez Soriano, 2014) El presente trabajo de investigación se centra en la identificación y valoración de los riesgos de taludes y terraplenes en obras lineales de la República Dominicana. De esta manera se busca disminuir el riesgo de la rotura o fallo de un talud o terraplén, y el número de víctimas que puedan ser afectadas por el mismo, concluyendo con un catálogo de riesgos para taludes y terraplenes en obras lineales. Based on the idea of the occurrence of a phenomenon and the actions to be carried out, Risk Management should be understood as a determination of relationships between what is assumed to be vulnerable, and the way in which the probability of a particular event will be estimated. The issue of vulnerability and risk becomes more important worldwide as time goes on. Vulnerability becomes more evident in the presence of certain natural hazards such as: landslides, earthquakes, overflowing rivers and flooding. (Pérez Soriano, 2014) It has become evident that vulnerability is mostly affected by human actions, such as the construction of housing in high-risk locations, for example, man-made slopes. These properties are conditioned by the location, land use, infrastructure, distribution and density of the population, organizational skills, etc. Risk management is becoming more discerning about the quality of the services offered, in addition to the compliance of the legal responsibility that the conception, design and construction of slopes and embankments in linear works require. (Pérez Soriano, 2014) This research project focuses on the identification and valuation of the risks derived from slopes and embankments in linear works of the Dominican Republic, in order to reduce the risks of failure of a slope or embankment and the number of victims who may be affected from it. Concluding with a risk catalogue for slopes and embankments in linear works.

Desarrollo de un Modelo de Gestión de Riesgos Operacionales basado en las Directrices del Comité de Basilea

La presente investigación tiene como objetivo principal diseñar un Modelo de Gestión de Riesgos Operacionales (MGRO) según las Directrices de los Acuerdos II y III del Comité de Supervisión Bancaria de Basilea del Banco de Pagos Internacionales (CSBB-BPI). Se considera importante realizar un estudio sobre este tema dado que son los riesgos operacionales (OpR) los responsables en gran medida de las últimas crisis financieras mundiales y por la dificultad para detectarlos en las organizaciones. Se ha planteado un modelo de gestión subdividido en dos vías de influencias. La primera acoge el paradigma holístico en el que se considera que hay múltiples maneras de percibir un proceso cíclico, así como las herramientas para observar, conocer y entender el objeto o sujeto percibido. La segunda vía la representa el paradigma totalizante, en el que se obtienen datos tanto cualitativos como cuantitativos, los cuales son complementarios entre si. Por otra parte, este trabajo plantea el diseño de un programa informático de OpR Cualitativo, que ha sido diseñado para determinar la raíz de los riesgos en las organizaciones y su Valor en Riesgo Operacional (OpVaR) basado en el método del indicador básico. Aplicando el ciclo holístico al caso de estudio, se obtuvo el siguiente diseño de investigación: no experimental, univariable, transversal descriptiva, contemporánea, retrospectiva, de fuente mixta, cualitativa (fenomenológica y etnográfica) y cuantitativa (descriptiva y analítica). La toma de decisiones y recolección de información se realizó en dos fases en la unidad de estudio. En la primera se tomó en cuenta la totalidad de la empresa Corpoelec-EDELCA, en la que se presentó un universo estadístico de 4271 personas, una población de 2390 personas y una unidad de muestreo de 87 personas. Se repitió el proceso en una segunda fase, para la Central Hidroeléctrica Simón Bolívar, y se determinó un segundo universo estadístico de 300 trabajadores, una población de 191 personas y una muestra de 58 profesionales. Como fuentes de recolección de información se utilizaron fuentes primarias y secundarias. Para recabar la información primaria se realizaron observaciones directas, dos encuestas para detectar las áreas y procesos con mayor nivel de riesgos y se diseñó un cuestionario combinado con otra encuesta (ad hoc) para establecer las estimaciones de frecuencia y severidad de pérdidas operacionales. La información de fuentes secundarias se extrajo de las bases de datos de Corpoelec-EDELCA, de la IEA, del Banco Mundial, del CSBB-BPI, de la UPM y de la UC at Berkeley, entre otras. Se establecieron las distribuciones de frecuencia y de severidad de pérdidas operacionales como las variables independientes y el OpVaR como la variable dependiente. No se realizó ningún tipo de seguimiento o control a las variables bajo análisis, ya que se consideraron estas para un instante especifico y solo se determinan con la finalidad de establecer la existencia y valoración puntual de los OpR en la unidad de estudio. El análisis cualitativo planteado en el MGRO, permitió detectar que en la unidad de investigación, el 67% de los OpR detectados provienen de dos fuentes principales: procesos (32%) y eventos externos (35%). Adicionalmente, la validación del MGRO en Corpoelec-EDELCA, permitió detectar que el 63% de los OpR en la organización provienen de tres categorías principales, siendo los fraudes externos los presentes con mayor regularidad y severidad de pérdidas en la organización. La exposición al riesgo se determinó fundamentándose en la adaptación del concepto de OpVaR que generalmente se utiliza para series temporales y que en el caso de estudio presenta la primicia de aplicarlo a datos cualitativos transformados con la escala Likert. La posibilidad de utilizar distribuciones de probabilidad típicas para datos cuantitativos en distribuciones de frecuencia y severidad de pérdidas con datos de origen cualitativo fueron analizadas. Para el 64% de los OpR estudiados se obtuvo que la frecuencia tiene un comportamiento semejante al de la distribución de probabilidad de Poisson y en un 55% de los casos para la severidad de pérdidas se obtuvo a las log-normal como las distribuciones de probabilidad más comunes, con lo que se concluyó que los enfoques sugeridos por el BCBS-BIS para series de tiempo son aplicables a los datos cualitativos. Obtenidas las distribuciones de frecuencia y severidad de pérdidas, se convolucionaron estas implementando el método de Montecarlo, con lo que se obtuvieron los enfoques de distribuciones de pérdidas (LDA) para cada uno de los OpR. El OpVaR se dedujo como lo sugiere el CSBB-BPI del percentil 99,9 o 99% de cada una de las LDA, obteniéndose que los OpR presentan un comportamiento similar al sistema financiero, resultando como los de mayor peligrosidad los que se ubican con baja frecuencia y alto impacto, por su dificultad para ser detectados y monitoreados. Finalmente, se considera que el MGRO permitirá a los agentes del mercado y sus grupos de interés conocer con efectividad, fiabilidad y eficiencia el status de sus entidades, lo que reducirá la incertidumbre de sus inversiones y les permitirá establecer una nueva cultura de gestión en sus organizaciones. ABSTRACT This research has as main objective the design of a Model for Operational Risk Management (MORM) according to the guidelines of Accords II and III of the Basel Committee on Banking Supervision of the Bank for International Settlements (BCBS- BIS). It is considered important to conduct a study on this issue since operational risks (OpR) are largely responsible for the recent world financial crisis and due to the difficulty in detecting them in organizations. A management model has been designed which is divided into two way of influences. The first supports the holistic paradigm in which it is considered that there are multiple ways of perceiving a cyclical process and contains the tools to observe, know and understand the subject or object perceived. The second way is the totalizing paradigm, in which both qualitative and quantitative data are obtained, which are complementary to each other. Moreover, this paper presents the design of qualitative OpR software which is designed to determine the root of risks in organizations and their Operational Value at Risk (OpVaR) based on the basic indicator approach. Applying the holistic cycle to the case study, the following research design was obtained: non- experimental, univariate, descriptive cross-sectional, contemporary, retrospective, mixed-source, qualitative (phenomenological and ethnographic) and quantitative (descriptive and analytical). Decision making and data collection was conducted in two phases in the study unit. The first took into account the totality of the Corpoelec-EDELCA company, which presented a statistical universe of 4271 individuals, a population of 2390 individuals and a sampling unit of 87 individuals. The process was repeated in a second phase to the Simon Bolivar Hydroelectric Power Plant, and a second statistical universe of 300 workers, a population of 191 people and a sample of 58 professionals was determined. As sources of information gathering primary and secondary sources were used. To obtain the primary information direct observations were conducted and two surveys to identify the areas and processes with higher risks were designed. A questionnaire was combined with an ad hoc survey to establish estimates of frequency and severity of operational losses was also considered. The secondary information was extracted from the databases of Corpoelec-EDELCA, IEA, the World Bank, the BCBS-BIS, UPM and UC at Berkeley, among others. The operational loss frequency distributions and the operational loss severity distributions were established as the independent variables and OpVaR as the dependent variable. No monitoring or control of the variables under analysis was performed, as these were considered for a specific time and are determined only for the purpose of establishing the existence and timely assessment of the OpR in the study unit. Qualitative analysis raised in the MORM made it possible to detect that in the research unit, 67% of detected OpR come from two main sources: external processes (32%) and external events (35%). Additionally, validation of the MORM in Corpoelec-EDELCA, enabled to estimate that 63% of OpR in the organization come from three main categories, with external fraud being present more regularly and greater severity of losses in the organization. Risk exposure is determined basing on adapting the concept of OpVaR generally used for time series and in the case study it presents the advantage of applying it to qualitative data transformed with the Likert scale. The possibility of using typical probability distributions for quantitative data in loss frequency and loss severity distributions with data of qualitative origin were analyzed. For the 64% of OpR studied it was found that the frequency has a similar behavior to that of the Poisson probability distribution and 55% of the cases for loss severity it was found that the log-normal were the most common probability distributions. It was concluded that the approach suggested by the BCBS-BIS for time series can be applied to qualitative data. Once obtained the distributions of loss frequency and severity have been obtained they were subjected to convolution implementing the Monte Carlo method. Thus the loss distribution approaches (LDA) were obtained for each of the OpR. The OpVaR was derived as suggested by the BCBS-BIS 99.9 percentile or 99% of each of the LDA. It was determined that the OpR exhibits a similar behavior to the financial system, being the most dangerous those with low frequency and high impact for their difficulty in being detected and monitored. Finally, it is considered that the MORM will allows market players and their stakeholders to know with effectiveness, efficiency and reliability the status of their entities, which will reduce the uncertainty of their investments and enable them to establish a new management culture in their organizations.

Modelo de estimación del contenido de carbono sobre el suelo en bosque húmedo tropical de la Amazonía Occidental usando tecnología LiDAR aerotransportada.

El modelo aplicado en el estudio, se basa en una metodología publicada en trabajos anteriores en bosques tropicales. El análisis se apoya en información obtenida con tecnología LiDAR aerotransportada y datos de censo de una mega-parcela forestal de 50 ha, en el Parque Nacional ?Yasuní? (PNY) en Ecuador. La metodología considera los hábitats topográficos, la distribución espacial de las parcelas de ajuste y validación en el diseño de muestreo y un modelo general para estimación de densidad de carbono. El modelo se ajusta usando técnicas de análisis de regresión lineal (R2=0,89, RMSE=8,67 MgC ha-1). Este enfoque metodológico hace posible calibraciones más precisas y permite cartografiar la distribución espacial de la biomasa sobre el suelo en alta resolución para una mejor estimación de las reservas de carbono en bosque tropical.

Respuesta fisiológica, productiva y cualitativa del cultivar Syrah (Vitis vinifera L) con diferentes frecuencias de riego y densidad de goteros en suelos arcillosos

La disponibilidad hídrica es uno de los principales factores que determinan el rendimiento del viñedo en muchas regiones vitícolas, por lo que sus consecuencias han sido ampliamente estudiadas. Sin embargo, para una cantidad de agua de riego determinada, otros aspectos como la frecuencia de aplicación, o la combinación entre el caudal de los goteros y la distancia entre los mismos (es decir, el patrón de distribución de agua en el suelo), pueden jugar un papel relevante, pero estos factores han sido poco estudiados. El objetivo de este trabajo ha sido evaluar las implicaciones agronómicas y fisiológicas de dos frecuencias de riego (IrrF, cada 2 y 4 días) y dos patrones de distribución de agua (DisP, goteros de 2 L h-1 separados 0,6 m vs. goteros de 4 L h-1 separados 1,2 m). El experimento se llevó a cabo durante cuatro temporadas consecutivas en un viñedo cv. Syrah con un suelo arcilloso en el centro de España, y los dos factores fueron evaluados bajo dos condiciones de disponibilidad hídrica (Baja: 20% de ETo y Media: 40% de ETo). El efecto de la frecuencia de riego y el patrón de distribución de agua en la respuesta agronómica del cv. Syrah se ha estudiado en el capítulo 1. La frecuencia de riego y el patrón de distribución de agua en el suelo afectaron a algunos aspectos de los componentes de rendimiento y desarrollo vegetativo en las dos condiciones de disponibilidad hídrica, aunque los efectos observados no fueron los mismos todos los años. Los efectos fueron más evidentes para IrrF en condiciones de baja disponibilidad hídrica y para DisP en condiciones de disponibilidad hídrica media. Dos de los cuatro años del experimento, el pasar de frecuencia de riego de 2 días a 4 días causó un incremento medio de rendimiento del 20% para la situación de baja disponibilidad hídrica. La textura del suelo, sin duda ha condicionado los resultados obtenidos en los tratamientos regados con el 20% de la ETo, ya que regar cada dos días implicaba la aplicación de pequeñas cantidades de agua y se formaban bulbos de riego superficiales, probablemente favoreciendo las pérdidas por evaporación. En el capítulo 2, se ha analizado el efecto de la frecuencia de riego y del patrón de distribución de agua en el estado hídrico de la planta y el intercambio gaseoso a nivel de hoja con el fin de explicar las diferencias observadas en la respuesta agronómica. En lo que respecta a la frecuencia de riego, en condiciones de baja disponibilidad hídrica, las plantas regadas cada 4 días (plantas 4d), mostraron mayores tasas de asimilación neta y conductancia estomática que las plantas regadas cada 2 días (plantas 2d), lo que es consistente con la hipótesis de que con la frecuencia de riego de 2 días se produjo una pérdida de eficiencia del uso del agua, probablemente debido a una mayor evaporación como consecuencia del hecho de que el volumen de suelo mojado creado era pequeño y cerca de la superficie. En condiciones de disponibilidad hídrica media, las diferencias en el intercambio gaseoso a nivel de hoja fueron mucho más pequeñas. Al comienzo del verano cada frecuencia de riego se comportó mejor uno de los días de medida, compensando al final del ciclo de riego de 4 días. Sin embargo, a medida que avanzó el verano y el déficit de agua se hizo más alto, las diferencias significativas aparecieron sólo en el 'día 4' del ciclo de riego, cuando las plantas 2d se comportaron mejor que las plantas regadas 4d que llevaban tres días sin regarse. Estas diferencias fisiológicas fueron menores que en condiciones de baja disponibilidad hídrica y al parecer no suficientes para afectar el comportamiento agronómico. En cuanto al patrón de distribución de agua, el efecto fue poco significativo, pero la mayor densidad de goteros tendió a presentar un mayor intercambio gaseoso a nivel de hoja, especialmente a media mañana. El efecto fue más importante para las condiciones de disponibilidad hídrica media. En el capítulo 3, se han comparado las relaciones entre el intercambio gaseoso a nivel de hoja, el estado hídrico y la demanda atmosférica, con el fin de explicar los cambios en la intensidad de la respuesta fisiológica observados en el Capítulo 2. No se han encontrado diferencias en dichas relaciones para el patrón de distribución de agua, por lo que sólo se ha analizado el efecto de la frecuencia de riego. El estudio se ha centrado fundamentalmente en si las plantas mostraron una respuesta fisiológica diferente a los cambios en el estado hídrico y en la demanda atmosférica según el tiempo transcurrido desde el último riego. Las diferencias observadas explican los resultados obtenidos en los capítulos anteriores, y sugieren la existencia de procesos de aclimatación vinculados a la frecuencia de riego y a la disponibilidad hídrica. Las plantas bajo condiciones de baja disponibilidad hídrica se mostraron más aclimatadas al estrés hídrico que aquellas en condiciones de disponibilidad hídrica media. La frecuencia de riego afectó claramente la relación entre los parámetros de intercambio gaseoso a nivel de hoja, el estado hídrico de la planta y las condiciones atmosféricas, y junto con la cantidad de agua aplicada tuvo implicaciones en el desarrollo de mecanismos de aclimatación que afectaron a la respuesta fisiológica de la planta, afectando a la eficiencia del riego. ABSTRACT Water availability is one of the major factors that determine vineyard performance in many grape growing regions, so its implications have been widely studied before. However, for a given irrigation water amount, other aspects such as application frequency, or emitter spacing and flow rate (i.e., distribution pattern), may play a relevant role, but these factors have been scarcely studied. The aim of this work was to evaluate the agronomic and physiological implications of two irrigation frequencies (IrrF, every 2 and 4 days) and two water distribution patterns (DisP, 2 L h−1 emitters every 0.6 m vs. 4 L h−1 emitters every 1.2 m). The experiment was carried out during four consecutive seasons in a cv. Syrah vineyard with a clay soil in central Spain, and the two factors were evaluated under two water availability conditions (LOW WA: 20% of ETo and MEDIUM WA: 40% of ETo). The effect of irrigation frequency and water distribution pattern on the agronomical response of cv. Syrah was studied in Chapter 1. IrrF and DisP affected some aspects of vegetative development and yield components under both water availability conditions, although the effects observed were not the same every year. The effects were more evident for IrrF under low water availability and for DisP under medium water availability. Two out of the four years of the experiment, the change of irrigation frequency from 2 days to 4 days promoted an average yield increase of 20% for the LOW WA situation. Soil texture certainly conditioned the results obtained under LOW WA conditions, since high frequency irrigation implied applying small amounts of water that resulted in limited superficial water bulbs, which probably favored water evaporation. In Chapter 2, the effect of irrigation frequency and water distribution pattern on plant water status and leaf gas exchange was analyzed to explain the differences observed in the agronomical response. Concerning irrigation frequency, under LOW WA conditions, applying irrigation every 4 days, resulted in higher net assimilation rates and stomatal conductance than doing it every 2 days, supporting the hypothesis that the latter frequency resulted in a water use efficiency loss, probably due to higher evaporation as a consequence of the fact the wetted soil volume created was small and close to the surface. Under MEDIUM WA conditions, differences in leaf gas exchange were much smaller. At the beginning of the summer each irrigation frequency behaved better one of the measurements days, compensating at the end of the 4-day irrigation cycle. However, as the summer progressed and water deficit became higher, significant differences appeared only on ‘day 4’ of the irrigation cycle, when 2d plants behaved better than 4d plants. These physiological differences were smaller than under LOW WA conditions and apparently not sufficient to affect agronomical performance. Regarding water distribution pattern, the effect was less significant but the closest emitter spacing resulted in general terms in a higher leaf gas exchange, especially at midmorning. The effect was more noticeable for MEDIUM WA conditions. In Chapter 3, the relationships between leaf gas exchange and leaf water status and atmospheric demand were compared to explain the changes in the intensity of the physiological response observed in Chapter 2. No differences were found in the relationships for water distribution pattern, so only the effect of irrigation frequency was analyzed focusing on whether the plants have a different physiological response to changes in water status and atmospheric demand according to the time elapsed since the last irrigation. Differences observed in the relationships explained the results obtained in the previous chapters, and point at the occurrence of acclimation processes linked to irrigation frequency and to water availability. Plants under LOW WATER AVAILABILITY conditions seemed to be more acclimated to water stress than those under MEDIUM WATER AVAILABILITY conditions. Irrigation frequency clearly affected the relationship between leaf gas exchange parameters, plant water status and atmospheric conditions, and together with the amount of water applied had implications in the development of acclimation mechanisms that affected plant physiological response, thus affecting irrigation efficiency.

Geografía del carbono en alta resolución en bosque tropical amazónico del Ecuador mediante sensores aerotransportados

La estimación de la biomasa de la vegetación terrestre en bosque tropical no sólo es un área de investigación en rápida expansión, sino también es un tema de gran interés para reducir las emisiones de carbono asociadas a la deforestación y la degradación forestal (REDD+). Las estimaciones de densidad de carbono sobre el suelo (ACD) en base a inventarios de campo y datos provenientes de sensores aerotransportados, en especial con sensores LiDAR, han conducido a un progreso sustancial en el cartografiado a gran escala de las reservas de carbono forestal. Sin embargo, estos mapas de carbono tienen incertidumbres considerables, asociadas generalmente al proceso de calibración del modelo de regresión utilizado para producir los mapas. En esta tesis se establece una metodología para la calibración y validación de un modelo general de estimación de ACD usando LiDAR en un sector del Parque Nacional Yasuní en Ecuador. En el proceso de calibración del modelo se considera el tamaño y la ubicación de las parcelas, la influencia de la topografía y la distribución espacial de la biomasa. Para el análisis de los datos se utilizan técnicas geoestadísticas en combinación con variables geomorfométricas derivadas de datos LiDAR, y se propone un esquema de muestreo estratificado por posiciones topográficas (valle, ladera y cima). La validación del modelo general para toda la zona de estudio presentó valores de RMSE = 5.81 Mg C ha-1, R2 = 0.94 y sesgo = 0.59, mientras que, al considerar las posiciones topográficas, el modelo presentó valores de RMSE = 1.67 Mg C ha-1, R2 = 0.98 y sesgo = 0.23 para el valle; RMSE = 3.13 Mg C ha-1, R2 = 0.98 y sesgo = - 0.34 para la ladera; y RMSE = 2.33 Mg C ha-1, R2 = 0.97 y sesgo = 0.74 para la cima. Los resultados obtenidos demuestran que la metodología de muestreo estratificado por posiciones topográficas propuesto, permite calibrar de manera efectiva el modelo general con las estimaciones de ACD en campo, logrando reducir el RMSE y el sesgo. Los resultados muestran el potencial de los datos LiDAR para caracterizar la estructura vertical de la vegetación en un bosque altamente diverso, permitiendo realizar estimaciones precisas de ACD, y conocer patrones espaciales continuos de la distribución de la biomasa aérea y del contenido de carbono en la zona de estudio. ABSTRACT Estimating biomass of terrestrial vegetation in tropical forest is not only a rapidly expanding research area, but also a subject of tremendous interest for reducing carbon emissions associated with deforestation and forest degradation (REDD+). The aboveground carbon density estimates (ACD) based on field inventories and airborne sensors, especially LiDAR sensors have led to a substantial progress in large-scale mapping of forest carbon stocks. However, these carbon maps have considerable uncertainties generally associated with the calibration of the regression model used to produce these maps. This thesis establishes a methodology for calibrating and validating a general ACD estimation model using LiDAR in Ecuador´s Yasuní National Park. The size and location of the plots are considered in the model calibration phase as well as the influence of topography and spatial distribution of biomass. Geostatistical analysis techniques are used in combination with geomorphometrics variables derived from LiDAR data, and then a stratified sampling scheme considering topographic positions (valley, slope and ridge) is proposed. The validation of the general model for the study area showed values of RMSE = 5.81 Mg C ha-1, R2 = 0.94 and bias = 0.59, while considering the topographical positions, the model showed values of RMSE = 1.67 Mg C ha-1, R2 = 0.98 and bias = 0.23 for the valley; RMSE = 3.13 Mg C ha-1, R2 = 0.98 and bias = - 0.34 for the slope; and RMSE = 2.33 Mg C ha-1, R2 = 0.97 and bias = 0.74 for the ridge. The results show that the stratified sampling methodology taking into account topographic positions, effectively calibrates the general model with field estimates of ACD, reducing RMSE and bias. The results show the potential of LiDAR data to characterize the vertical structure of vegetation in a highly diverse forest, allowing accurate estimates of ACD, and knowing continuous spatial patterns of biomass distribution and carbon stocks in the study area.

Fiabilidad estructural de coníferas españolas

Se ha procedido al análisis estadístico de un banco de datos de algo más de nueve mil probetas de tamaño estructural ensayadas a lo largo de la década final del s. XX y la inicial del s. XXI, provenientes de masas forestales de coníferas españolas. Las especies estudiadas (pinus sylsvestris, p. halepensis, p. pinaster, p. radiata, p. nigra) representan más del 80% de la superficie forestal española cubierta por coníferas, y ninguna especie de relevancia por volumen potencial de corta a medio plazo ha quedado excluida (quizá con la única excepción de p. pinea, y, con matices juniperus thurifera). Por lo tanto, puede considerarse una información razonablemente representativa de la población de coníferas españolas (en adelante, población objetivo), si bien los procedimientos de muestreo presentan marcadas diferencias en función de las especies. Los procedimientos de ensayo se han atenido, esencialmente, a la normativa europea de referencia, así como a la normativa española de clasificación de madera aserrada, básicamente UNE 56544, con los matices propios de las sucesivas ediciones de las indicadas normas durante un período de unos veinte años. El estudio se ha realizado con el objeto de caracterizar el estado de naturaleza de los principales parámetros resistentes la población objetivo (densidad, módulo de elasticidad longitudinal, y tensión de rotura en flexión pura), así como identificar posibles inconsistencias en la normativa vigente relacionada. Los resultados se han elaborado de forma operativa para su utilización, en una primera aproximación, en la ejecución de cálculos directos de fiabilidad estructural, tal como se plantean en las normativas ISO 2394, ISO 12491, ISO 13822, UNE EN 1990, Código Técnico de la Edificación español y Código Modelo Probabilístico. Inequívocamente, las variables resistentes de referencia de las coníferas españolas susceptibles de uso estructural presentan una variabilidad significativamente superior a la habitual en coníferas de origen en el Centro y Norte de Europa, así como en Norteamérica (en adelante, poblaciones de comparación). Esta diferencia es extrema en el caso de la resistencia, importante en el caso de la rigidez, y prácticamente testimonial en el caso de la densidad. Por otra parte, la rigidez y la densidad de la población objetivo es mucho mayor, para el mismo valor de la resistencia, que las de las poblaciones de comparación. Asimismo, las correlaciones entre las tres variables básicas son inferiores en la población objetivo respecto a la de comparación. Estos hechos llevan a que la aplicación del sistema de clases resistentes en uso conduzca a situaciones de dimensionado poco racionales, y marcadamente antieconómicas. Si bien el objeto del estudio no es el establecimiento de una “jerarquía” estructural de especies, sí deben subrayarse algunos aspectos particulares por especies. El p. nigra tiene los mayores valores centrales y los menores de dispersión de toda la población objetivo, y mayores valores centrales que cualquier referencia de la población de comparación: para la misma fiabilidad objetivo en supuestos habituales, en p.nigra puede ser necesario en torno a un 20% menos de volumen de madera que en p. sylvestris, y en torno a un 40% que en p. radiata . El p. radiata y el p. pinaster presentan los menores valores de centralidad y los mayores de dispersión, de toda la muestra objetivo en formatos aserrados (con la excepción de la densidad: el p. radiata presenta valores muy bajos, pero atípicamente baja variabilidad). Especies habitualmente “no consideradas” como el p. pinea, presentan valores notablemente altos de centralidad, especialmente en la variable densidad. Otro caso de interés es p. halepensis, con prestaciones estructurales muy elevadas, teniendo en cuenta la proveniencia de masas carentes de un manejo adecuado (a efectos de producción maderable). Al contrario que en las poblaciones de comparación, la densidad de las coníferas españolas se representa claramente por distribuciones lognormales, frente a las distribuciones normales que se consideran en la literatura para la misma variable básica. Esto no tiene consecuencias relevantes en la fiabilidad estructural de las uniones, debido a los bajos coeficientes de variación (entre el 10 y 13%, típicamente). Para la tensión de rotura en flexión, parece haber una marcada diferencia entre el tipo de distribución que representa a las poblaciones no clasificadas o rechazadas (que es claramente Lognormal) y la que representa a las poblaciones clasificadas (que parece ser de tipo Weibull de tres parámetros). La distinción va más allá de la sutileza matemática: si la distribución fuese lognormal para todos los casos, la ineficiencia de la aplicación del sistema de clases resistentes aún sería más acusada que lo anteriormente indicado. La distribución normal, sólo representa adecuadamente el estado de naturaleza del módulo de elasticidad y la tensión de rotura de maderas no escuadradas (cilindradas, descortezadas, o simplemente desramadas), y del módulo de elasticidad de las muestras rechazadas. Las coníferas no clasificadas, arrojan valores relativamente elevados de resistencia, lo que hace pensar en la conveniencia de la definición de un “protocolo de mínimos” para posibilitar su uso en el marco normativo actual. Se ha ejemplificado cómo de la aplicación de técnicas probabilísticas utilizando como fundamento los resultados del análisis realizado, resulta un ahorro material potencialmente notable en la población objetivo (en la horquilla del 30 al 60%, dependiendo del estado límite que gobierne el dimensionamiento), respecto a la aplicación directa de los protocolos normativos a partir de la aplicación de la norma de clasificación española. Complementariamente, se han elaborado modelos de variabilidad de los parámetros de las distribuciones propuestas, al objeto de facilitar la aplicación de técnicas de actualización bayesiana. Estas técnicas tienen particular relevancia en el entorno español, por dos razones. Una de ellas, es que suponen, de hecho, la única manera de tratar la cuestión de la seguridad de estructuras prexistentes en un marco matemático-formal consistente con la normativa. La otra, es que permiten la utilización de maderas procedentes de medios de producción con niveles de control de calidad mínimos, incluso casi inexistentes, con inversiones en coste de ensayos comparativamente irrelevantes. ABSTRACT A group of databases making a total of ca. 9.000 beams produced by sawing coniferous species from Spain, tested according EN 384 (1) and visually graded according Spanish standards were analyzed. The main goal was to reach a detailed comprehension about the true state of nature of the three reference structural properties (density, MOR and MOE), in the conditions the material reaches the real market. This goal has been addressed through the detailed adjustment of the proper probabilistic model to the data available. The most important outcome has been the proposal of a set of distribution parameters for the target species, ready to be used in direct probabilistic calculation. This use is exemplified through structural reliability calculations, establishing a priori data for Bayesian reliability updating, Spanish code calibration suggestions, and parameters proposal for the Probabilistic Model Code (2).

Deforestación y Pérdida de hábitat en Bosques de montaña en la Cuenca alta del Río Zamora (Loja, Ecuador)

Los bosques húmedos de montaña se encuentran reconocidos como uno de los ecosistemas más amenazados en el mundo, llegando inclusive a ser considerado como un “hotspot” por su alta diversidad y endemismo. La acelerada pérdida de cobertura vegetal de estos bosques ha ocasionado que, en la actualidad, se encuentren restringidos a una pequeña fracción de su área de distribución histórica. Pese a esto, los estudios realizados sobre cual es efecto de la deforestación, fragmentación, cambios de uso de suelo y su efecto en las comunidades de plantas presentes en este tipo de vegetación aún son muy escuetos, en comparación a los realizados con sus similares amazónicos. En este trabajo, el cual se encuentra dividido en seis capítulos, abordaremos los siguientes objetivos: a) Comprender cuál es la dinámica que han seguido los diferentes tipos de bosques montanos andinos de la cuenca del Rio Zamora, Sur de Ecuador durante entre 1976 y 2002. b) Proveer de evidencia de las tasas de deforestación y fragmentación de todos los tipos diferentes de bosques montanos andinos presentes en la cuenca del Rio Zamora, Sur de Ecuador entre 1976 y 2002. c) Determinar qué factores inducen a la fragmentación de bosques de montaña en la cuenca alta del río Zamora entre 1976 y 2002. d) Determinar cuáles son y cómo afectan los factores ambientales y socioeconómicos a la dinámica de la deforestación y regeneración (pérdida y recuperación del hábitat) sufrida por los bosques de montaña dentro de la zona de estudio y e) Determinar si la deforestación y fragmentación actúan sobre la diversidad y estructura de las comunidades de tres tipos de organismos (comunidades de árboles, comunidades de líquenes epífitos y comunidades de hepáticas epífitas). Este estudio se centró en el cuenca alta del río Zamora, localizada al sur de Ecuador entre las coordenadas 3º 00´ 53” a 4º 20´ 24.65” de latitud sur y 79º 49´58” a 78º 35´ 38” de longitud oeste, que cubre alrededor de 4300 km2 de territorio situado entre las capitales de las provincias de Loja y Zamora-Chinchipe. Con objeto de predecir la dinámica futura de la deforestación en la región de Loja y cómo se verán afectados los diferentes tipos de hábitat, así como para detectar los factores que más influyen en dicha dinámica, se han construido modelos basados en la historia de la deforestación derivados de fotografías aéreas e imágenes satelitales de tres fechas (1976, 1989 y 2002). La cuantificación de la deforestación se realizó mediante la tasa de interés compuesto y para la caracterización de la configuración espacial de los fragmentos de bosque nativo se calcularon índices de paisaje los cuales fueron calculados utilizando el programa Fragstats 3.3. Se ha clasificado el recubrimiento del terreno en forestal y no forestal y se ha modelado su evolución temporal con Modelos Lineales Generalizados Mixtos (GLMM), empleando como variables explicativas tanto variables ambientales espacialmente explícitas (altitud, orientación, pendiente, etc) como antrópicas (distancia a zonas urbanizadas, deforestadas, caminos, entre otras). Para medir el efecto de la deforestación sobre las comunidades modelo (de árboles, líquenes y hepáticas) se monitorearon 11 fragmentos de vegetación de distinto tamaño: dos fragmentos de más de cien hectáreas, tres fragmentos de entre diez y noventa ha y seis fragmentos de menos de diez hectáreas. En ellos se instalaron un total de 38 transectos y 113 cuadrantes de 20 x 20 m a distancias que se alejaban progresivamente del borde en 10, 40 y 80 m. Nuestros resultados muestran una tasa media anual de deforestación del 1,16% para todo el período de estudio, que el tipo de vegetación que más alta tasa de destrucción ha sufrido, es el páramo herbáceo, con un 2,45% anual. El análisis de los patrones de fragmentación determinó un aumento en 2002 de más del doble de fragmentos presentes en 1976, lo cual se repite en el análisis del índice de densidad promedio. El índice de proximidad media entre fragmentos muestra una reducción progresiva de la continuidad de las áreas forestadas. Si bien las formas de los fragmentos se han mantenido bastante similares a lo largo del período de estudio, la conectividad entre estos ha disminuido en un 84%. Por otro lado, de nuestros análisis se desprende que las zonas con mayor probabilidad de deforestarse son aquellas que están cercanas a zonas previamente deforestadas; la cercanía a las vías también influye significativamente en la deforestación, causando un efecto directo en la composición y estructura de las comunidades estudiadas, que en el caso de los árboles viene mediado por el tamaño del fragmento y en el caso del componente epífito (hepáticas y líquenes), viene mediado tanto por el tamaño del fragmento como por la distancia al borde del mismo. Se concluye la posibilidad de que, de mantenerse esta tendencia, este tipo de bosques desaparecerá en corto tiempo y los servicios ecosistémicos que prestan, se verán seriamente comprometidos. ABSTRACT Mountain rainforests are recognized as one of the most threatened ecosystems in the world, and have even come to be considered as a “hotspot” due to their high degree of diversity and endemism. The accelerated loss of plant cover of these forests has caused them to be restricted today to a small fraction of their area of historic distribution. In spite of this, studies done on the effect of deforestation, fragmentation, changes in soil use and their effect on the plant communities present in this type of vegetation are very brief compared to those done on their analogues in the Amazon region. In this study, which is divided into six chapters, we will address the following objectives: a) To understand what the dynamic followed by the different types of Andean mountain forests in the Zamora River watershed of southern Ecuador has been between 1976 and 2002. b) To provide evidence of the rates of deforestation and fragmentation of all the different types of Andean mountain forests existing in the upper watershed of the Zamora River between 1976 and 2002. c) To determine the factors that induces fragmentation of all different types of Andean mountain forests existing in the upper watershed of the Zamora River between 1976 and 2002. d) To determine what the environmental and anthropogenic factors are driving the dynamic of deforestation and regeneration (loss and recuperation of the habitat) suffered by the mountain forests in the area of the study and e) To determine if the deforestation and fragmentation act upon the diversity and structure of three model communities: trees, epiphytic lichens and epiphytic liverworts. This study is centered on the upper Zamora River watershed, located in southern Ecuador between 3º 00´ 53” and 4º 20´ 24.65 south latitude and 79º 49´ 58” to 78º 35´ 38” west longitude, and covers around 4,300 km2 of territory located between Loja and Zamora-Chinchipe provinces. For the purpose of predicting the future dynamic of deforestation in the Loja region and how different types of habitats will be affected, as well as detecting the environmental and socioeconomic factors that influence landscape dynamics, models were constructed based on deforestation history, derived from aerial photographs and satellite images for three dates (1976, 1989 and 2002). Quantifying the deforestation was done using the compound interest rate; to characterize the spatial configuration of fragments of native forest, landscape indices were calculated with Fragstats 3.3 program. Land cover was classified as forested and not forested and its evolution over time was modeled with Generalized Linear Mixed Models (GLMM), using spatially explicit environmental variables (altitude, orientation, slope, etc.) as well as anthropic variables (distance to urbanized, deforested areas and roads, among others) as explanatory variables. To measure the effects of fragmentation on three types of model communities (forest trees and epiphytic lichen and liverworts), 11 vegetation fragments of different sizes were monitored: two fragments of more than one hundred hectares, three fragments of between ten and ninety ha and six fragments of fewer than ten hectares . In these fragments, a total of 38 transects and 113 20 x 20 m quadrats were installed at distances that progressively moved away from the edge of the fragment by 10, 40 and 80 m. Our results show an average annual rate of deforestation of 1.16% for the entire period of the study, and that the type of vegetation that suffered the highest rate of destruction was grassy paramo, with an annual rate of 2.45%. The analysis of fragmentation patterns determined the number of fragments in 2002 more than doubled the number of fragments present in 1976, and the same occurred for the average density index. The variation of the average proximity index among fragments showed a progressive reduction of the continuity of forested areas. Although fragment shapes have remained quite similar over the period of the study, connectivity among them has diminished by 84%. On the other hand, it emerged from our analysis that the areas of greatest probability of deforestation were those that are close to previously deforested areas; proximity to roads also significantly favored the deforestation causing a direct effect on the composition of our model communities, that in the case of forest trees is determined by the size of the fragment, and in the case of the epiphyte communities (liverworts and lichens), is determined, by the size of the fragment as well as the distance to edge. A subject under discussion is the possibility that if this tendency continues, this type of forest will disappear in a short time, and the ecological services it provides, will be seriously endangered.

Optimización de operaciones unitarias utilizando líquidos iónicos como agentes materiales de separación

En este trabajo se estudia la modelización y optimización de procesos industriales de separación mediante el empleo de mezclas de líquidos iónicos como disolventes. Los disolventes habitualmente empleados en procesos de absorción o extracción suelen ser componentes orgánicos muy volátiles y dañinos para la salud humana. Las innovadoras propiedades que presentan los líquidos iónicos, los convierten en alternativas adecuadas para solucionar estos problemas. La presión de vapor de estos compuestos es muy baja y apenas varía con la temperatura. Por tanto, estos compuestos apenas se evaporan incluso a temperaturas altas. Esto supone una gran ventaja en cuanto al empleo de estos compuestos como disolventes industriales ya que permite el reciclaje continuo del disolvente al final del proceso sin necesidad de introducir disolvente fresco debido a la evaporación del mismo. Además, al no evaporarse, estos compuestos no suponen un peligro para la salud humana por inhalación; al contrario que otros disolventes como el benceno. El único peligro para la salud que tienen estos compuestos es por tanto el de contacto directo o ingesta, aunque de hecho muchos Líquidos Iónicos son inocuos con lo cual no existe peligro para la salud ni siquiera a través de estas vías. Los procesos de separación estudiados en este trabajo, se rigen por la termodinámica de fases, concretamente el equilibrio líquido-vapor. Para la predicción de los equilibrios se ha optado por el empleo de modelos COSMO (COnductor-like Screening MOdel). Estos modelos tienen su origen en el empleo de la termodinámica de solvatación y en la mecánica cuántica. En el desarrollo de procesos y productos, químicos e ingenieros frecuentemente precisan de la realización de cálculos de predicción de equilibrios de fase. Previamente al desarrollo de los modelos COSMO, se usaban métodos de contribución de grupos como UNIFAC o modelos de coeficientes de actividad como NRTL.La desventaja de estos métodos, es que requieren parámetros de interacción binaria que únicamente pueden obtenerse mediante ajustes por regresión a partir de resultados experimentales. Debido a esto, estos métodos apenas tienen aplicabilidad para compuestos con grupos funcionales novedosos debido a que no se dispone de datos experimentales para llevar a cabo los ajustes por regresión correspondientes. Una alternativa a estos métodos, es el empleo de modelos de solvatación basados en la química cuántica para caracterizar las interacciones moleculares y tener en cuenta la no idealidad de la fase líquida. Los modelos COSMO, permiten la predicción de equilibrios sin la necesidad de ajustes por regresión a partir de resultados experimentales. Debido a la falta de resultados experimentales de equilibrios líquido-vapor de mezclas en las que se ven involucrados los líquidos iónicos, el empleo de modelos COSMO es una buena alternativa para la predicción de equilibrios de mezclas con este tipo de materiales. Los modelos COSMO emplean las distribuciones superficiales de carga polarizada (sigma profiles) de los compuestos involucrados en la mezcla estudiada para la predicción de los coeficientes de actividad de la misma, definiéndose el sigma profile de una molécula como la distribución de probabilidad de densidad de carga superficial de dicha molécula. Dos de estos modelos son COSMO-RS (Realistic Solvation) y COSMO-SAC (Segment Activity Coefficient). El modelo COSMO-RS fue la primera extensión de los modelos de solvatación basados en continuos dieléctricos a la termodinámica de fases líquidas mientras que el modelo COSMO-SAC es una variación de este modelo, tal y como se explicará posteriormente. Concretamente en este trabajo se ha empleado el modelo COSMO-SAC para el cálculo de los coeficientes de actividad de las mezclas estudiadas. Los sigma profiles de los líquidos iónicos se han obtenido mediante el empleo del software de química computacional Turbomole y el paquete químico-cuántico COSMOtherm. El software Turbomole permite optimizar la geometría de la molécula para hallar la configuración más estable mientras que el paquete COSMOtherm permite la obtención del perfil sigma del compuesto mediante el empleo de los datos proporcionados por Turbomole. Por otra parte, los sigma profiles del resto de componentes se han obtenido de la base de datos Virginia Tech-2005 Sigma Profile Database. Para la predicción del equilibrio a partir de los coeficientes de actividad se ha empleado la Ley de Raoult modificada. Se ha supuesto por tanto que la fracción de cada componente en el vapor es proporcional a la fracción del mismo componente en el líquido, dónde la constante de proporcionalidad es el coeficiente de actividad del componente en la mezcla multiplicado por la presión de vapor del componente y dividido por la presión del sistema. Las presiones de vapor de los componentes se han obtenido aplicando la Ley de Antoine. Esta ecuación describe la relación entre la temperatura y la presión de vapor y se deduce a partir de la ecuación de Clausius-Clapeyron. Todos estos datos se han empleado para la modelización de una separación flash usando el algoritmo de Rachford-Rice. El valor de este modelo reside en la deducción de una función que relaciona las constantes de equilibrio, composición total y fracción de vapor. Para llevar a cabo la implementación del modelado matemático descrito, se ha programado un código empleando el software MATLAB de análisis numérico. Para comprobar la fiabilidad del código programado, se compararon los resultados obtenidos en la predicción de equilibrios de mezclas mediante el código con los resultados obtenidos mediante el simulador ASPEN PLUS de procesos químicos. Debido a la falta de datos relativos a líquidos iónicos en la base de datos de ASPEN PLUS, se han introducido estos componentes como pseudocomponentes, de manera que se han introducido únicamente los datos necesarios de estos componentes para realizar las simulaciones. El modelo COSMO-SAC se encuentra implementado en ASPEN PLUS, de manera que introduciendo los sigma profiles, los volúmenes de la cavidad y las presiones de vapor de los líquidos iónicos, es posible predecir equilibrios líquido-vapor en los que se ven implicados este tipo de materiales. De esta manera pueden compararse los resultados obtenidos con ASPEN PLUS y como el código programado en MATLAB y comprobar la fiabilidad del mismo. El objetivo principal del presente Trabajo Fin de Máster es la optimización de mezclas multicomponente de líquidos iónicos para maximizar la eficiencia de procesos de separación y minimizar los costes de los mismos. La estructura de este problema es la de un problema de optimización no lineal con variables discretas y continuas, es decir, un problema de optimización MINLP (Mixed Integer Non-Linear Programming). Tal y como se verá posteriormente, el modelo matemático de este problema es no lineal. Por otra parte, las variables del mismo son tanto continuas como binarias. Las variables continuas se corresponden con las fracciones molares de los líquidos iónicos presentes en las mezclas y con el caudal de la mezcla de líquidos iónicos. Por otra parte, también se ha introducido un número de variables binarias igual al número de líquidos iónicos presentes en la mezcla. Cada una de estas variables multiplican a las fracciones molares de sus correspondientes líquidos iónicos, de manera que cuando dicha variable es igual a 1, el líquido se encuentra en la mezcla mientras que cuando dicha variable es igual a 0, el líquido iónico no se encuentra presente en dicha mezcla. El empleo de este tipo de variables obliga por tanto a emplear algoritmos para la resolución de problemas de optimización MINLP ya que si todas las variables fueran continuas, bastaría con el empleo de algoritmos para la resolución de problemas de optimización NLP (Non-Linear Programming). Se han probado por tanto diversos algoritmos presentes en el paquete OPTI Toolbox de MATLAB para comprobar cuál es el más adecuado para abordar este problema. Finalmente, una vez validado el código programado, se han optimizado diversas mezclas de líquidos iónicos para lograr la máxima recuperación de compuestos aromáticos en un proceso de absorción de mezclas orgánicas. También se ha usado este código para la minimización del coste correspondiente a la compra de los líquidos iónicos de la mezcla de disolventes empleada en la operación de absorción. En este caso ha sido necesaria la introducción de restricciones relativas a la recuperación de aromáticos en la fase líquida o a la pureza de la mezcla obtenida una vez separada la mezcla de líquidos iónicos. Se han modelizado los dos problemas descritos previamente (maximización de la recuperación de Benceno y minimización del coste de operación) empleando tanto únicamente variables continuas (correspondientes a las fracciones o cantidades molares de los líquidos iónicos) como variables continuas y binarias (correspondientes a cada uno de los líquidos iónicos implicados en las mezclas).

New systems of management to improve the health and welfare of rabbits in farm conditions = Nuevos sistemas de cría y manejo para mejorar la salud y bienestar de conejos de granja

El objetivo general de esta Tesis Doctoral fue evaluar nuevos sistemas de alojamiento y cría de conejos de granja, estudiando tanto parámetros comportamentales (experimento 1) como productivos y reproductivos (experimento 3). Además, se evaluaron diferentes técnicas de muestreo con el fin de optimizar el tiempo empleado para el estudio del comportamiento animal (experimento 2). En el experimento 1, se estudió el comportamiento de conejas alojadas en dos tipos de jaulas (TJ), convencionales vs. alternativas con una plataforma elevada, en distintos estados fisiológicos (EF), lactantes y gestantes. Se observó el comportamiento de 12 conejas reproductoras con grabaciones de una duración de 24 h continuas. Independientemente del EF y TJ, las conejas pasaron gran parte de su tiempo sobre el reposapatas (57,7 %, de media). Sin embargo, debido al uso de la plataforma (23,0% del tiempo, de media), las conejas lactantes permanecieron un 36,6 % menos de tiempo (P<0,001) sobre el reposapatas y las gestantes un 27,0% menos (P<0,001) sobre el enrejillado en jaulas alternativas que en convencionales. En las jaulas alternativas, las conejas podían adoptar la postura “levantada”, sin embargo ésta fue observada solamente en conejas gestantes una media de 4,6 veces al día. Las conejas bebieron con mas frecuencia en jaulas convencionales que en alternativas (24,6 vs 19,1 veces al día; P<0,05). Se observó una mayor duración y frecuencia del comportamiento “interactuando con compañeras” en conejas gestantes alojadas en jaulas convencionales (276 s/d y 4,6 veces/d; P<0,05). La frecuencia de “interactuando con gazapos” fue menor en jaulas alternativas que en convencionales (2,4 vs 8,6 veces al día; P<0,01). La hora del día afectó al comportamiento de las conejas, teniendo un comportamiento menos activo durante las horas centrales del día. Durante las horas de oscuridad, las conejas estuvieron más inquietas realizando comportamientos como ‘encabritarse’ o amamantar, coincidiendo éstos en el tiempo en el cual las conejas pasaron más tiempo en la plataforma. Las conejas utilizaron frecuentemente la plataforma, independientemente del estado fisiológico. En la fase de lactación, las conejas utilizaron la plataforma para huir de los intentos de mamar por parte de los gazapos cuando éstas no estaban receptivas. El uso de la plataforma puede dar lugar a problemas higiénicos debidos tanto por la acumulación de heces sobre ella como por la caída de heces y orina sobre los animales que están en la parte inferior. La ausencia de estereotipias por parte de las conejas tanto en jaulas alternativas como en convencionales no sugiere una falta de bienestar debida al sistema de alojamiento. En el experimento 2, se compararon distintos métodos de observación simplificada con respecto un método de referencia usando grabaciones continuas de 24 h para la evaluación del comportamiento de conejas en distintos estados fisiológicos (gestantes y lactantes) alojadas en dos tipos de jaulas (convencionales y alternativas). Se analizaron un total de 576 h de grabaciones continuas de 24 h en 12 conejas reproductoras al final del periodo de lactación y en las mismas conejas después del destete. Los comportamientos observados se clasificaron en tres categorías independientes (localización en la jaula, postura y comportamientos funcionales). Se utilizaron grabaciones continuas de 24 h como método de referencia para validar otros cuatro métodos de observación simplificados, utilizando grabaciones de distinta duración y frecuencia a lo largo del día. Métodos regulares: corto y largo con 2.4 y 8 h de observación respectivamente, y métodos irregulares: corto y largo con 6 y 8 h de observación, respectivamente. Como resultado, se observó que independientemente del sistema de alojamiento, el mejor método para reducir el tiempo de observación necesario para evaluar el comportamiento de conejas reproductoras depende del tipo de variable a estudiar y del estado fisiológico de las conejas. En gestantes, los métodos irregulares no fueron adecuados para estimar comportamientos de larga duración tales como tumbada, sentada, descansando y acicalándose. Sin embargo, en ambos estados fisiológicos, los métodos regulares fueron precisos para los comportamientos de los grupos localización y postura y para comportamientos funcionales de larga duración. Por otro lado, los coeficientes de variación de los comportamientos poco frecuentes realizados principalmente durante el periodo de oscuridad fueron muy altos, y el método irregular largo obtuvo los menores errores de estimación para éstos comportamientos. En el experimento 3, se estudió el efecto de un uso combinado de lactaciones largas (hasta 46 días) con jaulas alternativas sobre los parámetros productivos y reproductivos de 104 conejas y sus camadas durante cinco ciclos reproductivos. La mitad de las conejas fueron alojadas en jaulas polivalentes convencionales (39 cm x 100 cm x 30 cm) y la otra mitad en jaulas polivalentes alternativas (39 cm x 100 cm x 60 cm), con una plataforma elevada. Dentro de cada grupo de alojamiento, la mitad de las conejas se destetaron a 32 días y la otra mitad a 46 días tras el parto. Las lactaciones más largas afectaron negativamente al peso (P<0,001), contenido en grasa y energía (P<0,05) de las conejas al final del periodo de lactación, pero éste efecto disminuyó con el número de partos. La fertilidad, prolificidad y la mortalidad de las conejas no fue afectada por la duración de la lactación. El destete tardío dio lugar a un mayor tamaño y peso de la camada al final del periodo de crecimiento (8,9 y 11,3 %, respectivamente) y a un menor índice de conversión por jaula durante el todo el periodo experimental (13,5 %) con respecto al destete convencional (P<0,001). Éstos resultados fueron paralelos a la menor mortalidad global (12,6 vs 17,6 %; P<0,05) observada en gazapos con destete tardío. Las diferencias en los parámetros productivos con las distintas edades al destete sólo fueron observadas en los ciclos con peor estado sanitario (tercer y quinto ciclo), en los cuales el destete tardío redujo la mortalidad. El tipo de jaula no afectó al peso de la coneja, condición corporal, mortalidad, fertilidad ni tamaño de camada durante los cinco primeros ciclos reproductivos. Sin embargo, el peso de la camada y el índice de conversión a los 21 días de edad fueron 4,2% mayor (P<0,001) y 5,0% menor (P<0,005) en animales alojados en jaulas alternativas que en jaulas convencionales. A día 59 las jaulas alternativas dieron lugar a camadas más pesadas (P<0,01); sin embargo, éste efecto fue influenciado por la densidad alcanzada en cada ciclo, ya que cuando la densidad de los animales fue menor que 40kg/m2 (tercer y quinto ciclo), el efecto del tipo de jaula sobre el peso de la camada no fue significativo. De los resultados obtenidos se puede concluir que el uso combinado de lactaciones más largas y jaulas con mayor superficie disponible con una plataforma elevada podría ser una alternativa para mejorar el bienestar animal en determinadas situaciones productivas. ABSTRACT The general aim of this PhD Thesis was to evaluate new housing and husbandry systems of farmed rabbits, studying behavioral (experiment 1), productive and reproductive (experiment 3) parameters. Moreover, different sampling techniques were evaluated in order to optimize the assessment of rabbit behaviour (experiment 2). In experiment 1, the behaviour of rabbit does housed in two different types of cage (TC), conventional vs. alternative with an elevated platform, at different physiological stages (PS), lactation and gestation was to study. Behavioural observations were carried out on 12 commercial rabbit does using continuous 24 hour video recording. Independently of PS and TC, rabbit does spent most of their time on foot mats (57.7 %, as average). However, due to the use of platforms (23.0% of time, as average), lactating does spent 36.6% less time (P<0.001) on foot mats and gestating does spent 27.0% less (P<0.001) time on wire mesh in alternative cages than in conventional cages. Alternative cages allowed for standing posture but this behaviour was only observed in gestating does (4.6 times a day, as average). Frequency of drinking was higher in conventional than in alternative cages (24.6 vs. 19.1 times a day; P<0.05). Gestating does housed in conventional cages reached the highest duration and frequency of interacting with neighbours (276 s/d and 4.6 times/d; P<0.05). The frequency of interacting with kits was lower in alternative than in conventional cages (2.4 vs. 8.6 times a day; P<0.01). Does’ behaviour was influenced by hour of day, being less active at the midday hours. During dark hours rabbit does more frequently performed restless behaviour such as hyperactivity or nursing, matching the time at which rabbit does spent more time on the platform. The platform was frequently used by rabbit does, independent of their physiological state, and during late lactation phase, when mothers were not receptive to nursing, does housed in alternative cages used the platform as a mean to flee from kids trying to suckle. The use of the platform might lead to hygienic problems due to retained faeces on the platform and faeces and urine falling onto animals located in the lower part of the cage. Stereotypies were not observed in any housing system, therefore conventional cages do not suggest lack of animal welfare. In experiment 2, it was compared the results of different simplified sampling methods of behavioural data with respect to reference records of 24-h in order to assess rabbit does behaviours at different physiological stages (gestation and lactation) in animals housed in two types of cages (conventional and alternative). A total of 576 h of continuous video of 12 rabbit does at the end of lactation and on the same females after weaning were analysed. The behavioural observations were studied using three independent categories of classification (location in the cage, posture and functional behaviours). Continuous behavioural recordings of 24 h were considered as the reference method to validate another 4 sampling methods of data collection by aggregated video recordings of different frequency and duration (regular short and long methods with 2.4 and 8 h of observation respectively, and irregular short and long methods with 6 and 8 h of observation, respectively). The current results showed that, independently of housing system, the best method to reduce the total observation time required to assess rabbit does behaviour depends on the trait studied and physiological stage of does. In gestating does, irregular methods were not suitable to estimate behaviours of long duration such as lying, sitting, resting and grooming. However, in both physiological stages, regular methods were accurate for location behaviours, postures and functional behaviours of long duration. Instead, for the study of infrequent behaviours performed mainly during dark period, where coefficients of variation were high, the irregular long method led to the lowest mean estimation errors. In experiment 3, the effects of the combined use of long lactation periods (46 days) with alternative cages on the reproductive and growth performance of 104 rabbit does and their litters during five consecutive reproductive cycles were studied. Half of does were housed in conventional polyvalent cages (39 cm x 100 cm x 30 cm) and the other half in alternative polyvalent cages (39 cm x 100 cm x 60 cm), with a raised platform. Half of the rabbit does in each type of cages were weaned at 32 and the other half at 46 days after parturition. Longer lactations affected negatively to body weight (P<0.001), fat and energy content (P<0.05) of rabbit does at the end of the lactation period, but this effect decreased with the number of parturitions. Fertility, prolificacy and doe mortality were not affected by lactation length. Late weaning led to higher litter size (by 8.9 %) and litter weight (by 11.3 %) at the end of growing period and lower feed conversion ratio per cage during the overall experimental period (13.5 %) than standard weaning (P<0.001). These results were parallels to a lower mortality (12.6 vs 17.6 %; P<0.05) of young rabbit weaned later during the overall experimental period. Differences in performances at different weaning ages were only observed during cycles with worst health status (third and fifth cycles) in which late weaning decreased mortality. Type of cage did not affect doe body weight and body condition, mortality, fertility, prolificacy and litter size during the five firsts reproductive cycles. Nevertheless, at day 21 litter weight and feed conversion ratio were 4.2 % higher (P<0.001) and 5.0 % lower (P<0.005) in animals housed in alternative than in conventional cages. Alternative cages also led to heavier litters at 59 days (P<0.01); however, this effect was influenced by density reached in each cycle, as when the density of animals was lower than 40 kg/m2 (cycles three and five), the difference of litter weight between alternative and conventional cages was not significant. From the results obtained it can be concluded that the combined use of longer lactations and cages with higher available surface with a raised platform could be an alternative to improve animal welfare in some productive situations.

Tamaño y densidad urbana: análisis de la ocupación de suelo por las áreas urbanas españolas

Tamaño y densidad son atributos imprescindibles para la caracterización de cualquier ciudad, sin embargo no siempre se les presta la misma atención. Mientras la densidad ocupa un lugar destacado en cualquier manual de planificación urbana, como variable crítica para modular la intervención urbanística, el tamaño de la ciudad tiende a tomarse como un imponderable más allá del control de la planificación. Sin embargo, no está claro que ambas variables se puedan analizar de forma autónoma, ya que las interrelaciones entre ambas son significativas. Con el fin de evaluar cómo se influyen mutuamente tamaño y densidad se plantea una metodología novedosa adaptada a las nuevas formas contemporáneas de urbanización que combina una delimitación funcional de las áreas urbanas con un análisis de su morfología en términos de densidad y compacidad. Los resultados obtenidos a partir de una muestra de 47 áreas urbanas españolas permiten clarificar algunos aspectos de la relación entre tamaño y densidad. Las grandes metrópolis españolas presentan valores relativamente altos de densidad y compacidad, pero los valores más elevados de ambos valores están presentes en áreas urbanas de tamaño intermedio o pequeño. En realidad se puede apreciar una enorme diversidad entre las áreas urbanas intermedias y pequeñas que desaparece casi por completo en las más grandes. Otro aspecto reseñable es la constatación de que densidad y compacidad son variables absolutamente independientes: existen áreas urbanas densas y compactas, pero también densas y poco compactas, o compactas pero poco densas. Sobre esta base, se propone una caracterización de la tipología tradicional de las ciudades españolas. Las áreas urbanas del Cantábrico presentan una alta densidad pero baja compacidad, en gran medida por la difícil orografía; por el contrario la áreas urbanas del interior presentan baja densidad y valores muy diversos de compacidad. Finalmente las áreas urbanas mediterráneas se caracterizan más por su alta compacidad que por su densidad, que es muy variable. En general los resultados obtenidos reflejan adecuadamente aspectos conocidos de la diversidad urbana española, pero también otros menos evidentes. Las conclusiones del estudio están muy mediatizadas por el contexto español, pero la metodología empleada ha demostrado ser de gran utilidad para el análisis de los territorios urbanos contemporáneos. ABSTRACT Size and density are essential attributes for characterizing cities, however they don't always deserve the same attention. While density is considered a critical variable in every urban planning guidebook, city size tends to be taken as an imponderable beyond the control of planning. However, it is unclear whether both variables can be analyzed independently, because their reciprocal influences are evident. To assess how city size and urban density influence each other, we propose a new methodology adapted to contemporary forms of urbanization, bringing together urban areas functional delimitation and morphological analysis in terms of density and compactness. Results obtained from a sample of 47 Spanish urban areas clarify some aspects of the size-density link. The largest Spanish urban areas are relatively dense and compact, but the highest values of both density and compactness are found in intermediate or small urban areas. Actually there is a huge diversity among intermediate and small urban areas which almost disappears in larger ones. Another noteworthy aspect is the realization that density and compactness are absolutely independent variables: there are dense and compact urban areas, but also low-density compact, and low-compactness dense areas. On this basis, we make a revision of the tradicional classification of Spanish urban areas. Cantabrian urban areas display high density and low compactness, largely because of the difficult terrain; on the contrary, inland urban areas exhibit low density and very different values of compactness. Finally, Mediterranean urban areas are characterized more by its high compactness than by its density, which is remarkably variable. In general results accurately reflect known aspects of Spanish urban diversity, but they also expose less obvious attributes. These findings may be closedly related to Spanish context, but the methodology has proven a useful tool for analyzing contemporary urban areas.

Evaluación de los efectos derivados de usos recreativos sobre las comunidades vegetales próximas a cursos de agua en el P.N. de Tijuca (Rio de Janeiro, Brasil)

La particular ubicación del Parque Nacional da Tijuca en el entorno urbano de la ciudad de Rio de Janeiro justifica la estrecha relación que sus habitantes han tenido históricamente con las montañas, sobre las que se asienta el P.N. Su población de más de 6 millones de habitantes, las altas temperaturas de su clima tropical húmedo, y el creciente interés por el turismo ecológico, hacen del baño en las cascadas de la Floresta da Tijuca una actividad habitual entre los cariocas. El posible impacto de los usos recreativos sobre las comunidades biológicas del entorno de las cascadas, en ocasiones localizadas en zonas de alto valor ecológico, aún no es suficientemente conocido, no ha sido evaluado, ni está contemplado en los planes oficiales del P.N. La decisión de permitir oficialmente el baño en una sucesión de cascadas concretas, situadas en un sector del Parque donde el baño no está permitido por el Plan de Gestión, pero donde se ha comprobado que este uso se produce con regularidad, y el potencial incremento en el número de visitantes que puede derivar de la divulgación al público de este permiso justifican el presente estudio. El trabajo, realizado bajo la co-dirección del profesor responsable del Laboratório de Ecologia de Florestal de la Universidad Federal del Estado de Rio de Janeiro, responde a la solicitud del P.N. de un estudio que evalúe el estado de conservación de la vegetación en el entorno de dichas cascadas, en el momento previo a la divulgación del permiso. Se pretende comprobar si existe degradación en el estrato arbóreo-arbustivo que pueda estar asociada a los usos recreativos que se dan en estas cascadas, bajo los niveles de visitación actuales. Se ha realizado una comparación del estado de conservación de dicho estrato, entre el entorno de las cascadas -representado por el muestreo de 10 parcelas que componen el conjunto CP-, y una zona alejada que se considera “libre” de presión antrópica debido a su menor accesibilidad -representada por el muestreo de 10 parcelas que componen el conjunto CL-. Dichos conjuntos se han comparado en lo referente a su composición florística (diversidad, en términos de riqueza observada, riqueza esperada, y equidad; especies predominantes y especies exóticas y amenazadas), atributos de desarrollo morfo-estructural (altura media, área basimétrica, proporción de individuos multicaules, densidad de individuos y proporción de individuos pertenecientes a especies exóticas y amenazadas) y atributos del medio físico que diferencian, fundamentalmente, las dos zonas comparadas (pendiente y exposición del suelo). Se han realizado comparaciones directas, en el caso de la composición florística, y comparaciones estadísticas básicas, de análisis de varianzas (ANOVA) y de análisis multivariante gráficos (de agrupamiento y NMDS) y numéricos (ANOSIM, NPMANOVA) para identificar diferencias significativas entre el entorno de las cascadas y la zona considerada “libre de impacto”. Los resultados muestran que en el momento inicial no es posible afirmar que exista degradación patente asociada a los usos recreativos en las cascadas, al no haber sido encontradas diferencias significativas en la mayor parte de las comparaciones, ni en la evaluación conjunta de todos los indicadores. Sin embargo, se han detectado importantes indicios de simplificación de la vegetación en el entorno del río, tanto en los indicadores de desarrollo estructural, como en la composición florística; y se han detectado diferencias relevantes en determinados indicadores que señalan la existencia de perturbación derivada del tránsito de visitantes, en especial sobre la exposición del suelo. Estos resultados podrían deberse a un estado incipiente de degradación, aún no claramente perceptible, y recalcan la importancia de estudiar, en los años sucesivos, la evolución del estado de conservación de la vegetación de esta zona, correlativamente a la evolución del volumen de visitantes, para comprobar y cuantificar el posible impacto antrópico derivado de los usos recreativos en estas cascadas.