Measurements and modelling of radio propagation in subway tunnels

The work presented in this article is focused on the RF measurement campaign carried out in several subway tunnels in Metro Madrid (Spain). Most common segments such as straight lines, curves and passing through station as well as other unique scenarios in metropolitan lines were the selected locations during this campaign. Measurements were conducted in tunnels of diverse cross section shapes and taken at three frequency bands: 900 MHz, 2.4GHz and 5.7 GHz for both horizontal and vertical polarization.

Propagation losses by tunnelling in barrier optical waveguides

An analytical model is proposed in order to estimate the optical propagation losses due to tunnelling in barrier waveguides. The results are validated by means of a beam propagation method (BPM) simulations for different waveguides conditions.

Propagation of nuclear data uncertainties for PWR core analysis

An uncertainty propagation methodology based on the Monte Carlo method is applied to PWR nuclear design analysis to assess the impact of nuclear data uncertainties. The importance of the nuclear data uncertainties for 235,238 U, 239 Pu, and the thermal scattering library for hydrogen in water is analyzed. This uncertainty analysis is compared with the design and acceptance criteria to assure the adequacy of bounding estimates in safety margins.

Experimental investigation of radio signal propagation in scientific facilities for telerobotic applications

Understanding the radio signal transmission characteristics in the environment where the telerobotic application is sought is a key part of achieving a reliable wireless communication link between a telerobot and a control station. In this paper, wireless communication requirements and a case study of a typical telerobotic application in an underground facility at CERN are presented. Then, the theoretical and experimental characteristics of radio propagation are investigated with respect to time, distance, location and surrounding objects. Based on analysis of the experimental findings, we show how a commercial wireless system, such as Wi-Fi, can be made suitable for a case study application at CERN.

Development of a methodology for nuclear data uncertainty propagation on isotopic evolution calculations for advanced nuclear systems

Una apropiada evaluación de los márgenes de seguridad de una instalación nuclear, por ejemplo, una central nuclear, tiene en cuenta todas las incertidumbres que afectan a los cálculos de diseño, funcionanmiento y respuesta ante accidentes de dicha instalación. Una fuente de incertidumbre son los datos nucleares, que afectan a los cálculos neutrónicos, de quemado de combustible o activación de materiales. Estos cálculos permiten la evaluación de las funciones respuesta esenciales para el funcionamiento correcto durante operación, y también durante accidente. Ejemplos de esas respuestas son el factor de multiplicación neutrónica o el calor residual después del disparo del reactor. Por tanto, es necesario evaluar el impacto de dichas incertidumbres en estos cálculos. Para poder realizar los cálculos de propagación de incertidumbres, es necesario implementar metodologías que sean capaces de evaluar el impacto de las incertidumbres de estos datos nucleares. Pero también es necesario conocer los datos de incertidumbres disponibles para ser capaces de manejarlos. Actualmente, se están invirtiendo grandes esfuerzos en mejorar la capacidad de analizar, manejar y producir datos de incertidumbres, en especial para isótopos importantes en reactores avanzados. A su vez, nuevos programas/códigos están siendo desarrollados e implementados para poder usar dichos datos y analizar su impacto. Todos estos puntos son parte de los objetivos del proyecto europeo ANDES, el cual ha dado el marco de trabajo para el desarrollo de esta tesis doctoral. Por tanto, primero se ha llevado a cabo una revisión del estado del arte de los datos nucleares y sus incertidumbres, centrándose en los tres tipos de datos: de decaimiento, de rendimientos de fisión y de secciones eficaces. A su vez, se ha realizado una revisión del estado del arte de las metodologías para la propagación de incertidumbre de estos datos nucleares. Dentro del Departamento de Ingeniería Nuclear (DIN) se propuso una metodología para la propagación de incertidumbres en cálculos de evolución isotópica, el Método Híbrido. Esta metodología se ha tomado como punto de partida para esta tesis, implementando y desarrollando dicha metodología, así como extendiendo sus capacidades. Se han analizado sus ventajas, inconvenientes y limitaciones. El Método Híbrido se utiliza en conjunto con el código de evolución isotópica ACAB, y se basa en el muestreo por Monte Carlo de los datos nucleares con incertidumbre. En esta metodología, se presentan diferentes aproximaciones según la estructura de grupos de energía de las secciones eficaces: en un grupo, en un grupo con muestreo correlacionado y en multigrupos. Se han desarrollado diferentes secuencias para usar distintas librerías de datos nucleares almacenadas en diferentes formatos: ENDF-6 (para las librerías evaluadas), COVERX (para las librerías en multigrupos de SCALE) y EAF (para las librerías de activación). Gracias a la revisión del estado del arte de los datos nucleares de los rendimientos de fisión se ha identificado la falta de una información sobre sus incertidumbres, en concreto, de matrices de covarianza completas. Además, visto el renovado interés por parte de la comunidad internacional, a través del grupo de trabajo internacional de cooperación para evaluación de datos nucleares (WPEC) dedicado a la evaluación de las necesidades de mejora de datos nucleares mediante el subgrupo 37 (SG37), se ha llevado a cabo una revisión de las metodologías para generar datos de covarianza. Se ha seleccionando la actualización Bayesiana/GLS para su implementación, y de esta forma, dar una respuesta a dicha falta de matrices completas para rendimientos de fisión. Una vez que el Método Híbrido ha sido implementado, desarrollado y extendido, junto con la capacidad de generar matrices de covarianza completas para los rendimientos de fisión, se han estudiado diferentes aplicaciones nucleares. Primero, se estudia el calor residual tras un pulso de fisión, debido a su importancia para cualquier evento después de la parada/disparo del reactor. Además, se trata de un ejercicio claro para ver la importancia de las incertidumbres de datos de decaimiento y de rendimientos de fisión junto con las nuevas matrices completas de covarianza. Se han estudiado dos ciclos de combustible de reactores avanzados: el de la instalación europea para transmutación industrial (EFIT) y el del reactor rápido de sodio europeo (ESFR), en los cuales se han analizado el impacto de las incertidumbres de los datos nucleares en la composición isotópica, calor residual y radiotoxicidad. Se han utilizado diferentes librerías de datos nucleares en los estudios antreriores, comparando de esta forma el impacto de sus incertidumbres. A su vez, mediante dichos estudios, se han comparando las distintas aproximaciones del Método Híbrido y otras metodologías para la porpagación de incertidumbres de datos nucleares: Total Monte Carlo (TMC), desarrollada en NRG por A.J. Koning y D. Rochman, y NUDUNA, desarrollada en AREVA GmbH por O. Buss y A. Hoefer. Estas comparaciones demostrarán las ventajas del Método Híbrido, además de revelar sus limitaciones y su rango de aplicación. ABSTRACT For an adequate assessment of safety margins of nuclear facilities, e.g. nuclear power plants, it is necessary to consider all possible uncertainties that affect their design, performance and possible accidents. Nuclear data are a source of uncertainty that are involved in neutronics, fuel depletion and activation calculations. These calculations can predict critical response functions during operation and in the event of accident, such as decay heat and neutron multiplication factor. Thus, the impact of nuclear data uncertainties on these response functions needs to be addressed for a proper evaluation of the safety margins. Methodologies for performing uncertainty propagation calculations need to be implemented in order to analyse the impact of nuclear data uncertainties. Nevertheless, it is necessary to understand the current status of nuclear data and their uncertainties, in order to be able to handle this type of data. Great eórts are underway to enhance the European capability to analyse/process/produce covariance data, especially for isotopes which are of importance for advanced reactors. At the same time, new methodologies/codes are being developed and implemented for using and evaluating the impact of uncertainty data. These were the objectives of the European ANDES (Accurate Nuclear Data for nuclear Energy Sustainability) project, which provided a framework for the development of this PhD Thesis. Accordingly, first a review of the state-of-the-art of nuclear data and their uncertainties is conducted, focusing on the three kinds of data: decay, fission yields and cross sections. A review of the current methodologies for propagating nuclear data uncertainties is also performed. The Nuclear Engineering Department of UPM has proposed a methodology for propagating uncertainties in depletion calculations, the Hybrid Method, which has been taken as the starting point of this thesis. This methodology has been implemented, developed and extended, and its advantages, drawbacks and limitations have been analysed. It is used in conjunction with the ACAB depletion code, and is based on Monte Carlo sampling of variables with uncertainties. Different approaches are presented depending on cross section energy-structure: one-group, one-group with correlated sampling and multi-group. Differences and applicability criteria are presented. Sequences have been developed for using different nuclear data libraries in different storing-formats: ENDF-6 (for evaluated libraries) and COVERX (for multi-group libraries of SCALE), as well as EAF format (for activation libraries). A revision of the state-of-the-art of fission yield data shows inconsistencies in uncertainty data, specifically with regard to complete covariance matrices. Furthermore, the international community has expressed a renewed interest in the issue through the Working Party on International Nuclear Data Evaluation Co-operation (WPEC) with the Subgroup (SG37), which is dedicated to assessing the need to have complete nuclear data. This gives rise to this review of the state-of-the-art of methodologies for generating covariance data for fission yields. Bayesian/generalised least square (GLS) updating sequence has been selected and implemented to answer to this need. Once the Hybrid Method has been implemented, developed and extended, along with fission yield covariance generation capability, different applications are studied. The Fission Pulse Decay Heat problem is tackled first because of its importance during events after shutdown and because it is a clean exercise for showing the impact and importance of decay and fission yield data uncertainties in conjunction with the new covariance data. Two fuel cycles of advanced reactors are studied: the European Facility for Industrial Transmutation (EFIT) and the European Sodium Fast Reactor (ESFR), and response function uncertainties such as isotopic composition, decay heat and radiotoxicity are addressed. Different nuclear data libraries are used and compared. These applications serve as frameworks for comparing the different approaches of the Hybrid Method, and also for comparing with other methodologies: Total Monte Carlo (TMC), developed at NRG by A.J. Koning and D. Rochman, and NUDUNA, developed at AREVA GmbH by O. Buss and A. Hoefer. These comparisons reveal the advantages, limitations and the range of application of the Hybrid Method.

Modelling of landslides propagation with SPH : effects of rheology and pore water pressure

El estudio desarrollado en este trabajo de tesis se centra en la modelización numérica de la fase de propagación de los deslizamientos rápidos de ladera a través del método sin malla Smoothed Particle Hydrodynamics (SPH). Este método tiene la gran ventaja de permitir el análisis de problemas de grandes deformaciones evitando operaciones costosas de remallado como en el caso de métodos numéricos con mallas tal como el método de los Elementos Finitos. En esta tesis, particular atención viene dada al rol que la reología y la presión de poros desempeñan durante estos eventos. El modelo matemático utilizado se basa en la formulación de Biot-Zienkiewicz v - pw, que representa el comportamiento, expresado en términos de velocidad del esqueleto sólido y presiones de poros, de la mezcla de partículas sólidas en un medio saturado. Las ecuaciones que gobiernan el problema son: • la ecuación de balance de masa de la fase del fluido intersticial, • la ecuación de balance de momento de la fase del fluido intersticial y de la mezcla, • la ecuación constitutiva y • una ecuación cinemática. Debido a sus propiedades geométricas, los deslizamientos de ladera se caracterizan por tener una profundidad muy pequeña frente a su longitud y a su anchura, y, consecuentemente, el modelo matemático mencionado anteriormente se puede simplificar integrando en profundidad las ecuaciones, pasando de un modelo 3D a 2D, el cual presenta una combinación excelente de precisión, sencillez y costes computacionales. El modelo propuesto en este trabajo se diferencia de los modelos integrados en profundidad existentes por incorporar un ulterior modelo capaz de proveer información sobre la presión del fluido intersticial a cada paso computacional de la propagación del deslizamiento. En una manera muy eficaz, la evolución de los perfiles de la presión de poros está numéricamente resuelta a través de un esquema explicito de Diferencias Finitas a cada nodo SPH. Este nuevo enfoque es capaz de tener en cuenta la variación de presión de poros debida a cambios de altura, de consolidación vertical o de cambios en las tensiones totales. Con respecto al comportamiento constitutivo, uno de los problemas principales al modelizar numéricamente deslizamientos rápidos de ladera está en la dificultad de simular con la misma ley constitutiva o reológica la transición de la fase de iniciación, donde el material se comporta como un sólido, a la fase de propagación donde el material se comporta como un fluido. En este trabajo de tesis, se propone un nuevo modelo reológico basado en el modelo viscoplástico de Perzyna, pensando a la viscoplasticidad como a la llave para poder simular tanto la fase de iniciación como la de propagación con el mismo modelo constitutivo. Con el fin de validar el modelo matemático y numérico se reproducen tanto ejemplos de referencia con solución analítica como experimentos de laboratorio. Finalmente, el modelo se aplica a casos reales, con especial atención al caso del deslizamiento de 1966 en Aberfan, mostrando como los resultados obtenidos simulan con éxito estos tipos de riesgos naturales. The study developed in this thesis focuses on the modelling of landslides propagation with the Smoothed Particle Hydrodynamics (SPH) meshless method which has the great advantage of allowing to deal with large deformation problems by avoiding expensive remeshing operations as happens for mesh methods such as, for example, the Finite Element Method. In this thesis, special attention is given to the role played by rheology and pore water pressure during these natural hazards. The mathematical framework used is based on the v - pw Biot-Zienkiewicz formulation, which represents the behaviour, formulated in terms of soil skeleton velocity and pore water pressure, of the mixture of solid particles and pore water in a saturated media. The governing equations are: • the mass balance equation for the pore water phase, • the momentum balance equation for the pore water phase and the mixture, • the constitutive equation and • a kinematic equation. Landslides, due to their shape and geometrical properties, have small depths in comparison with their length or width, therefore, the mathematical model aforementioned can then be simplified by depth integrating the equations, switching from a 3D to a 2D model, which presents an excellent combination of accuracy, computational costs and simplicity. The proposed model differs from previous depth integrated models by including a sub-model able to provide information on pore water pressure profiles at each computational step of the landslide's propagation. In an effective way, the evolution of the pore water pressure profiles is numerically solved through a set of 1D Finite Differences explicit scheme at each SPH node. This new approach is able to take into account the variation of the pore water pressure due to changes of height, vertical consolidation or changes of total stress. Concerning the constitutive behaviour, one of the main issues when modelling fast landslides is the difficulty to simulate with the same constitutive or rheological model the transition from the triggering phase, where the landslide behaves like a solid, to the propagation phase, where the landslide behaves in a fluid-like manner. In this work thesis, a new rheological model is proposed, based on the Perzyna viscoplastic model, thinking of viscoplasticity as the key to close the gap between the triggering and the propagation phase. In order to validate the mathematical model and the numerical approach, benchmarks and laboratory experiments are reproduced and compared to analytical solutions when possible. Finally, applications to real cases are studied, with particular attention paid to the Aberfan flowslide of 1966, showing how the mathematical model accurately and successfully simulate these kind of natural hazards.

Modelling of short runout propagation landslides and debris flows

This paper proposes an extension of methods used to predict the propagation of landslides having a long runout to smaller landslides with much shorter propagation distances. The method is based on: (1) a depth-integrated mathematical model including the coupling between the soil skeleton and the pore fluids, (2) suitable rheological models describing the relation between the stress and the rate of deformation tensors for fluidised soils and (3) a meshless numerical method, Smooth Particle Hydrodynamics, which separates the computational mesh (or set of computational nodes) from the mesh describing the terrain topography, which is of structured type – thus accelerating search operations. The proposed model is validated using two examples for which there are analytical solutions, and then it is applied to two short runout landslides which happened in Hong Kong in 1995, for which there is available information.

Neurite, a finite difference large scale parallel program for the simulation of the electrical signal propagation in neurites under mechanical loading

With the growing body of research on traumatic brain injury and spinal cord injury, computational neuroscience has recently focused its modeling efforts on neuronal functional deficits following mechanical loading. However, in most of these efforts, cell damage is generally only characterized by purely mechanistic criteria, function of quantities such as stress, strain or their corresponding rates. The modeling of functional deficits in neurites as a consequence of macroscopic mechanical insults has been rarely explored. In particular, a quantitative mechanically based model of electrophysiological impairment in neuronal cells has only very recently been proposed (Jerusalem et al., 2013). In this paper, we present the implementation details of Neurite: the finite difference parallel program used in this reference. Following the application of a macroscopic strain at a given strain rate produced by a mechanical insult, Neurite is able to simulate the resulting neuronal electrical signal propagation, and thus the corresponding functional deficits. The simulation of the coupled mechanical and electrophysiological behaviors requires computational expensive calculations that increase in complexity as the network of the simulated cells grows. The solvers implemented in Neurite-explicit and implicit-were therefore parallelized using graphics processing units in order to reduce the burden of the simulation costs of large scale scenarios. Cable Theory and Hodgkin-Huxley models were implemented to account for the electrophysiological passive and active regions of a neurite, respectively, whereas a coupled mechanical model accounting for the neurite mechanical behavior within its surrounding medium was adopted as a link between lectrophysiology and mechanics (Jerusalem et al., 2013). This paper provides the details of the parallel implementation of Neurite, along with three different application examples: a long myelinated axon, a segmented dendritic tree, and a damaged axon. The capabilities of the program to deal with large scale scenarios, segmented neuronal structures, and functional deficits under mechanical loading are specifically highlighted.

Propagation of DEM Varying Accuracy into Terrain-Based Analysis.

Análisis de precisión en modelos digitales de elevación globales. ABSTRACT: Terrain-Based Analysis results in derived products from an input DEM and these products are needed to perform various analyses. To efficiently use these products in decision-making, their accuracies must be estimated systematically. This paper proposes a procedure to assess the accuracy of these derived products, by calculating the accuracy of the slope dataset and its significance, taking as an input the accuracy of the DEM. Based on the output of previously published research on modeling the relative accuracy of a DEM, specifically ASTER and SRTM DEMs with Lebanon coverage as the area of study, analysis have showed that ASTER has a low significance in the majority of the area where only 2% of the modeled terrain has 50% or more significance. On the other hand, SRTM showed a better significance, where 37% of the modeled terrain has 50% or more significance. Statistical analysis deduced that the accuracy of the slope dataset, calculated on a cell-by-cell basis, is highly correlated to the accuracy of the input DEM. However, this correlation becomes lower between the slope accuracy and the slope significance, whereas it becomes much higher between the modeled slope and the slope significance.

Development of somatic embryogenesis for cloning and conservation of mature holm oak trees (Quercus ilex L.)

La encina (Quercus ilex L.) es una de las especies forestales mediterráneas más importantes. Constituye gran parte del estrato arbóreo de dehesas o montados, produce bellota como alimento del ganado y establece simbiosis con hongos micorrizógenos de gran valor económico. La encina está considerada como una especie recalcitrante en términos de conservación de semillas y capacidad morfogénica, lo que dificulta los programas de conservación de recursos genéticos y la mejora de la especie. La propagación vegetativa es una potente herramienta de los programas de mejora, por lo que es preciso desarrollar protocolos de regeneración somática en encina. La embriogénesis somática está considerada como la modalidad más adecuada de regeneración basada en técnicas de cultivo de tejidos vegetales utilizada en biotecnología forestal. Este trabajo se centra en el estudio de determinados aspectos de la embriogénesis somática para la regeneración clonal de encinas adultas. La memoria de esta tesis se ha dividido en capítulos que se corresponden con diferentes aspectos del sistema embriogénico. La embriogénesis somática se indujo en tegumentos maternos de óvulos en desarrollo procedentes de bellotas inmaduras de encinas adultas. A pesar de las bajas frecuencias de inducción, las líneas embriogénicas generadas se amplificaron mediante embriogénesis secundaria observándose cierta pérdida de la capacidad de diferenciación con el tiempo. Tanto el genotipo como la formulación del medio de cultivo influyeron en la respuesta embriogénica, concluyendo que la formulación de macronutrientes de Schenk y Hildebrant del medio sin reguladores de crecimiento fue la combinación más efectiva en la inducción. Los resultados sugirieron la existencia de una ventana en el desarrollo del óvulo más sensible a la inducción. El genotipo in[luyó en la capacidad proliferativa de los cultivos y en la conversión de los embriones somáticos, que se incrementó suplementando el medio con ácido indol-3-butírico y 6-benciladenina. El cultivo en medio líquido de líneas embriogénicas en condiciones de inmersión transitoria incrementó el crecimiento, dependiendo del genotipo, con respecto al cultivo en medio semisólido. Sin embargo, no mejoró la capacidad de diferenciar embriones cotiledonares aislados. Se estableció un protocolo de inicio y mantenimiento de cultivos en suspensión para varias líneas embriogénicas mediante inoculación en alta densidad de agregados embrionarios procedentes del medio semisólido. Para evitar la pérdida de vigor y la capacidad morfogénica debida al cultivo prolongado se desarrolló un protocolo de crioconservación de líneas embriogénicas mediante vitrificación. Al determinar la influencia de los agentes crioprotectores antes y después de su inmersión en nitrógeno líquido se concluyó que las respuestas de capacidad de crecimiento y de diferenciación del material embriogénico son independientes, además de estar bajo influencia del genotipo y el tipo de material crioconservado. La combinación de sacarosa y PVS2 previa a la inmersión en nitrógeno líquido proporcionó la mayor tasa de recuperación. Cuando las líneas fueron crioconservadas 30 días la capacidad de diferenciación se perdió en todas ellas. El análisis de SSR detectó variación somaclonal en el material crioconservado a corto plazo. SSR y RAPD mostraron importantes diferencias genéticas entre los árboles donantes y el material embriogénico que dependieron del genotipo. El grado de detección dependió del marcador empleado. Ambos marcadores revelaron baja inestabilidad intraclonal. Los RAPD revelaron variación genética intra-individuo en las encinas donantes. Se discuten la variación genética pre-existente en encina, su aparición durante las primeras fases de la inducción de embriogénesis, y la presencia de tejidos provenientes de la fertilización en el explanto materno. Esto hace preciso definir la identidad genética del material donante y acometer ensayos de detección precoz de variación somaclonal. ABSTRACT Holm oak (Quercus ilex L.) is one of the most important Mediterranean forest species. It conforms the tree layer of dehesas or montados, it produces acorns to feed the livestock and it establishes symbiosis with profitable mycorrhizal fungi. Holm oak is considered as recalcitrant species in terms of seed conservation and morphogenic capacities, which complicates the development of genetic conservation and improvement programs. Vegetative propagation is one of the mightiest tools for breeding programs therefore; developing protocols for clonal regeneration of holm oak is essential. Somatic embryogenesis is considered the best tissue culture-based way of plant regeneration in forest biotechnology. The present study is focused on the study of certain aspects of somatic embryogenesis for clonal regeneration of mature holm oak. This thesis manuscript is divided into several chapters that match with different aspects of the embryogenic system. Somatic embryogenesis induction was achieved on maternal teguments of developing ovules from immature acorns of adult holm oak trees. Despite the low induction frequencies, the generated embryogenic lines were amplified by secondary embryogenesis. A decline in the differentiation capacity over time was also observed. It was concluded that both genotype and culture media formulation influenced the embryogenic response, being the Schenk and Hildebrandt´s macronutrients formulation from culture medium and the lack of plant growth regulators the most effective combination for the induction of the embryogenic response. It has been suggested the existence of a developmental window in which ovules are prone to induction. Genotype influenced the proliferation capacity and the plant conversion of somatic embryos, which was also favoured by the presence of indol-3-butyric acid and 6-bencyladenine. The use of temporary immersion systems as proliferation in liquid culture of the embryogenic lines increased the growth depending on genotype, when compared to semisolid cultures. However, it did not improve the differentiation of single cotyledonary embryos. A protocol for the initiation and maintenance of embryogenic suspension cultures was established for several embryogenic lines with highly dense inoculi of embryogenic clusters from proliferating semisolid cultures. In order to avoid the loss of vigour and morphogenic ability of embryogenic lines due to prolonged cultures, a cryopreservation protocol for embryogenic lines of holm oak has been developed. During the determination of the influence of cryoprotective agents on the growth and differentiation capacities before and after liquid nitrogen immersion, it was concluded that both responses were independent from each other and also under the influence of genotype and the type of cryopreserved material. The combination of sucrose and PVS2 prior liquid nitrogen immersion provided higher recovery rates. When the same embryogenic lines were cryopreserved for 30 days, none was able to differentiate. The SSRs analysis of the short-term cryopreserved material detected somaclonal variation. Both SSR and RAPD markers showed high sensitivity to detect genetic differences between the donor trees and the generated embryogenic material. Nevertheless, the degree of instability detection depended on the marker. The SSR analysis indicated a relationship between genotype, the studied loci and the located polymorphisms. Also, both markers revealed low intraclonal genetic variation. The RAPD detected genetic variation within the donor trees. The presence of pre-existent genetic variation within mature trees, in addition to its occurrence during the early stages of the embryogenic induction, and the presence of tissues of fertilisation origin within the maternal explants are all discussed. Nonetheless, the determination of the genetic identity of donor material is required, in addition to early detection methods of somaclonal variation.

Development of somatic embryogenesis for cloning and conservation of mature holm oak trees (Quercus ilex L.)

La encina (Quercus ilex L.) es una de las especies forestales mediterráneas más importantes. Constituye gran parte del estrato arbóreo de dehesas o montados, produce bellota como alimento del ganado y establece simbiosis con hongos micorrizógenos de gran valor económico. La encina está considerada como una especie recalcitrante en términos de conservación de semillas y capacidad morfogénica, lo que dificulta los programas de conservación de recursos genéticos y la mejora de la especie. La propagación vegetativa es una potente herramienta de los programas de mejora, por lo que es preciso desarrollar protocolos de regeneración somática en encina. La embriogénesis somática está considerada como la modalidad más adecuada de regeneración basada en técnicas de cultivo de tejidos vegetales utilizada en biotecnología forestal. Este trabajo se centra en el estudio de determinados aspectos de la embriogénesis somática para la regeneración clonal de encinas adultas. La memoria de esta tesis se ha dividido en capítulos que se corresponden con diferentes aspectos del sistema embriogénico. La embriogénesis somática se indujo en tegumentos maternos de óvulos en desarrollo procedentes de bellotas inmaduras de encinas adultas. A pesar de las bajas frecuencias de inducción, las líneas embriogénicas generadas se amplificaron mediante embriogénesis secundaria observándose cierta pérdida de la capacidad de diferenciación con el tiempo. Tanto el genotipo como la formulación del medio de cultivo influyeron en la respuesta embriogénica, concluyendo que la formulación de macronutrientes de Schenk y Hildebrant del medio sin reguladores de crecimiento fue la combinación más efectiva en la inducción. Los resultados sugirieron la existencia de una ventana en el desarrollo del óvulo más sensible a la inducción. El genotipo in[luyó en la capacidad proliferativa de los cultivos y en la conversión de los embriones somáticos, que se incrementó suplementando el medio con ácido indol-3-butírico y 6-benciladenina. El cultivo en medio líquido de líneas embriogénicas en condiciones de inmersión transitoria incrementó el crecimiento, dependiendo del genotipo, con respecto al cultivo en medio semisólido. Sin embargo, no mejoró la capacidad de diferenciar embriones cotiledonares aislados. Se estableció un protocolo de inicio y mantenimiento de cultivos en suspensión para varias líneas embriogénicas mediante inoculación en alta densidad de agregados embrionarios procedentes del medio semisólido. Para evitar la pérdida de vigor y la capacidad morfogénica debida al cultivo prolongado se desarrolló un protocolo de crioconservación de líneas embriogénicas mediante vitrificación. Al determinar la influencia de los agentes crioprotectores antes y después de su inmersión en nitrógeno líquido se concluyó que las respuestas de capacidad de crecimiento y de diferenciación del material embriogénico son independientes, además de estar bajo influencia del genotipo y el tipo de material crioconservado. La combinación de sacarosa y PVS2 previa a la inmersión en nitrógeno líquido proporcionó la mayor tasa de recuperación. Cuando las líneas fueron crioconservadas 30 días la capacidad de diferenciación se perdió en todas ellas. El análisis de SSR detectó variación somaclonal en el material crioconservado a corto plazo. SSR y RAPD mostraron importantes diferencias genéticas entre los árboles donantes y el material embriogénico que dependieron del genotipo. El grado de detección dependió del marcador empleado. Ambos marcadores revelaron baja inestabilidad intraclonal. Los RAPD revelaron variación genética intra-individuo en las encinas donantes. Se discuten la variación genética pre-existente en encina, su aparición durante las primeras fases de la inducción de embriogénesis, y la presencia de tejidos provenientes de la fertilización en el explanto materno. Esto hace preciso definir la identidad genética del material donante y acometer ensayos de detección precoz de variación somaclonal. ABSTRACT Holm oak (Quercus ilex L.) is one of the most important Mediterranean forest species. It conforms the tree layer of dehesas or montados, it produces acorns to feed the livestock and it establishes symbiosis with profitable mycorrhizal fungi. Holm oak is considered as recalcitrant species in terms of seed conservation and morphogenic capacities, which complicates the development of genetic conservation and improvement programs. Vegetative propagation is one of the mightiest tools for breeding programs therefore; developing protocols for clonal regeneration of holm oak is essential. Somatic embryogenesis is considered the best tissue culture-based way of plant regeneration in forest biotechnology. The present study is focused on the study of certain aspects of somatic embryogenesis for clonal regeneration of mature holm oak. This thesis manuscript is divided into several chapters that match with different aspects of the embryogenic system. Somatic embryogenesis induction was achieved on maternal teguments of developing ovules from immature acorns of adult holm oak trees. Despite the low induction frequencies, the generated embryogenic lines were amplified by secondary embryogenesis. A decline in the differentiation capacity over time was also observed. It was concluded that both genotype and culture media formulation influenced the embryogenic response, being the Schenk and Hildebrandt´s macronutrients formulation from culture medium and the lack of plant growth regulators the most effective combination for the induction of the embryogenic response. It has been suggested the existence of a developmental window in which ovules are prone to induction. Genotype influenced the proliferation capacity and the plant conversion of somatic embryos, which was also favoured by the presence of indol-3-butyric acid and 6-bencyladenine. The use of temporary immersion systems as proliferation in liquid culture of the embryogenic lines increased the growth depending on genotype, when compared to semisolid cultures. However, it did not improve the differentiation of single cotyledonary embryos. A protocol for the initiation and maintenance of embryogenic suspension cultures was established for several embryogenic lines with highly dense inoculi of embryogenic clusters from proliferating semisolid cultures. In order to avoid the loss of vigour and morphogenic ability of embryogenic lines due to prolonged cultures, a cryopreservation protocol for embryogenic lines of holm oak has been developed. During the determination of the influence of cryoprotective agents on the growth and differentiation capacities before and after liquid nitrogen immersion, it was concluded that both responses were independent from each other and also under the influence of genotype and the type of cryopreserved material. The combination of sucrose and PVS2 prior liquid nitrogen immersion provided higher recovery rates. When the same embryogenic lines were cryopreserved for 30 days, none was able to differentiate. The SSRs analysis of the short-term cryopreserved material detected somaclonal variation. Both SSR and RAPD markers showed high sensitivity to detect genetic differences between the donor trees and the generated embryogenic material. Nevertheless, the degree of instability detection depended on the marker. The SSR analysis indicated a relationship between genotype, the studied loci and the located polymorphisms. Also, both markers revealed low intraclonal genetic variation. The RAPD detected genetic variation within the donor trees. The presence of pre-existent genetic variation within mature trees, in addition to its occurrence during the early stages of the embryogenic induction, and the presence of tissues of fertilisation origin within the maternal explants are all discussed. Nonetheless, the determination of the genetic identity of donor material is required, in addition to early detection methods of somaclonal variation.

A new regularization method for fast and accurate propagation of roto-translationally coupled asteroids

There is a well-distinguished group of asteroids for which the roto-translational cou-pling is known to have a non-negligible e�ect in the long-term. The study of such asteroids suggests the use of specialized propagation techniques, where perturbation methods make their best. The techniques from which the special regularization method DROMO is derived, have now been extended to the attitude dynamics, with equally remarkable results in terms of speed and accuracy, thus making the combination of these algorithms specially. well-suited to deal with the propagation of bodies with strong attitude coupling.

Unearthing the roots of degradation of Quercus pyrenaica coppices: an integrative perspective from clonal structure to carbon budgets

Quercus pyrenaica es una especie rebrotadora de raíz intensa e históricamente aprovechada en monte bajo para la obtención de leñas, carbón y pastos. Debido al éxodo rural y a la aparición de nuevas fuentes energéticas, este aprovechamiento fue abandonado en la década de 1970. Desde entonces, las bajas producciones de madera y bellota y el puntisecado de los pies evidencian el generalizado estancamiento de estas masas. Uno de los mayores retos actuales de la selvicultura en el ámbito mediterráneo es encontrar usos alternativos para estos montes abandonados, siendo la conversión a monte alto una de las alternativas preferidas. Se han realizado resalveos de conversión, sin embrago, éstos se aplican sin un conocimiento integral de las causas de la degradación. En esta tesis doctoral, estudiamos un hipotético desequilibrio entre la parte radical y la parte aérea (R:S) de las cepas de rebollo como causa subyacente de su decaimiento. En una parcela experimental, aprovechada al menos desde el siglo XII, se realizaron análisis genéticos a priori para elucidar la estructura genética del rodal, y así estudiar la influencia del tamaño clonal en el funcionamiento de las cepas. Las cepas de mayor tamaño presentaron un menor crecimiento diametral de sus pies, así como mayores tasas de respiración radical, estimadas a partir de flujos internos de CO2 a través del xilema (FT) y de los flujos de CO2 del suelo. Estos resultados sugieren que el desequilibrio R:S aumenta con el tamaño clonal, dado que la eliminación periódica de órganos aéreos, al mismo tiempo que las raíces permanecen intactas, da lugar a un gran desarrollo del sistema radical que consume gran parte de los carbohidratos no estructurales (NSC) en respiración de mantenimiento, comprometiendo así el desarrollo de órganos aéreos. Se excavaron y pesaron dos cepas compuestas por cuatro y ocho pies, las cuales mostraron ratios R:S (0.5 y 1, respectivamente) superiores a los registrados en pies de origen sexual. Al igual que en otras especies rebrotadoras de raíz, se observaron altas concentraciones de NSC en las raíces (> 20% en primavera) y una gran proporción de albura en el sistema radical (52%) que alberga una notable reserva de NSC (87 kg en la cepa de mayor tamaño). En el sistema radical de dicha cepa, estimada mediante dataciones radiocarbónicas en 550 años de edad, se contaron 248 uniones radicales. La persistencia de sistemas radicales grandes, viejos, y altamente interconectados sugiere que la gran cantidad de recursos almacenados y consumidos en las raíces compensan un pobre desarrollo aéreo con una alta resiliencia vegetativa. Para un mejor entendimiento de los balances de carbono y del agotamiento de NSC en las cepas de rebollo, se midieron los flujos internos y externos de CO2 en troncos y los flujos de CO2 del suelo, y se estimó la respiración de órganos aéreos (RS) y subterráneos (RR). Estacionalmente, RS y RR reflejaron las dinámicas de flujo de savia y de crecimiento del tronco, y estuvieron determinadas principalmente por los flujos externos de CO2, dada la escasa contribución de FT a RS y RR (< 10% y < 2%, respectivamente). En una escala circadiana, la contribución de FT a RS aumentó hasta un 25% en momentos de alta transpiración. Las bajas concentraciones de CO2 en el xilema ([CO2] hasta un 0.11%) determinaron comparativamente unos bajos FT, probablemente causados por una limitada respiración del xilema y una baja resistencia a la difusión radial del CO2 impuestos por la sequía estival. Los pulsos de [CO2] observados tras las primeras lluvias de otoño apoyan esta idea. A lo largo del periodo vegetativo, el flujo medio de CO2 procedente del suelo (39 mol CO2 day-1) fue el mayor flujo respiratorio, tres y cuatro veces superior a RS (12 mol CO2 day-1) y RR (8-9 mol CO2 day-1), respectivamente. Ratios RR/RS menores que la unidad evidencian un importante peso de la respiración aérea como sumidero de carbono adicional. Finalmente, se ensayó el zanjado de raíces y el anillamiento de troncos como tratamientos selvícolas alternativos con el objetivo de aumentar las reservas de NSC en los troncos de las cepas. Los resultados preliminares desaconsejan el zanjado de raíces por el alto coste derivado posiblemente de la cicatrización de las heridas. El anillado de troncos imposibilitó el transporte de NSC a las raíces y aumentó la concentración de almidón por encima de la zona anillada, mientras que sistema radical se mantiene por los pies no anillados de la cepa. Son necesarias más mediciones y datos adicionales para comprobar el mantenimiento de esta respuesta positiva a largo plazo. Para concluir, destacamos la necesidad de estudios multidisciplinares que permitan una comprensión integral de la degradación de los rebollares ibéricos para poder aplicar a posteriori una gestión adecuada en estos montes bajos abandonados. ABSTRACT Quercus pyrenaica is a vigorous root-resprouting species intensively and historically coppiced for firewood, charcoal and woody pastures. Due to the rural exodus and the appearance of new energy sources, coppicing was abandoned towards 1970. Since then, tree overaging has resulted in stand stagnation displayed by slow stem growth, branch dieback, and scarce acorn production. The urgent need to find new alternative uses for abandoned coppices is recognized as one of the biggest challenges which currently faces Mediterranean silviculture; conversion into high forest by thinning is one of the preferred alternatives. For this aim, thinning has been broadly applied and seldom tested, although without a comprehensive understanding of the causes of stand stagnation. In this PhD study, we test the hypothesis of an imbalance between above- and below-ground organs, result of long term coppicing, as the underlying cause of Q. pyrenaica decay. In an experimental plot coppiced since at least the 12th century, genetic analyses were performed a priori to elucidate inconspicuous clonal structure of Q. pyrenaica to evaluate how clonal size affects the functioning of these multi-stemmed trees. Clonal size negatively affected diametric stem growth, whereas root respiration rates, measured by internal fluxes of CO2 through xylem (FT) and soil CO2 efflux, increased with clonal size. These results suggest root-to-shoot (R:S) imbalance intensifying with clonal size: periodic removal of aboveground organs whilst belowground organs remain undisturbed may have led to massive root systems which consume a great proportion of non-structural carbohydrates (NSC) for maintenance respiration, thus constraining aboveground performance. Furthermore, excavation of two multi-stemmed trees, composed by four and eight stems, revealed R:S ratios (0.5 and 1, respectively) greater than those reported for sexually regenerated trees. Moreover, as similarly observed in several root-resprouting species, NSC allocation to roots was favored ([NSC] > 20% in spring): a large proportion of sapwood maintained throughout the root system (52%) stored a remarkable NSC pool of 87 kg in the case of the largest clone. In this root system of the eight-stemmed tree, 248 root connections were counted and, by radiocarbon dating, its age was estimated to be 550-years-old. Persistence of massive, old and highly interconnected root systems suggests that enhanced belowground NSC storage and consumption reflects a trade-off between vegetative resilience and aboveground development. For a better understanding of tree carbon budget and the potential role of carbon starvation in Q. pyrenaica decay, internal and external stem CO2 fluxes and soil CO2 effluxes were monitored to evaluate respiratory costs above- and below-ground. On a seasonal scale, stem and root respiration (RS and RR) mirrored sap flow and stem growth dynamics. Respiration was determined to the greatest extent by external fluxes of CO2 to the atmosphere or soil, since FT accounted for a low proportion of RS and RR (< 10% and < 2%, respectively). On a diel scale, the contribution of FT to RS increased up to 25% at high transpiration rates. Comparatively low FT was determined by the low concentration of xylem CO2 registered ([CO2] as low as 0.11%), likely as a consequence of constrained xylem respiration and reduced resistance to CO2 radial diffusion imposed by summer drought. Xylem [CO2] pulses following first autumn rains support this idea. Averaged over the growing season, soil CO2 efflux was the greatest respiratory flux (39 mol CO2 day-1), three and four times greater than RS (12 mol CO2 day-1) and RR (8-9 mol CO2 day-1), respectively. Ratios of RR/RS below one evidence an additional and important weight of aboveground respiration as a tree carbon sink. Finally, root trenching and stem girdling were tested as complimentary treatments to thinning as a means to improve carbon reserves in stems of clonal trees. Preliminary results discouraged root trenching due to the high cost likely incurred for wound closure. Stem girdling successfully blocked NSC translocation downward, increasing starch concentrations above the girdled zone whilst the root system is fed by non-girdled stems within the clone. Further measurements and ancillary data are necessary to verify that this positive effect hold over time. To conclude, the need of multidisciplinary approaches for an integrative understanding on the functioning of abandoned Q pyrenaica coppices is highlighted for an appropriate management of these stands.

Depth Averaged Models for Fast Landslide Propagation: Mathematical, Rheological and Numerical Aspects

This paper presents an overview of depth averaged modelling of fast catastrophic landslides where coupling of solid skeleton and pore fluid (air and water) is important. The first goal is to show how Biot-Zienkiewicz models can be applied to develop depth integrated, coupled models. The second objective of the paper is to consider a link which can be established between rheological and constitutive models. Perzyna´s viscoplasticity can be considered a general framework within which rheological models such as Bingham and cohesive frictional fluids can be derived. Among the several alternative numerical models, we will focus here on SPH which has not been widely applied by engineers to model landslide propagation. We propose an improvement, based on combining Finite Difference meshes associated to SPH nodes to describe pore pressure evolution inside the landslide mass. We devote a Section to analyze the performance of the models, considering three sets of tests and examples which allows to assess the model performance and limitations: (i) Problems having an analytical solution, (ii) Small scale laboratory tests, and (iii) Real cases for which we have had access to reliable information