Assessment of Forest Biomass Estimation from Dry and Wet SAR Acquisitions Collected during the 2019 UAVSAR AM-PM Campaign in Southeastern United States

: Forest above-ground biomass (AGB) estimation from SAR backscatter is affected by varying imaging and environmental conditions. This paper quantifies and compares the performance of forest biomass estimation from L-band SAR backscatter measured selectively under dry and wet conditions during the 2019 AM-PM NASA airborne campaign. Seven Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) images acquired between June and October 2019 over a temperate deciduous forest in Southeastern United States with varying moisture and precipitation conditions are examined in conjunction with LIDAR and field measurements. Biomass is estimated by fitting a 3-parameter modified Water Cloud Model (WCM) to radiometric terrain corrected SAR backscatter. Our experiment is designed to quantify the biomass estimation errors when biomass models are calibrated and validated on varying acquisition conditions (dry or wet). Multi-temporal estimation strategies are also evaluated and compared with single-acquisition estimation approaches. As an outcome, the experiment shows that the WCM model calibrated and validated on single acquisitions adapts to different soil moisture conditions with RMSD up to 18.7 Mg/ha. The AGB estimation performance, however, decreases with RMSD upwards of 30 Mg/ha when the model is cross-validated on moisture and precipitation conditions different than the calibration conditions. Results confirm that calibrating the model over the multi-temporal data using averaged backscatter or weighted combinations of individual AGB estimates, improves the biomass estimation accuracy up to about 20% at L-band. This study helps design biomass cal/val procedures and biomass estimation algorithms for dense time-series to be collected by low-frequency radar missions such as NASA-ISRO SAR (NISAR) and BIOMASS.


Introduction
Synthetic aperture radar (SAR) backscatter has been extensively explored for above-ground biomass (AGB) estimation across various forest ecosystems ranging from tropical to boreal [1][2][3][4][5][6][7][8][9]. The sensitivity of SAR backscatter to AGB has been modeled in many studies to estimate forest AGB using empirical [10][11][12], semi-empirical [13] and numerical [14,15] models, and machine learning algorithms [16][17][18]. SAR backscatter is sensitive to AGB up to a saturation level [19], although, for dense forests, a decrease in backscatter above saturation level may be observed [20]. The saturation at a specific biomass is determined by radar frequency along with other factors including polarization and forest conditions [2,4]. L-band SAR backscatter has been extensively used for forest AGB retrieval forest. The site has deciduous hardwood trees and woody wetlands prone to seasonal flooding each spring. The mean elevation of the site is 20 m above mean sea level. The dominant tree species are American sweetgum (Liquidambar styraciflua), American hornbeam (Carpinus caroliniana), and loblolly pine (Pinus taeda). Extensive field survey was carried out under the National Ecological Observatory Network (NEON) [33] to collect parameters such as tree DBH (diameter at breast height or 1.3 m above ground), species, and tree height. 41 plots measuring 20 m × 20 m (0.04 ha) have been surveyed between 2015 to 2019. Additionally, ten plots measuring 40 m × 40 m are surveyed in November 2019. The 5484 trees surveyed in these plots can be grouped according to species-groups listed by Jenkins et al. [34] as hardwood, woodland, and softwood. The distribution of these tree species-groups in the surveyed plots is hardwood species (71%), woodland species (11%), softwood species including conifers (2%) while the remaining are shrubs or are unclassified. The national land-cover database (NLCD) 2016 released by multi-resolution land characteristics consortium (MRLC) is also used to identify the major vegetation type for the plot. 54.5% of plots are in woody wetlands, 29% are in deciduous forests, 12% are in mixed forest while the remaining plots are in shrub/scrub forested areas. The geolocation of the plot corners is accurately determined using a dual-frequency GPS receiver. The DBH is measured for all trees with DBH greater than 2.5 cm. The height of a few trees in each plot is also measured. The DBH measured for each tree is converted to growing stock volume (GSV) and biomass using species-specific allometric equations and wood-specific density [34]. The aggregated plot-level AGB ranges between 0.05 Mg/ha to 373 Mg/ha with a mean AGB of 157 Mg/ha.
A permanent flux tower and network of soil moisture sensors near the flux tower measure precipitation, temperature, and volumetric soil water content or soil moisture (cm 3 /cm 3 ) at 0.01 m below the ground. The location of the field plots and the flux tower is shown in Figure 1, and the measured soil moisture and cumulative precipitation before each acquisition is shown in Table 1. The flux tower also has a phenology camera that captures images of the canopy at 15 min intervals. Figure 2 shows the images captured around one hour before the UAVSAR acquisition.   15 October 2019 (Acq7) Figure 2. Lenoir Landing site forest canopy images captured from the phenology camera at the flux tower. The images are captured approximately one hour before the acquisition of UAVSAR AM-PM campaign data. The precipitation is visible for the last acquisition on 15 October 2019 (g).

Lidar Data and Reference AGB Map
A reference AGB map is generated from small-footprint LIDAR-derived canopy height to mitigate errors that would result from calibrating SAR-AGB relationships directly with small plots [4]. Canopy height in meters above the bare Earth is generated at 1 m resolution using small footprint LIDAR data collected from NEON Airborne Observation Platform (AOP) LIDAR-Optech LIDAR Gemini instrument [33]. The plot-level AGB is related to LIDAR measured canopy height using a least-squares minimization regression fit using a power function [35]. This provides a reference AGB map at 1 m posting resampled to the UAVSAR SLC data in slant range. The reference AGB map generated using LIDAR data has RMSE of 17.5 Mg/ha for AGB less than 100 Mg/ha. The RMSE for higher AGB (up to 335 Mg/ha) regions is 48 Mg/ha. The extrapolation of field AGB to LIDAR provides wider spatial coverage and higher training samples [3]. 311 random samples (referred to as polygons) of approximately 1 ha are selected in the test site. Figure 1 shows the LIDAR height map and the surveyed 0.04 ha plots and LIDAR-upscaled 1 ha polygons. The AGB for these 311 polygons varies from 1.8 Mg/ha to 335 Mg/ha with a mean of 125 Mg/ha. The major NLCD forest-types represented by these 1 ha polygons are woody wetlands (43%), evergreen forests (26%), and mixed and deciduous forests (15%) while the remaining polygons are in scrub and herbaceous forests. Among these 311 polygons, 24 polygons (8%) are excluded from the analysis as these partially cover non-forest regions including roads and water-bodies.

AM-PM Campaign SAR Data
A total of seven L-band SAR acquisitions were collected over the Lenoir Landing forest site between 6 June to 15 October 2019 as part of the UAVSAR AM-PM campaign. The SAR backscatter for these data are collected under different acquisition conditions (soil moisture and precipitation) as detailed in Table 1. The data on October, 1 (Acq6) is acquired under dry conditions with no precipitation in the preceding 72 h (hours) and a very low soil moisture of 0.08 cm 3 /cm 3 or 8%. On the other hand, the data on October, 15 (Acq7) can be considered as a case of very wet acquisition as it was acquired during a precipitation event and a very high soil moisture of 0.48 cm 3 /cm 3 . Note that 0 h precipitation denotes rainfall measured during the approximately 5 min of acquisition time of the UAVSAR flight.
The UAVSAR data are acquired at 6-pm local time to match the NISAR ascending pass. The look angle of UAVSAR varies from 21 • in near range to 65 • in far range. However, the field plots and 1 ha polygons cover a smaller incidence angle range of 34 • to 50 • . The data are provided as SLC with azimuth resolution of 0.6 m and range resolution coarser than 1.87 m. Here, we examine the cross-polarimetric (HV) channel that is relevant to biomass estimation from systematic dual-polarimetric SAR missions like ALOS-2 and NISAR.
The UAVSAR data is radiometrically calibrated and co-registered using the 1-arc second SRTM DEM to form a stack of seven SLC images. The SLC images store radar backscatter as σ 0 (referenced to the flat ellipsoid), which is modulated by topography. An improved radiometric terrain correction (RTC) algorithm [36] is applied to flatten the backscatter and generate γ 0 . Our approach for RTC includes compensation of backscatter variation due to changing radar look angle across the imaged swath. The algorithm generates a flattened backscatter image ideally free of radiometric modulations caused by topography and incidence angle variations. In the RTC process, the plot corners are projected to radar geometry, and sub-pixel aggregation of backscatter is carried out. Further, this aggregation approach ensures that there is no mismatch between pixel size and plot size, as even partially covered pixels are weighted according to their areal coverage. The 1-ha backscatter averaging reduces noise and speckle from the full resolution SAR data.

Methodology
The 3-parameter Water Cloud Model (WCM) [37] is adopted in this research to estimate the forest AGB. The WCM model is a physically-based model that presents the SAR backscatter as a function of forest biophysical parameters. For modeling L-band SAR backscatter, a version of WCM is used which accounts for vertical and horizontal discontinuities or gaps in the forest canopy [38,39]. The WCM model is selected as it has a strong physical foundation as shown in [37][38][39] and has been used to map forest biomass from L-band backscatter in numerous studies [2,22,24,40]. In addition, the NISAR and BIOMASS missions are planning their cal/val and biomass mapping activities on models based on WCM given its simplistic form and wide applicability.
The WCM model expresses backscatter as the incoherent sum of ground (parameter a) and vegetation (parameter b) scattering contributions weighted by the attenuation due to a vegetation layer (factor e −c β ) where β is the above-ground biomass in Mg/ha. The derivation of (1) is not discussed here as it has been extensively discussed by [38,39]. The model parameters a, b, and c are estimated from SAR and lidar-based AGB map data after applying a χ-squared minimization approach using 137 1-ha polygons with AGB up to 100 Mg/ha.

Model Calibration and Validation
To analyze the impact of varying acquisition conditions, and the utility of multi-temporal data to reduce the AGB estimation errors, we consider the following distinct cases: • Single acquisition backscatter: The model is trained and validated using the seven L-band cross-polarized backscatter with different acquisition conditions (case 1). • Multi-temporal averaged backscatter: The temporal mean of HV backscatter over all the seven acquisitions (MT-all, case 2) and the temporal mean of HV backscatter over acquisitions without rainfall in the 24 h preceding the acquisition time (MT-24, case 3) • The multi-temporal weighted average (WA) of AGB is estimated from the seven acquisitions (case 4) with the weights explained below.
Case 1 addresses the first objective (see Section 1) of this study, while the other three cases (case 2-case 4) address the second objective and explore the ability of multi-temporal data to improve the AGB estimation performance. For the first three cases, the backscatter is modeled to estimate AGB, while for case 4, the modeled AGB from case 1 is the input to estimate the weighted average AGB. Hence, the WCM model parameters are reported for the first three cases.
To map AGB for the first three cases, the WCM model is calibrated using 60% randomly-selected polygons and validated over the remaining 40% polygons. The calibration and validation process is iterated 30 times with the data split randomly into training and test sets to allow robust estimation of the error. The WCM model is calibrated for each case by estimating the model parameters a, b and c. The root mean squared deviation (RMSD), coefficient of determination (R 2 ), saturation level, and number of valid solutions are used to compare the performance obtained for each case. The backscatter saturation level β sat for a given WCM model is obtained by deriving the value of biomass that gives a change in backscatter due to speckle equal to the backscatter associated with the desired biomass error (20 Mg/ha according to the NISAR requirements) [19]. For a WCM model characterized by parameters a, b and c, solving (1) for β = β sat provides the backscatter (γ 0 sat ) at which the WCM curve saturates. The valid solution space is the region between the modeled backscatter values ranging from a to γ 0 sat . In case a polygon has backscatter outside this valid range [39] a biomass value is assigned. Polygons with backscatter lower than the ground scattering component a are assumed to have 0 Mg/ha biomass and polygons with backscatter higher than γ 0 sat are assigned the maximum AGB. Multi-temporal data stack can improve the AGB retrieval results using a weighted combination of the AGB estimated from individual acquisitions as shown by [2,13,23,24,39]. The multi-temporal biomass, β mt is computed for a particular polygon from a weighted combination of biomass estimates (β i ) from the seven acquisitions as: The weights w i are calculated using two approaches (1) the difference between the parameters a and b, i.e., the dynamic range between backscatter from the ground and from the forest with infinite biomass (WA-DR); (2) the inverse of the soil moisture measured for each acquisition (WA-SM).

Temporal Cross-Validation of the Model
The robustness of the model trained and validated under different acquisition conditions is analyzed by considering the scenario of calibrating the model from each of the first three cases and validating it with all the other cases. The WCM model parameters (a, b, and c) obtained for each of the first three cases listed above are applied to randomly-selected 80% polygons in the remaining cases to generate all cross-cases combinations of calibration and validation of the AGB. This means, for instance, that the model generated for one date (Acq 1) is evaluated on other dates (Acq2-Acq7) and the multi-temporal mean backscatter (MT-all and MT-24). This step is iterated 30 times with randomly selected validation polygons to provide a robust estimate of cross-validation and model stability.

Analysis of Backscatter versus Biomass
The variation of backscatter with respect to forest AGB is analyzed as a function of forest type and with acquisition conditions such as precipitation and soil moisture. Figure 3 shows the HV backscatter and the WCM model curves as a function of AGB for all 7 individual acquisitions (case 1 in Section 3.1) along with the mean backscatter for all the seven acquisitions (case 2 in Section 3.1). The SAR backscatter increases with increasing forest AGB, especially for low forest biomass. Overall, for the seven acquisitions, the SAR backscatter-AGB trends agree with previous L-band analysis reported in temperate as well as other forest environments [21,23,[41][42][43][44][45][46][47]. Observation taken during significant rainfall event (Acq7) showed increased backscatter and lower dynamic range.
SAR backscatter magnitude changes with the dielectric constant as well as the surface roughness [48][49][50][51]. For a polygon representing a forested region, the backscatter varies over the seven acquisitions due to changes in soil moisture and precipitation. Overall, the backscatter increases for acquisitions during precipitation (Acq7) or with higher soil moisture (Acq1). Similarly, acquisitions with dry soil conditions (Acq6) have lower backscatter. Mean backscatter for Acq7 with wet soil conditions (0.48 cm 3 /cm 3 ) is −9.2 dB ± 1.5 dB, whereas for Acq6 with dry soil conditions (0.08 cm 3 /cm 3 ) it is about 2 dB lower (−11.1 dB ± 2.3 dB). The backscatter variation is also dependent on the forest AGB. The impact of precipitation and soil moisture is particularly visible in the low AGB range, where the WCM curve of Acq7 lies about 2 to 5 dB higher than the WCM curve of Acq6. The variation in backscatter is caused by the change in soil permittivity as a result of increased soil moisture. This is consistent with higher sensitivity for low AGB regions reported for tropical [27] and boreal forests [8].
The WCM model curves show that the model adapts to the variation in the backscatter in response to soil moisture or precipitation changes. We attempt to evaluate whether this backscatter change is influenced by the major forest-type of the polygons. The polygons covering herbaceous forests (see Figure 3) have low AGB (up to 20 Mg/ha) and the backscatter is more sensitive to soil moisture and precipitation due to sparser canopy cover. For Acq7, the backscatter for these herbaceous forests is about 2.23 dB lower. Polygons dominated by deciduous and mixed forests have AGB varying between 20 Mg/ha and 100 Mg/ha. These polygons show a uniform 1 dB decrease in backscatter for Acq7 and this decrease is uniform across the AGB range. The polygons dominated by woody wetlands also cover an AGB range between 20 Mg/ha and 100 Mg/ha. However, for Acq7 the backscatter decreases by almost 2 dB for low AGB (<20 Mg/ha) polygons and by 0.75 dB for polygons with AGB above 40 Mg/ha. This higher sensitivity of backscatter to lower AGB polygons is explained by the presence of wet soil or water in the woody wetlands and lower canopy density.

Estimates of WCM Parameters
Model training on the first three cases explained in Section 3.1 was carried out by estimating the WCM parameters a, b and c with a χ-squared minimization approach. Figure 4 shows the variation of WCM parameters with soil moisture for case 1. Parameters a, b and c represent the ground and canopy scattering contribution and the microwave extinction due to vegetation, respectively. Overall, with the increase in soil moisture, parameters a and b increase while c decreases. The parameters a and b are positively correlated with soil moisture having a correlation of 0.88 (p < 0.05) and 0.96 (p < 0.05), respectively. The parameter c is negatively correlated with soil moisture (r = −0.71), however, the correlation is not statistically significant within the 95% confidence interval with p = 0.07. Compared to other acquisitions, an additional increase in parameters a and b and a decrease in parameter c is observed for Acq7. This is attributed to high precipitation measured during the acquisition. The increase is highest in a with a 5 dB excursion from dry to wet conditions. This is expected as a represents the scattering contribution from the ground, and therefore it is more sensitive to changes in soil water content. The dynamic range of the WCM model decreases from 6 dB to 5 dB with an increase in soil moisture from 0.08 cm 3 /cm 3 to 0.40 cm 3 /cm 3 . This reduction in dynamic range translates into smaller signal attenuation in the canopy modeled by the WCM model (parameter c). The dynamic range is a measure of the sensitivity of the backscatter to forest AGB [52]. Interestingly, for Acq7, the dynamic range of the WCM model reduces to 3.5 dB which shows the impact of precipitation in addition to soil moisture increase. As seen in Figure 3, the impact of precipitation on Acq7 backscatter is higher for low AGB polygons and reduces as AGB increases. Therefore, precipitation appears to reduce the range of the backscatter measured over the polygons and the WCM model's dynamic range. The effect of these variations in WCM parameters on the estimated AGB is analyzed in the next section.

AGB Retrieval Performance
The estimated WCM parameters are used to invert the WCM model and estimate the forest AGB. The section first analyzes the AGB retrieval performance using single acquisitions (case 1) followed by multi-temporal average backscatter (case 2 and case 3) and multi-temporal weighted average of AGB (case 4). The AGB co-validation of models calibrated and validated for each backscatter case of Section 3.1 is shown in Table 2. The corresponding validation scatter-plots between the LIDAR AGB and WCM modeled AGB are shown in Figure 5. The figure clearly shows that the AGB retrieval performance varies for each acquisition. At 1 ha scale, AGB estimation from single acquisitions (case 1) leads to RMS deviation with respect to LIDAR AGB between 14.5 Mg/ha and 18.7 Mg/ha. Data acquired during precipitation and wet soil conditions (Acq7) has the highest error with RMSD of 18.7 Mg/ha. With multi-temporal average backscatter, the RMSD is 16.1 Mg/ha (MT-all, case 2) and further reduces to 14.4 Mg/ha (MT-24, case 3). From Figures 4 and 5, it is observed that the dynamic range of the WCM model is related to the soil moisture and a lower dynamic range results in higher RMSD. These results are consistent with previous findings [2,47]. Excluding acquisitions with precipitation (MT-24) reduces the backscatter variations and improves the AGB estimate by 12% compared to the mean RMSD of individual acquisitions.    Figure 6 illustrates the advantage of having multi-temporal AGB estimates. The dashed line shows the RMSD for all seven single acquisitions (case 1) while the solid lines represent the improvement in RMSD when each acquisition is successively added to the multi-temporal weighted average. The multi-temporal weighted average of the AGB estimates from all seven acquisitions in case 1, weighted by the inverse of the soil moisture (WA-SM) or by the dynamic range (WA-DR) have higher accuracy with RMSD of 14.13 Mg/ha and 14.17 Mg/ha respectively. Although both the parameters used as weights for averaging result in the similar performance, it is important to note the advantage of using dynamic range over regions where soil moisture products are not available at sufficient temporal or spatial resolutions.
It is important to also discuss the impact of multi-temporal weighted average AGB estimates on valid solution range. As discussed in Section 3, the polygons with backscatter outside the range defined by parameter a and b are assigned an AGB value of either 0 Mg/ha or maximum biomass. However, this adds further uncertainty in AGB estimates as the true value of AGB is not known. The percentage of polygons with backscatter outside the dynamic range for each of the above cases is computed. It is observed that the weighted average AGB estimates have the least number of polygons outside the dynamic range (7%) while the corresponding values for single acquisition ranges between 24% to 16%. This demonstrates the impact of varying acquisition conditions on AGB estimation and the advantage of multi-temporal acquisitions.

Factors Influencing AGB Retrieval
SAR backscatter is influenced by (1) forest structural parameters including forest density, tree species, shape of canopy, and (2) environmental parameters including soil moisture and precipitation [22].
Conifers and deciduous species have distinct canopy shape. However, the Lenoir Landing site is dominated by hardwood and woodland species with a small proportion (2%) of the 0.04 Ha plots dominated with coniferous species (see Section 2.1). The effect of acquisition conditions (soil moisture and precipitation) on WCM modeling parameters and AGB retrieval are discussed earlier. The accuracy of AGB retrieval across different AGB regions and for different NLCD forest types are analyzed in this section. Figure 7 shows the AGB inversion RMSD box-plots (mean, quantiles, and minimum and maximum RMSD values) for the seven single acquisitions over five AGB ranges. The AGB class boundary is selected such that equal sampling of polygons is possible across all regions. The dashed and dash-dotted lines show the RMSD for the weighted average RMSD. The advantage of having multi-temporal acquisitions and using a weighted average RMSD can be appreciated. The RMSD increases with AGB and stabilizes around a mean of 22.5 Mg/ha for case 1. With WA-SM and WA-DR, the RMSD is much lower at 17 Mg/ha which shows a significant improvement in AGB retrieval performance over single-acquisitions. The advantage of using multi-temporal weighted average AGB is clearly observed, especially for higher AGB regions.

NLCD Forest Type
The NLCD provides a broad perspective regarding the major forest canopy type within the study area. Figure 8 shows the box plot of the RMSD for the seven acquisitions over five major NLCD classes. The colors represent the mean AGB for each NLCD class. Forests in all the major NLCD classes have RMSD below 25 Mg/ha with the woody wetlands having a higher deviation in AGB retrieval. This might be due to the presence of under-canopy wet soil not accurately modeled by the WCM model. The low-AGB forests (grassland NLCD class) have a relatively higher RMSD compared to their mean AGB. In these low-AGB regions, the changes in ground characteristics have a higher impact on SAR backscatter leading to higher RMSD.

Temporal AGB Cross-Validation
The RMSD from all cross-validation cases are computed and shown by the box-plot in Figure 9. For Acq7, cross-validation leads to the highest RMSD as backscatter is affected by precipitation.
The cross-validation RMSD ranges between 22 and 33 Mg/ha in contrast to the co-validation RMSD of 18.70 Mg/ha. This shows that due caution is needed when applying models trained under wet acquisition conditions across multiple acquisitions. The model trained with Acq6 resulted in invalid solutions for about 60% of cross-validation samples due to saturation of the backscatter-AGB curve. This suggests that the valid inversion backscatter interval needs to be considered for cross-validation. Backscatter does not increase significantly for biomass values beyond the saturation level β sat (see Figure 3). γ 0 sat is the backscatter at which the WCM curve saturates (see Section 3.1) and provides the upper limit of backscatter that can be inverted within a pre-defined error limit using the WCM model. This is important for cross-validation as the backscatter can vary between acquisitions due to changes in environmental conditions as shown in this analysis. The WCM curve saturates at −8.0 dB for acquisition with wet soil conditions (Acq7, see Figure 3) and at −9.8 dB for acquisition with dry soil conditions (Acq6). This implies that when the model calibrated using Acq6 is validated over other backscatter cases, all polygons with backscatter higher than −9.8 dB will not result in a valid solution. The percentage of polygons with valid solutions is shown by the color gradient in Figure 9.

Limitations of the Study
The results presented in this section show the impact of dry and wet acquisition conditions on AGB estimation performance and the advantage of multi-temporal data to mitigate this impact. It is essential to discuss a few limitations related to the data and methodology.
First, as mentioned in Section 3.1, we use data acquired during dry and wet imaging conditions. These conditions are assessed from precipitation and soil moisture measurements available from one location (flux tower). The conditions can deviate locally from those observed at the flux tower. However, as these measuring sensors are within the forest site, the spatial variability is reduced.
When using multiple SAR acquisitions, bias can be introduced due to temporal changes in forest or soil. In this study, this is minimized as all the acquisitions are collected within the summer and beginning of fall 2019. As no acquisitions are collected during the spring season when flooding occurs, the impact of standing water on the backscatter-AGB relationship is not analyzed in this paper. Extending this analysis with multi-temporal SAR data acquired through all seasons would provide a more extensive data-set to capture the combined impact of weather and seasonal changes (leaf on/off) on backscatter-AGB relation. The NISAR mission will provide dense time-series similar to the one used in this analysis, capable of generating multiple acquisitions over different seasons or growth cycles unique to each forest eco-region. This analysis is limited to forests with AGB less than 100 Mg/ha in agreement with the upper limit established by the NISAR mission science requirements.
The reported AGB accuracy does not include errors from field-estimated and LIDAR-derived AGB, which may be caused by uncertainties in allometry, geolocation, and relation between LIDAR Canopy Height Model (CHM) and field AGB. Chave [29] reported errors of 5% for allometric equations over tropical forests with field plots of 1 ha. This error is expected to be higher for smaller field plots. The GPS geolocation uncertainty for NEON field plots is below 1 m. For the small plot size (0.04 ha), this can lead to errors in highly heterogeneous forests when the field AGB is used to generate LIDAR AGB maps. For the Lenoir Landing site, the forest in most field plots is homogeneous, so we expect this error to be minimal.
Furthermore, the inherent difference in viewing geometry of LIDAR and SAR will impact the accuracy of LIDAR AGB-SAR backscatter relationship when LIDAR AGB is projected to slant-range geometry. The airborne geometry, motion compensation and processing may introduce small artifacts in the data. The accuracy of reference LIDAR AGB depends on the ability of field plots to represent the species diversity within the study site. For forests with higher heterogeneity within plots, the intra-plot tree-density distribution would lead to a higher variance of observed backscatter for similar AGB values. This leads to a larger variance in model coefficients and higher error in AGB estimates. These errors, however, do not affect the conclusions drawn from this analysis. To obtain an absolute AGB accuracy, propagation of errors, which is out of the scope of this paper, should be carried out [4,19].

Conclusions
This study reports quantitatively the impact of backscatter perturbations due to changes in acquisition conditions on AGB estimation performance over a temperate broadleaf forest using new airborne campaign data acquired during the 2019 UAVSAR AM-PM Campaign. The advantage of using multi-temporal data to reduce the AGB estimation errors is also assessed and confirmed in the context of upcoming SAR missions. Precipitation affects AGB retrieval accuracy with up to 20% increase in AGB estimation errors. The WCM model parameters are highly correlated with measured soil moisture with correlation (r) of 0.88 and 0.96 for parameters a and b, respectively. Furthermore, soil moisture changes do not significantly affect the co-validation AGB estimation accuracy. Models calibrated and validated using acquisitions with significantly different meteorological conditions provide erroneous AGB estimates, with the RMSD increasing by up to a factor of 2 compared with co-validation RMSD. This emphasizes the challenges of calibrating backscatter-AGB models coping with soil moisture changes and precipitation and expected to be valid over multiple acquisitions in a time-series. This is evident in the temporal cross-validation analysis shown in this study. Capitalizing on the time-series data, we used the weighted average of the AGB estimated from multiple acquisitions weighted by the dynamic range and soil moisture to reduce the impact of changing acquisition conditions. The analysis can be further expanded over other forest eco-regions and with the use of a denser SAR time-series data-set.