Novel parameters for evaluating the Spatial and Thematic accuracy of land cover maps.

Se proponen novedosas fórmulas para evaluar la certeza de la cartografía

Land cover dynamics and climate change implications on water resources in South Pacific Costa Rica

Water is fundamental to human life and the availability of freshwater is often a constraint on human welfare and economic development. Consequently, the potential effects of global changes on hydrology and water resources are considered among the most severe and vital ones. Water scarcity is one of the main problems in the rural communities of Central America, as a result of an important degradation of catchment areas and the over-exploitation of aquifers. The present Thesis is focused on two critical aspects of global changes over water resources: (1) the potential effects of climate change on water quantity and (2) the impacts of land cover and land use changes on the hydrological processes and water cycle. Costa Rica is among the few developing countries that have recently achieved a land use transition with a net increase in forest cover. Osa Region in South Pacific Costa Rica is an appealing study site to assess water supply management plans and to measure the effects of deforestation, forest transitions and climate change projections reported in the region. Rural Community Water Supply systems (ASADAS) in Osa are dealing with an increasing demand of freshwater due to the growing population and the change in the way of life in the rural livelihoods. Land cover mosaics which have resulted from the above mentioned processes are characterized by the abandonment of marginal farmland with the spread over these former grasslands of high return crops and the expansion of secondary forests due to reforestation initiatives. These land use changes have a significant impact on runoff generation in priority water-supply catchments in the humid tropics, as evidenced by the analysis of the Tinoco Experimental Catchment in the Southern Pacific area of Costa Rica. The monitoring system assesses the effects of the different land uses on the runoff responses and on the general water cycle of the basin. Runoff responses at plot scale are analyzed for secondary forests, oil palm plantations, forest plantations and grasslands. The Oil palm plantation plot presented the highest runoff coefficient (mean RC=32.6%), twice that measured under grasslands (mean RC=15.3%) and 20-fold greater than in secondary forest (mean RC=1.7%). A Thornthwaite-type water balance is proposed to assess the impact of land cover and climate change scenarios over water availability for rural communities in Osa Region. Climate change projections were obtained by the downscaling of BCM2, CNCM3 and ECHAM5 models. Precipitation and temperature were averaged and conveyed by the A1B, A2 and B1 IPCC climate scenario for 2030, 2060 and 2080. Precipitation simulations exhibit a positive increase during the dry season for the three scenarios and a decrease during the rainy season, with the highest magnitude (up to 25%) by the end of the 21st century under scenario B1. Monthly mean temperature simulations increase for the three scenarios throughout the year with a maximum increase during the dry season of 5% under A1B and A2 scenarios and 4% under B1 scenario. The Thornthwaite-type Water Balance model indicates important decreases of water surplus for the three climate scenarios during the rainy season, with a maximum decrease on May, which under A1B scenario drop up to 20%, under A2 up to 40% and under B1 scenario drop up to almost 60%. Land cover scenarios were created taking into account current land cover dynamics of the region. Land cover scenario 1 projects a deforestation situation, with forests decreasing up to 15% due to urbanization of the upper catchment areas; land cover scenario 2 projects a forest recovery situation where forested areas increase due to grassland abandonment on areas with more than 30% of slope. Deforestation scenario projects an annual water surplus decrease of 15% while the reforestation scenario projects a water surplus increase of almost 25%. This water balance analysis indicates that climate scenarios are equal contributors as land cover scenarios to future water resource estimations.

Land Cover Mapping by an Optimized Object-Oriented Approach. Case of Study Mediterranean Landscapes

Remote sensing information from spaceborne and airborne platforms continues to provide valuable data for different environmental monitoring applications. In this sense, high spatial resolution im-agery is an important source of information for land cover mapping. For the processing of high spa-tial resolution images, the object-based methodology is one of the most commonly used strategies. However, conventional pixel-based methods, which only use spectral information for land cover classification, are inadequate for classifying this type of images. This research presents a method-ology to characterise Mediterranean land covers in high resolution aerial images by means of an object-oriented approach. It uses a self-calibrating multi-band region growing approach optimised by pre-processing the image with a bilateral filtering. The obtained results show promise in terms of both segmentation quality and computational efficiency.

Linking land cover changes in the subalpine and montane belts to changes in a torrential river.

Channel cbanges are the consequence of cbanges in sediment yield from the slopes and in the connectivity between slopes and channels bccause of distinct land use and climate inJpacts. In Ibis study, we investigated the characteristics and evolution of a short reach in the headwater of the !juez River, central-soutbem Pyrenees. Assessment of a series of sedimentar)' and geomorphic structures confirtned major cbanges to the valley boUom, mainly related to changes in the intensity of human activity. 'The oldest sedimentar)' structure is a terrace leve! located 3 10 4 m above the current alluvial plain. General deforestation, overgrazing and recwring tires in the montane belt ( 1100-1600 m a.sJ.) have led 10 increased soil erosioo and connectivity, and to tbe triggering of debris llows thal have been deposited on the fluvial tmace. Woody fragments from within the debris llows were dated using acceleraror mass spectrOmetry '"e radiocaroon tcchniques (AMS), yielding ages between 1 00 and 115 cal years BP, whicb coincides with tbe period of maximum deforestalion and human density in the Pyrenees. Depopulation and fannland abaodonment since tbe beginning of tbe 2001 oenrury has resulled in generalliJcod natural and artificial reforestation, a shrinkage of the eroded arcas aod a decline in connectivi¡y bdween slopes and the channel. 1be rnost impor1an1 consequence has been cbannel incision and oarrowing, and the development of a sedimed annour !ayer. Active sedimenl b8llSpOI1 is continui.Qg, although there has been a decrease in sed.iment yield from the slopes. Copyright O 2014 John Wiley & Soos, Ltd.

Integrating periurban agrarian ecosystem services into spatial planning to cope with urban pressure

According to Corine Land Cover databases, in Europe between 1990 and 2000,77% of new artificial surfaces were built on previous agrarian areas. Urban sprawl ¡s far from being under control, between 2000 and 2006 new artificial land has grown in larger proportion than the decade before. In Spain, like in most countries, the impact of urban sprawl during the last decades has been especially significant in periurban agrarian spaces: between 2000 and 2006, 73% of new artificial surfaces were built on previous agrarian areas. The indirect impact of this trend has been even more relevant, as the expectations of appreciation in the value of land after new urban developments reinforce the ongoing trend of abandonment of agricultural land. In Madrid between 1980 and 2000 the loss of agricultural land due to abandonment of exploitation was 2-fold that due to transformation into urban areas. By comparing four case studies: Valladolild, Montpellier.Florence and Den Haag, this paper explores if urban and territorial planning may contribute to reduce urban pressure on the hinterland. In spite of their diversity, these regions have in common a relative prosperity arising from their territorial endowments, though their landscapes are still under pressure. The three last ones have been working for years on mainstream concepts like multifunctional agriculture. The systematic comparison and the analysis of successful approaches provide some clues on how to reconsider urban planning in order to preserve agricultural land. The final remarks highlight the context in which public commitment, legal protection instruments and financial strategies may contribute to the goals of urban, peri-urban or regional planning about fostering agrarian ecosystem services

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.

Análisis de edificación y sostenibilidad

El desarrollo sostenible, es un concepto amplio que afecta a todos los ámbitos de nuestro entorno natural, social y económico. Existen multitud de interrelaciones entre los diferentes ámbitos de éste que nos permitirán realizar un análisis tangible y apoyado en datos sobre las interacciones que se dan. La intención es realizar un análisis de los datos más relevantes de los indicadores del desarrollo sostenible que establecen relaciones con las edificaciones, a nivel del sector de la construcción y del parque residencial y terciario existente en España. No debemos olvidar el análisis de las leyes de la edificación, y de como se ha ido incorporando el desarrollo sostenible en los últimos años y que grado de repercusión ha alcanzado. Para analizar la edificación y sostenibilidad nos debemos respaldar en datos, los cuales se recogen en los indicadores más destacados en materia de desarrollo sostenible como son los: de la EEDS, la UE-EDS, el Banco Público de Indicadores Ambientales (BPIA) y los elaborados por el OSE. El carácter tan singular de la relación que se quiere establecer, provoca la búsqueda de otros muchos datos que establezcan relaciones más íntimas entre las edificaciones o la construcción. Estos datos se han sacado de las fuentes institucionales más relevantes, siempre con la finalidad de dar luz a las relaciones que se establecen entre los edificios y el sector de la construcción con los diferentes aspectos más notables con los que interactúan, en términos de desarrollo sostenible. Las relaciones se han establecido unas veces con la construcción como sector económico, otras atendiendo simplemente a las edificaciones como punto de consumo o desecho (hogares) de ciertos ítems. Las fuentes usadas para la elaboración de este análisis, aparte de los indicadores de sostenibilidad han sido: el Eurostat, proyecto Europeo CORINE Land Cover (CLC), el Instituto Nacional de Estadística (INE), Banco Público de Indicadores Ambientales (BPIA) del MARM, datos de documentos de diversos Ministerios, organismos e instituciones.

Space-time analysis of the urban-rural grafient in the metropolitan area of Madrid and ist surroundings -Guadalajara and Toledo.

According to UN provisions in the period from 2007 to 2050 world population will grow up to 9200 million people. In fact, for the first time in history, in the year 2008 world urban population became higher than rural population. The increase of urban areas and their transport infrastructures has influenced agricultural land use due to their irreversible change, especially when they remain as periurban vacant land, losing their character and identity. In the Europe of the nineties, the traditional urban-rural gradient, characterized by a neat contact between both land types, has become so complex that it has change to a gradient in which it is difficult to separate urban and rural land uses. [Antrop 2004]. A literature review has been made on methodologies used for the urban-rural gradient analysis. One of these methodologies was selected that integrates ecological characterization based on the use of spatial metrics and geographical characterization based on spatial components. Cartographical sources used were Corine Land Cover at 1: 100000 scale and the Spanish Land Use Information System at 1:25000 scale. Urban-rural gradient paradigm is an analysis methodology, coming from landscape ecology, which enables to investigate how urbanization provokes changes in ecological patterns and processes into landscape. [Hahs and McDonnell 2006].The present research adapt this methodology to study the urban-rural gradient in the outskirts of Madrid, Toledo and Guadalajara. Both scales (1:25000 and 1:100000) were simultaneously used to reach the next objectives: 1) Analysis of landscape pattern dynamics in relation to distance to the town centre and major infrastructures. 2) Analysis of landscape pattern dynamics in the fringe of protected areas. The paper presents a new approach to the urban-rural relationship which allows better planning and management of urban áreas.

Análisis espacio-temporal del gradiente urbano-rural del sur de la Región Metropolitana de Madrid y su entorno. Caracterización de los procesos y patrones paisajísticos acaecidos en el período 1990-2006

El objetivo de esta tesis es proponer una metodología capaz de cuantificar la dinámica paisajística a lo largo del gradiente urbano – rural perteneciente al sur de la Región Metropolitana de Madrid y su entorno. Esta tesis se estructura en ocho capítulos, planos y anejos. El primero se refiere a los antecedentes tanto conceptuales como metodológicos. Los primeros se centran en los diversos enfoques existentes en relación al estudio de la dinámica paisajística, con el objetivo de encontrar los puntos en común existentes entre el enfoque del metabolismo social, el de la ecología del paisaje y el humanista, para obtener un diagnóstico que permita entender la complejidad de la realidad a la que esta Tesis se enfrenta. Los segundos se centran en los antecedentes de carácter metodológico que también desde diversos enfoques han abordado el análisis de la dinámica paisajística. El segundo capítulo se centra en los objetivos concretos derivados del objetivo general ya expresado, la tesis considera que para comprender y cuantificar la dinámica hay que identificar en primer lugar los procesos de transformación, como manifestación espacial de los factores socioeconómicos y naturales responsables en última instancia de la variación de los patrones paisajísticos existentes. En segundo lugar se identifican los patrones paisajísticos con el objetivo de analizar sus características espaciales y su evolución en el período analizado. Por último se identifican los procesos paisajísticos, es decir qué tipos de variaciones espaciales se producen en los patrones paisajísticos como consecuencia de los procesos de transformación identificados así como su pauta de distribución a lo largo del gradiente urbano ‐ rural. El tercer capítulo se dedica a la caracterización del ámbito de estudio, ésta se extiende al sur del límite del suelo urbano de la capital madrileña en el año 1990, comprende la totalidad de los municipios madrileños que contactan con los municipios castellano – manchegos que se encuentran en el área de influencia de la capital, abarcando el área 9.968 km2. El cuarto capítulo se centra en la metodología. Como material de partida se ha utilizado en la cartografía del Corine Land Cover y como herramienta de análisis se ha utilizado los Sistemas de Información Geográfica. En primer lugar se identifican los procesos de transformación, acaecidos en los períodos 1990 – 2000 y 2000 – 2006, mediante la aplicación de matrices de transición. Se han identificado cuatro tipos de procesos dinámicos: Urbanización, abandono, renaturalización y agrarización. Se ha realizado un análisis de indicadores compuestos lo que ha permitido identificar los tipos de patrones paisajísticos existentes a lo largo del gradiente urbano – rural. Del mismo modo se ha calculado la variación de los indicadores individuales para identificar los procesos paisajísticos mediante el análisis de indicadores compuestos que se produjeron en el período 1990 – 2000 y 2000 – 2006. En el quinto capítulo se aportan los resultados tanto de carácter cuantitativo como gráfico de los tres componentes analizados tanto de forma independiente como integrada. En el sexto capítulo se describen las conclusiones producto de la investigación realizada. En el séptimo capítulo se identifican qué líneas de investigación podrían desarrollarse en el futuro para continuar la línea de investigación iniciada con esta tesis. ABSTRACT The aim of this thesis was to propose a methodology to characterize landscape dynamics along the urban – rural gradient in the south Madrid area. It´s structured in eight chapters, planes and annexes: the first one describes previous research. Firstly to make a diagnosis of the effects of landscape dynamic we have performed an integrated analysis from social metabolism, landscape ecology and the humanistic point of view. Secondly we have focused on previous methodological research mainly developped by landscape ecology. The second chapter focuses on specific objectives derived from the general objective. The thesis considers that to understand and quantify landscape dynamics must first identify the transformation processes: spatial manifestation of natural and socioeconomics factors that induce the change of landscape patterns. Secondly landscape patterns have been identified in order to analyze their spatial characteristics and evolution. Finally the landscape processes have been identified, i.e. what kind of spatial variations cause changes in landscape patterns along urban – rural gradient. The third chapter describes the study area. The study area occupies 9968 km2. It covers the area to the south of Madrid’s 1990 urban land area, and takes in the southeast of the Madrid Autonomous Region plus all the municipal areas of the Castilla–La Mancha Autonomous Region directly influenced by the expansion of Madrid. The fourth chapter contains the methodology. To identify the changes in the landscape of the study area, the land cover data for the area held in the CORINE LandCover Project Database was examined. To characterize the transformations processes in the period 1990 ‐ 2000 and 2000 – 2006, transition matrices were constructed. We have identified four clear changes: Urbanization, renaturalization, abandonment and agrarianization. We have characterized landscape patterns using composite indicators by integrating individual spatial metrics. Similarly we have characterized landscape processes using composite indicators by integrating the variation of individual spatial metrics. Chapter fifth includes the results, both for each component and its final integration. The conclusions of this research have been described in the sixth chapter. The seventh chapter describes what kind of investigations could be done in the future.

Análisis espacio-temporal del gradiente urbano-rural del sur de la Región Metropolitana de Madrid y su entorno. Caracterización de los procesos y patrones paisajísticos acaecidos en el período 1990-2006 (versión revisada y corregida)

El objetivo de esta tesis es proponer una metodología capaz de cuantificar la dinámica paisajística a lo largo del gradiente urbano – rural perteneciente al sur de la Región Metropolitana de Madrid y su entorno. Esta tesis se estructura en ocho capítulos, planos y anejos. El primero se refiere a los antecedentes tanto conceptuales como metodológicos. Los primeros se centran en los diversos enfoques existentes en relación al estudio de la dinámica paisajística, con el objetivo de encontrar los puntos en común existentes entre el enfoque del metabolismo social, el de la ecología del paisaje y el humanista, para obtener un diagnóstico que permita entender la complejidad de la realidad a la que esta Tesis se enfrenta. Los segundos se centran en los antecedentes de carácter metodológico que también desde diversos enfoques han abordado el análisis de la dinámica paisajística. El segundo capítulo se centra en los objetivos concretos derivados del objetivo general ya expresado, la tesis considera que para comprender y cuantificar la dinámica hay que identificar en primer lugar los procesos de transformación, como manifestación espacial de los factores socioeconómicos y naturales responsables en última instancia de la variación de los patrones paisajísticos existentes. En segundo lugar se identifican los patrones paisajísticos con el objetivo de analizar sus características espaciales y su evolución en el período analizado. Por último se identifican los procesos paisajísticos, es decir qué tipos de variaciones espaciales se producen en los patrones paisajísticos como consecuencia de los procesos de transformación identificados así como su pauta de distribución a lo largo del gradiente urbano ‐ rural. El tercer capítulo se dedica a la caracterización del ámbito de estudio, ésta se extiende al sur del límite del suelo urbano de la capital madrileña en el año 1990, comprende la totalidad de los municipios madrileños que contactan con los municipios castellano – manchegos que se encuentran en el área de influencia de la capital, abarcando el área 9.968 km2. El cuarto capítulo se centra en la metodología. Como material de partida se ha utilizado en la cartografía del Corine Land Cover y como herramienta de análisis se ha utilizado los Sistemas de Información Geográfica. En primer lugar se identifican los procesos de transformación, acaecidos en los períodos 1990 – 2000 y 2000 – 2006, mediante la aplicación de matrices de transición. Se han identificado cuatro tipos de procesos dinámicos: Urbanización, abandono, renaturalización y agrarización. Se ha realizado un análisis de indicadores compuestos lo que ha permitido identificar los tipos de patrones paisajísticos existentes a lo largo del gradiente urbano – rural. Del mismo modo se ha calculado la variación de los indicadores individuales para identificar los procesos paisajísticos mediante el análisis de indicadores compuestos que se produjeron en el período 1990 – 2000 y 2000 – 2006. En el quinto capítulo se aportan los resultados tanto de carácter cuantitativo como gráfico de los tres componentes analizados tanto de forma independiente como integrada. En el sexto capítulo se describen las conclusiones producto de la investigación realizada. En el séptimo capítulo se identifican qué líneas de investigación podrían desarrollarse en el futuro para continuar la línea de investigación iniciada con esta tesis. ABSTRACT The aim of this thesis was to propose a methodology to characterize landscape dynamics along the urban – rural gradient in the south Madrid area. It´s structured in eight chapters, planes and annexes: the first one describes previous research. Firstly to make a diagnosis of the effects of landscape dynamic we have performed an integrated analysis from social metabolism, landscape ecology and the humanistic point of view. Secondly we have focused on previous methodological research mainly developped by landscape ecology. The second chapter focuses on specific objectives derived from the general objective. The thesis considers that to understand and quantify landscape dynamics must first identify the transformation processes: spatial manifestation of natural and socioeconomics factors that induce the change of landscape patterns. Secondly landscape patterns have been identified in order to analyze their spatial characteristics and evolution. Finally the landscape processes have been identified, i.e. what kind of spatial variations cause changes in landscape patterns along urban – rural gradient. The third chapter describes the study area. The study area occupies 9968 km2. It covers the area to the south of Madrid’s 1990 urban land area, and takes in the southeast of the Madrid Autonomous Region plus all the municipal areas of the Castilla–La Mancha Autonomous Region directly influenced by the expansion of Madrid. The fourth chapter contains the methodology. To identify the changes in the landscape of the study area, the land cover data for the area held in the CORINE LandCover Project Database was examined. To characterize the transformations processes in the period 1990 ‐ 2000 and 2000 – 2006, transition matrices were constructed. We have identified four clear changes: Urbanization, renaturalization, abandonment and agrarianization. We have characterized landscape patterns using composite indicators by integrating individual spatial metrics. Similarly we have characterized landscape processes using composite indicators by integrating the variation of individual spatial metrics. Chapter fifth includes the results, both for each component and its final integration. The conclusions of this research have been described in the sixth chapter. The seventh chapter describes what kind of investigations could be done in the future.

Implications of agricultural production, policy and land use changes on water resource assessment

Una gestión más eficiente y equitativa del agua a escala de cuenca no se puede centrar exclusivamente en el recurso hídrico en sí, sino también en otras políticas y disciplinas científicas. Existe un consenso creciente de que, además de la consideración de las cambiantes condiciones climáticas, es necesaria una integración de ámbitos de investigación tales como la agronomía, planificación del territorio y ciencias políticas y económicas a fin de satisfacer de manera sostenible las demandas de agua por parte de la sociedad y del medio natural. La Política Agrícola Común (PAC) es el principal motor de cambio en las tendencias de paisajes rurales y sistemas agrícolas, pero el deterioro del medio ambiente es ahora una de las principales preocupaciones. Uno de los cambios más relevantes se ha producido con la expansión e intensificación del olivar en España, principalmente con nuevas zonas de regadío o la conversión de olivares de secano a sistemas en regadío. Por otra parte, el cambio de las condiciones climáticas podría ejercer un papel importante en las tendencias negativas de las aportaciones a los ríos, pero no queda claro el papel que podrían estar jugando los cambios de uso de suelo y cobertura vegetal sobre las tendencias negativas de caudal observadas. Esta tesis tiene como objetivo mejorar el conocimiento de los efectos de la producción agrícola, política agraria y cambios de uso de suelo y cobertura vegetal sobre las condiciones de calidad del agua, respuesta hidrológica y apropiación del agua por parte de la sociedad. En primer lugar, el estudio determina las tendencias existentes de nitratos y sólidos en suspensión en las aguas superficiales de la cuenca del río Guadalquivir durante el periodo de 1998 a 2009. Desde una perspectiva de política agraria, la investigación trata de evaluar mediante un análisis de datos de panel las principales variables, incluyendo la reforma de la PAC de 2003, que están teniendo una influencia en ambos indicadores de calidad. En segundo lugar, la apropiación del agua y el nivel de contaminación por nitratos debido a la producción del aceite de oliva en España se determinan con una evaluación de la huella hídrica (HH), teniendo en cuenta una variabilidad espacial y temporal a largo de las provincias españolas y entre 1997 y 2008. Por último, la tesis analiza los efectos de los cambios de uso de suelo y cobertura vegetal sobre las tendencias negativas observadas en la zona alta del Turia, cabecera de la cuenca del río Júcar, durante el periodo 1973-2008 mediante una modelización ecohidrológica. En la cuenca del Guadalquivir cerca del 20% de las estaciones de monitoreo muestran tendencias significativas, lineales o cuadráticas, para cada indicador de calidad de agua. La mayoría de las tendencias significativas en nitratos están aumentando, y la mayoría de tendencias cuadráticas muestran un patrón en forma de U. Los modelos de regresión de datos de panel muestran que las variables más importantes que empeoran ambos indicadores de calidad del agua son la intensificación de biomasa y las exportaciones de ambos indicadores de calidad procedentes de aguas arriba. En regiones en las que el abandono agrícola y/o desintensificación han tenido lugar han mejorado las condiciones de calidad del agua. Para los nitratos, el desacoplamiento de las subvenciones a la agricultura y la reducción de la cuantía de las subvenciones a tierras de regadío subyacen en la reducción observada de la concentración de nitratos. Las medidas de modernización de regadíos y el establecimiento de zonas vulnerables a nitratos reducen la concentración en subcuencas que muestran una tendencia creciente de nitratos. Sin embargo, el efecto de las exportaciones de nitratos procedente de aguas arriba, la intensificación de la biomasa y los precios de los cultivos presentan un mayor peso, explicando la tendencia creciente observada de nitratos. Para los sólidos en suspensión, no queda de forma evidente si el proceso de desacoplamiento ha influido negativa o positivamente. Sin embargo, los mayores valores de las ayudas agrarias aún ligadas a la producción, en particular en zonas de regadío, conllevan un aumento de las tasas de erosión. Aunque la cuenca del Guadalquivir ha aumentado la producción agrícola y la eficiencia del uso del agua, el problema de las altas tasas de erosión aún no ha sido mitigado adecuadamente. El estudio de la huella hídrica (HH) revela que en 1 L de aceite de oliva español más del 99,5% de la HH está relacionado con la producción de la aceituna, mientras que menos del 0,5% se debe a otros componentes, es decir, a la botella, tapón y etiqueta. Durante el período estudiado, la HH verde en secano y en regadío representa alrededor del 72% y 12%, respectivamente, del total de la HH. Las HHs azul y gris representan 6% y 10%, respectivamente. La producción de aceitunas se concentra en regiones con una HH menor por unidad de producto. La producción de aceite de oliva ha aumentado su productividad del agua durante 1997-2008, incentivado por los crecientes precios del aceite, como también lo ha hecho la cantidad de exportaciones de agua virtual. De hecho, las mayores zonas productoras presentan una eficiencia alta del uso y de productividad del agua, así como un menor potencial de contaminación por nitratos. Pero en estas zonas se ve a la vez reflejado un aumento de presión sobre los recursos hídricos locales. El aumento de extracciones de agua subterránea relacionadas con las exportaciones de aceite de oliva podría añadir una mayor presión a la ya estresada cuenca del Guadalquivir, mostrando la necesidad de equilibrar las fuerzas del mercado con los recursos locales disponibles. Los cambios de uso de suelo y cobertura vegetal juegan un papel importante en el balance del agua de la cuenca alta del Turia, pero no son el principal motor que sustenta la reducción observada de caudal. El aumento de la temperatura es el principal factor que explica las mayores tasas de evapotranspiración y la reducción de caudales. Sin embargo, los cambios de uso de suelo y el cambio climático han tenido un efecto compensatorio en la respuesta hidrológica. Por un lado, el caudal se ha visto afectado negativamente por el aumento de la temperatura, mientras que los cambios de uso de suelo y cobertura vegetal han compensado positivamente con una reducción de las tasas de evapotranspiración, gracias a los procesos de disminución de la densidad de matorral y de degradación forestal. El estudio proporciona una visión que fortalece la interdisciplinariedad entre la planificación hidrológica y territorial, destacando la necesidad de incluir las implicaciones de los cambios de uso de suelo y cobertura vegetal en futuros planes hidrológicos. Estos hallazgos son valiosos para la gestión de la cuenca del río Turia, y el enfoque empleado es útil para la determinación del peso de los cambios de uso de suelo y cobertura vegetal en la respuesta hidrológica en otras regiones. ABSTRACT Achieving a more efficient and equitable water management at catchment scale does not only rely on the water resource itself, but also on other policies and scientific knowledge. There is a growing consensus that, in addition to consideration of changing climate conditions, integration with research areas such as agronomy, land use planning and economics and political science is required to meet sustainably the societal and environmental water demands. The Common Agricultural Policy (CAP) is a main driver for trends in rural landscapes and agricultural systems, but environmental deterioration is now a principal concern. One of the most relevant changes has occurred with the expansion and intensification of olive orchards in Spain, taking place mainly with new irrigated areas or with the conversion from rainfed to irrigated systems. Moreover, changing climate conditions might exert a major role on water yield trends, but it remains unclear the role that ongoing land use and land cover changes (LULCC) might have on observed river flow trends. This thesis aims to improve the understanding of the effects of agricultural production, policies and LULCC on water quality conditions, hydrological response and human water appropriation. Firstly, the study determines the existing trends for nitrates and suspended solids in the Guadalquivir river basin’s surface waters (south Spain) during the period from 1998 to 2009. From a policy perspective, the research tries to assess with panel data analysis the main drivers, including the 2003 CAP reform, which are having an influence on both water quality indicators. Secondly, water appropriation and nitrate pollution level originating from the production of olive oil in Spain is determined with a water footprint (WF) assessment, considering a spatial temporal variability across the Spanish provinces and from 1997 to 2008 years. Finally, the thesis analyzes the effects of the LULCC on the observed negative trends over the period 1973-2008 in the Upper Turia basin, headwaters of the Júcar river demarcation (east Spain), with ecohydrological modeling. In the Guadalquivir river basin about 20% of monitoring stations show significant trends, linear or quadratic, for each water quality indicator. Most significant trends of nitrates are augmenting than decreasing, and most significant quadratic terms of both indicators exhibit U-shaped patterns. The panel data models show that the most important drivers that are worsening nitrates and suspended solids in the basin are biomass intensification and exports of both water quality indicators from upland regions. In regions that agricultural abandonment and/or de-intensification have taken place the water quality conditions have improved. For nitrates, the decoupling of agricultural subsidies and the reduction of the amount of subsidies to irrigated land underlie the observed reduction of nitrates concentration. Measures of irrigation modernization and establishment of vulnerable zones to nitrates ameliorate the concentration of nitrates in subbasins showing an increasing trend. However, the effect of nitrates load from upland areas, intensification of biomass and crop prices present a greater weight leading to the final increasing trend in this subbasins group, where annual crops dominate. For suspended solids, there is no clear evidence that decoupling process have influenced negatively or positively. Nevertheless, greater values of subsidies still linked to production, particularly in irrigated regions, lead to increasing erosion rates. Although agricultural production has augmented in the basin and water efficiency in the agricultural sector has improved, the issue of high erosion rates has not yet been properly faced. The water footprint (WF) assessment reveals that for 1 L Spanish olive oil more than 99.5% of the WF is related to the olive fruit production, whereas less than 0.5% is due to other components i.e. bottle, cap and label. Over the studied period, the green WF in rainfed and irrigated systems represents about 72% and 12%, respectively, of the total WF. Blue and grey WFs represent 6% and 10%, respectively. The olive production is concentrated in regions with the smallest WF per unit of product. The olive oil production has increased its apparent water productivity from 1997 to 2008 incentivized by growing trade prices, but also did the amount of virtual water exports. In fact, the largest producing areas present high water use efficiency per product and apparent water productivity as well as less nitrates pollution potential, but this enhances the pressure on the available water resources. Increasing groundwater abstractions related to olive oil exports may add further pressure to the already stressed Guadalquivir basin. This shows the need to balance the market forces with the available local resources. Concerning the effects of LULCC on the Upper Turia basin’s streamflow, LULCC play a significant role on the water balance, but it is not the main driver underpinning the observed reduction on Turia's streamflow. Increasing mean temperature is the main factor supporting larger evapotranspiration rates and streamflow reduction. In fact, LULCC and climate change have had an offsetting effect on the streamflow generation during the study period. While streamflow has been negatively affected by increasing temperature, ongoing LULCC have positively compensated with reduced evapotranspiration rates, thanks to mainly shrubland clearing and forest degradation processes. The research provides insight for strengthening the interdisciplinarity between hydrological and spatial planning, highlighting the need to include the implications of LULCC in future hydrological plans. These findings are valuable for the management of the Turia river basin, as well as a useful approach for the determination of the weight of LULCC on the hydrological response in other regions.

Two-dimensional direct rainfall hydraulic models applied to direct calculation of hydrologic events in dams to prevent overtopping

One of the main concerns when conducting a dam test is the acute determination of the hydrograph for a specific flood event. The use of 2D direct rainfall hydraulic mathematical models on a finite elements mesh, combined with the efficiency of vector calculus that provides CUDA (Compute Unified Device Architecture) technology, enables nowadays the simulation of complex hydrological models without the need for terrain subbasin and transit splitting (as in HEC-HMS). Both the Spanish PNOA (National Plan of Aereal Orthophotography) Digital Terrain Model GRID with a 5 x 5 m accuracy and the CORINE GIS Land Cover (Coordination of INformation of the Environment) that allows assessment of the ground roughness, provide enough data to easily build these kind of models

Classification of Satellite Images by means of Fuzzy Rules generated by a Genetic Algorithm

The data acquired by Remote Sensing systems allow obtaining thematic maps of the earth's surface, by means of the registered image classification. This implies the identification and categorization of all pixels into land cover classes. Traditionally, methods based on statistical parameters have been widely used, although they show some disadvantages. Nevertheless, some authors indicate that those methods based on artificial intelligence, may be a good alternative. Thus, fuzzy classifiers, which are based on Fuzzy Logic, include additional information in the classification process through based-rule systems. In this work, we propose the use of a genetic algorithm (GA) to select the optimal and minimum set of fuzzy rules to classify remotely sensed images. Input information of GA has been obtained through the training space determined by two uncorrelated spectral bands (2D scatter diagrams), which has been irregularly divided by five linguistic terms defined in each band. The proposed methodology has been applied to Landsat-TM images and it has showed that this set of rules provides a higher accuracy level in the classification process

Modelización del cálculo de Indicadores Ambientales, a partir de sensores ópticos de alta resolución temporal a bordo de satélite

La presente tesis doctoral tiene por objeto el estudio y análisis de técnicas y modelos de obtención de parámetros biofísicos e indicadores ambientales, de manera automatizada a partir de imágenes procedentes de satélite de alta resolución temporal. En primer lugar se revisan los diferentes programas espaciales de observación del territorio, con especial atención a los que proporcionan dicha resolución. También se han revisado las metodologías y procesos que permiten la obtención de diferentes parámetros cuantitativos y documentos cualitativos, relacionados con diversos aspectos de las cubiertas terrestres, atendiendo a su adaptabilidad a las particularidades de los datos. En segundo lugar se propone un modelo de obtención de parámetros ambientales, que integra información proveniente de sensores espaciales y de otras fuentes auxiliares utilizando, en cierta medida, las metodologías presentadas en apartados anteriores y optimizando algunas de las referidas o proponiendo otras nuevas, de manera que se permita dicha obtención de manera eficiente, a partir de los datos disponibles y de forma sistemática. Tras esta revisión de metodologías y propuesta del modelo, se ha procedido a la realización de experimentos, con la finalidad de comprobar su comportamiento en diferentes casos prácticos, depurar los flujos de datos y procesos, así como establecer las situaciones que pueden afectar a los resultados. De todo ello se deducirá la evaluación del referido modelo. Los sensores considerados en este trabajo han sido MODIS, de alta resolución temporal y Thematic Mapper (TM), de media resolución espacial, por tratarse de instrumentos de referencia en la realización de estudios ambientales. También por la duración de sus correspondientes misiones de registro de datos, lo que permite realizar estudios de evolución temporal de ciertos parámetros biofísicos, durante amplios periodos de tiempo. Así mismo. es de destacar que la continuidad de los correspondientes programas parece estar asegurada. Entre los experimentos realizados, se ha ensayado una metodología para la integración de datos procedentes de ambos sensores. También se ha analizado un método de interpolación temporal que permite obtener imágenes sintéticas con la resolución espacial de TM (30 m) y la temporal de MODIS (1 día), ampliando el rango de aplicación de este último sensor. Asimismo, se han analizado algunos de los factores que afectan a los datos registrados, tal como la geometría de la toma de los mismos y los episodios de precipitación, los cuales alteran los resultados obtenidos. Por otro lado, se ha comprobado la validez del modelo propuesto en el estudio de fenómenos ambientales dinámicos, en concreto la contaminación orgánica de aguas embalsadas. Finalmente, se ha demostrado un buen comportamiento del modelo en todos los casos ensayados, así como su flexibilidad, lo que le permite adaptarse a nuevos orígenes de datos, o nuevas metodologías de cálculo. Abstract This thesis aims to the study and analysis of techniques and models, in order to obtain biophysical parameters and environmental indicators in an automated way, using high temporal resolution satellite data. Firstly we have reviewed the main Earth Observation Programs, paying attention to those that provide high temporal resolution. Also have reviewed the methodologies and process flow diagrams in order to obtain quantitative parameters and qualitative documents, relating to various aspects of land cover, according to their adaptability to the peculiarities of the data. In the next stage, a model which allows obtaining environmental parameters, has been proposed. This structure integrates information from space sensors and ancillary data sources, using the methodologies presented in previous sections that permits the parameters calculation in an efficient and automated way. After this review of methodologies and the proposal of the model, we proceeded to carry out experiments, in order to check the behavior of the structure in real situations. From this, we derive the accuracy of the model. The sensors used in this work have been MODIS, which is a high temporal resolution sensor, and Thematic Mapper (TM), which is a medium spatial resolution instrument. This choice was motivated because they are reference sensors in environmental studies, as well as for the duration of their corresponding missions of data logging, and whose continuity seems assured. Among the experiments, we tested a methodology that allows the integration of data from cited sensors, we discussed a proposal for a temporal interpolation method for obtaining synthetic images with spatial resolution of TM (30 m) and temporal of MODIS (1 day), extending the application range of this one. Furthermore, we have analyzed some of the factors that affect the recorded data, such as the relative position of the satellite with the ground point, and the rainfall events, which alter the obtained results. On the other hand, we have proven the validity of the proposed model in the study of the organic contamination in inland water bodies. Finally, we have demonstrated a good performance of the proposed model in all cases tested, as well as its flexibility and adaptability.

Comparison of fractal dimensions based on segmented NDVI fields obtained from different remote sensors.

Satellite image data have become an important source of information for monitoring vegetation and mapping land cover at several scales. Beside this, the distribution and phenology of vegetation is largely associated with climate, terrain characteristics and human activity. Various vegetation indices have been developed for qualitative and quantitative assessment of vegetation using remote spectral measurements. In particular, sensors with spectral bands in the red (RED) and near-infrared (NIR) lend themselves well to vegetation monitoring and based on them [(NIR - RED) / (NIR + RED)] Normalized Difference Vegetation Index (NDVI) has been widespread used. Given that the characteristics of spectral bands in RED and NIR vary distinctly from sensor to sensor, NDVI values based on data from different instruments will not be directly comparable. The spatial resolution also varies significantly between sensors, as well as within a given scene in the case of wide-angle and oblique sensors. As a result, NDVI values will vary according to combinations of the heterogeneity and scale of terrestrial surfaces and pixel footprint sizes. Therefore, the question arises as to the impact of differences in spectral and spatial resolutions on vegetation indices like the NDVI. The aim of this study is to establish a comparison between two different sensors in their NDVI values at different spatial resolutions.