Spatial and temporal variations of trace element distribution in soils and street dust of an industrial town in NW Spain: 15 years of study

Extensive spatial and temporal surveys, over 15 years, have been conducted in soil in urban parks and street dusts in one of the most polluted cities in western Europe, Avilés (NW Spain). The first survey was carried out in 1996, and since then monitoring has been undertaken every five years. Whilst the sampling site is a relatively small town, industrial activities (mainly the steel industry and Zn and Al metallurgy) and other less significant urban sources, such as traffic, strongly affect the load of heavy metals in the urban aerosol. Elemental tracers have been used to characterise the influence of these sources on the composition of soil and dust. Although PM10 has decreased over these years as a result of environmental measures undertaken in the city, some of the “industrial” elements still remain in concentrations of concern for example, up to 4.6% and 0.5% of Zn in dust and soil, respectively. Spatial trends in metals such as Zn and Cd clearly reflect sources from the processing industries. The concentrations of these elements across Europe have reduced over time, however the most recent results from Avilés revealed an upward trend in concentration for Zn, Cd, Hg and As. A risk assessment of the soil highlighted As as an element of concern since its cancer risk in adults was more than double the value above which regulatory agencies deem it to be unacceptable. If children were considered to be the receptors, then the risk nearly doubles from this element.

On the charged particle drag acting on LAGEOS

A recent study by the authors points to Charged Particle Drag (CPD) as a contributor to revisit in the LAGEOS non-gravitational perturbations problem. Such perturbations must account for dynamical contributions in the order of pms−2 . The simulated effect takes into account: (i) spatial and temporal variations of the plasmatic parameters (temperature and concentration of the species), (ii) spacecraft potential variations caused by both the eclipse passages and variations in the parameters mentioned above, and (iii) solar and geomagnetic conditions. Furthermore, recent theoretical improvements concerning scattering drag overcome previous limitations allowing for a complete formulation of this effect. For each satellite the lifetime CPD instantaneous acceleration is computed. The plasmatic parameters have been obtained fromthe Sheffield Coupled Thermosphere-Ionosphere-Plasmasphere (SCTIP) semi-empirical model (up to the polar region), as well as alytical/empirical approximations based on spacecraft measurements for the auroral and polar regions. Results show that maximum amplitudes for LAGEOSI are larger than those for LAGEOS-II: −85 pms−2 and −70 pms−2 respectively. This is due to the almost (magnetically) polar orbit configuration of the first, producing larger combinations of plasmatic parameter values. High solar activity has a huge impact in the resulting LAGEOS accelerations: it yields a perfect modulation of the resulting acceleration with maximum amplitudes up to a factor of 10 when comparing low and high activity periods. On the other hand, the impact of the geomagnetic activity results into a reduction of the effect itself, probably due to a decrease in the hydrogen concentration for high energy input periods. The acceleration results will be used in a refined orbit computation in a subsequent investigation.

Farm nitrogen balances in six European landscapes as an indicator for nitrogen losses and basis for improved management.

Improved management of nitrogen (N) in agriculture is necessary to achieve a sustainable balance between the production of food and other biomass, and the unwanted effects of N on water pollution, greenhouse gas emissions, biodiversity deterioration and human health. To analyse farm N-losses and the complex interactions within farming systems, efficient methods for identifying emissions hotspots and evaluating mitigation measures are therefore needed. The present paper aims to fill this gap at the farm and landscape scales. Six agricultural landscapes in Poland (PL), the Netherlands (NL), France (FR), Italy (IT), Scotland (UK) and Denmark (DK) were studied, and a common method was developed for undertaking farm inventories and the derivation of farm N balances, N surpluses and for evaluating uncertainty for the 222 farms and 11 440 ha of farmland included in the study. In all landscapes, a large variation in the farm N surplus was found, and thereby a large potential for reductions. The highest average N surpluses were found in the most livestock-intensive landscapes of IT, FR, and NL; on average 202 ± 28, 179 ± 63 and 178 ± 20 kg N ha−1 yr−1, respectively. All landscapes showed hotspots, especially from livestock farms, including a special UK case with large-scale landless poultry farming. Overall, the average N surplus from the land-based UK farms dominated by extensive sheep and cattle grazing was only 31 ± 10 kg N ha−1 yr−1, but was similar to the N surplus of PL and DK (122 ± 20 and 146 ± 55 kg N ha−1 yr−1, respectively) when landless poultry farming was included. We found farm N balances to be a useful indicator for N losses and the potential for improving N management. Significant correlations to N surplus were found, both with ammonia air concentrations and nitrate concentrations in soils and groundwater, measured during the period of N management data collection in the landscapes from 2007–2009. This indicates that farm N surpluses may be used as an independent dataset for validation of measured and modelled N emissions in agricultural landscapes. No significant correlation was found with N measured in surface waters, probably because of spatial and temporal variations in groundwater buffering and biogeochemical reactions affecting N flows from farm to surface waters. A case study of the development in N surplus from the landscape in DK from 1998–2008 showed a 22% reduction related to measures targeted at N emissions from livestock farms. Based on the large differences in N surplus between average N management farms and the most modern and N-efficient farms, it was concluded that additional N-surplus reductions of 25–50%, as compared to the present level, were realistic in all landscapes. The implemented N-surplus method was thus effective for comparing and synthesizing results on farm N emissions and the potentials of mitigation options. It is recommended for use in combination with other methods for the assessment of landscape N emissions and farm N efficiency, including more detailed N source and N sink hotspot mapping, measurements and modelling.

Desarrollo, implementación y análisis de un modelo cinemático local de velocidades para Costa Rica basado en la integración y homogenización de datos de estaciones GNSS continuas

El primer procesamiento estricto realizado con el software científico Bernese y contemplando las más estrictas normas de cálculo recomendadas internacionalmente, permitió obtener un campo puntual de alta exactitud, basado en la integración y estandarización de los datos de una red GPS ubicada en Costa Rica. Este procesamiento contempló un total de 119 semanas de datos diarios, es decir unos 2,3 años, desde enero del año 2009 hasta abril del año 2011, para un total de 30 estaciones GPS, de las cuales 22 están ubicadas en el territorio nacional de Costa Rica y 8 internaciones pertenecientes a la red del Sistema Geocéntrico para las Américas (SIRGAS). Las denominadas soluciones semilibres generaron, semana a semana, una red GPS con una alta exactitud interna definida por medio de los vectores entre las estaciones y las coordenadas finales de la constelación satelital. La evaluación semanal dada por la repetibilidad de las soluciones brindó en promedio errores de 1,7 mm, 1,4 mm y 5,1 mm en las componentes [n e u], confirmando una alta consistencia en estas soluciones. Aunque las soluciones semilibres poseen una alta exactitud interna, las mismas no son utilizables para fines de análisis cinemático, pues carecen de un marco de referencia. En Latinoamérica, la densificación del Marco Internacional Terrestre de Referencia (ITRF), está representado por la red de estaciones de operación continua GNSS de SIRGAS, denominada como SIRGAS-CON. Por medio de las denominadas coordenadas semanales finales de las 8 estaciones consideradas como vínculo, se refirió cada una de las 119 soluciones al marco SIRGAS. La introducción del marco de referencia SIRGAS a las soluciones semilibres produce deformaciones en estas soluciones. Las deformaciones de las soluciones semilibres son producto de las cinemática de cada una de las placas en las que se ubican las estaciones de vínculo. Luego de efectuado el amarre semanal a las coordenadas SIRGAS, se hizo una estimación de los vectores de velocidad de cada una de las estaciones, incluyendo las de amarre, cuyos valores de velocidad se conocen con una alta exactitud. Para la determinación de las velocidades de las estaciones costarricenses, se programó una rutina en ambiente MatLab, basada en una ajuste por mínimos cuadrados. Los valores obtenidos en el marco de este proyecto en comparación con los valores oficiales, brindaron diferencias promedio del orden de los 0,06 cm/a, -0,08 cm/a y -0,10 cm/a respectivamente para las coordenadas [X Y Z]. De esta manera se logró determinar las coordenadas geocéntricas [X Y Z]T y sus variaciones temporales [vX vY vZ]T para el conjunto de 22 estaciones GPS de Costa Rica, dentro del datum IGS05, época de referencia 2010,5. Aunque se logró una alta exactitud en los vectores de coordenadas geocéntricas de las 22 estaciones, para algunas de las estaciones el cálculo de las velocidades no fue representativo debido al relativo corto tiempo (menos de un año) de archivos de datos. Bajo esta premisa, se excluyeron las ocho estaciones ubicadas al sur de país. Esto implicó hacer una estimación del campo local de velocidades con solamente veinte estaciones nacionales más tres estaciones en Panamá y una en Nicaragua. El algoritmo usado fue el denominado Colocación por Mínimos Cuadrados, el cual permite la estimación o interpolación de datos a partir de datos efectivamente conocidos, el cual fue programado mediante una rutina en ambiente MatLab. El campo resultante se estimó con una resolución de 30' X 30' y es altamente constante, con una velocidad resultante promedio de 2,58 cm/a en una dirección de 40,8° en dirección noreste. Este campo fue validado con base en los datos del modelo VEMOS2009, recomendado por SIRGAS. Las diferencias de velocidad promedio para las estaciones usadas como insumo para el cálculo del campo fueron del orden los +0,63 cm/a y +0,22 cm/a para los valores de velocidad en latitud y longitud, lo que supone una buena determinación de los valores de velocidad y de la estimación de la función de covarianza empírica, necesaria para la aplicación del método de colocación. Además, la grilla usada como base para la interpolación brindó diferencias del orden de -0,62 cm/a y -0,12 cm/a para latitud y longitud. Adicionalmente los resultados de este trabajo fueron usados como insumo para hacer una aproximación en la definición del límite del llamado Bloque de Panamá dentro del territorio nacional de Costa Rica. El cálculo de las componentes del Polo de Euler por medio de una rutina programa en ambiente MatLab y aplicado a diferentes combinaciones de puntos no brindó mayores aportes a la definición física de este límite. La estrategia lo que confirmó fue simplemente la diferencia en la dirección de todos los vectores velocidad y no permitió reveló revelar con mayor detalle una ubicación de esta zona dentro del territorio nacional de Costa Rica. ABSTRACT The first strict processing performed with the Bernese scientific software and contemplating the highest standards internationally recommended calculation, yielded a precise field of high accuracy, based on the integration and standardization of data from a GPS network located in Costa Rica. This processing watched a total of 119 weeks of daily data, is about 2.3 years from January 2009 to April 2011, for a total of 30 GPS stations, of which 22 are located in the country of Costa Rica and 8 hospitalizations within the network of Geocentric System for the Americas (SIRGAS). The semi-free solutions generated, every week a GPS network with high internal accuracy defined by vectors between stations and the final coordinates of the satellite constellation. The weekly evaluation given by repeatability of the solutions provided in average errors of 1.7 mm 1.4 mm and 5.1 mm in the components [n e u], confirming a high consistency in these solutions. Although semi-free solutions have a high internal accuracy, they are not used for purposes of kinematic analysis, because they lack a reference frame. In Latin America, the densification of the International Terrestrial Reference Frame (ITRF), is represented by a network of continuously operating GNSS stations SIRGAS, known as SIRGAS-CON. Through weekly final coordinates of the 8 stations considered as a link, described each of the solutions to the frame 119 SIRGAS. The introduction of the frame SIRGAS to semi-free solutions generates deformations. The deformations of the semi-free solutions are products of the kinematics of each of the plates in which link stations are located. After SIRGAS weekly link to SIRGAS frame, an estimate of the velocity vectors of each of the stations was done. The velocity vectors for each SIRGAS stations are known with high accuracy. For this calculation routine in MatLab environment, based on a least squares fit was scheduled. The values obtained compared to the official values, gave average differences of the order of 0.06 cm/yr, -0.08 cm/yr and -0.10 cm/yr respectively for the coordinates [XYZ]. Thus was possible to determine the geocentric coordinates [XYZ]T and its temporal variations [vX vY vZ]T for the set of 22 GPS stations of Costa Rica, within IGS05 datum, reference epoch 2010.5. The high accuracy vector for geocentric coordinates was obtained, however for some stations the velocity vectors was not representative because of the relatively short time (less than one year) of data files. Under this premise, the eight stations located in the south of the country were excluded. This involved an estimate of the local velocity field with only twenty national stations plus three stations in Panama and Nicaragua. The algorithm used was Least Squares Collocation, which allows the estimation and interpolation of data from known data effectively. The algorithm was programmed with MatLab. The resulting field was estimated with a resolution of 30' X 30' and is highly consistent with a resulting average speed of 2.58 cm/y in a direction of 40.8° to the northeast. This field was validated based on the model data VEMOS2009 recommended by SIRGAS. The differences in average velocity for the stations used as input for the calculation of the field were of the order of +0.63 cm/yr, +0.22 cm/yr for the velocity values in latitude and longitude, which is a good determination velocity values and estimating the empirical covariance function necessary for implementing the method of application. Furthermore, the grid used as the basis for interpolation provided differences of about -0.62 cm/yr, -0.12 cm/yr to latitude and longitude. Additionally, the results of this investigation were used as input to an approach in defining the boundary of Panama called block within the country of Costa Rica. The calculation of the components of the Euler pole through a routine program in MatLab and applied to different combinations of points gave no further contributions to the physical definition of this limit. The strategy was simply confirming the difference in the direction of all the velocity vectors and not allowed to reveal more detail revealed a location of this area within the country of Costa Rica.

Automatic physical activity pattern recognition and meal intake modelling in type 1 diabetes for their inclusion in artificial pancreas systems

La diabetes comprende un conjunto de enfermedades metabólicas que se caracterizan por concentraciones de glucosa en sangre anormalmente altas. En el caso de la diabetes tipo 1 (T1D, por sus siglas en inglés), esta situación es debida a una ausencia total de secreción endógena de insulina, lo que impide a la mayoría de tejidos usar la glucosa. En tales circunstancias, se hace necesario el suministro exógeno de insulina para preservar la vida del paciente; no obstante, siempre con la precaución de evitar caídas agudas de la glucemia por debajo de los niveles recomendados de seguridad. Además de la administración de insulina, las ingestas y la actividad física son factores fundamentales que influyen en la homeostasis de la glucosa. En consecuencia, una gestión apropiada de la T1D debería incorporar estos dos fenómenos fisiológicos, en base a una identificación y un modelado apropiado de los mismos y de sus sorrespondientes efectos en el balance glucosa-insulina. En particular, los sistemas de páncreas artificial –ideados para llevar a cabo un control automático de los niveles de glucemia del paciente– podrían beneficiarse de la integración de esta clase de información. La primera parte de esta tesis doctoral cubre la caracterización del efecto agudo de la actividad física en los perfiles de glucosa. Con este objetivo se ha llevado a cabo una revisión sistemática de la literatura y meta-análisis que determinen las respuestas ante varias modalidades de ejercicio para pacientes con T1D, abordando esta caracterización mediante unas magnitudes que cuantifican las tasas de cambio en la glucemia a lo largo del tiempo. Por otro lado, una identificación fiable de los periodos con actividad física es un requisito imprescindible para poder proveer de esa información a los sistemas de páncreas artificial en condiciones libres y ambulatorias. Por esta razón, la segunda parte de esta tesis está enfocada a la propuesta y evaluación de un sistema automático diseñado para reconocer periodos de actividad física, clasificando su nivel de intensidad (ligera, moderada o vigorosa); así como, en el caso de periodos vigorosos, identificando también la modalidad de ejercicio (aeróbica, mixta o de fuerza). En este sentido, ambos aspectos tienen una influencia específica en el mecanismo metabólico que suministra la energía para llevar a cabo el ejercicio y, por tanto, en las respuestas glucémicas en T1D. En este trabajo se aplican varias combinaciones de técnicas de aprendizaje máquina y reconocimiento de patrones sobre la fusión multimodal de señales de acelerometría y ritmo cardíaco, las cuales describen tanto aspectos mecánicos del movimiento como la respuesta fisiológica del sistema cardiovascular ante el ejercicio. Después del reconocimiento de patrones se incorpora también un módulo de filtrado temporal para sacar partido a la considerable coherencia temporal presente en los datos, una redundancia que se origina en el hecho de que en la práctica, las tendencias en cuanto a actividad física suelen mantenerse estables a lo largo de cierto tiempo, sin fluctuaciones rápidas y repetitivas. El tercer bloque de esta tesis doctoral aborda el tema de las ingestas en el ámbito de la T1D. En concreto, se propone una serie de modelos compartimentales y se evalúan éstos en función de su capacidad para describir matemáticamente el efecto remoto de las concetraciones plasmáticas de insulina exógena sobre las tasas de eleiminación de la glucosa atribuible a la ingesta; un aspecto hasta ahora no incorporado en los principales modelos de paciente para T1D existentes en la literatura. Los datos aquí utilizados se obtuvieron gracias a un experimento realizado por el Institute of Metabolic Science (Universidad de Cambridge, Reino Unido) con 16 pacientes jóvenes. En el experimento, de tipo ‘clamp’ con objetivo variable, se replicaron los perfiles individuales de glucosa, según lo observado durante una visita preliminar tras la ingesta de una cena con o bien alta carga glucémica, o bien baja. Los seis modelos mecanísticos evaluados constaban de: a) submodelos de doble compartimento para las masas de trazadores de glucosa, b) un submodelo de único compartimento para reflejar el efecto remoto de la insulina, c) dos tipos de activación de este mismo efecto remoto (bien lineal, bien con un punto de corte), y d) diversas condiciones iniciales. ABSTRACT Diabetes encompasses a series of metabolic diseases characterized by abnormally high blood glucose concentrations. In the case of type 1 diabetes (T1D), this situation is caused by a total absence of endogenous insulin secretion, which impedes the use of glucose by most tissues. In these circumstances, exogenous insulin supplies are necessary to maintain patient’s life; although caution is always needed to avoid acute decays in glycaemia below safe levels. In addition to insulin administrations, meal intakes and physical activity are fundamental factors influencing glucose homoeostasis. Consequently, a successful management of T1D should incorporate these two physiological phenomena, based on an appropriate identification and modelling of these events and their corresponding effect on the glucose-insulin balance. In particular, artificial pancreas systems –designed to perform an automated control of patient’s glycaemia levels– may benefit from the integration of this type of information. The first part of this PhD thesis covers the characterization of the acute effect of physical activity on glucose profiles. With this aim, a systematic review of literature and metaanalyses are conduced to determine responses to various exercise modalities in patients with T1D, assessed via rates-of-change magnitudes to quantify temporal variations in glycaemia. On the other hand, a reliable identification of physical activity periods is an essential prerequisite to feed artificial pancreas systems with information concerning exercise in ambulatory, free-living conditions. For this reason, the second part of this thesis focuses on the proposal and evaluation of an automatic system devised to recognize physical activity, classifying its intensity level (light, moderate or vigorous) and for vigorous periods, identifying also its exercise modality (aerobic, mixed or resistance); since both aspects have a distinctive influence on the predominant metabolic pathway involved in fuelling exercise, and therefore, in the glycaemic responses in T1D. Various combinations of machine learning and pattern recognition techniques are applied on the fusion of multi-modal signal sources, namely: accelerometry and heart rate measurements, which describe both mechanical aspects of movement and the physiological response of the cardiovascular system to exercise. An additional temporal filtering module is incorporated after recognition in order to exploit the considerable temporal coherence (i.e. redundancy) present in data, which stems from the fact that in practice, physical activity trends are often maintained stable along time, instead of fluctuating rapid and repeatedly. The third block of this PhD thesis addresses meal intakes in the context of T1D. In particular, a number of compartmental models are proposed and compared in terms of their ability to describe mathematically the remote effect of exogenous plasma insulin concentrations on the disposal rates of meal-attributable glucose, an aspect which had not yet been incorporated to the prevailing T1D patient models in literature. Data were acquired in an experiment conduced at the Institute of Metabolic Science (University of Cambridge, UK) on 16 young patients. A variable-target glucose clamp replicated their individual glucose profiles, observed during a preliminary visit after ingesting either a high glycaemic-load or a low glycaemic-load evening meal. The six mechanistic models under evaluation here comprised: a) two-compartmental submodels for glucose tracer masses, b) a single-compartmental submodel for insulin’s remote effect, c) two types of activations for this remote effect (either linear or with a ‘cut-off’ point), and d) diverse forms of initial conditions.

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.