Assessing L-Band GNSS-Reflectometry and Imaging Radar for Detecting Sub-Canopy Inundation Dynamics in a Tropical Wetlands Complex

Despite the growing number of remote-sensing products from satellite sensors, mapping of the combined spatial distribution and temporal variability of inundation in tropical wetlands remains challenging. An emerging innovative approach is offered by Global Navigation Satellite System reflectometry (GNSS-R), a concept that takes advantage of GNSS-transmitting satellites and independent radar receivers to provide bistatic radar observations of Earth’s surface with large-scale coverage. The objective of this paper is to assess the capability of spaceborne GNSS reflections to characterize surface inundation dynamics in a complex wetlands environment in the Peruvian Amazon with respect to current state-of-the-art methods. This study examines contemporaneous ALOS2 PALSAR-2 L-band imaging radar, CYGNSS GNSS reflections, and ground measurements to assess associated advantages and challenges to mapping inundation dynamics, particularly in regions under dense tropical forest canopies. Three derivatives of CYGNSS Delay-Doppler maps (1) peak signal-to-noise ratio (SNR), (2) leading edge slope, and (3) trailing edge slope, demonstrated statistically significant logarithmic relationships with estimated flooded area percentages determined from SAR, with SNR exhibiting the strongest association. Aggregated Delay-Doppler maps SNR time series data examined for inundated regions undetected by imaging radar suggests GNSS-R exhibits a potentially greater sensitivity to inundation state beneath dense forest canopies relative to SAR. Results demonstrate the capability for mapping extent and dynamic wetlands ecosystems in complex tropical landscapes, alone or in combination with other remote-sensing techniques such as those based on imaging radar, contributing to enhanced mapping of these regions. However, several aspects of GNSS-R observations such as noise level, spatial resolution, and signal coherence need to be further examined.


Introduction
Freshwater ecosystems, such as streams, lakes, floodplains, marshes, and swamps, are estimated to cover between 14-29% of the Amazon basin area [1].The seasonal flood pulse and runoff from heavy seasonal rains cause extensive and prolonged flooding of many vegetated regions of the Amazon.These dynamically inundated areas, in addition to forming biologically diverse habitats, also regulate biogeochemical processes such as the generation of atmospheric methane and the outgassing of carbon dioxide [2,3].These flooding regimes can be impacted by future changes in rainfall, evapotranspiration and land use.The ability to accurately monitor the current state and changes in inundation extent would enable further examination of potential climatological and anthropogenic tipping points in these wetlands regions.
Satellite remote sensing is the only practical approach that can provide insight into the spatial and temporal dynamics of wetlands on a large scale on a continuous basis.Although they can offer high spatial resolution and sensitivity to photochemical properties of vegetation, optical sensors are severely limited in capturing inundation dynamics in tropical wetlands due to frequent cloud cover and difficulties detecting sub-canopy inundation [4,5].Microwave sensors, on the other hand, are not significantly affected by clouds or changing solar illumination and are able to observe processes below vegetation canopy.However, tradeoffs exist in the use of microwave remote sensing with respect to characterizing dynamic tropical wetlands.Low spatial resolution (e.g., 20-40 km) radiometers and scatterometers have the advantage of frequent, often daily, coverage but limit detection to large wetlands or to regions where the cumulative area of small wetlands is a substantial portion of the instrument field of view.High spatial resolution (<100 m) synthetic aperture radar (SAR) provides more detailed spatial information but suffers from lower temporal coverage (6-12 days, but often more infrequent).
The nature of backscatter responses to surface water depends on the presence and character of vegetation cover.The complex assemblies of vegetation and dynamic surface hydrologic conditions often found in wetlands contribute to a diverse variety of backscatter signatures observed by SAR in these regions.Specular surfaces such as open water return low backscatter, while flooded vegetation can return a very strong backscatter response because of so-called double-bounce scattering.Depending on the density and structure of overlying vegetation, the radar energy can be substantially attenuated by the canopy.Apart from vegetation conditions, this attenuation depends largely on the sensor frequency (wavelength).With higher frequency, there is less radar energy penetration through the canopy and more contribution to backscatter from within-canopy volume scattering.L-band is the most common SAR frequency of choice used for remote sensing of wetland extent and dynamics over large regions because it provides superior canopy penetration and water surface discrimination relative to higher frequency microwave sensors [6,7].
Preliminary results based on data collections from the NASA's aircraft-mounted UAVSAR imaging radar over the Napo River in Ecuador and findings from Amazon-scale inundation dynamics mapped using PALSAR show notable limitations of L-band SAR-based inundation extent maps, particularly with systemic underestimation of inundated areas in dense forest [6,8,9].Yuan et al. [10] found that increases in double-bounce backscatter with rising water levels in flooded vegetation was significantly lower in areas characterized by high amounts of woody biomass.Although SAR offers the current best sensitivity to inundation dynamics in wetlands, there remains a need for focused efforts to assess uncertainty in and improve the quality of inundation extent mappings provided by microwave remote sensing.
Global Navigation Satellite System reflectometry (GNSS-R) is a bistatic radar concept that takes advantage of L-band signals transmitted by navigation constellations such as the Global Positioning System (GPS).The GPS signals reflect from the Earth's surface with opportunity for global-scale bistatic radar coverage by appropriately positioned Earth-orbiting receivers.GNSS-R offers a new, complementary approach to monitoring inundation with Earth-orbiting radar.Although GNSS-R was first developed for measuring sea surface parameters [11,12], more recent studies have focused on the monitoring of terrestrial surface characteristics including soil moisture [13], vegetation phenology [14,15], mapping of flooded areas [16,17], and the detection of snow and ice [18,19].Experimental campaigns from fixed [20,21], airborne [22][23][24], and Low-Earth Orbit platforms have confirmed the sensitivity of GNSS receivers to soil and biomass characteristics [25].Parallel to these experimental developments, progress has been made with theoretical bistatic models that have improved the understanding of the forward-scattering effects taking place in GNSS-R [26][27][28].This technology is gaining momentum as demonstrated by the launch of the UK TechDemoSat-1 (TDS-1) in July 2014 and the launch of the NASA Cyclone GNSS (CYGNSS) mission in December 2016, both of which carry GNSS-R receivers.
GNSS-R scattering and associated surface characteristics are determined from Delay-Doppler maps (DDM) of the bistatically-scattered signal.A DDM is produced by cross-correlating the recorded signal from a downward-looking antenna with a replica of the pseudorandom noise of the GNSS satellite over a range of time lags and carrier frequency offsets [29].Various aspects of the DDM can be used for analysis (e.g., peak power, leading or trailing slope of the waveform, volume).Studies from a 2005 aircraft field campaign and the more recent TDS-1 satellite observations over known wetland areas have demonstrated the capability of GNSS-R DDM peak power measurements to detect standing water beneath low vegetation [30].The forward scatter from GNSS reflected that signals strengthen in the presence of standing surface water in spite of vegetation attenuation.There have yet to be studies that quantify the amount of vegetation that GNSS-R can penetrate or examine the relative sensitivities of GNSS-R and L-band SAR backscatter when observing inundation below vegetation canopies.
The goal of this research is to assess the capability of spaceborne GNSS reflections to characterize surface inundation dynamics in a tropical wetlands complex with respect to current state-of-the-art SAR-based methods.In an integrative analysis of contemporaneous L-band SAR, GNSS-R, and ground measurements, this study examines the associated advantages and obstacles of mapping inundation dynamics under the canopy with these datasets.

Study Area
This study focuses on Peru's Ucayali-Marañon region of the Western Amazon (Figure 1).This region is home to the Pacaya Samiria National Reserve, a protected area bordered by the Marañon and Ucayali rivers.Spanning over 20,000 km 2 , this reserve encompasses the largest tropical flooded forest in the Peruvian Amazon and is the most expansive protected wetlands area in entire Amazon Basin.Topographic and climatic conditions, along with water availability, support the Pacaya Samiria's high ecological richness, biodiversity, and prevalent wetlands ecosystems.The seasonally dynamic flooding of this area with a diverse mosaic of vegetation and open water makes this an ideal study site for evaluating remote sensing assessments of inundation dynamics across a range of tropical vegetation conditions.
The dominant sources of water to these ecosystems are rainfall and upland run-off [1].There is a pronounced rainy season which lasts approximately from December to April, with the seasonal flood pulse raising river water levels in the order of 10 m (Figure 2).As a result, a vast amount of the region experiences enhanced inundation for months every year.The duration of this flooding can vary considerably based on micro-topography and can drive the local vegetation composition.

ALOS2 PALSAR-2 ScanSAR Multi-Tmporal Mosaics
The Advanced Land Observing Satellite 2 (ALOS2) Phased Array L-band Synthetic Aperture Radar 2 (PALSAR-2) is an orbiting L-band SAR that was launched by the Japanese Aerospace Exploration Agency (JAXA) in May 2014.This study uses images acquired while the sensor was operating in ScanSAR mode, which is characterized by a wide swath up to 350 km.This mode allows frequent data acquisition over large regions such as the Amazon [34].Through the ALOS Kyoto and Carbon (K&C) Initiative, processed dual polarized (HH and HV) mosaics have been assembled from ScanSAR data that have undergone radiometric calibration and terrain correction.These data have been formatted in 16-bit amplitude mosaics of topographic normalized backscatter coefficient (γ 0 ) posted at 50 m resolution in 1° × 1° tiles for each 14-day ScanSAR data acquisition cycle.Cycles between May 2014 through February 2018 were available for this analysis, 29 of which were found to either fully or partially observe the study area.Acquisition dates are summarized in Table A1 (see

ALOS2 PALSAR-2 ScanSAR Multi-Tmporal Mosaics
The Advanced Land Observing Satellite 2 (ALOS2) Phased Array L-band Synthetic Aperture Radar 2 (PALSAR-2) is an orbiting L-band SAR that was launched by the Japanese Aerospace Exploration Agency (JAXA) in May 2014.This study uses images acquired while the sensor was operating in ScanSAR mode, which is characterized by a wide swath up to 350 km.This mode allows frequent data acquisition over large regions such as the Amazon [34].Through the ALOS Kyoto and Carbon (K&C) Initiative, processed dual polarized (HH and HV) mosaics have been assembled from ScanSAR data that have undergone radiometric calibration and terrain correction.These data have been formatted in 16-bit amplitude mosaics of topographic normalized backscatter coefficient (γ 0 ) posted at 50 m resolution in 1° × 1° tiles for each 14-day ScanSAR data acquisition cycle.Cycles between May 2014 through February 2018 were available for this analysis, 29 of which were found to either fully or partially observe the study area.Acquisition dates are summarized in

ALOS2 PALSAR-2 ScanSAR Multi-Tmporal Mosaics
The Advanced Land Observing Satellite 2 (ALOS2) Phased Array L-band Synthetic Aperture Radar 2 (PALSAR-2) is an orbiting L-band SAR that was launched by the Japanese Aerospace Exploration Agency (JAXA) in May 2014.This study uses images acquired while the sensor was operating in ScanSAR mode, which is characterized by a wide swath up to 350 km.This mode allows frequent data acquisition over large regions such as the Amazon [34].Through the ALOS Kyoto and Carbon (K&C) Initiative, processed dual polarized (HH and HV) mosaics have been assembled from ScanSAR data that have undergone radiometric calibration and terrain correction.These data have been formatted in 16-bit amplitude mosaics of topographic normalized backscatter coefficient (γ 0 ) posted at 50 m resolution in    tiles for each 14-day ScanSAR data acquisition cycle.Cycles between May 2014 through February 2018 were available for this analysis, 29 of which were found to either fully or partially observe the study area.Acquisition dates are summarized in Table A1 (see Appendix A).Statistical analysis of backscatter was performed using amplitude values (digital numbers, DN).Final results were converted from DN to normalized backscatter coefficients expressed in dB using the following equation: where CF is a calibration factor (−83) used for standard PALSAR products [35].

Image Classification
SAR images were classified for open water and inundated vegetation to support evaluation of contemporaneous GNSS reflected signals.A two-step, object-oriented classification scheme was developed using temporal backscatter patterns (Figure 3), building on methods presented by Reference [36].Thresholds applied in the decision tree classifier were derived from analysis of a number of small regions of interest (ROIs).Twelve ROIs were selected pertaining to different homogeneous land cover classes based on a priori knowledge of the reserve, making a distinction between forest and non-forest (NF) vegetation (Figure 4 and Table 1).Multi-temporal statistics of the HH-and HV-polarized backscatter time series were computed Appendix A).Statistical analysis of backscatter was performed using amplitude values (digital numbers, DN).Final results were converted from DN to normalized backscatter coefficients expressed in dB using the following equation: where CF is a calibration factor (−83) used for standard PALSAR products [35].

Image Classification
SAR images were classified for open water and inundated vegetation to support evaluation of contemporaneous GNSS reflected signals.A two-step, object-oriented classification scheme was developed using temporal backscatter patterns (Figure 3), building on methods presented by Reference [36].Thresholds applied in the decision tree classifier were derived from analysis of a number of small regions of interest (ROIs).Twelve ROIs were selected pertaining to different homogeneous land cover classes based on a priori knowledge of the reserve, making a distinction between forest and non-forest (NF) vegetation (Figure 4 and Table 1).Multi-temporal statistics of the HH-and HV-polarized backscatter time series were computed, including: minimum ( γ , γ ), maximum ( γ , γ ), mean ( γ , γ ), and standard deviation (γ , γ ).Additional dry season information was collected including: (1) mean dry season HV-polarized backscatter (γ ) calculated from scenes observed in July-September (Cycles 27, 30, 53, 56, 79, 82), and (2) median dry season normalized difference vegetation index (NDVI-DS) from Landsat 8 derived from cloud-free observations made in July-September of 2015-2017.A general assumption was made that flooding occurs only in regions not significantly elevated above the river network, as supported by past literature e.g., [37,38].A potential flood area mask was created based on the Height Above Nearest Drainage (HAND) index [39,40], which is effectively the elevation difference between a DEM pixel and the nearest pixel that is part of the drainage network.The HAND index was calculated from flow direction and a void-filled DEM from the HydroSHEDS database [41], that is derived from the Shuttle Radar Topography Mission (SRTM) [42] at 3 arc-second resolution.Similar to what was done by Reference [38], the stream network was defined by deriving flow accumulation for each pixel and identifying drainage pixels as those with 1000 or greater contributing pixels.The resulting stream network was further reduced by retaining segments with a Strahler stream order of three or greater.The height differences were calculated from the void-filled DEM provided by HydroSHEDS.A threshold of 20 m above the drainage network was used to mask out regions not prone to flooding.This approach supports a more appropriate estimate of the floodable area in regions of low topographic relief compared to estimates based on elevation alone, as regions of slowly increasing elevation may be accounted for as part of the inundatable wetlands complex.In Step 1A (Figure 3), the maximum flooded extent observed in the entire ScanSAR time series (November 2014-February 2018) was delineated, distinguishing between open water, inundated vegetation and non-inundated areas.In Step 1B, masked regions and those found to never flood were merged into a single class named "Non-Inundated."The Inundated Vegetation class is split into Inundated Forest and Inundated Non-

Classification Evaluation
Classification accuracy was evaluated for both Step 1 and 2 using independent approaches.The maximum flooded extent (Step 1) was evaluated against a static land cover map of the Pacaya Samiria National Reserve produced by the Instituto de Investigaciones de la Amazonia Peruana (IIAP) based on aerial surveys, ground observations, and Landsat imagery [43].Pixels classified as intermittent open water in Step 1 were excluded from this assessment since only permanent water extent is mapped by IIAP.Furthermore, pixels classified as "Intervened Forest" and "Urban Population Centers" by IIAP were also excluded from the assessment as it was unclear if these were considered floodable by IIAP.This afforded 5,217,938 of testing pixels at 50m grid posting to assess against the maximum extent map (~16% of the entire study area extent).The IIAP land cover classes were generalized to match the Step 1 classes, summarized in Table A3 (Appendix B).
Individual cycle classifications generated in Step 2 were more difficult to robustly assess due to lack of a time series of ground reference data.River level data from four monitoring stations in the reserve (locations shown in Figure 1b) were compared against the time series of flooded area percentage in a 5 km radius of each station.

CYGNSS
The Cyclone GNSS (CYGNSS) mission was launched in December 2016.CYGNSS consists of a constellation of eight satellites in a low inclination orbit, each carrying a specialized GPS navigation receiver integrated with a reflections processor, with each receiver capable of receiving four GNSS reflection signals simultaneously [44].With each satellite producing four DDMs per second, the CYGNSS constellation receives 32 concurrent measurements.Although the primary focus of the CYGNSS mission is to retrieve ocean surface winds, it also records measurements over land surfaces.Daily observations from 18 March 2017 through 28 February 2018 of CYGNSS Level 1 (L1) version 2.1 were provided by the CYGNSS Science Team, offering nearly a year of overlap with the PALSAR-2 ScanSAR time series (Cycles 71, 74, 77, 79, 82, 85, 88, 91, 93).

Data Processing
This analysis focuses on three aspects of the DDM: The peak signal to noise ratio (SNR), leading edge slope (LES) and trailing edge slope (TES).SNR, as defined in Equation (2), is provided as an observable corrected for antenna gain and range in the L1 dataset, while the LES and TES were derived from the two-dimensional (Doppler × delay) bistatic radar cross section DDM (Figure 5).There are 27 quality flags associated with each DDM observation, some of which are critical for ocean applications but not necessary for terrestrial studies.A summary of the quality flags that were used in this analysis is presented in Table A2 (Appendix A).As observed with traditional imaging radar, the incidence angle (θ inc ) affects the signal path length and can noticeably impact the signal particularly at both very steep and very shallow angles.For these reasons, observations with incidence angles of <15 • and >60 • were excluded from the analysis.
where P is maximum power of the DDM, and N is the background noise level.The baseline Level 2 CYGNSS wind speed product has a nominal spatial resolution of 25 km × 25 km (centered on the specular reflection point), but the native resolution is much greater.CYGNSS averages multiple DDM bins together to yield the 25 km × 25 km footprint.However, the L1 data products can be exploited to optimize resolution over wetlands.Previous analyses have shown that over land surface, particularly wetlands, there is a dominant coherent component to the reflected signal.Thus, the active scattering area is defined by the first Fresnel zone, ~650 m at smaller incidence angles (<40 • ) and reaches up to ~1 km at larger angles (>50 • ) [30].The effective cross-track resolution is assumed as this first Fresnel zone size, but the along-track resolution is elongated since CYGNSS measurements are averaged incoherently over 1 s and the distance traveled by the specular point during this time is ~6 km.CYGNSS DDM SNR, LES, and TES values were mapped to a fixed 50 m grid identical to the one at which the ALOS2 PALSAR-2 ScanSAR mosaics are posted, adopting a method first presented by Reference [45].The size of the semi-major and semi-minor axes of the ellipses corresponding to individual DDM footprints were defined using Equations ( 3) and ( 4), respectively: where ε is the elevation angle of the reflected signal with respect to the local tangent plane at the specular point, or 90° − θinc; Rrx is the distance between the receiver antenna and the specular point; Rtx is the distance between the transmitter antenna to the specular point; λ is the wavelength of the signal (0.19 m).Derivation of these equations is described in more detail by Reference [46].The footprint ellipse is rotated by the angle of the track direction around the specular point, and all grid cells found to be in the ellipse are assigned the corresponding DDM values.DDM footprints were examined using the data posted at this 50 m grid (Section 3.4.1).However, a pixel-wise time series analysis cannot be done properly using these data since in this case individual DDM SNR, LES, and TES values are not representative of any single 50 m pixel, but rather the entire footprint.Footprints were aggregated to a 0.01° (~1 km) grid, where grid cells are assigned DDM values when there is 50% or greater coverage by a DDM footprint.Multi-temporal observations were assessed in Sections 3.4.2and 3.4.3using DDM values at this grid posting.

Ground Reference Data
Field data were collected in July-August 2017 during low-water conditions in the Pacaya Samiria National Reserve to assess the impact of biomass on the ability to detect inundation with SAR and GNSS-R observations.During the 18-day excursion, ground measurements were collected on vegetation structure and past inundation depth history at four primary floodplain sites along Canal de Puinahua (locations shown in Figure 1b; photographs shown in Figure 6).CYGNSS DDM SNR, LES, and TES values were mapped to a fixed 50 m grid identical to the one at which the ALOS2 PALSAR-2 ScanSAR mosaics are posted, adopting a method first presented by Reference [45].The size of the semi-major and semi-minor axes of the ellipses corresponding to individual DDM footprints were defined using Equations ( 3) and ( 4), respectively: where ε is the elevation angle of the reflected signal with respect to the local tangent plane at the specular point, or 90 • − θ inc ; R rx is the distance between the receiver antenna and the specular point; R tx is the distance between the transmitter antenna to the specular point; λ is the wavelength of the signal (0.19 m).Derivation of these equations is described in more detail by Reference [46].The footprint ellipse is rotated by the angle of the track direction around the specular point, and all grid cells found to be in the ellipse are assigned the corresponding DDM values.DDM footprints were examined using the data posted at this 50 m grid (Section 3.4.1).However, a pixel-wise time series analysis cannot be done properly using these data since in this case individual DDM SNR, LES, and TES values are not representative of any single 50 m pixel, but rather the entire footprint.Footprints were aggregated to a 0.01 • (~1 km) grid, where grid cells are assigned DDM values when there is 50% or greater coverage by a DDM footprint.Multi-temporal observations were assessed in Sections 3.4.2and 3.4.3using DDM values at this grid posting.

Ground Reference Data
Field data were collected in July-August 2017 during low-water conditions in the Pacaya Samiria National Reserve to assess the impact of biomass on the ability to detect inundation with SAR and GNSS-R observations.During the 18-day excursion, ground measurements were collected on vegetation structure and past inundation depth history at four primary floodplain sites along Canal de Puinahua (locations shown in Figure 1b; photographs shown in Figure 6).1b) with sparse tree cover dominated by Pseudobombax munguba; (c) Dense canopy cover found in Site #2 (in Figure 1b); (d) inundated aquatic shrub species (Montrichardia arborescens) found along margins of lake at Site #3(in Figure 1b); (e) denser tree cover and understory found in Site #3 (in Figure 1b); (f) inundated swamp found at Site #4 (in Figure 1b).Examples of high flood marks on trees (shown red rectangles) indicate the level of previous inundation.

Circular Plots and Biomass Estimation
Biometric measurements were made at seven circular plots (15 m radius each) hosting varying levels of biomass density.Stand leaf area index (LAI) (m 2 m −2 ) was estimated using the LAI-2200C plant canopy analyzer (Li-Cor, Inc., Lincoln, Nebraska, NE, USA) following the procedure detailed in Appendix C. Canopy closure (%) was estimated using a densiometer.Four readings were taken in each of the cardinal directions in each quadrant and then averaged.In addition, evidence of past flooding stage was recorded to the nearest 10 cm.Individual tree measurements were recorded for each living tree ≥5 cm in diameter found in each plot.Diameter at breast height (D) was measured with a diameter tape at~1.3 m and recorded to the nearest 0.1 cm.Height to the base of the live crown (HC) and total height (HT) were estimated using a distance meter; and crown depth (CD) was calculated as the difference between HT and HC.All trees were identified as their species or genus when possible and assigned a wood density value (ρ, oven-dry mass over green volume, g cm −3 ) from the Global Wood Density Database [47,48].Unidentified trees were assigned the average ρ for their containing plot.Aboveground biomass (AGB) was estimated for individual trees using the Equation ( 5), an allometric equation for moist pantropical forests developed by Reference [49].AGB density (Mg ha −1 ) was estimated for each plot by summing the individual tree AGBs per plot. = 0.0673 ×    . (5)

Line Transects
In addition to the detailed measurements recorded at circular plots, a limited set of biometry measurements were collected along six line transects (~0.5-1 km in length) in 25-50 m intervals, depending on canopy conditions.Measurements taken at each point included: LAI, canopy closure, and any past inundation depth-all collected in a similar fashion as those recorded in the circular plots (detailed in Section 2.4.1).One longer transect (~3 km in length) was made via canoe through an inundated swamp.Although no biometry measurements were collected along this transect, notes on vegetation character, inundation, and geolocated photographs were recorded approximately every 250 m.  1b) with sparse tree cover dominated by Pseudobombax munguba; (c) Dense canopy cover found in Site #2 (in Figure 1b); (d) inundated aquatic shrub species (Montrichardia arborescens) found along margins of lake at Site #3(in Figure 1b); (e) denser tree cover and understory found in Site #3 (in Figure 1b); (f) inundated swamp found at Site #4 (in Figure 1b).Examples of high flood marks on trees (shown red rectangles) indicate the level of previous inundation.

Circular Plots and Biomass Estimation
Biometric measurements were made at seven circular plots (15 m radius each) hosting varying levels of biomass density.Stand leaf area index (LAI) (m 2 m −2 ) was estimated using the LAI-2200C plant canopy analyzer (Li-Cor, Inc., Lincoln, Nebraska, NE, USA) following the procedure detailed in Appendix C. Canopy closure (%) was estimated using a densiometer.Four readings were taken in each of the cardinal directions in each quadrant and then averaged.In addition, evidence of past flooding stage was recorded to the nearest 10 cm.Individual tree measurements were recorded for each living tree ≥5 cm in diameter found in each plot.Diameter at breast height (D) was measured with a diameter tape at~1.3 m and recorded to the nearest 0.1 cm.Height to the base of the live crown (HC) and total height (HT) were estimated using a distance meter; and crown depth (CD) was calculated as the difference between HT and HC.All trees were identified as their species or genus when possible and assigned a wood density value (ρ, oven-dry mass over green volume, g cm −3 ) from the Global Wood Density Database [47,48].Unidentified trees were assigned the average ρ for their containing plot.Aboveground biomass (AGB) was estimated for individual trees using the Equation ( 5), an allometric equation for moist pantropical forests developed by Reference [49].AGB density (Mg ha −1 ) was estimated for each plot by summing the individual tree AGBs per plot.

Line Transects
In addition to the detailed measurements recorded at circular plots, a limited set of biometry measurements were collected along six line transects (~0.5-1 km in length) in 25-50 m intervals, depending on canopy conditions.Measurements taken at each point included: LAI, canopy closure, and any past inundation depth-all collected in a similar fashion as those recorded in the circular plots (detailed in Section 2.4.1).One longer transect (~3 km in length) was made via canoe through an inundated swamp.Although no biometry measurements were collected along this transect, notes on vegetation character, inundation, and geolocated photographs were recorded approximately every 250 m.

Sensitivity Analysis
The SAR-based classifications were utilized in two general ways: (1) To assist in interpreting CYGNSS DDM observations, and (2) to evaluate against CYGNSS DDMs with respect to sensitivity to surface inundation.First, the footprints of CYGNSS DDMs were examined with contemporaneous flooded state, estimated using the derived PALSAR-2 classifications.Secondly, the limitations of the SAR-based classification were assessed with ground information.Time series of both the SAR-based classifications and CYGNSS DDMs of known areas were examined.
The Global Forest Watch developed a pantropical 30 m resolution dataset of estimated woody AGB density of the year 2000 [50].This dataset was derived from a combination of ground measurements, Geoscience Laser Altimeter System (GLAS) LiDAR, Landsat 7 ETM+ surface reflectance, elevation, and biophysical variables.Although this is a dated estimate, a large portion of the study area is protected and hence, large amounts of drastic changes in vegetation over the past 18 years are assumed unlikely.Plot level AGB density values were compared against corresponding values in GFW-AGB to assess the suitability of large scale analysis using the GFW-AGB dataset as a proxy for biomass density across the study area.The GFW-AGB data were co-registered to the 50 m grid at which the PALSAR-2 SAR data were posted using nearest-neighbor interpolation and used to assess the impact of biomass on inundation detection.

Multi-Temporal SAR Signatures
Multi-temporal HH-and HV-polarized radar backscatter from PALSAR-2 and dry-season NDVI from Landsat 8 were examined in 12 ROIs (representing six land cover types).Figure 7 shows the histograms of these classes (ROIs of the same class were merged) for eight multi-temporal images.A thresholding approach was developed to delineate flooded extent by land cover type, as demonstrated in previous studies of classifying wetland types with L-band SAR e.g., [6,31,36,51].Open water, which presents a highly reflective surface for which specular radar scatter dominates, exhibits low backscatter values.Thus, the maximum extent of open water regions can be delineated with γ 0 HH−Min and γ 0 HV−Min (Figure 7a,c).Regions of temporally intermittent open water can be differentiated from permanent open water areas with γ 0 HV−Max .(Figure 7d).High biomass vegetation where any inundation occurs can be identified by high γ 0 HH−Max (Figure 7c) due to plant-water double bounce scattering.Low vegetation (e.g., pioneer formations, such as grasses and other herbaceous species) where inundation may occur poses a more complex backscatter dynamic with greater temporal variation (Figure 7e,f) due to factors like soil moisture, vegetation structure, and varying inundation depth.If partially flooded, double bounce interactions may yield high HH-polarized backscatter, but if submerged, these regions will appear dark with low backscatter.Forested areas observe both greater γ 0 HV−DS (Figure 7g) and NDVI-DS (Figure 7h) than lower density, non-forest vegetation.

SAR Classification
The classification algorithm and decision tree derived for Step 1 are illustrated in Figure 8a, with the resulting map shown in Figure 8b.This maximum flooded extent classification showed generally high agreement with the IIAP Pacaya Samiria land cover map, depicting inundated regions with an overall accuracy of 91.8% and kappa index value of 0.84 (Table 2).

SAR Classification
The classification algorithm and decision tree derived for Step 1 are illustrated in Figure 8a, with the resulting map shown in Figure 8b.This maximum flooded extent classification showed generally high agreement with the IIAP Pacaya Samiria land cover map, depicting inundated regions with an overall accuracy of 91.8% and kappa index value of 0.84 (Table 2).The highest commission error was observed for the Non-Flooded class (20.2%), which suggests notable underestimation of flooded maximum extent.Open Water, Inundated Non-Forest, and Inundated Forest each had comparable omission error rates (8-11%).Floodable regions of low vegetation areas may potentially have been classified as Open Water at maximum extent if they were completely submerged, but may have been identified as Inundated Non-Forest by the IIAP map.
Remote Sens. 2018, 10, 1431 13 of 29 Some regions classified as Non-Inundated noticeably follow along river banks.Many of these regions exhibit variations in micro-topography with elevated embankments that have shallow or no inundation depths as detected by the SAR observations.Furthermore, it is possible that regions of particularly dense canopy conditions limit the ability for radar to detect inundation.Other potential error sources that are not related to the backscatter retrieval include: (1) Disagreement in definitions of Forest and Non-Forest vegetation, (2) co-registration issues between IIAP map and PALSAR-2 ScanSAR mosaics and (3) changes in river extent and vegetation regimes since the IIAP map was derived in 2002.The classification of individual ScanSAR cycles in Step 2 was performed using the scheme summarized in Figure 9. Daily river level at four field stations (locations shown in Figure 1b) were compared with the corresponding percentage of local classified flood extent (considering a radius of 5 km around each station) (Figure 10).Since the relationship between river level and flood extent is not linear, observations were fit to a logistic curve and demonstrated reasonable goodness of fit (R 2 ranging from 0.92-0.95).

Biometry Measurements and Biomass Estimation
A variety of primary and secondary forest stands were surveyed, providing a range of biomass conditions for inundation detection assessment.In Appendix D, biometry data of circular plot characteristics are summarized in Table A4 and descriptive statistics derived from line transects are summarized in Table A5.Each stand at which circular plots were made showed evidence of past flooding, while 78 out of the 81 transect points surveyed did so.Circular plot-level estimated AGB density ranged from 42.2-278.3Mg ha −1 .Although this only provides seven data points, these values were compared to the GFW-AGB dataset and illustrated in Figure 11.The GFW dataset appears to over-estimate AGB density compared to our biometry collections.The one site where the GFW-AGB density was found to be less than the plot-measured AGB had considerable biomass contributed to from several ficus trees (Ficus trigonata) which exhibited large networks of exposed roots.It is possible that these unique features could lead to an underestimation of standing biomass in the GFW data.The strong, positive correlation (Pearson's correlation coefficient = 0.835, two-tailed p-value = 0.019) is encouraging, and the GFW-AGB dataset was deemed suitable for large-scale analysis.that these unique features could lead to an underestimation of standing biomass in the GFW data.The strong, positive correlation (Pearson's correlation coefficient = 0.835, two-tailed p-value = 0.019) is encouraging, and the GFW-AGB dataset was deemed suitable for large-scale analysis.Transect measurements demonstrated positive relationships among vegetation structure metrics, with observed LAI and canopy closure shown in Figure 12a, and each of these compared to GFW-AGB in Figure 12b,c respectively.Figure 13a shows the long transect made through a swamp with a vegetation gradient.It leads from a small lake (point 0) through tall grasses and other herbaceous species (points 1, 7, 8), mixed herbaceous and sparse trees (points 2-6), mixed trees and herbaceous communities (points 9-13), and forest (points 14-18), ending near the river bank.At the time of observation (30 July 2017) in the middle of low-water season, points 1-15 were inundated while 16-18 were dry but showed evidence of past seasonal flooding to depths >0.5 m.This transect sampled three 0.01° × 0.01° grid cells (labelled A, B and C in Figure 18a) that displayed a general increase in biomass density.A and B are both primarily inundated swamp, with B hosting greater amounts of trees, and C primarily observed as a seasonally inundated forest.The GWF-AGB density and γ (Figure 13b) both generally increase across the transect and appear strongly correlated, with the last three points (16)(17)(18) exhibiting an indication of high biomass.Transect measurements demonstrated positive relationships among vegetation structure metrics, with observed LAI and canopy closure shown in Figure 12a, and each of these compared to GFW-AGB in Figure 12b,c respectively.Figure 13a shows the long transect made through a swamp with a vegetation gradient.It leads from a small lake (point 0) through tall grasses and other herbaceous species (points 1, 7, 8), mixed herbaceous and sparse trees (points 2-6), mixed trees and herbaceous communities (points 9-13), and forest (points 14-18), ending near the river bank.At the time of observation (30 July 2017) in the middle of low-water season, points 1-15 were inundated while 16-18 were dry but showed evidence of past seasonal flooding to depths >0.5 m.This transect sampled three 0.01 • × 0.01 • grid cells (labelled A, B and C in Figure 13a) that displayed a general increase in biomass density.A and B are both primarily inundated swamp, with B hosting greater amounts of trees, and C primarily observed as a seasonally inundated forest.The GWF-AGB density and γ 0 HV−DS (Figure 13b) both generally increase across the transect and appear strongly correlated, with the last three points (16)(17)(18) exhibiting an indication of high biomass.

Footprint Analysis of SAR Classification vs. CYGNSS DDMs
For each CYGNSS footprint observed during each of the 14-day ALOS2 PALSAR-2 ScanSAR cycles, the percentage of flooded area in each footprint was calculated based on the PALSAR-2 dynamics inundation classification and plotted against corresponding CYGNSS DDM observations.Summary plots are shown in Figure 14, with the top row (Figure 14a-d

Footprint Analysis of SAR Classification vs. CYGNSS DDMs
For each CYGNSS footprint observed during each of the 14-day ALOS2 PALSAR-2 ScanSAR cycles, the percentage of flooded area in each footprint was calculated based on the PALSAR-2 dynamics inundation classification and plotted against corresponding CYGNSS DDM observations.Summary plots are shown in Figure 14, with the top row (Figure 14a-d) presenting SNR, center row (Figure 14e-h) presenting LES, and bottom row (Figure 14i-l) presenting TES.The leftmost column displays the relationship observed between the percentage of all flooded area (Open Water, Inundated Non-Forest, and Inundated Forest) of all CYGNSS footprints and the respective DDM values.The center left column illustrates a subset of these footprints, showing only those with negligible amount of classified flooded vegetation (<5%) and focusing on the dependence of DDM measurements on the flooded area percentage of Open Water.The center right column features a different subset, containing footprints with a negligible amount of Open Water (<1%) and illustrates the relationship between the area percentage of flooded vegetation (Inundated Non-Forest and Inundated Forest) and DDM values.Lastly, the rightmost column is an even smaller subset of those footprints in the third column, consisting of those with negligible amounts of Inundated Non-Forest and Open Water (<1%) and thereby isolating the relationship between the area of Inundated Forest and DDM measurements.The CYGNSS footprint data exhibited a high level of noise so were binned in equal intervals (n = 1000) by classified flood percentage in each plot.All DDM observations exhibited a generally logarithmic relationship with flooded area percentages regardless of vegetation conditions, with DDM SNR and LES increasing with greater flooded footprint area and TES decreasing with flooded area.SNR and TES observed greater goodness of fit with a logarithmic curve (R 2 ranging from 0.88 to 0.96) than LES.TES appeared to exhibit greater variability than SNR at higher flooded area percentages.The GNSS-reflected signal approaches saturation at variable percentages of flooding, depending on vegetation cover and the type of observable.The SNR in particular appears to saturate at lower flooded percentages, showing near saturation at flooded area percentages as low as 20%.The maximum value at which the DDM values saturated varied across the footprint subsets, illustrating the sensitivity of the GNSS-R signal to canopy attenuation.The largest ranges of binned DDM measurements are observed in the center left column which focuses on Open Water alone (SNR: 0 to ~11, LES: 0 to ~1500, TES: 0 to −70), and the narrowest range observed in the rightmost column focusing on Inundated Forest percentages alone (SNR: 0 to 4.5, LES: 0 to ~70, TES: 0 to −3).
The influence of canopy attenuation was examined further by binning the DDM observations by overall flooded area percentage (Open Water, Inundated Non-Forest, Inundated Forest) in four distinct biomass density categories, as informed by the GFW-AGB dataset (Figure 15).The best-fit logarithmic curves derived from the lowest AGB density category (<100 Mg ha −1 ) observed the largest ranges in DDM values, with signal range decreasing as biomass density increases.The notable differences in the GNSS-R signatures with respect to flooded area percentage observed by these four biomass categories underscores the need for land cover information to interpret the GNSS-R signal.differences in the GNSS-R signatures with respect to flooded area percentage observed by these four biomass categories underscores the need for land cover information to interpret the GNSS-R signal.

Comparison of SAR Classification and Peak DDM SNR at 1 km Scale
Analysis of multi-temporal DDM observations are focused on SNR since it exhibits the least noise of the three observations examined.DDM SNR values are posted to a 0.01° (~1 km) grid, with the time series maximum displayed over the study area in Figure 16a.The distribution of temporal maximum SNR in grid cells dominantly (>99%) classified as a single maximum flooded extent class with the PALSAR-2 data shown in Figure 16b.The median peak SNR value observed in Open Water grid cells is highest, at 13.97 dB, second highest in Inundated Non-Forest at 12.70 dB, and third highest in Inundated Forest at 8.54 dB, following the inverse relationship observed between peak SNR and amount of canopy attenuation above the flooded surface.The Non-Inundated class observed the lowest maximum SNR, with a median value among all dominantly Non-Inundated grid cells of 2.24 dB.There is a degree of overlap between the distributions of the Non-Inundated and all of the floodable classes, particularly between the top 75% of the Non-Inundated and the bottom 50% of the Inundated Forest, indicating that it may be difficult to differentiate flooded dense forest from non-inundated forest with GNSS-R signals alone.However, it also suggests there are potential regions of seasonally flooded vegetation that failed to be classified as flooded using the PALSAR-2 data but demonstrates increased DDM SNR that may be associated with surface inundation.There is a degree of overlap between the distributions of the Non-Inundated and all of the floodable classes, particularly between the top 75% of the Non-Inundated and the bottom 50% of the Inundated Forest, indicating that it may be difficult to differentiate flooded dense forest from non-inundated forest with GNSS-R signals alone.However, it also suggests there are potential regions of seasonally flooded vegetation that failed to be classified as flooded using the PALSAR-2 data but demonstrates increased DDM SNR that may be associated with surface inundation.
Remote Sens. 2018, 10, x FOR PEER REVIEW 19 of 29 differences in the GNSS-R signatures with respect to flooded area percentage observed by these four biomass categories underscores the need for land cover information to interpret the GNSS-R signal.

Comparison of SAR Classification and Peak DDM SNR at 1 km Scale
Analysis of multi-temporal DDM observations are focused on SNR since it exhibits the least noise of the three observations examined.DDM SNR values are posted to a 0.01° (~1 km) grid, with the time series maximum displayed over the study area in Figure 16a.The distribution of temporal maximum SNR in grid cells dominantly (>99%) classified as a single maximum flooded extent class with the PALSAR-2 data shown in Figure 16b.The median peak SNR value observed in Open Water grid cells is highest, at 13.97 dB, second highest in Inundated Non-Forest at 12.70 dB, and third highest in Inundated Forest at 8.54 dB, following the inverse relationship observed between peak SNR and amount of canopy attenuation above the flooded surface.The Non-Inundated class observed the lowest maximum SNR, with a median value among all dominantly Non-Inundated grid cells of 2.24 dB.There is a degree of overlap between the distributions of the Non-Inundated and all of the floodable classes, particularly between the top 75% of the Non-Inundated and the bottom 50% of the Inundated Forest, indicating that it may be difficult to differentiate flooded dense forest from non-inundated forest with GNSS-R signals alone.However, it also suggests there are potential regions of seasonally flooded vegetation that failed to be classified as flooded using the PALSAR-2 data but demonstrates increased DDM SNR that may be associated with surface inundation.interpretation of the impact of vegetative cover on inundation detection with GNSS-R.However, as L-band SAR may have limited sensitivity to inundation dynamics under dense, high biomass forest canopies, PALSAR-2 does not provide an ideal validation dataset.Ground data collected along the line transects made in Pacaya Samiria were used to offer insight into the SAR-based classification errors.Seventy-eight out of 81 transect points showed clear evidence of seasonal flooding.Of these 78 points, 31 (39.7%) are contained in pixels classified as "Non-Inundated" in the maximum flooded extent map derived from the PALSAR-2 ScanSAR mosaics.All three of the points that did not show evidence of past flooding were in pixels classified as Non-Inundated by the static PALSAR-based classification.Measured LAI and canopy cover of the 78 transect points (that showed evidence of past flooding) are illustrated in Figure 17a,b, along with the corresponding biomass density estimated by GFW-AGB shown in Figure 17c.In these boxplots, these vegetation structure metrics are grouped by the PALSAR-2 maximum flooded extent class.The 31 points where no flooding was ever detected by the PALSAR-2 classification time series but exhibited evidence of past inundation (shown in white) featured the highest values in all three vegetation metric categories (LAI, canopy cover, and AGB), further confirming the role of dense vegetation cover in under-classification of wetland area.
Another parameter to consider is maximum flood depth.Although it is expected that the height of the markings on trees left behind by past inundation events is generally representative of a past minimum inundation level, the depth observed by sustained inundation can be generally inferred from this value (although tree type can also have notable influence on this).Also, the presence of these inundation markings do not provide a direct indication of inundation timing, within or between years.The past inundation depths noted for the misclassified "Non-Inundated" transect points are markedly lower (median value of 0.5 m) than those observed by the classified flooded classes (median value of 1.2 m).A combination of dense canopy and comparatively shallow inundation depths may have contributed to under-classification of inundation in these areas by PALSAR-2.Another parameter to consider is maximum flood depth.Although it is expected that the height of the markings on trees left behind by past inundation events is generally representative of a past minimum inundation level, the depth observed by sustained inundation can be generally inferred from this value (although tree type can also have notable influence on this).Also, the presence of these inundation markings do not provide a direct indication of inundation timing, within or between years.The past inundation depths noted for the misclassified "Non-Inundated" transect points are markedly lower (median value of 0.5 m) than those observed by the classified flooded classes (median value of 1.2 m).A combination of dense canopy and comparatively shallow inundation depths may have contributed to under-classification of inundation in these areas by PALSAR-2.
Figure 18 illustrates a comparison of seasonal dynamics at four different 0.01 • (~1 km) grid cells with varying vegetative cover and inundation state."Inundated Swamp I and II" and "Seasonally Inundated Forest" represent regions A, B, and C, respectively, from the large transect made through a swamp shown in Figure 13."Terra Firme" represents an equivalently-sized grid cell identified in the non-inundated forest ROI K in Figure 4.As Figure 18a illustrates, a gradient of biomass density is observed across these four regions.All observed transect points in A and B were flooded during our visit in the dry season and are likely flooded throughout most of the year.However, C was observed to be partially non-inundated during the dry season but showed evidence of seasonal variation in inundation, as well as hosts much denser vegetation, and is thus referred to as "Seasonally Inundated Forest".The fourth region, referred to as "Terra Firme", represents an area known to rarely, if ever, flood and consists of high biomass density vegetation (GFW estimates ranging from 160-306 Mg ha −1 ). Figure 18b,c illustrate the mean DDM SNR and PALSAR-2-derived flooded area percentage estimates observed during the high water (March, April, May) and low water (July, August, September) seasons of 2017 for each of the four regions.Inundated Swamp I and II both demonstrate high flood signals in both seasons, consistent with a perennially flooded swamp.A small decrease in DDM SNR and SAR-derived flooded percentage is observed in the low water season.Terra Firme, on the other hand, demonstrates a low flood signal in both datasets.The Seasonally Inundated Forest exhibited a lower DDM SNR and SAR-derived flooded area percentage than both Inundated Swamp areas, with a notable decrease in signal observed in the low water season by both datasets, likely corresponding with a seasonal retreat of inundation.The higher biomass density found in the Inundated Forest has presumably contributed to the overall lower flooded signal compared to the signals observed in the inundated swamps, but it is difficult to separate the influence of biomass and existent inundation.The Seasonally Inundated Forest DDM SNR signal appears to be higher than that of the SAR-derived flooded area percentages relative to the inundated swamps' measurements.This difference may be indicative of greater sensitivity to flooded vegetation by GNSS-R, but also may be due to the non-linear response of GNSS-R to the presence of surface water.

Backscatter vs. Forward-Scatter
As demonstrated in numerous studies employing SAR over the past 30+ years, surface interactions between radar backscatter and inundated surfaces are generally well understood.Classification of L-band PALSAR-2 ScanSAR mosaics using established techniques were found to be generally very successful in detecting flooded surfaces, but did fail to detect inundation in several points in our ground data collection in transects of highly dense vegetation.One aspect not addressed in this study is the influence of the local incidence angle on backscatter.Because the SAR imagery used in this study was acquired using ScanSAR technology, the images cover a large swath width observed at a wide range of viewing angles (~18 • -70 • ).The imaging angle influences the radar scattering physics and may in turn impact the land surface classification.In this study, PALSAR-2 data were provided pre-assembled as multi-temporal mosaics with little flexibility in image selection.Classification of these SAR images can be potentially improved with custom assembly of ScanSAR datasets, as was used by Reference [6], taking advantage of optimal viewing geometries supporting landcover classification.
GNSS-R, on the other hand, provides a new mapping approach.The potential of using reflected GNSS signals for studying terrestrial surface properties has been established, but the exploitation of such signals is still in its early stages.Our results have found that three aspects of the DDM (peak SNR, LES, and TES) are all sensitive to both inundation state and vegetative cover in a tropical wetlands complex, with SNR exhibiting the lowest noise of the three.The SNR signal approaches saturation at relatively low percentages of flooded area in a given footprint (~20%), which strongly supports the sensitivity to the presence of surface water but presents a challenge in quantifying the amount of surface water.The influence of canopy attenuation can be delineated using ancillary information on above-ground biomass.Derivation of DDM LES and TES relies on the estimated scattering area of each DDM and the estimated location of each DDM's specular point, which depends highly on signal calibration.The CYGNSS Science Team is still assessing GNSS-R calibration techniques over land, so an element of the uncertainty found in these current measurements may be improved with future updates to the CYGNSS L1 dataset [52].The effects of incidence angle on these measurements have yet to be quantified and should be assessed in future work.Furthermore, in this study all surface reflections are assumed to be coherent over the wetland region.A calibration more accurately reflecting land cover composition would also account for contribution from the incoherent reflection.Radiometric modelling of reflected signals that allows assessment and interpretation of bistatic scattering physics should be conducted to further examine the degree of signal coherence in wetlands regions in addition to the influence of canopy attenuation.The role of biomass variability within a footprint on the DDM signal should also be explored in a modeling framework.
Validation of remote sensing-based, dynamic inundation datasets is immensely challenging due to a lack of large-scale in situ spatio-temporal information of inundation state.The very reason that makes these products so useful is what makes them so difficult to assess.Since CYGNSS footprints are much larger than imaging radar pixels, assessment of GNSS-R must rely on comparison of disparate remote sensing-based data at various spatial resolutions and with different temporal samplings.This study relied on generally successful, yet imperfect, classification of L-band SAR mosaics to assess GNSS-R.Comparison of SAR backscatter with the GNSS-R data at a radiometric level directly is problematic for complex wetland environments as radar backscatter can increase or decrease with inundation state depending on vegetation cover.Improvements can be made with access to greater amounts of ground information.Focused experiments that take advantage of controlled data acquisitions, such as one designed with aircraft instruments, combined with ground data collections and radiometric modeling, would limit sampling ambiguity and support direct comparison and assessment of SAR and GNSS-R technologies for mapping wetlands in complex environments.

Potential for a Merged Approach
Results from this study support the potential for a combined backscatter and forward-scatter approach to mapping wetlands that can leverage the combined benefits offered by SAR and GNSS-R.CYGNSS reflections resolved to a grid over wetland regions of interest can support estimates of the fractional coverage of open water and inundated vegetation at sub-weekly intervals (repeat time is inversely related to the grid size).Sub-grid cell information from contemporaneous SAR images can inform on vegetation cover and areas of open water and inundated vegetation, while reflected signals observed by CYGNSS can be used to estimate the total fraction of surface inundation, potentially using empirically-or radiometrically-derived relationships driven by biomass metrics.Such a technique may benefit from the improved capability of bistatic-scattered signal from GNSS-R by enhanced sensitivity to inundation under dense vegetation canopies.
Calculations of the temporal repeat frequency of CYGNSS data depend on the assumed footprint size.For example, the CYGNSS mission estimates that the mean revisit time over the ocean, where the spatial footprint is ~25 × 25 km, is approximately 7 h [44].Due to the finer spatial resolution over the land surface, the temporal repeat time is longer.This study focused only on CYGNSS DDMs processed at the first Fresnel zone, but a larger footprint area can be achieved with greater understanding of the role of signal coherence on the DDM processing.The coarser the spatial grid in use, the greater the temporal coverage.For example, a grid cell of ~9 × 9 km in the tropics is sampled on average every ~3 days.Further examination is needed to focus on determining the optimal spatial resolution at which to resolve the reflected CYGNSS measurements for terrestrial applications.

Conclusions
This study documents the comparative capabilities of L-band SAR and GNSS-R for mapping tropical wetlands extent and inundation over a large tropical wetlands region with complex landcover.Dual-polarized multi-temporal L-band SAR mosaics provided by ALOS2 PALSAR-2 offer a rich time series (repeat time of ~1-2 months) of interannual inundation dynamics in the Peruvian Amazon.Consistent with previous studies, ground measurements in areas of high biomass and canopy cover confirmed inundatable areas where inundation was not identified with SAR.
GNSS-R data from NASA's CYGNSS mission were analyzed in combination with the PALSAR-2-based wetlands assessments.Findings show that GNSS-R signals demonstrate a comparable sensitivity to the presence of surface water as standard techniques of classifying L-band SAR, with greater sensitivity possible but unconfirmed.A number of aspects of GNSS-R observations remain unanswered, including the roles of biomass variability and signal coherence on DDM observables.
One of the compelling characteristics of surface-reflected GNSS signals is their sensitivity to relatively high resolution (~1 km) changes in inundation extent with potentially high temporal coverage.Harmonization of GNSS-R and L-band SAR has the potential for broader application to further the goal of quantifying wetlands extent and inundation patterns at pan-tropical scales.This offers an opportunity for research in time-sensitive hydrodynamic systems, such as assessing inter-annual variations in local flood pulse intensity, understanding water exchanges between floodplains and rivers, or validating estimates of carbon exchange in wetlands.

Figure 1 .
Figure 1.(a) The Amazon Basin, illustrating primary river sub-basins and freshwater ecosystem extent, with the study region boxed in red.Data sources: open water and wetland delineation from Reference [31]; sub-basin boundaries and GTOPO30 digital elevation model (DEM) from Reference [32]; (b) Extent of the study area, bounded by 3-7°S, 73-77°W, containing the Pacaya Samiria National Reserve.Locations of biometry measurement sites 1-4 along Canal de Puinahua are labelled in black.Field station location where river gauge data are available are shown in red.

Figure 1 . 29 Figure 1 .
Figure 1.(a) The Amazon Basin, illustrating primary river sub-basins and freshwater ecosystem extent, with the study region boxed in red.Data sources: open water and wetland delineation from Reference [31]; sub-basin boundaries and GTOPO30 digital elevation model (DEM) from Reference [32]; (b) Extent of the study area, bounded by 3-7 • S, 73-77 • W, containing the Pacaya Samiria National Reserve.Locations of biometry measurement sites 1-4 along Canal de Puinahua are labelled in black.Field station location where river gauge data are available are shown in red.

Figure 3 .
Figure 3. Two-step classification scheme used to classify static maximum flooded extent and dynamic land cover types and flooding status.The decision tree was implemented with thresholds derived from ROIs corresponding to associated homogenous land cover classes.

Figure 3 .
Figure 3. Two-step classification scheme used to classify static maximum flooded extent and dynamic land cover types and flooding status.The decision tree was implemented with thresholds derived from ROIs corresponding to associated homogenous land cover classes.

Figure 3 .
Figure 3. Two-step classification scheme used to classify static maximum flooded extent and dynamic land cover types and flooding status.The decision tree was implemented with thresholds derived from ROIs corresponding to associated homogenous land cover classes.
Forest classes, separated by vegetation density informed by γ 0 HV−DS .The Open Water class is divided into Permanent Open Water and Intermittent Open Water classes based on minimum temporal open water extent detected.The five static classes derived by the end of Step 1 each had its own decision rule-set developed for dynamic classification in Step 2. Here, classification was performed on individual scenes in the time series describing the flooding status at a given point in time.Image segmentation was performed at both steps using the open-source RSGISLib Python module.Segmentation in Step 1 considered the following images: γ 0 HH−Min , γ 0 HH−Max , γ 0 HV−Min , γ 0 HV−Max , γ 0 HV−DS , and NDVI-DS.In Step 2, the images used for segmentation include an individual cycle's HHand HV-polarized backscatter, and the ratio of these two (HH/HV).

29 Figure 5 .
Figure 5. (a) Example of a bistatic radar cross section Delay Doppler map (DDM); (b) integrated delay waveform (IDW) derived from (a).The leading edge slope (LES) is estimated from the peak in the IDW and the two preceding points, while trailing edge slope (TES) is derived from the IDW peak and three following points.

Figure 5 .
Figure 5. (a) Example of a bistatic radar cross section Delay Doppler map (DDM); (b) integrated delay waveform (IDW) derived from (a).The leading edge slope (LES) is estimated from the peak in the IDW and the two preceding points, while trailing edge slope (TES) is derived from the IDW peak and three following points.

Figure 6 .
Figure 6.Variety of floodplains observed at study sites (shown in Figure 1b) in Pacaya Samiria at low water season.(a) View of river bank near Site #1, (b) open area in Site #1(in Figure1b) with sparse tree cover dominated by Pseudobombax munguba; (c) Dense canopy cover found in Site #2 (in Figure1b); (d) inundated aquatic shrub species (Montrichardia arborescens) found along margins of lake at Site #3(in Figure1b); (e) denser tree cover and understory found in Site #3 (in Figure1b); (f) inundated swamp found at Site #4 (in Figure1b).Examples of high flood marks on trees (shown red rectangles) indicate the level of previous inundation.

Figure 6 .
Figure 6.Variety of floodplains observed at study sites (shown in Figure 1b) in Pacaya Samiria at low water season.(a) View of river bank near Site #1, (b) open area in Site #1(in Figure1b) with sparse tree cover dominated by Pseudobombax munguba; (c) Dense canopy cover found in Site #2 (in Figure1b); (d) inundated aquatic shrub species (Montrichardia arborescens) found along margins of lake at Site #3(in Figure1b); (e) denser tree cover and understory found in Site #3 (in Figure1b); (f) inundated swamp found at Site #4 (in Figure1b).Examples of high flood marks on trees (shown red rectangles) indicate the level of previous inundation.

Figure 8 .
Figure 8.(a) Step 1, object-based classification algorithm for maximum flooded extent.Boxes with rounded edges indicate inputs; boxes with sharp corners indicate a process/decision; ovals indicate assigned classes.In the decision tree, bold red arrows correspond with "True" and bold black arrows with "False".(b) Maximum flooded extent classification map.

Figure 8 .
Figure 8.(a) Step 1, object-based classification algorithm for maximum flooded extent.Boxes with rounded edges indicate inputs; boxes with sharp corners indicate a process/decision; ovals indicate assigned classes.In the decision tree, bold red arrows correspond with "True" and bold black arrows with "False".(b) Maximum flooded extent classification map.

Figure 9 .
Figure 9. Step 2 object-based classification algorithm for delineating time series inundation state.Boxes with rounded edges indicate inputs; boxes with sharp corners indicate a process/decision; ovals indicate assigned classes.In the decision tree, bold red arrows correspond with "True" and bold black arrows with "False".

Figure 9 .
Figure 9. Step 2 object-based classification algorithm for delineating time series inundation state.Boxes with rounded edges indicate inputs; boxes with sharp corners indicate a process/decision; ovals indicate assigned classes.In the decision tree, bold red arrows correspond with "True" and bold black arrows with "False".

Figure 10 .
Figure 10.Time series of PALSAR-2 classified flooded area percentage (open water, inundated nonforest and inundated forest) in the surrounding area within a 5 km radius and river level measurements (left) and scatter plots of corresponding values fitted to a logistic model (right) of the four field stations with locations shown in Figure 1: (a) Flor de Punga, (b) Requena, (c) Nauta, (d) San Regis.

Figure 10 .
Figure 10.Time series of PALSAR-2 classified flooded area percentage (open water, inundated non-forest and inundated forest) in the surrounding area within a 5 km radius and river level measurements (left) and scatter plots of corresponding values fitted to a logistic model (right) of the four field stations with locations shown in Figure 1: (a) Flor de Punga, (b) Requena, (c) Nauta, (d) San Regis.

Figure 11 .
Figure 11.Plot-level AGB density (x-axis) versus GFW-AGB density (y-axis).1:1 line shown in solid black.Linear fit shown in dashed black, with equation and r-squared reported.
) presenting SNR, center row (Figure14e-h) presenting LES, and bottom row (Figure14i-l) presenting TES.The leftmost column displays the relationship observed between the percentage of all flooded area (Open Water, Inundated Non-Forest, and Inundated Forest) of all CYGNSS footprints and the respective DDM values.The center left column illustrates a subset of these footprints, showing only those with negligible amount of classified flooded vegetation (<5%) and focusing on the dependence of DDM measurements on the flooded area percentage of Open Water.The center right column features a different subset, containing footprints with a negligible amount of Open Water (<1%) and illustrates the relationship between the area percentage of flooded vegetation (Inundated Non-Forest and Inundated Forest) and DDM values.Lastly, the rightmost column is an even smaller subset of those footprints in the third column, consisting of those with negligible amounts of Inundated Non-Forest and Open Water (<1%) and thereby isolating the relationship between the area of Inundated Forest and DDM measurements.The CYGNSS footprint data exhibited a high level of noise so were binned in equal intervals (n = 1000) by classified flood percentage in each plot.

Figure 14 .
Figure 14.CYGNSS footprint flooded area percentage estimated with classified PALSAR-2 mosaics vs. CYGNSS DDM observations, including SNR in top row (a-d), LES in center row (e-h), and TES in bottom row (i-l).Leftmost column displays all footprints, where % flooded counts flooded area regardless of vegetation (Open Water, Inundated Non-Forest, and Inundated Forest).The center left column presents the subset of the CYGNSS footprints with negligible inundated vegetation (Inundated Non-Forest and Inundated Forest of less than 5%), focusing only on the percentage of Open Water in each footprint.The center right column displays the subset of footprints with negligible Open Water (<1%) and displays the percentage of inundated vegetation (either Inundated Non-Forest or Inundated Forest).The rightmost column shows the subset of footprints with negligible Open Water or Inundated Non-Forest (<1%) and displays the percentage of Inundated Forest.Median values of DDM values are displayed in black, derived from bins of equal size (n = 1000).

Figure 17 .
Figure 17.Distribution of vegetation structure metrics of transect points that showed evidence of seasonal flooding, grouped by PALSAR-2 maximum flooded extent class: (a) LAI; (b) canopy closure; (c) GFW-AGB density.

Figure 18 Figure 17 .
Figure18illustrates a comparison of seasonal dynamics at four different 0.01° (~1 km) grid cells with varying vegetative cover and inundation state."Inundated Swamp I and II" and "Seasonally Inundated Forest" represent regions A, B, and C, respectively, from the large transect made through a swamp shown in Figure13."Terra Firme" represents an equivalently-sized grid cell identified in the non-inundated forest ROI K in Figure4.As Figure18aillustrates, a gradient of biomass density is observed across these four regions.All observed transect points in A and B were flooded during our visit in the dry season and are likely flooded throughout most of the year.However, C was observed to be partially non-inundated during the dry season but showed evidence of seasonal

Figure 18 .
Figure 18.Comparison of high vs. low water seasons in 2017 at four 0.01 • grid cells.(a) Distribution of estimated GFW-AGB density observed in each grid cell; (b) observed CYGNSS DDM SNR during high vs. low water seasons in each region (number of observations ranged from 3-8 per season for each region); (c) estimated flooded area percentages derived from PALSAR-2 during high vs. low water seasons in each region (number of observations was 2 or 3 per season for each region).Bar charts (b,c) show mean values observed in seasonal window, with standard deviation indicated by error bars.

Table 1 .
Characteristics of selected ROIs and descriptive summaries of dry season NDVI and multi-temporal γ 0 statistics observed.

Table 2 .
Accuracy assessment of Step 1 maximum flooded extent classification compared with IIAP reference map of Pacaya Samiria National Reserve.

Table A2 .
CYGNSS L1 data quality flags used in this study.