Structural Characteristics of the Tallest Mangrove Forests of the American Continent: A Comparison of Ground-Based, Drone and Radar Measurements

The Panama Bight eco-region along the Pacific coast of central and South America is considered to have one of the best-preserved mangrove ecosystems in the American continent. The regional climate, with rainfall easily reaching 5–8 m every year and weak wind conditions, contribute to the exceptionally tall mangroves along the southern Colombian and northern Ecuadorian Pacific coasts (Nariño Department and Esmeraldas Province areas). Here we evaluate the use of different methods (ground-based measurements, drone imagery and radar data [Shuttle Radar Topography mission-SRTM and TanDEM-X]) to characterize the structure of the tallest of these forests. In November 2019, three mangrove sites with canopy heights between 50 and 60 m, previously identified with SRTM data, were sampled close to the town of Guapi, Colombia. In addition to in situ field measurements of trees, we conducted airborne drone surveys in order to generate georeferenced orthomosaics and digital surface models (DSMs). We found that the extensive mangrove forests in this area of the Colombian Pacific are almost entirely composed of Rhizophora spp. trees. The tallest mangrove tree measured in the three plots was 57 m. With ca. 900 drone photographs, three orthomosaics (2 cm pixel–1 resolution) and digital surface models (3.5 cm pixel–1) with average area of 4,0 ha were generated. The field-measured canopy heights were used to validate the drone-derived and radar-derived data, confirming these mangrove forests as the tallest in the Americas. The drone-derived orthomosaics showed significant patches of the Golden Leather Fern, Acrostichum aureum, an opportunistic species that can be associated to mangrove degradation, indicating that the mangrove forests investigated here may be threatened from increased selective logging requiring improvements and effective implementation of the current mangrove management plans in Colombia. The techniques used here are highly complementary and may represent the three tiers for carbon reporting, whereby the drone-derived canopy height maps, calibrated with local in situ measurements, provides cheap but reliable Tier 3 estimates of carbon stocks at the project level.

The Panama Bight eco-region along the Pacific coast of central and South America is considered to have one of the best-preserved mangrove ecosystems in the American continent. The regional climate, with rainfall easily reaching 5-8 m every year and weak wind conditions, contribute to the exceptionally tall mangroves along the southern Colombian and northern Ecuadorian Pacific coasts (Nariño Department and Esmeraldas Province areas). Here we evaluate the use of different methods (groundbased measurements, drone imagery and radar data [Shuttle Radar Topography mission-SRTM and TanDEM-X]) to characterize the structure of the tallest of these forests. In November 2019, three mangrove sites with canopy heights between 50 and 60 m, previously identified with SRTM data, were sampled close to the town of Guapi, Colombia. In addition to in situ field measurements of trees, we conducted airborne drone surveys in order to generate georeferenced orthomosaics and digital surface models (DSMs). We found that the extensive mangrove forests in this area of the Colombian Pacific are almost entirely composed of Rhizophora spp. trees. The tallest mangrove tree measured in the three plots was 57 m. With ca. 900 drone photographs, three orthomosaics (2 cm pixel −1 resolution) and digital surface models (3.5 cm pixel −1 ) with average area of 4,0 ha were generated. The field-measured canopy heights were used to validate the drone-derived and radar-derived data, confirming these mangrove forests as the tallest in the Americas. The drone-derived orthomosaics showed significant patches of the Golden Leather Fern, Acrostichum aureum, an opportunistic species that can be associated to mangrove degradation, indicating that the mangrove forests investigated here may be threatened from increased selective logging requiring improvements and effective implementation of the current mangrove INTRODUCTION Spaceborne remote sensing (SRS) is nowadays a key tool to understand Earth phenomena including the monitoring of biodiversity on our planet (Pettorelli et al., 2016). This is the case for the use of SRS in monitoring mangrove ecosystems. Accelerating since 2000 and facilitated by the free availability of satellite products, different global datasets of mangrove cover and loss are now available (Worthington et al., 2020). These baseline products (e.g., Giri et al., 2011;Bunting et al., 2018) have helped to derive global datasets of important properties of mangrove such as above-ground biomass (AGB) and carbon (AGC) after extrapolating or averaging site-specific values (Rovai et al., 2016;Rovai et al., 2018;Hamilton and Friess, 2018). Other SRS products like Lidar and interferometric synthetic aperture radar (InSAR) have been recently used to estimate mangrove heights at national, regional and global scales greatly improving previous AGB estimates (Fatoyinbo and Simard, 2013;Shapiro et al., 2015;Lagomasino et al., 2019;Simard et al., 2019). Nevertheless, there are still limitations faced when using some of these SRS global products with respect to their resolution (usually > 25 m per pixel) and the climatic conditions under which images can be taken (e.g., cloud free conditions).
Ideally, local and national mangrove mapping datasets that use remote sensing products with higher resolutions and that incorporate ground-truth measurements can play an important role in providing training data and validation for many of the on-going global mapping exercises (e.g., Global Mangrove Watch initiative; Worthington et al., 2020). Hence coordination and integration of global, national, and local datasets are needed to overcome many of the limitations of individual datasets.
Consumer-grade drones, belonging to the unoccupied aircraft systems (UAS) sector, appear to be an appropriate tool to overcome some of the limitations of SRS products in mangrove areas that are < 100 ha (Castellanos-Galindo et al., 2019). Examples of mapping mangrove forest with drones at spatial resolutions with centimeter scale are rapidly appearing in the literature (Otero et al., 2018;Li et al., 2019;Navarro et al., 2020). These finer resolution products provide unique opportunities to understand ecological processes and derive essential biodiversity variables that would not be possible to obtain with SRS at the moment. The use of UAS in mangrove mapping is therefore likely to simplify and complement traditional field inventories and additionally increase the accuracy of AGB and AGC estimates at the local level (e.g., Lucas et al., 2020). Collecting drone-derived information can therefore greatly benefit the calibration of the regional and global estimates that have been produced so far using SRS products of coarser resolution.
The Panama Bight eco-region in the west coast of Central and South America is considered to have one of the bestpreserved mangrove ecosystems in the American continent and is thus a widely recognized conservation hotspot (Spalding et al., 2007). Highly developed mangrove forests are found along the southern Colombian and northern Ecuadorian Pacific coasts (Nariño and Esmeraldas areas; Castellanos-Galindo et al., 2015;Hamilton et al., 2017). Already in the mid-twentieth century, the American Geographer Robert West referred to these intertidal forests as the "most luxuriant mangroves of the world" (West, 1956). A more recent study (Simard et al., 2019) has identified the southern Colombian Pacific coast (Cauca and Nariño Departments), together with Gabon and Equatorial Guinea in Africa, as the regions containing the tallest mangroves in the world. A combination of high precipitation (i.e., potentially lower salinity), high temperature and low cyclone landfall frequency plus local geomorphological factors are associated to the observed exceptionally large areas and high canopy heights observed in those regions (Simard et al., 2019).
Mangroves in the southern Colombian Pacific coast occur in areas with geomorphological features like alluvial plains that form dynamic barrier islands surrounded by tidal channels (Martínez et al., 1995). These mangroves may be the wettest in the world with rainfall easily reaching 5 to ca. 8 m every year. The remoteness of this vast area (with no coastal road and only accessible by long boat journeys or small airplanes) contributes to the relatively pristine nature of some of these forests, but also to the difficulty in obtaining scientific ecological information (Castellanos-Galindo et al., 2021). This difficulty has prevented so far corroborating most of the information from modeling studies that have highlighted the unique characteristics of the mangrove areas in this coast (e.g., Hutchison et al., 2014;Rovai et al., 2018).
To validate previous remote sensed measurements (i.e., Simard et al., 2019) and to assess the complementarity of different techniques, we compare here mangrove tree height and AGB measurements from ground-based inventories, drone imagery and radar data (from the Shuttle Radar Topography mission-SRTM and the TerraSAR-X add-on for Digital Elevation Measurement-TanDEM-X) taken from a remote area in the Colombian Pacific coast recognized to contain the tallest mangroves in the American continent. The combination of fieldbased methods (tree DBH and height measurements and drone photogrammetry) with existing SRS products in the most extreme range of mangrove canopy height (∼50 m) allows for the first time: (1) the validation of previous models that highlighted the ecological value of mangroves in this region and (2) the recognition that these areas are in need of urgent protection due on-going localized alterations.

Study Site
The Colombian Pacific coast encompasses ca. 1,500 km in the tropical eastern Pacific biogeographic region that extends from the Gulf of California, Mexico to northern Peru. Almost 2/3 of this coastal region is dominated by mangroves, with ∼80% of the total mangrove area of Colombia (Mejía-Rentería et al., 2018). The whole coast presents annual precipitations of > 2,000 mm, reaching in some areas 8,000 mm yr −1 . This extreme precipitation translates into year-long low water salinities (<30) and permanent brackish areas inside bays, small deltas and tidal channels in what is identified as the Eastern Pacific fresh pool (Alory et al., 2012). The tidal regime is semidiurnal, meso-to macro-tidal with amplitudes during spring tides greater than 4 m (Correa and Morton, 2010). These tidal conditions indicate that mangroves and mudflats are completely exposed at low tides and completely inundated at high tide (during both neap and spring tides).
Nariño, the southernmost coastal Department along the Pacific coast of Colombia, borders the neighboring Esmeraldas province of Ecuador to the South, and the Colombian Cauca Department in Colombia to the North. Nariño alone contains 46% of the total mangrove forests of the country (Mejía-Rentería et al., 2018). This tectonically active coast is characterized by the presence of two high water discharge and dynamic deltas (Mira and Patía; Restrepo and Cantera, 2013) and an almost uninterrupted mangrove belt. The coast is also underdeveloped in terms of infrastructure with the presence of small villages (<5,000 inhabitants) scattered in a mangrovedominated landscape. Nariño contains the largest mangrove marine protected area (MPA) of the whole tropical eastern Pacific (the Sanquianga National Park with ∼80,000 ha) and a more recently created mangrove MPA bordering Ecuador (Distrito Nacional de Manejo Integrado Cabo Manglares Bajo Mira y Frontera). North of Sanquianga MPA lies the mouth of the Iscuandé River, an estuarine area dominated by extensive mangrove forests and small human settlements that constitute the Afro-Colombian community council of Esfuerzo Pescador. The water discharge of the Iscuandé River is 6.71 km 3 yr −1 (Restrepo and Kjerfve, 2004).

Field Sampling
In November 2019 (6th-8th), three different mangrove sites in the northern coast of the Nariño Department of the Colombian Pacific coast were visited (Figure 1). The sites are located near the mouth of the Iscuandé River, an estuarine area completely dominated by mangrove forests, tidal creeks and mud-sand flats that emerge during low tides. Circular plots of 15 m-diameter were established at each site. Ground measurements within these plots included the identification of tree species, diameter at breast height (DBH) and top height of selected trees. While DBH measurements were generally made around the trunk of each tree at 1.3 m above the ground surface, several measurements were made above the high prop roots of Rhizophora spp. (e.g., Figure 2). To measure DBH, a diameter tape was wrapped around the girth of each tree and a diameter measurement read was made to the nearest 0.1 cm. Heights of selected trees within each plot were collected with the aid of a Nikon Forestry Pro II laser hypsometer. There trees were then geolocated using a handheld GPS (Bad Elf GNSS Surveyor).

Allometric Equations
For each of the plots surveyed, we used the DBH values and two allometric equations to determine tree biomass per hectare. First, we used the allometric equation for Rhizophora mangle developed by Fromard et al. (1998) to estimate tree biomass that only takes into account DBH: We also used the pantropical tree allometric model to estimate AGB proposed by Chave et al. (2014) which takes into account tree height, DBH and wood specific gravity: where D is DBH in cm, H is tree height in m and ρ is wood specific gravity in g cm −3 . ρ was computed as the mean value of the wood specific gravity values available for Rhizophora spp. in the Central and South America from the Global wood density database (Chave et al., 2009;Zanne et al., 2009).
To implement the pantropical tree allometric model of Chave et al. (2014), we back-calculated the heights of the trees for which we only had DBH. For that we fitted, a Log-Log model of DBH:H with 26 trees for which we had both DBH and Heights (see Supplementary Figure 1). Total living carbon was then calculated using the relationship total living biomass:living carbon of 1:0.464 (Kauffman et al., 2011(Kauffman et al., , 2020.

Drone Surveys
Drone surveys in the direct vicinity of each plot were carried out at low tide with calm winds. The drone patches surveyed were ∼3 ha in size at each site and always covered the plot areas where ground measurements were taken. Low-altitude (100 m) aerial images were acquired using two UAS, DJI Mavic Pro and Mavic 2 Pro. These drones are equipped with integrated photo cameras. DJI Mavic Pro has a FC220 camera, with a sensor 1/2.3" CMOS (effective pixel:12,35 Megapixel), focal length of 5 mm, a pixel size of 1.7 × 1.7 µm and a resolution of 4,000 × 2,250 pixels. The Mavic 2 Pro has a L1D-20c camera with a sensor 1" CMOS camera (effective pixel:20 Megapixel), focal length of 10 mm, a pixel size of 2.41 × 2.41 µm and a resolution of 5,472 × 3,648 pixels. Photos are geotagged by the DJI integrated GPS/GLONASS location system with a horizontal and vertical precision of up to ± 1.5 and ± 0.5 m, respectively. 1 Flights were performed in automated mode and were programmed using the commercial web app "Drone Deploy." Due to the nature of the environment, characterized by dense vegetation with lack of large gaps exposing the terrain within the forest, it was not possible to easily place and collect ground control points (GCPs). Although few small gaps are present, they are not easily accessible and lack FIGURE 1 | Location of the mangrove sites visited in the southern Colombian Pacific coast (Iscuandé River mouth, Nariño Department; Imagery ©2021 Google, Maxar Technologies, Data SIO, NOAA, U.S. Navy, NGA, GEBCO, Landsat/Copernicus, Imagery ©2021 Terrametrics, Map data ©2021). Three mangrove patches were surveyed in the estuarine complex formed close to the mouth of the Iscuandé River: from left to right, sites are known as: "La Punta del Codo," "La Rotura," and "Madrid." the necessary sky view to operate GNSS and collect GCPs with acceptable accuracy.
The acquired images were analyzed and processed using the software Agisoft Metashape. 2 Metashape is based on Structure from Motion (SfM) (Ullman, 1979) and Multi-View Stereo reconstruction (MVS) methods (Scharstein and Szeliski, 2002;Seitz et al., 2006). For a comprehensive description of the SfM method implemented in Metashape, the reader 2 www.agisoft.com is referred to Westoby et al. (2012). Metashape gives as outputs an orthorectified photomosaic (orthomosaic) and a Digital Surface Model (DSM) from nadir photos collected during the flight. The output models were georeferenced in Metashape to the WGS84 datum using camera positions. The results of the photogrammetric suite were used in the three locations to estimate different environmental variables. Orthomosaics with 2 cm/pixel resolution allowed detecting the heterogeneity of vegetation and the presence of opportunistic species associated to mangrove degradation. DSMs were used to estimate the heights of trees within the in situ sampling areas. Since the ground was not visible within the forest, we used the elevation of mud-sand flats that emerge during low tides as a reference. The observed DSM was vertically shifted to match the elevation of the mud-sand flats. Consequently, the height of the trees was extracted from the vertically shifted DSM.

Satellite Data
To compare our in situ and drone generated data, we used two global RS products generated by Simard et al. (2019). One of maximum mangrove canopy height for the year 2000 with a resolution of ∼30 m that used the Shuttle Radar Topography Mission (SRTM) global digital elevation model (DEM), and the Geoscience Laser Altimeter System (GLAS) global Lidar altimetry products. The second RS product was an aboveground mangrove biomass map that was generated by Simard et al. (2019) linking field-measured biomass-height allometry with SRTM estimates of basal area weighted height (all RS products available at https: //daac.ornl.gov/cgi-bin/dsviewer.pl?ds_id=1665).
We additionally compared our mangrove height fieldgenerated data against the high-resolution data product (12 × 12 m) derived from the TerraSAR-X add-on for Digital Elevation Measurement-TanDEM-X mission (Simard unpublished data). We performed a similar procedure as with the drone DSMs and adjusted the TanDEM-X DEM shifting it vertically to match the elevation of the adjacent water. For details about use of the TanDEM-X in mangrove canopy height estimations see Lee and Fatoyinbo (2015) and Lagomasino et al. (2016).
The tree height comparisons between the three methods used was done using the ground height measurements taken at the three sites surveyed (n = 21 trees). With the georeferenced location of each of these trees, a vector (∼5 m radius circle) was created in QGIS and using the zonal statistics, height information (mean, maximum and minimum) from the drone, the SRTM and the TanDEM-X rasters were extracted. We fitted linear regressions between in situ maximum tree heights and calculated drone, SRTM and TanDEM-X tree heights. In addition, Pearson, determination coefficients and RMSE (root mean squared error) were computed to assess the fit of the regressions.

RESULTS
All surveyed sites were dominated by the red mangrove Rhizophora spp. In one of the sites (Madrid; Figure 3), the golden leather fern Acrostrichum aureum, a mangrove associate, dominated the understory which in many cases reached up to 4 m (see Figures 2, 4). While still difficult to distinguish in the drone orthomosaic, we estimated that at least 5% of the orthomosaic obtained for the Madrid site was completely covereded by A. aurem (i.e., gaps with no Rhizophora spp. canopies). Moreover, our ground surveys at this site showed that the understory in the remaining 95% of the orthomosaic was completely dominated by this species, which was not captured by the drone orthomosaic (Figure 4).
Mean DBH values were very high at the Madrid site in comparison with the other two surveyed sites. The Madrid site was characterized by the largest trees with a mean DBH of 52.4 cm reaching a DBH and height of 96.7 cm and 57 m, respectively. In contrast, the mean DBH at La Rotura and La Punta del Codo sites, were considerably lower (around 20 cm) with maximum tree heights of 50.8 and 53.6 m, respectively. The low mean DBH values indicate the presence of comparatively smaller trees in those two latter sites. Nonetheless, the total AGB was, as expected, dominated by the presence of the largest trees in all sites. Depending on the allometric equation used (Fromard et al., 1998;Chave et al., 2014), the mean AGB across all sites was 862.2 or 626.4 Mg ha −1 , respectively. Similarly, mean above-ground mangrove carbon stocks were estimated in 400.1 or 290.7 Mg C ha −1 ( Table 1). These values, especially in the Madrid site do not account for the contribution of A. aurem to the total biomass and carbon. The mangrove AGB values derived from the ground data were two to three times higher than those estimated with the SRTM data in Simard et al. (2019). AGB values for the examined plots derived from SRTM, and based on continental scale allometry, ranged from 296.2 to 413.3 to Mg ha −1 .

DISCUSSION
We conducted a field campaign specifically to validate and study the regions with the tallest mangrove forests in the Americas (Simard et al., 2019). Field measurements of tree DBH and height were used to validate drone-and SRTM-derived maps of mangrove canopy height, confirming the findings of Simard et al. (2019). The tallest field-measured tree was 57 m, the dronederived canopy height revealed maximum heights of 57.6 m and the SRTM data for these sites showed maximum canopy heights of 54.3 m. While our analysis is performed in forest stands that may not be representative of the entire region, we show that in situ, drone and spaceborne data provide accurate results in the most extreme range of mangrove canopy height. All the three data sources corroborate the existence of exceptionally tall mangrove forests in this region of the Western American tropics. However, each of the methods used to estimate tree height data presented challenges and associated errors ( Figure 5 and Table 2). As revealed in previous modeling studies (e.g., Rovai et al., 2016;Hamilton and Friess, 2018;Rovai et al., 2018), in situ measurements confirmed that these sites are on the upper range of mangrove biomass and carbon per area in the world. Field tree height data was collected with a laser range finder (Nikon Forestry Pro II laser hypsometer). In a recent study, Saliu et al. (2021) indicated that this method presents the lowest amount of error (8%) among a variety of methods used to measure mangrove tree heights in Malaysia. Larjavaara and Muller-Landau (2013) have also identified that both methods (tangent and sine) available when using the laser range finder can incur in high random errors and underestimations in moist tropical forests. A few aspects complicate mangrove tree height calculations with this instrument. For example, in highly dense forests, the identification of the ground and the top of the tree crown can be challenging and thus is not exempt from errors. For a mangrove system in the Colombian Caribbean, Simard et al. (2008) calculated that the random error of tree height measurements in the field can be ∼10%. In the forests examined here, the presence of a dense understory in some sites plus the stilt root nature of Rhizophora trees also complicate the identification of the ground with the laser range finder. While we consider these in situ measurements to be the most reliable estimates of tree height, we acknowledge the potentially significant uncertainty associated with this method. Mangrove Above-Ground Biomass of the Iscuandé Region in the Global Context The mean mangrove AGB values obtained here using the in-situ data are extremely high varying between 862.2 and 626.4 Mg ha −1 (depending on which allometric equation is used). These values are located in the upper range of the recent global synthesis by Rovai et al. (2021) and also correspond to the upper range values recently estimated by Trettin et al. (2021) for a mangrove forest in Gabon where mangrove trees of up to 60 m can be found. Our estimated mean AGB values are generally higher than those obtained using the AGB values calculated by Simard et al. (2019) for the same areas. This result warrants a re-examination of the relationship between SRTM data and in situ inventories considering a larger dataset than the one currently considered in our study. Likewise, it is important to consider that our plots represent very tall mangrove forests within the Iscuandé area. A more representative AGB (and carbon) estimate would need a sampling design that considers the whole range of mangrove type mosaics and patch heterogeneity that are present in this large area (e.g., like the recent work of Trettin et al., 2021 in Gabon). Such a study could yield slightly lower mean AGB values for these mangrove forests but would also contribute to more accurate numbers on carbon stocks of these areas.

Advantages of Drone-Derived Data to Detect Mangrove Degradation
The drone derived data (orthomosaics and DSMs) of these three mangrove sites provided sufficiently detailed information (<4 cm pixel −1 resolution) to detect landscape features that are not captured with the ground surveys nor with the SRTM or the TanDEM-X products. Especially important are (1) the heterogeneity in canopy heights observed and also (2) the detection of canopy gaps that are the result of possible deforestation (i.e., selective logging). On one hand, the drone-derived DSMs (3.5 cm pixel −1 resolution) recorded highly variable tree heights within each of the 3 ha patches. On the FIGURE 6 | Comparison of altimetry products derived from surveys in mangroves of Iscuandé (Colombian Pacific coast). Mangrove patches correspond to "La Punta del Codo," "La Rotura" and "Madrid" from left to right. The upper row shows the vertically shifted DSMs obtained from the drone data. The middle row depicts canopy heights from the TanDEM-X product (Simard unpublished;© DLR 2020) The lower row corresponds to canopy heights from the Simard et al. (2019) dataset.
other hand, the pixels SRTM and TanDEM-X products may include tree crowns of short and tall trees or ground areas, thus these elevation of these mixed pixels represent overall canopy height rather than the height of the tallest tree. This effect may be particularly important in this study (Figure 6), with tree heights measured near forest edges along the shore or around canopy gaps.
On the other hand, the drone orthomosaics were sufficiently detailed to detect the presence of gaps within one of the mangrove patches that were dominated by the golden leather fern A. aureum, an opportunistic species associated to mangrove degradation (Madrid patch in Figure 4). Acrostichum aureum is a highly resistant fern able to colonize the understory of mangrove areas when facilitated either by the presence of natural canopy gaps (see Amir, 2012;Amir and Duke, 2019) or by man-made gaps produced after forest logging. The fern is commonly associated to deforestation and is able to rapidly colonize mangrove areas that have been disturbed, a process that has been defined as cryptic (mangrove) ecological degradation (Dahdouh-Guebas et al., 2005;Blanco et al., 2012). Once established, A. aureum modifies the geomorphology and hydrology, impeding natural regeneration and hindering animal movement (Medina et al., 1990;Biswas et al., 2018). Other aspects of the effects of this fern on ecosystem processes within the mangrove forests are poorly understood. The information provided by the drone products on A. aureum can be valuable as it provides evidence at the scale of 10 s of hectares on loss of mangrove ecosystem health that would otherwise be unnoticed when using coarser SRS products. These highly detailed drone products can also aid in the development of complexity indexes, such as those proposed in Blanco et al. (2001), that could be used to characterize health status and degradation in mangrove forests. Being able to identify the understory vegetation (e.g., A. aureum in this case) of a mangrove area with the drone products implies also that the contribution of this vegetation to the ecosystem biomass and carbon storage, which could become very high, could be accounted for. However, it is unclear how persistent in time A. aureum can become in this region. Surprisingly, allometric equations and information on the carbon content and establishment dynamics of this opportunistic fern are scarce and needed.
The mangrove areas examined here, still contain well developed forests. However, historical mangrove timber exploitation has occurred and continues to occur. Although we do not have data on the magnitude of recent and on-going logging, the local communities identify this activity as an increasingly pressing factor on the surrounding mangrove ecosystems. Moreover, unaccounted as a mangrove degradation source remains the effects that the anti-narcotics national program for eradication of coca crops that the Colombian government undertook in the 2000s with aerial aspersion of glyphosate (and currently seeks to continue). One of the few accounts on the use of herbicides in mangrove forests is that of the severe effects that its use caused on Vietnam forests during the second Indochina war (Westing, 1983). Anecdotical information from locals in Iscuandé reported on the effects that glyphosate had on native flora and fauna around mangroves in the past. The localized effects that both timber logging and glyphosate could have on these exceptionally well-developed mangroves could be studied with the finer resolution products from drones.

Implications for Mangrove Carbon Monitoring in Remote Areas
Our study reveals significant challenges in collecting field and drone data in these mangrove areas. These can be categorized as logistical and geometrical. The access to these sites requires careful travel planning that includes site accessibility by plane and boat in a highly dynamic setting driven by the large tidal range, and close collaboration with local authorities and organizations. The latter facilitates, not only safe passage to sites, but gathering of a wealth of knowledge about local landscapes, cultural practices and customs related to mangroves and water management (see Oslender, 2016). On the other hand, the geometrical challenges are related to the complex terrain with mangroves roots, litter, and mud that constrains the possibility to collect tree heights in a systematic and time-efficient way. To accurately process the drone data, collecting ground control points-GCPs (ideally with a Real-time kinematic-RTK GPS) is desirable. However, identifying sufficiently clear gaps to place those GCPs inside mangrove forests can prove an unachievable task. Ways to overcome this limitation could include using drones that have incorporated an RTK GPS (although this technique does not completely exclude the need of GCPs) and/or using adjacent ecosystems within the mangrove forest (e.g., mudflats) where placing GCPs is feasible. Zeng et al. (2021) suggests that mangrove blue carbon projects in Colombia may only be profitable in the southern Colombian Pacific coast, the area where this study was conducted. Our results, therefore provide valuable assessment methodologies and an initial in situ assessment of carbon stocks for potential blue carbon projects in an area that has been highlighted by modeling studies as a blue carbon hotspot in the world with much need of effective conservation actions. The techniques are complementary and may represent the three tiers for carbon reporting, whereby the drone-derived canopy height maps, calibrated with local in situ measurements, provides affordable and reliable Tier 3 estimates of carbon stocks at the project level. These detailed estimates can in turn help calibrating the coarser Tier 1 and Tier 2 assessments. Additionally, the combination of drone and ground survey data can provide a way to detect and monitor small-scale ecosystem degradation symptoms that so far are not possible to detect with the SRS products currently used for mangrove mapping. This can bring a more detailed understanding of landscape processes occurring within these vulnerable ecosystems that so far have been excessively focused on coarse area measurements.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS
All authors contributed to the article and approved the submitted version.

FUNDING
Funding for fieldwork was provided by WWF Colombia and WWF US (GLOBAL MANGROVE ALLIANCE 510-51000000-4350). Bezos Earth Fund (agreement CN11408) supported article publication costs. GC-G was funded by the Smithsonian Tropical Research Institution postdoctoral fellowship. MS was funded by NASA's LCLUC program.

ACKNOWLEDGMENTS
We thank the members of the community council Esfuerzo Pescador and its Legal representative (Kennedy Caicedo) in Iscuandé for their hospitality and facilitating logistical resources during our field campaign. We also thank the local experts (Carmelo Castillo and others) that guided us through the marvelous mangrove forests of Iscuandé. Part of the research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration (LCLUC program). We thank DLR for providing TanDEM-X DEM data through their science program. This project (Manglares y productividad biológica en la región de Iscuandé-Tapaje, Pacífico colombiano) was part of those approved within the frame of the II Scientific Expedition to the Colombian Pacific -Bocas de Sanquianga (Plan Nacional de Expediciones Científicas Pacífico) organized by the Colombian Commission for the Oceans (Comisión Colombiana del Océano-CCO). We thank the comments made by the two reviewers and the guest associate editor that helped to improve the clarity of this manuscript.