Impact of Atmospheric Correction on Classiﬁcation and Quantiﬁcation of Seagrass Density from WorldView-2 Imagery

: Mapping the seagrass distribution and density in the underwater landscape can improve global Blue Carbon estimates. However, atmospheric absorption and scattering introduce errors in space-based sensors’ retrieval of sea surface reﬂectance, affecting seagrass presence, density, and above-ground carbon (AGC seagrass ) estimates. This study assessed atmospheric correction’s impact on mapping seagrass using WorldView-2 satellite imagery from Saint Joseph Bay, Saint George Sound, and Keaton Beach in Florida, USA. Coincident in situ measurements of water-leaving radiance ( L w ), optical properties, and seagrass leaf area index (LAI) were collected. Seagrass classiﬁcation and the retrieval of LAI were compared after empirical line height (ELH) and dark-object subtraction (DOS) methods were used for atmospheric correction. DOS left residual brightness in the blue and green bands but had minimal impact on the seagrass classiﬁcation accuracy. However, the brighter reﬂectance values reduced LAI retrievals by up to 50% compared to ELH-corrected images and ground-based observations. This study offers a potential correction for LAI underestimation due to incomplete atmospheric correction, enhancing the retrieval of seagrass density and above-ground Blue Carbon from WorldView-2 imagery without in situ observations for accurate atmospheric interference correction.


Introduction
Seagrasses represent the greatest data gap in Blue Carbon habitat mapping, with their spatial extent largely under-represented in the Blue Carbon inventories for Canada, Mexico, and the United States [1].This under-representation is due to a combination of fundamental uncertainties about the true extent and density of seagrass meadows, as well as the logistical and labor challenges involved in the ground-based mapping of submerged seagrass populations, many of which are located in remote coastal areas [2].Remote sensing imagery obtained from air-and space-borne sensors can quantify seagrass and other benthic ecosystems, and the ability to remotely quantify seagrass meadows has been greatly enhanced by high-spatial-resolution (≤30 m) space-based sensors such as Landsat, IKONOS, Quickbird-2, RapidEye, WorldView-2 and -3, and Sentinel-2 [3][4][5][6][7][8][9][10][11][12].
Here, we explored the extent to which accuracy and precision in AC impact the ability to retrieve the seagrass distribution and density from WorldView-2 multispectral images from the Gulf Coast of Florida.Our primary goal was to develop a robust algorithm for seagrass classification and density determination that provides a path for the routine, accurate mapping of submerged aquatic vegetation (SAV) across the submarine landscape.The future ability to exploit the growing archive of images from orbiting sensors to quantify important changes in coastal environments requires us to assess the impacts of incomplete atmospheric correction on our classification and quantification algorithms.

Study Sites
We focused on three sites located along the Florida Gulf Coast, USA (Figure 1).Full characterization of the water column optical properties and quantification of seagrass abundance were undertaken at stations across all study sites (Figure 2) between 18 and 26 May 2010 near Keaton Beach, 1 and 9 November 2010 in Saint Joseph Bay, and 24 April and 5 May 2012 in Saint George Sound, coincident with the scheduled image acquisition at each site.
Here, we explored the extent to which accuracy and precision in AC impact the ability to retrieve the seagrass distribution and density from WorldView-2 multispectral images from the Gulf Coast of Florida.Our primary goal was to develop a robust algorithm for seagrass classification and density determination that provides a path for the routine, accurate mapping of submerged aquatic vegetation (SAV) across the submarine landscape.The future ability to exploit the growing archive of images from orbiting sensors to quantify important changes in coastal environments requires us to assess the impacts of incomplete atmospheric correction on our classification and quantification algorithms.

Study Sites
We focused on three sites located along the Florida Gulf Coast, USA (Figure 1).Full characterization of the water column optical properties and quantification of seagrass abundance were undertaken at stations across all study sites (Figure 2) between 18 and 26 May 2010 near Keaton Beach, 1 and 9 November 2010 in Saint Joseph Bay, and 24 April and 5 May 2012 in Saint George Sound, coincident with the scheduled image acquisition at each site.
Saint Joseph Bay is the only embayment in the eastern Gulf of Mexico not directly influenced by riverine input.It does, however, receive water rich in colored dissolved organic ma er (CDOM) from the Gulf County Canal, constructed in 1938, which links the Intracoastal Waterway to the bay [35].At Port Saint Joe, seagrass meadows are extensively distributed in shallow nearshore waters (<2 m depth) throughout the bay, especially at the southern (closed) end, and their distribution has been temporally stable from 1990 to 2020 [5].Saint Joseph Bay is the only embayment in the eastern Gulf of Mexico not directly influenced by riverine input.It does, however, receive water rich in colored dissolved organic matter (CDOM) from the Gulf County Canal, constructed in 1938, which links the Intracoastal Waterway to the bay [35].At Port Saint Joe, seagrass meadows are extensively distributed in shallow nearshore waters (<2 m depth) throughout the bay, especially at the southern (closed) end, and their distribution has been temporally stable from 1990 to 2020 [5].Saint George Sound forms the eastern portion of the Apalachicola Bay estuary/lagoon system, which is heavily influenced by the Apalachicola River, receiving discharge from agricultural, industrial, and municipal activities in the 50,000 km 2 Apalachicola-Flint-Cha ahoochee River Basin.Seagrasses are abundant in the shallow, quiescent waters along the mainland coast of the eastern portion of Saint George Sound between Carrabelle and St. Teresa and in the lee of Dog Island Reef (Figure 2B) that separates the southern boundary of St. George Sound from the northern Gulf of Mexico.Seagrasses are absent from the highly turbid waters in the western portion of Apalachicola Bay and in the much smaller Alligator Harbor at the east end.
Keaton Beach is located at the northern end of the Big Bend Seagrass Aquatic Preserve, which encompasses the second-largest contiguous area of seagrass habitat in the eastern Gulf of Mexico [36].Unlike St. Joseph's Bay and St. George Sound, in which seagrass meadows are confined to shallow embayments or protected lagoon waters, the seagrass meadows along Keaton Beach inhabit the open Gulf Coast to depths extending to the 3 m isobath.The coastal waters of Keaton Beach are influenced by CDOM-rich but oligotrophic discharge from the Steinhatchee River, a short river derived from the relatively undeveloped watershed dominated by cypress swamps, titi swamps, peat bogs [37], and other non-point sources along this section of the coast.

Measurement of Water Column Optical Properties and In Situ Rrs
All symbols and definitions described in this and following sections are summarized in Table 1.The inherent optical properties (IOPs) of spectral absorption (a, m −1 ) and beam a enuation (c, m −1 ) were measured at nine channels (412, 440, 488, 510, 532, 555, 650, 676, 715 nm) using an ac-9 Plus (Sea-Bird Instruments) deployed 1 m below the surface at each station.The ac-9 data were corrected for ambient temperature, salinity, and sca ering according to Method 3 [38] using a Conductivity, Temperature, and Depth (CTD) sensor incorporated into the ac-9 package.Sca ering correction was performed using data from a Hydroscat-6P series 300 (HOBI Labs, Inc) deployed alongside the ac-9.All data processing of the Hydroscat-6 data was undertaken in accordance with manufacturers' protocols, including the σ correction for highly sca ering waters.A enuation and absorption by particles plus CDOM (cpg, apg) were calculated by subtracting the absorption and Saint George Sound forms the eastern portion of the Apalachicola Bay estuary/lagoon system, which is heavily influenced by the Apalachicola River, receiving discharge from agricultural, industrial, and municipal activities in the 50,000 km 2 Apalachicola-Flint-Chattahoochee River Basin.Seagrasses are abundant in the shallow, quiescent waters along the mainland coast of the eastern portion of Saint George Sound between Carrabelle and St. Teresa and in the lee of Dog Island Reef (Figure 2B) that separates the southern boundary of St. George Sound from the northern Gulf of Mexico.Seagrasses are absent from the highly turbid waters in the western portion of Apalachicola Bay and in the much smaller Alligator Harbor at the east end.
Keaton Beach is located at the northern end of the Big Bend Seagrass Aquatic Preserve, which encompasses the second-largest contiguous area of seagrass habitat in the eastern Gulf of Mexico [36].Unlike St. Joseph's Bay and St. George Sound, in which seagrass meadows are confined to shallow embayments or protected lagoon waters, the seagrass meadows along Keaton Beach inhabit the open Gulf Coast to depths extending to the 3 m isobath.The coastal waters of Keaton Beach are influenced by CDOM-rich but oligotrophic discharge from the Steinhatchee River, a short river derived from the relatively undeveloped watershed dominated by cypress swamps, titi swamps, peat bogs [37], and other non-point sources along this section of the coast.

Measurement of Water Column Optical Properties and In Situ R rs
All symbols and definitions described in this and following sections are summarized in Table 1.The inherent optical properties (IOPs) of spectral absorption (a, m −1 ) and beam attenuation (c, m −1 ) were measured at nine channels (412, 440, 488, 510, 532, 555, 650, 676, 715 nm) using an ac-9 Plus (Sea-Bird Instruments) deployed 1 m below the surface at each station.The ac-9 data were corrected for ambient temperature, salinity, and scattering according to Method 3 [38] using a Conductivity, Temperature, and Depth (CTD) sensor incorporated into the ac-9 package.Scattering correction was performed using data from a Hydroscat-6P series 300 (HOBI Labs, Inc.) deployed alongside the ac-9.All data processing of the Hydroscat-6 data was undertaken in accordance with manufacturers' protocols, including the σ correction for highly scattering waters.Attenuation and absorption by particles plus CDOM (c pg , a pg ) were calculated by subtracting the absorption and scattering of pure water [39] from the temperature-and salinity-corrected data.The total spectral scattering coefficient (b p (λ),m −1 ) was calculated as: under the assumption that CDOM did not contribute to scattering.Remote sensing reflectance from ELH atmospheric correction sr A HyperPro profiling radiometer (SatLantic Instruments) and hyperspectral tethered spectroradiometer buoy (HTSRB, SatLantic Instruments) were floated in tandem at each station to measure downwelling spectral irradiance (E s (λ,0 + )) above the sea surface and upwelling spectral radiances (L u (λ)) at 0.21 m and 0.65 m, respectively, beneath the sea surface.The diffuse attenuation coefficient for upwelling radiance was calculated as: where ∆z was the difference in depth (0.44 m) between the sensors placed at 0.21 and 0.65 m.The water-leaving radiance (L w (λ)) was calculated by propagating L u (λ) to the surface using K Lu (λ) and transmitting it across the air-water interface: where τ u (=0.54,) represented the fraction of upwelling radiance transmitted upward across the water-air interface of a flat ocean [40,41].In situ remotely sensed reflectances (R rs (λ), sr −1 ) were then calculated as: R rs (λ) was then averaged to match the multispectral wavelength bins of the WorldView-2 sensor, using the response functions provided by DigitalGlobe (Table 2).

In Situ Seagrass Density
Direct counts of Thalassia testudinum (turtle grass), Halodule wrightii (shoal grass), and Syringodium filiforme (manatee grass) shoots were performed by SCUBA divers using 20 to 30 quadrats (0.04 m 2 ) randomly located within a 20 m radius at optically shallow water stations at each site (Figure 2).Halophila spp.(paddlegrass) and Ruppia maritima (widgeon grass) were not present in any of the quadrats.One shoot of each species present was collected from each quadrat for subsequent measurement of total leaf length, leaf width, and number of leaves per shoot and determination of leaf area index (LAI = shoots m 2 seafloor × m 2 leaf area shoot ).The LAI of Syringodium filiforme was calculated using the diameter of the cylindrical leaves for leaf width.

Top-of-Canopy Reflectance Measurements
Downwelling (E d (λ,z)) and upwelling (E u (λ,z)) spectral irradiances were measured at the top of the seagrass canopy and at 15 cm above the bare sediment at ~10 randomly selected spots at each station using our Diver-Operated Benthic Bio-Optical Spectrometer (DOBBS), a three-channel radiometer (HR-3, HOBI Labs, Inc.) configured for portable in situ operation by a diver.Instrument spectra (nominally 0.3 nm resolution) were interpolated to 1 nm with a cubic spline and smoothed using a 21 nm running average.Spectral reflectance (R b (λ), dimensionless) of the seagrass canopy or bare substrate was calculated as (λ,z) .R b (green) was then calculated by averaging R b (λ) across the WorldView-2 green band (511 to 581 nm).The mean and standard error of R b (green) were calculated for each station.
The top-of-canopy reflectance measurements were paired with coincident LAI measurements (Section 2.3) to generate a regression model that could be used with top-ofcanopy reflectance retrieved from R rs (Section 2.6.3) to estimate LAI from the images.

Satellite Data
The WorldView-2 satellite imaging system operated by Maxar Technologies, Inc. (formerly DigitalGlobe) was tasked to obtain multispectral images of Keaton Beach, Saint Joseph Bay, and Saint George Sound.Tasking information included the optimal geometry relative to the target to minimize direct solar reflectance (sunglint) and local cloud cover.Image collection took place on 20 May 2010 (Keaton Beach), 14 November 2010 (Saint Joseph Bay), and 27 April 2012 (Saint George Sound, Table 3 and Figure 2).Ground resolution, which depended on the view angle required to minimize glint, varied from 2 m at Saint Joseph Bay to 3 m at both Keaton Beach and Saint George Sound (Table 3).Raw, top-of-atmosphere digital values (Level 1b data product) were converted to sea surface reflectance (R rs ) for each image using two atmospheric correction methods (Figure 3).The ELH method incorporates both radiometric calibration and atmospheric correction into a single step [27,42] using in situ observations of R rs collected coincidently with the satellite overpass (Section 2.2) at "optically deep" stations where the geometric water depth exceeded the depth of the 10% isolume, defined by the diffuse attenuation coefficient K d .These in situ reflectance values were regressed against the average digital value obtained from a 4 × 4 pixel array extracted from the image corresponding to the location of each station to create R rs gains (slopes) and offsets (intercepts) for all bands that were subsequently applied to the entire image.Accuracy of the in situ L w (λ) measurements was validated against HydroLight [43] simulations parameterized with IOP data from the ac-9 and Hydroscat-6 measurements at those stations (Section 2.2).R rs generated from this method will be referred to in the text as R rs(ELH) (Figure 3).Raw, top-of-atmosphere digital values (Level 1b data product) were converted to sea surface reflectance (Rrs) for each image using two atmospheric correction methods (Figure 3).The ELH method incorporates both radiometric calibration and atmospheric correction into a single step [27,42] using in situ observations of Rrs collected coincidently with the satellite overpass (Section 2.2) at "optically deep" stations where the geometric water depth exceeded the depth of the 10% isolume, defined by the diffuse a enuation coefficient Kd.These in situ reflectance values were regressed against the average digital value obtained from a 4 × 4 pixel array extracted from the image corresponding to the location of each station to create Rrs gains (slopes) and offsets (intercepts) for all bands that were subsequently applied to the entire image.Accuracy of the in situ Lw(λ) measurements was validated against HydroLight [43] simulations parameterized with IOP data from the ac-9 and Hydroscat-6 measurements at those stations (Section 2.2).Rrs generated from this method will be referred to in the text as Rrs(ELH) (Figure 3).Atmospheric correction was also undertaken using the DOS approach [44].Radiometric calibration and DOS were performed on each image in separate steps, as described in Coffer et al. [3], using radiometric calibration coefficients from Kuester [45].The red-edge band (MS5, Table 2) was used as the DOS reference band because radiometric calibration of the near-infrared (NIR) bands using coefficients provided by Maxar frequently yielded negative radiances over dark water targets before atmospheric correction, preventing their use for DOS [3].Inland water pixels were excluded from the red-edge analysis by masking the land area using a shapefile of the coastline [46].A normalized difference water index (NDWI) was then applied to remove any remaining shoreline pixels [47].The median of the red-edge R rs values was then computed from the optically darkest 5% of the remaining water pixels.We assumed that half of the red-edge radiance emanating from the dark pixels represented atmospheric contamination [3], while the remainder represented true water-leaving radiances.The scattering factor (s) was computed as σ = λ 4.75  red edge × α red edge (5) where α is anchored to the correction at the red edge at 0.5.The wavelength-dependent Rayleigh exponent was set to 4.75 [3].The scattering factor was then applied to the center wavelength of each band to compute the correction factor (η band ), which was subtracted from the band-specific R rs for all pixels in the scene:

Date
The radiometrically calibrated and atmospherically corrected R rs values generated from this method are identified below as R rs(DOS) (Figure 3).

Machine Learning Supervised Classification
A deep convolutional neural network (DCNN, [17]) previously used to separate seagrass from other targets in imagery from WorldView-2 and RapidEye [3] sensors was used to classify the R rs(ELH) and R rs(DOS) images.Four classes were identified at all sites: (i) optically deep water, in which the bottom type was unidentifiable, (ii) optically shallow submerged sand, (iii) submerged seagrass, and (iv) dry land consisting of natural vegetation and human development.A fifth class of emergent intertidal sand covered with benthic algal film was identified only in Saint Joseph Bay.The DCNN was trained using representative patches (rectangular polygons) of each class, referred to as regions of interest (ROIs), identified in each image through a combination of local knowledge and visual confirmation (Figure 2).ROIs were split into training and testing subsets for the DCNN model, and a subsequent patch size of 3 × 3 pixels was extracted from the ROIs to use as training samples.The DCNN consisted of a sequential model, or linear stack, of six layers.The first layer was a convolutional layer consisting of 32 filters with a kernel size of 1 × 1 × 8 (containing 1 pixel and 8 channels), followed by a rectified linear activation function.Next, a dropout layer randomly set 1% of the outputs from the first layer to zero.The second convolutional layer consisted of 16 filters with a kernel size of 3 × 3 × 32 pixels, again, followed by a rectified linear activation function.Then, another 1% dropout layer was added before flattening the model.Finally, a dense layer with a SoftMax activation function was used to compute pixel-based probabilities for each class [48].This DCNN model was trained for 500 epochs, where an epoch represented the point at which all training data have been processed one time, using a batch size of 256.Once trained, the model was applied across the entire image to classify all pixels, resulting in separately classified images for R rs(ELH) and R rs(DOS) values.

Calculation of Seagrass Density
R rs(ELH) and R rs(DOS) values were used to generate benthic reflectance at the top of the seagrass canopy in the green band (R b (green)) [22,23]: where R rs (green) represents either R rs(ELH) or R rs(DOS) values from band MS3 (green, centered at 545 nm).Water depth (z b ) was extracted from NOAA's Continuously Updated Digital Elevation Model [49] at 1 m horizontal resolution and resampled to match the grid dimensions and spatial resolution of the satellite imagery (2 to 3 m) using ArcGIS Pro.The resulting digital elevation models (DEM) for all sites were then used along with local tide predictions (Table 1) to calculate water depth across each scene at the times of WorldView-2 image acquisition.Downwelling diffuse attenuation (K d ) and depth-averaged upwelling radiance attenuation −K lu (z b ) coefficients were modeled using Hydrolight (© C. Mob- ley, Sequoia Scientific) parameterized with inherent optical properties (IOPs) measured in situ.Q b represents the ratio of E u to L u at the base of the water column produced by the reflectance of E d from the sea floor, which was assumed to be π for a Lambertian bottom boundary.The relative transmittance of L u through the water-air interface (τ u ) was approximated as 0.54 [40].Leaf area index was determined by linear regression of R b (green) against direct shoot counts made at in situ stations distributed across each scene.AGC seagrass was estimated from LAI using a series of transfer coefficients for fresh weight (Equation ( 8)), dry weight (Equation ( 9) [25,26]), and organic carbon (Equation (10) [24]).
Above-ground Fresh Wt (g FW) = 500 Above-ground Dry Wt (g DW) = fresh wt × 0.2 (g DW g −1 FW) Above-ground seagrass carbon stock (g C) = dry wt (g DW) × 0.35 (g C g −1 DW) ( 10) R rs(ELH) and R rs(DOS) values at 400 randomly generated points were selected across each study site for comparison.Points were generated using the "create random points" tool in ArcGIS Pro across all marine targets to collect a range of R rs values.R rs(ELH) and R rs(DOS) for each band were extracted from these points and subjected to linear regression.Differences in the band-specific relationships between R rs(ELH) and R rs(DOS) among the three sites were evaluated using Analysis of Covariance (ANCOVA).

Comparison of In Situ and Retrieved LAI
Direct comparisons of in situ and satellite-retrieved LAI are challenging due to the small-scale spatial heterogeneity found in seagrass meadows, where quadrat-scale LAI measured using a 0.04 m 2 quadrat can range from 0 to 4 m 2 m −2 within the 3 m footprint of the WV-2 pixel.In order to determine whether in situ and satellite-derived LAI values were statistically similar, we extracted LAI from image pixels within a 20 m radius of our in situ seagrass stations using the "create buffer tool" in ArcGIS Pro (Figure 2).Box-andwhisker plots were generated to illustrate similarities in mean, median, and data ranges.Statistical similarity between in situ LAI calculated from direct quadrat counts (Section 2.3) and remotely detected LAI was determined via a non-parametric Kruskal-Wallis test [50] followed by a pairwise Mann-Whitney U-test [51], run using the SciPy statistical package [52].

Water Column Optical Properties
The water column at all three sites consisted of optically complex Case 2 waters, with high absorption and/or scattering dominated by a combination of non-algal particles and CDOM.Keaton Beach had the highest a pg ; values at 412 and 443 nm were twice those observed at the other sites (Figure 4A).The comparatively low contribution of b p to c pg (Figure 4B,C) at Keaton Beach indicates that the a pg was dominated by CDOM and not high levels of phytoplankton or other particles.Saint Joseph Bay and Saint George Sound had similar values of a pg (Figure 4A).However, the b p in Saint George Sound was three times higher than that observed at Saint Joseph Bay and accounted for 75% of the c pg , suggesting a high load of non-algal particles in this lagoon fed by the Apalachicola River (Figure 5B).There was a high degree of variation in the b p across Saint George Sound, with lower b p and c pg observed offshore near Dog Island Reef.In Saint Joseph Bay, both a pg and b p contributed equally to light attenuation.Based on the optical properties, the depth of the 10% isolume, which we used to define the limit for optically shallow water and which often corresponds to the lower limit of seagrass colonization [53,54], was 0.9 m in Saint George Sound, 2.1 m at Keaton Beach, and 2.3 m in Saint Joseph Bay.

Saint Joseph Bay
Despite differences observed in Rrs derived by the two atmospheric correction methods, there was broad agreement in the seagrass classification, with extensive continuous seagrass beds identified covering the shallow areas of less than 2 m depth along the southern edge of the bay (Figure 6A,B).The total seagrass area was classified as 25.2 km 2 in the ELH image and 27.0 km 2 in the DOS image (Table 6), with an area of 23 km 2 commonly

Saint Joseph Bay
Despite differences observed in Rrs derived by the two atmospheric correction methods, there was broad agreement in the seagrass classification, with extensive continuous seagrass beds identified covering the shallow areas of less than 2 m depth along the southern edge of the bay (Figure 6A,B).The total seagrass area was classified as 25.2 km 2 in the ELH image and 27.0 km 2 in the DOS image (Table 6), with an area of 23 km 2 commonly

Comparison of Remote Sensing Reflectance between Atmospheric Correction Methods
R rs(DOS) was higher than R rs(ELH) across all bands at all three sites, especially in bands MS7 (425 nm) and MS4 (480 nm), where R rs(DOS) was an order of magnitude higher than R rs(ELH) (Figure 5A,B).The relationship between R rs(DOS) and R rs(ELH) was positive and linear in all bands and sites, except for Keaton Beach, where we found a negative, albeit linear, relationship between R rs(ELH) and R rs(DOS) (Figure 5A 425 nm, Table 4).Overall, the Keaton Beach image was darker than both Saint Joseph's Bay and Saint George Sound images in bands MS7, 4, and 3 in both ELH and DOS corrections (Figure 5A-C).R rs(DOS) was more divergent from R rs(ELH) (slopes > 1) in the yellow (MS6) and red (MS2) bands (Table 4 and Figure 5D,E) at Keaton Beach and Saint George Sound relative to St. Josephs Bay, where R rs(DOS) converged with R rs(ELH) at higher values.The NIR bands MS1 and 8 (not shown) had an insufficient signal-to-noise ratio to be useful, regardless of the atmospheric correction method.The regression slopes between R rs(DOS) and R rs(ELH) were statistically different between each of the three sites for all bands (ANCOVA p < 0.001, Table 5).In the green band (MS3, 545 nm), the relationship between R rs(DOS) and R rs(ELH) appears to converge between all the sites, with slopes that were close but still did not overlap within the 95% confidence limits, rendering them significantly different (Figure 5C and Table 5).Despite differences observed in R rs derived by the two atmospheric correction methods, there was broad agreement in the seagrass classification, with extensive continuous seagrass beds identified covering the shallow areas of less than 2 m depth along the southern edge of the bay (Figure 6A,B).The total seagrass area was classified as 25.2 km 2 in the ELH image and 27.0 km 2 in the DOS image (Table 6), with an area of 23 km 2 commonly identified as seagrass in both images (Figure 6C).Differences in seagrass classification largely occurred along the transition zones between submerged seagrass meadows and the intertidal zone, which was dominated by filamentous red and brown algae, as well as along the deeper edge of the meadows, where sparse vegetation was transitioning to bare sand (Figure 6C).The greater seagrass area recovered from the DOS image was located in the intertidal and in shallow inlets that were classified as intertidal in the ELH image (Figure 6C).In situ observations identified these areas as benthic algae and not seagrass.The ELH image also retrieved benthic information at greater depths, resolving 42% more optically shallow sand area along the 2 m isobath that was classified as optically deep water by the DOS image (Figure 6).

Keaton Beach
Seagrass was distributed at Keaton Beach in a more-or-less continuous meadow from the shoreline out to optically deep water (~2 m isobath) (Figure 8).The total cover

Saint George Sound
In Saint George Sound, seagrass was identified along the shallow areas (<2 m depth) close to the mainland shore and along the lee side of Dog Island Reef that forms the seaward barrier between the Sound and the Gulf of Mexico (Figure 7A,B).A total of 17.9 km 2 of seagrass was classified from the ELH image, and 17.3 km 2 was identified from the DOS image (Table 6).Although the total area identified as seagrass differed by only 3% between the images, there was 3 km 2 of seagrass in the ELH-corrected image in a large patch to the west of Alligator Harbor (Figure 7C, indicated by the circle) and a patch in deeper water along the coastal margin between Lanark Reef and Turkey Point that was not resolved in the DOS-corrected image.The presence of seagrass in these disputed areas was verified from in situ diver surveys conducted by the authors in 2017.Visual in situ surveys also confirmed the sandbank south of Turkey Point to consist of benthic algae, reef invertebrates, and shelly aggregates that were classified as submerged sand in the ELH image but incorrectly classified as seagrass in the DOS image (Figure 7C, indicated by the rectangle).Twice as many pixels of optically shallow sand were identified in the ELH image (16.7 km 2 ) compared to the DOS image (7.8 km 2 ), which classified these disputed pixels as optically deep water.In particular, the ELH image resolved more of the optically shallow sand around Dog Island Reef and the optically shallow sand between the shore and Dog Island Reef.These discrepancies have no in situ verification but visually appear to be dominated by submerged sand, as classified in the ELH image.

Keaton Beach
Seagrass was distributed at Keaton Beach in a more-or-less continuous meadow from the shoreline out to optically deep water (~2 m isobath) (Figure 8).The total cover

Keaton Beach
Seagrass was distributed at Keaton Beach in a more-or-less continuous meadow from the shoreline out to optically deep water (~2 m isobath) (Figure 8).The total cover encompassed 70.5 km 2 in the ELH image and 72.3 km 2 in the DOS image, equating to a 2% difference between calibration/atmospheric correction methods (Table 6).A narrow strip of seagrass of approximately 3.2 km 2 identified along the deeper edge of the seagrass beds in the ELH image and positively identified as sparse seagrass by diver surveys (Figure 8A,C) was classified as optically deep water in the DOS image (Figure 8B).Conversely, 5.1 km 2 of seagrass along the shallow intertidal shoreline was classified as seagrass in the DOS image and as land in the ELH image, including Hagens Cove, which is an optically shallow area with some seagrass, as seen in the DOS image (Figure 8C).There may be more seagrass at the deeper edge of identified meadows on the south end of the image; however, the verification of seagrass in water with high CDOM loads required in situ observations that we did not make.The strong negative relationship between in situ log-transformed R b (green) and LAI for all stations across the three study sites yielded the predictive relationship shown in Figure 9 and expressed in Equation ( 11) (R 2 = 0.75, slope SE = 0.31, intercept SE = 0.39).The negative log-transformed R b values in the green band represent lower (darker) reflectance at the top of the canopy and are correlated with a higher seagrass density (LAI).This relationship provided the mathematical basis for predicting LAI on a pixel-by-pixel basis from R b(green) retrieved from the WorldView-2 imagery using Equation (7).LAI = −2.98logR b(green) − 2.17 (11) It should be noted when using this model that, once densities increase past 3 m 2 leaf m −2 ground, the relationship is expected to no longer be linear, with little additional reduction in canopy reflectance in the green occurring once LAI becomes dense enough to occlude any reflectance from the brighter sand at the base of the vegetation.

Saint Joseph Bay
LAI was spatially variable throughout Saint Joseph Bay (Figure 10A).The densest areas identified in the ELH image were in the shallow waters at the southern end of the bay, reaching 2.6 m 2 leaf m −2 ground.The median LAI of the whole bay was 1.89 m 2 leaf m −2 ground (Table 7).The seagrass density decreased sharply as the deeper edges of the meadows transitioned to bare sand.AGC seagrass in Saint Joseph Bay estimated from LAI (ELH) was 1.66 Gg, with an area-specific carbon density of 66 g m −2 (Table 7).Although the spatial pattern of seagrass distribution retrieved from R rs(ELH) was similar to that produced using R rs(DOS) (Figure 10A,B), LAI (DOS) was approximately 50% lower than LAI (ELH) across the entire image (median LAI 0.98 m 2 leaf m −2 ground, Table 7 and Figure 10B).Correspondingly, the underestimation of LAI reduced the bay-wide AGC seagrass by 47% to 0.88 Gg (Table 7).LAIs obtained through in situ measurements (LAI (IN SITU) ) were not significantly different from the retrieved LAI (ELH) from pixels within our station locations (Table 8 and Figure 11A).There were, however, notable differences in the range of values between these two datasets, with the LAI (IN SI) exhibiting a substantially broader span in comparison to the values retrieved via satellite.The distribution of LAI (ELH) values exclusively fell within the interquartile range of the in situ dataset but failed to capture LAI values surpassing 2.5 m 2 m −2 , densities that were observed in situ.LAI (DOS) values were significantly lower than observed in situ (Table 8), with the majority of LAI (DOS) values falling below the lower quartile of the LAI (IN SITU) and entirely below the lowest LAI (ELH) values (Figure 11A).3.4.3.Saint George Sound LAI (ELH) values were high along the mainland shore of Saint George Sound behind the complex of shallow sandbars created by Lanark Reef and Turkey Point Reef and along the leeward edge of Dog Island Reef that separates the Sound from the Gulf of Mexico (Figure 12A).The median LAI (ELH) across the scene was 1.96 m 2 leaf m −2 seafloor, and the highest LAI (ELH) of ~2.6 m 2 leaf m −2 seafloor was found in the relatively turbid shallow water (~0.5 m deep) close to shore (Figure 12A).Similar densities (~2.6 m 2 leaf m −2 ground) were found in the meadows occupying the deeper (~2 m) but clearer waters along the lee of Dog Island Reef.The total AGC seagrass derived from LAI (ELH) was 1.24 Gg, with an area-specific carbon of 68 g m −2 (Table 7).As with Saint Joseph Bay, LAI (DOS) values were ~50% lower than LAI (ELH) (Figure 12B), reducing the estimation of total AGC seagrass to 0.62 Gg and the area-specific carbon to 36g m −2 (Table 7).from pixels collected from a 20 m radius around the in situ station locations (Figure 1).(A).Saint Joseph Bay, (B).Saint George Sound, (C).Keaton Beach.

Saint George Sound
LAI(ELH) values were high along the mainland shore of Saint George Sound behind the complex of shallow sandbars created by Lanark Reef and Turkey Point Reef and along the leeward edge of Dog Island Reef that separates the Sound from the Gulf of Mexico (Figure 12A).The median LAI(ELH) across the scene was 1.96 m 2 leaf m −2 seafloor, and the highest LAI(ELH) of ~2.6 m 2 leaf m −2 seafloor was found in the relatively turbid shallow water (~0.5 m deep) close to shore (Figure 12A).Similar densities (~2.6 m 2 leaf m −2 ground) were found in the meadows occupying the deeper (~2 m) but clearer waters along the lee of Dog Island Reef.The total AGCseagrass derived from LAI(ELH) was 1.24 Gg, with an areaspecific carbon of 68 g m −2 (Table 7).As with Saint Joseph Bay, LAI(DOS) values were ~50% lower than LAI(ELH) (Figure 12B), reducing the estimation of total AGCseagrass to 0.62 Gg and the area-specific carbon to 36g m −2 (Table 7).While the median values of LAI(IN SITU) and LAI(ELH) appear to be similar (Figure 11B), the two datasets were found to be significantly different (Table 8) due to the small dynamic range of LAI(ELH).The middle 50% of LAI(ELH) encompassed values of 1.98 to 2.10 m 2 leaf m −2 ground (Figure 11B), in contrast to the larger range observed in LAI(IN SITU), from 1.7 to While the median values of LAI (IN SITU) and LAI (ELH) appear to be similar (Figure 11B), the two datasets were found to be significantly different (Table 8) due to the small dynamic range of LAI (ELH) .The middle 50% of LAI (ELH) encompassed values of 1.98 to 2.10 m 2 leaf m −2 ground (Figure 11B), in contrast to the larger range observed in LAI (IN SITU), from 1.7 to 2.8 m 2 leaf m −2 ground.With such a narrow range of distribution of LAI (ELH), the higher and lower range of values measured in situ was not resolved.As in Saint Joseph Bay, both the mean and median of LAI (DOS) fell below the lower quartile of the in situ dataset and entirely below the lowest LAI (ELH) retrieval (Figure 11B), making these retrievals significantly different from both LAI (ELH) and LAI (IN SITU) (Table 8) .

Keaton Beach
The seagrass at Keaton Beach occupied a 3 km wide band of homogeneous density ~3 m 2 leaf m −2 ground parallel to the shore (Figure 13A).The 2 m isobath marked the edge of the continuous seagrass meadow, with densities dropping sharply to bare sand, transitioning to optically deep water.The median LAI (ELH) was the highest of all three regions at 3.33 m 2 leaf m −2 ground (Figure 13A and Table 7).Both the higher total area and densities in Keaton Beach were reflected in the highest AGC seagrass and area-specific carbon of all three study sites at 8.34 Gg, with an area-specific carbon of 118 g m 2 (Table 7).As with the other two sites, LAI (DOS) was lower than LAI (ELH) (Figure 13B), which translated into a 58% underestimation of AGC seagrass calculated from LAI (DOS) compared to that from LAI (ELH) (Table 7).When comparing co-located LAI(IN SITU) and LAI(ELH), the satellite retrievals were found to be significantly higher (Table 8), exceeding the upper quartile of the in situ data while remaining within the range of the maximum values (Figure 11C).Although the mean LAI(ELH) was greater, the retrievals effectively capture the elevated in situ densities observed at this site when compared to Saint Joseph Bay and Saint George Sound.LAI(DOS) was significantly lower than LAI(IN SITU) (Table 8), exhibited a significantly narrower range of values when compared to LAI(ELH), and had the greatest divergence from LAI(ELH) of all three sites (Figure 11C).When comparing co-located LAI (IN SITU) and LAI (ELH ), the satellite retrievals were found to be significantly higher (Table 8), exceeding the upper quartile of the in situ data while remaining within the range of the maximum values (Figure 11C).Although the mean LAI (ELH ) was greater, the retrievals effectively capture the elevated in situ densities observed at this site when compared to Saint Joseph Bay and Saint George Sound.LAI (DOS) was significantly lower than LAI (IN SITU) (Table 8), exhibited a significantly narrower range of values when compared to LAI (ELH) , and had the greatest divergence from LAI (ELH) of all three sites (Figure 11C).The underestimation of LAI from R rs(DOS) relative to both R rs(ELH) was a consistent pattern across all three sites (Figure 14).The highly predictive relationship between LAI (ELH) and LAI (DOS) (r 2 0.94) can be approximated across the entire dataset with an exponential function (Figure 14A and Table 9).When analyzed separately, the relationships at both Saint Joseph Bay and Saint George Sound were best represented by linear functions (Figure 14B,C and Table 9), while the Keaton Beach fit remained slightly exponential (Figure 14D and Table 9).
in RMSE.Using either the single exponential function that encompassed all the sites or the individual linear/exponential models specific to each site, the accuracy of correcting LAI (DOS) to values equivalent to LAI (ELH) was high (Table 10).Using algorithm 1, the mean absolute error (MAE) at all three sites was less than 0.2 m 2 leaf m −2 ground, which is approximately 10% of the mean LAI values for Saint Joseph Bay and Saint George and 6% for Keaton Beach.The root-mean-square error (RMSE) was similar in magnitude to the MAE (Table 10), with the greatest errors occurring for LAI predicted at Saint George Sound.Using algorithms tuned for individual study sites improved both the MAE and RMSE for Saint Joseph Bay and Saint George Sound by approximately 50%, reducing the error in the predicted LAI to less than 5% of the mean values (Table 10).For Keaton Beach, an individual exponential fit increased MAE over algorithm 1 and produced only a small improvement in RMSE.

Discussion
This study revealed that the use of published radiometric calibration factors combined with DOS atmospheric correction produced images with much higher reflectances than those produced by empirical calibration of the raw satellite images against direct ground observations.However, the residual brightness, which we assume to be primarily atmospheric contamination but could include radiometric calibration inaccuracy, did not dramatically affect the supervised classification of seagrass presence in these optically complex coastal environments.Classification differences between the ELH and DOS methods were most apparent at the transition zones defining the shallow and deep edges of the meadows, but these transition zones represented a relatively small fraction of the total area.Thus, the choice of atmospheric correction method may not be critical for detecting the presence of seagrass beds, but it may be important for accurately identifying the boundaries of the meadows.Furthermore, the residual brightness in the green band resulting from the DOS approach caused the retrievals of seagrass LAI, and therefore, the carbon biomass, to be underestimated by approximately 50% compared to the ELH-processed image, which showed a stronger concordance with in situ measurements.Nonetheless, our comparison provides a potential path for rectifying DOS-based LAI retrievals that may permit seagrass LAI and above-ground carbon to be estimated from the existing collection of WorldView-2 images that lack in situ observations needed to accurately correct atmospheric interference in the green region of the spectrum.
Residual atmospheric contamination leading to higher radiances, especially in the blue and green bands, relative to in situ observations has been observed in both WorldView-2 and other high-spatial-resolution sensors, such as Pleiades, Sentinel-2, and Landsat-5 and 7 [30,33,34], and may be a widespread problem in atmospheric correction methods generally (B.Collister, NASA, pers comm).Coffer et al. [3] compared WorldView-2 R rs after DOS correction to in situ observations and concluded that it may not be just an atmospheric correction problem, as uncalibrated drift in post-launch sensor responsivity may also play a role.
Despite the presence of high residual radiances in the blue and green bands (bands 1 and 2), the classification of the DOS-processed image successfully separated seagrass pixels from optically shallow sand and optically deep water in three distinct WorldView-2 images that represented a range of optically complex waters characteristic of coastal Florida.This indicates that the WorldView-2 archive can be used to produce reliable supervised classification maps of seagrass distributions, even when atmospheric correction is imperfect.
The differences in the seagrass area between ELH-and DOS-corrected images were found to account for less than 8% of the total coverage, primarily resulting from the overclassification of (i) optically shallow bare sand along the deep margins and (ii) intertidal environments vegetated by benthic algae as seagrass from R rs(DOS) .These misclassifications were caused not only by the overestimation of R rs(DOS) but also by the narrow dynamic range seen in the blue band, which limits the ability to resolve targets with small spectral differences, such as the shallow edges of the seagrass meadows and green intertidal environments.The ELH method better replicated the true spectral shape of sea surface R rs , preserving the larger dynamic range in all bands, but especially in the blue band, which improved the resolution of some boundaries between the deeper seagrass meadows and bare sand that were spectrally indistinguishable in the DOS image.
A few estimates of the seagrass area for these sites are available for comparison.In September 2006, airborne hyperspectral imagery was collected at Saint Joseph Bay, from which seagrass was estimated to cover 18.3 km 2 [23].The flight lines did not cover the northernmost section of the bay, which may account for the slightly lower estimate.A 30-year analysis of seagrass populations in Saint Joseph Bay between 1990 and 2021 using Landsat 5 through 8 imagery returned areal estimates varying from 17.5 to 27 km 2 , with no discernable temporal trend [5], placing our 2010 estimate at the higher end of seagrass coverage estimates.Variations in the total seagrass coverage in Saint Joseph Bay have been linked to hurricanes, which can cause short-term (6-month) declines in seagrass [5].
Importantly, our DCNN analysis of this multispectral imagery generated a classification map of Saint Joseph Bay that was nearly identical to that achieved with hyperspectral imagery collected over the same site in 2006 [23] that included an intertidal class dominated by filamentous red algae.However, Hill et al. (2014) also demonstrated how the delineation of the seagrass habitat deteriorated when the spatial resolution was coarsened beyond 10 m.Thus, the ability to generate similar classification results with multispectral and hyperspectral sensors suggests that the greatest need in categorical classification coastal remote sensing may be spatial rather than spectral resolution.

Seagrass Density Quantification
Although the ELH imagery produced seagrass densities that were comparable to in situ observations, direct comparisons to the latter present challenges due to the small-scale spatial heterogeneity found in seagrass meadows.This heterogeneity results in larger variances in shoot counts measured by small quadrats (e.g., 0.04 m 2 ) within a 20 m radius than in corresponding satellite pixels at the same locations.At our sites, in situ seagrass densities exhibited high variability on small scales, ranging from sparse to dense within a few meters.This variability is reflected in the large range of LAI observed in the in situ datasets.In contrast, the coarser spatial resolution of the WorldView-2 sensor (4 to 9 m 2 ) resulted in lower variability in retrieved LAI due to the integration of the bottom reflectance signal over the larger footprint of each pixel.Despite the limitations in assessing match-ups between remotely sensed and in situ seagrass densities, LAI (ELH) in Saint Joseph Bay and Saint George Sound from pixels within our in situ station locations did show similarity to the mean densities measured in situ.At Keaton Beach, the mean LAI (ELH) from our station locations was double the mean from the collocated in situ observations.However, high seagrass densities were prevalent in the in situ counts conducted at Keaton Beach, as seen in the presented box plots, which show a skew toward higher LAI, with several high-density quadrats more than 5 m 2 leaf m −2 ground.LAI (ELH) retrievals at Keaton Beach did resolve lower densities closer to the coastline, which is seen when considering the whole Keaton Beach image, which had a mean of 3.38 m 2 leaf m −2 ground.At all sites, the imagery effectively captured density variations across the underwater landscape, which can provide crucial insights into the spatial ecology of seagrass meadows.The ability to differentiate between sparse and dense seagrass can be critical in the man-agement of seagrass resources, as sparse areas are more susceptible to disturbances, and therefore, the retrieval of density from satellites provides a tool to aid in the early detection of changes.
The underestimation of LAI (DOS) was caused by the residual reflectance in band MS3 (green) remaining after calibration/atmospheric correction.The strong correlation between LAI (DOS) and LAI (ELH) provides a mechanism for correcting LAI (DOS ) to values equivalent to those retrieved from an image with accurate atmospheric correction (i.e., ELH).Both the comprehensive regression encompassing all datasets and the individual regressions for each dataset yielded high accuracy in correcting LAI (DOS) .This study therefore indicates that there exists a model to correct the underestimation of LAI (DOS) that is stable across a 2-year gap and 220 km of the Gulf of Mexico coastline, which supports a considerable range of water column optical properties.LAI (DOS) was also underestimated in a study utilizing Landsat 5 images of Saint Joseph Bay [5].The offset in that case was a simple upward correction of LAI (DOS) of 0.98, indicating that while the underestimation of LAI from incomplete atmospheric correction is prevalent, the LAI correction may be sensor-specific.
WorldView-2 was tasked to collect the images used in this study, but an archive of images extending back to 2010 is available over other coastal areas with seagrass.Although it appears that the use of published radiometric calibration values and DOS atmospheric correction are sufficient for the reliable classification of areal distribution, the retrieval of biogeochemically relevant products such as LAI and, therefore, AGC seagras) may be underestimated without in situ measurements to apply accurate calibration/atmospheric correction.Ultimately, the robustness of this correction needs to be further studied if it is to be a pathway for processing the WorldView-2 archive for density and Blue Carbon estimates.Further studies of LAI retrieved from WorldView-2 using available atmospheric correction methods will be needed to determine whether underestimation is consistent across geographical regions and sensors.

Conclusions
WorldView-2 provides sufficient spectral resolution to separate spectral characteristics of seagrass from the surrounding benthic targets in optically complex waters.Accuracy in atmospheric correction was found to be less important than precision in enabling the identification of seagrass.Therefore, images could be classified routinely even if atmospheric contaminants remain in the data.Incomplete atmospheric correction appears to be universal in the blue and green regions of the spectrum and not just in WorldView-2 images; therefore, we expect that without accounting for this omission, the seagrass density will be consistently underestimated from coastal imagery.As coincident in situ R rs measurements are not available with most image collections, the development and testing of an empirical correction to LAI is a critical step for carbon retrieval from WorldView-2 and other sensors.
WorldView-2 and other high-spatial-resolution satellites provide the ability to retrieve the absolute density on a per-pixel basis, as opposed to classifying polygons using relatively broad categories of percent cover.The benefit of retrieving LAI lies in the strong correlation between the density and above-ground seagrass carbon, which in turn can be used to estimate below-ground carbon.The percent seagrass cover is a commonly retrieved attribute that can be estimated via manual delineation of aerial imagery or through sea surface reflectance, but unlike LAI, it is not strongly predictive of biomass [19].For this reason, we suggest that finding an appropriate correction for LAI when residual brightness remains in the spectra after atmospheric correction or finding a way to improve our atmospheric correction is critical in moving our estimates forward for the seagrass contribution to Blue Carbon.

Figure 1 .
Figure 1.Mercator projection showing the three study regions located along the Florida Gulf Coast, USA, and the three WorldView-2 images analyzed in this study.From left to right, Saint Joseph Bay, Saint George Sound, and Keaton Beach.Underlying imagery sources: CONANP, Esri, Garmin, FAO, NOAA, USGS, EPA, Earthstar Geographics.

Figure 1 .
Figure 1.Mercator projection showing the three study regions located along the Florida Gulf Coast, USA, and the three WorldView-2 images analyzed in this study.From left to right, Saint Joseph Bay, Saint George Sound, and Keaton Beach.Underlying imagery sources: CONANP, Esri, Garmin, FAO, NOAA, USGS, EPA, Earthstar Geographics.

Figure 2 .
Figure 2. RGB images of all three study sites, stretched using histogram equalization (γ   0.5) to assist visualization of spectral differences.(A).St.Joseph Bay. (B).St. George Sound.(C).Keaton Beach.Mapped are the locations of stations used for atmospheric correction (white) and those that were fully characterized with respect to water column optics and seagrass abundance (purple) and Regions of Interest (ROIs) used for DCNN classification.

Figure 2 .
Figure 2. RGB images of all three study sites, stretched using histogram equalization (γ = 0.5) to assist visualization of spectral differences.(A).St.Joseph Bay. (B).St. George Sound.(C).Keaton Beach.Mapped are the locations of stations used for atmospheric correction (white) and those that were fully characterized with respect to water column optics and seagrass abundance (purple) and Regions of Interest (ROIs) used for DCNN classification.

Figure 4 .
Figure 4. Mean spectra for (A) absorption (apg), (B) sca ering (bp), and (C) beam a enuation (cpg) coefficients measured using Sea-Bird/WetLabs ac-9 instrument in Saint Joseph Bay, Saint George Sound, and Keaton Beach.Error bars indicate ±1 standard error of the mean for each wavelength bin.

Figure 5 .
Figure 5.Comparison between Rrs retrieval for ELH and DOS atmospheric correction for WV-2 bands MS2 through MS7.Data for each site are indicated by symbol color.Do ed line represents the 1:1 line of perfect agreement.Data consist of 300 randomly chosen points from each site, 100 from seagrass, 100 from submerged sand, and 100 from deep water.

Figure 4 .
Figure 4. Mean spectra for (A) absorption (a pg ), (B) scattering (b p ), and (C) beam attenuation (c pg ) coefficients measured using Sea-Bird/WetLabs ac-9 instrument in Saint Joseph Bay, Saint George Sound, and Keaton Beach.Error bars indicate ±1 standard error of the mean for each wavelength bin.

Figure 4 .
Figure 4. Mean spectra for (A) absorption (apg), (B) sca ering (bp), and (C) beam a enuation (cpg) coefficients measured using Sea-Bird/WetLabs ac-9 instrument in Saint Joseph Bay, Saint George Sound, and Keaton Beach.Error bars indicate ±1 standard error of the mean for each wavelength bin.

Figure 5 .
Figure 5.Comparison between Rrs retrieval for ELH and DOS atmospheric correction for WV-2 bands MS2 through MS7.Data for each site are indicated by symbol color.Do ed line represents the 1:1 line of perfect agreement.Data consist of 300 randomly chosen points from each site, 100 from seagrass, 100 from submerged sand, and 100 from deep water.

Figure 5 .
Figure 5.Comparison between R rs retrieval for ELH and DOS atmospheric correction for WV-2 bands MS2 through MS7.Data for each site are indicated by symbol color.Dotted line represents the 1:1 line of perfect agreement.Data consist of 300 randomly chosen points from each site, 100 from seagrass, 100 from submerged sand, and 100 from deep water.

Figure 6 .
Figure 6.Habitat classification maps for Saint Joseph Bay using (A) ELH atmospheric correction and (B) DOS atmospheric correction.Habitat class is indicated by color identified in the legend in (A).(C) Comparison of pixels identified as seagrass using ELH and DOS atmospheric correction.Pixels mapped in green were identified by both.Pixels mapped in blue and red were identified only by DOS and ELH, respectively.

Figure 7 .
Figure 7. Habitat classification maps for Saint George Sound using (A) ELH atmospheric correction and (B) DOS atmospheric correction.(C) Comparison of pixels identified as seagrass using ELH and DOS atmospheric correction.Pixels mapped in green were identified by both.Pixels mapped in blue and red were identified only by DOS and ELH, respectively.The rectangular and circular shapes highlight areas that are mentioned in the text.

Figure 6 .
Figure 6.Habitat classification maps for Saint Joseph Bay using (A) ELH atmospheric correction and (B) DOS atmospheric correction.Habitat class is indicated by color identified in the legend in (A).(C) Comparison of pixels identified as seagrass using ELH and DOS atmospheric correction.Pixels mapped in green were identified by both.Pixels mapped in blue and red were identified only by DOS and ELH, respectively.

Figure 6 .
Figure 6.Habitat classification maps for Saint Joseph Bay using (A) ELH atmospheric correction and (B) DOS atmospheric correction.Habitat class is indicated by color identified in the legend in (A).(C) Comparison of pixels identified as seagrass using ELH and DOS atmospheric correction.Pixels mapped in green were identified by both.Pixels mapped in blue and red were identified only by DOS and ELH, respectively.

Figure 7 .
Figure 7. Habitat classification maps for Saint George Sound using (A) ELH atmospheric correction and (B) DOS atmospheric correction.(C) Comparison of pixels identified as seagrass using ELH and DOS atmospheric correction.Pixels mapped in green were identified by both.Pixels mapped in blue and red were identified only by DOS and ELH, respectively.The rectangular and circular shapes highlight areas that are mentioned in the text.

Figure 7 .
Figure 7. Habitat classification maps for Saint George Sound using (A) ELH atmospheric correction and (B) DOS atmospheric correction.(C) Comparison of pixels identified as seagrass using ELH and DOS atmospheric correction.Pixels mapped in green were identified by both.Pixels mapped in blue and red were identified only by DOS and ELH, respectively.The rectangular and circular shapes highlight areas that are mentioned in the text.

Figure 8 .
Figure 8. Habitat classification maps for Keaton Beach using (A) ELH atmospheric correction and (B) DOS atmospheric correction.Habitat class is indicated by color identified in the legend in Figure 6A.(C) Comparison of pixels identified as seagrass using ELH and DOS atmospheric correction.Pixels mapped in green were identified by both.Pixels mapped in blue and red were identified only by DOS and ELH, respectively.

Figure 9 .
Figure 9.In situ measurements of Rb(green) and leaf area index from all three study sites.LAI values are averaged from randomized sampling within 20 m of a central location and are shown with ±1 standard error bars.

Figure 8 .
Figure 8. Habitat classification maps for Keaton Beach using (A) ELH atmospheric correction and (B) DOS atmospheric correction.Habitat class is indicated by color identified in the legend in Figure 6A.(C) Comparison of pixels identified as seagrass using ELH and DOS atmospheric correction.Pixels mapped in green were identified by both.Pixels mapped in blue and red were identified only by DOS and ELH, respectively.

Figure 8 .
Figure 8. Habitat classification maps for Keaton Beach using (A) ELH atmospheric correction and (B) DOS atmospheric correction.Habitat class is indicated by color identified in the legend in Figure 6A.(C) Comparison of pixels identified as seagrass using ELH and DOS atmospheric correction.Pixels mapped in green were identified by both.Pixels mapped in blue and red were identified only by DOS and ELH, respectively.

Figure 9 .
Figure 9.In situ measurements of Rb(green) and leaf area index from all three study sites.LAI values are averaged from randomized sampling within 20 m of a central location and are shown with ±1 standard error bars.

Figure 9 .
Figure 9.In situ measurements of R b (green) and leaf area index from all three study sites.LAI values are averaged from randomized sampling within 20 m of a central location and are shown with ±1 standard error bars.

Figure 10 .
Figure 10.Maps of leaf area index (LAI) retrieval for Saint Joseph Bay calculated from (A) Rrs(ELH) and (B) Rrs(DOS).

Figure 10 .Table 8 .
Figure 10.Maps of leaf area index (LAI) retrieval for Saint Joseph Bay calculated from (A) R rs(ELH) and (B) R rs(DOS) .Table 8. Results of pairwise Mann-Whitney tests to compare satellite-based LAI values with in situ observations for each study site.Raster values were collected from a 20 m radius around the in situ station locations (Figure 1).Symbol * indicates a significant difference.Saint Joseph Bay Saint George Sound Keaton Beach LAI (INS ITU) LAI (ELH) LAI (DOS) LAI (INS ITU) LAI (ELH) LAI (DOS) LAI (INS ITU) LAI (ELH) LAI (DOS) LAI (INS ITU) U = 280,464 p = 0.118 U = 423,261 p = <0.00*

Figure 11 .
Figure 11.Box-and-whisker plots for LAI measured in situ and retrieved from LAI(ELH) and LAI(DOS) from pixels collected from a 20 m radius around the in situ station locations (Figure 1).(A).Saint Joseph Bay, (B).Saint George Sound, (C).Keaton Beach.

Figure 11 . 26 Figure 12 .
Figure 11.Box-and-whisker plots for LAI measured in situ and retrieved from LAI (ELH) and LAI (DOS) from pixels collected from a 20 m radius around the in situ station locations (Figure 1).(A).Saint Joseph Bay, (B).Saint George Sound, (C).Keaton Beach.Remote Sens. 2023, 15, x FOR PEER REVIEW 18 of 26

Figure 12 .
Figure 12.Maps of leaf area index (LAI) retrieved for Saint George Sound calculated from (A) R rs(ELH) and (B) R rs(DOS) .

Figure 13 .
Figure 13.Maps of leaf area index (LAI) retrieval for Keaton Beach calculated from (A) R rs(ELH) and (B) R rs(DOS) .

Figure 14 .
Figure 14.LAI(DOS) versus LAI(ELH) from 1000 randomly selected pixels from each site.(A).All sites combined.(B).Saint Joseph Bay. (C).Saint George Sound.(D).Keaton Beach Do ed line represents the 1:1 line, and points above this line are higher in the LAI(ELH) retrieval.

Table 1 .
Summary of symbols, their definitions, and their dimensions.

Table 3 .
WorldView-2 image information for all three sites.Data from NOAA tide stations were used to correct mean lower low water levels to water depth at the time of the image collection.

Table 4 .
Slope (95% confidence intervals) of the regression between R rs(ELH) and R rs(DOS) for each band and site.

Table 5 .
Results of the ANCOVA statistical test to determine differences in slopes between R rs(DOS) and R rs(ELH) between study sites.Only results for comparison of slopes are shown.F ratio is the ANCOVA test statistic.

Table 6 .
Retrieved area of benthic targets for empirical line height (ELH) and dark-object subtraction (DOS) atmospherically corrected images for Saint Joseph Bay, Saint George Sound, and Keaton Beach.

Table 7 .
Retrieved leaf area index (LAI) and derived above-ground carbon from Saint Joseph Bay, Saint George Sound, and Keaton Beach.Results are shown for LAI (ELH) and LAI (DOS).Statistics are calculated from all pixels classified as seagrass with LAI greater than 0 m 2 m −2 in each raster.

Table 9 .
Regression parameters for data presented in Figure13.Values in parentheses are standard errors of the slope and intercept.

Table 10 .
Statistical output comparing accuracy of regression models for predicting LAI from LAI(DOS) when using regression relationships detailed in Table9.

Table 9 .
Regression parameters for data presented in Figure13.Values in parentheses are standard errors of the slope and intercept.

Table 10 .
Statistical output comparing accuracy of regression models for predicting LAI from LAI (DOS) when using regression relationships detailed in Table9.