Retrieval of Salt Marsh Above-Ground Biomass from High-Spatial Resolution Hyperspectral Imagery Using PROSAIL

: Salt marsh vegetation density varies considerably on short spatial scales, complicating attempts to evaluate plant characteristics using airborne remote sensing approaches. In this study, we used a mast-mounted hyperspectral imaging system to obtain cm-scale imagery of a salt marsh chronosequence on Hog Island, VA, where the morphology and biomass of the dominant plant species, Spartina alterniﬂora , varies widely. The high-resolution hyperspectral imagery allowed the detailed delineation of variations in above-ground biomass, which we retrieved from the imagery using the PROSAIL radiative transfer model. The retrieved biomass estimates correlated well with contemporaneously collected in situ biomass ground truth data ( R 2 = 0.73). In this study, we also rescaled our hyperspectral imagery and retrieved PROSAIL salt marsh biomass to determine the applicability of the method across spatial scales. Histograms of retrieved biomass changed considerably in characteristic marsh regions as the spatial scale of the imagery was progressively degraded. This rescaling revealed a loss of spatial detail and a shift in the mean retrieved biomass. This shift is indicative of the loss of accuracy that may occur when scaling up through a simple averaging approach that does not account for the detail found in the landscape at the natural scale of variation of the salt marsh system. This illustrated the importance of developing methodologies to appropriately scale results from very ﬁne scale resolution up to the more coarse-scale resolutions commonly obtained in airborne and satellite remote sensing.


Introduction
Coastal salt marshes form a critical transitional zone between land and sea and provide numerous ecological benefits including the removal and transformation of nutrients, essential habitat, protection from storm surges, support for coastal fisheries, and carbon sequestration [1][2][3][4][5][6]. The cumulative carbon stored in the biomass and deep sediments of vegetated coastal systems, including salt marshes-"Blue Carbon"-plays a critical role in global carbon sequestration [7,8]. However, primary production and carbon storage in marshes are highly variable over small spatial scales, leading to substantial uncertainty in total carbon sequestration [9]. Further, even in remote areas, these ecosystems are vulnerable to anthropogenic stressors, with rising sea-levels and shifting climates acting as a chief contributor to change in primary production and in some cases total marsh loss [8,[10][11][12][13][14][15]. Coupled with other environmental stressors, such as coastal eutrophication, significant marsh loss is occurring [16][17][18]. Hence, it is crucial to have tools to assess both spatial and temporal variability in plant biomass and marsh health to predict future loss scenarios and guide conservation decisions.
Vegetation community structure and aboveground plant biomass (AGB) can be a good indicator of marsh health [15,19,20]. Although field data and in situ measurements are essential for monitoring change, it is not feasible to cover large spatial extents using these methods alone. Given this limitation, hyperspectral remote sensing can play an important role in mapping ecosystems over large spatial areas, and advances in remote sensing have significantly enhanced our ability to monitor wetland environments. Remote sensing data have been used, for example, to map species composition [21][22][23][24], classify vegetation [22][23][24][25][26][27][28][29][30], and monitor marsh dieback [31][32][33][34]. Mapping of change in vegetation communities over time has been useful to predict future prospects for marsh survival, even when conducted at the coarse scale of satellite imagery [35][36][37][38]. However, in marshes of the Western Atlantic, where the low marsh zone is dominated by Spartina alterniflora, biomass can change by an order of magnitude over short distances due to differences in invertebrate densities [39,40], elevation [41,42], or die off [17,18]. This is particularly true in the narrow fringing back-barrier marshes that line the Atlantic Coast of North America [42]. New methods are needed to capture the variability in salt marsh structure on the scale at which it occurs and provide accurate assessment of marsh condition. The critical scale on which change occurs requires centimeter-scale accuracy, and, in this study, we used a hyperspectral imaging system that provides imagery of a salt marsh system at resolutions varying from sub-centimeter-to centimeter-scale spatial resolution.
Statistical approaches and vegetation indices have been used to characterize carbon storage within different ecosystems (including wetlands) [43][44][45]. In some cases, estimates of plant biomass have been made [46], but with spatial resolution on the order of 1 m pixel size. In spatially heterogeneous marshes, even this scale may be too coarse to capture the detailed variability inherent in these systems and derive accurate models of biomass. Physics-based radiative transfer models, such as PROSAIL, have also been widely used to characterize biophysical properties of vegetation [47][48][49][50][51][52][53][54][55][56][57][58]. Physics-based models are capable of simulating the propagation of light and interaction with plant canopies, and interpreting vegetation reflectance in terms of plant structure and biophysical parameters [48]. The inversion of these physics-based models allows for targeted retrieval of different biophysical parameters. In this study, we inverted the radiative transfer model PROSAIL using optimization techniques to retrieve the aboveground leaf biomass from salt marsh vegetation, typically dominated by Spartina alterniflora.
For this study, we collected data using a high-efficiency hyperspectral imaging (HSI) sensor mounted on a telescopic mast [59]. The HSI system is capable of collecting imagery at very fine spatial resolutions (∼0.5-3 cm), which can provide the level of detail required for more accurate mapping of wetland ecosystems where the above-ground biomass varies substantially over small spatial scales due to variability in age, elevation, and proximity to die-off zones. In contrast to drone-based imaging systems, the oblique imaging geometry of our mast-mounted imaging system provides a unique view of the wetland canopy allowing in many instances access to views of both lower and upper canopy elements. No view of the canopy is without the potential for occlusion, thus, in the oblique view, parts of one plant may block another. Similarly, for an erectophile canopy, nadir views from a drone may be limited by the fact that upper canopy elements may obscure lower canopy elements and sediment may become a large part of the imagery with little of the plant surface directly exposed to the imaging system. Drone imagery is, of course, directly comparable to the nadir views from satellite systems, however, in some satellite imaging systems, off-nadir geometries are possible offering the possibility of oblique views of vegetation canopies. Our mast geometry for this study was limited to ±34 • degrees from the horizon, however, in the future, we plan to modify the configuration to obtain a broader range of viewing geometries that can be compared more directly with satellite oblique views.
The specific objectives of our study were to: (1) retrieve biophysical parameters from inversion of the PROSAIL model; (2) compare retrieved results with ground-based biophysical measurements; (3) evaluate how the model performs across gradients of marsh elevation, age, and heterogeneity; and (4) consider the effect of spatial scale on retrieved estimates of biomass.

Study Area
The work was conducted at the the Virginia Coast Reserve Long Term Ecological Research site on the Atlantic coast of the Delmarva Peninsula [14,41,42,60]. Hog Island (37 • 25 5.91 N, 75 • 41 36.71 W), a 10 km long by 2.5 km wide barrier island roughly 14 km from the Peninsula, is typical of Western Atlantic barrier islands ( Figure 1). The island, predominantly uninhabited since the mid-1900s and owned by The Nature Conservancy, is bounded on the eastern shore by the Atlantic Ocean and on the western shore by Hog Island Bay, a shallow coastal lagoon [61]. S. alterniflora is the dominant macrophyte in the marsh that fringes the lagoonal shore of the island, where it forms a classic zonation pattern of tall form plants in the low marsh, and short form in the mid and high zones [60]. In 1962, a strong nor'easter deposited nearly 1 m of sand on the southern end of the island, burying the existing marsh and elevating the intertidal platform [62]. Vegetation recolonized at different rates across the overwash fan, creating a chronosequence of different aged marshes. We define marsh age as the time at which S. alterniflora was first detected in aerial imagery. At the time of this study the younger marshes varied in age from 6 (established in 2011) to 43 years (established in 1974) old. The oldest marsh in our study was present on navigational maps dating back to at least 1850, and was thus at least 170 years old in 2017 [42,63]. The chronosequence marshes are experiencing both erosion [12,14] and expansion, and are punctuated by patches of marsh die off (Tyler pers. obs). The resulting heterogeneity in marsh age and health makes this site ideal for the evaluation of biomass, as the aerial biomass of S. alterniflora varies by more than two orders of magnitude [41,42,60,63]. The research presented here primarily focused on the low to mid marsh zones, and on the retrieval of S. alterniflora biomass. A major field campaign was conducted in the summer of 2017 to characterize and map the salt marsh ecosystem in southern Hog Island, as shown in Figure 1. The long-term objective of the study is to accurately model carbon storage within salt marsh environments across various spatial scales using very fine-scale hyperspectral imaging. Within the context of the larger team effort, the focus of the work presented herein was to estimate the aboveground biomass of S. alterniflora using PROSAIL across the different aged marshes in the southern part of Hog Island. The sampling sites during this first of several campaigns appear in Figure 1b. These spanned all three marsh ages shown, but within the available time and manpower of the first campaign, our sampling initially emphasized the marshes from 1850 and 1974. Sampling within individual marshes was designed to cover the spatial variations in biomass observed in the marsh system, and, where possible, we aimed to have 10-15 field ground truth points per hyperspectral scene. Future work will emphasize data taken from multiple years and probe the relationship between biomass and marsh age more completely.

In Situ Measurements
We collected ground truth biomass data in 48 plots spanning the different aged marshes in the field of view of our hyperspectral imaging system. In all ground-truth plots, the total number of culms and height to the tallest leaf of the first 10 culms along a diagonal transect through the plot were measured within a 0.25 m 2 or 0.125 m 2 quadrat. The plant height (mean ± standard deviation (minimum-maximum)) for the 1850, 1974 and 2011 marshes was 42.7 ± 18.6 (5-96), 30.6 ± 11.8 , and 55.8 ± 10.0 (18-86), respectively. A height/weight regression was applied to these measurements to obtain total aboveground biomass per plot, as described in [42,64]. Briefly, in three 0.25 m 2 randomly selected plots per marsh site, all S. alterniflora culms were clipped and cleaned of sediment. The height to the tallest leaf of 10-15 randomly selected culms was measured, and individual plants were dried at 60 • C for 24 h to obtain individual dry weights that were then used to develop allometric regression equations for each marsh. The remaining culms were dried in bulk to verify the regression. The regressions relating culm heights to the dry weights for the three different marsh ages are defined by the following equations: (1) where y represents the individual culm biomass and x is the culm height.

Design and Instrumentation
A major part of the 2017 field campaign focused on imaging the salt marsh vegetation using a state-of-the-art hyperspectral imaging system [59]. The hyperspectral imager, a Headwall VNIR Micro-hyperspec High Efficiency E-series [65], is a pushbroom system providing spectral measurements with 16-bit dynamic range from 400 to 1000 nm with 371 spectral bands and 1600 across-track spatial pixels. Figure 2 illustrates the different components and imagery collected by the system. The HSI sensor resides within a maritime-rated General Dynamics Vector 20 pan-tilt unit [66]. The pan-tilt unit can cover the angular range from −34 • to 34 • in pitch, and −175 • to 175 • in yaw. In practice, any azimuthal view orientation can be obtained simply by physically rotating the mast; however, full 360 • coverage of the landscape was possible without rotating the mast because the field of view of our imaging system was 57 • with the 12 mm lens that we typically used. The imaging system also hosts a Vectornav VN-300 IMU-GPS [67] to provide pointing information and GPS time-stamps. In the field, the HSI system is placed on top of a BlueSky AL-3 (15 m) telescopic mast [68], which is designed to be raised 1-15 m above the ground [59].
From the pan-tilt unit mounted on the telescopic mast, the hyperspectral imaging system can capture scenes from multiple viewing geometries, which allowed us to perform hyperspectral Bidirectional Reflectance Distribution Function (BRDF) measurements from the imaging sensor [59]. The ability to perform measurements at different heights, using the telescopic mast, also allowed us to observe the effect that spatial resolution has on hyperspectral imagery of a particular scene. The sensor can also be used as a low-rate video system, providing hyperpsectral measurements at frame-rates (rate of acquisition of a line of hyperspectral pixels) at up to 250 Hz, leading to full hyperspectral scenes of size 1600 across-track × 212 along track spatial pixels and 371 spectral bands at low video rates of 1.5 Hz [59]. The hyperspectral scenes described in this paper were acquired at slower scan rates and longer integration times and typically had larger along-track spatial pixel dimensions.

Image Acquisition and Processing
We acquired hyperspectral imagery using the Headwall system described in Section 3.1 during the 2017 field survey contemporaneously with the in situ measurements. The portable configuration of our system enabled us to collect data over a large spatial extent in the southern part of Hog Island. The system captured imagery from areas with three different marsh ages: ≥170, 56, and 7 years. The HSI system, which is inside a pan-tilt unit and mounted upon a telescopic mast, acquired imagery from multiple viewing geometries and at various spatial resolutions. The hyperspectral imagery captured had a ground sampling distance (GSD) typically between 0.5 and 3 cm, which provided a level of detail not usually available with most traditional remote sensing systems. The instrument design and collection method of the system produced varying GSD for each along-track pixel captured. As the line scanning unit was swept vertically across a scene, pixels capturing the foreground (bottom of image) had smaller GSD (0.5 cm) than those capturing information in the background (top of image), where GSD was 3 cm. This high spatial resolution is crucial for studying wetland ecosystems, which have such high heterogeneity.
Our HSI system was characterized in our lab at RIT using our LabSphere Helios Integrating Sphere, which is coupled to an on-board calibrated spectrometer to acquire calibrated spectral radiance (Wm −2 sr −1 nm −1 ) from the raw digital number (DN) collected by the Headwall imaging system [59]. The placement of Spectralon ® panels within the scene enabled the conversion of the imagery into surface reflectance. Using a Garmin RTK GPS TRM55971, we also conducted GPS surveys of fiducials placed in each scene, allowing us to geo-register the imagery. The number of ground control points (GCPs) used to geo-register the images varied between 4-12. The HSI images were registered using either first-, second-or third-order polynomial equations. The registration accuracy of the geo-referenced images was within approximately 1-10 cm, depending on the number of fiducials present within the scene.

PROSAIL: Model Overview
The radiative transfer model PROSAIL combines the canopy reflectance model SAIL with the leaf optical properties model PROSPECT [48,[69][70][71]. The model simulates the hemispherical reflectance of plant canopies from the visible near-infrared (VNIR) to the short wave-infrared (SWIR), 400-2500 nm, with a spectral resolution of 1 nm.
The PROSPECT model, based on the plate model, simulates canopy spectra as a function of leaf mesophyll structure (N), chlorophyll a+b concentration (C ab ), leaf water content (C w ), leaf mass area (LMA), brown pigments (C bp ), and total carotenoid content (C cx ) [69,72,73]. SAIL simulates the bidirectional reflectance of the canopy using the following inputs: leaf reflectance obtained from PROSPECT, leaf inclination distribution function (LIDF), leaf area index (LAI), the sensor and solar geometries (zenith and azimuth angles), the hotspot effect (hspot), and the soil brightness factor (ρsoil) [48,70,71]. The factor, ρsoil, is a scaling parameter that accounts for changes in the brightness of the soil due to the presence of moisture and soil roughness [74,75]. The LIDF is a key parameter describing the angular orientation of leaves in vegetation, and is an important variable for characterizing the plant structure. The LIDF generally describes six different types of plant structure: planophile, erectophile, plagiophile, extremophile, spherical, and uniform. In the model, the LIDF is defined by two different parameters, LIDFa and LIDFb, which describe the average leaf slope and the distribution bimodality, respectively [76].

PROSAIL: Inversion Methodology
In this study, we inverted the PROSAIL model to retrieve the above ground biomass from hyperspectral data collected during the 2017 field campaign. Our approach for the inversion of PROSAIL used a simplex-based search method, the Nelder-Mead Simplex Method [77]. This optimization technique minimizes the residual between the estimated reflectance from the PROSAIL model and the measured reflectance from our mast-based HSI system. The ranges (minimum and maximum) of the input parameters required for the inversion of the PROSAIL model are summarized in Table 1. The optimization technique attempts to find the best combination of input values for the PROSAIL model that minimizes the residual between the estimated and measured reflectance. Based on assessments from our field campaign, we found that the leaf area distribution of S. alterniflora has an erectophile structure. This was modeled in PROSAIL by setting LIDFa = 1 and LIDFb = 0. We set the C bp to zero in our inversion, since the data were collected during the growing season. The hspot was defined in terms of the LAI (hspot = 0.5/LAI) [78,79]. The range of other variables in the PROSAIL model were selected based on previous studies and field measurements. We estimated the AGB as the product of the leaf mass area and leaf area index (LMA × LAI) [80]. Over the years, many studies have demonstrated the use of only a limited number of spectral bands to efficiently describe the reflectance spectra of various materials, e.g., soil, water, vegetation, and snow [81][82][83][84]. Optimal spectral sampling has also been widely used in various studies involving PROSAIL [81,[85][86][87][88]. In our work, we had a similar objective of finding optimal wavelengths to perform our inversion to retrieve above ground biomass. We investigated the sensitivity of the PROSAIL model to variations of the different input parameters reported in Table 1. The use of sensitivity analysis (SA) provides us the opportunity to observe the dependence of models on each individual input parameter. SA helps us to understand model behavior and its fidelity in describing real-world phenomenology [89]. We simulated PROSAIL a total of 70,000 times for the range (minimum-maximum) of its model parameters. The SA used in the study was the Extended Fourier Amplitude Sensitivity Test (EFAST), which has been applied to remote sensing data including in PROSAIL studies [89][90][91][92]. The EFAST computes the total contribution of each model parameter by considering the associated variance of the model output (in our case this is the variability in the reflectance) [89]. The total contribution includes both the effect of the individual input parameter as well as the interactions between all of the parameters.
The results from the EFAST are shown in Figure 3. We report the SA only for wavelengths in the VNIR (400-1000 nm), since this was the operational wavelength range of our hyperspectral imaging system. The two model parameters of interest in our analysis were the LMA and LAI. The impact of LAI was most noticeable in the wavelength range of 400-500 nm, accounting for 65-75% of the variation. In this wavelength range, the impact of the other variables is significantly reduced (typically less than 25%). Similarly, in the wavelength range of 720-1000 nm, the LMA was expressed most strongly compared to other model inputs, accounting for 50-75% or more of the variability, while the contribution due to other inputs were again typically well below 25%. This sensitivity analysis provided us with a range of wavelengths (400-500 nm and 720-1000 nm) to perform our inversion, which integrated 235 out of the 371 bands of the mast-based HSI system used in this study. By incorporating the optimal wavelength ranges into the inversion methodology, we avoided common pitfalls typically observed with optimization techniques, such as the risks of local "false" minima in the parameter space. The inversion scheme to retrieve aboveground biomass as well as our methodology to collect/process hyperspectral imagery is illustrated in Figure 4.  Overview of the inversion scheme to retrieve aboveground biomass of Spartina alterniflora from imagery collected by a ground-based HSI system. After data preprocessing, the NDVI threshold mask is used to remove soil pixels from the image in the analysis. EFAST provides band selection, while optimization of the PROSAIL model to obtain LAI and LMA uses the Nelder-Mead algorithm [77].

PROSAIL: Spatial Resampling
In this work, we also wanted to illustrate the effects of spatial resolution on the retrieved PROSAIL salt marsh biomass. To assess the effects of spatial scaling, we binned collected imagery to reduce image resolution by factors of 4 and 100. These scale factors allowed us to consider, respectively, the effects of both relatively moderate vs. large changes in spatial scale. To preserve global scaling across an image, each row was binned independently based upon a predefined window size. To reduce spatial resolution by a factor of 4, for example, we constructed an averaging filter window 4 times the original GSD. Thus, a 2 cm wide window was constructed for a row containing pixels of 0.5 cm GSD. Neighboring pixels that fell within the window bounds were averaged to obtain a single reflectance.

Results
Biomass across the three sites varied substantially from roughly 30 to 1000 g/m 2 , with the greatest range in the oldest marsh, and the lowest values in the 35-year-old marsh ( Figure 5). Biomass results are generally consistent with those observed previously in these marshes [42,63] and are representative of the high variability that occurs in fringing barrier island marshes over relatively short spatial scales (1-10 m). Figure 5 demonstrates a relatively good relationship between in situ measured above-ground biomass and estimates from the inversion of PROSAIL (the measured and estimated AGB of Spartina alterniflora: R 2 = 0.73; normalized root mean square error (NRMSE) = 23.40%). In Figure 6, we show examples of the retrieved biomass estimate of the PROSAIL model biomass for individual hyperspectral pixels in the scene. However, to compare against our ground truth quadrats of size 0.25 m 2 , we averaged the reflectance in our hyperspectral imagery over the region of the quadrat and compared the retrieved biomass from PROSAIL against that determined directly from the quadrat through laboratory analysis. This comparison provided the results shown in Figure 5. Where the data in our present analysis ( Figure 5) show greater scatter, noise is likely due to registration errors of the ground truth quadrat within the hyperspectral scene, mixed pixels (vegetation/soil) or soil pixels not fully removed by the NDVI spectral test (Figure 4), obscuration, and shadows. The scatterplot in Figure 5 does not show particular bias for one age fitting better than another to the regression line, but a more comprehensive analysis as a function of age will require samples from subsequent campaigns to fill in data gaps, for example, for the 2011 marsh. At the pixel level, as noted above, we applied the inversion methodology to hyperspectral images collected during our 2017 field survey. The inversion was applied to images from three different marsh ages: the marshes originating in 2011, 1974, and 1850. AGB maps from six selected images collected from our HSI system are shown in Figure 6. These six HSI images were collected simultaneously with the in situ measurements, which were used for the validation of the PROSAIL inversion methodology shown in Figure 5. The HSI images had a GSD ranging from 0.5 to 3 cm. The very high resolution scale hyperspectral imagery obtained in this study, at the sub-centimeter-to-centimeter scale, addresses a critical need in that the imagery is precisely at the scale where the variation in vegetation structure and plant biomass can best be characterized, that is, by looking directly at components of the individual plants and how these components vary across the landscape. The retrieved AGB from these images shows good correlation with past and present field surveys. For example, the AGB maps indicate higher marsh densities on the berm, while die-off areas show lower biomass. The 2011 marsh site generally had taller and denser Spartina alterniflora compared to the other two aged marsh sites, as indicated by the higher AGB observed in our inversion results in Figure 6. To evaluate how the biomass estimates shift as the spatial resolution changes, and make comparisons with the resolution that may be more characteristic of that obtained from aerial imagery, we downsampled the imagery to evaluate how resolution impacts the accuracy of assessment. We selected certain region of interests (ROIs) from the AGB maps in Figure 6 (the ROIs are shown in white) to observe the influence of spatial heterogeneity and resolution on the overall assessment of the biomass. These results are shown in the form of histograms in Figure 7. The first column in Figure 7 shows the histogram at the original spatial resolution, in the second column the spatial resolution was changed by a factor of four, and in column three by a factor of a hundred. The statistics for the ROIs are detailed in Table 2. This table and the histograms show that the standard deviation represents a large fraction of the mean biomass retrieved in the original image resolution. As the imagery spatial resolution is progressively degraded, the mean retrieved biomass can shift significantly, for example, in the 2011 marsh site for the ROI labeled "medium zone with die-off". As spatial resolution of the input imagery decreased, the trend was that these shifts in the mean appeared to depend on the relative density of the marsh site. The mean of the lower biomass sites did not change substantially, however, larger changes in the mean biomass occurred for the higher biomass regions. ; and (f) 2011 "medium zone with die-off". The first column shows the histogram at the original spatial resolution, in the second column the spatial resolution has been changed by a factor of four, and in column three by a factor of one hundred.  Figures 6 and 7. The spatial resolution "1" is the original resolution, "2" the resolution was changed by a factor of four, and "3" by a factor of hundred.

Discussion and Conclusions
We demonstrated a new approach to quantifying biomass in salt marsh systems using a novel mast-mounted hyperspectral imaging system [59]. This system provides pixels at sub-centimeter-scale spatial resolution, allowing the characterization of the marsh system at significantly higher spatial resolutions than typical airborne hyperspectral imaging systems. Imaging at this very fine spatial scale provides the opportunity to remotely retrieve and analyze small-scale details of coastal marsh ecosystems, which exhibit high spatial heterogeneity.
Using this fine-scale hyperspectral imagery at a range of salt marsh sites, which vary in biomass across an order of magnitude, as the input to a workflow that incorporates the PROSAIL radiative transfer model, we were able to successfully retrieve above-ground biomass at a resolution substantially greater than previously published work [22][23][24][25][26][27][28][29][30]. Our model output captured the high spatial heterogeneity of the system and retrieved AGB that correlated moderately well with ground-truth measurements (R 2 = 0.73). During PROSAIL model inversion, the relative accuracy of the biomass retrieval result relied heavily on using optimal wavelengths derived from EFAST. By incorporating the optimal wavelength ranges into the inversion methodology, we avoided common pitfalls typically observed with optimization techniques, such as the risks of local "false" minima in the parameter space. We saw significant improvement in our retrieval of AGB by performing the inversion within the wavelength ranges of 400-500 nm and 720-1000 nm, compared with using the complete wavelength range of our HSI system. The R 2 value improved by a factor greater than 2, increasing from 0.35 to 0.73 due to the incorporation of optimal wavelengths found by EFAST into the inversion methodology.
Inspection of retrieved biomass products showed that variations followed observed patterns in the field setting, with, for example, high biomass values obtained in the tall S. alterniflora near creekbank edges and lower biomass values found in die-off patches and in the short zone near the salt panne ( Figure 6). Although our HSI system was capable of collecting imagery at very high spatial resolution, we still have mixed pixels of vegetation, sediment and water in our data, which impacted the overall accuracy of our model in retrieving AGB. The accuracy of the retrieved biomass was also affected by the geo-registration method utilized for the mast-mounted HSI system. We conducted GPS surveys of fiducials within our scene to geo-register the images. The registration accuracy of the geo-referenced images were within approximately 1-10 cm, depending on the number of fiducials present within the scene. This affected our ability to precisely select ROIs where we collected ground-truth biomass data.
In this work, we also rescaled the retrieved PROSAIL salt marsh biomass to two different spatial scales (×4 and ×100) for selected regions within the marsh system in order to determine the applicability of the method across spatial scales. This rescaling revealed a loss of spatial detail ( Figure 7) and a shift in the mean biomass in the region of interest (Table 2). For low biomass regions (<500 g/m 2 ), the retrievals were similar across the three scales. However, for higher biomass, the variation across scales increased, with underestimates at coarser scales in the mid-range of biomass, and a 34% overestimate at the coarsest scale for the highest biomass region. This shift is indicative of the loss of accuracy that may occur when scaling up through a simple averaging approach that does not account for the detail found in the landscape at the natural scale of variation of the salt marsh system. This suggests that coarse-scale imagery, such as that obtained from aerial systems, may not provide the accuracy needed to assess biomass and thereby proportional carbon sequestration potential in systems with high spatial heterogeneity, such as salt marshes, and that additional work to verify our ability to use such imagery is required. Similar discrepancies in scaling across platforms have been found in a variety of other remote sensing applications (e.g., [93][94][95]). Additional work to compare imagery collected simultaneously at different scales may help to reconcile accuracy issues associated with the retrieval of biomass across platforms. Our scaling analysis illustrates that simple averaging may introduce some error in the process of retrieving biomass likely because a simple average does not account for the structure found at the very fine scales of our imagery. Indeed, to be successful in scaling, modeling will have to account more completely for the spatial relationships present in the plant structures and how they vary across the landscape. Our very fine scale imagery gives us access to this information, and long-term will be useful in improving models.
Biomass products such as those that we have obtained can support other types of analyses-for example, the development of models related to plant stress in the salt marsh landscape. Plant stress can directly influence biomass and may be an important factor in determining ecosystem services such as carbon storage. While there have been prior studies examining the effects of various stressors on S. alterniflora biomass [31,37,96], these models may be limited by the broad spatial scales of data acquisition when applied at the landscape level. Using our imaging system, we now have access to imagery at very high spatial resolutions ( sub-centimeter-to centimeter-scale). Spectral data have long been used to analyze the effects of stress on vegetation [96,97]; however, this enhanced spatial scale would allow us to observe the details of stress across the plant structure as opposed to over averages of multiple plants, which would normally be the scale observed in past remote sensing studies. The level of high spatial detail applied in this study may allow for better detection of spatially heterogeneous parameters such as salinity stress, oxidation-reduction potential, or leaf nitrogen content, which may be influencing biomass production. The level of high spatial detail applied in this study makes possible correlations with other plant biophysical measures as well as biophysical measures of the sediment surface, such as soil organic carbon, which are likely to be correlated with the above-ground biomass in salt marsh ecosystems [42]. In the long term, understanding these patterns of salt marsh stress and C sequestration are important for predicting marsh resilience. The authors are also grateful for support provided by NSF DEB Grants #1237733 #1832221 (VCR-LTER), which contributed to the fieldwork described in this paper.

Conflicts of Interest:
The authors declare no conflict of interest.