Environmental variables modeling based on remote sensing data using time series analysis: forestry and agricultural applications

En la actualidad, el seguimiento de la dinámica de los procesos medio ambientales está considerado como un punto de gran interés en el campo medioambiental. La cobertura espacio temporal de los datos de teledetección proporciona información continua con una alta frecuencia temporal, permitiendo el análisis de la evolución de los ecosistemas desde diferentes escalas espacio-temporales. Aunque el valor de la teledetección ha sido ampliamente probado, en la actualidad solo existe un número reducido de metodologías que permiten su análisis de una forma cuantitativa. En la presente tesis se propone un esquema de trabajo para explotar las series temporales de datos de teledetección, basado en la combinación del análisis estadístico de series de tiempo y la fenometría. El objetivo principal es demostrar el uso de las series temporales de datos de teledetección para analizar la dinámica de variables medio ambientales de una forma cuantitativa. Los objetivos específicos son: (1) evaluar dichas variables medio ambientales y (2) desarrollar modelos empíricos para predecir su comportamiento futuro. Estos objetivos se materializan en cuatro aplicaciones cuyos objetivos específicos son: (1) evaluar y cartografiar estados fenológicos del cultivo del algodón mediante análisis espectral y fenometría, (2) evaluar y modelizar la estacionalidad de incendios forestales en dos regiones bioclimáticas mediante modelos dinámicos, (3) predecir el riesgo de incendios forestales a nivel pixel utilizando modelos dinámicos y (4) evaluar el funcionamiento de la vegetación en base a la autocorrelación temporal y la fenometría. Los resultados de esta tesis muestran la utilidad del ajuste de funciones para modelizar los índices espectrales AS1 y AS2. Los parámetros fenológicos derivados del ajuste de funciones permiten la identificación de distintos estados fenológicos del cultivo del algodón. El análisis espectral ha demostrado, de una forma cuantitativa, la presencia de un ciclo en el índice AS2 y de dos ciclos en el AS1 así como el comportamiento unimodal y bimodal de la estacionalidad de incendios en las regiones mediterránea y templada respectivamente. Modelos autorregresivos han sido utilizados para caracterizar la dinámica de la estacionalidad de incendios y para predecir de una forma muy precisa el riesgo de incendios forestales a nivel pixel. Ha sido demostrada la utilidad de la autocorrelación temporal para definir y caracterizar el funcionamiento de la vegetación a nivel pixel. Finalmente el concepto “Optical Functional Type” ha sido definido, donde se propone que los pixeles deberían ser considerados como unidades temporales y analizados en función de su dinámica temporal. ix SUMMARY A good understanding of land surface processes is considered as a key subject in environmental sciences. The spatial-temporal coverage of remote sensing data provides continuous observations with a high temporal frequency allowing the assessment of ecosystem evolution at different temporal and spatial scales. Although the value of remote sensing time series has been firmly proved, only few time series methods have been developed for analyzing this data in a quantitative and continuous manner. In the present dissertation a working framework to exploit Remote Sensing time series is proposed based on the combination of Time Series Analysis and phenometric approach. The main goal is to demonstrate the use of remote sensing time series to analyze quantitatively environmental variable dynamics. The specific objectives are (1) to assess environmental variables based on remote sensing time series and (2) to develop empirical models to forecast environmental variables. These objectives have been achieved in four applications which specific objectives are (1) assessing and mapping cotton crop phenological stages using spectral and phenometric analyses, (2) assessing and modeling fire seasonality in two different ecoregions by dynamic models, (3) forecasting forest fire risk on a pixel basis by dynamic models, and (4) assessing vegetation functioning based on temporal autocorrelation and phenometric analysis. The results of this dissertation show the usefulness of function fitting procedures to model AS1 and AS2. Phenometrics derived from function fitting procedure makes it possible to identify cotton crop phenological stages. Spectral analysis has demonstrated quantitatively the presence of one cycle in AS2 and two in AS1 and the unimodal and bimodal behaviour of fire seasonality in the Mediterranean and temperate ecoregions respectively. Autoregressive models has been used to characterize the dynamics of fire seasonality in two ecoregions and to forecasts accurately fire risk on a pixel basis. The usefulness of temporal autocorrelation to define and characterized land surface functioning has been demonstrated. And finally the “Optical Functional Types” concept has been proposed, in this approach pixels could be as temporal unities based on its temporal dynamics or functioning.

Joint interpretation of magnetotelluric, seismic, and well-log data in Hontomín (Spain)

Acknowledgements. This work is dedicated to the memory of Andrés Pérez-Estaún, brilliant scientist, colleague, and friend. The authors sincerely thank Ian Ferguson and an anonymous reviewer for their useful comments on the manuscript. Xènia Ogaya is currently supported in the Dublin Institute for Advanced Studies by a Science Foundation Ireland grant IRECCSEM (SFI grant 12/IP/1313). Juan Alcalde is funded by NERC grant NE/M007251/1, on interpretational uncertainty. Juanjo Ledo, Pilar Queralt and Alex Marcuello thank Ministerio de Economía y Competitividad and EU Feder Funds through grant CGL2014- 54118-C2-1-R. Funding for this Project has been partially provided by the Spanish Ministry of Industry, Tourism and Trade, through the CIUDEN-CSIC-Inst. Jaume Almera agreement (ALM-09-027: Characterization, Development and Validation of Seismic Techniques applied to CO2 Geological Storage Sites), the CIUDEN-Fundació Bosch i Gimpera agreement (ALM-09-009 Development and Adaptation of Electromagnetic techniques: Characterisation of Storage Sites) and the project PIERCO2 (Progress In Electromagnetic Research for CO2 geological reservoirs CGL2009-07604). The CIUDEN project is co-financed by the European Union through the Technological Development Plant of Compostilla OXYCFB300 Project (European Energy Programme for Recovery).

Biogenic opal estimation of deep-sea sediment cores from the Scotia Sea

We present biogenic opal flux records from two deep-sea sites in the Scotia Sea (MD07-3133 and MD07-3134) at decadal-scale resolution, covering the last glacial cycle. Besides conventional and time-consuming biogenic opal measuring methods, we introduce new biogenic opal estimation methods derived from sediment colour b*, wet bulk density, Si/Ti-count ratio, and Fourier transform infrared spectroscopy (FTIRS). All methods capture the biogenic opal amplitude, however, FTIRS - a novel method for marine sediment - yields the most reliable results. 230Th normalization data show strong differences in sediment focusing with intensified sediment focusing during glacial times. At MD07-3134 230Th normalized biogenic opal fluxes vary between 0.2 and 2.5 g/cm2/kyr. Our biogenic opal flux records indicate bioproductivity changes in the Southern Ocean, strongly influenced by sea ice distribution and also summer sea surface temperature changes. South of the Antarctic Polar Front, lowest bioproductivity occurred during the Last Glacial Maximum when upwelling of mid-depth water was reduced and sea ice cover intensified. Around 17 ka, bioproductivity increased abruptly, corresponding to rising atmospheric CO2 contents and decreasing seasonal sea ice coverage.

River and Hillslope Response to Localized Uplift Along a Left Bend of the San Andreas Fault, Santa Cruz Mountains, CA

Studying landscape evolution of the Earthís surface is difficult because both tectonic forces and surface processes control its response to perturbation, and ultimately, its shape and form. Researchers often use numerical models to study erosional response to deformation because there are rarely natural settings in which we can evaluate both tectonic activity and topographic response over appropriate time scales (103-105 years). In certain locations, however, geologic conditions afford the unique opportunity to study the relationship between tectonics and topography. One such location is along the Dragonís Back Pressure Ridge in California, where the landscape moves over a structural discontinuity along the San Andreas Fault and landscape response to both the initiation and cessation of uplift can be observed. In their landmark study, Hilley and Arrowsmith (2008) found that geomorphic metrics such as channel steepness tracked uplift and that hillslope response lagged behind that of rivers. Ideal conditions such as uniform vegetation density and similar lithology allowed them to view each basin as a developmental stage of response to uplift only. Although this work represents a significant step forward in understanding landscape response to deformation, it remains unclear how these results translate to more geologically complex settings. In this study, I apply similar methodology to a left bend along the San Andreas Fault in the Santa Cruz Mountains, California. At this location, the landscape is translated through a zone of localized uplift caused by the bend, but vegetation, lithology, and structure vary. I examine the geomorphic response to uplift along the San Andreas Fault bend in order to determine whether predicted landscape patterns can be observed in a larger, more geologically complex setting than the Dragonís Back Pressure Ridge. I find that even with a larger-scale and a more complex setting, geomorphic metrics such as channel steepness index remain useful tools for evaluating landscape evolution through time. Steepness indices in selected streams of study record localized uplift caused by the restraining bend, while hillslope adjustment in the form of landsliding occurs over longer time scales. This project illustrates that it is possible to apply concepts of landscape evolution models to complex settings and is an important contribution to the body of geomorphological study.

A generalized adaptive mathematical morphological filter for LIDAR data

Airborne Light Detection and Ranging (LIDAR) technology has become the primary method to derive high-resolution Digital Terrain Models (DTMs), which are essential for studying Earth's surface processes, such as flooding and landslides. The critical step in generating a DTM is to separate ground and non-ground measurements in a voluminous point LIDAR dataset, using a filter, because the DTM is created by interpolating ground points. As one of widely used filtering methods, the progressive morphological (PM) filter has the advantages of classifying the LIDAR data at the point level, a linear computational complexity, and preserving the geometric shapes of terrain features. The filter works well in an urban setting with a gentle slope and a mixture of vegetation and buildings. However, the PM filter often removes ground measurements incorrectly at the topographic high area, along with large sizes of non-ground objects, because it uses a constant threshold slope, resulting in "cut-off" errors. A novel cluster analysis method was developed in this study and incorporated into the PM filter to prevent the removal of the ground measurements at topographic highs. Furthermore, to obtain the optimal filtering results for an area with undulating terrain, a trend analysis method was developed to adaptively estimate the slope-related thresholds of the PM filter based on changes of topographic slopes and the characteristics of non-terrain objects. The comparison of the PM and generalized adaptive PM (GAPM) filters for selected study areas indicates that the GAPM filter preserves the most "cut-off" points removed incorrectly by the PM filter. The application of the GAPM filter to seven ISPRS benchmark datasets shows that the GAPM filter reduces the filtering error by 20% on average, compared with the method used by the popular commercial software TerraScan. The combination of the cluster method, adaptive trend analysis, and the PM filter allows users without much experience in processing LIDAR data to effectively and efficiently identify ground measurements for the complex terrains in a large LIDAR data set. The GAPM filter is highly automatic and requires little human input. Therefore, it can significantly reduce the effort of manually processing voluminous LIDAR measurements.

Increase in soil organic carbon by agricultural intensification in northern China

Acknowledgements. This research was supported by National Natural Science Foundation of China (no. 31370527 and 31261140367) and the National Science and Technology Support Program of China (no. 2012BAD14B01-2). The authors gratefully thank the Huantai Agricultural Station for providing of the Soil Fertility Survey data. We also thank Zheng Liang from China Agricultural University for the soil sampling and analysis in 2011. Thanks are extended to Jessica Bellarby for helpful discussion and suggestions.

Soil water dynamics and availability for citrus and peanut along a hillslope at the Sunjia Red Soil Critical Zone Observatory (CZO)

Acknowledgments This work was granted by the China-UK jointed Red Soil Critical Zone project from National Natural Science Foundation of China (NSFC: 41571130053; 41301233) and from Natural Environmental Research Council (NERC: Code: NE/N007611/1), and by the National Key Technology R&D Program of China (2011BAD31B04). We thank two anonymous reviewers for their constructive comments.

Patterns of carbon processing at the seafloor : The role of faunal and microbial communities in moderating carbon flows

Peer reviewed

Sediment budget and geomorphic process units of the Swedish Kärkevagge, link to ArcInfo E00 files

Geomorphic process units have been derived in order to allow quantification via GIS techniques at a catchment scale. Mass movement rates based on existing field measurements are employed in the budget calculations. In the Kärkevagge catchment, Northern Sweden, 80% of the area can be identified either as a source area for sediments or as a zone where sediments are deposited. The overall budget for the slopes beneath the rockwalls in the Kärkevagge is approximately 680 t/a whilst about 150 t a-1 are transported into the fluvial system.

Incisión fluvial del río Darro durante el periodo 1890-2010 en el sector urbano de Granada (Alambra-Valparaíso, España) determinada a partir de fotografías históricas

En este trabajo se calcula la tasa media de incisión fluvial del río Darro (Granada, España) durante el periodo 1890-2010 en su tramo urbano (sector Alhambra-Valparaíso). Para ello se han utilizado fotografías históricas en las que aparece dicho río, a partir de las cuales se ha podido determinar la posición del cauce en el momento en el que se realizaron las fotografías. La comparación con los escenarios actuales de tales imágenes ha permitido determinar la diferencia de altura del cauce a través de medidas de cotas absolutas realizadas mediante teodolito. Esta metodología ha permitido estimar de modo cuantitativo un índice de encajamiento vertical medio del río de 1,05 cm/año para el periodo histórico considerado.