Stand scale variability of topsoil organic matter composition in a high-elevation Norway spruce forest ecosystem

Our knowledge about the effect of single-tree influence areas on the physicochemical properties of the underlying mineral soil in forest ecosystems is still limited. This restricts our ability to adequately estimate future changes in soil functioning due to forest management practices. We studied the stand scale spatial variation of different soil organic matter species investigated by 13C NMR spectroscopy, lignin phenol and neutral sugar analysis under an unmanaged mountainous high-elevation Norway spruce (Picea abies L.) forest in central Europe. Multivariate geostatistical approaches were applied to relate the spatial patterns of the different soil organic matter species to topographic parameters, bulk density, oxalate- and dithionite-extractable iron, pH, and the impact of tree distribution. Soil samples were taken from the mineral top soil. Generally, the stand scale distribution patterns of different soil organic matter compounds could be divided into two groups: Those compounds, which were significantly spatially correlated with topography/altitude and those with small scale spatial pattern (range ≤ 10 m) that was closely related to tree distribution. The concentration of plant-derived soil organic matter components, such as lignin, at a given sampling point was significantly spatially related to the distance of the nearest tree (p ≤ 0.05). In contrast, the spatial distribution of mainly microbial-derived compounds (e.g. galactose and mannose) could be attributed to the dominating impact of small-scale topography and the contribution of poorly crystalline iron oxides that were significantly larger in the central depression of the study site compared to crest and slope positions. Our results demonstrate that topographic parameters dominate the distribution of overall topsoil organic carbon (OC) stocks at temperate high-elevation forest ecosystems, particularly in sloped terrain. However, trees superimpose topography-controlled OC biogeochemistry beneath their crown by releasing litter and changing soil conditions in comparison to open areas. This may lead to distinct zones with different mechanisms of soil organic matter degradation and also stabilization in forest stands.

Modelización espacial de huevos de Thrips (Thysanoptera: Frankliniella occidentalis) en tomate de cáscara por medio de técnicas geoestadísticas

El tomate de cáscara (Physalis ixocarpa Brot.) es un cultivo alimenticio de gran importancia económica en México. Sin embargo, es afectado por diversas plagas y enfermedades tales como los Thrips (Thysanoptera: Frankliniella occidentalis) y el virus de la marchitez manchada del tomate (TSWV) que llegan a causar hasta un 80% de pérdidas. El objetivo del presente trabajo fue modelizar la distribución espacial de huevos de Thrips mediante técnicas geoestadísticas y obtener, en consecuencia, mapas de incidencia por medio del Kriging. Se georreferenciaron 121 puntos de muestreo en cada una de las parcelas comerciales de los municipios de Luvianos, Jocotitlán e Ixtlahuaca, a través del método de transectos en tres etapas fenológicas del cultivo. Se contabilizó el número de huevos de Thrips en cada punto de muestreo. Los resultados mostraron que las poblaciones de huevos deThrips presentan una distribución agregada, identificándose varios centros de conglomeración a través de los mapas obtenidos. Los semivariogramas obtenidos de la distribución espacial se ajustaron principalmente a los modelos gaussianos y esféricos. La distribución de huevos de Thrips se presentó en centros de agregación dentro de las parcelas estudiadas, lo cual permitirá establecer estrategias y medidas de control o mitigación en términos de sitios específicos de infestación de huevos de Thrips.

Probability inundation map for a scenario combining storm surge and sea-level rise

Coastal managers require reliable spatial data on the extent and timing of potential coastal inundation, particularly in a changing climate. Most sea level rise (SLR) vulnerability assessments are undertaken using the easily implemented bathtub approach, where areas adjacent to the sea and below a given elevation are mapped using a deterministic line dividing potentially inundated from dry areas. This method only requires elevation data usually in the form of a digital elevation model (DEM). However, inherent errors in the DEM and spatial analysis of the bathtub model propagate into the inundation mapping. The aim of this study was to assess the impacts of spatially variable and spatially correlated elevation errors in high-spatial resolution DEMs for mapping coastal inundation. Elevation errors were best modelled using regression-kriging. This geostatistical model takes the spatial correlation in elevation errors into account, which has a significant impact on analyses that include spatial interactions, such as inundation modelling. The spatial variability of elevation errors was partially explained by land cover and terrain variables. Elevation errors were simulated using sequential Gaussian simulation, a Monte Carlo probabilistic approach. 1,000 error simulations were added to the original DEM and reclassified using a hydrologically correct bathtub method. The probability of inundation to a scenario combining a 1 in 100 year storm event over a 1 m SLR was calculated by counting the proportion of times from the 1,000 simulations that a location was inundated. This probabilistic approach can be used in a risk-aversive decision making process by planning for scenarios with different probabilities of occurrence. For example, results showed that when considering a 1% probability exceedance, the inundated area was approximately 11% larger than mapped using the deterministic bathtub approach. The probabilistic approach provides visually intuitive maps that convey uncertainties inherent to spatial data and analysis.

Spatially explicit methane inventory for Switzerland with 500 m x 500 m resolution for 2011 (NetCDF 25 MB)

We present the first high-resolution (500 m × 500 m) gridded methane (CH4) emission inventory for Switzerland, which integrates the national emission totals reported to the United Nations Framework Convention on Climate Change (UNFCCC) and recent CH4 flux studies conducted by research groups across Switzerland. In addition to anthropogenic emissions, we also include natural and semi-natural CH4 fluxes, i.e., emissions from lakes and reservoirs, wetlands, wild animals as well as uptake by forest soils. National CH4 emissions were disaggregated using detailed geostatistical information on source locations and their spatial extent and process- or area-specific emission factors. In Switzerland, the highest CH4 emissions in 2011 originated from the agricultural sector (150 Gg CH4/yr), mainly produced by ruminants and manure management, followed by emissions from waste management (15 Gg CH4/yr) mainly from landfills and the energy sector (12 Gg CH4/yr), which was dominated by emissions from natural gas distribution. Compared to the anthropogenic sources, emissions from natural and semi-natural sources were relatively small (6 Gg CH4/yr), making up only 3 % of the total emissions in Switzerland. CH4 fluxes from agricultural soils were estimated to be not significantly different from zero (between -1.5 and 0 Gg CH4/yr), while forest soils are a CH4 sink (approx. -2.8 Gg CH4/yr), partially offsetting other natural emissions. Estimates of uncertainties are provided for the different sources, including an estimate of spatial disaggregation errors deduced from a comparison with a global (EDGAR v4.2) and a European CH4 inventory (TNO/MACC). This new spatially-explicit emission inventory for Switzerland will provide valuable input for regional scale atmospheric modeling and inverse source estimation.

Small-scale variation of vegetation in a mixed forest understorey is partly controlled by the effect of overstory composition on litter accumulation

We investigated how richness and composition of vascular plant species in the understory of a mixed hardwood forest stand varied with respect to the abundance and composition of the overstory. The stand is in central Spain and represents the southernmost range of distribution of several tree and herbaceous species in Europe. Understory species were identified in 46 quadrats (0.25 m2) where variables litter depth and light availability were measured. In addition, we estimated tree density, basal area, and percent basal area by tree species within 6-m-radius areas around each plot. Species richness and composition were studied using path analysis and scale-dependent geostatistical methods, respectively. We found that the relative abundance of certain trees species in the overstory was more important than total overstory abundance in explaining understory species richness. Richness decreased as soil litter depth increased, and soil litter increased as the relative proportion of Fagus sylvatica in the overstory increased, which accounted for a negative, indirect effect of Fagus sylvatica on richness. Regarding understory species composition, we found that some species distributed preferentially below certain tree species. For example, Melica uniflora was most frequent below Fagus sylvatica and Quercus petraea while the increasing proportion of Q. pyrenaica in the overstory favored the presence of Cruciata glabra, Arenaria montana, Prunus avium, Conopodium bourgaei, Holcus mollis, Stellaria media and Galium aparine in the understory. Overall, these results emphasize the importance of individual tree species in controlling the assemblage and richness of understory species in mixed stands. We conclude that soil litter accumulation is one way through which overstory composition shapes the understory community.

Sampling strategies using the "accumulation chamber" for monitoring geological

Fundación Ciudad de la Energía (CIUDEN) is carrying out a project of geological storage of CO2, where CO2 injection tests are planned in saline aquifers at a depth of 1500 m for scientific objectives and project demonstration. Before any CO2 is stored, it is necessary to determine the baseline flux of CO2 in order to detect potential leakage during injection and post-injection monitoring. In November 2009 diffuse flux measurements of CO2 using an accumulationchamber were made in the area selected by CIUDEN for geological storage, located in Hontomin province of Burgos (Spain). This paper presents the tests carried out in order to establish the optimum sampling methodology and the geostatistical analyses performed to determine the range, with which future field campaigns will be planned.

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

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

Comparative assessment of global irradiation from a satellite estimate model (CM SAF) and on-ground measurements (SIAR): a Spanish case study

An analysis and comparison of daily and yearly solar irradiation from the satellite CM SAF database and a set of 301 stations from the Spanish SIAR network is performed using data of 2010 and 2011. This analysis is completed with the comparison of the estimations of effective irradiation incident on three different tilted planes (fixed, two axis tracking, north-south hori- zontal axis) using irradiation from these two data sources. Finally, a new map of yearly values of irradiation both on the horizontal plane and on inclined planes is produced mixing both sources with geostatistical techniques (kriging with external drift, KED) The Mean Absolute Difference (MAD) between CM SAF and SIAR is approximately 4% for the irradiation on the horizontal plane and is comprised between 5% and 6% for the irradiation incident on the inclined planes. The MAD between KED and SIAR, and KED and CM SAF is approximately 3% for the irradiation on the horizontal plane and is comprised between 3% and 4% for the irradiation incident on the inclined planes. The methods have been implemented using free software, available as supplementary ma- terial, and the data sources are freely available without restrictions.

Comparative assessment of global irradiation from a satellite estimate model (CM SAF) and on-ground measurements (SIAR): a Spanish case study

An analysis and comparison of daily and yearly solar irradiation from the satellite CM SAF database and a set of 301 stations from the Spanish SIAR network is performed using data of 2010 and 2011. This analysis is completed with the comparison of the estimations of effective irradiation incident on three different tilted planes (fixed, two axis tracking, north-south hori- zontal axis) using irradiation from these two data sources. Finally, a new map of yearly values of irradiation both on the horizontal plane and on inclined planes is produced mixing both sources with geostatistical techniques (kriging with external drift, KED) The Mean Absolute Difference (MAD) between CM SAF and SIAR is approximately 4% for the irradiation on the horizontal plane and is comprised between 5% and 6% for the irradiation incident on the inclined planes. The MAD between KED and SIAR, and KED and CM SAF is approximately 3% for the irradiation on the horizontal plane and is comprised between 3% and 4% for the irradiation incident on the inclined planes. The methods have been implemented using free software, available as supplementary ma- terial, and the data sources are freely available without restrictions.

Integración de la calidad del aire en la planificación física

En esta tesis se presenta una nueva aproximación para la realización de mapas de calidad del aire, con objeto de que esta variable del medio físico pueda ser tenida en cuenta en los procesos de planificación física o territorial. La calidad del aire no se considera normalmente en estos procesos debido a su composición y a la complejidad de su comportamiento, así como a la dificultad de contar con información fiable y contrastada. Además, la variabilidad espacial y temporal de las medidas de calidad del aire hace que sea difícil su consideración territorial y exige la georeferenciación de la información. Ello implica la predicción de medidas para lugares del territorio donde no existen datos. Esta tesis desarrolla un modelo geoestadístico para la predicción de valores de calidad del aire en un territorio. El modelo propuesto se basa en la interpolación de las medidas de concentración de contaminantes registradas en las estaciones de monitorización, mediante kriging ordinario, previa homogeneización de estos datos para eliminar su carácter local. Con el proceso de eliminación del carácter local, desaparecen las tendencias de las series muestrales de datos debidas a las variaciones temporales y espaciales de la calidad del aire. La transformación de los valores de calidad del aire en cantidades independientes del lugar de muestreo, se realiza a través de parámetros de uso del suelo y de otras variables características de la escala local. Como resultado, se obtienen unos datos de entrada espacialmente homogéneos, que es un requisito fundamental para la utilización de cualquier algoritmo de interpolación, en concreto, del kriging ordinario. Después de la interpolación, se aplica una retransformación de los datos para devolver el carácter local al mapa final. Para el desarrollo del modelo, se ha elegido como área de estudio la Comunidad de Madrid, por la disponibilidad de datos reales. Estos datos, valores de calidad del aire y variables territoriales, se utilizan en dos momentos. Un momento inicial, donde se optimiza la selección de los parámetros más adecuados para la eliminación del carácter local de las medidas y se desarrolla cada una de las etapas del modelo. Y un segundo momento, en el que se aplica en su totalidad el modelo desarrollado y se contrasta su eficacia predictiva. El modelo se aplica para la estimación de los valores medios y máximos de NO2 del territorio de estudio. Con la implementación del modelo propuesto se acomete la territorialización de los datos de calidad del aire con la reducción de tres factores clave para su efectiva integración en la planificación territorial o en el proceso de toma de decisiones asociado: incertidumbre, tiempo empleado para generar la predicción y recursos (datos y costes) asociados. El modelo permite obtener una predicción de valores del contaminante objeto de análisis en unas horas, frente a los periodos de modelización o análisis requeridos por otras metodologías. Los recursos necesarios son mínimos, únicamente contar con los datos de las estaciones de monitorización del territorio que, normalmente, están disponibles en las páginas web viii institucionales de los organismos gestores de las redes de medida de la calidad del aire. Por lo que respecta a las incertidumbres de la predicción, puede decirse que los resultados del modelo propuesto en esta tesis son estadísticamente muy correctos y que los errores medios son, en general, similares o menores que los encontrados con la aplicación de las metodologías existentes. ABSTRACT This thesis presents a new approach for mapping air quality, so that this variable of physical environment can be taken into account in physical or territorial planning. Ambient air quality is not normally considered in territorial planning mainly due to the complexity of its composition and behavior and the difficulty of counting with reliable and contrasted information. In addition, the wide spatial and temporal variability of the measurements of air quality makes his territorial consideration difficult and requires georeferenced information. This involves predicting measurements in the places of the territory where there are no data. This thesis develops a geostatistical model for predicting air quality values in a territory. The proposed model is based on the interpolation of measurements of pollutants from the monitoring stations, using ordinary kriging, after a detrending or removal of the local character of sampling values process. With the detrending process, the local character of the time series of sampling data, due to temporal and spatial variations of air quality, is removed. The transformation of the air quality values into site-independent quantities is performed using land use parameters and other characteristic parameters of local scale. This detrending of the monitoring data process results in a spatial homogeneous input set which is a prerequisite for a correct use of any interpolation algorithm, particularly, ordinary kriging. After the interpolation step, a retrending or retransformation is applied in order to incorporate the local character in the final map at places where no monitoring data is available. For the development of this model, the Community of Madrid is chosen as study area, because of the availability of actual data. These data, air quality values and local parameters, are used in two moments. A starting point, to optimize the selection of the most suitable indicators for the detrending process and to develop each one of the model stages. And a second moment, to fully implement the developed model and to evaluate its predictive power. The model is applied to estimate the average and maximum values of NO2 in the study territory. With the implementation of the proposed model, the territorialization of air quality data is undertaken with the reduction in three key factors for the effective integration of this parameter in territorial planning or in the associated decision making process: uncertainty, time taken to generate the prediction and associated resources (data and costs). This model allows the prediction of pollutant values in hours, compared to the implementation time periods required for other modeling or analysis methodologies. The required resources are also minimal, only having data from monitoring stations in the territory, that are normally available on institutional websites of the authorities responsible for air quality networks control and management. With regard to the prediction uncertainties, it can be concluded that the results of the proposed model are statistically very accurate and the mean errors are generally similar to or lower than those found with the application of existing methodologies.

Sampling strategies using the "accumulation chamber" for monitoring geological storage of CO2

Fundación Ciudad de la Energía (CIUDEN) is carrying out a project of geological storage of CO2, where CO2 injection tests are planned in saline aquifers at a depth of 1500 m for scientific objectives and project demonstration. Before any CO2 is stored, it is necessary to determine the baseline flux of CO2 in order to detect potential leakage during injection and post-injection monitoring. In November 2009 diffuse flux measurements of CO2 using an accumulation chamber were made in the area selected by CIUDEN for geological storage, located in Hontomin province of Burgos (Spain). This paper presents the tests carried out in order to establish the optimum sampling methodology and the geostatistical analyses performed to determine the range, with which future field campaigns will be planned.

Análisis de la distribución espacial de los parámetros acústicos en recintos

El principal objetivo de este estudio es la evaluación de la distribución espacial de los parámetros acústicos en un recinto a través de la aplicación de técnicas geoestadísticas como el método Kriging. Mediante el uso de la herramienta de análisis espacial ArcMap, perteneciente a la plataforma ArcGIS, se ha analizado el comportamiento acústico del Salón de Actos común a la Escuela Técnica Superior de Ingeniería y Sistemas de Telecomunicación (ETSIST) y la Escuela Técnica Superior de Ingeniería de Sistemas Informáticos (ETSISI), ambas situadas en el Campus Sur de la Universidad Politécnica de Madrid. Se han realizado mediciones in-situ del recinto no ocupado utilizando la herramienta de medición DIRAC Room Acoustics y el método de la respuesta impulsiva integrada, extrayéndose los parámetros acústicos de tiempo de reverberación (RT), tiempo de reverberación inicial (EDT), fuerza sonora relativa (Grel), claridad (C80), tiempo central (Ts), definición (D50) e Índice de Transmisión Rápida de la Palabra (RASTI). Se ha analizado la adecuación de los valores observados de Grel, C80 y Ts al modelo teórico de Barron y estimado, mediante el método Kriging Ordinario, el conjunto de parámetros medidos en el recinto, obteniéndose los semivariograma y mapas de estimación correspondientes. Además, se ha evaluado la calidad de la estimación en base a un número de puntos de medición reducido. A la vista de los resultados obtenidos, en general, el método Kriging puede considerarse un buen interpolador de los parámetros acústicos en un recinto, observándose que los parámetros que evalúan relaciones energéticas, especialmente la fuerza sonora relativa (Grel) proporcionan mejores estimaciones en comparación con aquellos relacionados con la reverberación y la inteligibilidad del habla. El coeficiente de determinación (R2) constituye una medida útil para evaluar la precisión de la estimación. Además, la entropía de los datos observados puede ser un buen indicador a priori de la precisión de la estimación. Asimismo, se ha demostrado que, basándose en un reducido número de puntos de medición, es posible obtener una estimación precisa de los parámetros acústicos de fuerza sonora relativa (Grel) y tiempo central (Ts). ABSTRACT. This project aims to evaluate the feasibility of using geostatistical techniques such as Kriging on the analysis of the spatial distribution of the acoustic parameters in rooms. The acoustic behaviour of the Assembly Hall of the ETSIST and ETSISI (Universidad Politécnica de Madrid) is investigated using ArcMap, which is the main component of ArcGIS suite of geospatial processing programs. For this purpose, in-situ acoustic measurements are carried out in the unoccupied room using DIRAC Room Acoustics software. The following acoustic parameters are measured by means of the integrated impulse response method for further examination: Reverberation Time (RT), Early Decay Time (EDT), Relative Strength (Grel), Clarity (C80), Centre Time (Ts), Definition (D50) and Rapid Speech Transmission Index (RASTI). Goodness-of-fit of measured Grel, C80 and Ts values to Barron’s theory is determined and Ordinary Kriging is applied to all the measured parameters in order to calculate the semivariogram and prediction surfaces. The prediction performance is also analysed when significantly fewer receiver positions are used for the prediction. The experimental results obtained lead to conclude that Kriging can be successfully applied to room acoustics. Energy\based acoustic parameters can be estimated with higher accuracy compared to those related to reverberation and speech intelligibility. Coefficient of determination (R2) is a reliable statistic for assessing the prediction accuracy, for which measured data entropy can also be a good a priori indicator. Furthermore, based on fewer receiver positions, it is demonstrated that accurate predictions of Grel and Ts can be achieved.

Predicción espacio-temporal de la irradiancia solar global a corto plazo en España mediante geoestadística y redes neuronales artificiales

El enriquecimiento del conocimiento sobre la Irradiancia Solar (IS) a nivel de superficie terrestre, así como su predicción, cobran gran interés para las Energías Renovables (ER) - Energía Solar (ES)-, y para distintas aplicaciones industriales o ecológicas. En el ámbito de las ER, el uso óptimo de la ES implica contar con datos de la IS en superficie que ayuden tanto, en la selección de emplazamientos para instalaciones de ES, como en su etapa de diseño (dimensionar la producción) y, finalmente, en su explotación. En este último caso, la observación y la predicción es útil para el mercado energético, la planificación y gestión de la energía (generadoras y operadoras del sistema eléctrico), especialmente en los nuevos contextos de las redes inteligentes de transporte. A pesar de la importancia estratégica de contar con datos de la IS, especialmente los observados por sensores de IS en superficie (los que mejor captan esta variable), estos no siempre están disponibles para los lugares de interés ni con la resolución espacial y temporal deseada. Esta limitación se une a la necesidad de disponer de predicciones a corto plazo de la IS que ayuden a la planificación y gestión de la energía. Se ha indagado y caracterizado las Redes de Estaciones Meteorológicas (REM) existentes en España que publican en internet sus observaciones, focalizando en la IS. Se han identificado 24 REM (16 gubernamentales y 8 redes voluntarios) que aglutinan 3492 estaciones, convirtiéndose éstas en las fuentes de datos meteorológicos utilizados en la tesis. Se han investigado cinco técnicas de estimación espacial de la IS en intervalos de 15 minutos para el territorio peninsular (3 técnicas geoestadísticas, una determinística y el método HelioSat2 basado en imágenes satelitales) con distintas configuraciones espaciales. Cuando el área de estudio tiene una adecuada densidad de observaciones, el mejor método identificado para estimar la IS es el Kriging con Regresión usando variables auxiliares -una de ellas la IS estimada a partir de imágenes satelitales-. De este modo es posible estimar espacialmente la IS más allá de los 25 km identificados en la bibliografía. En caso contrario, se corrobora la idoneidad de utilizar estimaciones a partir de sensores remotos cuando la densidad de observaciones no es adecuada. Se ha experimentado con el modelado de Redes Neuronales Artificiales (RNA) para la predicción a corto plazo de la IS utilizando observaciones próximas (componentes espaciales) en sus entradas y, los resultados son prometedores. Así los niveles de errores disminuyen bajo las siguientes condiciones: (1) cuando el horizonte temporal de predicción es inferior o igual a 3 horas, las estaciones vecinas que se incluyen en el modelo deben encentrarse a una distancia máxima aproximada de 55 km. Esto permite concluir que las RNA son capaces de aprender cómo afectan las condiciones meteorológicas vecinas a la predicción de la IS. ABSTRACT ABSTRACT The enrichment of knowledge about the Solar Irradiance (SI) at Earth's surface and its prediction, have a high interest for Renewable Energy (RE) - Solar Energy (SE) - and for various industrial and environmental applications. In the field of the RE, the optimal use of the SE involves having SI surface to help in the selection of sites for facilities ES, in the design stage (sizing energy production), and finally on their production. In the latter case, the observation and prediction is useful for the market, planning and management of the energy (generators and electrical system operators), especially in new contexts of smart transport networks (smartgrid). Despite the strategic importance of SI data, especially those observed by sensors of SI at surface (the ones that best measure this environmental variable), these are not always available to the sights and the spatial and temporal resolution desired. This limitation is bound to the need for short-term predictions of the SI to help planning and energy management. It has been investigated and characterized existing Networks of Weather Stations (NWS) in Spain that share its observations online, focusing on SI. 24 NWS have been identified (16 government and 8 volunteer networks) that implies 3492 stations, turning it into the sources of meteorological data used in the thesis. We have investigated five technical of spatial estimation of SI in 15 minutes to the mainland (3 geostatistical techniques and HelioSat2 a deterministic method based on satellite images) with different spatial configurations. When the study area has an adequate density of observations we identified the best method to estimate the SI is the regression kriging with auxiliary variables (one of them is the SI estimated from satellite images. Thus it is possible to spatially estimate the SI beyond the 25 km identified in the literature. Otherwise, when the density of observations is inadequate the appropriateness is using the estimates values from remote sensing. It has been experimented with Artificial Neural Networks (ANN) modeling for predicting the short-term future of the SI using observations from neighbor’s weather stations (spatial components) in their inputs, and the results are promising. The error levels decrease under the following conditions: (1) when the prediction horizon is less or equal than 3 hours the best models are the ones that include data from the neighboring stations (at a maximum distance of 55 km). It is concluded that the ANN is able to learn how weather conditions affect neighboring prediction of IS at such Spatio-temporal horizons.

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.

Estudio para un Modelo de Gestión de un Acuífero Costero, mediante Metodologías Participativas y Análisis Geoestadístico en el marco del Desarrollo Local. Manglaralto, Ecuador

El presente estudio se fundamenta en la investigación-acción-participativa (IAP), para buscar alternativas que tiendan al desarrollo local de un territorio. Se centra en la cuenca hidrográfica del rio Manglaralto-Santa Elena-Ecuador, aplicando un sistema metodológico participativo que considera las características peculiares del territorio, que se analizan geoespacialmente reconociendo la influencia de la dinámica de sus cambios y observando los móviles que la propiciaban. A través de mecanismos participativos, se conectan los aspectos técnicos para el conocimiento y el aprovechamiento racional del acuífero costero, con los valores de los habitantes del territorio, para mejorar su abastecimiento de agua y crear nuevas condiciones y oportunidades en el camino del desarrollo local, vislumbrando la sostenibilidad. Cabe indicar que el ente administrativo y propulsor es la Junta de Agua Potable Regional Manglaralto (JAPRM). La hipótesis del estudio considera, que los métodos participativos generan en la comunidad una respuesta basada en su identidad y sus deseos de mejorar, que propiciará una gestión del acuífero costero que conlleve al desarrollo local. Otra hipótesis complementaria estipula que las estrategias del gobierno respecto al turismo propicia un crecimiento en la demanda del agua del acuífero. En Manglaralto-Ecuador, una parroquia de 30.000 habitantes aproximadamente, donde la JAPRM, administra y suministra agua a 23.586 habitantes que cuenta en su organización, llevada por 6 representantes de las comunidades rurales que la conforman, empezaron hace 7 años a buscar una forma de lograr un cambio, de tener agua para el desarrollo de la comunidad. Buscaron ayuda por diferentes medios, políticos, económicas, sociales y encontraron como base fundamental a la cooperación con el Organismo Internacional de Energía Atómica (OIEA) y la Escuela Superior Politécnica del Litoral (ESPOL) para entrelazar aspectos técnicos, ambientales, sociales y culturales. La gestión del acuífero costero, desde la perspectiva del IAP repercute en el desarrollo de Manglaralto. También se realiza un análisis geoespacial-geoestadístico, para vislumbrar aspectos de cambios en el territorio ligados al crecimiento turístico, que afectan a la demanda del recurso agua proveniente del acuífero costero bajo la administración de la JAPRM. La tesis presenta el modelo integral y propio de la comunidad de Manglaralto, que refleja una evolución que alcanzó un apogeo en 2011 y parte del 2012, con 9 pozos de agua que daban servicio los 365 días del año, 24 horas al día ininterrumpidamente. Las condiciones externas (promociones turísticas de la ruta del Spondylus) han repercutido en nuevas problemáticas (crecimiento elevado de la demanda del agua). El acuífero costero se convierte en el emblema y móvil de solución, gracias a la gestión integral y a la interacción IAP que se amolda a la evolución de las condiciones, buscando soluciones para la comunidad y su entorno. El modelo integral del territorio con la participación de sus pobladores, considera el aspecto turístico, como un agente que propicia la mayor demanda del agua. Situación a la que hay que dar respuesta mediante la observación-reflexión en el ciclo del IAP para generar nuevas directrices estratégicas y gestionar el desarrollo local. ABSTRACT The present study is based on the participatory action research (PAR) methodology in order to look for alternatives which tend to the local development of a territory. It focuses on the Manglaralto hydrographic river basin located in Santa Elena-Ecuador through the application of the participatory methodology which considers the peculiar characteristics of the territory. These are geospatially analyzed recognizing the influence of its dynamic of changes and observing the causes that originated them. Through the use of participatory mechanisms, technical aspects are connected for stimulating knowledge and rational use of the coastal aquifer with the values of inhabitants of the territory to improve the water supply and create new conditions of sustainability. It is important to point out that the administrative organism and promoter is the Manglaralto Regional Fresh Water Board (JAPRM). In Manglaralto-Ecuador, a parish of approximately 30,000 inhabitants, the MRFWB manages and supplies water to 23.586 inhabitants. This organization is composed by 6 representatives of rural communities. It started 7 years ago looking for a way to achieve a change, from obtaining water to developing the community. They seeked for help in different fields such as: political, economic and social and they found International Atomic Energy Agency (IAEA) and Escuela Superior Politécnica del Litoral (ESPOL) as a fundamental basis for cooperation to bond technical, environmental, social and cultural aspects. Management of coastal aquifer, from the PAR perspective affects the development of Manglaralto. Also, a geospatial and geostatistical analysis is carried out to distinguish change aspects in territories related to touristy growth which affects the demand of water obtained from the coastal aquifer under the management of the MRFWB. The thesis presents a comprehensive model that belongs to the Manglaralto community and reveals an evolution that reached a peak in 2011 and part of 2012, with 9 water wells that operated the 365 days of the year 24 hours a day without interruption. The external conditions (touristic packages of Spondylus route) have created new problems (higher demand of water). The coastal aquifer is a symbol and solution, thanks to the comprehensive management and PAR interaction which fits the evolution of conditions, looking for solutions for the community and its surroundings. The comprehensive model of territory with the participation of inhabitants considers the touristic aspect as an agent which brings about a higher demand of water. This situation requests a response through the observation-reflection in the PAR cycle to generate new strategic guidelines and promote the local development.