The Extended Water Footprint and illegal groundwater use in the Upper Guadiana Basin (Spain)

Tablas de Daimiel National Park is located in the Upper Guadiana Basin and represents one of the largest and most important wetlands in Europe. The long term ecological integrity of this wetland is inherently associated with the maintenance of a shallow groundwater table, namely the Western Mancha aquifer (WMA) or Aquifer 23. The intensive use of groundwater, mainly for irrigation, has led over the last decades to deep socio‐economic changes. Such intensive use has also lowered the water table of Aquifer 23, drastically reducing the flooded area of the wetland and threatening its ecological integrity. A number of plans and measures have been developed and implemented since the declaration of overexploitation of Aquifer 23 in the year 1987. The most recent one is the Special Plan for the Upper Guadiana (SPUG), approved in 2008. This Plan is the main measure to comply with achieving the objective of good quantitative and qualitative status required under the Water Framework Directive (2000). This paper offers a new type of integrated analysis which allows assessing under a common lens the physical, economic and social dimensions of groundwater use in the area. The first objective is to calculate the groundwater footprint of agricultural production in the Upper Guadiana basin and its evolution during 2000‐2008. For this purpose, we have applied the Extended Water Footprint (EWF) methodology ‐a novel approach based on the classical Water Footprint (WF) approach‐ that includes an assessment of the water productivity from an economic and social perspective. Compared to the classical WF, the EWF allows for a more complete overview of the sector, providing new insights for policy decisions (e.g. to define options and possibilities on water re‐allocation in order to achieve both better ecosystem conservation and social equity). The second objective is to use the EWF to compare the existing authorized and non‐authorized or illegal use of water. This allows us to discuss current initiatives by public authorities in relation to the existing frame of water rights

Investigaciones dendrogeomorfológicas aplicadas al estudio de la peligrosidad por avenidas e inundaciones en el parque nacional de la Caldera de Taburiente (La Palma, Islas Canarias).

Las Angustias torrent is an ungauged stream located in the Caldera de Taburiente National Park (Canary Islands), where frequent and intense flood events occur (even causing fatalities, such as in November 2001). The aim of this research is to analyse the flood hazard in one of the most visited areas of this protected area. To carry out this objective, during two dendrogeomorphological campaigns, all pine trees located on the stream bed and banks showing external evidence of flash flood damage were sampled. In addition, a detailed topographical survey using GPS and total station allowed us to obtain the inputs for hydraulic modelling. In the preliminary results, two flood events have been identified (1967-1968 and 1999- 2000), and evidence of other dates has been observed (1976, 1983, 1993 and 2001). Magnitude reconstruction of these events will improve flood hazard and risk analyses and will be useful for National Park managers.

Geomorphological analysis and evolution of an altered floodplain drainage system. The case of the Partido Stream (Spain)

The Partido Stream is a small torrential course that flows into the marsh of the Doñana National Park, an area that was declared a World Heritage Site in 1994. Before 1981, floods occurred, and the stream overflowed onto a floodplain. As an old alluvial fan, the floodplain has its singular orography and functionality. Fromthe floodplain, several drainage channels, locally called caño, discharged into themarsh. The Partido Streamhad themorphology of a caño and covered approximately 8 km from the old fan to the marsh. The stream was straightened and channelised in 1981 to cultivate the old fan. This resulted in floods that were concentrated between the banks in the following years, which caused the depth of water and the shear stress to increase, thus, scouring the river bed and river banks. In this case, the eroded materials were carried towards the marsh where a new alluvial fan evolved. Control measures on the old fan were implemented in 2006 to stop the development of the new alluvial fan downstream over the marsh. Thus, the stream would partially recover its original behaviour that it had before channelisation, moving forwards in a new, balanced state. The present study describes the geomorphological evolution that channelisation has caused since 1981 and the later slow process of recovery of the original hydraulic-sedimentation regime since 2006. Additionally, it deepens the understanding of the original hydraulic behaviour of the stream, combining field data and 2D simulations.

Evaluación del ruido ambiental en espacios naturales protegidos: implicaciones para su gestión

Esta Tesis Doctoral trata sobre la caracterización acústica de los ecosistemas naturales y la evaluación del impacto ambiental del ruido antropogénico sobre sus potenciales receptores en estos lugares, incluidos los receptores no humanos y sus efectos ecológicos, además, analiza las implicaciones para su gestión a distintas escalas y se lleva a cabo una valoración económica. Este trabajo ofrece soluciones para caracterizar los paisajes sonoros de forma compatible con distintas escalas de trabajo, nivel de esfuerzo técnico y en contextos de recursos limitados que haga viable su tratamiento como cualquier otra variable ambiental en el ámbito de la conservación y gestión del medio natural. Se han adaptado herramientas y metodologías propias de disciplinas como la acústica ambiental, bioacústica y ecología del paisaje, para servir a los objetivos específicos de la evaluación y gestión de los paisajes sonoros y el ruido ambiental en amplias extensiones geográficas. Se ha establecido un método general de muestreo sistemático para trabajo de campo y también se han adaptado métodos de modelización informática, que permiten analizar escenarios sonoros dinámicos en el tiempo y en el espacio, desde localizaciones puntuales hasta la escala del paisaje. Es posible elaborar cartografía ambiental con esta información y se ha representado gráficamente la zona de influencia de distintas fuentes de ruido sobre la calidad de distintos hábitats faunísticos. Se recomienda el uso del indicador del nivel de presión sonora equivalente (Leq) por su operatividad en medición y modelización, y su adaptabilidad a cualquier dimensión espacial y temporal que se requiera, por ejemplo en función del paisaje, actividades o especies que se establezcan como objeto de análisis. Se ha comprobado que las voces y conversaciones de parte de los excursionistas en zonas de reposo, observación y descanso (Laguna Grande de Peñalara) es la fuente de ruido que con mayor frecuencia identifican los propios visitantes (51%) y causa un incremento del nivel de presión sonora equivalente de unos 4,5 dBA sobre el nivel correspondiente al ambiente natural (Lnat). También se ha comprobado que carreteras con bajo nivel de tráfico (IMD<1000) pueden causar estrés fisiológico sobre la fauna y afectar a la calidad de sus hábitats. La isófona de 30 dBA del índice Leq (24h) permite dividir a los corzos de la zona de estudio en dos grupos con diferente nivel de estrés fisiológico, más elevado en los que se sitúan más cerca de la carretera con mayor volumen de tráfico y se expone a mayores niveles de ruido. Por otro lado, ha sido posible delimitar una zona de exclusión para la nidificación de buitre negro alrededor de las carreteras, coincidente con la isófona Leq (24h) de 40 dBA que afecta al 11% de su hábitat potencial. Además se ha llevado a cabo una novedosa valoración económica de la contaminación acústica en espacios naturales protegidos, mediante el análisis de la experiencia sonora de los visitantes del antiguo Parque Natural de Peñalara, y se ha constatado su disposición al pago de una entrada de acceso a estos lugares (aproximadamente 1 euro) si redundara en una mejora de su estado de conservación. En conclusión, los espacios naturales protegidos pueden sufrir un impacto ambiental significativo causado por fuentes de ruido localizadas en su interior pero también lejanas a ellos, que se sitúan fuera del ámbito de competencias de sus gestores. Sucesos sonoros como el sobrevuelo de aviones pueden incrementar en aproximadamente 8 dBA el nivel de referencia Lnat en las zonas tranquilas del parque. Se recomienda llevar a cabo una gestión activa del medio ambiente sonoro y se considera necesario extender la investigación sobre los efectos ecológicos del ruido ambiental a otros lugares y especies animales. ABSTRACT This PhD Thesis deals with acoustic characterization of natural ecosystems and anthropogenic noise impact assessment on potential receivers, including non-human receivers and their ecological effects. Besides, its management implications at different scales are analyzed and an economic valuation is performed. This study provides solutions for characterizing soundscapes in a compatible way with different working scales, level of technical effort and in a context of limited resources, so its treatment becomes feasible as for any other environmental variable in conservation and environmental management. Several tools and methodologies have been adapted from a variety of disciplines such as environmental acoustics, bioacoustics and landscape ecology, to better serve the specific goals of assessing and managing soundscapes and environmental noise in large areas. A procedure has been established for systematic field measurement surveys and noise common computer modelling methods have also been adapted in order to analyze dynamic soundscapes across time and space, from local to landscape scales. It is possible to create specific thematic cartography as for instance delimiting potential influence zone from different noise sources on animal habitats quality. Use of equivalent continuous sound pressure level index (Leq) is recommended because it provides great flexibility in operation for noise measurement and modelling, and because of its adaptability to any required temporal and spatial dimension, for instance landscape, activities or the target species established as study subjects. It has been found that human voices and conversations in a resting and contemplation area (Laguna Grande de Peñalara) is the most frequently referred noise source by national park visitors (51 %) when asked. Human voices alter this recreational area by increasing the sound pressure level approximately 4.5 dBA over the natural ambient level (Lnat). It has also been found that low traffic roads (AADT<1000 ) may cause physiological stress on wildlife and affect the quality of their habitats. It has also been possible to define a road-effect zone by noise mapping, which suggests an effective habitat loss within the Leq (24h) 30 dBA isophone in case of Roe deer and also divide the study area in two groups with different physiological stress level, higher for those exposed to higher noise levels and traffic volume. On the other hand, it has been possible to determine an exclusion area for Cinereous vulture nesting surrounding roads which is coincident with the Leq (24h) 40 dBA isophone and affects 11 % of the vulture potential habitat. It has also been performed an economic estimation of noise pollution impact on visitors’ perception and results showed that visitors would be willing to pay an entrance fee of approximately 1 euro if such payment is really bringing an improvement of the conservation status. In conclusion, protected areas may be significantly affected by anthropogenic noise sources located within the park borders but perturbations may also be caused by large-distance noise sources outside the park managers’ jurisdiction. Aircraft overflight events disrupted quietness and caused Leq increases of almost 8 dBA during a monitoring period with respect to Lnat reference levels in the park quiet areas. It is recommended to actively manage the acoustic environment. Finally, further research on ecological impacts of environmental noise needs to be extended to other species and places.

Alteración del comportamiento hidráulico natural de la vertiente noreste de la marisma de El Rocío por efectos antrópicos

El geo-dinamismo torrencial es resultado de la incidencia de precipitaciones torrenciales en la cuenca hidrográfica, y se manifiesta por una intensificación del ciclo del agua (con la presencia de mayores volúmenes de escorrentía y la consiguiente aparición de repentinos y elevados caudales de avenida) e intensificación del ciclo de los sedimentos (mayor erosión del suelo, incremento del transporte de sedimentos en los cauces de drenaje y la formación de grades depósitos de sedimentos en el área dominada de la cuenca). Estos procesos naturales se pueden agravar si se realizan en la cuenca actuaciones que alteren sustancialmente su morfología; es el caso del arroyo del Partido, tributario a la marisma de El Rocío en el Parque Nacional de Doñana, cuyo encauzamiento del tramo final sobre su antiguo cono de sedimentación fue el causante de la formación de un nuevo cono de sedimentación sobre la marisma. Pero las intervenciones antrópicas también pueden dirigir el sistema hacia la estabilidad sedimentaria y el tramo final del arroyo del Partido ha experimentado asimismo actuaciones de carácter restaurador. El antiguo cono de sedimentación del arroyo del Partido es un lugar complejo en el que drenan también otros tres cursos que intercambian sus flujos con dicho arroyo durante las avenidas. Esta Memoria analiza cuatro situaciones del paraje objeto de estudio: la primera responde al estado original del lugar en cuestión (situación A); la segunda a la creada tras el encauzamiento del arroyo del Partido en 1981 (situación B); la tercera a la formada pasado un tiempo desde el encauzamiento, cuando el geo-dinamismo del lugar hizo sus efectos sobre el paraje en cuestión (situación C); y la cuarta la generada tras la restauración en el verano de 2006 del régimen hidráulico-sedimentario del arroyo del Partido con la Actuación núm. 3 del Proyecto Doñana 2005. Se han investigado los efectos que la alteración de la morfología del sistema de drenaje del antiguo cono de sedimentación del arroyo del Partido causa en los parámetros hidráulicos que desencadenan el geo-dinamismo torrencial. Para ello se han reproducido las condiciones morfológicas de las cuatro situaciones más representativas ya comentadas que ha experimentado el paraje en los últimos 30 años y que han condicionado el comportamiento de las avenidas torrenciales ocurridas en dicho periodo. En este contexto se ha elaborado un modelo hidrológico con el programa HECHMS a partir de la información disponible y un modelo hidrodinámico en 2D con el programa IBER; con este último se han comparado ocho escenarios característicos derivados de la incidencia de dos avenidas representativas de los eventos torrenciales de la zona sobre las morfologías contempladas. También se ha estudiado la tensión de arrastre de las avenidas dada su incidencia en el terreno sobre el que transitan, analizando de manera separada los efectos en los cuatro cauces de drenaje del antiguo cono de sedimentación y del resto de la superficie del mismo. La investigación realizada ha permitido describir minuciosamente el comportamiento de los drenajes locales del antiguo y complejo cono de sedimentación que bordea el noroeste de la marisma del Parque Nacional de Doñana, desde la marisma de El Rocío hasta la marisma gallega, lo que redundará en beneficio de la gestión del Parque Nacional de Doñana y de su entorno, y especialmente en la seguridad del perímetro oriental de la aldea de El Rocío. ABSTRACT The torrential geo-dynamism results from the effect of torrential precipitations on the hydrological watershed, and it shows ups with the intensification of the water cycle (with presence of higher volumes of run-off and the consequent appearance of sudden and high food flows) and intensification of the sediment cycle (higher soil erosion, increase of sediment transport in the drainage courses and the formation of big sediment deposits in the area downstream the watershed). These natural processes can worsen if actions that change the morphology of the watershed are carried out. It is the case of the Partido stream, tributary of the marsh of El Rocío in the Doñana National Park, whose channelization of the final stretch over its old alluvial fan was the cause of the formation of a new alluvial fan over the marsh. But human actions can also go toward the sediment stabilization, and the final stretch of the Partido stream has experienced restoration actions as well. The old alluvial fan of the Partido stream is a complex terrain where also other three courses exchange their flows during floods. This document analyses four situations of the area studied: the first corresponds to the original state of the place (situation A); the second is the state created after the channel of the Partido stream was built in 1981 (situation B); the third is the state formed after some time passed and the geo-dynamism caused effect in the channelization (situation C); and the forth situation is the one created after the restoration in summer of 2006 of the hydraulic-sedimentary regime of the Partido stream with de Action number 3 of the Doñana 2005 Project. It has been researched the effects that the alteration of the morphology of the drainage system of the old alluvial fan of the Partido stream causes to the hydraulic parameters that trigger the torrential geo-dynamism. For that, the morphologic conditions of the four more representative situations were reproduced. Those situations experienced in the area in the last 30 years have conditioned the behaviour of torrential floods in that period. In this context, a hydrologic model with HEC-HMS and a 2D hydro-dynamic model with IBER have been set up. With IBER, eight scenes have been compared resulting from the impact of two representative floods from torrential events in the region against the morphologies considered. Also, the shear stress of the floods on the terrain has been studied, analysing separately the effect in the four drainage courses of the old alluvial fan and in the rest of its surfaces. The research carried out has allowed describing in minute detail the behavior of the local drainages of the old and complex alluvial fan that border the northwest part of the marsh of Doñana National Park, from the marsh of El Rocío to the gallega marsh. This will result in benefit for the management of Doñana National Park and its surroundings, and especially for the safety of the eastern limit of El Rocío village.

Unravelling past flash flood activity in a forested mountain catchment of the Spanish Central System

Flash floods represent one of the most common natural hazards in mountain catchments, and are frequent in Mediterranean environments. As a result of the widespread lack of reliable data on past events, the understanding of their spatio-temporal occurrence and their climatic triggers remains rather limited. Here, we present a dendrogeomorphic reconstruction of past flash flood activity in the Arroyo de los Puentes stream (Sierra de Guadarrama, Spanish Central System). We analyze a total of 287 increment cores from 178 disturbed Scots pine trees (Pinus sylvestris L.) which yielded indications on 212 growth disturbances related to past flash flood impact. In combination with local archives, meteorological data, annual forest management records and highly-resolved terrestrial data (i.e., LiDAR data and aerial imagery), the dendrogeomorphic time series allowed dating 25 flash floods over the last three centuries, with a major event leaving an intense geomorphic footprint throughout the catchment in 1936. The analysis of meteorological records suggests that the rainfall thresholds of flash floods vary with the seasonality of events. Dated flash floods in the 20th century were primarily related with synoptic troughs owing to the arrival of air masses from north and west on the Iberian Peninsula during negative indices of the North Atlantic Oscillation. The results of this study contribute considerably to a better understanding of hazards related with hydrogeomorphic processes in central Spain in general and in the Sierra de Guadarrama National Park in particular.

Planning and management of mobility in natural protected areas

The final purpose is the identification of: a) the effects of various choices in transport planning, both at long term and strategic level; b) the most effective policies of mobility management. The preliminary work was articulated in the following steps: 1) definition of protected area on the basis of ecological and socio-economic criteria and legislative constraints; 2) analysis of mobility needs in the protected areas; 3) reconstruction of the state of the art of mobility management in natural parks at European level; 4) analysis of used traffic flows measurement methods; 5) analysis of environmental impacts due to transport systems modelling (limited to air pollution and noise); 6) identification of mitigation measures to the potentially applied. The whole methodology has been firstly tested on the case study of the National Park of ?Gran Sasso and Monti della Laga? and further validated on the National Park of ?Gargano?, both located Italy: i) the concerned area has been zoned according to the land-use peculiarities; ii) the local situations of transport infrastructure (roads and parking), services (public transport systems) and rules (traffic regulations) have been mapped with references to physical and functional attributes; iii) the mobility, both systematic and touristic, has been synthetically represented in an origin-destination matrix. By means of an assignment model it has been determined the distribution of flows and the corresponding average speeds to quantify gaseous and noise emissions. On this basis the environmental criticalities in the reference scenario have been highlighted, as well as some alternative scenarios including both operational and infrastructural measures have been identified. The comparison between the projects and the reference scenario allowed the quantification of the effects (variation of emissions) for each scenario and a selection of the most effective management actions to be taken.

Spatially-Explicit Testing of a General Aboveground Carbon Density Estimation Model in aWestern Amazonian Forest Using Airborne LiDAR

Mapping aboveground carbon density in tropical forests can support CO2 emissionmonitoring and provide benefits for national resource management. Although LiDAR technology has been shown to be useful for assessing carbon density patterns, the accuracy and generality of calibrations of LiDAR-based aboveground carbon density (ACD) predictions with those obtained from field inventory techniques should be intensified in order to advance tropical forest carbon mapping. Here we present results from the application of a general ACD estimation model applied with small-footprint LiDAR data and field-based estimates of a 50-ha forest plot in Ecuador?s Yasuní National Park. Subplots used for calibration and validation of the general LiDAR equation were selected based on analysis of topographic position and spatial distribution of aboveground carbon stocks. The results showed that stratification of plot locations based on topography can improve the calibration and application of ACD estimation using airborne LiDAR (R2 = 0.94, RMSE = 5.81 Mg?C? ha?1, BIAS = 0.59). These results strongly suggest that a general LiDAR-based approach can be used for mapping aboveground carbon stocks in western lowland Amazonian forests.

LiDAR remote sensing applied to forest resources assessment

Disponer de información precisa y actualizada de inventario forestal es una pieza clave para mejorar la gestión forestal sostenible y para proponer y evaluar políticas de conservación de bosques que permitan la reducción de emisiones de carbono debidas a la deforestación y degradación forestal (REDD). En este sentido, la tecnología LiDAR ha demostrado ser una herramienta perfecta para caracterizar y estimar de forma continua y en áreas extensas la estructura del bosque y las principales variables de inventario forestal. Variables como la biomasa, el número de pies, el volumen de madera, la altura dominante, el diámetro o la altura media son estimadas con una calidad comparable a los inventarios tradicionales de campo. La presente tesis se centra en analizar la aplicación de los denominados métodos de masa de inventario forestal con datos LIDAR bajo diferentes condiciones y características de masa forestal (bosque templados puros y mixtos) y utilizando diferentes bases de datos LiDAR (información proveniente de vuelo nacionales e información capturada de forma específica). Como consecuencia de lo anterior, se profundiza en la generación de inventarios forestales continuos con LiDAR en grandes áreas. Los métodos de masa se basan en la búsqueda de relaciones estadísticas entre variables predictoras derivadas de la nube de puntos LiDAR y las variables de inventario forestal medidas en campo con el objeto de generar una cartografía continua de inventario forestal. El rápido desarrollo de esta tecnología en los últimos años ha llevado a muchos países a implantar programas nacionales de captura de información LiDAR aerotransportada. Estos vuelos nacionales no están pensados ni diseñados para fines forestales por lo que es necesaria la evaluación de la validez de esta información LiDAR para la descripción de la estructura del bosque y la medición de variables forestales. Esta información podría suponer una drástica reducción de costes en la generación de información continua de alta resolución de inventario forestal. En el capítulo 2 se evalúa la estimación de variables forestales a partir de la información LiDAR capturada en el marco del Plan Nacional de Ortofotografía Aérea (PNOA-LiDAR) en España. Para ello se compara un vuelo específico diseñado para inventario forestal con la información de la misma zona capturada dentro del PNOA-LiDAR. El caso de estudio muestra cómo el ángulo de escaneo, la pendiente y orientación del terreno afectan de forma estadísticamente significativa, aunque con pequeñas diferencias, a la estimación de biomasa y variables de estructura forestal derivadas del LiDAR. La cobertura de copas resultó más afectada por estos factores que los percentiles de alturas. Considerando toda la zona de estudio, la estimación de la biomasa con ambas bases de datos no presentó diferencias estadísticamente significativas. Las simulaciones realizadas muestran que las diferencias medias en la estimación de biomasa entre un vuelo específico y el vuelo nacional podrán superar el 4% en áreas abruptas, con ángulos de escaneo altos y cuando la pendiente de la ladera no esté orientada hacia la línea de escaneo. En el capítulo 3 se desarrolla un estudio en masas mixtas y puras de pino silvestre y haya, con un enfoque multi-fuente empleando toda la información disponible (vuelos LiDAR nacionales de baja densidad de puntos, imágenes satelitales Landsat y parcelas permanentes del inventario forestal nacional español). Se concluye que este enfoque multi-fuente es adecuado para realizar inventarios forestales continuos de alta resolución en grandes superficies. Los errores obtenidos en la fase de ajuste y de validación de los modelos de área basimétrica y volumen son similares a los registrados por otros autores (usando un vuelo específico y parcelas de campo específicas). Se observan errores mayores en la variable número de pies que los encontrados en la literatura, que pueden ser explicados por la influencia de la metodología de parcelas de radio variable en esta variable. En los capítulos 4 y 5 se evalúan los métodos de masa para estimar biomasa y densidad de carbono en bosques tropicales. Para ello se trabaja con datos del Parque Nacional Volcán Poás (Costa Rica) en dos situaciones diferentes: i) se dispone de una cobertura completa LiDAR del área de estudio (capitulo 4) y ii) la cobertura LiDAR completa no es técnica o económicamente posible y se combina una cobertura incompleta de LiDAR con imágenes Landsat e información auxiliar para la estimación de biomasa y carbono (capitulo 5). En el capítulo 4 se valida un modelo LiDAR general de estimación de biomasa aérea en bosques tropicales y se compara con los resultados obtenidos con un modelo ajustado de forma específica para el área de estudio. Ambos modelos están basados en la variable altura media de copas (TCH por sus siglas en inglés) derivada del modelo digital LiDAR de altura de la vegetación. Los resultados en el área de estudio muestran que el modelo general es una alternativa fiable al ajuste de modelos específicos y que la biomasa aérea puede ser estimada en una nueva zona midiendo en campo únicamente la variable área basimétrica (BA). Para mejorar la aplicación de esta metodología es necesario definir en futuros trabajos procedimientos adecuados de medición de la variable área basimétrica en campo (localización, tamaño y forma de las parcelas de campo). La relación entre la altura media de copas del LiDAR y el área basimétrica (Coeficiente de Stock) obtenida en el área de estudio varía localmente. Por tanto es necesario contar con más información de campo para caracterizar la variabilidad del Coeficiente de Stock entre zonas de vida y si estrategias como la estratificación pueden reducir los errores en la estimación de biomasa y carbono en bosques tropicales. En el capítulo 5 se concluye que la combinación de una muestra sistemática de información LiDAR con una cobertura completa de imagen satelital de moderada resolución (e información auxiliar) es una alternativa efectiva para la realización de inventarios continuos en bosques tropicales. Esta metodología permite estimar altura de la vegetación, biomasa y carbono en grandes zonas donde la captura de una cobertura completa de LiDAR y la realización de un gran volumen de trabajo de campo es económica o/y técnicamente inviable. Las alternativas examinadas para la predicción de biomasa a partir de imágenes Landsat muestran una ligera disminución del coeficiente de determinación y un pequeño aumento del RMSE cuando la cobertura de LiDAR es reducida de forma considerable. Los resultados indican que la altura de la vegetación, la biomasa y la densidad de carbono pueden ser estimadas en bosques tropicales de forma adecuada usando coberturas de LIDAR bajas (entre el 5% y el 20% del área de estudio). ABSTRACT The availability of accurate and updated forest data is essential for improving sustainable forest management, promoting forest conservation policies and reducing carbon emissions from deforestation and forest degradation (REDD). In this sense, LiDAR technology proves to be a clear-cut tool for characterizing forest structure in large areas and assessing main forest-stand variables. Forest variables such as biomass, stem volume, basal area, mean diameter, mean height, dominant height, and stem number can be thus predicted with better or comparable quality than with costly traditional field inventories. In this thesis, it is analysed the potential of LiDAR technology for the estimation of plot-level forest variables under a range of conditions (conifer & broadleaf temperate forests and tropical forests) and different LiDAR capture characteristics (nationwide LiDAR information vs. specific forest LiDAR data). This study evaluates the application of LiDAR-based plot-level methods in large areas. These methods are based on statistical relationships between predictor variables (derived from airborne data) and field-measured variables to generate wall to wall forest inventories. The fast development of this technology in recent years has led to an increasing availability of national LiDAR datasets, usually developed for multiple purposes throughout an expanding number of countries and regions. The evaluation of the validity of nationwide LiDAR databases (not designed specifically for forest purposes) is needed and presents a great opportunity for substantially reducing the costs of forest inventories. In chapter 2, the suitability of Spanish nationwide LiDAR flight (PNOA) to estimate forest variables is analyzed and compared to a specifically forest designed LiDAR flight. This study case shows that scan angle, terrain slope and aspect significantly affect the assessment of most of the LiDAR-derived forest variables and biomass estimation. Especially, the estimation of canopy cover is more affected than height percentiles. Considering the entire study area, biomass estimations from both databases do not show significant differences. Simulations show that differences in biomass could be larger (more than 4%) only in particular situations, such as steep areas when the slopes are non-oriented towards the scan lines and the scan angles are larger than 15º. In chapter 3, a multi-source approach is developed, integrating available databases such as nationwide LiDAR flights, Landsat imagery and permanent field plots from SNFI, with good resultos in the generation of wall to wall forest inventories. Volume and basal area errors are similar to those obtained by other authors (using specific LiDAR flights and field plots) for the same species. Errors in the estimation of stem number are larger than literature values as a consequence of the great influence that variable-radius plots, as used in SNFI, have on this variable. In chapters 4 and 5 wall to wall plot-level methodologies to estimate aboveground biomass and carbon density in tropical forest are evaluated. The study area is located in the Poas Volcano National Park (Costa Rica) and two different situations are analyzed: i) available complete LiDAR coverage (chapter 4) and ii) a complete LiDAR coverage is not available and wall to wall estimation is carried out combining LiDAR, Landsat and ancillary data (chapter 5). In chapter 4, a general aboveground biomass plot-level LiDAR model for tropical forest (Asner & Mascaro, 2014) is validated and a specific model for the study area is fitted. Both LiDAR plot-level models are based on the top-of-canopy height (TCH) variable that is derived from the LiDAR digital canopy model. Results show that the pantropical plot-level LiDAR methodology is a reliable alternative to the development of specific models for tropical forests and thus, aboveground biomass in a new study area could be estimated by only measuring basal area (BA). Applying this methodology, the definition of precise BA field measurement procedures (e.g. location, size and shape of the field plots) is decisive to achieve reliable results in future studies. The relation between BA and TCH (Stocking Coefficient) obtained in our study area in Costa Rica varied locally. Therefore, more field work is needed for assessing Stocking Coefficient variations between different life zones and the influence of the stratification of the study areas in tropical forests on the reduction of uncertainty. In chapter 5, the combination of systematic LiDAR information sampling and full coverage Landsat imagery (and ancillary data) prove to be an effective alternative for forest inventories in tropical areas. This methodology allows estimating wall to wall vegetation height, biomass and carbon density in large areas where full LiDAR coverage and traditional field work are technically and/or economically unfeasible. Carbon density prediction using Landsat imaginery shows a slight decrease in the determination coefficient and an increase in RMSE when harshly decreasing LiDAR coverage area. Results indicate that feasible estimates of vegetation height, biomass and carbon density can be accomplished using low LiDAR coverage areas (between 5% and 20% of the total area) in tropical locations.

Geografía del carbono en alta resolución en bosque tropical amazónico del Ecuador mediante sensores aerotransportados

La estimación de la biomasa de la vegetación terrestre en bosque tropical no sólo es un área de investigación en rápida expansión, sino también es un tema de gran interés para reducir las emisiones de carbono asociadas a la deforestación y la degradación forestal (REDD+). Las estimaciones de densidad de carbono sobre el suelo (ACD) en base a inventarios de campo y datos provenientes de sensores aerotransportados, en especial con sensores LiDAR, han conducido a un progreso sustancial en el cartografiado a gran escala de las reservas de carbono forestal. Sin embargo, estos mapas de carbono tienen incertidumbres considerables, asociadas generalmente al proceso de calibración del modelo de regresión utilizado para producir los mapas. En esta tesis se establece una metodología para la calibración y validación de un modelo general de estimación de ACD usando LiDAR en un sector del Parque Nacional Yasuní en Ecuador. En el proceso de calibración del modelo se considera el tamaño y la ubicación de las parcelas, la influencia de la topografía y la distribución espacial de la biomasa. Para el análisis de los datos se utilizan técnicas geoestadísticas en combinación con variables geomorfométricas derivadas de datos LiDAR, y se propone un esquema de muestreo estratificado por posiciones topográficas (valle, ladera y cima). La validación del modelo general para toda la zona de estudio presentó valores de RMSE = 5.81 Mg C ha-1, R2 = 0.94 y sesgo = 0.59, mientras que, al considerar las posiciones topográficas, el modelo presentó valores de RMSE = 1.67 Mg C ha-1, R2 = 0.98 y sesgo = 0.23 para el valle; RMSE = 3.13 Mg C ha-1, R2 = 0.98 y sesgo = - 0.34 para la ladera; y RMSE = 2.33 Mg C ha-1, R2 = 0.97 y sesgo = 0.74 para la cima. Los resultados obtenidos demuestran que la metodología de muestreo estratificado por posiciones topográficas propuesto, permite calibrar de manera efectiva el modelo general con las estimaciones de ACD en campo, logrando reducir el RMSE y el sesgo. Los resultados muestran el potencial de los datos LiDAR para caracterizar la estructura vertical de la vegetación en un bosque altamente diverso, permitiendo realizar estimaciones precisas de ACD, y conocer patrones espaciales continuos de la distribución de la biomasa aérea y del contenido de carbono en la zona de estudio. ABSTRACT Estimating biomass of terrestrial vegetation in tropical forest is not only a rapidly expanding research area, but also a subject of tremendous interest for reducing carbon emissions associated with deforestation and forest degradation (REDD+). The aboveground carbon density estimates (ACD) based on field inventories and airborne sensors, especially LiDAR sensors have led to a substantial progress in large-scale mapping of forest carbon stocks. However, these carbon maps have considerable uncertainties generally associated with the calibration of the regression model used to produce these maps. This thesis establishes a methodology for calibrating and validating a general ACD estimation model using LiDAR in Ecuador´s Yasuní National Park. The size and location of the plots are considered in the model calibration phase as well as the influence of topography and spatial distribution of biomass. Geostatistical analysis techniques are used in combination with geomorphometrics variables derived from LiDAR data, and then a stratified sampling scheme considering topographic positions (valley, slope and ridge) is proposed. The validation of the general model for the study area showed values of RMSE = 5.81 Mg C ha-1, R2 = 0.94 and bias = 0.59, while considering the topographical positions, the model showed values of RMSE = 1.67 Mg C ha-1, R2 = 0.98 and bias = 0.23 for the valley; RMSE = 3.13 Mg C ha-1, R2 = 0.98 and bias = - 0.34 for the slope; and RMSE = 2.33 Mg C ha-1, R2 = 0.97 and bias = 0.74 for the ridge. The results show that the stratified sampling methodology taking into account topographic positions, effectively calibrates the general model with field estimates of ACD, reducing RMSE and bias. The results show the potential of LiDAR data to characterize the vertical structure of vegetation in a highly diverse forest, allowing accurate estimates of ACD, and knowing continuous spatial patterns of biomass distribution and carbon stocks in the study area.

Crecimiento forestal en el bosque tropical de montaña: efectos de la diversidad florística y de la manipulación de nutrientes

El bosque tropical de montaña, es considerado zona de mega diversidad y de alto grado de endemismo, por las diferentes zonas ecológicas que presenta. Durante las últimas décadas estos bosques han recibido mayor atención por parte de investigadores, sin embargo, aún existe poca investigación en determinar cuáles son las respuestas de los bosques a los cambios ambientales a los que son sometidos. Estos bosques están sufriendo serias amenazas como pérdida de cobertura vegetal y cambios en los ciclos de nutrientes. El trabajo se dividió en cuatro objetivos específicos. i) Caracterización y análisis de patrones altitudinales de la riqueza de especies en el bosque tropical de montaña en el sur del Ecuador; con el fin de conocer cómo varía la diversidad de especies riqueza lo largo de un gradiente altitudinal. ii) Conocer los patrones espaciales del crecimiento en tres remanentes boscosos de un bosque tropical de montaña para determinar cómo la vecindad y la semejanza funcional de ésta influyen en el crecimiento forestal. iii) Conocer los efectos de la fertilización en el crecimiento diamétrico de especies arbóreas, en el bosque tropical de montaña; se analizó cómo reaccionan los árboles a la adición de nutrientes N y P en tres tipos de bosque. iv) Saber la respuesta de la comunidad de árboles a la adición de nutrientes en el bosque montano andino; este objetivo se basó con el supuesto de la deficiencia de tres tipos de nutrientes N, P y Ca, en esta formación boscosa y cómo reaccionan los árboles a la adición de nutrientes. El presente trabajo se llevó a cabo, en el bosque tropical de montaña que se encuentra localizada en la parte adyacente del Parque Nacional Podocarpus (PNP) en la cordillera del Consuelo, forma parte de la cadena oriental de los Andes del sur del Ecuador El trabajo de desarrollo entre los años 2008 y 2014. Para abordar el primer objetivo se establecieron 54 parcelas ubicadas aleatoriamente a lo largo de un gradiente altitudinal (3 niveles de altitud) y se e midieron e identificaron todos los individuos mayores a 5cm de DAP. Se construyó una filogenia con Phylocom y se calcularon diferentes componentes de diversidad para cada parcela ( riqueza taxonómica, diversidad filogenética y edad media de las especies). Ajustando modelos lineares se contrastó el efecto de la altitud sobre dichos componentes y se vio que la riqueza taxonómica y la edad media de las especies aumentaron con la altitud, en sentido contrario a las predicciones de la "hipótesis del conservadurismo tropical" (Tropical Conservatism Hypothesis). Para abordar el segundo objetivo se realizó una remedición de todos los árboles cartografiados en tres parcelas permanentes de alrededor de 5000 m2 cada una, representativas de tres estados diferentes de la sucesión del bosque montano. A partir de las coordenadas y de los datos de registrados, y empleando diferentes funciones de correlación de marca se analizó la distribución espacial del tamaño y del crecimiento relativo y del tamaño. Se constató que mientras que el tamaño de los árboles presentó una correlación espacial negativa, el crecimiento presentó correlación espacial positiva, en ambos casos a distancias cortas. El rango y la magnitud de ambas correlaciones aumentaron al avanzar la sucesión. La distribución espacial del crecimiento mostró una correlación negativa con la distribución espacial de tamaños. Por otro lado, la distribución espacial del crecimiento mostró una correlación negativa para árboles semejantes funcionalmente y positiva cuando se calculó entre árboles con diferente estrategia funcional. En conjunto, los resultados obtenidos señalan un aumento de la importancia de procesos competitivos y una mayor estructuración espacial del crecimiento y de la distribución de tamaños al avanzar la sucesión. Para el tercer y cuarto objetivo se instalaron 52 parcelas distribuidas en bloques donde se fertilizaron dos veces al año durante 6,4 años, se identificaron todos los individuos mayores a 10 cm de DAP, y se midió el crecimiento diamétrico durante estos años Con la adición de nutrientes realizada a los diferentes tipos de bosque en la gradiente altitudinal, encontramos que el efecto sobre el crecimiento diamétrico en la comunidad varia con el rango altitudinal, y el tipo de nutriente, analizando a nivel de las especies, en la mayoría de los casos las especies comunes no tuvieron cambios significativos a la adición de nutrientes. Los resultados de este estudio aportan nuevas evidencias para el entendimiento de la diversidad, estructura y dinámica de los bosques tropicales de montaña. ABSTRACT The montane tropical forest is considered a megadiverse habitat that harbor an enormous degree of endemism. This is mainly due to the high degree of environmental heterogeneity found and the presence of different well defined ecological areas. These forests have received more attention during the last decades, however, the information regarding the responses of these forests to environmental change, is still scarce. These forests are seriously endangered and are suffering serious threats, such as loss of vegetative cover, changes in the nutrient cycles. The work was divided in four specific objectives: i) Characterization and analysis of the species richness altitudinal patterns in the montane tropical forest of south Ecuador. Specifically, how species diversity changes along altitudinal gradients. ii) Exploring the spatial patterns of tree growth in three remnants of a montane tropical forest, and analyze how tree neighborhood and functional similarity among trees influence tree growth. Tropical Conservatism Hypothesis iii) Understanding the effects of fertilization in arboreal species growth (increase in diameter) of the montane tropical forest. Specifically we studied the effects of P and N addition on three different forests across an altitudinal gradient. iv) Know the response of the community of trees to the addition of nutrients in the Andean montane forest; this objective was based on the supposition of deficiency of three types of nutrients: P, N and Ca in this forest all formation and how the trees react to the addition of these nutrients. The present work was carried out in the montane tropical forest located in Bombuscaro, San Francisco and Cajanuma close to Podocarpus National Park (PNP) on Consuelo mountain range (Andean oriental range) at South of Ecuador. Field work was carried out during 2008 and 2014. To address the first objective, we randomly placed 54 plots along an altitudinal gradient. In these plots, every individual larger than 5 cm of DBH was measured and identified. A phylogeny was build with Phylocom and different diversity components (taxonomic richness, phylogenetic diversity and average species age) were computed for each plot. Linear models were used to test the effects of altitude on the diversity components. Our results showed that, contrary to the Tropical Conservatism Hypothesis, both taxonomic richness and average species age increased with altitude. To address our second objective, all mapped trees in three successional permanent plots (around ~5000 m2 each) were re-measured. Using different mark correlation functions, we analyzed the spatial distribution of tree-size and tree relative growth rate. Whereas tree size showed negative spatial correlation at fine spatial scales, relative growth rate showed positive correlation at the same scales. The range and magnitude of those correlations increased along successional stage. The spatial distribution of the relative growth rate was negatively correlated with the spatial distribution of tree sizes. Additionally, we found that the spatial correlation of the relative growth rate was negative for functionally similar trees and positive when computed for functionally different trees. In synthesis, our results point to an increase of competitive processes and strong spatial structure of relative growth rate and tree size along succession. For the third and fourth objectives, 52 plots were placed in a block design and were fertilized twice a year for 6,4 years. In these plots all the individuals with DBH > 10 cm were identified, and the diametrical growth was measured during these years. The nutrient addition at the three different altitude forests, revealed that the effect on the diametrical growth in the community varied with the altitudinal range. When analyzed at species level, the addition of nutrients was no significant in most cases. These results represent new evidences that will improved our understanding of diversity patterns and structure, and the dynamics of tropical montane forests.

Lipid biomarkers in the Enol Lake (Asturias, Northern Spain): coupled natural and human induced environmental history

The lipid content of three cores from Lake Enol (Picos de Europa National Park, Asturias, Northern Spain) was studied. The n-alkane profiles indicated a major input from terrigenous plants [predominance of high molecular weight (HMW) alkanes] since ca. 1695 AD to the water body, although the uppermost cm revealed a predominance of organic matter (OM) derived from algae, as the most abundant alkane was C17. Three units revealing different environmental conditions were defined. Unit A (ca. 1695–1860 AD) in the lowermost parts of ENO13-10 (< 12 cm) and ENO13-15 (< 28 cm) was identified and was characterized by higher OM input and evidence of minimal degradation (high CPI values, predominance of HMW n-alkanoic acids and good correspondence between the predominant n-alkane and n-alkanoic acid chains). These findings could be linked to the Little Ice Age, when cold and humid conditions may have favored an increase in total organic carbon (TOC) and n-alkane and n-alkanoic acid content (greater terrigenous OM in-wash), and may have also reduced bacterial activity. In Unit B (ca. 1860–1980 AD) the lack of correspondence between the n-alkane and n-alkanoic acid profiles of ENO13-10 (12–4 cm) and ENO13-15 (28–8 cm) suggested a certain preferential microbial synthesis of long chain saturated fatty acids from primary OM and/or bacterial activity, coinciding with a decrease in OM input, which could be linked to the global warming that started in the second half of the 19th century. In ENO13-7 the low OM input (low TOC) was accompanied by some bacterial degradation (predominance ofHMWn-alkanoic acids but with a bimodal distribution) in the lowermost 16–5 cm. Evidence of considerable phytoplankton productivity and microbial activity was especially significant in Unit C (ca. 1980–2013 AD) identified in the uppermost part of all three cores (5 cm in ENO13-7, 4 cm in ENO13-10 and 8 cm in ENO13-15), coinciding with higher concentrations of n-alkanes and n-alkanoic acids, which were considered to be linked to warmer and drier conditions, as well as to greater anthropogenic influence in modern times. Plant sterols, such as b-sitosterol, campesterol and stigmasterol, were significantly present in the cores. In addition, fecal stanols, such as 24-ethylcoprostanol from herbivores, were present, thereby indicating a continuous and significant pollution input derived from these animals since the 17th century, being more important in the last 20 years.

A genetic record of population isolation in pocket gophers during Holocene climatic change

A long-standing question in Quaternary paleontology is whether climate-induced, population-level phenotypic change is a result of large-scale migration or evolution in isolation. To directly measure genetic variation through time, ancient DNA and morphologic variation was measured over 2,400 years in a Holocene sequence of pocket gophers (Thomomys talpoides) from Lamar Cave, Yellowstone National Park, Wyoming. Ancient specimens and modern samples collected near Lamar Cave share mitochondrial cytochrome b sequences that are absent from adjacent localities, suggesting that the population was isolated for the entire period. In contrast, diastemal length, a morphologic character correlated with body size and nutritional level, changed predictably in response to climatic change. Our results demonstrate that small mammal populations can experience the long-term isolation assumed by many theoretical models of microevolutionary change.

The compensatory potential of increased immigration following intensive American mink population control is diluted by male-biased dispersal

Acknowledgments This research was funded by NERC Grants NE/J01396X/1; NE/E006434/1 to XL. Y. M. was funded by a Marie Curie FP7-PEOPLE-2011-IEF 300288. We profusely thank the Scottish Mink Initiative, staff, funders and multiple mink volunteers for the continued effort, samples and enthusiasm. The Scottish Water Vole Conservation Project was funded by The Tubney Charitable trust, Scottish Natural Heritage, the Cairngorms National Park Authority and the People’s Trust for Endangered Species. We also thank the Associate Editor and two anonymous referees for their thoughtful reviews.

Perspectives on archaeal diversity, thermophily and monophyly from environmental rRNA sequences.

Phylogenetic analysis of ribosomal RNA sequences obtained from uncultivated organisms of a hot spring in Yellowstone National Park reveals several novel groups of Archaea, many of which diverged from the crenarchaeal line of descent prior to previously characterized members of that kingdom. Universal phylogenetic trees constructed with the addition of these sequences indicate monophyly of Archaea, with modest bootstrap support. The data also show a specific relationship between low-temperature marine Archaea and some hot spring Archaea. Two of the environmental sequences are enigmatic: depending upon the data set and analytical method used, these sequences branch deeply within the Crenarchaeota, below the bifurcation between Crenarchaeota and Euryarchaeota, or even as the sister group to Eukaryotes. If additional data confirm either of the latter two placements, then the organisms represented by these ribosomal RNA sequences would merit recognition as a new kingdom, provisionally named "Korarchaeota."