Intra-annual patterns of tracheid size in the Mediterranean tree Juniperus thurifera as an indicator of seasonal water stress

Because climate can affect xylem cell anatomy, series of intra-annual cell anatomical features have the potential to retrospectively supply seasonal climatic information. In this study, we explored the ability to extract information about water stress conditions from tracheid features of the Mediterranean conifer Juniperus thurifera L. Tracheidograms of four climatic years from two drought-sensitive sites in Spain were compared to evaluate whether it is possible to link intra-annual cell size patterns to seasonal climatic conditions. Results indicated site-specific anatomical adjustment such as smaller and thicker tracheids at the dryer site but also showed a strong climatic imprint on the intra-annual pattern of tracheid size. Site differences in cell size reflected expected structural adjustments against cavitation failures. Differences between intra-annual patterns, however, indicated a response to seasonal changes in water availability whereby cells formed under drought conditions were smaller and thicker, and vice versa. This relationship was more manifest and stable at the dryer site

Spatial and temporal precipitation patterns of the Ebro River Basin, Spain

The Ebro River Basin, with around 85 000 km2 and located in NE Spain, is characterized by the high spatial heterogeneity of its geology, topography, climatology and land use. Rainfall is one of the most important climatic variables studied owing to its non-homogenous behaviour in event and intensity, which creates drought, water runoff and soil erosion with negative environmental and social consequences. In this work we characterized the rainfall variability pattern in the Ebro River Basin using universal multifractal (UM) analysis, which estimates the concentration of the data around the precipitation average (C1, codimension average), the degree of multiscaling behaviour in time (? index) and the maximum probable singularity in the rainfall distribution ( s). A spatial and temporal analysis of the UM parameters is applied to study the possible changes. With this porpoise, 60 daily rainfall series were selected from 132 synthetic series generated by Luna and Balairón (AEMet). These daily rainfall series present a length of 60 years, from 1950 to 2009. Each one of them was subdivided (1950?1970 and 1980?2009) to analyse the difference between the two periods. The range of variation of precipitation amounts and the frequency of dry events between both periods are discussed, as well as the evolution of the UM parameters through the years.

Origin and patterns of distribution of trace elements in street dust. Unleaded petrol and urban lead

The elemental composition, patterns of distribution and possible sources of street dust are not common to all urban environments, but vary according to the peculiarities of each city. The common features and dissimilarities in the origin and nature of street dust were investigated through a series of studies in two widely different cities, Madrid (Spain) and Oslo (Norway), between 1990 and 1994. The most comprehensive sampling campaign was carried out in the Norwegian capital during the summer of 1994. An area of 14 km2, covering most of downtown Oslo and some residential districts to the north of the city, was divided into 1 km2 mapping units, and 16 sampling increments of approximately 150 g were collected from streets and roads in each of them. The fraction below 100 μm was acid-digested and analysed by ICP-MS. Statistical analyses of the results suggest that chemical elements in street dust can be classified into three groups: “urban” elements (Ba, Cd, Co, Cu, Mg, Pb, Sb, Ti, Zn), “natural” elements (Al, Ga, La, Mn, Na, Sr, Th, Y) and elements of a mixed origin or which have undergone geochemical changes from their original sources (Ca, Cs, Fe, Mo, Ni, Rb, Sr, U). Soil resuspension and/or mobilisation appears to be the most important source of “natural” elements, while “urban” elements originate primarily from traffic and from the weathering and corrosion of building materials. The data for Pb seem to prove that the gradual shift from leaded to unleaded petrol as fuel for automobiles has resulted in an almost proportional reduction in the concentration of Pb in dust particles under 100 μm. This fact and the spatial distribution of Pb in the city strongly suggest that lead sources other than traffic (i.e. lead accumulated in urban soil over the years) may contribute as much lead, if not more, to urban street dust.

Spatial variability patterns of some Vertisol properties at a field scale using standardized data

Spatial variability of Vertisol properties is relevant for identifying those zones with physical degradation. In this sense, one has to face the problem of identifying the origin and distribution of spatial variability patterns. The objectives of the present work were (i) to quantify the spatial structure of different physical properties collected from a Vertisol, (ii) to search for potential correlations between different spatial patterns and (iii) to identify relevant components through multivariate spatial analysis. The study was conducted on a Vertisol (Typic Hapludert) dedicated to sugarcane (Saccharum officinarum L.) production during the last sixty years. We used six soil properties collected from a squared grid (225 points) (penetrometer resistance (PR), total porosity, fragmentation dimension (Df), vertical electrical conductivity (ECv), horizontal electrical conductivity (ECh) and soil water content (WC)). All the original data sets were z-transformed before geostatistical analysis. Three different types of semivariogram models were necessary for fitting individual experimental semivariograms. This suggests the different natures of spatial variability patterns. Soil water content rendered the largest nugget effect (C0 = 0.933) while soil total porosity showed the largest range of spatial correlation (A = 43.92 m). The bivariate geostatistical analysis also rendered significant cross-semivariance between different paired soil properties. However, four different semivariogram models were required in that case. This indicates an underlying co-regionalization between different soil properties, which is of interest for delineating management zones within sugarcane fields. Cross-semivariograms showed larger correlation ranges than individual, univariate, semivariograms (A ≥ 29 m). All the findings were supported by multivariate spatial analysis, which showed the influence of soil tillage operations, harvesting machinery and irrigation water distribution on the status of the investigated area.

Secundary school profile in science and technology versus the first year at the Technical University of Madrid. Cause of failure or quitting?

Se analiza el absentismo, el fallo y el abandono de los estudiantes en los primeros semestres del grado sobre la base de su formación en la educación secundaria.

Regional and seasonal patterns of litterfall in tropical South America

The production of aboveground soft tissue represents an important share of total net primary production in tropical rain forests. Here we draw from a large number of published and unpublished datasets (n = 81 sites) to assess the determinants of litterfall variation across South American tropical forests. We show that across old-growth tropical rainforests, litterfall averages 8.61±1.91Mgha?1 yr?1 (mean±standard deviation, in dry mass units). Secondary forests have a lower annual litterfall than old-growth tropical forests with a mean of 8.01±3.41Mgha?1 yr?1. Annual litterfall shows no significant variation with total annual rainfall, either globally or within forest types. It does not vary consistently with soil type, except in the poorest soils (white sand soils), where litterfall is significantly lower than in other soil types (5.42±1.91Mgha?1 yr?1). We also study the determinants of litterfall seasonality, and find that it does not depend on annual rainfall or on soil type. However, litterfall seasonality is significantly positively correlated with rainfall seasonality. Finally, we assess how much carbon is stored in reproductive organs relative to photosynthetic organs. Mean leaf fall is 5.74±1.83Mgha?1 yr?1 (71% of total litterfall). Mean allocation into reproductive organs is 0.69±0.40Mgha?1 yr?1 (9% of total litterfall). The investment into reproductive organs divided by leaf litterfall increases with soil fertility, suggesting that on poor soils, the allocation to photosynthetic organs is prioritized over that to reproduction. Finally, we discuss the ecological and biogeochemical implications of these results.

High fidelity simulations of failure in fiber-reinforced composites

Los ensayos virtuales de materiales compuestos han aparecido como un nuevo concepto dentro de la industria aeroespacial, y disponen de un vasto potencial para reducir los enormes costes de certificación y desarrollo asociados con las tediosas campañas experimentales, que incluyen un gran número de paneles, subcomponentes y componentes. El objetivo de los ensayos virtuales es sustituir algunos ensayos por simulaciones computacionales con alta fidelidad. Esta tesis es una contribución a la aproximación multiescala desarrollada en el Instituto IMDEA Materiales para predecir el comportamiento mecánico de un laminado de material compuesto dadas las propiedades de la lámina y la intercara. La mecánica de daño continuo (CDM) formula el daño intralaminar a nivel constitutivo de material. El modelo de daño intralaminar se combina con elementos cohesivos para representar daño interlaminar. Se desarrolló e implementó un modelo de daño continuo, y se aplicó a configuraciones simples de ensayos en laminados: impactos de baja y alta velocidad, ensayos de tracción, tests a cortadura. El análisis del método y la correlación con experimentos sugiere que los métodos son razonablemente adecuados para los test de impacto, pero insuficientes para el resto de ensayos. Para superar estas limitaciones de CDM, se ha mejorado la aproximación discreta de elementos finitos enriqueciendo la cinemática para incluir discontinuidades embebidas: el método extendido de los elementos finitos (X-FEM). Se adaptó X-FEM para un esquema explícito de integración temporal. El método es capaz de representar cualitativamente los mecanismos de fallo detallados en laminados. Sin embargo, los resultados muestran inconsistencias en la formulación que producen resultados cuantitativos erróneos. Por último, se ha revisado el método tradicional de X-FEM, y se ha desarrollado un nuevo método para superar sus limitaciones: el método cohesivo X-FEM estable. Las propiedades del nuevo método se estudiaron en detalle, y se concluyó que el método es robusto para implementación en códigos explícitos dinámicos escalables, resultando una nueva herramienta útil para la simulación de daño en composites. Virtual testing of composite materials has emerged as a new concept within the aerospace industry. It presents a very large potential to reduce the large certification costs and the long development times associated with the experimental campaigns, involving the testing of a large number of panels, sub-components and components. The aim of virtual testing is to replace some experimental tests by high-fidelity numerical simulations. This work is a contribution to the multiscale approach developed in Institute IMDEA Materials to predict the mechanical behavior of a composite laminate from the properties of the ply and the interply. Continuum Damage Mechanics (CDM) formulates intraply damage at the the material constitutive level. Intraply CDM is combined with cohesive elements to model interply damage. A CDM model was developed, implemented, and applied to simple mechanical tests of laminates: low and high velocity impact, tension of coupons, and shear deformation. The analysis of the results and the comparison with experiments indicated that the performance was reasonably good for the impact tests, but insuficient in the other cases. To overcome the limitations of CDM, the kinematics of the discrete finite element approximation was enhanced to include mesh embedded discontinuities, the eXtended Finite Element Method (X-FEM). The X-FEM was adapted to an explicit time integration scheme and was able to reproduce qualitatively the physical failure mechanisms in a composite laminate. However, the results revealed an inconsistency in the formulation that leads to erroneous quantitative results. Finally, the traditional X-FEM was reviewed, and a new method was developed to overcome its limitations, the stable cohesive X-FEM. The properties of the new method were studied in detail, and it was demonstrated that the new method was robust and can be implemented in a explicit finite element formulation, providing a new tool for damage simulation in composite materials.

Tailoring the magnetic anisotropy and domain patterns of sputtered TbFeGa alloys

We report the magnetic anisotropy and domain configuration of cosputtered TbFeGa alloys. The layers were deposited from two targets with compositions TbFe2 and Fe3Ga, respectively. The structural and magnetic properties do not only depend on the composition but also on the growth conditions. Alloys with the same composition but deposited using a DC or a pulsed power source in the TbFe2 target exhibit a different magnetic anisotropy. The perpendicular magnetic anisotropy, the size and topology of domain patterns can be tailored by changing the evaporation parameters of TbFe2. The width of the stripe domain increases from 235 to 835 nm when using the DC source in the TbFe2. We correlate this effect with Tb enrichment of the TbxFe1−x phases present in the samples.

Changes in spatial point patterns of pioneer woody plants across a large tropical landslide

We assessed whether the relative importance of positive and negative interactions in early successional communities varied across a large landslide on Casita Volcano (Nicaragua). We tested several hypotheses concerning the signatures of these processes in the spatial patterns of woody pioneer plants, as well as those of mortality and recruitment events, in several zones of the landslide differing in substrate stability and fertility, over a period of two years (2001 and 2002). We identified all woody individuals with a diameter >1 cm and mapped them in 28 plots measuring 10 × 10-m. On these maps, we performed a spatial point pattern analysis using univariate and bivariate pair-correlation functions; g (r) and g12 (r), and pairwise differences of univariate and bivariate functions. Spatial signatures of positive and negative interactions among woody plants were more prevalent in the most and least stressful zones of the landslide, respectively. Natural and human-induced disturbances such as the occurrence of fire, removal of newly colonizing plants through erosion and clearcutting of pioneer trees were also identified as potentially important pattern-creating processes. These results are in agreement with the stress-gradient hypothesis, which states that the relative importance of facilitation and competition varies inversely across gradients of abiotic stress. Our findings also indicate that the assembly of early successional plant communities in large heterogeneous landslides might be driven by a much larger array of processes than previously thought.

Maintaining distances with the engineer: patterns of coexistence in plant communities beyond the patch-bare dichotomy

Two-phase plant communities with an engineer conforming conspicuous patches and affecting the performance and patterns of coexisting species are the norm under stressful conditions. To unveil the mechanisms governing coexistence in these communities at multiple spatial scales, we have developed a new point-raster approach of spatial pattern analysis, which was applied to a Mediterranean high mountain grassland to show how Festuca curvifolia patches affect the local distribution of coexisting species. We recorded 22 111 individuals of 17 plant perennial species. Most coexisting species were negatively associated with F. curvifolia clumps. Nevertheless, bivariate nearest-neighbor analyses revealed that the majority of coexisting species were confined at relatively short distances from F. curvifolia borders (between 0-2 cm and up to 8 cm in some cases). Our study suggests the existence of a fine-scale effect of F. curvifolia for most species promoting coexistence through a mechanism we call 'facilitation in the halo'. Most coexisting species are displaced to an interphase area between patches, where two opposite forces reach equilibrium: attenuated severe conditions by proximity to the F. curvifolia canopy (nutrient-rich islands) and competitive exclusion mitigated by avoiding direct contact with F. curvifolia.

The gait patterns of an elite athlete with low vision: a case study

Presentación Póster en ESMAC 2012

New methods for the synthesis of radiation patterns of coupled antenna arrays = Nuevos métodos de síntesis de diagramas de radiación de agrupaciones de antenas acopladas

El principal objetivo de esta tesis es el desarrollo de métodos de síntesis de diagramas de radiación de agrupaciones de antenas, en donde se realiza una caracterización electromagnética rigurosa de los elementos radiantes y de los acoplos mutuos existentes. Esta caracterización no se realiza habitualmente en la gran mayoría de métodos de síntesis encontrados en la literatura, debido fundamentalmente a dos razones. Por un lado, se considera que el diagrama de radiación de un array de antenas se puede aproximar con el factor de array que únicamente tiene en cuenta la posición de los elementos y las excitaciones aplicadas a los mismos. Sin embargo, como se mostrará en esta tesis, en múltiples ocasiones un riguroso análisis de los elementos radiantes y del acoplo mutuo entre ellos es importante ya que los resultados obtenidos pueden ser notablemente diferentes. Por otro lado, no es sencillo combinar un método de análisis electromagnético con un proceso de síntesis de diagramas de radiación. Los métodos de análisis de agrupaciones de antenas suelen ser costosos computacionalmente, ya que son estructuras grandes en términos de longitudes de onda. Generalmente, un diseño de un problema electromagnético suele comprender varios análisis de la estructura, dependiendo de las variaciones de las características, lo que hace este proceso muy costoso. Dos métodos se utilizan en esta tesis para el análisis de los arrays acoplados. Ambos están basados en el método de los elementos finitos, la descomposición de dominio y el análisis modal para analizar la estructura radiante y han sido desarrollados en el grupo de investigación donde se engloba esta tesis. El primero de ellos es una técnica de análisis de arrays finitos basado en la aproximación de array infinito. Su uso es indicado para arrays planos de grandes dimensiones con elementos equiespaciados. El segundo caracteriza el array y el acoplo mutuo entre elementos a partir de una expansión en modos esféricos del campo radiado por cada uno de los elementos. Este método calcula los acoplos entre los diferentes elementos del array usando las propiedades de traslación y rotación de los modos esféricos. Es capaz de analizar agrupaciones de elementos distribuidos de forma arbitraria. Ambas técnicas utilizan una formulación matricial que caracteriza de forma rigurosa el campo radiado por el array. Esto las hace muy apropiadas para su posterior uso en una herramienta de diseño, como los métodos de síntesis desarrollados en esta tesis. Los resultados obtenidos por estas técnicas de síntesis, que incluyen métodos rigurosos de análisis, son consecuentemente más precisos. La síntesis de arrays consiste en modificar uno o varios parámetros de las agrupaciones de antenas buscando unas determinadas especificaciones de las características de radiación. Los parámetros utilizados como variables de optimización pueden ser varios. Los más utilizados son las excitaciones aplicadas a los elementos, pero también es posible modificar otros parámetros de diseño como son las posiciones de los elementos o las rotaciones de estos. Los objetivos de las síntesis pueden ser dirigir el haz o haces en una determinada dirección o conformar el haz con formas arbitrarias. Además, es posible minimizar el nivel de los lóbulos secundarios o del rizado en las regiones deseadas, imponer nulos que evitan posibles interferencias o reducir el nivel de la componente contrapolar. El método para el análisis de arrays finitos basado en la aproximación de array infinito considera un array finito como un array infinito con un número finito de elementos excitados. Los elementos no excitados están físicamente presentes y pueden presentar tres diferentes terminaciones, corto-circuito, circuito abierto y adaptados. Cada una de estas terminaciones simulará mejor el entorno real en el que el array se encuentre. Este método de análisis se integra en la tesis con dos métodos diferentes de síntesis de diagramas de radiación. En el primero de ellos se presenta un método basado en programación lineal en donde es posible dirigir el haz o haces, en la dirección deseada, además de ejercer un control sobre los lóbulos secundarios o imponer nulos. Este método es muy eficiente y obtiene soluciones óptimas. El mismo método de análisis es también aplicado a un método de conformación de haz, en donde un problema originalmente no convexo (y de difícil solución) es transformado en un problema convexo imponiendo restricciones de simetría, resolviendo de este modo eficientemente un problema complejo. Con este método es posible diseñar diagramas de radiación con haces de forma arbitraria, ejerciendo un control en el rizado del lóbulo principal, así como en el nivel de los lóbulos secundarios. El método de análisis de arrays basado en la expansión en modos esféricos se integra en la tesis con tres técnicas de síntesis de diagramas de radiación. Se propone inicialmente una síntesis de conformación del haz basado en el método de la recuperación de fase resuelta de forma iterativa mediante métodos convexos, en donde relajando las restricciones del problema original se consiguen unas soluciones cercanas a las óptimas de manera eficiente. Dos métodos de síntesis se han propuesto, donde las variables de optimización son las posiciones y las rotaciones de los elementos respectivamente. Se define una función de coste basada en la intensidad de radiación, la cual es minimizada de forma iterativa con el método del gradiente. Ambos métodos reducen el nivel de los lóbulos secundarios minimizando una función de coste. El gradiente de la función de coste es obtenido en términos de la variable de optimización en cada método. Esta función de coste está formada por la expresión rigurosa de la intensidad de radiación y por una función de peso definida por el usuario para imponer prioridades sobre las diferentes regiones de radiación, si así se desea. Por último, se presenta un método en el cual, mediante técnicas de programación entera, se buscan las fases discretas que generan un diagrama de radiación lo más cercano posible al deseado. Con este método se obtienen diseños que minimizan el coste de fabricación. En cada uno de las diferentes técnicas propuestas en la tesis, se presentan resultados con elementos reales que muestran las capacidades y posibilidades que los métodos ofrecen. Se comparan los resultados con otros métodos disponibles en la literatura. Se muestra la importancia de tener en cuenta los diagramas de los elementos reales y los acoplos mutuos en el proceso de síntesis y se comparan los resultados obtenidos con herramientas de software comerciales. ABSTRACT The main objective of this thesis is the development of optimization methods for the radiation pattern synthesis of array antennas in which a rigorous electromagnetic characterization of the radiators and the mutual coupling between them is performed. The electromagnetic characterization is usually overlooked in most of the available synthesis methods in the literature, this is mainly due to two reasons. On the one hand, it is argued that the radiation pattern of an array is mainly influenced by the array factor and that the mutual coupling plays a minor role. As it is shown in this thesis, the mutual coupling and the rigorous characterization of the array antenna influences significantly in the array performance and its computation leads to differences in the results obtained. On the other hand, it is difficult to introduce an analysis procedure into a synthesis technique. The analysis of array antennas is generally expensive computationally as the structure to analyze is large in terms of wavelengths. A synthesis method requires to carry out a large number of analysis, this makes the synthesis problem very expensive computationally or intractable in some cases. Two methods have been used in this thesis for the analysis of coupled antenna arrays, both of them have been developed in the research group in which this thesis is involved. They are based on the finite element method (FEM), the domain decomposition and the modal analysis. The first one obtains a finite array characterization with the results obtained from the infinite array approach. It is specially indicated for the analysis of large arrays with equispaced elements. The second one characterizes the array elements and the mutual coupling between them with a spherical wave expansion of the radiated field by each element. The mutual coupling is computed using the properties of translation and rotation of spherical waves. This method is able to analyze arrays with elements placed on an arbitrary distribution. Both techniques provide a matrix formulation that makes them very suitable for being integrated in synthesis techniques, the results obtained from these synthesis methods will be very accurate. The array synthesis stands for the modification of one or several array parameters looking for some desired specifications of the radiation pattern. The array parameters used as optimization variables are usually the excitation weights applied to the array elements, but some other array characteristics can be used as well, such as the array elements positions or rotations. The desired specifications may be to steer the beam towards any specific direction or to generate shaped beams with arbitrary geometry. Further characteristics can be handled as well, such as minimize the side lobe level in some other radiating regions, to minimize the ripple of the shaped beam, to take control over the cross-polar component or to impose nulls on the radiation pattern to avoid possible interferences from specific directions. The analysis method based on the infinite array approach considers an infinite array with a finite number of excited elements. The infinite non-excited elements are physically present and may have three different terminations, short-circuit, open circuit and match terminated. Each of this terminations is a better simulation for the real environment of the array. This method is used in this thesis for the development of two synthesis methods. In the first one, a multi-objective radiation pattern synthesis is presented, in which it is possible to steer the beam or beams in desired directions, minimizing the side lobe level and with the possibility of imposing nulls in the radiation pattern. This method is very efficient and obtains optimal solutions as it is based on convex programming. The same analysis method is used in a shaped beam technique in which an originally non-convex problem is transformed into a convex one applying symmetry restrictions, thus solving a complex problem in an efficient way. This method allows the synthesis of shaped beam radiation patterns controlling the ripple in the mainlobe and the side lobe level. The analysis method based on the spherical wave expansion is applied for different synthesis techniques of the radiation pattern of coupled arrays. A shaped beam synthesis is presented, in which a convex formulation is proposed based on the phase retrieval method. In this technique, an originally non-convex problem is solved using a relaxation and solving a convex problems iteratively. Two methods are proposed based on the gradient method. A cost function is defined involving the radiation intensity of the coupled array and a weighting function that provides more degrees of freedom to the designer. The gradient of the cost function is computed with respect to the positions in one of them and the rotations of the elements in the second one. The elements are moved or rotated iteratively following the results of the gradient. A highly non-convex problem is solved very efficiently, obtaining very good results that are dependent on the starting point. Finally, an optimization method is presented where discrete digital phases are synthesized providing a radiation pattern as close as possible to the desired one. The problem is solved using linear integer programming procedures obtaining array designs that greatly reduce the fabrication costs. Results are provided for every method showing the capabilities that the above mentioned methods offer. The results obtained are compared with available methods in the literature. The importance of introducing a rigorous analysis into the synthesis method is emphasized and the results obtained are compared with a commercial software, showing good agreement.