Short communication. Collection and characterisation of a population of Triticum boeoticum Boiss., a wild wheat species not previously found in the Mediterranean western region

En julio del 2010 se recolectó una población de trigo silvestre en una zona abandonada cerca de Madrid, España. Esta zona posee una biodiversidad botánica elevada y un tipo de suelo muy peculiar denominado “arcillas verdes”. Se recogió una muestra de trigo y se multiplicó y caracterizó para varios caracteres agro-morfológicos y subunidades de gluteninas. El número cromosómico 2n de las semillas demostró que es una especie diploide de trigo y los datos de caracterización indicaron que es Triticum boeoticum Boiss. Esta especie llegó probablemente como mala hierba del cultivo de escaña que se producía en la zona hasta al menos la primera mitad del s. xix. Las características edáficas y climáticas del lugar y el hecho de que no haya referencias hasta ahora de esta especie en la zona oeste de la región Mediterránea aumentan el valor de esta adquisición para la mejora del trigo. La nueva accesión se conserva en el Centro Nacional de Recursos Fitogenéticos y se ha depositado una hoja de herbario en el Real Jardín Botánico de Madrid

Desarrollo de un modelo para la caracterización termoeconómica de ciclos combinados de turbinas de gas y de vapor en condiciones de carga variable

En la actualidad, el interés por las plantas de potencia de ciclo combinado de gas y vapor ha experimentado un notable aumento debido a su alto rendimiento, bajo coste de generación y rápida construcción. El objetivo fundamental de la tesis es profundizar en el conocimiento de esta tecnología, insuficientemente conocida hasta el momento debido al gran número de grados de libertad que existen en el diseño de este tipo de instalaciones. El estudio se realizó en varias fases. La primera consistió en analizar y estudiar las distintas tecnologías que se pueden emplear en este tipo de centrales, algunas muy recientes o en fase de investigación, como las turbinas de gas de geometría variable, las turbinas de gas refrigeradas con agua o vapor del ciclo de vapor o las calderas de paso único que trabajan con agua en condiciones supercríticas. Posteriormente se elaboraron los modelos matemáticos que permiten la simulación termodinámica de cada uno de los componentes que integran las plantas, tanto en el punto de diseño como a cargas parciales. Al mismo tiempo, se desarrolló una metodología novedosa que permite resolver el sistema de ecuaciones que resulta de la simulación de cualquier configuración posible de ciclo combinado. De esa forma se puede conocer el comportamiento de cualquier planta en cualquier punto de funcionamiento. Por último se desarrolló un modelo de atribución de costes para este tipo de centrales. Con dicho modelo, los estudios se pueden realizar no sólo desde un punto de vista termodinámico sino también termoeconómico, con lo que se pueden encontrar soluciones de compromiso entre rendimiento y coste, asignar costes de producción, determinar curvas de oferta, beneficios económicos de la planta y delimitar el rango de potencias donde la planta es rentable. El programa informático, desarrollado en paralelo con los modelos de simulación, se ha empleado para obtener resultados de forma intensiva. El estudio de los resultados permite profundizar ampliamente en el conocimiento de la tecnología y, así, desarrollar una metodología de diseño de este tipo de plantas bajo un criterio termoeconómico. ABSTRACT The growing energy demand and the need of shrinking costs have led to the design of high efficiency and quick installation power plants. The success of combined cycle gas turbine power plants lies on their high efficiency, low cost and short construction lead time. The main objective of the work is to study in detail this technology, which is not thoroughly known owing to the great number of degrees of freedom that exist in the design of this kind of power plants. The study is divided into three parts. Firstly, the different technologies and components that could be used in any configuration of a combined cycle gas turbine power plant are studied. Some of them could be of recent technology, such as the variable inlet guide vane compressors, the H-technology for gas turbine cooling or the once-through heat recovery steam generators, used with water at supercritical conditions. Secondly, a mathematical model has been developed to simulate at full and part load the components of the power plant. At the same time, a new methodology is proposed in order to solve the equation system resulting for any possible power plant configuration. Therefore, any combined cycle gas turbine could be simulated at any part load condition. Finally a themoeconomic model is proposed. This model allows studying the power plant not only from a thermodynamic point of view but also from a thermoeconomic one. Likewise, it allows determining the generating costs or the cash flow, thus achieving a trade off between efficiency and cost. Likewise, the model calculates the part load range where the power plant is profitable. Once the thermodynamic and thermoeconomic models are developed, they are intensively used in order to gain knowledge in the combined cycle gas turbine technology and, in this way, to propose a methodology aimed at the design of this kind of power plants from a thermoeconomic point of view.

Influence of local wind speed and direction on wind power dynamics - Application to offshore very short-term forecasting

Wind power time series usually show complex dynamics mainly due to non-linearities related to the wind physics and the power transformation process in wind farms. This article provides an approach to the incorporation of observed local variables (wind speed and direction) to model some of these effects by means of statistical models. To this end, a benchmarking between two different families of varying-coefficient models (regime-switching and conditional parametric models) is carried out. The case of the offshore wind farm of Horns Rev in Denmark has been considered. The analysis is focused on one-step ahead forecasting and a time series resolution of 10 min. It has been found that the local wind direction contributes to model some features of the prevailing winds, such as the impact of the wind direction on the wind variability, whereas the non-linearities related to the power transformation process can be introduced by considering the local wind speed. In both cases, conditional parametric models showed a better performance than the one achieved by the regime-switching strategy. The results attained reinforce the idea that each explanatory variable allows the modelling of different underlying effects in the dynamics of wind power time series.

Evolution and treatment of vitamin B12 deficiency as a risk factor for (cognitive and functional) neurodegenerative diseases in institutionalized elderly = Evolución y tratamiento de la deficiencia de vitamina B12 como factor de riesgo de enfermedades neurodegenerativas (cognitivas y funcionales) en las personas mayores institucionalizadas.

Prevalence of vitamin B12 deficiency is very common in elderly people and can reach values as high as 40.5% of the population. It can be the result of the interaction among several factors. Vitamin B12 deficiencies have been associated with neurological, cognitive deterioration, haematological abnormalities and cardiovascular diseases that have an important influence on the health of the elderly and their quality of life. It is necessary to approach the problems arisen from the lack of data relative to them. The main objective of this thesis was to analyse the evolution of vitamin B12 status and related parameters, lipid and haematological profiles and their relationship to health risk factors, and to functional and cognitive status over one year and to determine the effect of an oral supplementation of 500 μg of cyanocobalamin for a short period of 28 days. An additional objective was to analyze the possible effects of medicine intakes on vitamin B status. Three studies were performed: a) a one year longitudinal follow-up with four measure points; b) an intervention study providing an oral liquid supplement of 500 μg of cyanocobalamin for a 28 days period; and c) analysis of the possible effect of medication intake on vitamin B status using the ATC classification of medicines. The participants for these studies were recruited from nursing homes for the elderly in the Region of Madrid. Sixty elders (mean age 84 _ 7y, 19 men and 41 women) were recruited for Study I and 64 elders (mean age 82 _ 7y, 24 men and 40 women) for Study II. For Study III, baseline data from the initially recruited participants of the first two studies were used. An informed consent was obtained from all participants or their mentors. The studies were approved by the Ethical Committee of the University of Granada. Blood samples were obtained at each examination date and were analyzed for serum cobalamin, holoTC, serum and RBC folate and total homocysteine according to laboratory standard procedures. The haematological parameters analyzed were haematocrit, haemoglobin and MCV. For the lipid profile TG, total cholesterol, LDL- and HDLcholesterol were analyzed. Anthropometric measures (BMI, skinfolds [triceps and subscapular], waist girth and waist to hip ratio), functional tests (hand grip, arm and leg strength tests, static balance) and MMSE were obtained or administered by trained personal. The vitamin B12 supplement of Study II was administered with breakfast and the medication intake was taken from the residents’ anamnesis. Data were analyzed by parametric and non-parametric statistics depending on the obtained data. Comparisons were done using the appropriate ANOVAs or non-parametric tests. Pearsons’ partial correlations with the variable “time” as control were used to define the association of the analyzed parameters. XIII The results showed that: A) Over one year, in relationship to vitamin B status, serum cobalamin decreased, serum folate and mean corpuscular volumen increased significantly and total homocysteine concentrations were stable. Regarding blood lipid profile, triglycerides increased and HDL-cholesterol decreased significantly. Regarding selected anthropometric measurements, waist circumference increased significantly. No significant changes were observed for the rest of parameters. B) Prevalence of hyperhomocysteinemia was high in the elderly studied, ranging from 60% to 90 % over the year depending on the cut-off used for the classification. LDL-cholesterol values were high, especially among women, and showed a tendency to increase over the year. Results of the balance test showed a deficiency and a tendency to decrease; this indicates that the population studied is at high risk for falls. Lower extremity muscular function was deficient and showed a tendency to decrease. A highly significant relationship was observed between the skinfold of the triceps and blood lipid profile. C) Low cobalamin concentrations correlated significantly with low MMSE scores in the elderly studied. No correlations were observed between vitamin B12 status and functional parameters. D) Regarding vitamin B12 status, holo-transcobalamin seems to be more sensitive for diagnosis; 5-10% of the elderly had a deficiency using serum cobalamin as a criterion, and 45-52% had a deficiency when using serum holotranscobalamin as a criterion. E) 500 μg of cyanocobalamin administered orally during 28 days significantly improved vitamin B12 status and significantly decreased total homocysteine concentrations in institutionalized elderly. No effect of the intervention was observed on functional and cognitive parameters. F) The relative change (%) of improvement of vitamin B12 status was higher when using serum holo-transcobalamin as a criterion than serum cobalamin. G) Antiaenemic drug intake normalized cobalamin, urologic drugs and corticosteroids serum folate, and psychoanaleptics holo-transcobalamin levels. Drugs treating pulmonary obstruction increased total homocysteine concentration significantly. H) The daily mean drug intake was 5.1. Fiftynine percent of the elderly took medication belonging to 5 or more different ATC groups. The most prevalent were psycholeptic (53%), antiacid (53%) and antithrombotic (47%) drugs.

Propagation mechanism modelling in the near region of circular tunnels

Artículo sobre comunicaciones ferroviarias. Abstract: Along with the increase in operating frequencies in advanced radio communication systems utilised inside tunnels, the location of the break point is further and further away from the transmitter. This means that the near region lengthens considerably and even occupies the whole propagation cell or the entire length of some short tunnels. To begin with, this study analyses the propagation loss resulting from the free-space mechanism and the multi-mode waveguide mechanism in the near region of circular tunnels, respectively. Then, by conjunctive employing the propagation theory and the three-dimensional solid geometry, a general analytical model of the dividing point between two propagation mechanisms is presented for the first time. Moreover, the model is validated by a wide range of measurement campaigns in different tunnels at different frequencies. Finally, discussions on the simplified formulae of the dividing point in some application situations are made. The results in this study can be helpful to grasp the essence of the propagation mechanism inside tunnels.

Propagation Mechanism modeling in the Near-Region of Arbitrary Cross-Sectional Tunnels.

Along with the increase of the use of working frequencies in advanced radio communication systems, the near-region inside tunnels lengthens considerably and even occupies the whole propagation cell or the entire length of some short tunnels. This paper analytically models the propagation mechanisms and their dividing point in the near-region of arbitrary cross-sectional tunnels for the first time. To begin with, the propagation losses owing to the free space mechanism and the multimode waveguide mechanism are modeled, respectively. Then, by conjunctively employing the propagation theory and the three-dimensional solid geometry, the paper presents a general model for the dividing point between two propagation mechanisms. It is worthy to mention that this model can be applied in arbitrary cross-sectional tunnels. Furthermore, the general dividing point model is specified in rectangular, circular, and arched tunnels, respectively. Five groups of measurements are used to justify the model in different tunnels at different frequencies. Finally, in order to facilitate the use of the model, simplified analytical solutions for the dividing point in five specific application situations are derived. The results in this paper could help deepen the insight into the propagation mechanisms in tunnels.

Sustainability of irrigated farming systems in a Tunisian region: A recursive stochastic programming analysis

The aim of this study was to evaluate the sustainability of farm irrigation systems in the Cébalat district in northern Tunisia. It addressed the challenging topic of sustainable agriculture through a bio-economic approach linking a biophysical model to an economic optimisation model. A crop growth simulation model (CropSyst) was used to build a database to determine the relationships between agricultural practices, crop yields and environmental effects (salt accumulation in soil and leaching of nitrates) in a context of high climatic variability. The database was then fed into a recursive stochastic model set for a 10-year plan that allowed analysing the effects of cropping patterns on farm income, salt accumulation and nitrate leaching. We assumed that the long-term sustainability of soil productivity might be in conflict with farm profitability in the short-term. Assuming a discount rate of 10% (for the base scenario), the model closely reproduced the current system and allowed to predict the degradation of soil quality due to long-term salt accumulation. The results showed that there was more accumulation of salt in the soil for the base scenario than for the alternative scenario (discount rate of 0%). This result was induced by applying a higher quantity of water per hectare for the alternative as compared to a base scenario. The results also showed that nitrogen leaching is very low for the two discount rates and all climate scenarios. In conclusion, the results show that the difference in farm income between the alternative and base scenarios increases over time to attain 45% after 10 years.

Evaluation and improvement of empirical models of global solar irradiation: Case study northern Spain

This paper presents a new methodology to build parametric models to estimate global solar irradiation adjusted to specific on-site characteristics based on the evaluation of variable im- portance. Thus, those variables higly correlated to solar irradiation on a site are implemented in the model and therefore, different models might be proposed under different climates. This methodology is applied in a study case in La Rioja region (northern Spain). A new model is proposed and evaluated on stability and accuracy against a review of twenty-two already exist- ing parametric models based on temperatures and rainfall in seventeen meteorological stations in La Rioja. The methodology of model evaluation is based on bootstrapping, which leads to achieve a high level of confidence in model calibration and validation from short time series (in this case five years, from 2007 to 2011). The model proposed improves the estimates of the other twenty-two models with average mean absolute error (MAE) of 2.195 MJ/m2 day and average confidence interval width (95% C.I., n=100) of 0.261 MJ/m2 day. 41.65% of the daily residuals in the case of SIAR and 20.12% in that of SOS Rioja fall within the uncertainty tolerance of the pyranometers of the two networks (10% and 5%, respectively). Relative differences between measured and estimated irradiation on an annual cumulative basis are below 4.82%. Thus, the proposed model might be useful to estimate annual sums of global solar irradiation, reaching insignificant differences between measurements from pyranometers.

Structure Variation And Dynamics of the Nuclear Ribosomal Intergenic Spacer Region in Key Gibberella Fujikuroi Complex Species

The intergenic spacer (IGS) region of the ribosomal DNA was cloned and sequenced in eight species within the Gibberella fujikuroi species complex with anamorphs in the genus Fusarium , a group that includes the most relevant toxigenic species. DNA sequence analyses revealed two categories of repeated elements: long repeats and short repeats of 125 and 8 bp, respectively. Long repeats were present in two copies and were conserved in all the species analyzed, whereas different numbers of short repeat elements were observed, leading to species-specific IGS sequences with different length. In Fusarium subglutinans and Fusarium nygamai , these differences seemed to be the result of duplication and deletion events. Here, we propose a model based on unequal crossing over that can explain these processes. The partial IGS sequence of 22 Fusarium proliferatum isolates was also obtained to study variation at the intraspecific level. The results revealed no differences in terms of number or pattern of repeated elements and detected frequent gene conversion events. These results suggest that the homogenization observed at the intraspecific level might not be achieved primarily by unequal crossing-over events but rather by processes associated with recombination such as gene conversion events.

Short-term effects of four tillage practices on soil physical properties, soil water potential, and maize yield

The area cultivated using conservation tillage has recently increased in central Spain. However, soil compaction and water retention with conservation tillage still remains a genuine concern for landowners in this region be- cause of its potential effect on the crop growth and yield. The aim of this research is to determine the short- term influences of four tillage treatments on soil physical properties. In the experiment, bulk density, cone index, soil water potential, soil temperature and maize (Zea mays L.) productivity have been measured. A field experiment was established in spring of 2013 on a loamy soil. The experiment compared four tillage methods (zero tillage, ZT; reservoir tillage, RT; minimum tillage, MT; and conventional tillage, CT). Soil bulk density and soil cone index were measured during maize growing season and at harvesting time. Furthermore, the soil water potential was monitored by using a wireless sensors network with sensors at 20 and 40 cm depths. Also, soil temperatures were registered at depths of 5 and 12 cm. Results indicated that there were significant differ- ences between soil bulk density and cone index of ZT method and those of RT, MT, and CT, during the growing season; although, this difference was not significant at the time of harvesting in some soil layers. Overall, in most soil layers, tillage practice affected bulk density and cone index in the order: ZT N RT N MT N CT. Regardless oftheentireobservationperiod,results exhibited that soils under ZT and RT treatments usually resulted in higher water potential and lower soil temperature than the other two treatments at both soil depths. In addition, clear differences in maize grain yield were observed between ZT and CT treatments, with a grain yield (up to 15.4%) increase with the CT treatment. On the other hand, no significant differences among (RT, MT, and CT) on maizeyieldwerefound.Inconclusion,the impact of soil compaction increase and soil temperature decrease,pro- duced by ZT treatment is a potential reason for maize yield reduction in this tillage method. We found that RT could be certainly a viable option for farmers incentral Spain,particularly when switching to conservation tillage from conventional tillage. This technique showed a moderate and positive effect on soil physical properties and increased maize yields compared to ZT and MT, and provides an opportunity to stabilize maize yields compared to CT.

Variable-density jet flows induced by concentrated sources of momentum and energy

The planar and axisymmetric variable-density flows induced in a quiescent gas by a concentrated source of momentum that is simultaneously either a source or a sink of energy are investigated for application to the description of the velocity and temperature far fields in laminar gaseous jets with either large or small values of the initial jet-to-ambient temperature ratio. The source fluxes of momentum and heat are used to construct the characteristic scales of velocity and length in the region where the density differences are of the order of the ambient density, which is slender for the large values of the Reynolds number considered herein. The problem reduces to the integration of the dimensionless boundary-layer conservation equations, giving a solution that depends on the gas transport properties but is otherwise free of parameters. The boundary conditions at the jet exit for integration are obtained by analysing the self-similar flow that appears near the heat source in planar and axisymmetric configurations and also near the heat sink in the planar case. Numerical integrations of the boundary-layer equations with these conditions give solutions that describe accurately the velocity and temperature fields of very hot planar and round jets and also of very cold plane jets in the far field region where the density and temperature differences are comparable to the ambient values. Simple scaling arguments indicate that the point source description does not apply, however, to cold round jets, whose far field region is not large compared with the jet development region, as verified by numerical integrations

Railway Rapid Transit timetables with variable and elastic demand.

This paper focuses on the design of railway timetables considering a variable elastic demand profile along a whole design day. Timetabling is the third stage in the classical hierarchical railway planning process. Most of previous works on this topic consider a uniform demand behavior for short planning intervals. In this paper, we propose a MINLP model for designing non-periodic timetables on a railway corridor where demand is dependent on waiting times. In the elastic demand case, long waiting times lead to a loss of passengers, who may select an alternative transportation mode. The mode choice is modeled using two alternative methods. The first one is based on a sigmoid function and can be used in case of absence of information for competitor modes. In the second one, the mode choice probability is obtained using a Logit model that explicitly considers the existence of a main alternative mode. With the purpose of obtaining optimal departure times, in both cases, a minimization of the loss of passengers is used as objective function. Finally, as illustration, the timetabling MINLP model with both mode choice methods is applied to a real case and computational results are shown.

Temporally flat top pulse generation from gain switched semiconductor lasers based on a polarization interferometer with variable transfer function

We propose the use of a polarization based interferometer with variable transfer function for the generation of temporally flat top pulses from gain switched single mode semiconductor lasers. The main advantage of the presented technique is its flexibility in terms of input pulse characteristics, as pulse duration, spectral bandwidth and operating wavelength. Theoretical predictions and experimental demonstrations are presented and the proposed technique is applied to two different semiconductor laser sources emitting in the 1550 nm region. Flat top pulses are successfully obtained with input seed pulses with duration ranging from 40 ps to 100 ps.

Methodologies for Collision Risk Computation = Metodologías para Cálculo de Riesgo de Colisión

Esta tesis aborda metodologías para el cálculo de riesgo de colisión de satélites. La minimización del riesgo de colisión se debe abordar desde dos puntos de vista distintos. Desde el punto de vista operacional, es necesario filtrar los objetos que pueden presentar un encuentro entre todos los objetos que comparten el espacio con un satélite operacional. Puesto que las órbitas, del objeto operacional y del objeto envuelto en la colisión, no se conocen perfectamente, la geometría del encuentro y el riesgo de colisión deben ser evaluados. De acuerdo con dicha geometría o riesgo, una maniobra evasiva puede ser necesaria para evitar la colisión. Dichas maniobras implican un consumo de combustible que impacta en la capacidad de mantenimiento orbital y por tanto de la visa útil del satélite. Por tanto, el combustible necesario a lo largo de la vida útil de un satélite debe ser estimado en fase de diseño de la misión para una correcta definición de su vida útil, especialmente para satélites orbitando en regímenes orbitales muy poblados. Los dos aspectos, diseño de misión y aspectos operacionales en relación con el riesgo de colisión están abordados en esta tesis y se resumen en la Figura 3. En relación con los aspectos relacionados con el diseño de misión (parte inferior de la figura), es necesario evaluar estadísticamente las características de de la población espacial y las teorías que permiten calcular el número medio de eventos encontrados por una misión y su capacidad de reducir riesgo de colisión. Estos dos aspectos definen los procedimientos más apropiados para reducir el riesgo de colisión en fase operacional. Este aspecto es abordado, comenzando por la teoría descrita en [Sánchez-Ortiz, 2006]T.14 e implementada por el autor de esta tesis en la herramienta ARES [Sánchez-Ortiz, 2004b]T.15 proporcionada por ESA para la evaluación de estrategias de evitación de colisión. Esta teoría es extendida en esta tesis para considerar las características de los datos orbitales disponibles en las fases operacionales de un satélite (sección 4.3.3). Además, esta teoría se ha extendido para considerar riesgo máximo de colisión cuando la incertidumbre de las órbitas de objetos catalogados no es conocida (como se da el caso para los TLE), y en el caso de querer sólo considerar riesgo de colisión catastrófico (sección 4.3.2.3). Dichas mejoras se han incluido en la nueva versión de ARES [Domínguez-González and Sánchez-Ortiz, 2012b]T.12 puesta a disposición a través de [SDUP,2014]R.60. En fase operacional, los catálogos que proporcionan datos orbitales de los objetos espaciales, son procesados rutinariamente, para identificar posibles encuentros que se analizan en base a algoritmos de cálculo de riesgo de colisión para proponer maniobras de evasión. Actualmente existe una única fuente de datos públicos, el catálogo TLE (de sus siglas en inglés, Two Line Elements). Además, el Joint Space Operation Center (JSpOC) Americano proporciona mensajes con alertas de colisión (CSM) cuando el sistema de vigilancia americano identifica un posible encuentro. En función de los datos usados en fase operacional (TLE o CSM), la estrategia de evitación puede ser diferente debido a las características de dicha información. Es preciso conocer las principales características de los datos disponibles (respecto a la precisión de los datos orbitales) para estimar los posibles eventos de colisión encontrados por un satélite a lo largo de su vida útil. En caso de los TLE, cuya precisión orbital no es proporcionada, la información de precisión orbital derivada de un análisis estadístico se puede usar también en el proceso operacional así como en el diseño de la misión. En caso de utilizar CSM como base de las operaciones de evitación de colisiones, se conoce la precisión orbital de los dos objetos involucrados. Estas características se han analizado en detalle, evaluando estadísticamente las características de ambos tipos de datos. Una vez concluido dicho análisis, se ha analizado el impacto de utilizar TLE o CSM en las operaciones del satélite (sección 5.1). Este análisis se ha publicado en una revista especializada [Sánchez-Ortiz, 2015b]T.3. En dicho análisis, se proporcionan recomendaciones para distintas misiones (tamaño del satélite y régimen orbital) en relación con las estrategias de evitación de colisión para reducir el riesgo de colisión de manera significativa. Por ejemplo, en el caso de un satélite en órbita heliosíncrona en régimen orbital LEO, el valor típico del ACPL que se usa de manera extendida es 10-4. Este valor no es adecuado cuando los esquemas de evitación de colisión se realizan sobre datos TLE. En este caso, la capacidad de reducción de riesgo es prácticamente nula (debido a las grandes incertidumbres de los datos TLE) incluso para tiempos cortos de predicción. Para conseguir una reducción significativa del riesgo, sería necesario usar un ACPL en torno a 10-6 o inferior, produciendo unas 10 alarmas al año por satélite (considerando predicciones a un día) o 100 alarmas al año (con predicciones a tres días). Por tanto, la principal conclusión es la falta de idoneidad de los datos TLE para el cálculo de eventos de colisión. Al contrario, usando los datos CSM, debido a su mejor precisión orbital, se puede obtener una reducción significativa del riesgo con ACPL en torno a 10-4 (considerando 3 días de predicción). Incluso 5 días de predicción pueden ser considerados con ACPL en torno a 10-5. Incluso tiempos de predicción más largos se pueden usar (7 días) con reducción del 90% del riesgo y unas 5 alarmas al año (en caso de predicciones de 5 días, el número de maniobras se mantiene en unas 2 al año). La dinámica en GEO es diferente al caso LEO y hace que el crecimiento de las incertidumbres orbitales con el tiempo de propagación sea menor. Por el contrario, las incertidumbres derivadas de la determinación orbital son peores que en LEO por las diferencias en las capacidades de observación de uno y otro régimen orbital. Además, se debe considerar que los tiempos de predicción considerados para LEO pueden no ser apropiados para el caso de un satélite GEO (puesto que tiene un periodo orbital mayor). En este caso usando datos TLE, una reducción significativa del riesgo sólo se consigue con valores pequeños de ACPL, produciendo una alarma por año cuando los eventos de colisión se predicen a un día vista (tiempo muy corto para implementar maniobras de evitación de colisión).Valores más adecuados de ACPL se encuentran entre 5•10-8 y 10-7, muy por debajo de los valores usados en las operaciones actuales de la mayoría de las misiones GEO (de nuevo, no se recomienda en este régimen orbital basar las estrategias de evitación de colisión en TLE). Los datos CSM permiten una reducción de riesgo apropiada con ACPL entre 10-5 y 10-4 con tiempos de predicción cortos y medios (10-5 se recomienda para predicciones a 5 o 7 días). El número de maniobras realizadas sería una en 10 años de misión. Se debe notar que estos cálculos están realizados para un satélite de unos 2 metros de radio. En el futuro, otros sistemas de vigilancia espacial (como el programa SSA de la ESA), proporcionarán catálogos adicionales de objetos espaciales con el objetivo de reducir el riesgo de colisión de los satélites. Para definir dichos sistemas de vigilancia, es necesario identificar las prestaciones del catalogo en función de la reducción de riesgo que se pretende conseguir. Las características del catálogo que afectan principalmente a dicha capacidad son la cobertura (número de objetos incluidos en el catalogo, limitado principalmente por el tamaño mínimo de los objetos en función de las limitaciones de los sensores utilizados) y la precisión de los datos orbitales (derivada de las prestaciones de los sensores en relación con la precisión de las medidas y la capacidad de re-observación de los objetos). El resultado de dicho análisis (sección 5.2) se ha publicado en una revista especializada [Sánchez-Ortiz, 2015a]T.2. Este análisis no estaba inicialmente previsto durante la tesis, y permite mostrar como la teoría descrita en esta tesis, inicialmente definida para facilitar el diseño de misiones (parte superior de la figura 1) se ha extendido y se puede aplicar para otros propósitos como el dimensionado de un sistema de vigilancia espacial (parte inferior de la figura 1). La principal diferencia de los dos análisis se basa en considerar las capacidades de catalogación (precisión y tamaño de objetos observados) como una variable a modificar en el caso de un diseño de un sistema de vigilancia), siendo fijas en el caso de un diseño de misión. En el caso de las salidas generadas en el análisis, todos los aspectos calculados en un análisis estadístico de riesgo de colisión son importantes para diseño de misión (con el objetivo de calcular la estrategia de evitación y la cantidad de combustible a utilizar), mientras que en el caso de un diseño de un sistema de vigilancia, los aspectos más importantes son el número de maniobras y falsas alarmas (fiabilidad del sistema) y la capacidad de reducción de riesgo (efectividad del sistema). Adicionalmente, un sistema de vigilancia espacial debe ser caracterizado por su capacidad de evitar colisiones catastróficas (evitando así in incremento dramático de la población de basura espacial), mientras que el diseño de una misión debe considerar todo tipo de encuentros, puesto que un operador está interesado en evitar tanto las colisiones catastróficas como las letales. Del análisis de las prestaciones (tamaño de objetos a catalogar y precisión orbital) requeridas a un sistema de vigilancia espacial se concluye que ambos aspectos han de ser fijados de manera diferente para los distintos regímenes orbitales. En el caso de LEO se hace necesario observar objetos de hasta 5cm de radio, mientras que en GEO se rebaja este requisito hasta los 100 cm para cubrir las colisiones catastróficas. La razón principal para esta diferencia viene de las diferentes velocidades relativas entre los objetos en ambos regímenes orbitales. En relación con la precisión orbital, ésta ha de ser muy buena en LEO para poder reducir el número de falsas alarmas, mientras que en regímenes orbitales más altos se pueden considerar precisiones medias. En relación con los aspectos operaciones de la determinación de riesgo de colisión, existen varios algoritmos de cálculo de riesgo entre dos objetos espaciales. La Figura 2 proporciona un resumen de los casos en cuanto a algoritmos de cálculo de riesgo de colisión y como se abordan en esta tesis. Normalmente se consideran objetos esféricos para simplificar el cálculo de riesgo (caso A). Este caso está ampliamente abordado en la literatura y no se analiza en detalle en esta tesis. Un caso de ejemplo se proporciona en la sección 4.2. Considerar la forma real de los objetos (caso B) permite calcular el riesgo de una manera más precisa. Un nuevo algoritmo es definido en esta tesis para calcular el riesgo de colisión cuando al menos uno de los objetos se considera complejo (sección 4.4.2). Dicho algoritmo permite calcular el riesgo de colisión para objetos formados por un conjunto de cajas, y se ha presentado en varias conferencias internacionales. Para evaluar las prestaciones de dicho algoritmo, sus resultados se han comparado con un análisis de Monte Carlo que se ha definido para considerar colisiones entre cajas de manera adecuada (sección 4.1.2.3), pues la búsqueda de colisiones simples aplicables para objetos esféricos no es aplicable a este caso. Este análisis de Monte Carlo se considera la verdad a la hora de calcular los resultados del algoritmos, dicha comparativa se presenta en la sección 4.4.4. En el caso de satélites que no se pueden considerar esféricos, el uso de un modelo de la geometría del satélite permite descartar eventos que no son colisiones reales o estimar con mayor precisión el riesgo asociado a un evento. El uso de estos algoritmos con geometrías complejas es más relevante para objetos de dimensiones grandes debido a las prestaciones de precisión orbital actuales. En el futuro, si los sistemas de vigilancia mejoran y las órbitas son conocidas con mayor precisión, la importancia de considerar la geometría real de los satélites será cada vez más relevante. La sección 5.4 presenta un ejemplo para un sistema de grandes dimensiones (satélite con un tether). Adicionalmente, si los dos objetos involucrados en la colisión tienen velocidad relativa baja (y geometría simple, Caso C en la Figura 2), la mayor parte de los algoritmos no son aplicables requiriendo implementaciones dedicadas para este caso particular. En esta tesis, uno de estos algoritmos presentado en la literatura [Patera, 2001]R.26 se ha analizado para determinar su idoneidad en distintos tipos de eventos (sección 4.5). La evaluación frete a un análisis de Monte Carlo se proporciona en la sección 4.5.2. Tras este análisis, se ha considerado adecuado para abordar las colisiones de baja velocidad. En particular, se ha concluido que el uso de algoritmos dedicados para baja velocidad son necesarios en función del tamaño del volumen de colisión proyectado en el plano de encuentro (B-plane) y del tamaño de la incertidumbre asociada al vector posición entre los dos objetos. Para incertidumbres grandes, estos algoritmos se hacen más necesarios pues la duración del intervalo en que los elipsoides de error de los dos objetos pueden intersecar es mayor. Dicho algoritmo se ha probado integrando el algoritmo de colisión para objetos con geometrías complejas. El resultado de dicho análisis muestra que este algoritmo puede ser extendido fácilmente para considerar diferentes tipos de algoritmos de cálculo de riesgo de colisión (sección 4.5.3). Ambos algoritmos, junto con el método Monte Carlo para geometrías complejas, se han implementado en la herramienta operacional de la ESA CORAM, que es utilizada para evaluar el riesgo de colisión en las actividades rutinarias de los satélites operados por ESA [Sánchez-Ortiz, 2013a]T.11. Este hecho muestra el interés y relevancia de los algoritmos desarrollados para la mejora de las operaciones de los satélites. Dichos algoritmos han sido presentados en varias conferencias internacionales [Sánchez-Ortiz, 2013b]T.9, [Pulido, 2014]T.7,[Grande-Olalla, 2013]T.10, [Pulido, 2014]T.5, [Sánchez-Ortiz, 2015c]T.1. ABSTRACT This document addresses methodologies for computation of the collision risk of a satellite. Two different approaches need to be considered for collision risk minimisation. On an operational basis, it is needed to perform a sieve of possible objects approaching the satellite, among all objects sharing the space with an operational satellite. As the orbits of both, satellite and the eventual collider, are not perfectly known but only estimated, the miss-encounter geometry and the actual risk of collision shall be evaluated. In the basis of the encounter geometry or the risk, an eventual manoeuvre may be required to avoid the conjunction. Those manoeuvres will be associated to a reduction in the fuel for the mission orbit maintenance, and thus, may reduce the satellite operational lifetime. Thus, avoidance manoeuvre fuel budget shall be estimated, at mission design phase, for a better estimation of mission lifetime, especially for those satellites orbiting in very populated orbital regimes. These two aspects, mission design and operational collision risk aspects, are summarised in Figure 3, and covered along this thesis. Bottom part of the figure identifies the aspects to be consider for the mission design phase (statistical characterisation of the space object population data and theory computing the mean number of events and risk reduction capability) which will define the most appropriate collision avoidance approach at mission operational phase. This part is covered in this work by starting from the theory described in [Sánchez-Ortiz, 2006]T.14 and implemented by this author in ARES tool [Sánchez-Ortiz, 2004b]T.15 provided by ESA for evaluation of collision avoidance approaches. This methodology has been now extended to account for the particular features of the available data sets in operational environment (section 4.3.3). Additionally, the formulation has been extended to allow evaluating risk computation approached when orbital uncertainty is not available (like the TLE case) and when only catastrophic collisions are subject to study (section 4.3.2.3). These improvements to the theory have been included in the new version of ESA ARES tool [Domínguez-González and Sánchez-Ortiz, 2012b]T.12 and available through [SDUP,2014]R.60. At the operation phase, the real catalogue data will be processed on a routine basis, with adequate collision risk computation algorithms to propose conjunction avoidance manoeuvre optimised for every event. The optimisation of manoeuvres in an operational basis is not approached along this document. Currently, American Two Line Element (TLE) catalogue is the only public source of data providing orbits of objects in space to identify eventual conjunction events. Additionally, Conjunction Summary Message (CSM) is provided by Joint Space Operation Center (JSpOC) when the American system identifies a possible collision among satellites and debris. Depending on the data used for collision avoidance evaluation, the conjunction avoidance approach may be different. The main features of currently available data need to be analysed (in regards to accuracy) in order to perform estimation of eventual encounters to be found along the mission lifetime. In the case of TLE, as these data is not provided with accuracy information, operational collision avoidance may be also based on statistical accuracy information as the one used in the mission design approach. This is not the case for CSM data, which includes the state vector and orbital accuracy of the two involved objects. This aspect has been analysed in detail and is depicted in the document, evaluating in statistical way the characteristics of both data sets in regards to the main aspects related to collision avoidance. Once the analysis of data set was completed, investigations on the impact of those features in the most convenient avoidance approaches have been addressed (section 5.1). This analysis is published in a peer-reviewed journal [Sánchez-Ortiz, 2015b]T.3. The analysis provides recommendations for different mission types (satellite size and orbital regime) in regards to the most appropriate collision avoidance approach for relevant risk reduction. The risk reduction capability is very much dependent on the accuracy of the catalogue utilized to identify eventual collisions. Approaches based on CSM data are recommended against the TLE based approach. Some approaches based on the maximum risk associated to envisaged encounters are demonstrated to report a very large number of events, which makes the approach not suitable for operational activities. Accepted Collision Probability Levels are recommended for the definition of the avoidance strategies for different mission types. For example for the case of a LEO satellite in the Sun-synchronous regime, the typically used ACPL value of 10-4 is not a suitable value for collision avoidance schemes based on TLE data. In this case the risk reduction capacity is almost null (due to the large uncertainties associated to TLE data sets, even for short time-to-event values). For significant reduction of risk when using TLE data, ACPL on the order of 10-6 (or lower) seems to be required, producing about 10 warnings per year and mission (if one-day ahead events are considered) or 100 warnings per year (for three-days ahead estimations). Thus, the main conclusion from these results is the lack of feasibility of TLE for a proper collision avoidance approach. On the contrary, for CSM data, and due to the better accuracy of the orbital information when compared with TLE, ACPL on the order of 10-4 allows to significantly reduce the risk. This is true for events estimated up to 3 days ahead. Even 5 days ahead events can be considered, but ACPL values down to 10-5 should be considered in such case. Even larger prediction times can be considered (7 days) for risk reduction about 90%, at the cost of larger number of warnings up to 5 events per year, when 5 days prediction allows to keep the manoeuvre rate in 2 manoeuvres per year. Dynamics of the GEO orbits is different to that in LEO, impacting on a lower increase of orbits uncertainty along time. On the contrary, uncertainties at short prediction times at this orbital regime are larger than those at LEO due to the differences in observation capabilities. Additionally, it has to be accounted that short prediction times feasible at LEO may not be appropriate for a GEO mission due to the orbital period being much larger at this regime. In the case of TLE data sets, significant reduction of risk is only achieved for small ACPL values, producing about a warning event per year if warnings are raised one day in advance to the event (too short for any reaction to be considered). Suitable ACPL values would lay in between 5•10-8 and 10-7, well below the normal values used in current operations for most of the GEO missions (TLE-based strategies for collision avoidance at this regime are not recommended). On the contrary, CSM data allows a good reduction of risk with ACPL in between 10-5 and 10-4 for short and medium prediction times. 10-5 is recommended for prediction times of five or seven days. The number of events raised for a suitable warning time of seven days would be about one in a 10-year mission. It must be noted, that these results are associated to a 2 m radius spacecraft, impact of the satellite size are also analysed within the thesis. In the future, other Space Situational Awareness Systems (SSA, ESA program) may provide additional catalogues of objects in space with the aim of reducing the risk. It is needed to investigate which are the required performances of those catalogues for allowing such risk reduction. The main performance aspects are coverage (objects included in the catalogue, mainly limited by a minimum object size derived from sensor performances) and the accuracy of the orbital data to accurately evaluate the conjunctions (derived from sensor performance in regards to object observation frequency and accuracy). The results of these investigations (section 5.2) are published in a peer-reviewed journal [Sánchez-Ortiz, 2015a]T.2. This aspect was not initially foreseen as objective of the thesis, but it shows how the theory described in the thesis, initially defined for mission design in regards to avoidance manoeuvre fuel allocation (upper part of figure 1), is extended and serves for additional purposes as dimensioning a Space Surveillance and Tracking (SST) system (bottom part of figure below). The main difference between the two approaches is the consideration of the catalogue features as part of the theory which are not modified (for the satellite mission design case) instead of being an input for the analysis (in the case of the SST design). In regards to the outputs, all the features computed by the statistical conjunction analysis are of importance for mission design (with the objective of proper global avoidance strategy definition and fuel allocation), whereas for the case of SST design, the most relevant aspects are the manoeuvre and false alarm rates (defining a reliable system) and the Risk Reduction capability (driving the effectiveness of the system). In regards to the methodology for computing the risk, the SST system shall be driven by the capacity of providing the means to avoid catastrophic conjunction events (avoiding the dramatic increase of the population), whereas the satellite mission design should consider all type of encounters, as the operator is interested on avoiding both lethal and catastrophic collisions. From the analysis of the SST features (object coverage and orbital uncertainty) for a reliable system, it is concluded that those two characteristics are to be imposed differently for the different orbital regimes, as the population level is different depending on the orbit type. Coverage values range from 5 cm for very populated LEO regime up to 100 cm in the case of GEO region. The difference on this requirement derives mainly from the relative velocity of the encounters at those regimes. Regarding the orbital knowledge of the catalogues, very accurate information is required for objects in the LEO region in order to limit the number of false alarms, whereas intermediate orbital accuracy can be considered for higher orbital regimes. In regards to the operational collision avoidance approaches, several collision risk algorithms are used for evaluation of collision risk of two pair of objects. Figure 2 provides a summary of the different collision risk algorithm cases and indicates how they are covered along this document. The typical case with high relative velocity is well covered in literature for the case of spherical objects (case A), with a large number of available algorithms, that are not analysed in detailed in this work. Only a sample case is provided in section 4.2. If complex geometries are considered (Case B), a more realistic risk evaluation can be computed. New approach for the evaluation of risk in the case of complex geometries is presented in this thesis (section 4.4.2), and it has been presented in several international conferences. The developed algorithm allows evaluating the risk for complex objects formed by a set of boxes. A dedicated Monte Carlo method has also been described (section 4.1.2.3) and implemented to allow the evaluation of the actual collisions among a large number of simulation shots. This Monte Carlo runs are considered the truth for comparison of the algorithm results (section 4.4.4). For spacecrafts that cannot be considered as spheres, the consideration of the real geometry of the objects may allow to discard events which are not real conjunctions, or estimate with larger reliability the risk associated to the event. This is of particular importance for the case of large spacecrafts as the uncertainty in positions of actual catalogues does not reach small values to make a difference for the case of objects below meter size. As the tracking systems improve and the orbits of catalogued objects are known more precisely, the importance of considering actual shapes of the objects will become more relevant. The particular case of a very large system (as a tethered satellite) is analysed in section 5.4. Additionally, if the two colliding objects have low relative velocity (and simple geometries, case C in figure above), the most common collision risk algorithms fail and adequate theories need to be applied. In this document, a low relative velocity algorithm presented in the literature [Patera, 2001]R.26 is described and evaluated (section 4.5). Evaluation through comparison with Monte Carlo approach is provided in section 4.5.2. The main conclusion of this analysis is the suitability of this algorithm for the most common encounter characteristics, and thus it is selected as adequate for collision risk estimation. Its performances are evaluated in order to characterise when it can be safely used for a large variety of encounter characteristics. In particular, it is found that the need of using dedicated algorithms depend on both the size of collision volume in the B-plane and the miss-distance uncertainty. For large uncertainties, the need of such algorithms is more relevant since for small uncertainties the encounter duration where the covariance ellipsoids intersect is smaller. Additionally, its application for the case of complex satellite geometries is assessed (case D in figure above) by integrating the developed algorithm in this thesis with Patera’s formulation for low relative velocity encounters. The results of this analysis show that the algorithm can be easily extended for collision risk estimation process suitable for complex geometry objects (section 4.5.3). The two algorithms, together with the Monte Carlo method, have been implemented in the operational tool CORAM for ESA which is used for the evaluation of collision risk of ESA operated missions, [Sánchez-Ortiz, 2013a]T.11. This fact shows the interest and relevance of the developed algorithms for improvement of satellite operations. The algorithms have been presented in several international conferences, [Sánchez-Ortiz, 2013b]T.9, [Pulido, 2014]T.7,[Grande-Olalla, 2013]T.10, [Pulido, 2014]T.5, [Sánchez-Ortiz, 2015c]T.1.

Compatibilidad de Eretmocerus mundus Mercet (Hymenoptera: Aphelinidae) y Amblyseius swirskii Athias-Henriot (Acari: Phytoseiidae), importantes enemigos naturales de la mosca blanca Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) en cultivos hortícolas, con nuevas barreras físicas selectivas y modernos plaguicidas

Los programas de Gestión Integrada de Plagas (GIP) promueven el uso de estrategias de control que sean respetuosas con el medio ambiente, sin embargo el uso de insecticidas en los cultivos hortícolas sigue siendo necesario para el control de determinadas plagas, como es el caso de la mosca blanca Bemisia tabaci (Gennadius). Por ello, el objetivo de esta tesis es el estudio de la integración de las tres estrategias de control más empleadas hoy en día para el control de plagas: el control biológico, el físico y el químico. Una primera parte de este trabajo ha consistido en el estudio de los efectos letales y subletales de once insecticidas, aplicados a la dosis máxima de campo, sobre los enemigos naturales Eretmocerus mundus Mercet y Amblyseius swirskii Athias-Henriot, mediante ensayos de laboratorio y persistencia (laboratorio extendido). Para la evaluación de la toxicidad de los insecticidas sobre los estados de vida más protegidos de estos enemigos naturales, se trataron bajo la Torre de Potter las pupas de E. mundus y los huevos de A. swirskii. Además, se llevaron a cabo ensayos de contacto residual para determinar los efectos letales y subletales de estos insecticidas sobre el estado adulto de ambas especies de enemigos naturales. Para ello, los pesticidas se aplicaron sobre placas de cristal (laboratorio) o sobre plantas (laboratorio extendido: persistencia). Los resultados mostraron que los insecticidas flonicamida, flubendiamida, metaflumizona, metoxifenocida, spiromesifen y spirotetramat eran compatibles con el estado de pupa de E. mundus (OILB 1: Inocuos). Sin embargo, abamectina, deltametrina y emamectina fueron categorizadas como ligeramente tóxicas (OILB 2) al causar efectos deletéreos. Los dos pesticidas más tóxicos fueron spinosad y sulfoxaflor, los cuales redujeron significativamente la emergencia de las pupas tratadas (OILB 4: Tóxicos). Flonicamida, flubendiamida, metoxifenocida y spiromesifen fueron compatibles con el estado adulto de E. mundus (OILB 1: Inocuos). Abamectina, deltametrina, emamectina, metaflumizona y spiromesifen pueden ser recomendados para su uso en programas de GIP, si se usan los plazos de seguridad apropiados, de acuerdo con la persistencia de cada uno de estos insecticidas, antes de la liberación del enemigo natural. Al contrario, spinosad y sulfoxaflor no resultaron ser compatibles (OILB D: Persistentes), aunque la realización de ensayos adicionales es necesaria para ver los efectos de los mismos en campo. Todos los insecticidas estudiados, excepto el spirotetramat (OILB 2: Ligeramente tóxico), fueron selectivos para el estado de huevo de A. swirskii (OILB 1: Inocuos). Flonicamida, flubendiamida, metaflumizona, metoxifenocida, spiromesifen, spirotetramat y sulfoxaflor, fueron compatibles con el estado adulto de A. swirskii (OILB 1: Inocuos). Abamectina, deltametrina, emamectina y spinosad pueden ser recomendados para su uso en programas de GIP, si se usan los plazos de seguridad apropiados, de acuerdo con la persistencia de cada uno de estos insecticidas, antes de la liberación del enemigo natural. Entre las nuevas estrategias de la GIP, los plásticos y mallas fotoselectivas han demostrado ser una herramienta importante para el control de plagas y enfermedades en cultivos hortícolas protegidos. Por ello, en una segunda parte de este trabajo, se estudiaron tanto los efectos directos, como la combinación de efectos directos y mediados por planta y plaga de ambientes pobres en luz UV, en presencia o ausencia del Virus del rizado amarillo del tomate (TYLCV), sobre E. mundus. En primer lugar, se realizó un ensayo al aire libre para la evaluación de la capacidad de vuelo de E. mundus en cajas tipo túnel (1 x 0,6 x 0,6 m) cubiertas con distintas barreras absorbentes de luz UV. Se detectó un efecto directo en la capacidad de orientación de E. mundus, debido a que este parasitoide utiliza estímulos visuales para localizar a sus huéspedes, únicamente en las barreras que bloqueaban más del 65% de la luz UV (malla G). En segundo lugar, bajo condiciones de invernadero, se evaluó la combinación de efectos directos y mediados por planta y plaga sobre E. mundus, usando plantas de tomate sanas o infectadas con el TYLCV y cajas (30 x 30 x 60 cm) cubiertas con los distintos plásticos fotoselectivos. En este caso, no se observó ningún efecto en la capacidad benéfica del parasitoide cuando este estaba en contacto con plantas de tomate infestadas con ninfas de B. tabaci, lo que demuestra que este insecto usa estímulos táctiles para encontrar a sus huéspedes a cortas distancias. Además, las diferentes condiciones de radiación UV estudiadas tuvieron cierto impacto en la morfología, fisiología y bioquímica de las plantas de tomate, infestadas o no con el virus de la cuchara, detectándose pequeñas alteraciones en alguno de los parámetros estudiados, como el peso fresco y seco, el contenido en H y el espesor de las cutículas y de las paredes celulares de la epidermis foliar. Por último, no se observaron efectos de la radiación UV mediados por planta, ni en B. tabaci ni en su parasitoide, E. mundus. En una tercera parte, se evaluaron los efectos de una malla tratada con bifentrin sobre ambos enemigos naturales, en ensayos de laboratorio, semicampo y campo. Las mallas tratadas fueron diseñadas originariamente para el control de mosquitos vectores de la malaria, y actualmente se está trabajando para su uso en agricultura, como una nueva estrategia de control de plagas. En ensayos de laboratorio, cuando adultos de E. mundus y A. swirskii se expusieron por contacto durante 72 horas con la malla tratada (cajas de 6 cm diámetro), se registró una alta mortalidad. Sin embargo, en el ensayo de preferencia, estos enemigos naturales no fueron capaces de detectar la presencia de bifentrin y, en aquellos individuos forzados a atravesar la malla tratada, no se observó mortalidad a corto plazo (72 horas). En estudios de semicampo, llevados a cabo bajo condiciones de invernadero en cajas de 25 x 25 x 60 cm de altura, la capacidad benéfica de E. mundus no se vio afectada. Finalmente, en ensayos de campo llevados a cabo en invernaderos comerciales (4000m2) en Almería, A. swirskii no se vio afectado por la presencia en el cultivo de la malla tratada con bifentrin y los niveles de infestación de B. tabaci y F. occidentalis detectados bajo dicha malla, fueron inferiores a los del control. Por último, se ha evaluado la composición de la microflora bacteriana de tres especies de parasitoides, E. mundus, Eretmocerus eremicus Rose & Zolnerowich y Encarsia formosa Gahan, y la influencia de la misma en su susceptibilidad a insecticidas. Se llevó a cabo una extracción total de ADN de los insectos y la región variable V4 del ARNr se amplificó usando cebadores universales bacterianos. Para identificar las secuencias de los géneros bacterianos presentes en los parasitoides, se realizó una Next Generation sequencing (Illumina sequencing). Una vez identificados los géneros bacterianos, el gen ADNr 16S de las Actinobacterias se amplificó del ADN extraído de los insectos, usando cebadores universales bacterianos y específicos de Actinobacterias, y los productos de la Nested PCR fueron clonados para identificar todas las especies del género Arthrobacter. Tres bacterias (A. aurescens Phillips, A. nicotinovarans Kodama, Yamamoto, Amano and Amichi y A. uratoxydans Stackebrandt, Fowler, Fiedler and Seiler), próximas a las especies de Arthrobacter presentes en los parasitoides, se obtuvieron de la colección bacteriana del BCCMTM/LMG y se midió su actividad esterasa. Finalmente, se realizaron ensayos con antibióticos (tetraciclina) y de contacto residual con insecticidas (abamectina) para determinar la influencia de las especies de Arthrobacter en la susceptibilidad de E. mundus a insecticidas. Los resultados muestran que este género bacteriano puede afectar a la toxicidad de E. mundus a abamectina, mostrando la importancia de la comunidad microbiana en enemigos naturales, factor que debe ser considerado en los estudios de evaluación de los riesgos de los insecticidas. ABSTRACT Integrated Pest Management (IPM) programs promote the use of control strategies more respectful with the environment; however the use of insecticides in vegetable crops is still needed to control certain pests, such as the whitefly Bemisia tabaci (Gennadius). Therefore, the objective of this work is to study the integration of the three most commonly used pest control strategies nowadays: biological, physical and chemical control. Firstly, the lethal and sublethal effects of eleven insecticides, applied at their maximum field recommended concentration, on the parasitic wasp Eretmocerus mundus Mercet and the predator Amblyseius swirskii Athias-Henriot has been assessed in the laboratory and in persistence tests (extended laboratory). To test the effects of pesticides on the most protected life stage of these natural enemies, E. mundus pupae and A. swirskii eggs were sprayed under a Potter precision spray tower. Laboratory contact tests were therefore conducted to determine the lethal and sublethal effects of these pesticides on the adult stage of these natural enemies. In the residual contact tests the pesticides were applied on glass plates (laboratory) or plants (extended laboratory: persistence). The study showed that the insecticides flonicamid, flubendiamide, metaflumizone, methoxyfenozide, spiromesifen and spirotetramat were selective for E. mundus pupae (IOBC 1: Harmless). Nevertheless, abamectin, deltamethrin and emamectin were categorized as slightly harmful (IOBC 2) due to the deleterious effects caused. The two most harmful pesticides were spinosad and sulfoxaflor, which significantly reduced the adult emergence from treated pupae (IOBC 4: Harmful). Flonicamid, flubendiamide, methoxyfenozide and spiromesifen were compatible with E. mundus adults (IOBC 1: Harmless). Base on the duration of the harmful activity, abamectin, deltamethrin, emamectin, metaflumizone and spirotetramat could be recommended for use in IPM programs if appropriate safety deadlines are used before the natural enemy release. On the contrary, spinosad and sulfoxaflor were not compatible (IOBC D: persistent), although additional studies are required to determine their effects under field conditions. All the pesticides tested, except spirotetramat (IOBC 2: Slightly harmful), were selective for A. swirskii eggs (IOBC 1: Harmless). Flonicamid, flubendiamide, metaflumizone, methoxyfenozide, spiromesifen, spirotetramat and sulfoxaflor were compatible with A. swirskii adults (IOBC 1: Harmless). However, abamectin, deltamethrin, emamectin and spinosad could be recommended for use in IPM programs if appropriate safety deadlines are used before the natural enemy release. Among new IPM strategies, UV-absorbing photoselective plastic films and nets have been shown to be an important tool for the control of pests and diseases in horticultural protected crops. Because of that, we secondly studied the plant and pest insect-mediated and/or the direct effects on E. mundus under different UV radiation conditions, in presence or absence of the Tomato Yellow Leaf Curl Virus (TYLCV). In the first experiment, performed outdoors, the flight activity of E. mundus was studied in one-chamber tunnels (1 x 0.6 x 0.6 m) covered with different photoselective barriers. Because E. mundus uses visual cues for host location at a long distance, a direct effect on its host location ability was detected, but only in the UV-absorbing barriers blocking more than 65% of the UV light (G net). In a second experiment, the direct and plant and pest insect-mediated effects of different UV radiation conditions on E. mundus were studied, inside cages (30 x 30 x 60 cm) covered with the different UVplastic films and under greenhouse conditions, using healthy or TYLCV-virus infected tomato plants. In this case, not any effect on the beneficial capacity of this parasitoid was detected, proving that he uses tactile cues at a short distance of the host. Moreover, the different UV radiation conditions studied had a certain direct impact in the morphology, physiology and biochemistry of tomato plants infested or not with the TYLCV, and small alterations in some parameters such as fresh and dry weight, H percentage and cuticle and cell wall thickness of epidermal cells of the leaves, were detected. Finally, none plant-mediated UV effects neither in the whitefly B. tabaci nor in their parasitic wasp were found. Thirdly, the effects of a bifenthrin treated net were evaluated in different laboratory, semi-field and field experiments on the natural enemies studied. Treated nets were developed long time ago aiming at the control of the mosquitoes vectors of malaria, and nowadays, there is a great interest on assessing the possibility of their use in agriculture. In laboratory assays, a high mortality was recorded when E. mundus and A. swirskii adults were exposed by contact to the bifenthrin treated net for 72 hours in small cages (12 cm diameter). However, these natural enemies were not able to detect the presence of bifenthrin in a dual-choice test and no short-term mortality (72 hours) was recorded in those individuals that went through the treated net. In semi-field assays, performed under greenhouse conditions with cages of 25 x 25 x 60 cm high, the beneficial capacity of E. mundus was not affected. Finally, in field assays carried out in commercial multispan greenhouses (4000 m2) in Almería, A. swirskii was not affected by the presence of the bifenthrin treated net in the crop and the B. tabaci and F. occidentalis infestation levels were significantly lower than in the control. Finally, the composition of the microflora present in three species of parasitoids, E. mundus, Eretmocerus eremicus Rose & Zolnerowich and Encarsia formosa Gahan, and its influence in their susceptibility to insecticides, have been assessed. A total DNA extraction was performed on insects and universal bacterial primers were used to amplify the variable V4 region of the rRNA. A Next Generation sequencing (Illumina sequencing) was performed to identify the sequences of the bacterial genera present in the parasitic wasps. Once, the bacterial genera were identified, 16S rDNA gene of Actinobacteria were amplified from insects DNA extracts using the universal bacterial and actinobacterial primers, and the nested PCR products, were cloned to identify the Arthrobacter species. Three bacteria (A. aurescens Phillips, A. nicotinovarans Kodama, Yamamoto, Amano and Amichi and A. uratoxydans Stackebrandt, Fowler, Fiedler and Seiler), having the closest match with the Arthrobacter species present in the parasitic wasps, were obtained from the BCCMTM/LMG bacteria collection and its esterase activity was measured. Finally, antibiotic and residual contact tests were done to determine the influence of Arthrobacter species in the susceptibility of E. mundus to pesticides (abamectin). The results suggest that this bacterial genus can affect the toxicity of E. mundus to abamectin, which in turn supports the importance of the microbial community in natural enemies that it should be considered as a factor in risk assessment tests of pesticides.