Open Access
Identifying forest harvesting practices: clear-cutting and thinning in diverse tree species using dense Landsat time series
(Elsevier, 2024-12-07) Giambelluca, Ana Laura; Hermosilla, Txomin; Álvarez-Mozos, Jesús; González de Audícana Amenábar, María; Ingeniería; Ingeniaritza; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD; Gobierno de Navarra / Nafarroako Gobernua
Forest monitoring plays a critical role in achieving sustainable forest management practices. The ability to identify ongoing harvesting activities is crucial for developing targeted strategies to maintain forest health. Traditional monitoring methods, which rely on field inventories, are often expensive and time-consuming. Remote sensing offers an interesting alternative, leveraging dense time series of satellite imagery and various algorithms for disturbance detection. This study presents and assesses a novel methodology for identifying forest harvesting practices (clear-cutting and thinning) using Continuous Change Detection and Classification (CCDC) algorithm, available in Google Earth Engine. The methodology comprises two steps. In the first step, performed at the pixel level, the CCDC algorithm was used to detect changes in the vegetation cover by considering Landsat 8 spectral bands, vegetation indices, and different combinations thereof. In the second step, two optimal thresholds were determined to identify forest harvesting practices based on the proportion of pixels flagged as change. This study was conducted in forest stands consisting of different conifer and broadleaf species. Accuracy was assessed using an independent set of photo-interpreted samples. The results indicated that the short-wave infrared 2 was the best individual band for forest harvesting practices identification, with an average F-score of 0.77 ± 0.06, overperforming vegetation indices. The combination of all spectral bands was the most effective to identify both clear-cuts and thinning (F-score = 0.85 ± 0.05). This combination was used to evaluate the accuracy of this approach for identifying harvesting practices over different tree species. Poplar (Populus sp.) had the highest identification rate (F-score = 0.99 ± 0.02), while black pine (Pinus nigra J.F. Arnold) stands had the lowest F-score (0.74 ± 0.05). These results highlight the ability to accurately identify forest harvesting practices even in heterogeneous forests with a high diversity of tree species using dense time series of Landsat imagery.
Open Access
Evaluation of 2D models for the prediction of surface depression storage using realistic reference values
(Wiley, 2016) Giménez Díaz, Rafael; Mezkiritz Barberena, Irantzu; Campo-Bescós, Miguel; Álvarez-Mozos, Jesús; González de Audícana Amenábar, María; Martínez de Aguirre Escobar, Alejandro; Casalí Sarasíbar, Javier; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
Depression storage (DS) is the maximum storage of precipitation and runoff in the soil surface at a given slope. The DS is determined by soil roughness that in agricultural soils is largely affected by tillage. The direct measurement of DS is not straightforward because of the natural permeability of the soil. Therefore, DS has generally been estimated from 2D/3D empirical relationships and numerical algorithms based on roughness indexes and height measurements of the soil surface, respectively. The objective of this work was to evaluate the performance of some 2D models for DS, using direct and reliable measurements of DS in an agricultural soil as reference values. The study was carried out in experimental microplots where DS was measured in six situations resulting from the combination of three types of tillage carried out parallel and perpendicular to the main slope. Those data were used as reference to evaluate four empirical models and a numerical method. Longitudinal altitudinal profiles of the relief were obtained by a laser profilometer. Infiltration measurements were carried out before and after tillage. The DS was largely affected by tillage and its direction. Highest values of DS are found on rougher surfaces mainly when macroforms cut off the dominant slope. The empirical models had a limited performance while the numerical method was the most effective, even so, with an important variability. In addition, a correct hydrological management should take into account that each type of soil tillage affects infiltration rate differently.
Open Access
The added value of stratified topographic correction of multispectral images
(MDPI, 2016) Sola Torralba, Ion; González de Audícana Amenábar, María; Álvarez-Mozos, Jesús; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
Satellite images in mountainous areas are strongly affected by topography. Different studies demonstrated that the results of semi-empirical topographic correction algorithms improved when a stratification of land covers was carried out first. However, differences in the stratification strategies proposed and also in the evaluation of the results obtained make it unclear how to implement them. The objective of this study was to compare different stratification strategies with a non-stratified approach using several evaluation criteria. For that purpose, Statistic-Empirical and Sun-Canopy-Sensor + C algorithms were applied and six different stratification approaches, based on vegetation indices and land cover maps, were implemented and compared with the non-stratified traditional option. Overall, this study demonstrates that for this particular case study the six stratification approaches can give results similar to applying a traditional topographic correction with no previous stratification. Therefore, the non-stratified correction approach could potentially aid in removing the topographic effect, because it does not require any ancillary information and it is easier to implement in automatic image processing chains. The findings also suggest that the Statistic-Empirical method performs slightly better than the Sun-Canopy-Sensor + C correction, regardless of the stratification approach. In any case, further research is necessary to evaluate other stratification strategies and confirm these results.
Open Access
Evaluación multitemporal de métodos de corrección topográfica mediante el uso de imágenes sintéticas multiespectrales
(Asociación Española de Teledetección, 2014) Sola Torralba, Ion; Álvarez-Mozos, Jesús; González de Audícana Amenábar, María; Torres Escribano, José Luis; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
En este trabajo se presentan los resultados de la evaluación multitemporal de varios métodos de corrección topográfica (TOC), cuya bondad se determina de forma cuantitativa mediante el uso de imágenes sintéticas multiespectrales simuladas para diferentes fechas de adquisición a lo largo del año. Para cada fecha se generan dos imágenes sintéticas, una considerando el relieve real (imagen SR), y otra el relieve horizontal (imagen SH). Las imágenes SR se corrigen utilizando distintos TOC y estas imágenes corregidas se comparan con la corrección ideal (imagen SH) mediante el índice de similitud estructural (SSIM). Los valores de SSIM nos permiten evaluar la eficacia de cada corrección para distintas fechas, es decir, para distintos ángulos de elevación solar.
Open Access
Validation of a simplified model to generate multispectral synthetic images
(MDPI, 2015) Sola Torralba, Ion; González de Audícana Amenábar, María; Álvarez-Mozos, Jesús; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
A new procedure to assess the quality of topographic correction (TOC) algorithms applied to remote sensing imagery was previously proposed by the authors. This procedure was based on a model that simulated synthetic scenes, representing the radiance an optical sensor would receive from an area under some specific conditions. TOC algorithms were then applied to synthetic scenes and the resulting corrected scenes were compared with a horizontal synthetic scene free of topographic effect. This comparison enabled an objective and quantitative evaluation of TOC algorithms. This approach showed promising results but had some shortcomings that are addressed herein. First, the model, originally built to simulate only broadband panchromatic scenes, is extended to multispectral scenes in the visible, near infrared (NIR), and short wave infrared (SWIR) bands. Next, the model is validated by comparing synthetic scenes with four Satellite pour l'Observation de la Terre 5 (SPOT5) real scenes acquired on different dates and different test areas along the Pyrenees mountain range (Spain). The results obtained show a successful simulation of all the spectral bands. Therefore, the model is deemed accurate enough for its purpose of evaluating TOC algorithms.
Open Access
Synthetic images for evaluating topographic correction algorithm
(IEEE, 2013) Sola Torralba, Ion; González de Audícana Amenábar, María; Álvarez-Mozos, Jesús; Torres Escribano, José Luis; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak; Gobierno de Navarra / Nafarroako Gobernua
In the last years, many topographic correction (TOC) methods have been proposed to correct the illumination differences between the areas observed by optical remote sensors. Although the available number of TOC methods is high, the evaluation of their performance generally relies on the existence of precise land-cover information, and a standardized and objective evaluation procedure has not been proposed yet. In this paper, we propose an objective procedure to assess the accuracy of these TOC methods on the basis of simulated scenes, i.e., synthetically generated images. These images represent the radiance an optical sensor would receive under specific geometric and temporal acquisition conditions and assuming a certain land-cover type. A simplified method for creating synthetic images using the stateof- the-art irradiance models is proposed, both considering the real topography of a certain area [synthetic real (SR) image] or considering the relief of this area as being completely flat [synthetic horizontal image (SH)]. The comparison between the corrected image obtained by applying a TOC method to the SR and SH images of the same area, allows assessing the performance of each TOC algorithm. This comparison is quantitatively carried out using the structural similarity index. The proposed TOC evaluation procedure is applied to a specific case study in northern Spain to explain its implementation and demonstrate its potential. The procedure proposed in this paper could be also used to assess the behavior of TOC methods operating under different scenarios considering diverse topographic, geometrical, and temporal acquisition configurations.
Open Access
Multitemporal evaluation of topographic correction algorithms using synthetic images
(SPIE, 2012) Sola Torralba, Ion; Álvarez-Mozos, Jesús; González de Audícana Amenábar, María; Torres Escribano, José Luis; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
Land cover classification and quantitative analysis of multispectral data in mountainous regions is considerably hampered by the influence of topography on the spectral response pattern. In the last years, different topographic correction (TOC) algorithms have been proposed to correct illumination differences between sunny and shaded areas observed by optical remote sensors. Although the available number of TOC methods is high, the evaluation of their performance usually relies on the existence of precise land cover information, and a standardised and objective evaluation procedure has not been proposed yet. Besides, previous TOC assessment studies only considered a limited set of illumination conditions, normally assuming favourable illumination conditions. This paper presents a multitemporal evaluation of TOC methods based on synthetically generated images in order to evaluate the influence of solar angles on the performance of TOC methods. These synthetic images represent the radiance an optical sensor would receive under specific geometric and temporal acquisition conditions and assuming a certain land-cover type. A method for creating synthetic images using state-of-the-art irradiance models has been tested for different periods of the year, which entails a variety of solar angles. Considering the real topography of a specific area a Synthetic Real image (SR) is obtained, and considering the relief of this area as being completely flat a Synthetic Horizontal image (SH) is obtained. The comparison between corrected image obtained applying a TOC method to SR image and SH image of the same area, i.e. considered the ideal correction, allows assessing the performance of each TOC algorithm.
Open Access
Inter-comparison of atmospheric correction methods on Sentinel-2 images applied to croplands
(IEEE, 2018) Sola Torralba, Ion; Álvarez-Mozos, Jesús; González de Audícana Amenábar, María; Ingeniería; Ingeniaritza
Atmospheric correction of high resolution satellite scenery is a necessary preprocessing step for applications where bottom of atmosphere (BOA) reflectances are needed. The selection of the best atmospheric correction method to use on images acquired from new platforms, such as Sentinel-2, is essential to provide accurate BOA reflectances. In this work the performance of three atmospheric correction methods (6S, MAJA and SEN2COR) applied to Sentinel-2 scenes are compared by evaluating the resultant spectral signatures of six crop types on two specific dates, and their NDVI time series along a complete year. Although SEN2COR introduced greater corrections, especially in the infrared bands, the results suggest a varying performance of the methods depending on the land cover and the atmospheric conditions. Further research, particularly incorporating ground truth data, is recommended to rigorously validate the different atmospheric methods.
Open Access
Influence of surface roughness spatial variability and temporal dynamics on the retrieval of soil moisture from SAR observations
(MDPI, 2009) Álvarez-Mozos, Jesús; Verhoest, Niko E. C.; Larrañaga Urien, Arantzazu; Casalí Sarasíbar, Javier; González de Audícana Amenábar, María; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
Radar-based surface soil moisture retrieval has been subject of intense research during the last decades. However, several difficulties hamper the operational estimation of soil moisture based on currently available spaceborne sensors. The main difficulty experienced so far results from the strong influence of other surface characteristics, mainly roughness, on the backscattering coefficient, which hinders the soil moisture inversion. This is especially true for single configuration observations where the solution to the surface backscattering problem is ill-posed. Over agricultural areas cultivated with winter cereal crops, roughness can be assumed to remain constant along the growing cycle allowing the use of simplified approaches that facilitate the estimation of the moisture content of soils. However, the field scale spatial variability and temporal variations of roughness can introduce errors in the estimation of soil moisture that are difficult to evaluate. The objective of this study is to assess the impact of roughness spatial variability and roughness temporal variations on the retrieval of soil moisture from radar observations. A series of laser profilometer measurements were performed over several fields in an experimental watershed from September 2004 to March 2005. The influence of the observed roughness variability and its temporal variations on the retrieval of soil moisture is studied using simulations performed with the Integral Equation Model, considering different sensor configurations. Results show that both field scale roughness spatial variability and its temporal variations are aspects that need to be taken into account, since they can introduce large errors on the retrieved soil moisture values.
Open Access
Assessment of atmospheric correction methods for Sentinel-2 images in Mediterranean landscapes
(Elsevier, 2018) Sola Torralba, Ion; García-Martín, Alberto; Sandonís Pozo, Leire; Álvarez-Mozos, Jesús; González de Audícana Amenábar, María; Ingeniería; Ingeniaritza
Atmospheric correction of optical satellite imagery is an essential pre-processing for modelling biophysical variables, multi-temporal analysis, and digital classification processes. Sentinel-2 products available for users are distributed by the European Space Agency (ESA) as Top Of Atmosphere reflectance values in cartographic geometry (Level-1C product). In order to obtain Bottom Of Atmosphere reflectance images (Level-2A product) derived from this Level-1C products, ESA provides the SEN2COR module, which is implemented in the Sentinel Application Platform. Alternatively, ESA recently distributes Level-2A products processed by SEN2COR with a default configuration. On the other hand, the conversion from Level-1C to Level-2A product can be generated using alternative atmospheric correction methods, such as MAJA, 6S, or iCOR. In this context, this paper aims to evaluate the quality of Level-2A products obtained through different methods in Mediterranean shrub and grasslands by comparing data obtained from Sentinel-2 imagery with field spectrometry data. For that purpose, six plots with different land covers (asphalt, grass, shrub, pasture, and bare soil) were analyzed, by using synchronous imagery to fieldwork (from July to September 2016). The results suggest the suitability of the applied atmospheric corrections, with coefficients of determination higher than 0.90 and root mean square error lower than 0.04 achieving a relative error in bottom of atmosphere reflectance of only 2–3%. Nevertheless, minor differences were observed between the four tested methods, with slightly varying results depending on the spectral band and land cover.