Source: Applied GeoSolutions, LLC submitted to NRP
OPERATIONALIZING MULTISCALE SAR METRICS FOR RAPID FOREST DISTURBANCE ASSESSMENT
Sponsoring Institution
National Institute of Food and Agriculture
Project Status
COMPLETE
Funding Source
Reporting Frequency
Annual
Accession No.
1011291
Grant No.
2016-33610-25687
Cumulative Award Amt.
$600,000.00
Proposal No.
2016-03968
Multistate No.
(N/A)
Project Start Date
Sep 1, 2016
Project End Date
May 31, 2019
Grant Year
2016
Program Code
[8.1]- Forests & Related Resources
Recipient Organization
Applied GeoSolutions, LLC
87 Packers Falls Road
Durham,NH 03824
Performing Department
(N/A)
Non Technical Summary
Demand for mapping forest disturbance and structural characteristics exists from USDA and Federal programs, industry, and international programs such as REDD. Applications range from operational forestry, assessing carbon sequestration, monitoring wildlife habitat, weather and disaster response, infestations, gauging forest productivity, and evaluating programs such as Environmental Quality Incentives Program (EQIP), Easements, FIA, and state management plans. There is an opportunity to operationalize multiscale Synthetic Aperture Radar (SAR) data to generate metrics of stand structure and spatiotemporal changes. Specifically, this SBIR Phase II is to prototype an automated system and scale products for above ground biomass, forest stand height, crown canopy cover, and disturbance detection using multiscale SAR. During Phase II we build on previous work and leverage ongoing partnerships with USDA USFS and space agencies (NASA, USGS, JAXA, ISRO). Our focused Phase II builds on Phase I case studies by conducting a coordinated campaign of near-simultaneous collection of field data, SAR, and Lidar in partnership with USFS; adds the recently launched operational Sentinel C-band satellites; and scales products to larger areas in a robust and automated approach. We considered recommendations from Phase 1 reviews, Science Advisory Panel feedback, and the most promising outcomes from Phase 1. The long-term (Phase III) vision is to 1.) provide disturbance mapping services, 2.) build and offer robust Monitoring, Reporting, and Verification (MRV) forest metrics derived from multiscale SAR, and 3.) develop Public Private Partnerships (PPP) to support regional, Federal (USFS, NASA, EPA), and international programs centered on monitoring forest disturbances and characterizing ecosystems services.
Animal Health Component
33%
Research Effort Categories
Basic
33%
Applied
33%
Developmental
34%
Classification

Knowledge Area (KA)Subject of Investigation (SOI)Field of Science (FOS)Percent
12306202060100%
Goals / Objectives
The goal of the SBIR is to develop innovative technology to assess and monitor our valuable forest resources and their characteristics. This SBIR also addresses NIFA Societal Challenge Area Climate Change by developing standardized metrics and easy-to-execute methods that will contribute to carbon science and reducing greenhouse gas emissions through programs such as Forest Inventory and Analysis (FIA) National Program, national USFS management plan for sustainable utilization, and international programs such as Reducing Emissions from Deforestation and Degradation (REDD), SilvaCarbon, Global Forest Observations Initiative (GFOI), and Global Observation of Forest and Land Cover Dynamics (GOFC-GOLD).Specifically, this SBIR Phase II is to prototype an automated system and scale products for above ground biomass, forest stand height, crown canopy cover, and disturbance detection using multiscale SAR.
Project Methods
Synthetic Aperture Radar are active sensors that measure echo through clouds with sensitivity to vegetation structure, geometry, and moisture content (dielectric properties). Wavelength frequency, polarization, and viewing geometry are key factors determining choice of SAR data for a given application. In short, longer wavelengths (smaller frequency) have greater ability to 'penetrate' a given target (fig 2). SAR signal can transmit horizontal (H) or vertical (V) field vectors, and receive either horizontal (H) or vertical (V) return signals, or both. The incidence angle (θ) can be defined as the angle formed by the radar beam and a line perpendicular to the surface. By modeling the given responses of a target with known imaging parameters a detailed description of scattering can be generated. At L-band (1.25 GHz; 23cm wavelength), radar waves penetrate into the forest canopy through gaps and areas with relatively low density of scatterers (sparse medium effect) and scatter back from various components of the forest such as the crown (leaves and branches) and stems (Saatchi and McDonald 1997). The co- and cross-polarized radar measurements contain different information about the orientation and structure of forest canopy and tree stems within the resolution cell of radar images (Ulaby and Elachi 1990). Thus, SAR has particular advantages over traditional optical satellite imaging at the same time creates barriers to use due to relatively more complexities and requirements for data handling.There is tremendous opportunity to take advantage of SAR archives and new SAR imagery. SAR observations have been used in many forestry research applications and academic institutions. The most common applications identify relationships between SAR backscatter and the vegetation structural attributes in one of two ways, via empirical models or mechanistic models. For example, Saatchi et al. (2007) used an empirical approach to map biomass and fuel loads in Yellowstone using L- and C- band radar. Mitchard et al (2009) used co-polarized L-band data to map above ground biomass across several sites of woody vegetation up to 150 Mg ha-1 before relationships became noisy from saturation. At high biomass L-band can become saturated (Saatchi et al 2011); Zheng et al (2008) showed that the average biomass in the northeast around 120 Mg-ha.As part of our Phase II we will continue to gauge limitations to mapping Above Ground Biomass (AGB) at multiple sites including the entire northeast USA with L-band from ALOS-1 and ALOS-2 and fusion with C-band Sentinel-1A. Hame et al (2013) found L-band effective in mapping AGB in Laos with dense tropical forest, but noted many of those sites had undergone degradation. To estimate foliar biomass in conifer forests, Moghaddam et al. (2002) found that L-band data from AIRSAR in combination with C-band data and multispectral Landsat data improved results over any of the three sensors alone. Wolter and Townsend (2011) generated a suite of metrics from SAR and optical data and applied partial least squares regression to map forest species and abundance in Minnesota. Antropov et al (2013) developed a semi-empirical approach based on dual pol L-band to estimate stand-level stem volume for boreal sites in Finland. Recently, Cartus et al (2012) applied a water-cloud approach to L-band PALSAR for the northeast USA with success for AGB and noted moisture as an issue. During Phase II we will operationalize a flow building off all these research efforts.Other efforts (e.g. Treuhaft and Siqueira 2000, Papathanassiou and Cloude 2001) have used the full interferometric signature for estimating vegetation height. One advantage of these techniques, and others like it, is that they do not require an external measure of the bare ground topography. The magnitude, also known as the coherence, is related to the volume scattering within a resolution element and therefore can be used for determining height using simple volume-scattering models (Treuhaft and Siqueira 2004), or models that take in to account the motion of the scatterers between observing passes of the interferometric instrument (Lavalle et al. 2012). As part of Phase II we will build an automated routine to scale repeat-pass InSAR maps of Forest Stand Height.The results from all these studies indicate that 1.) SAR data significantly improve estimates of structure over those using optical data alone, and 2.) expert and intensive approaches are required. While the use of L-band data itself is not new, none of the metrics we propose are currently systematic products available to USDA/USFS or other end users. This Phase II will allow us to develop the science and prototype to operationalize products to improve mapping of above ground biomass, forest stand height, canopy cover, and automated degradation and disturbance detection with L-band PALSAR and C-band Sentinel.

Progress 09/01/16 to 05/31/19

Outputs
Target Audience: Nothing Reported Changes/Problems:Our contractor (ArgenTech) had technical trouble with their main UAS; Vapor Pulse, which reduced our field sampling scheme. To augment we integrated additional plots from colleagues at USFS and UNH and shifted our sampling period. What opportunities for training and professional development has the project provided? Open access software released ongithubfor Forest Stand Heighthttps://github.com/leiyangleon/FSHThis software can be used for training. How have the results been disseminated to communities of interest?We had discussions with White Mountain National Forest and University of New Hampshire. We presented our work at the NESAF Annual Winter Meeting New England Forest Stories: the People, the Management, the Technical Knowledge, March, 2018. Additional outreach efforts: Project development with Winrock on MRV tools Member of science team on NASA's NISAR mission Project development with Conservation International on radar MRV tools Promotion at NASA and USDA science team meetings and programs such as SERVIR, GOFC-GOLD, SilvaCarbon, LEAF, IWMI Promotion with TNC on forest monitoring tools Discussions with Prime on Remote Sensing and Geospatial Technology Support Services" contract (Solicitation Number: 1284JC18R0001) under USFS solicitation What do you plan to do during the next reporting period to accomplish the goals? Nothing Reported

Impacts
What was accomplished under these goals? Open access software released on github for Forest Stand Height https://github.com/leiyangleon/FSH Tested ALOS-2 currently for updated baselines and working with NASA JPL and ISCE Executed two main field campaigns focused on assessing Gypsy Moth infestation and another for fused algorithm development with measurements of Above Ground Biomass (AGB), species, TLS, crown canopy cover, DBH, Lorey's Height, and related metrics. Worked with WMNF and USFS community to assess management plans, harvest activities, and disturbances (infestation, ice damage, fire) across WMNF and regional hot spots; 3b. conduct tech transfer to transition tools to WMNF as part of Public Private Partnership. A multitemporal analysis routine (MT) was developed. The MT routine is modular, allowing time series operations to be chained together in potentially complex arrangements, starting with one or more data sources, followed by processing steps, each step having one or more inputs, and an output that can take several different structural forms. MT includes modules for smoothing, recompositing, statistical summaries, trend calculation (e.g., decreasing biomass), seasonal phenological metrics calculation (e.g., Start Of Season), anomaly detection, gap filling, correlation analysis between two time series, and fusion of two time series based on correlation. Taken together, these infrastructure developments have enabled and improved our ability to create fully automated agricultural metrics derived from SAR. During this period we continued to utilize and extend our Geospatial Information Processing System (GIPS) data infrastructure for retrieval, management, and analysis of data from Sentinel-1 and PALAR-2 ScanSAR. GIPS is a modular open source framework that is based on "driver" technology, which enables extension to include data from essentially any source, abstracting away the details of data formats, projections, source, retrieval mechanisms, structure, or units. GIPS drivers include the definition of processing blocks, so we can create custom products on the fly. We've built around the ESA S1Toolbox v5.3 to generate despeckled and terrain corrected gamma and / or sigma naught stacks along with ancillary information such as layover shadow masks and incidence angle maps depending on custom options via XML. This allows for automated data feeds into DSTs or further ingestion into work flows. Annotated data sets hold metadata on the main characteristics including acquisition, image properties, polarization, Doppler information, swath merging, calibration, and geographic location. Our driver is designed to operate with any combination of specifications (e.g., Single Look Complex (SLC), Ground Range Detected (GRD), etc...) depending on application objectives. Worked with ISCE and WINSAR on SW development but determined ISCE is not optimal for commercial applications Deployed GIPS on cloud computing environment

Publications

  • Type: Journal Articles Status: Published Year Published: 2018 Citation: Huang et al., 2018, Assessment of Forest above Ground Biomass Estimation Using Multi-Temporal C-band Sentinel-1and Polarimetric L-band PALSAR-2 Data, Remote Sens.2018,10, 1424; doi:10.3390/rs10091424
  • Type: Conference Papers and Presentations Status: Published Year Published: 2017 Citation: Y. Lei, P. Siqueira, N. Torbick, D. Chowdhury, W. Salas and R. Treuhaft, "Large-scale product of forest height using a new approach from spacborne repeat-pass sar interferometry and lidar," 2017 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Fort Worth, TX, 2017, pp. 2875-2878. doi: 10.1109/IGARSS.2017.8127598;


Progress 09/01/17 to 08/31/18

Outputs
Target Audience:Worked with WMNF and USFS community to assess management plans, harvest activities, and disturbances (infestation, ice damage, fire) across WMNF and regional hot spots; conducted tech transfer to transition tools to WMNF as part of Public Private Partnership. Changes/Problems: Nothing Reported What opportunities for training and professional development has the project provided? We released an open access software released on github for Forest Stand Heigh Working with NASA JPL and ISCE Staff training at Small Business in radar procesing and forest carbon modeling How have the results been disseminated to communities of interest? Presentation at the NESAF Annual Winter Meeting New England Forest Stories: the People, the Management, the Technical Knowledge, March, 2018 Open access software released on github for Forest Stand Height and SAR drivers, ie. https://github.com/leiyangleon/FSH and https://github.com/Applied-GeoSolutions more coming soon Discussions with end users in WMNF and UNH Discussion of project interest and development with Forwarn Project development with Winrock on MRV tools Member of science team on NASA's NISAR mission Project development with Conservation International on radar MRV tools Promotion at NASA and USDA science team meetings and programs such as SERVIR, GOFC-GOLD, SilvaCarbon, LEAF, IWMI Promotion with TNC on forest monitoring tools Discussions with Prime on Remote Sensing and Geospatial Technology Support Services" contract (Solicitation Number: 1284JC18R0001) under USFS solicitation Publication in IEEE, Y. Lei, P. Siqueira, N. Torbick, D. Chowdhury, W. Salas and R. Treuhaft, "Large-scale product of forest height using a new approach from spacborne repeat-pass sar interferometry and lidar," 2017 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Fort Worth, TX, 2017, pp. 2875-2878. doi: 10.1109/IGARSS.2017.8127598; http://ieeexplore.ieee.org/document/8127598/?anchor=authors What do you plan to do during the next reporting period to accomplish the goals? Execute commercialization plan.

Impacts
What was accomplished under these goals? Mapped Above Ground Biomass (AGB) and Crown Canopy Cover (CCC) across White Mountain National Forest and entire northeast USA for ~2008 and ~2016. Example structural mapping of Kingman Farm (NH) using UAS lidar and time series Sentinel-1. "C1" is the control site vs the "T1" is the treatment site. Treatment site contains a more shorter scrubby/shrubby vegetation (regrowth or understory) as compared with the control site. Evaluated ABG modeling. Theassessment scatterplot illustrated random forest model of AGB with strong outcome (R2 0.8), middle highlights variable importance within model, and right shows AGB for development site. Coefficient of variation for cell intensity above 2m had the strong importance in model runs. Initial mapping of gypsy moth disturbance mapping. Results shown below. Operationalized Interferometric Synthetic Aperture Radar (InSAR) Forest Stand Height (FSH) maps across northeast USA using PALSAR-2 (2015-2017 data); 4a. work with USFS, NASA JPL, and EPA on using FSH to improve carbon models

Publications


    Progress 09/01/16 to 08/31/17

    Outputs
    Target Audience:United States Department of Agriculture (USDA) United States Forest Service (USFS)- Forest lands comprise 304million hectares of the total land area in the US with 44% controlled by Federal, State, and local governments. Characterizing forest structure, quantifying forest disturbance, and monitoring forest utilization on these lands is a primary and continuous need by multiple end user communities, including the USDA USFS (see letters). Both natural and anthropogenic stressors are influencing forest disturbance processes and patterns. Climate change, extreme weather, renewable energy demands, timber, land conversion, and invasive species are several key actors influencing forest function and services. These can range in frequency, intensity, and magnitude. Timely and accurate assessment of disturbance events and condition are critical to respond effectively and appropriately at local, regional, and national scales. In addition, assessing the severity and understanding the scales of disturbance are critical for broad-scale management, policy, and economic decision-making. While optical-based mapping techniques have made significant contributions toward assessment, they are limited by cloud coverage and only provided limited information on stand structure. The wide adaptation of operational SAR derived metrics into USFS programs and disturbance protocol has remained elusive due to the technical barriers. In this SBIR we intend to bridge this gap and provide enhanced products and services (metrics, maps, techniques, and expertise) that are easy to apply and integrate into decision support tools. Monitoring, Reporting, and Verification (MRV) community- The rapid growth of carbon science and characterization of ecosystem services has created a huge market for tools, products, and services. International monitoring frameworks, such as Reducing Emissions from Deforestation and Degradation (REDD), have created the need for robust Monitoring, Reporting, and Verification (MRV) platforms that can serve programs and initiatives aimed at carbon cycle science, conservation, and conventions. Maps of deforestation and degradation at multiple scales have taken on increased importance with the recognition of the value of forests and the rates of forest loss. International programs, like REDD, are being constructed to help preserve the world's forests, improve livelihoods, and promote sustainable practices. Under REDD, emissions and removals of carbon due to changes in forests are estimated from two components: activity data (e.g. deforestation, degradation) and forest carbon stock data (e.g. carbon in biomass, soil c). As systems are built in which some people pay others to slow the rate of carbon emissions from forest loss, the need for accurate, spatially explicit, and timely information on forests and forest change becomes critical. Satellite observations have been used to map forests since the 1970s. While optical reflectance data are the most abundant and common type of observations available, radar and Lidar data can be used to effectively map forest characteristics and changes. Optical reflectance observations passively measure the sun's radiation off of the land surface. Trees, forested areas, and changes in these areas are identified based on the amount and change in reflectance. SAR and Lidar are active systems that send out energy to the Earth's surface and measure the subsequent energy return. These active systems often provide more information regarding forest structure and, therefore, are useful for biomass and carbon content. Additionally, SAR systems are much less sensitive to cloud and haze cover than are optical systems, making them useful for forest monitoring in cloudy areas. However, many MRV/REDD platforms do not have the capacity to integrate SAR forest metrics. In this SBIR, we have strategically aligned the goals to serve this community (see letters). Our small business (AGS) has several high profile example MRV or "Decision Support Tools" (DST) currently operating that blend satellite maps, geospatial information, graphing and tabular data, mobile Apps, and end user options within web-GIS. Changes/Problems:Lidar contractor had technical problems. Workign with UNH partners on work around and leveraging GRANIT lidar as solution. What opportunities for training and professional development has the project provided?Conducted workshops and training at NASA Science team meetings (LCLUC) and at academic partbner institutions Serving on PhD commitee at partner institution (UNH) How have the results been disseminated to communities of interest?In part, working on final products now. What do you plan to do during the next reporting period to accomplish the goals?Extensive marketign and promotional campaign with strategic partners

    Impacts
    What was accomplished under these goals? Working with regional USFS in White Mountian National Forest & New England Society of American Foresters http://www.nesaf.org/society-american-foresters-annual-winter-meeting.asp

    Publications