Search NASA⌕ Search

Engineering topics

Bruce D. Cook

Publications and source records attributed to Bruce D. Cook.

9 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Three-dimensional estimation of deciduous forest canopy structure and leaf area using multi-directional, leaf-on and leaf-off airborne lidar data

Airborne laser scanning (ALS) has been widely used to map gap probability and leaf area index (LAI) distribution at plot and landscape scales. As an indirect measurement, most ALS methods to estimate LAI combine waveform or point density information with supporting field measurements such as the leaf angle distribution, gap probability, or direct LAI measures. The development of a more independent estimation approach would facilitate more widespread use of existing ALS data to investigate patterns of forest structure and build realistic 3-D vegetation scenes to simulate remote sensing imagery and energy balance. Here, we develop a data processing workflow (named PVlad) using ALS point cloud apparent reflectance to estimate LAI and voxel-based leaf area density (LAD), aiming to reduce the need for associated field measurements such as the gap probability. The adaptation of the path volume (PV) concept derived from apparent reflectance integrates information from multi-directional ALS pulses, and quantifies the percentage exploration of each voxel for classification and occlusion correction, such that rigorous volumetric sampling approaches can be developed to derive LAI and LAD. The PVlad workflow was applied to discrete-return lidar data (Riegl VQ480i) acquired by NASA Goddard's LiDAR, Hyperspectral and Thermal Imager (G-LiHT) Airborne Imager during leaf-on (summer) and leaf-off (spring) conditions at the Smithsonian Environmental Research Center (SERC). The estimates of LAI and LAD captured structural differences between mature, logged, and intermediate-aged stands over eight deciduous forest plots. The derived LAI values were compared to field litter collection measurements, and the derived LAD vertical distribution was compared to the output of the VoxLAD model using terrestrial laser scan (TLS) field survey data. Using voxel sizes ranging from 0.5 m to 5 m, overall LAI estimation showed linear fitting coefficient bias and for 1 and 2 m voxel sizes, and vertical LAD distribution showed strong correlation with and for 0.5 and 1m voxel sizes. For every forest stand, upper-canopy LAD had a low variance for voxel sizes of ≤ . Application of PVlad to the G-LiHT and other similar ALS data archives enables the development of fine-resolution LAI map products, including voxelization of LAD for ecosystem science and radiative transfer simulations of remote sensing imagery or surface energy balance.

Tiangang Yin↗
NASA NTRS2021

Tracking the Rates and Mechanisms of Canopy Damage and Recovery Following Hurricane Maria Using Multitemporal Lidar Data

Hurricane Maria, a Category 4 storm, snapped and uprooted canopy trees, removed large branches, and defoliated vegetation across Puerto Rico. The magnitude of forest damages and the rates and mechanisms of forest recovery following Maria provide important benchmarks for understanding the ecology of extreme events. We used airborne Lidar data acquired before (2017) and after Maria (2018, 2020) to quantify landscape-scale changes in forest structure along a 439-ha elevational gradient (100–800 m) in the Luquillo Experimental Forest. Damages from Maria were widespread, with 73% of the study area losing ≥ 1 m in canopy height (mean = −7.1 m). Taller forests at lower elevations suffered more damage than shorter forests above 600 m. Yet only 13.5% of the study area had canopy heights ≤ 2 m in 2018, a typical threshold for forest gaps, highlighting the importance of damaged trees and advanced regeneration on post-storm forest structure. Heterogeneous patterns of regrowth and recruitment yielded shorter and more open forests by 2020. Nearly 45% of forests experienced initial height loss > 1 m (2017–2018) followed by rapid height gain > 1 m (2018–2020), whereas 21.6% of forests with initial height losses showed little or no height gain, and 17.8% of forests exhibited no height changes larger than ± 1 m in either period. Canopy layers < 10 m tall accounted for most increases in canopy height and fractional cover between 2018 and 2020, with gains split evenly between height growth and lateral crown expansion by surviving individuals. These findings benchmark rates of gap formation, crown expansion, and canopy closure following hurricane damage and highlight the diversity of ecosystem impacts from heterogeneous spatial patterns and vertical stratification of forest regrowth following a major disturbance event.

forest structure↗
NASA NTRS2021

Tracking the Rates and Mechanisms of Canopy Damage and Recovery Following Hurricane Maria Using Multitemporal Lidar Data

Hurricane Maria, a Category 4 storm, snapped and uprooted canopy trees, removed large branches, and defoliated vegetation across Puerto Rico. The magnitude of forest damages and the rates and mechanisms of forest recovery following Maria provide important benchmarks for understanding the ecology of extreme events. We used airborne Lidar data acquired before (2017) and after Maria (2018, 2020) to quantify landscape-scale changes in forest structure along a 439-ha elevational gradient (100–800 m) in the Luquillo Experimental Forest. Damages from Maria were widespread, with 73% of the study area losing ≥ 1 m in canopy height (mean = −7.1 m). Taller forests at lower elevations suffered more damage than shorter forests above 600 m. Yet only 13.5% of the study area had canopy heights ≤ 2 m in 2018, a typical threshold for forest gaps, highlighting the importance of damaged trees and advanced regeneration on post-storm forest structure. Heterogeneous patterns of regrowth and recruitment yielded shorter and more open forests by 2020. Nearly 45% of forests experienced initial height loss > 1 m (2017–2018) followed by rapid height gain > 1 m (2018–2020), whereas 21.6% of forests with initial height losses showed little or no height gain, and 17.8% of forests exhibited no height changes larger than ± 1 m in either period. Canopy layers < 10 m tall accounted for most increases in canopy height and fractional cover between 2018 and 2020, with gains split evenly between height growth and lateral crown expansion by surviving individuals. These findings benchmark rates of gap formation, crown expansion, and canopy closure following hurricane damage and highlight the diversity of ecosystem impacts from heterogeneous spatial patterns and vertical stratification of forest regrowth following a major disturbance event.

forest structure↗
NASA NTRS2021

Assessing impacts of canopy 3D structure on chlorophyll fluorescence radiance and radiative budget of deciduous forest stands using DART

Although remote sensing (RS) of solar-induced chlorophyll fluorescence (SIF) is increasingly used as a valuable source of information about vegetation photosynthetic activity, the RS SIF observations are significantly influenced by canopy-specific structural features (i.e., canopy architecture including leaf area index and presence of woody components), atmospheric conditions during their acquisition (e.g., proportion of direct and diffuse irradiance) and observational geometric configurations (e.g., sun and viewing directions). Radiative transfer (RT) models have the potential to provide a better understanding of the canopy structural effects on the SIF emission and RS signals. Here, we used the DART model to assess the daily influence, from morning to evening, of forest 3D architecture on SIF nadir radiance, emission, escape factor and nadir yield of eight 100 m × 100 m forest study plots established in a temperate deciduous forest of the Smithsonian Environmental Research Center (Edgewater, MD, USA). The 3D architecture of each plot was derived from airborne LiDAR. DART simulations of these 3D forest plots and their 1D (i.e., vertical profile of sun-adapted and shade-adapted leaves) and 0D (i.e., homogeneous layer of sun-adapted leaves above an homogeneous layer of shade-adapted leaves) abstractions were compared to assess the relative errors (ε1D−3D and ε0D−3D) associated with horizontal and vertical structural heterogeneity, respectively. Forest 3D structure, especially horizontal heterogeneity, had a great influence on forest nadir SIF radiance, resulting in ε1D−3D up to 55% at 8:00 and 18:00 (i.e., for oblique sun directions). The key indicators of this impact, in the descending order of importance, were the SIF escape factor (ε1D−3D up to 40%), the attenuation of incident photosynthetically active radiation (ε1D−3D less than 5%), and the SIF emission yield (ε1D−3D less than 2%). The influence of forest architecture on the nadir SIF escape factor and SIF yield (ε1D−3D up to 40%) varied over time, with differences in forest stand structure, and per spectral domain, being always larger between 640 and 700 nm than between 700 and 850 nm. In addition, woody elements demonstrated a large influence on forest SIF radiance due to their “shading” effect (ε up to 17%) and their “blocking” effect (ε ≈ 10%), both of them higher for far-red than for red SIF. These results underline the importance of 3D forest canopy architecture, especially 2D heterogeneity, and inclusion of woody elements in RT modeling used for interpretation of the RS SIF signal, and subsequently for the estimation of gross primary production and detection of vegetation stress.

Omar Regaieg↗
NASA NTRS2021

Mapping Boreal Forest Spruce Beetle Health Status at the Individual Crown Scale Using Fused Spectral and Structural Data

The frequency and severity of spruce bark beetle outbreaks are increasing in boreal forests leading to widespread tree mortality and fuel conditions promoting extreme wildfire. Detection of beetle infestation is a forest health monitoring (FHM) priority but is hampered by the challenges of detecting early stage (“green”) attack from the air. There is indication that green stage might be detected from vertical gradients of spectral data or from shortwave infrared information distributed within a single crown. To evaluate the efficacy of discriminating “non-infested”, “green”, and “dead” health statuses at the landscape scale in Alaska, USA, this study conducted spectral and structural fusion of data from: (1) Unoccupied aerial vehicle (UAV) multispectral (6 cm) + structure from motion point clouds (~700 pts per sq. m); and (2) Goddard Lidar Hyperspectral Thermal (G-LiHT) hyperspectral (400 to 1000 nm, 0.5 m) + SWIR-band lidar (~32 pts per sq.m). We achieved 78% accuracy for all three health statuses using spectral + structural fusion from either UAV or G-LiHT and 97% accuracy for non-infested/dead using G-LiHT. We confirm that UAV 3D spectral (e.g., greenness above versus below median height in crown) and lidar apparent reflectance metrics (e.g., mean reflectance at 99th percentile height in crown), are of high value, perhaps capturing the vertical gradient of needle degradation. In most classification exercises, UAV accuracy was lower than G-LiHT indicating that collecting ultra-high spatial resolution data might be less important than high spectral resolution information. While the value of passive optical spectral information was largely confined to the discrimination of non-infested versus dead crowns, G-LiHT hyperspectral band selection (~400, 675, 755, and 940 nm) could inform future FHM mission planning regarding optimal wavelengths for this task. Interestingly, the selected regions mostly did not align with the band designations for our UAV multispectral data but do correspond to, e.g., Sentinel-2 red edge bands, suggesting a path forward for moderate scale bark beetle detection when paired with suitable structural data.

Janice Cessna↗
NASA NTRS2021

Classifying Forest Type in the National Forest Inventory Context with Airborne Hyperspectral and Lidar Data

Forest structure and composition regulate a range of ecosystem services, including biodiversity, water and nutrient cycling, and wood volume for resource extraction. Forest type is an important metric measured in the US Forest Service Forest Inventory and Analysis (FIA) program, the national forest inventory of the USA. Forest type information can be used to quantify carbon and other forest resources within specific domains to support ecological analysis and forest management decisions, such as managing for disease and pests. In this study, we developed a methodology that uses a combination of airborne hyperspectral and lidar data to map FIA-defined forest type between sparsely sampled FIA plot data collected in interior Alaska. To determine the best classification algorithm and remote sensing data for this task, five classification algorithms were tested with six different combinations of raw hyperspectral data, hyperspectral vegetation indices, and lidar-derived canopy and topography metrics. Models were trained using forest type information from 632 FIA subplots collected in interior Alaska. Of the thirty model and input combinations tested, the random forest classification algorithm with hyperspectral vegetation indices and lidar-derived topography and canopy height metrics had the highest accuracy (78% overall accuracy). This study supports random forest as a powerful classifier for natural resource data. It also demonstrates the benefits from combining both structural (lidar) and spectral (imagery) data for forest type classification.

random forest↗
NASA NTRS2020

Mapping tall shrub biomass in Alaska at landscape scale using structure-from-motion photogrammetry and lidar

Warming in arctic and boreal regions is increasing shrub cover and biomass. In southcentral Alaska, willow (Salix spp.) and alder (Alnus spp.) shrubs grow taller than many tree species and account for a substantial proportion of aboveground biomass, yet they are not individually measured as part of the operational Forest Inventory and Analysis (FIA) Program. The goal of this research was to test methods for landscape-scale mapping of tall shrub biomass in upper montane and subalpine environments using FIA-type plot measurements (n = 51) and predictor variables from imagery-based structure-from-motion (SfM) and airborne lidar. Specifically, we compared biomass models constructed from imagery acquired by unmanned aerial vehicle (UAV; ~1.7 cm pixels), imagery from the NASA Goddard's Lidar, Hyperspectral, and Thermal Airborne Imager (G-LiHT; ~3.1 cm pixels), and concomitant G-LiHT small-footprint lidar. Tall shrub biomass was most accurately predicted at 5 m resolution (R^2 = 0.81, RMSE = 1.09 kg/sq. m) using G-LiHT SfM color and structure variables. Lidar-only models had lower precision (R^2 = 0.74, RMSE = 1.26 kg/sq. m), possibly due to reduced model information content from variable multicollinearity or lower data density. Separate models for upper montane zones with trees and shrubs and subalpine zones with only shrubs were always chosen over single models based on minimization of Akaike's Information Criterion, indicating the need for variable sets robust to overhanging tree canopy. Decreasing point density from UAV (5000–8000 pts./sq. m) to the G-LiHT SfM point cloud (500–2000 pts./sq. m) had little impact on model fit, suggesting that high-resolution airborne imagery can extend SfM approaches well beyond line-of-sight restrictions for UAV platforms. Overall, our results confirmed that SfM from high-resolution imagery is a viable approach to estimate shrub biomass in the boreal region, especially when an existing lidar terrain model and local field calibration data are available to quantify uncertainty in the SfM point cloud and landscape-scale estimates of shrub biomass.

Michael Alonzo↗
NASA NTRS2019

Boreal canopy surfaces from spaceborne stereogrammetry

Surface elevation estimates from high resolution spaceborne image (HRSI) stereogrammetry are used to examine fine-scaled structure of boreal forest canopies. These data can depict detailed spatial patterns of vertical forest structure at remote sites across the circumpolar domain where these estimates would otherwise be unavailable. This work examines where these estimates are most effective at describing vertical forest structure to explain which canopy surfaces they represent. We evaluated the variation in canopy surface estimates captured from four general types of HRSI digital surface models (DSMs) across the full range of boreal canopy cover. These DSMs, classified into 4 types by grouping them according to the acquisition's (1) sun elevation angle (low or high) and (2) seasonality-driven ground surface condition (snow presence/absence), vary with acquisition characteristics and the details of this variation continues to be studied. We explored some of this variation by comparing the distributions of differences in boreal canopy percentile heights derived from reference small footprint lidar in Tanana Valley, Alaska with canopy surface elevations derived from these 4 types of HRSI DSMs. We examined how canopy surface estimates from HRSI DSMs differ according to acquisition characteristics and canopy cover, and ultimately which canopy surfaces are represented in these DSMs. Our results help clarify which boreal canopy surfaces are representative of those captured with HRSI DSMs. They show that in the Tanana Valley (1) DSMs grouped by sun elevation angle and ground surface condition provide different surface estimates of boreal canopies; (2) the two DSM types that appear to most differently capture boreal forest canopy surfaces are DSMs from snow-free images acquired at sun elevation angles <30° (Low sun elev. & snow-free) and those with snow-cover at sun elevation angles ≥30° (High sun elev. & snow-free); (3) DSMs with snow most often do not capture upper canopy surfaces; (4) the “Low sun elev. & snow-free” DSMs resolve surfaces that are most representative of upper canopy surfaces (dense forests >60% cover, 70th–80th percentile heights); and (5) in the most dense forests (>80% cover) where canopy gaps are least likely to bias downward the average surface estimates, the snow-free DSM types are representative of 70th - 80th percentile heights (“Low sun elev. & snow-free”) and 60th–70th percentile heights (“High sun elev. & snow-free”). The combination of horizontal structure (canopy cover) and acquisition characteristics affect the boreal vertical structure (canopy surface height) estimates from spaceborne stereogrammetry. These effects should be considered when analyzing products derived from HRSI DSMs, and as part of a comprehensive approach to spaceborne remote sensing of circumpolar boreal forests.

forest structure↗