Characteriation of Mediterranean Aleppo pine forest using low-density ALS data

Los espacios forestales son una fuente de servicios, tanto ambientales como económicos, de gran importancia para la sociedad. La caracterización de estos ambientes ha requerido tradicionalmente de un laborioso trabajo de campo. La aplicación de técnicas de teledetección ha proporcionado una visión más amplia a escala espacial y temporal, a la par que ha generado una reducción de los costes. La utilización de sensores óptico-pasivo multiespectrales y de sensores radar posibilita la estimación de parámetros forestales, si bien el desarrollo de sensores LiDAR, como el caso de los escáneres láser aeroportados (ALS), ha mejorado la caracterización tridimensional de la estructura de los bosques. La disponibilidad pública de dos coberturas LiDAR, generadas en el marco del Plan Nacional de Ortofotografía Aérea (PNOA), ha abierto nuevas líneas de investigación que permiten proporcionar información útil para la gestión forestal. La presente tesis utiliza datos LiDAR aeroportados de baja densidad para estimar diversas variables forestales, con ayuda de trabajo de campo, en masas forestales de Pino carrasco (Pinus halepensis Miller) en Aragón. La investigación aborda dos cuestiones relevantes como son la exploración de las metodologías más adecuadas para estimar variables forestales considerando escalas locales y regionales, teniendo en cuenta las posibles fuentes de error en el modelado; y, además, analiza la potencialidad de los datos LiDAR del PNOA para el desarrollo de aplicaciones forestales que valoricen las áreas forestales como recursos socio-económicos. La tesis se ha desarrollado según la modalidad de compendio de publicaciones, incluyendo cuatro trabajos que dan respuesta a los objetivos planteados. En primer lugar, se realiza un análisis comparativo de distintos modelos de regresión, paramétricos y no paramétricos, para estimar la pérdida de biomasa y las emisiones de CO2 en un incendio, mediante la utilización de datos LiDAR-PNOA y datos ópticos del satélite Landsat 8. En segundo lugar, se explora la idoneidad de distintos métodos de selección de variables para estimar biomasa total en masas de Pino carrasco utilizando datos LiDAR de baja densidad. En tercer lugar, se cuantificó y cartografió la biomasa residual forestal en el conjunto de masas de Pino carrasco de Aragón y se evaluó el efecto de diversas características de la tecnología LiDAR y de las variables ambientales en la precisión de los modelos. Finalmente, se analiza la transferibilidad temporal de modelos para estimar a escala regional siete variables forestales, utilizando datos LiDAR-PNOA multi-temporales. A este respecto, se compararon dos enfoques que permiten analizar la transferibilidad temporal: en primer lugar, el método directo ajusta un modelo para un determinado punto en el tiempo y estima las variables forestales para otra fecha; por otra parte, el método indirecto ajusta dos modelos diferentes para cada momento en el tiempo, estimando las variables forestales en dos fechas distintas. Los resultados obtenidos y las conclusiones derivadas de la investigación indican que la técnica basada en coeficientes de correlación de Spearman y el método de selección por todos los subconjuntos constituyen los métodos de selección de métricas LiDAR más apropiados para la modelización. El análisis de métodos de regresión para la estimación de variables forestales indicó que su idoneidad variaba de acuerdo con el tamaño y complejidad de la muestra. El método de regresión linear multivariante arrojó mejores resultados que los métodos no-paramétricos en el caso de muestras pequeñas. Por el contrario, el método Support Vector Machine produjo los mejores resultados con muestras grandes. El incremento de la densidad de puntos y de los valores de penetración de los pulsos LiDAR en el dosel, así como la presencia de ángulos de escaneo pequeños, incrementó la exactitud de los modelos. De forma similar, el incremento de la pendiente y la presencia de arbustos en el sotobosque implican una reducción en la exactitud de los modelos. En la estimación de variables forestales utilizando datos LiDAR multi-temporales, aunque la utilización del enfoque indirecto arrojó generalmente una mayor precisión en los modelos, se obtuvieron resultados similares con el enfoque directo, el cual constituye una alternativa óptima para reducir el tiempo de modelado y los costes de realización de trabajo de campo. La fusión de datos LiDAR y datos óptico-pasivos ha evidenciado la conveniencia de los métodos aplicados para cuantificar las emisiones de CO2 a la atmósfera generadas por un incendio. Esta metodología constituye una alternativa adecuada cuando no existen datos multi-temporales LiDAR. La estimación de variables de inventario forestal, así como de diversas fracciones de biomasa, como la biomasa total y la biomasa residual forestal, proporciona información valiosa para caracterizar las masas forestales mediterráneas de Pino carrasco y mejorar la gestión forestalForest ecosystems provide environmental and economic services of great importance to the society. The characterization of these environments has been traditionally accomplished with intense field work. In comparison, the application of remote sensing tools provides a greater overview over large spatial and temporal scales while minimizing costs. Although optical data and Synthetic Aperture Radar (SAR) allow estimating forest stand variables, the development of LiDAR sensors such as Airborne Laser Scanner (ALS) have improved three-dimensional characterization of forest structure. The availability of two ALS public data coverages for the Spanish territory, provided by the National Plan for Aerial Ortophotography (PNOA), opens new research opportunities to generate useful information for forest management. This PhD Thesis used low-density ALS-PNOA data to estimate different forest variables, with support in fieldwork, in the Aleppo pine (Pinus halepensis Miller) forests of Aragón region. The addressed research is relevant mainly for two reasons: first, the examination of suitable methodologies and error sources in forest stand variables prediction at local (small area) and regional scales (large area), and second, the application of ALS data to the characterization of forest areas as a socio-economic reservoir. This PhD Thesis is a compendium of four scientific papers, which sequentially answer the objectives established. Firstly, a comparative analysis of different parametric and non-parametric models was performed to estimate biomass losses and CO2 emissions using low-density ALS and Landsat 8 data in a burnt Aleppo pine forest. Secondly, we assess the suitability of variable selection methods when estimating total biomass in Aleppo pine forest stands using low-density ALS data. In the third manuscript, the quantification and mapping of forest residual biomass in Aleppo pine forest of Aragón region and the assessment of the effect of ALS and environmental variables in model accuracy were accomplished. Finally, the temporal transferability of seven forest stands attributes modelling using multi-temporal ALS-PNOA data in Aleppo pine forest at regional scale was explored. In this case, the temporal transferability was assessed comparing two methodologies; the direct and indirect approach. The first one fits a model for one point in time and estimates the forest variable for another point in time. The indirect approach adjusts two models in different points in time to estimate the forest variables in two different dates. The results derived from this research indicated that Spearman’s rank and All Subset Selection are the most appropriate methods in the ALS metrics selection step commonly applied in modelling. The suitability of the regression methods depends on the sample size and complexity. Thus, multivariate linear regression outperformed non-parametric methods with small samples while support vector machine was the most accurate method with larger samples. Model accuracy increased with higher point density and canopy pulse penetration, while decreasing with wider scan angles. Furthermore, the presence of steep slopes and shrub reduced model performance. In the case of forest stand variables prediction using multi-temporal ALS data, although the indirect approach produced generally a higher precision, the direct approach provided similar results, constituting a suitable alternative to reduce modelling time and fieldwork costs. The fusion of ALS and passive optical data have evidenced the suitability of this information for quantifying wildfire CO2 emissions to atmosphere, constituting a good alternative when multi-temporal ALS data is not available. The estimation of forest inventory variables as well as different biomass fractions, such as total biomass and forest residual biomass, provided valuable information to characterize Mediterranean Aleppo pine forests and improve forest management.<br /

Aboveground forest biomass estimation with Landsat and LiDAR data and uncertainty analysis of the estimates.

Landsat Thematic mapper (TM) image has long been the dominate data source, and recently LiDAR has offered an important new structural data stream for forest biomass estimations. On the other hand, forest biomass uncertainty analysis research has only recently obtained sufficient attention due to the difficulty in collecting reference data. This paper provides a brief overview of current forest biomass estimation methods using both TM and LiDAR data. A case study is then presented that demonstrates the forest biomass estimation methods and uncertainty analysis. Results indicate that Landsat TM data can provide adequate biomass estimates for secondary succession but are not suitable for mature forest biomass estimates due to data saturation problems. LiDAR can overcome TM?s shortcoming providing better biomass estimation performance but has not been extensively applied in practice due to data availability constraints. The uncertainty analysis indicates that various sources affect the performance of forest biomass/carbon estimation. With that said, the clear dominate sources of uncertainty are the variation of input sample plot data and data saturation problem related to optical sensors. A possible solution to increasing the confidence in forest biomass estimates is to integrate the strengths of multisensor data

The BIOMASS level 2 prototype processor : design and experimental results of above-ground biomass estimation

BIOMASS is ESA’s seventh Earth Explorer mission, scheduled for launch in 2022. The satellite will be the first P-band SAR sensor in space and will be operated in fully polarimetric interferometric and tomographic modes. The mission aim is to map forest above-ground biomass (AGB), forest height (FH) and severe forest disturbance (FD) globally with a particular focus on tropical forests. This paper presents the algorithms developed to estimate these biophysical parameters from the BIOMASS level 1 SAR measurements and their implementation in the BIOMASS level 2 prototype processor with a focus on the AGB product. The AGB product retrieval uses a physically-based inversion model, using ground-canceled level 1 data as input. The FH product retrieval applies a classical PolInSAR inversion, based on the Random Volume over Ground Model (RVOG). The FD product will provide an indication of where significant changes occurred within the forest, based on the statistical properties of SAR data. We test the AGB retrieval using modified airborne P-Band data from the AfriSAR and TropiSAR campaigns together with reference data from LiDAR-based AGB maps and plot-based ground measurements. For AGB estimation based on data from a single heading, comparison with reference data yields relative Root Mean Square Difference (RMSD) values mostly between 20% and 30%. Combining different headings in the estimation process significantly improves the AGB retrieval to slightly less than 20%. The experimental results indicate that the implemented retrieval scheme provides robust results that are within mission requirements

Aboveground forest biomass derived using multiple dates of WorldView-2 stereo-imagery : quantifying the improvement in estimation accuracy

The aim of this study was to investigate the capabilities of two date satellite-derived image-based point clouds (IPCs) to estimate forest aboveground biomass (AGB). The data sets used include panchromatic WorldView-2 stereo-imagery with 0.46 m spatial resolution representing 2014 and 2016 and a detailed digital elevation model derived from airborne laser scanning data. Altogether, 332 field sample plots with an area of 256 m(2) were used for model development and validation. Predictors describing forest height, density, and variation in height were extracted from the IPC 2014 and 2016 and used in k-nearest neighbour imputation models developed with sample plot data for predicting AGB. AGB predictions for 2014 (AGB(2014)) were projected to 2016 using growth models (AGB(Projected_2016)) and combined with the AGB estimates derived from the 2016 data (AGB(2016)). AGB prediction model developed with 2014 data was also applied to 2016 data (AGB(2016_pred2014)). Based on our results, the change in the 90(th) percentile of height derived from the WorldView-2 IPC was able to characterize forest height growth between 2014 and 2016 with an average growth of 0.9 m. Features describing canopy cover and variation in height derived from the IPC were not as consistent. The AGB(2016) had a bias of -7.5% (-10.6 Mg ha(-1)) and root mean square error (RMSE) of 26.0% (36.7 Mg ha(-1)) as the respective values for AGB(Projected_2016) were 7.0% (9.9 Mg ha(-1)) and 21.5% (30.8 Mg ha(-1)). AGB(2016_pred2014) had a bias of -19.6% (-27.7 Mg ha(-1)) and RMSE of 33.2% (46.9 Mg ha(-1)). By combining predictions of AGB(2016) and AGB(Projected_2016) at sample plot level as a weighted average, we were able to decrease the bias notably compared to estimates made on any single date. The lowest bias of -0.25% (-0.4 Mg ha(-1)) was obtained when equal weights of 0.5 were given to the AGB(Projected_2016) and AGB(2016) estimates. Respectively, RMSE of 20.9% (29.5 Mg ha(-1)) was obtained using equal weights. Thus, we conclude that combination of two date WorldView-2 stereo-imagery improved the reliability of AGB estimates on sample plots where forest growth was the only change between the two dates.Peer reviewe

Aboveground Forest Biomass Estimation with Landsat and LiDAR Data and Uncertainty Analysis of the Estimates

Landsat Thematic mapper (TM) image has long been the dominate data source, and recently LiDAR has offered an important new structural data stream for forest biomass estimations. On the other hand, forest biomass uncertainty analysis research has only recently obtained sufficient attention due to the difficulty in collecting reference data. This paper provides a brief overview of current forest biomass estimation methods using both TM and LiDAR data. A case study is then presented that demonstrates the forest biomass estimation methods and uncertainty analysis. Results indicate that Landsat TM data can provide adequate biomass estimates for secondary succession but are not suitable for mature forest biomass estimates due to data saturation problems. LiDAR can overcome TM’s shortcoming providing better biomass estimation performance but has not been extensively applied in practice due to data availability constraints. The uncertainty analysis indicates that various sources affect the performance of forest biomass/carbon estimation. With that said, the clear dominate sources of uncertainty are the variation of input sample plot data and data saturation problem related to optical sensors. A possible solution to increasing the confidence in forest biomass estimates is to integrate the strengths of multisensor data

Airborne Lidar Sampling Strategies to Enhance Forest Aboveground Biomass Estimation from Landsat Imagery

Accurately estimating aboveground biomass (AGB) is important in many applications, including monitoring carbon stocks, investigating deforestation and forest degradation, and designing sustainable forest management strategies. Although lidar provides critical three-dimensional forest structure information for estimating AGB, acquiring comprehensive lidar coverage is often cost prohibitive. This research focused on developing a lidar sampling framework to support AGB estimation from Landsat images. Two sampling strategies, systematic and classification-based, were tested and compared. The proposed strategies were implemented over a temperate forest study site in northern New York State and the processes were then validated at a similar site located in central New York State. Our results demonstrated that while the inclusion of lidar data using systematic or classification-based sampling supports AGB estimation, the systematic sampling selection method was highly dependent on site conditions and had higher accuracy variability. Of the 12 systematic sampling plans, R-2 values ranged from 0.14 to 0.41 and plot root mean square error (RMSE) ranged from 84.2 to 93.9 Mg ha(-1). The classification-based sampling outperformed 75% of the systematic sampling strategies at the primary site with R-2 of 0.26 and RMSE of 70.1 Mg ha(-1). The classification-based lidar sampling strategy was relatively easy to apply and was readily transferable to a new study site. Adopting this method at the validation site, the classification-based sampling also worked effectively, with an R-2 of 0.40 and an RMSE of 108.2 Mg ha(-1) compared to the full lidar coverage model with an R-2 of 0.58 and an RMSE of 96.0 Mg ha(-1). This study evaluated different lidar sample selection methods to identify an efficient and effective approach to reduce the volume and cost of lidar acquisitions. The forest type classification-based sampling method described in this study could facilitate cost-effective lidar data collection in future studies

Forest Aboveground Biomass Estimation Using Multi-Source Remote Sensing Data in Temperate Forests

Forests are a crucial part of global ecosystems. Accurately estimating aboveground biomass (AGB) is important in many applications including monitoring carbon stocks, investigating forest degradation, and designing sustainable forest management strategies. Remote sensing techniques have proved to be a cost-effective way to estimate forest AGB with timely and repeated observations. This dissertation investigated the use of multiple remotely sensed datasets for forest AGB estimation in temperate forests. We compared the performance of Landsat and lidar data—individually and fused—for estimating AGB using multiple regression models (MLR), Random Forest (RF) and Geographically Weight Regression (GWR). Our approach showed MLR performed similarly to GWR and both were better than RF. Integration of lidar and Landsat inputs outperformed either data source alone. However, although lidar provides valuable three-dimensional forest structure information, acquiring comprehensive lidar coverage is often cost prohibitive. Thus we developed a lidar sampling framework to support AGB estimation from Landsat images. We compared two sampling strategies—systematic and classification-based—and found that the systematic sampling selection method was highly dependent on site conditions and had higher model variability. The classification-based lidar sampling strategy was easy to apply and provides a framework that is readily transferable to new study sites. The performance of Sentinel-2 and Landsat 8 data for quantifying AGB in a temperate forest using RF regression was also tested. We modeled AGB using three datasets: Sentinel-2, Landsat 8, and a pseudo dataset that retained the spatial resolution of Sentinel-2 but only the spectral bands that matched those on Landsat 8. We found that while RF model parameters impact model outcomes, it is more important to focus attention on variable selection. Our results showed that the incorporation of red-edge information increased AGB estimation accuracy by approximately 6%. The additional spatial resolution improved accuracy by approximately 3%. The variable importance ranks in the RF regression model showed that in addition to the red- edge bands, the shortwave infrared bands were important either individually (in the Sentinel-2 model) or in band indices. With the growing availability of remote sensing datasets, developing tools to appropriately and efficiently apply remote sensing data is increasingly important