LiDAR simulation in modelled orchards to optimise the use of terretrial laser scanners and derived vegetative measures

Light detection and ranging (LiDAR) technology is beginning to have an impact on agriculture. Canopy volume and/or fruit tree leaf area can be estimated using terrestrial laser sensors based on this technology. However, the use of these devices may have different options depending on the resolution and scanning mode. As a consequence, data accuracy and LiDAR derived parameters are affected by sensor configuration, and may vary according to vegetative characteristics of tree crops. Given this scenario, users and suppliers of these devices need to know how to use the sensor in each case. This paper presents a computer program to determine the best configuration, allowing simulation and evaluation of different LiDAR configurations in various tree structures (or training systems). The ultimate goal is to optimise the use of laser scanners in field operations. The software presented generates a virtual orchard, and then allows the scanning simulation with a laser sensor. Trees are created using a hidden Markov tree (HMT) model. Varying the foliar structure of the orchard the LiDAR simulation was applied to twenty different artificially created orchards with or without leaves from two positions (lateral and zenith). To validate the laser sensor configuration, leaf surface of simulated trees was compared with the parameters obtained by LiDAR measurements: the impacted leaf area, the impacted total area (leaves and wood), and th impacted area in the three outer layers of leaves.

Aplicaciones de interés forestal de las imágenes obtenidas con el sensor Multi-angle Imaging SpectroRadiometer (MISR)

Las aplicaciones de la teledetección al seguimiento de lo que ocurre en la superficie terrestre se han ido multiplicando y afinando con el lanzamiento de nuevos sensores por parte de las diferentes agencias espaciales. La necesidad de tener información actualizada cada poco tiempo y espacialmente homogénea, ha provocado el desarrollo de nuevos programas como el Earth Observing System (EOS) de la National Aeronautics and Space Administration (NASA). Uno de los sensores que incorpora el buque insignia de ese programa, el satélite TERRA, es el Multi-angle Imaging SpectroRadiometer (MISR), diseñado para capturar información multiangular de la superficie terrestre. Ya desde los años 1970, se conocía que la reflectancia de las diversas ocupaciones y usos del suelo variaba en función del ángulo de observación y de iluminación, es decir, que eran anisotrópicas. Tal variación estaba además relacionada con la estructura tridimensional de tales ocupaciones, por lo que se podía aprovechar tal relación para obtener información de esa estructura, más allá de la que pudiera proporcionar la información meramente espectral. El sensor MISR incorpora 9 cámaras a diferentes ángulos para capturar 9 imágenes casi simultáneas del mismo punto, lo que permite estimar con relativa fiabilidad la respuesta anisotrópica de la superficie terrestre. Varios trabajos han demostrado que se pueden estimar variables relacionadas con la estructura de la vegetación con la información que proporciona MISR. En esta Tesis se ha realizado una primera aplicación a la Península Ibérica, para comprobar su utilidad a la hora de estimar variables de interés forestal. En un primer paso se ha analizado la variabilidad temporal que se produce en los datos, debido a los cambios en la geometría de captación, es decir, debido a la posición relativa de sensores y fuente de iluminación, que en este caso es el Sol. Se ha comprobado cómo la anisotropía es mayor desde finales de otoño hasta principios de primavera debido a que la posición del Sol es más cercana al plano de los sensores. También se ha comprobado que los valores máximo y mínimo se van desplazando temporalmente entre el centro y el extremo angular. En la caracterización multiangular de ocupaciones del suelo de CORINE Land Cover que se ha realizado, se puede observar cómo la forma predominante en las imágenes con el Sol más alto es convexa con un máximo en la cámara más cercana a la fuente de iluminación. Sin embargo, cuando el Sol se encuentra mucho más bajo, ese máximo es muy externo. Por otra parte, los datos obtenidos en verano son mucho más variables para cada ocupación que los de noviembre, posiblemente debido al aumento proporcional de las zonas en sombra. Para comprobar si la información multiangular tiene algún efecto en la obtención de imágenes clasificadas según ocupación y usos del suelo, se han realizado una serie de clasificaciones variando la información utilizada, desde sólo multiespectral, a multiangular y multiespectral. Los resultados muestran que, mientras para las clasificaciones más genéricas la información multiangular proporciona los peores resultados, a medida que se amplían el número de clases a obtener tal información mejora a lo obtenido únicamente con información multiespectral. Por otra parte, se ha realizado una estimación de variables cuantitativas como la fracción de cabida cubierta (Fcc) y la altura de la vegetación a partir de información proporcionada por MISR a diferentes resoluciones. En el valle de Alcudia (Ciudad Real) se ha estimado la fracción de cabida cubierta del arbolado para un píxel de 275 m utilizando redes neuronales. Los resultados muestran que utilizar información multiespectral y multiangular puede mejorar casi un 20% las estimaciones realizadas sólo con datos multiespectrales. Además, las relaciones obtenidas llegan al 0,7 de R con errores inferiores a un 10% en Fcc, siendo éstos mucho mejores que los obtenidos con el producto elaborado a partir de datos multiespectrales del sensor Moderate Resolution Imaging Spectroradiometer (MODIS), también a bordo de Terra, para la misma variable. Por último, se ha estimado la fracción de cabida cubierta y la altura efectiva de la vegetación para 700.000 ha de la provincia de Murcia, con una resolución de 1.100 m. Los resultados muestran la relación existente entre los datos espectrales y los multiangulares, obteniéndose coeficientes de Spearman del orden de 0,8 en el caso de la fracción de cabida cubierta de la vegetación, y de 0,4 en el caso de la altura efectiva. Las estimaciones de ambas variables con redes neuronales y diversas combinaciones de datos, arrojan resultados con R superiores a 0,85 para el caso del grado de cubierta vegetal, y 0,6 para la altura efectiva. Los parámetros multiangulares proporcionados en los productos elaborados con MISR a 1.100 m, no obtienen buenos resultados por sí mismos pero producen cierta mejora al incorporarlos a la información espectral. Los errores cuadráticos medios obtenidos son inferiores a 0,016 para la Fcc de la vegetación en tanto por uno, y 0,7 m para la altura efectiva de la misma. Regresiones geográficamente ponderadas muestran además que localmente se pueden obtener mejores resultados aún mejores, especialmente cuando hay una mayor variabilidad espacial de las variables estimadas. En resumen, la utilización de los datos proporcionados por MISR ofrece una prometedora vía de mejora de resultados en la media-baja resolución, tanto para la clasificación de imágenes como para la obtención de variables cuantitativas de la estructura de la vegetación. ABSTRACT Applications of remote sensing for monitoring what is happening on the land surface have been multiplied and refined with the launch of new sensors by different Space Agencies. The need of having up to date and spatially homogeneous data, has led to the development of new programs such as the Earth Observing System (EOS) of the National Aeronautics and Space Administration (NASA). One of the sensors incorporating the flagship of that program, the TERRA satellite, is Multi-angle Imaging Spectroradiometer (MISR), designed to capture the multi-angle information of the Earth's surface. Since the 1970s, it was known that the reflectance of various land covers and land uses varied depending on the viewing and ilumination angles, so they are anisotropic. Such variation was also related to the three dimensional structure of such covers, so that one could take advantage of such a relationship to obtain information from that structure, beyond which spectral information could provide. The MISR sensor incorporates 9 cameras at different angles to capture 9 almost simultaneous images of the same point, allowing relatively reliable estimates of the anisotropic response of the Earth's surface. Several studies have shown that we can estimate variables related to the vegetation structure with the information provided by this sensor, so this thesis has made an initial application to the Iberian Peninsula, to check their usefulness in estimating forest variables of interest. In a first step we analyzed the temporal variability that occurs in the data, due to the changes in the acquisition geometry, i.e. the relative position of sensor and light source, which in this case is the Sun. It has been found that the anisotropy is greater from late fall through early spring due to the Sun's position closer to the plane of the sensors. It was also found that the maximum and minimum values are displaced temporarily between the center and the ends. In characterizing CORINE Land Covers that has been done, one could see how the predominant form in the images with the highest sun is convex with a maximum in the camera closer to the light source. However, when the sun is much lower, the maximum is external. Moreover, the data obtained for each land cover are much more variable in summer that in November, possibly due to the proportional increase in shadow areas. To check whether the information has any effect on multi-angle imaging classification of land cover and land use, a series of classifications have been produced changing the data used, from only multispectrally, to multi-angle and multispectral. The results show that while for the most generic classifications multi-angle information is the worst, as there are extended the number of classes to obtain such information it improves the results. On the other hand, an estimate was made of quantitative variables such as canopy cover and vegetation height using information provided by MISR at different resolutions. In the valley of Alcudia (Ciudad Real), we estimated the canopy cover of trees for a pixel of 275 m by using neural networks. The results showed that using multispectral and multiangle information can improve by almost 20% the estimates that only used multispectral data. Furthermore, the relationships obtained reached an R coefficient of 0.7 with errors below 10% in canopy cover, which is much better result than the one obtained using data from the Moderate Resolution Imaging Spectroradiometer (MODIS), also onboard Terra, for the same variable. Finally we estimated the canopy cover and the effective height of the vegetation for 700,000 hectares in the province of Murcia, with a spatial resolution of 1,100 m. The results show a relationship between the spectral and the multi-angle data, and provide estimates of the canopy cover with a Spearman’s coefficient of 0.8 in the case of the vegetation canopy cover, and 0.4 in the case of the effective height. The estimates of both variables using neural networks and various combinations of data, yield results with an R coefficient greater than 0.85 for the case of the canopy cover, and 0.6 for the effective height. Multi-angle parameters provided in the products made from MISR at 1,100 m pixel size, did not produce good results from themselves but improved the results when included to the spectral information. The mean square errors were less than 0.016 for the canopy cover, and 0.7 m for the effective height. Geographically weighted regressions also showed that locally we can have even better results, especially when there is high spatial variability of estimated variables. In summary, the use of the data provided by MISR offers a promising way of improving remote sensing performance in the low-medium spatial resolution, both for image classification and for the estimation of quantitative variables of the vegetation structure.

Simulation of oil productivity and quality of N-S oriented olive hedgerow orchards in response to structure and interception of radiation

Simulations of oil yield and quality are presented for N–S oriented, hedgerow olive orchards of a range of structures (viz. canopy depth, canopy width, canopy slope and row spacing) using responses of yield and quality parameters to solar irradiance on canopy walls measured in a range of orchards, cv. Arbequina, in Spain. Results reveal that orchard yield of hedgerows of rectangular shape reaches a maximum when canopy depth equals alley width (row spacing−canopy width) and decreases at wider spacing, and/or with wider canopies, as the length of productive row decreases per unit area. Maximum yields for 4-m deep canopies were 2885 kg ha−1 at 1-m width and 5-m row spacing, 2400 kg ha−1 at 2-m width and 6-m spacing, and 2050 kg ha−1 at 3-m width and 7-m spacing. Illumination of canopies can be increased by applying slopes to form rhomboidal hedgerows. Substantial yield advantage can be achieved, especially for wide hedgerows, partly by closer row spacing that increases row length per unit area. By comparison, responses to latitude in the range 30–40◦ are small and do not warrant different row spacing. Oil quality parameters also respond to orchard structure. Responses are presented for oleic and palmitic acid, stability, and maturity index. Oleic acid content declines as alley spacing increases and is smaller, shallow than in wide, deep canopies. Palmitic acid content, stability, and maturity index increase with row alley spacing and are greater in narrow, shallow than in wide, deep canopies.

Ground cover and leaf area index relationship in a grass, legume and crucifer crop

Canopy characterization is essential for describing the interaction of a crop with its environment. The goal of this work was to determine the relationship between leaf area index (LAI) and ground cover (GC) in a grass, a legume and a crucifer crop, and to assess the feasibility of using these relationships as well as LAI-2000 readings to estimate LAI. Twelve plots were sown with either barley (Hordeum vulgare L.), vetch (Vicia sativa L.), or rape (Brassica napus L.). On 10 sampling dates the LAI (both direct and LAI-2000 estimations), fraction intercepted of photosynthetically active radiation (FIPAR) and GC were measured. Linear and quadratic models fitted to the relationship between the GC and LAI for all of the crops, but they reached a plateau in the grass when the LAI mayor que 4. Before reaching full cover, the slope of the linear relationship between both variables was within the range of 0.025 to 0.030. The LAI-2000 readings were linearly correlated with the LAI but they tended to overestimation. Corrections based on the clumping effect reduced the root mean square error of the estimated LAI from the LAI-2000 readings from 1.2 to less than 0.50 for the crucifer and the legume, but were not effective for barley.

Micro-scale habitat associations of woody plants in a neotropical cloud forest

Species?habitat associations may contribute to the maintenance of species richness in tropical forests, but previous research has been conducted almost exclusively in lowland forests and has emphasized the importance of topography and edaphic conditions. Is the distribution of woody plant species in a Peruvian cloud forest determined by microhabitat conditions? What is the role of environmental characteristics and forest structure in habitat partitioning in a tropical cloud forest? We examined species?habitat associations in three 1-ha plots using the torus-translation method. We used three different criteria to define habitats for habitat partitioning analyses, based on microtopography, forest structure and both sets of factors. The number of species associated either positively or negatively with each habitat was assessed. Habitats defined on the basis of environmental conditions and forest structure discriminated a greater number of positive and negative associations at the scale of our analyses in a tropical cloud forest. Both topographic conditions and forest structure contribute to small-scale microhabitat partitioning of woody plant species in a Peruvian tropical cloud forest. Nevertheless, canopy species were most correlated with the distribution of environmental variables, while understorey species displayed associations with forest structure.

Relationship between solar radiation on watercore on apple fruit assessed with MRI

This work is a preliminary studio of the possibility of assess a relationship between solar radiation and watercore development on apple fruit, during maturation, using a non destructive method such as Magnetic Resonance Imaging (MRI). For such purpose, several low cost solar radiation sensors were designed for the trial and placed at 2 different heights (1.5 and 2.5 m) on 6 adult ?Esperiega? apple trees, in a commercial orchard in Ademuz (Valencia). Sensors were connected along 27 days, during the end of the growth period and start of the fruit maturation process, and radiation measurements of the a-Si sensors were recorded every 1 minute. At the end of this period, fruits from the upper and the lower part of the canopy of each tree were harvested. In all, 152 apples were collected and images with MRI. A Principal Component Analysis, perfomed over the histograms of the images, as well as segmentation methods were performed on the MR images in order to find a pattern involving solar radiation and watercore incidence.

Optimización de la Exactitud Geométrica al Integrar Diversas Fuentes de Datos en Proyectos de Inventario Forestal Basados en Escaneo Laser Aerotransportado=Optimizing Geometric Accuracy in Data Source Integration for Airborne Laser Scanning-based Forest Inventory Projects

La mayoría de las aplicaciones forestales del escaneo laser aerotransportado (ALS, del inglés airborne laser scanning) requieren la integración y uso simultaneo de diversas fuentes de datos, con el propósito de conseguir diversos objetivos. Los proyectos basados en sensores remotos normalmente consisten en aumentar la escala de estudio progresivamente a lo largo de varias fases de fusión de datos: desde la información más detallada obtenida sobre un área limitada (la parcela de campo), hasta una respuesta general de la cubierta forestal detectada a distancia de forma más incierta pero cubriendo un área mucho más amplia (la extensión cubierta por el vuelo o el satélite). Todas las fuentes de datos necesitan en ultimo termino basarse en las tecnologías de sistemas de navegación global por satélite (GNSS, del inglés global navigation satellite systems), las cuales son especialmente erróneas al operar por debajo del dosel forestal. Otras etapas adicionales de procesamiento, como la ortorectificación, también pueden verse afectadas por la presencia de vegetación, deteriorando la exactitud de las coordenadas de referencia de las imágenes ópticas. Todos estos errores introducen ruido en los modelos, ya que los predictores se desplazan de la posición real donde se sitúa su variable respuesta. El grado por el que las estimaciones forestales se ven afectadas depende de la dispersión espacial de las variables involucradas, y también de la escala utilizada en cada caso. Esta tesis revisa las fuentes de error posicional que pueden afectar a los diversos datos de entrada involucrados en un proyecto de inventario forestal basado en teledetección ALS, y como las propiedades del dosel forestal en sí afecta a su magnitud, aconsejando en consecuencia métodos para su reducción. También se incluye una discusión sobre las formas más apropiadas de medir exactitud y precisión en cada caso, y como los errores de posicionamiento de hecho afectan a la calidad de las estimaciones, con vistas a una planificación eficiente de la adquisición de los datos. La optimización final en el posicionamiento GNSS y de la radiometría del sensor óptico permitió detectar la importancia de este ultimo en la predicción de la desidad relativa de un bosque monoespecífico de Pinus sylvestris L. ABSTRACT Most forestry applications of airborne laser scanning (ALS) require the integration and simultaneous use of various data sources, pursuing a variety of different objectives. Projects based on remotely-sensed data generally consist in upscaling data fusion stages: from the most detailed information obtained for a limited area (field plot) to a more uncertain forest response sensed over a larger extent (airborne and satellite swath). All data sources ultimately rely on global navigation satellite systems (GNSS), which are especially error-prone when operating under forest canopies. Other additional processing stages, such as orthorectification, may as well be affected by vegetation, hence deteriorating the accuracy of optical imagery’s reference coordinates. These errors introduce noise to the models, as predictors displace from their corresponding response. The degree to which forest estimations are affected depends on the spatial dispersion of the variables involved and the scale used. This thesis reviews the sources of positioning errors which may affect the different inputs involved in an ALS-assisted forest inventory project, and how the properties of the forest canopy itself affects their magnitude, advising on methods for diminishing them. It is also discussed how accuracy should be assessed, and how positioning errors actually affect forest estimation, toward a cost-efficient planning for data acquisition. The final optimization in positioning the GNSS and optical image allowed to detect the importance of the latter in predicting relative density in a monospecific Pinus sylvestris L. forest.

Flujos de gases de efecto invernadero en suelos forestales y agroforestales del centro de la Península Ibérica=Greenhouse gas flows in forest and agroforestry soils in the center of the Iberian Peninsula

Debido a la complejidad de los procesos que controlan el intercambio de gases de carbono (C) y nitrógeno (N) entre el suelo y la atmósfera, en los sistemas forestales y agroforestales, son comprensibles las incógnitas existentes respecto a la estimación de los flujos de los gases de efecto invernadero (GEI) y la capacidad como reservorios de carbono de los suelos, bajo diferentes formas de uso y regímenes de alteración a escala regional y global. Esta escasez de información justifica la necesidad de caracterizar la dinámica de intercambio de GEI en los ecosistemas Mediterráneos, en especial en el contexto actual de cambio climático, y el incremento asociado de temperatura y periodos de sequía, alteración de los patrones de precipitación, y el riesgo de incendios forestales; cuyas consecuencias afectarán tanto a los compartimentos de C y de N del suelo como a la capacidad de secuestro de C de estos ecosistemas. Dentro de este contexto se enmarca la presente tesis doctoral cuyo objetivo ha sido cuantificar y caracterizar los flujos de dióxido de carbono (CO2), de oxido nitroso (N2O) y de metano (CH4), junto con los stocks de C y N, en suelos forestales de Quercus ilex, Quercus pyrenaica y Pinus sylvestris afectados por incendios forestales; así como el estudiar el efecto de la gestión y la cubierta arbórea en la respiración del suelo y los stocks de C y N en una dehesa situada en el centro de la Península Ibérica. De manera que los flujos de CO2, N2O y CH4; y los parámetros físico-químicos y biológicos del suelo fueron estudiados en los diferentes tratamientos y ecosistemas a lo largo del trabajo que se presenta. Los resultados obtenidos muestran la existencia de variaciones temporales y espaciales de la respiración del suelo dentro de una escala geográfica pequeña, controladas principalmente por la temperatura y la humedad del suelo; y por los contenidos de C y N del suelo en un bosque de Pinus sylvestris en la vertiente norte de la Sierra de Guadarrama , en España. El análisis de los efectos de los incendios forestales a largo plazo (6-8 años) revela que las pérdidas anuales de C a través de la respiración del suelo en las zonas quemadas de Quercus ilex, Quercus pyrenaica y Pinus sylvestris fueron 450 gCm-2yr-1, 790 gCm-2yr-1 y 1220 gCm-2yr-1, respectivamente; lo que representa una reducción del 43%, 22% y 11% en comparación con las zonas no quemadas de dichas especies, debido a la destrucción de la masa arbórea. El efecto del fuego también alteró los flujos N2O y CH4 del suelo, de una forma diferente en los distintos ecosistemas y estacionalidades estudiadas. De tal modo, que los suelos quemados mostraron una mayor oxidación del CH4 en las masas de Q. ilex, y una menor oxidación en las de P. sylvestris; además de una disminución de los flujos de N2O en Q. pyrenaica. Los incendios también afectaron los parámetros microclimáticos de los suelos forestales, observándose un incremento de la temperatura del suelo y una disminución de la humedad en los emplazamientos quemados que en los no quemados. Los cationes intercambiables, el pH, el cociente C/N, el contenido en raicillas y la biomasa microbiana también disminuyeron en las zonas quemadas. Aunque el C orgánico del suelo no se alteró de manera significativa, si lo hizo la calidad de la materia orgánica, disminuyendo el carbono lábil y aumentando las formas recalcitrantes lo que se tradujo en menor sensibilidad de la respiración del suelo a la temperatura (valores de Q10) en las zonas quemadas. Los resultados del estudio realizado en la Dehesa muestran que las actividades silvopastorales estudiadas afectaron levemente y de forma no constante a la respiración del suelo y las condiciones microclimáticas del suelo. Se observó una reducción 12% de la respiración del suelo por efecto del pastoreo no intensivo. Sin embargo, se observaron incrementos de 3Mg/ha en los stocks de C y de 0.3 Mg/ha en los stocks de N en los suelos pastoreados en comparación con los no pastoreados. Aunque, no se observó un claro efecto de la labranza sobre la respiración del suelo en nuestro experimento, sin embargo si se observó una disminución de 3.5 Mg/ha en las reservas de C y de 0.3 Mg/ ha en las de N en los suelos labrados comparados con los no labrados. La copa del arbolado influyó de forma positiva tanto en la respiración del suelo, como en los stocks de C y N de los suelos. La humedad del suelo jugó un papel relevante en la sensibilidad de la respiración a la temperatura del suelo. Nuestros resultados ponen de manifiesto la sensibilidad de la respiración del suelo a cambios en la humedad y los parámetros edáficos, y sugieren que la aplicación de modelos estándar para estimar la respiración del suelo en áreas geográficas pequeñas puede no ser adecuada a menos que otros factores sean considerados en combinación con la temperatura del suelo. Además, las diferentes respuestas de los flujos de gases de efecto invernadero a los cambios, años después de la ocurrencia de incendios forestales, destaca la necesidad de incluir estos cambios en las futuras investigaciones de la dinámica del carbono en los ecosistemas mediterráneos. Por otra parte, las respuestas divergentes en los valores de respiración del suelo y en los contenidos de C y N del suelo observados en la dehesa, además de la contribución de la copa de los árboles en los nutrientes del suelo ilustran la importancia de mantener la gestión tradicional aplicada en beneficio de la capacidad de almacenar C en la dehesa estudiada. La información obtenida en este trabajo pretende contribuir a la mejora del conocimiento de la dinámica y el balance de C en los sistemas mediterráneos, además de ayudar a predecir el impacto del cambio climático en el intercambio de C entre los ecosistemas forestales y agroforestales y la atmósfera. ABSTRACT Due to the complexity of the processes that control the exchange of carbon (C) and nitrogen (N) gasses between soils and the atmosphere in forest and agroforestry ecosystems, understandable uncertainties exist as regards the estimation of greenhouse gas (GHG) fluxes and the soil sink capacity at regional and global scale under different forms of land use and disturbance regimes. These uncertainties justify the need to characterize the exchange dynamics of GHG between the atmosphere and soils in Mediterranean terrestrial ecosystems, particularly in the current context of climate change and the associated increase in temperature, drought periods, heavy rainfall events, and increased risk of wildfires, which affect not only the C and N pools but also the soil C sink capacity of these ecosystems. Within this context, the aims of the present thesis were, firstly, to quantify and characterize the fluxes of carbon dioxide (CO2), nitrous oxide (N2O) and methane (CH4) as well as the C and N stocks in Quercus ilex, Quercus pyrenaica and Pinus sylvestris stands affected by wildfires, and secondly, to study the effects of Quercus ilex canopy and management on both soil respiration and C and N pools in dehesa systems in the center of Iberian Peninsula. Soil CO2, N2O and CH4 fluxes, and soil physical-chemical and biological parameters were studied under the different treatments and ecosystems considered in this study. The results showed seasonal and spatial variations in soil respiration within small geographic areas, mainly controlled by soil temperature and moisture in addition to soil carbon and nitrogen stocks in mixed pine–oak forest ecosystems on the north facing slopes of the Sierra de Guadarrama in Spain. The analysis of long term effects of wildfires (6–8 years) revealed that annual carbon losses through soil respiration from burned sites in Quercus ilex, Quercus pyrenaica and Pinus sylvestris stands were 450 gCm-2yr-1, 790 gCm-2yr-1 and 1220 gCm-2yr-1, respectively; with burned sites emitting 43%, 22% and 11% less in burned as opposed to non-burned sites due the loss of trees. Fire may alter both N2O and CH4 fluxes although the magnitude of such variation depends on the site, soil characteristics and seasonal climatic conditions. The burned sites showed higher CH4 oxidation in Q.ilex stands, and lower oxidation rates in P. sylvestris stands. A reduction in N2O fluxes in Q. pyrenaica stands was detected at burned sites along with changes in soil microclimate; higher soil temperature and lower soil moisture content. Exchangeable cations, the C/N ratio, pH, fine root and microbial biomass were also found to decrease at burned sites. Although the soil organic carbon was not significantly altered, the quality of the organic matter changed, displaying a decrease in labile carbon and a relative increase in refractory forms, leading to lower sensitivity of soil respiration to temperature (Q10 values) at burned sites. The results from the dehesa study show that light grazing and superficial tilling practices used in the studied dehesa system in Spain had a slight but non-consistent impact on soil respiration and soil microclimate over the study period. The reduction in soil respiration in the dehesa system due to the effects of grazing was around 12 %. However, increments of 3Mg/ha in C stocks and 0.3 Mg/ha in N stocks in grazed soils were observed. Although no clear effect of tilling on soil respiration was found, a decrease of 3.5 Mg/ha in C stocks and 0.3 Mg/ha in N stocks was detected for tilled soils. The presence of a tree canopy induced increases in soil respiration, soil C and N stocks, while soil moisture was found to play an important role in soil respiration temperature response. Our results suggest that the use of standard models to estimate soil respiration in small geographical areas may not be adequate unless other factors are considered in addition to soil temperature. Furthermore, the different responses of GHG flux to climatic shifts, many years after the occurrence of wildfire, highlight the need to include these shifts in C dynamics in future research undertaken in Mediterranean ecosystems. Furthermore, divergent responses in soil respiration and soil C and N stocks to grazing or tilling practices in Dehesa systems, and the influence of tree canopy on soil respiration and soil nutrient content, illustrate the importance of maintaining beneficial management practices. Moreover, the carbon sequestration capacity of the Dehesa system studied may be enhanced through improvements in the management applied. It is hoped that the information obtained through this research will contribute towards improving our understanding of the dynamics and balance of C in Mediterranean systems, and help predict the impact of climate change on the exchange of C between forest and agroforestry ecosystems and the atmosphere.

Simplifying data acquisition in plant canopies- Measurements of leaf angles with a cell phone

1. Canopies are complex multilayered structures comprising individual plant crowns exposing a multifaceted surface area to sunlight. Foliage arrangement and properties are the main mediators of canopy functions. The leaves act as light traps whose exposure to sunlight varies with time of the day, date and latitude in a trade-off between photosynthetic light harvesting and excessive or photoinhibitory light avoidance. To date, ecological research based upon leaf sampling has been limited by the available echnology, with which data acquisition becomes labour intensive and time-consuming, given the verwhelming number of leaves involved. 2. In the present study, our goal involved developing a tool capable of easuring a sufficient number of leaves to enable analysis of leaf populations, tree crowns and canopies.We specifically tested whether a cell phone working as a 3Dpointer could yield reliable, repeatable and valid leaf anglemeasurements with a simple gesture. We evaluated the accuracy of this method under controlled conditions, using a 3D digitizer, and we compared performance in the field with the methods commonly used. We presented an equation to estimate the potential proportion of the leaf exposed to direct sunlight (SAL) at any given time and compared the results with those obtained bymeans of a graphicalmethod. 3. We found a strong and highly significant correlation between the graphical methods and the equation presented. The calibration process showed a strong correlation between the results derived from the two methods with amean relative difference below 10%. Themean relative difference in calculation of instantaneous exposure was below 5%. Our device performed equally well in diverse locations, in which we characterized over 700 leaves in a single day. 4. The newmethod, involving the use of a cell phone, ismuchmore effective than the traditionalmethods or digitizers when the goal is to scale up from leaf position to performance of leaf populations, tree crowns or canopies. Our methodology constitutes an affordable and valuable tool within which to frame a wide range of ecological hypotheses and to support canopy modelling approaches.

Low-cost irradiance sensors for irradiation assessments inside tree canopies.

The solar irradiation that a crop receives is directly related to the physical and biological processes that affect the crop. However, the assessment of solar irradiation poses certain problems when it must be measured through fruit inside the canopy of a tree. In such cases, it is necessary to check many test points, which usually requires an expensive data acquisition system. The use of conventional irradiance sensors increases the cost of the experiment, making them unsuitable. Nevertheless, it is still possible to perform a precise irradiance test with a reduced price by using low-cost sensors based on the photovoltaic effect. The aim of this work is to develop a low-cost sensor that permits the measurement of the irradiance inside the tree canopy. Two different technologies of solar cells were analyzed for their use in the measurement of solar irradiation levels inside tree canopies. Two data acquisition system setups were also tested and compared. Experiments were performed in Ademuz (Valencia, Spain) in September 2011 and September 2012 to check the validity of low-cost sensors based on solar cells and their associated data acquisition systems. The observed difference between solar irradiation at high and low positions was of 18.5% ± 2.58% at a 95% confidence interval. Large differences were observed between the operations of the two tested sensors. In the case of a-Si cells based mini-modules, an effect of partial shadowing was detected due to the larger size of the devices, the use of individual c-Si cells is recommended over a-Si cells based mini-modules.

Remote sensing of LAI, chlorophyll and leaf nitrogen pools of crop- and grasslands in five European landscapes.

Leaf nitrogen and leaf surface area influence the exchange of gases between terrestrial ecosystems and the atmosphere, and play a significant role in the global cycles of carbon, nitrogen and water. The purpose of this study is to use field-based and satellite remote-sensing-based methods to assess leaf nitrogen pools in five diverse European agricultural landscapes located in Denmark, Scotland (United Kingdom), Poland, the Netherlands and Italy. REGFLEC (REGularized canopy reFLECtance) is an advanced image-based inverse canopy radiative transfer modelling system which has shown proficiency for regional mapping of leaf area index (LAI) and leaf chlorophyll (CHLl) using remote sensing data. In this study, high spatial resolution (10–20 m) remote sensing images acquired from the multispectral sensors aboard the SPOT (Satellite For Observation of Earth) satellites were used to assess the capability of REGFLEC for mapping spatial variations in LAI, CHLland the relation to leaf nitrogen (Nl) data in five diverse European agricultural landscapes. REGFLEC is based on physical laws and includes an automatic model parameterization scheme which makes the tool independent of field data for model calibration. In this study, REGFLEC performance was evaluated using LAI measurements and non-destructive measurements (using a SPAD meter) of leaf-scale CHLl and Nl concentrations in 93 fields representing crop- and grasslands of the five landscapes. Furthermore, empirical relationships between field measurements (LAI, CHLl and Nl and five spectral vegetation indices (the Normalized Difference Vegetation Index, the Simple Ratio, the Enhanced Vegetation Index-2, the Green Normalized Difference Vegetation Index, and the green chlorophyll index) were used to assess field data coherence and to serve as a comparison basis for assessing REGFLEC model performance. The field measurements showed strong vertical CHLl gradient profiles in 26% of fields which affected REGFLEC performance as well as the relationships between spectral vegetation indices (SVIs) and field measurements. When the range of surface types increased, the REGFLEC results were in better agreement with field data than the empirical SVI regression models. Selecting only homogeneous canopies with uniform CHLl distributions as reference data for evaluation, REGFLEC was able to explain 69% of LAI observations (rmse = 0.76), 46% of measured canopy chlorophyll contents (rmse = 719 mg m−2) and 51% of measured canopy nitrogen contents (rmse = 2.7 g m−2). Better results were obtained for individual landscapes, except for Italy, where REGFLEC performed poorly due to a lack of dense vegetation canopies at the time of satellite recording. Presence of vegetation is needed to parameterize the REGFLEC model. Combining REGFLEC- and SVI-based model results to minimize errors for a "snap-shot" assessment of total leaf nitrogen pools in the five landscapes, results varied from 0.6 to 4.0 t km−2. Differences in leaf nitrogen pools between landscapes are attributed to seasonal variations, extents of agricultural area, species variations, and spatial variations in nutrient availability. In order to facilitate a substantial assessment of variations in Nl pools and their relation to landscape based nitrogen and carbon cycling processes, time series of satellite data are needed. The upcoming Sentinel-2 satellite mission will provide new multiple narrowband data opportunities at high spatio-temporal resolution which are expected to further improve remote sensing capabilities for mapping LAI, CHLl and Nl.

Evaluation of cover crops based on characteristics of nitrogen use, and aerial and root growth

La caracterización de los cultivos cubierta (cover crops) puede permitir comparar la idoneidad de diferentes especies para proporcionar servicios ecológicos como el control de la erosión, el reciclado de nutrientes o la producción de forrajes. En este trabajo se estudiaron bajo condiciones de campo diferentes técnicas para caracterizar el dosel vegetal con objeto de establecer una metodología para medir y comparar las arquitecturas de los cultivos cubierta más comunes. Se estableció un ensayo de campo en Madrid (España central) para determinar la relación entre el índice de área foliar (LAI) y la cobertura del suelo (GC) para un cultivo de gramínea, uno de leguminosa y uno de crucífera. Para ello se sembraron doce parcelas con cebada (Hordeum vulgare L.), veza (Vicia sativa L.), y colza (Brassica napus L.). En 10 fechas de muestreo se midieron el LAI (con estimaciones directas y del LAI-2000), la fracción interceptada de la radiación fotosintéticamente activa (FIPAR) y la GC. Un experimento de campo de dos años (Octubre-Abril) se estableció en la misma localización para evaluar diferentes especies (Hordeum vulgare L., Secale cereale L., x Triticosecale Whim, Sinapis alba L., Vicia sativa L.) y cultivares (20) en relación con su idoneidad para ser usadas como cultivos cubierta. La GC se monitorizó mediante análisis de imágenes digitales con 21 y 22 muestreos, y la biomasa se midió 8 y 10 veces, respectivamente para cada año. Un modelo de Gompertz caracterizó la cobertura del suelo hasta el decaimiento observado tras las heladas, mientras que la biomasa se ajustó a ecuaciones de Gompertz, logísticas y lineales-exponenciales. Al final del experimento se determinaron el C, el N y el contenido en fibra (neutrodetergente, ácidodetergente y lignina), así como el N fijado por las leguminosas. Se aplicó el análisis de decisión multicriterio (MCDA) con objeto de obtener un ranking de especies y cultivares de acuerdo con su idoneidad para actuar como cultivos cubierta en cuatro modalidades diferentes: cultivo de cobertura, cultivo captura, abono verde y forraje. Las asociaciones de cultivos leguminosas con no leguminosas pueden afectar al crecimiento radicular y a la absorción de N de ambos componentes de la mezcla. El conocimiento de cómo los sistemas radiculares específicos afectan al crecimiento individual de las especies es útil para entender las interacciones en las asociaciones, así como para planificar estrategias de cultivos cubierta. En un tercer ensayo se combinaron estudios en rhizotrones con extracción de raíces e identificación de especies por microscopía, así como con estudios de crecimiento, absorción de N y 15N en capas profundas del suelo. Las interacciones entre raíces en su crecimiento y en el aprovisionamiento de N se estudiaron para dos de los cultivares mejor valorados en el estudio previo: uno de cebada (Hordeum vulgare L. cv. Hispanic) y otro de veza (Vicia sativa L. cv. Aitana). Se añadió N en dosis de 0 (N0), 50 (N1) y 150 (N2) kg N ha-1. Como resultados del primer estudio, se ajustaron correctamente modelos lineales y cuadráticos a la relación entre la GC y el LAI para todos los cultivos, pero en la gramínea alcanzaron una meseta para un LAI>4. Antes de alcanzar la cobertura total, la pendiente de la relación lineal entre ambas variables se situó en un rango entre 0.025 y 0.030. Las lecturas del LAI-2000 estuvieron correlacionadas linealmente con el LAI, aunque con tendencia a la sobreestimación. Las correcciones basadas en el efecto de aglutinación redujeron el error cuadrático medio del LAI estimado por el LAI-2000 desde 1.2 hasta 0.5 para la crucífera y la leguminosa, no siendo efectivas para la cebada. Esto determinó que para los siguientes estudios se midieran únicamente la GC y la biomasa. En el segundo experimento, las gramíneas alcanzaron la mayor cobertura del suelo (83-99%) y la mayor biomasa (1226-1928 g m-2) al final del mismo. Con la mayor relación C/N (27-39) y contenido en fibra digestible (53-60%) y la menor calidad de residuo (~68%). La mostaza presentó elevadas GC, biomasa y absorción de N en el año más templado en similitud con las gramíneas, aunque escasa calidad como forraje en ambos años. La veza presentó la menor absorción de N (2.4-0.7 g N m-2) debido a la fijación de N (9.8-1.6 g N m-2) y escasa acumulación de N. El tiempo térmico hasta alcanzar el 30% de GC constituyó un buen indicador de especies de rápida cubrición. La cuantificación de las variables permitió hallar variabilidad entre las especies y proporcionó información para posteriores decisiones sobre la selección y manejo de los cultivos cubierta. La agregación de dichas variables a través de funciones de utilidad permitió confeccionar rankings de especies y cultivares para cada uso. Las gramíneas fueron las más indicadas para los usos de cultivo de cobertura, cultivo captura y forraje, mientras que las vezas fueron las mejor como abono verde. La mostaza alcanzó altos valores como cultivo de cobertura y captura en el primer año, pero el segundo decayó debido a su pobre actuación en los inviernos fríos. Hispanic fue el mejor cultivar de cebada como cultivo de cobertura y captura, mientras que Albacete como forraje. El triticale Titania alcanzó la posición más alta como cultiva de cobertura, captura y forraje. Las vezas Aitana y BGE014897 mostraron buenas aptitudes como abono verde y cultivo captura. El MCDA permitió la comparación entre especies y cultivares proporcionando información relevante para la selección y manejo de cultivos cubierta. En el estudio en rhizotrones tanto la mezcla de especies como la cebada alcanzaron mayor intensidad de raíces (RI) y profundidad (RD) que la veza, con valores alrededor de 150 cruces m-1 y 1.4 m respectivamente, comparados con 50 cruces m-1 y 0.9 m para la veza. En las capas más profundas del suelo, la asociación de cultivos mostró valores de RI ligeramente mayores que la cebada en monocultivo. La cebada y la asociación obtuvieron mayores valores de densidad de raíces (RLD) (200-600 m m-3) que la veza (25-130) entre 0.8 y 1.2 m de profundidad. Los niveles de N no mostraron efectos claros en RI, RD ó RLD, sin embargo, el incremento de N favoreció la proliferación de raíces de veza en la asociación en capas profundas del suelo, con un ratio cebada/veza situado entre 25 a N0 y 5 a N2. La absorción de N de la cebada se incrementó en la asociación a expensas de la veza (de ~100 a 200 mg planta-1). Las raíces de cebada en la asociación absorbieron también más nitrógeno marcado de las capas profundas del suelo (0.6 mg 15N planta-1) que en el monocultivo (0.3 mg 15N planta-1). ABSTRACT Cover crop characterization may allow comparing the suitability of different species to provide ecological services such as erosion control, nutrient recycling or fodder production. Different techniques to characterize plant canopy were studied under field conditions in order to establish a methodology for measuring and comparing cover crops canopies. A field trial was established in Madrid (central Spain) to determine the relationship between leaf area index (LAI) and ground cover (GC) in a grass, a legume and a crucifer crop. Twelve plots were sown with either barley (Hordeum vulgare L.), vetch (Vicia sativa L.), or rape (Brassica napus L.). On 10 sampling dates the LAI (both direct and LAI-2000 estimations), fraction intercepted of photosynthetically active radiation (FIPAR) and GC were measured. A two-year field experiment (October-April) was established in the same location to evaluate different species (Hordeum vulgare L., Secale cereale L., x Triticosecale Whim, Sinapis alba L., Vicia sativa L.) and cultivars (20) according to their suitability to be used as cover crops. GC was monitored through digital image analysis with 21 and 22 samples, and biomass measured 8 and 10 times, respectively for each season. A Gompertz model characterized ground cover until the decay observed after frosts, while biomass was fitted to Gompertz, logistic and linear-exponential equations. At the end of the experiment C, N, and fiber (neutral detergent, acid and lignin) contents, and the N fixed by the legumes were determined. Multicriteria decision analysis (MCDA) was applied in order to rank the species and cultivars according to their suitability to perform as cover crops in four different modalities: cover crop, catch crop, green manure and fodder. Intercropping legumes and non-legumes may affect the root growth and N uptake of both components in the mixture. The knowledge of how specific root systems affect the growth of the individual species is useful for understanding the interactions in intercrops as well as for planning cover cropping strategies. In a third trial rhizotron studies were combined with root extraction and species identification by microscopy and with studies of growth, N uptake and 15N uptake from deeper soil layers. The root interactions of root growth and N foraging were studied for two of the best ranked cultivars in the previous study: a barley (Hordeum vulgare L. cv. Hispanic) and a vetch (Vicia sativa L. cv. Aitana). N was added at 0 (N0), 50 (N1) and 150 (N2) kg N ha-1. As a result, linear and quadratic models fitted to the relationship between the GC and LAI for all of the crops, but they reached a plateau in the grass when the LAI > 4. Before reaching full cover, the slope of the linear relationship between both variables was within the range of 0.025 to 0.030. The LAI-2000 readings were linearly correlated with the LAI but they tended to overestimation. Corrections based on the clumping effect reduced the root mean square error of the estimated LAI from the LAI-2000 readings from 1.2 to less than 0.50 for the crucifer and the legume, but were not effective for barley. This determined that in the following studies only the GC and biomass were measured. In the second experiment, the grasses reached the highest ground cover (83- 99%) and biomass (1226-1928 g/m2) at the end of the experiment. The grasses had the highest C/N ratio (27-39) and dietary fiber (53-60%) and the lowest residue quality (~68%). The mustard presented high GC, biomass and N uptake in the warmer year with similarity to grasses, but low fodder capability in both years. The vetch presented the lowest N uptake (2.4-0.7 g N/m2) due to N fixation (9.8-1.6 g N/m2) and low biomass accumulation. The thermal time until reaching 30% ground cover was a good indicator of early coverage species. Variable quantification allowed finding variability among the species and provided information for further decisions involving cover crops selection and management. Aggregation of these variables through utility functions allowed ranking species and cultivars for each usage. Grasses were the most suitable for the cover crop, catch crop and fodder uses, while the vetches were the best as green manures. The mustard attained high ranks as cover and catch crop the first season, but the second decayed due to low performance in cold winters. Hispanic was the most suitable barley cultivar as cover and catch crop, and Albacete as fodder. The triticale Titania attained the highest rank as cover and catch crop and fodder. Vetches Aitana and BGE014897 showed good aptitudes as green manures and catch crops. MCDA allowed comparison among species and cultivars and might provide relevant information for cover crops selection and management. In the rhizotron study the intercrop and the barley attained slightly higher root intensity (RI) and root depth (RD) than the vetch, with values around 150 crosses m-1 and 1.4 m respectively, compared to 50 crosses m-1 and 0.9 m for the vetch. At deep soil layers, intercropping showed slightly larger RI values compared to the sole cropped barley. The barley and the intercropping had larger root length density (RLD) values (200-600 m m-3) than the vetch (25-130) at 0.8-1.2 m depth. The topsoil N supply did not show a clear effect on the RI, RD or RLD; however increasing topsoil N favored the proliferation of vetch roots in the intercropping at deep soil layers, with the barley/vetch root ratio ranging from 25 at N0 to 5 at N2. The N uptake of the barley was enhanced in the intercropping at the expense of the vetch (from ~100 mg plant-1 to 200). The intercropped barley roots took up more labeled nitrogen (0.6 mg 15N plant-1) than the sole-cropped barley roots (0.3 mg 15N plant-1) from deep layers.

Application of NMR techniques for the evaluation of fruit internal quality and physical properties of food

Los alimentos son sistemas complejos, formados por diversas estructuras a diferentes escalas: macroscópica y microscópica. Muchas propiedades de los alimentos, que son importantes para su procesamiento, calidad y tratamiento postcosecha, están relacionados con su microestructura. La presente tesis doctoral propone una metodología completa para la determinación de la estructura de alimentos desde un punto de vista multi-escala, basándose en métodos de Resonancia Magnética Nuclear (NMR). Las técnicas de NMR son no invasivas y no destructivas y permiten el estudio tanto de macro- como de microestructura. Se han utilizado distintos procedimientos de NMR dependiendo del nivel que se desea estudiar. Para el nivel macroestructural, la Imagen de Resonancia Magnética (MRI) ha resultado ser muy útil para la caracterización de alimentos. Para el estudio microestructural, la MRI requiere altos tiempos de adquisición, lo que hace muy difícil la transferencia de esta técnica a aplicaciones en industria. Por tanto, la optimización de procedimientos de NMR basados en secuencias relaxometría 2D T1/T2 ha resultado ser una estrategia primordial en esta tesis. Estos protocolos de NMR se han implementado satisfactoriamente por primera vez en alto campo magnético. Se ha caracterizado la microestructura de productos alimentarios enteros por primera vez utilizando este tipo de protocolos. Como muestras, se han utilizado dos tipos de productos: modelos de alimentos y alimentos reales (manzanas). Además, como primer paso para su posterior implementación en la industria agroalimentaria, se ha mejorado una línea transportadora, especialmente diseñada para trabajar bajo condiciones de NMR en trabajos anteriores del grupo LPF-TAGRALIA. Se han estudiado y seleccionado las secuencias más rápidas y óptimas para la detección de dos tipos de desórdenes internos en manzanas: vitrescencia y roturas internas. La corrección de las imágenes en movimiento se realiza en tiempo real. Asimismo, se han utilizado protocolos de visión artificial para la clasificación automática de manzanas potencialmente afectadas por vitrescencia. El presente documento está dividido en diferentes capítulos: el Capítulo 2 explica los antecedentes de la presente tesis y el marco del proyecto en el que se ha desarrollado. El Capítulo 3 recoge el estado del arte. El Capítulo 4 establece los objetivos de esta tesis doctoral. Los resultados se dividen en cinco sub-secciones (dentro del Capítulo 5) que corresponden con los trabajos publicados bien en revistas revisadas por pares, bien en congresos internacionales o bien como capítulos de libros revisados por pares. La Sección 5.1. es un estudio del desarrollo de la vitrescencia en manzanas mediante MRI y lo relaciona con la posición de la fruta dentro de la copa del árbol. La Sección 5.2 presenta un trabajo sobre macro- y microestructura en modelos de alimentos. La Sección 5.3 es un artículo en revisión en una revista revisada por pares, en el que se hace un estudio microestrcutural no destructivo mediante relaxometría 2D T1/T2. la Sección 5.4, hace una comparación entre manzanas afectadas por vitrescencia mediante dos técnicas: tomografía de rayos X e MRI, en manzana. Por último, en la Sección 5.5 se muestra un trabajo en el que se hace un estudio de secuencias de MRI en línea para la evaluación de calidad interna en manzanas. Los siguientes capítulos ofrecen una discusión y conclusiones (Capítulo 6 y 7 respectivamente) de todos los capítulos de esta tesis doctoral. Finalmente, se han añadido tres apéndices: el primero con una introducción de los principios básicos de resonancia magnética nuclear (NMR) y en los otros dos, se presentan sendos estudios sobre el efecto de las fibras en la rehidratación de cereales de desayuno extrusionados, mediante diversas técnicas. Ambos trabajos se presentaron en un congreso internacional. Los resultados más relevantes de la presente tesis doctoral, se pueden dividir en tres grandes bloques: resultados sobre macroestructura, resultados sobre microestructura y resultados sobre MRI en línea. Resultados sobre macroestructura: - La imagen de resonancia magnética (MRI) se aplicó satisfactoriamente para la caracterización de macroestructura. En particular, la reconstrucción 3D de imágenes de resonancia magnética permitió identificar y caracterizar dos tipos distintos de vitrescencia en manzanas: central y radial, que se caracterizan por el porcentaje de daño y la conectividad (número de Euler). - La MRI proveía un mejor contraste para manzanas afectadas por vitrescencia que las imágenes de tomografía de rayos X (X-Ray CT), como se pudo verificar en muestras idénticas de manzana. Además, el tiempo de adquisición de la tomografía de rayos X fue alrededor de 12 veces mayor (25 minutos) que la adquisición de las imágenes de resonancia magnética (2 minutos 2 segundos). Resultados sobre microestructura: - Para el estudio de microestructura (nivel subcelular) se utilizaron con éxito secuencias de relaxometría 2D T1/T2. Estas secuencias se usaron por primera vez en alto campo y sobre piezas de alimento completo, convirtiéndose en una forma no destructiva de llevar a cabo estudios de microestructura. - El uso de MRI junto con relaxometría 2D T1/T2 permite realizar estudios multiescala en alimentos de forma no destructiva. Resultados sobre MRI en línea: - El uso de imagen de resonancia magnética en línea fue factible para la identificación de dos tipos de desórdenes internos en manzanas: vitrescencia y podredumbre interna. Las secuencias de imagen tipo FLASH resultaron adecuadas para la identificación en línea de vitrescencia en manzanas. Se realizó sin selección de corte, debido a que la vitrescencia puede desarrollarse en cualquier punto del volumen de la manzana. Se consiguió reducir el tiempo de adquisición, de modo que se llegaron a adquirir 1.3 frutos por segundos (758 ms por fruto). Las secuencias de imagen tipo UFLARE fueron adecuadas para la detección en línea de la podredumbre interna en manzanas. En este caso, se utilizó selección de corte, ya que se trata de un desorden que se suele localizar en la parte central del volumen de la manzana. Se consiguió reducir el tiempo de adquisicón hasta 0.67 frutos por segundo (1475 ms por fruto). En ambos casos (FLASH y UFLARE) fueron necesarios algoritmos para la corrección del movimiento de las imágenes en tiempo real. ABSTRACT Food is a complex system formed by several structures at different scales: macroscopic and microscopic. Many properties of foods that are relevant to process engineering or quality and postharvest treatments are related to their microstructure. This Ph.D Thesis proposes a complete methodology for food structure determination, in a multiscale way, based on the Nuclear Magnetic Resonance (NMR) phenomenon since NMR techniques are non-invasive and non-destructive, and allow both, macro- and micro-structure study. Different NMR procedures are used depending on the structure level under study. For the macrostructure level, Magnetic Resonance Imaging (MRI) revealed its usefulness for food characterization. For microstructure insight, MRI required high acquisition times, which is a hindrance for transference to industry applications. Therefore, optimization of NMR procedures based on T1/T2 relaxometry sequences was a key strategy in this Thesis. These NMR relaxometry protocols, are successfully implemented in high magnetic field. Microstructure of entire food products have been characterized for the first time using these protocols. Two different types of food products have been studied: food models and actual food (apples). Furthermore, as a first step for the food industry implementation, a grading line system, specially designed for working under NMR conditions in previous works of the LPF-TAGRALIA group, is improved. The study and selection of the most suitable rapid sequence to detect two different types of disorders in apples (watercore and internal breakdown) is performed and the real time image motion correction is applied. In addition, artificial vision protocols for the automatic classification of apples potentially affected by watercore are applied. This document is divided into seven different chapters: Chapter 2 explains the thesis background and the framework of the project in which it has been worked. Chapter 3 comprises the state of the art. Chapter 4 establishes de objectives of this Ph.D thesis. The results are divided into five different sections (in Chapter 5) that correspond to published peered reviewed works. Section 5.1 assesses the watercore development in apples with MRI and studies the effect of fruit location in the canopy. Section 5.2 is an MRI and 2D relaxometry study for macro- and microstructure assessment in food models. Section 5.3 is a non-destructive microstructural study using 2D T1/T2 relaxometry on watercore affected apples. Section 5.4 makes a comparison of X-ray CT and MRI on watercore disorder of different apple cultivars. Section 5.5, that is a study of online MRI sequences for the evaluation of apple internal quality. The subsequent chapters offer a general discussion and conclusions (Chapter 6 and Chapter 7 respectively) of all the works performed in the frame of this Ph.D thesis (two peer reviewed journals, one book chapter and one international congress).Finally, three appendices are included in which an introduction to NMR principles is offered and two published proceedings regarding the effect of fiber on the rehydration of extruded breakfast cereal are displayed. The most relevant results can be summarized into three sections: results on macrostructure, results on microstructure and results on on-line MRI. Results on macrostructure: - MRI was successfully used for macrostructure characterization. Indeed, 3D reconstruction of MRI in apples allows to identify two different types of watercore (radial and block), which are characterized by the percentage of damage and the connectivity (Euler number). - MRI provides better contrast for watercore than X-Ray CT as verified on identical samples. Furthermore, X-Ray CT images acquisition time was around 12 times higher (25 minutes) than MRI acquisition time (2 minutes 2 seconds). Results on microstructure: - 2D T1/T2 relaxometry were successfully applied for microstructure (subcellular level) characterization. 2D T1/T2 relaxometry sequences have been applied for the first time on high field for entire food pieces, being a non-destructive way to achieve microstructure study. - The use of MRI together with 2D T1/T2 relaxometry sequences allows a non-destructive multiscale study of food. Results on on-line MRI: - The use of on-line MRI was successful for the identification of two different internal disorders in apples: watercore and internal breakdown. FLASH imaging was a suitable technique for the on-line detection of watercore disorder in apples, with no slice selection, since watercore is a physiological disorder that may be developed anywhere in the apple volume. 1.3 fruits were imaged per second (768 ms per fruit). UFLARE imaging is a suitable sequence for the on-line detection of internal breakdown disorder in apples. Slice selection was used, as internal breakdown is usually located in the central slice of the apple volume. 0.67 fruits were imaged per second (1475 ms per fruit). In both cases (FLASH and UFLARE) motion correction was performed in real time, during the acquisition of the images.

Conductancia del mesófilo en especies forestales: implicaciones en la respuesta fotosintética a la disponibilidad de agua y luz

La fotosíntesis es el proceso biológico que permite la producción primaria y, por tanto, la vida en nuestro planeta. La tasa fotosintética viene determinada por la ‘maquinaria’ bioquímica y las resistencias difusivas al paso del CO2 desde la atmósfera hasta su fijación en el interior de los cloroplastos. Históricamente la mayor resistencia difusiva se ha atribuido al cierre estomático, sin embargo ahora sabemos, debido a las mejoras en las técnicas experimentales, que existe también una resistencia grande que se opone a la difusión del CO2 desde los espacios intercelulares a los lugares de carboxilación. Esta resistencia, llamada normalmente por su inversa: la conductancia del mesófilo (gm), puede ser igual o incluso superior a la resistencia debida por el cierre estomático. En la presente tesis doctoral he caracterizado la limitación que ejerce la resistencia del mesófilo a la fijación de CO2 en diversas especies forestales y en distintos momentos de su ciclo biológico. En la fase de regenerado, hemos estudiado tres situaciones ambientales relevantes en el mayor éxito de su supervivencia, que son: el déficit hídrico, su interacción con la irradiancia y el paso del crecimiento en la sombra a mayor irradiancia, como puede suceder tras la apertura de un hueco en el dosel forestal. En la fase de arbolado adulto se ha caracterizado el estado hídrico y el intercambio gaseoso en hojas desarrolladas a distinta irradiancia dentro del dosel vegetal durante tres años contrastados en pluviometría. Para cada tipo de estudio se han empleado las técnicas ecofisiológicas más pertinentes para evaluar el estado hídrico y el intercambio gaseoso. Por su complejidad y la falta de un método que permita su cuantificación directa, la gm ha sido evaluada por los métodos más usados, que son: la discriminación isotópica del carbono 13, el método de la J variable, el método de la J constante y el método de la curvatura. Los resultados más significativos permiten concluir que la limitación relativa a la fotosíntesis por la conductancia estomática, del mesófilo y bioquímica es dependiente de la localización de la hoja en el dosel forestal. Por primera vez se ha documentado que bajo estrés hídrico las hojas desarrolladas a la sombra estuvieron más limitadas por una reducción en la gm, mientras que las hojas desarrolladas a pleno sol estuvieron más limitadas por reducción mayor de la conductancia estomática (gsw). Encontramos buena conexión entre el aparato fotosintético foliar y el sistema hídrico debido al alto grado de correlación entre la conductancia hidráulica foliar aparente y la concentración de CO2 en los cloroplastos en distintas especies forestales. Además, hemos mostrado diferentes pautas de regulación del intercambio gaseoso según las particularidades ecológicas de las especies estudiadas. Tanto en brinzales crecidos de forma natural y en el arbolado adulto como en plántulas cultivadas en el invernadero la ontogenia afectó a las limitaciones de la fotosíntesis producidas por estrés hídrico, resultando que las limitaciones estomáticas fueron dominantes en hojas más jóvenes mientras que las no estomáticas en hojas más maduras. La puesta en luz supuso un gran descenso en la gm durante los días siguientes a la transferencia, siendo este efecto mayor según el grado de sombreo previo en el que se han desarrollado las hojas. La aclimatación de las hojas a la alta irradiancia estuvo ligada a las modificaciones anatómicas foliares y al estado de desarrollo de la hoja. El ratio entre la gm/gsw determinó la mayor eficiencia en el uso del agua y un menor estado oxidativo durante la fase de estrés hídrico y su posterior rehidratación, lo cual sugiere el uso de este ratio en los programas de mejora genética frente al estrés hídrico. Debido a que la mayoría de modelos de estimación de la producción primaria bruta (GPP) de un ecosistema no incluye la gm, los mismos están incurriendo en una sobreestimación del GPP particularmente bajo condiciones de estrés hídrico, porque más de la mitad de la reducción en fotosíntesis en hojas desarrolladas a la sombra se debe a la reducción en gm. Finalmente se presenta un análisis de la importancia en las estimas de la gm bajo estrés hídrico de la refijación del CO2 emitido en la mitocondria a consecuencia de la fotorrespiración y la respiración mitocondrial en luz. ABSTRACT Photosynthesis is the biological process that supports primary production and, therefore, life on our planet. Rates of photosynthesis are determined by biochemical “machinery” and the diffusive resistance to the transfer of CO2 from the atmosphere to the place of fixation within the chloroplasts. Historically the largest diffusive resistance was attributed to the stomata, although we now know via improvements in experimental techniques that there is also a large resistance from sub-stomatal cavities to sites of carboxylation. This resistance, commonly quantified as mesophyll conductance (gm), can be as large or even larger than that due to stomatal resistance. In the present PhD I have characterized the limitation exerted by the mesophyll resistance to CO2 fixation in different forest species at different stages of their life cycle. In seedlings, we studied three environmental conditions that affect plant fitness, namely, water deficit, the interaction of water deficit with irradiance, and the transfer of plants grown in the shade to higher irradiance as can occur when a gap opens in the forest canopy. At the stage of mature trees we characterized water status and gas exchange in leaves developed at different irradiance within the canopy over the course of three years that had contrasting rainfall. For each study we used the most relevant ecophysiological techniques to quantify water relations and gas exchange. Due to its complexity and the lack of a method that allows direct quantification, gm was estimated by the most commonly used methods which are: carbon isotope discrimination, the J-variable, constant J and the curvature method The most significant results suggest that the relative limitation of photosynthesis by stomata, mesophyll and biochemistry depending on the position of the leaf within the canopy. For the first time it was documented that under water stress shaded leaves were more limited by a reduction in gm, while the sun-adapted leaves were more limited by stomatal conductance (gsw). The connection between leaf photosynthetic apparatus and the hydraulic system was shown by the good correlations found between the apparent leaf hydraulic conductance and the CO2 concentration in the chloroplasts in shade- and sun-adapted leaves of several tree species. In addition, we have revealed different patterns of gas exchange regulation according to the functional ecology of the species studied. In field grown trees and greenhouse-grown seedlings ontogeny affected limitations of photosynthesis due to water stress with stomatal limitations dominating in young leaves and nonstomatal limitations in older leaves. The transfer to high light resulted in major decrease of gm during the days following the transfer and this effect was greater as higher was the shade which leaves were developed. Acclimation to high light was linked to the leaf anatomical changes and the state of leaf development. The ratio between the gm/gsw determined the greater efficiency in water use and reduced the oxidative stress during the water stress and subsequent rehydration, suggesting the use of this ratio in breeding programs aiming to increase avoidance of water stress. Because most models to estimate gross primary production (GPP) of an ecosystem do not include gm, they are incurring an overestimation of GPP particularly under conditions of water stress because more than half of An decrease in shade-developed leaves may be due to reduction in gm. Finally, we present an analysis of the importance of how estimates of gm under water stress are affected by the refixation of CO2 that is emitted from mitochondria via photorespiration and mitochondrial respiration in light.

Respuesta fisiológica y agronómica del olivar en seto (cv. Arbequina) a la orientación de filas y al ancho de calle: relación entre la producción y calidad del aceite y la radiación interceptada

En las últimas dos décadas, los productores han plantado olivares en seto para lograr la mecanización de la poda y en especial de la cosecha, reducir los costes de mano de obra y permitir intervenciones de manejo rápidas y oportunas. Los olivares se desarrollaron en ausencia del conocimiento científico, sobre el diseño óptimo de la estructura de la copa, necesario para incrementar la producción y calidad del aceite. En contraste, con los árboles muy espaciados y distribuidos uniformemente de las plantaciones tradicionales, en el olivar en seto hay una marcada variabilidad espacial y temporal de la radiación disponible en función del diseño de la plantación. Así, conocer la respuesta fisiológica y productiva del olivo a la radiación resulta fundamental en el olivar en seto. La orientación de las filas y el ancho de calle son aspectos que se deciden en el diseño de las plantaciones en seto. Ambos aspectos modifican la radiación interceptada por la canopia y, por lo tanto, pueden incidir en la productividad y calidad del aceite. Una vez realizada la plantación no pueden ser modificados, y así las ventajas o desventajas permanecerán fijas durante toda la vida productiva del olivar. A pesar de esto, el impacto de la orientación de las filas y el ancho de calle han recibido poca atención en olivos y en la mayoría de los frutales conducidos en seto. Por todo ello, los objetivos principales de esta tesis fueron, (i) evaluar el efecto de la orientación del seto y del ancho de calle, sobre la productividad y calidad del aceite, (ii) evaluar un modelo que estime la radiación dentro de la canopia. Este modelo permitirá cuantificar las relaciones entre la radiación y los componentes del rendimiento y calidad del aceite de olivares en setos con un amplio rango de estructuras y (iii) conocer la variabilidad en las características de las hojas (morfológicas y fisiológicas) y de los tejidos del fruto (tamaño y composición) en diferentes posiciones de la copa de los setos. Para ello, se dispuso de 3 ensayos de olivar en seto (cv. Arbequina) implantados en 2008 en el municipio de La Puebla de Montalbán, Toledo. La primera cosecha fue en 2010 y a partir del 2012 los setos formaron una copa continua. A partir de ese año, los setos se mantuvieron mediante poda, con similar ancho (~1 m) y altura (~2,5 m), acordes a las dimensiones de la cosechadora vendimiadora. En los años 2012 y 2013 se estudió en profundidad la respuesta de las plantas de estos ensayos. En el ensayo 1, los setos fueron plantados con cuatro orientaciones de filas: N–S, NE–SO, NO–SE y E–O y el mismo ancho de calle (4 m). En los otros dos ensayos, los setos fueron plantados con tres anchos de calle (5,0, 4,0 y 2,5 m), y con filas orientadas N–S (ensayo 2) y E–O (ensayo 3). La respuesta de la orientación de las filas se evaluó a nivel de seto y de estratos del seto (alturas y caras), a través de mediciones del crecimiento de brotes, componentes reproductivos, características y temperatura del fruto, estado hídrico del suelo y de las plantas, fotosíntesis neta de las hojas y contenido de ácidos grasos. Los setos orientados NE–SO (2,7 t/ha) lograron la mayor producción de aceite, que fue significativamente más alta que la de los setos E–O (2,3 t/ha). La producción de aceite de los setos E–O no se diferenció estadísticamente de los setos N–S (2,5 t/ha). Las diferencias productivas entre orientaciones fueron explicadas por el número de frutos en cosecha, a su vez la variación en el número de frutos estuvo asociada al efecto de la orientación de las filas sobre el número de yemas desarrolladas y el porcentaje de inflorescencias fértiles. Las hojas en las caras iluminadas de los setos NE–SO y N–S presentaron mayor tasa fotosintética a la mañana (~10.0 h) que los setos E–O, en el año 2012, pero no en 2013. La orientación de las filas no tuvo un efecto significativo en el contenido de ácidos grasos de los aceites extraídos, esto ocurrió a pesar de variaciones en la temperatura interna de los frutos (3 °C) y de la radiación (40%) entre las distintas caras de los setos. La orientación del seto afectó significativamente al contenido relativo de agua del suelo, donde setos E–O presentaron valores más altos (12%) que setos N–S durante el verano y otoño. Sin embargo, el potencial hídrico de tallo fue similar entre orientaciones. En los ensayos 2 y 3, se evaluó el efecto que produce, a nivel de seto y de estratos (caras y alturas), reducir el ancho de calle de 5,0 a 4,0 y 2,5 m, en un seto orientado N–S y otro E–O, respectivamente. La relación entre altura/ancho de calle libre aumentó 0,6 a 0,8 y 1,6, al reducir 5,0, 4,0 y 2,5 m el ancho de calle, mientras la longitud de seto y el volumen de copa por hectárea incrementó 100% al reducir de 5,0 a 2,5 m, el ancho de calle. En los setos orientados N–S, la producción de aceite por ha acumulada en 4 campañas, incrementó significativamente un 52 %, al reducir de 5,0 a 2,5 m el ancho de calle. Los setos N–S con calle más estrecha (2,5 m) tuvieron un 19% menos frutos que los setos con calle más ancha (5,0 m) y a su vez el 60% de los mismos se localizaron los estratos altos de la canopia de los setos con calles estrecha en comparación al 40% en setos con calle de 5,0 m. En los estratos más bajos de los setos con calles de 2,5m hubo menor crecimiento de los brotes y los frutos tuvieron menor peso seco, contenido de aceite y madurez, que los frutos en los estratos bajos de los setos a 5,0 m. Los componentes del rendimiento y características de los frutos (agua y madurez) fueron similares entre la caras E y O, independientemente del ancho de calle. En los setos orientados E–O, la producción de aceite por ha acumulada en 4 campañas, no respondió significativamente al ancho de calle, debido a una disminución significativa en el número de frutos y producción de aceite por m de seto, al reducir de 5,0 a 2,5 m, el ancho de calle. En los setos orientados E–O, con calles de 5,0 m, los frutos presentaron similar peso seco, contenido de aceite y agua, en las caras S y N, sin embargo, cuando la calle fue reducida a 2,5, los frutos de la cara S fueron más pesado y maduros que en la cara N. Independientemente del ancho de calle y de la orientación del seto, el aceite presentó mayor contenido de ácidos palmitoleico, palmítico, esteárico y linoleico en los frutos del estrato más alto de la canopia disminuyendo hacia la base. En contraste, el contenido de ácido oleico aumentó desde el estrato más alto hacia la base de los setos. Las diferencias en el contenido de ácidos grasos entre la parte alta y baja de los setos, incrementó al reducir el ancho de calle en los setos N–S, pero no en los E-O. En conclusión, en olivares en seto, reducir el ancho de calle permite incrementar la producción de aceite, en setos orientados N–S, pero no en E–O. Un modelo que estima la cantidad y distribución de la radiación en toda la copa del seto, fue utilizado para estimar la radiación interceptada en distintos estratos del seto. El modelo requiere un valor del coeficiente de extinción (k) para estimar la transmisión de radiación a través de la copa, el cual fue obtenido experimentalmente (k=1,2). Utilizando los datos del ensayo 1, un único modelo lineal relacionó el peso seco y el rendimiento graso de setos con la radiación interceptada por los distintos estratos de setos con cuatro orientaciones de filas. La densidad de frutos fue también relacionada con la radiación, pero más débilmente. En los setos orientados N–S, plantados con tres anchos de calles, (ensayo 2) el contenido de ácidos palmitoleico y linoleico del aceite incrementó linealmente con el incremento de la radiación interceptada, mientras el contenido ácido oleico disminuyó linealmente con el incremento de la radiación. El contenido de ácidos grasos del aceite no estuvo relacionado con la radiación interceptada en setos orientados E–O (Ensayo 3). En los setos N–S y E–O, plantados con anchos de calle de 2,5 m, se estudiaron las interacciones entre la radiación y características de las hojas, número de fruto, tamaño y composición de los frutos a nivel de órgano, tejido y células. Independientemente de la orientación del seto, el área y el contenido de clorofila de las hojas incrementaron significativamente en los estratos más bajos de los setos. Mientras, las hojas de los estratos medios del seto presentaron mayor capacidad fotosintética que en los estratos bajos y alto de los setos. Los estratos del seto que interceptaron más radiación produjeron frutos con mayor tamaño y contenido de aceite en el mesocarpo, sin efectos sobre el tamaño y composición del endocarpo. A nivel celular, los frutos expuestos a mayor nivel de radiación desarrollaron en el mesocarpo células de mayor tamaño en comparación a frutos menos expuestos, mientras el número de células no fue afectado. Adicionalmente, el número y tamaño de las células estuvo relacionado con la composición del mesocarpo en términos de aceite, agua y peso seco menos aceite. Esta tesis, contribuye, desde una perspectiva integral del cultivo del olivo, a cuantificar el impacto de la orientación y ancho de calle sobre la producción y calidad del aceite en olivares conducidos en setos. El análisis y discusión de la relación entre la radiación y los componentes del rendimiento y calidad del aceite, puede ayudar a diseñar plantaciones en seto con dimensiones óptimas para la intercepción de la radiación. ABSTRACT In the last two decades, olive hedgerow system has been established by commercial growers to allow continuous mechanized pruning and especially harvest, reduce costs of manual labour and allow more rapid and timely management interventions. The adoption of hedgerow was done in the absence of adequate scientific knowledge of the impact of this orchard structure and associated mechanization on tree response, yield and quality, after centuries in low-density orchards and open-formed trees. The row orientation and width alley are fundamental aspects in the hedgerow design and have been scarcely studied in olive. Both aspects modify the radiation intercepted by the canopy, and consequently the productivity and oil quality, and once defined in orchard planting cannot be changed, so advantages and disadvantages remain fixed for the lifespan of the orchard. The main objectives of this thesis were to (i) evaluate the impact of the row orientation and width alley on productivity and oil quality by the measurements of profile of the determining processes of shoot growth, fruit temperature, yield components and fruit and oil characteristics on opposite sides of olive hedgerows. Additionally, the effect of row orientation on the plant water status was also evaluated; (ii) evaluate a mathematical model for estimating the radiation within the canopy and quantify the relationships between the radiation estimated and yield components and oil quality in olive hedgerows under wide range of structures and; (iii) determine the variability in the characteristics of the leaves (morphological and physiological) and fruit tissues (size and composition) in different positions of the hedgerows canopy. Three plots of olive hedgerows (cv. Arbequina) planted in 2008 in La Puebla de Montalbán, Toledo were evaluated during the 2012 and 2013 seasons. The hedgerows were maintained by lateral pruning and topping with the same width (1 m) and height (2.5 m) compatible with the intended harvester. In a plot (experiment 1), the hedgerows were planted with the same width alley (4 m) and four row orientations: N–S, NE–SW, NW–SE and E–W. Other two plots (Experiments 2 and 3) separated by approximately 100 m were planted with N–S and E–O oriented rows and three alley widths in each orientation: 5.0, 4.0 and 2.5 m. In the exp. 1, maximum fruit yield were achieved by NE–SW and NW–SW (15.7 t/ha). Of these, NE–SW achieved the highest oil yield (2.7 t/ha). There were no differences in fruit or oil yield between N–S (2.5 t oil/ha) and E–W (2.3 t oil/ha) orientations. Fruit number was the most important component to explain these differences, by previous influence on number of bud developed and percentage of fertile inflorescences. Fruit maturity and oil quality on both sides of the hedgerows were not affected by row orientation. This occurred despite significant variations in the internal fruit temperature, which was closely related to the irradiance received by the canopy and the time of day. Additionally, row orientation significantly affected the relative water content of the soil, where E–W oriented hedgerows showed consistently higher values than N–S during summer-autumn season. The stem water potential at midday, however, was similar between orientations, revealing possible lower water consumption of E–W than N–S oriented hedgerows. In the exp. 2, regardless of row orientation, reduction of row spacing from 5.0 to 4.0 and 2.5 m increases the ratio of canopy depth to free alley width (Al/An) from 0.6 to 0.8 and 1.6, respectively, and ads 25 and 100 % more hedgerow length per ha. In N–S oriented hedgerows, oil production per ha increased significantly by 14 and 52 % in 4.0 m and 2.5 m relative to 5.0 m row spacing, the effect being proportionally less than the increase in hedgerow length per ha. Hedgerows spaced 2.5 m with Al/An = 1.6 produced relatively fewer fruits per unit length than did wider spacings and were preferentially distributed in upper layers. Fruits located at the bottom of the canopy were smaller, with lower oil content and were less mature. In E–W oriented hedgerows, oil production per ha did not respond significantly to row spacing, despite the doubling of row length from the 5.0 to the 2.5 m row spacing. The explanation was found in fewer fruit per unit length of hedgerow and smaller oil content at 2.5 m than 5.0 m row spacing, averaged over the experimental period. In E–W hedgerows spaced at 5.0 m with Al/An = 0.6, the vertical profiles of fruit characteristics (mass, oil and water contents, and maturity) were similar between opposing sides, but at 4.0 m (Al/An= 0.8) and 2.5 m (Al/An=1.6) spacings, fruits on the S side were heavier and more mature than on N side. The oil extracted from fruits harvested at different heights of N–S and E–W oriented hedgerows showed higher palmitoleic, palmitic, stearic and linoleic contents at the canopy top decreasing toward base. The oleic content was reverse, increased from top to base. In N–S hedgerows, vertical gradients increased by reducing the alley width, but not in the E–W oriented hedgerows. The simulation of internal canopy irradiance was related in a single relationship (R2 = 0.63) to the vertical profiles of fruit weight and oil content of olive hedgerows with wide range of structures. The density of fruits was also associated with the irradiance but more weakly (R2 = 0.27), and revealed a more complex response involving changes in the vegetative structure by canopy management (topping) and the effect of radiation on the previous sequence that defines the number of fruits. The vertical profiles of oil quality traits were closely associated to canopy irradiance, but only when the N–S oriented hedgerows were considered. The contents of palmitoleic and linoleic acid in the oil increased linearly when intercepted irradiance increased from 9 to 19 mol PAR/m2. In contrast, oleic content decreased linearly in this irradiance range. Additionally, we advanced knowledge regarding the interactions among irradiance and leaf, fruit number, size and composition at organ-, tissue- and cellular- levels. The irradiance received at different positions in the canopy strongly affected the leaf area and chlorophyll content, and mesocarp size and composition (water and oil), without effects on endocarp size and composition. At the cellular level, light-exposed fruit developed larger mesocarp cells than shaded fruits, but cell number was not affected. Our results indicate that cell number and size are related to mesocarp composition in term of oil, water, and dry weight menus oil, although the specific manner in which they interact remains to be determined. This research contributes from an integral perspective of olive growing to quantify the impact of row orientation and width alley on productivity and oil quality in hedgerows systems. The analysis and discussion of the relationships between radiation and yield components and oil quality can help understand the impact of design olive hedgerows in general and in a wide range of environmental conditions.