Mangroves on the Edge: Anthrome-dependent Fragmentation Influences Ecological Condition (turbo, Colombia, Southern Caribbean)

Marine protected areas are commonly seen as the most effective strategy for protecting mangroves from external human pressures but little is known about the role of public land-tenure contexts (dense settlements, agricultural or range lands and wild anthromes) on clearing rates, patch properties, and ecological condition. We addressed the following questions using a peri-urban to wild gradient along the anthropogenic coastal-scape in Turbo Municipality (Colombia, Southern Caribbean): Do the different deforestation rates observed under peri-urban, rural, military-protected and wild land-use-and-tenure contexts, promote distinctive fragmentation patterns? Do these patterns influence loggers' access and ultimately ecosystem ecological condition? Loss rate (1938–2009) was the greatest in peri-urban mangroves and positively correlated with urban edge and patch density. Pasture edge was highest in rural mangroves while mean patch area was higher in protected and wild mangroves. An Anthropogenic Disturbance Index (ADI) was strongly correlated with reduced mean patch area and increased patch density, due to increased trampling and logging, that ultimately promoted high densities of thin (diameter: <5 cm) Laguncularia racemosa trees but had no significant effect on the presence of a dominant benthic gastropod. In conclusion, both protection and remoteness were effective in reducing anthropogenic edges and fragmentation, and thus contributed to a high ecological condition in mangroves at a major deforestation hotspot.


Introduction
The extensive human transformation of the biosphere imposes new challenges for the conservation of tropical forests [1] including mangroves [2].Firstly, the plethora of human-transformed biomes or anthromes (sensu [3,4]) suggests that conservation of the natural capital may take place under different spatial contexts observed across the continuum from peri-urban and rural landscapes to remote wilderness areas [1,5].Along this gradient, different land-tenure contexts are observed in public and private lands, thus imposing different threats to conservation targets.Human threats such as deforestation and hunting are increasingly observed around protected areas immersed in remote, rural, peri-urban and urban matrices.Therefore, new conservation targets such as "hybrid" and "novel" ecosystems are being selected, as a result of the discussions about long-term sustainability, due to the valuable services that they still provide across the anthropogenic biosphere [6][7][8][9].As a consequence, tropical conservation biologists must learn from the new opportunities offered by the anthropogenic landscapes, moving away from the paradigm of study sites and conservation reserves in isolation from human influences [9,10].
Preserving mangroves through protected areas, alternative types of management, and restoration have failed in many locations worldwide due to the incapability of such strategies for impeding access of illegal loggers from nearby populated areas (e.g.[11,12]).For instance, this situation has prevailed in some marine protected areas (MPAs) with limited personnel for law enforcement (i.e., paper MPAs) in the Caribbean and Central and South America [2,13].Illegal logging within MPAs and other land-based conservation initiatives are seemingly facilitated by deforestation in adjacent areas, as extensively documented for protected terrestrial forests [14].Therefore, mangrove conservation initiatives inside and outside of MPAs need to face a two-fold challenge: (1) preventing illegal extraction of woody products, and (2) reducing clearance in adjacent forests and buffer zones that eventually facilitates access to previously isolated and protected areas.Since land-grant-based mangrove conservation initiatives in public lands (other than MPAs) may take place in a range of landscape types (from urban to wild), law enforcement depends on the kind of tenure, and ultimately on the anthrome context.For example, local police, environmental authorities, community-based organizations and non-government organizations oversee natural resource protection across the urban-rural gradient.In addition to the above landscape contexts, conservation also takes place within restricted-access lands, such as military facilities [15].In military lands such as coastguard posts and military bases, unauthorized-personnel access is restricted or impeded by surveillance and fencing [15][16][17].However, enforcement of laws protecting mangrove ecosystems and species is absent in many urban to rural landscapes in tropical developing countries [2,8].Accordingly, mangrove conservation inside and outside MPAs largely relies upon the physical properties of the landscape context, particularly on isolation by distance from populated areas [18,19], similarly to those reported for terrestrial forests [20,21].Moreover, since mangroves are often found in highly developed coastal landscapes, they frequently experience strong direct and indirect influences from the neighboring human populations (e.g.[22,23]).
Remote and wild forests, as well as restricted-access public lands (e.g. military facilities), are potentially important scenarios for mangrove conservation but no study has tested this hypothesis.In terrestrial forests, isolated areas and restricted-access public lands exhibit lower clearing rates, and thus less anthropogenic fragmentation occurs [15,24].In isolation, forest fragmentation is lower and edge effects become less important drivers of negative ecological effects upon plant and animal biodiversity [20,21].Consequently, since edge-formation is expected to similarly affect mangrove vegetation and benthic fauna, remote or wild areas and restricted-access lands would serve as sources of propagules and larvae supporting mangrove biodiversity in degraded non-protected public lands (i.e., urban, peri-urban and rural) and become key areas for conservation in a regional or meta-community context.
Unfortunately, the landscape ecology approach has been scarcely applied in mangrove conservation science despite cartographic efforts and accounts of deforestation rates populating the literature [23,25,26].Little is known about how external drivers of mangrove degradation quantitatively interact with coastal-scape mosaic and mangrove-patch properties.Various reports have shown that peri-urban mangroves are structurally stressed due to the proximity to neighboring populated centers.For instance, studies in Kenya [18] have recently reported that peri-urban mangroves are characterized by canopy gaps, tree-diameter distributions biased towards small sizes, short stature, and highly clustered trees.In addition, the richness and abundance of mollusks on urban mangroves in Sydney (Australia) were negatively correlated to mangrove condition and particularly to proximity to residential areas, as opposed to the proximity to a national park [27].Nonetheless, many studies lacked comparisons with reference areas, either isolated or protected, due to reduced sampling effort at a landscape level, and therefore observations were not distributed across the urban-rural or urban-wild gradients as recommended in landscape and regional ecology [28].A few studies (other than global and national cartographic initiatives) have covered spatial extents large enough (>10 1 -10 2 km) to compare mangroves under contrasting landscape contexts (but see Kenya, Bosire et al. [18] and references therein).
According to the evidence on terrestrial forests, deforestation translates to ecological effects by increasing the number of fragments and the length of edges [20,21].For instance, the extensive research conducted in Amazonian forest fragments over the past three decades has documented edge effects on climatic and hydrologic disturbance regimes (particularly wind, fires and droughts), tree mortality, forest fauna, and ecosystem processes [29,30].Conversely, the existing literature on the effects of mangrove deforestation on vegetation and benthic fauna has only shown the localized effects of small-scale deforestations (gaps, walkways, and roads) (see [31,32] and references therein) and has scantly described land-cover/land-use transitions and the consequences on mangrove floristics in single locations (e.g., Sri Lanka, [33]).
Our recent explorations of the southern-most Caribbean mangroves located in the Turbo Municipality (Antioquia, Colombia) have demonstrated that they exist along a continuum from wild (or remote) to urban (encroached) landscape contexts along a >100 km coastline [34,35].A mangrove inventory found that tree diameter was thinner in the vicinity of Turbo City than in isolated, wild areas, due to the illegal selective logging of red mangroves (Rhizophora mangle).Further studies have shown that extensive mangrove areas persist under rural and peri-urban contexts surrounded by crops and pastures [34,36].As a consequence, the peri-urban mangroves of Turbo City are considered a deforestation hotspot, although official data are lacking [36].In these rural and peri-urban contexts, pasture expansion for cattle ranching formed sharp mangrove borders that have been responsible for small-scale ecological edge effects upon macro-benthic species (e.g., gastropods and crabs: [31,32,37]).These direct influences on mangrove ecosystems seem to be more important than indirect influences such as sea level rise and upland deforestation produced by human activities taking place at far distances, because they occur at a faster rate thus causing major negative ecological effects [36].Land-use transitions have been noticed, but not quantified, in other river deltas located in rural contexts to the South of Turbo City, and hence ecological edge effects seem to be occurring at a regional scale.The variety of land-tenure contexts (urban, peri-urban, rural, wild and military-protected) where mangroves persist in the Turbo Municipality provided a unique opportunity to address the following three questions: (1) How different are deforestation rates under different land-tenure contexts?(2) What are the landscape-scale consequences in terms of mangrove fragmentation (patch and edge formation)? and (3) What are the consequences on human access and ecological condition (mangrove structure, species composition, and presence of a dominant gastropods)?

Study Area
The study was undertaken in three mangrove neighborhoods in Turbo Municipality (Punta Yarumal, Punta Las Vacas and Punta Coquito), located along the Southeast coast of the Urabá Gulf (Antioquia State), Southern Caribbean of Colombia (Figure 1).These mangrove neighborhoods are smaller in area in comparison to those located on the Atrato River Delta (West coast), which comprise 78% of total mangrove extent within the Urabá Gulf.Mangroves along the East coast contain ca.11% of gulf's total, are highly patchy or fragmented, and are surrounded by anthropogenic land cover such as extensive pastures and banana crops.The Turbo Municipality exhibits the largest population density, mostly concentrated in the urban area (ca.47,000 inhabitants in 2005; [38]), with the largest along the Urabá Gulf coast and the Southern Caribbean of Colombia [39,40].
The species composition of mangroves in Turbo Municipality is similar to the Greater Caribbean, and particularly the Colombian Caribbean, where Rhizophora mangle, Avicennia germinans, Laguncularia racemosa, Pelliciera rhizophorae and Conocarpus erectus occur, the first two being the most important, forming nearly mono-specific stands in some locations [34,35].Physiographically, three kinds of mangroves are found: fringing, riverine and basin [35].Both fringing and basin mangroves tend to be mono-specific, dominated by R. mangle and A. germinans, respectively.In the study area, basin mangroves are more in contact with anthropogenic land covers such as pastures and crops that have reclaimed the former freshwater wetlands.Indeed, forest structure and area show clear signs of human influence in the vicinity of Turbo City [31,32,36,37].

Spatial Data, GIS Methods and Mangrove Change Calculations
High-resolution (0.3 m pixel) and color aerial photographs were taken for the study area in 2009, comprising a coastal fringe nearly 2 km of width containing mangrove neighborhoods.These aerial photos were obtained during the State-government-sponsored research expedition "Antioquia Expedition 2013" [41].We decided to use the neighborhood concept, defined as a spatial unit grouping mangrove patches found within a well-defined geo-form such as a river delta or a depositional intertidal area [42,43].This concept is synonymous to Eco-series, defined as a spatial level within the coastal landscape hierarchy comprising land cover and land use areas between 1.5 and 25 ha (sensu [42,43]).
These high-resolution photographs were ortho-rectified and assembled into a single block to depict the current mangrove extent and surrounding land cover.We used the ArcGIS 10.2 software to process the spatial information from this ortho-photo and we mapped 10 different types of land cover: natural types such as mangroves, alluvial forests and shrubs; anthropogenic types such as pasture, crops and built area; and geo-forms such as rivers, waterways, beaches and sea.Therefore, we described the full landscape complexity and either the natural or the anthropogenic domination around the mangroves.We also employed the categories recommended by the national geographic authority (IGAC: Spanish acronym for Agustín Codazzi Geographic Institute) in the spatial object catalog and the CORINE land-cover methodology adapted for Colombia [44,45].Due to the sub-metric resolution of the 2009-ortho-photo, the minimum mapped area was 0.01 ha and was allowed to discriminate between natural and anthropogenic coverage within small areas.Finally, we validated the photo-interpretation by ground-truthing over seven transects (3.5 km in total) that crossed all types of land cover.The field data was considered the observed data, while the cartographic data were considered the expected data.Observed and expected data for each polygon were contrasted using a confusion matrix to find overall classification accuracy and Kappa coefficients [46,47].The entire coastal landscape comprised 46.3 km 2 (a 79 km coastline).
Further analyzes were conducted only in four contrasting land-use and land-tenure contexts (Figure 2), in agreement with the anthrome classification [3]:  Peri-urban: Equivalent to "Dense settlements", where significant urban area is found.The landscape is a mixture of private and public areas, but "green areas", "lowlands" and "inter-tidal lands" are public domain. Rural: Equivalent to "Croplands or Rangelands", where dense agricultural or pastoral land cover are dominant, respectively.The "green elements" of the landscape are a mixture of private and public lands. Military-protected: Equivalent to "Semi-natural", where forests with minor human presence occur.It is mostly dominated by native tree species, and land-tenure is public. Wild: Equivalent to "Wildlands", where no land with human populations, agriculture or pastures occur.The land-tenure is public, mostly dominated by coastal freshwater wetlands and mangroves.These four land use and tenure contexts also differed in terms of the proximity to Turbo City as the largest populated coastal area, and being the wild context, the most remote, and the most isolated.Mangrove extent was computed for 1938 and 2009, under the four contexts, within observation windows of 1 km 2 for Rural and Peri-urban, and 2 km 2 for Military-protected and Wild contexts.In addition to the 2009 ortho-photo, color or black and white aerial photographs were obtained for each land use context in 1938, 1961, 1975 and 2004.Each photo was geo-referenced relative to the 2009-ortho-photo to compute either gain or loss of mangrove area change according to Puyrabaud's equation [48].

Mangrove Neighborhood Fragmentation
To describe the influence of land use context upon mangrove neighborhoods, we measured the length of the mangrove edge in contact with both natural and anthropogenic land cover surrounding each patch [49].We computed anthropogenic edge (or edge density) as the sum of all percentages of anthropogenic perimeter for all patches.In addition, the degree of fragmentation of each neighborhood was computed as patch density (the ratio between the number of patches and the sum of areas for all patches) [49].

Mangrove Ecological Condition: Vegetation, Anthropogenic Disturbance Index (ADI) and Dominant Snail Presence
The Point-Centered Quarter Method (PCQM, [50,51]) was employed to describe each patch and neighborhood in terms of species occurrence, diameter at breast height (DBH), density and basal area in 2011.Relative abundance, dominance and frequency were computed to obtain the Importance Value Index (IVI) for each species and transect per land use context.The sampling effort comprised a total of 12 patches, 239 points and 31 transects (oriented perpendicular to the tidal flooding gradient) (Figure 2).The number of patches, points and transects sampled per neighborhood was proportional to the total mangrove area, as follows: 126 points within 14 transects in four rural patches, 17 points within four transects in three peri-urban patches, 70 points within 10 transects in four military-protected patches, and 29 points within three transects in one wild patch.
Associated with the PCQM, we valued the anthropogenic disturbance within each sampling point and patch scoring the degree of influence according to the following categorical variables:  Trampling (T): Presence of human and livestock footprints. Logging (L): Evidence of selective logging on mangrove species, such as stumps and downed wood and logs. Wastes (W): Presence of solid wastes on the forest floor and entangled on the roots. Structures (S): Evidence of human modifications of mangrove hydrology and topography (diggings, infillings, cement canals and pipelines), and other structures related to pastoral activities (fences and troughs).
We combined the scores for each categorical variable in an Anthropogenic Disturbance Index (ADI).Each variable was scored from 0 to 3, where 0 means disturbance absence, 1 means little evidence, 2 means evident and 3 means very evident disturbance.Values for each variable were averaged for each land use context and summarized into the ADI as follows: The ADI ranged from 0 to 12.As a result, transects with values close to 1 exhibited the least anthropogenic disturbance, while those close to 12 exhibited the greatest.The stumps were identified when possible and the DBH measured, to provide further support to the logging component of ADI.
The dominant mangrove snail, Neritina virginea (Prosobranchia: Neritidae), a key component of faunal ecological condition, was selected as the benthic biological indicator of anthropogenic disturbance because it responds negatively on abundance and mean shell size to edge effects [31,32,52].At each sampling point, we recorded the presence of individuals of this species within a radius of 15 m (Figure 2).With this method, we favored the occurrence of the species over the abundance, and it is therefore conservative and less biased by observer's skill.Percent presence was computed for each transect.

Statistical Methods
The deforestation (loss) rates, total and anthropogenic edge length (including total, pasture and urban), total number of patches, patch density, total area and mean patch area, were altogether considered as estimates of mangrove fragmentation for each context.ADI and mean stump DBH for R. mangle were considered as estimates of human access to each transect within each context.Tree density, mean DBH, and IVI per species, and percent presence of gastropods were considered as indicators of ecological condition of each transect within each context.
A Non-Metric Multidimensional Scaling (NMDS) [53] was employed to extract the most important variables correlated with the anthrome classification.Because no replication was possible at the context level, only a single value was available for the 10 patch-based variables included in Tables 1 and 2 (e.g., deforestation rate and total mangrove area).Mean values of transects per context were employed for ecological condition variables measured using the PCQ method.A total of 22 variables were included in the ordination matrix.No significant difference in the ordination pattern was observed when values for each transect per context were included in the ordination, but the stress was improved.The NMDS was run using the vegan package in R [54].Data were square-root transformed and a Wisconsin double standardization was employed.After extracting the most important variables defining the anthrome context, variables were correlated using non-parametric (Spearman Rank) or parametric (Pearson) correlations depending on the adjustment to normality, homocedasticity, and sample size of the variables [53].ADI and ecological condition variables were considered as responses of mangrove fragmentation descriptors.All regression analyses were run in Minitab® v.16.

Deforestation Rates
The overall classification accuracy was 90.1% and the Kappa coefficient was 0.86 thus indicating that the cartographic products for the entire study landscape analyzed in 2009-ortho-photo were highly reliable.However, mangrove cover was fully distinctive from other covers and hence classification error was zero.Mangrove change rates depended on the context during the period 1938-2009 (Figure 3, Table 1).Loss rate was highest in the Peri-urban context (1.2%•year −1 ) and lowest in Military and Wild contexts (<1%•year −1 ).Mangrove net gain was only observed in the Rural context.In this context, despite of the net gain, mangroves have been historically converted to pastures as the deltaic fan has expanded due to increased sedimentation.Mangroves colonized the waterfront over the past six decades due to the relocation of the Turbo River mouth to the North of the city [36].In the Peri-urban context, mangroves were slowly converted to pastures during the first four decades, but mangroves and pastures have been rapidly converted to marginal human settlements since the late seventies (thus mangroves cleared at ca. 2%•year −1 , data not shown).In the Military context, coastal geomorphology was nearly steady, but mangroves were replaced by pastures during the first four decades with the establishment of the facility, and the trend has continued.In the Wild context, mangroves have been lost mostly due to coastline retreat, although interior mangroves were lost due to the expansion of banana croplands.The observed deforestation rates in the Peri-urban context exceeded the global average, while the observed rates in Military and Wild contexts were lower [23,55].Cumulative deforestation rate in Peri-urban mangroves in Turbo City was similar to other Peri-urban sites reported in the literature, considered as hotspots (e.g., Kenya: [18,56,57]).In addition, this figure exceeded that for many of the 10 countries with the largest and best-monitored mangrove areas [23].Globally, Singapore, both a country and a city, has experienced the largest deforestation rates due to the rapid urbanization process [22,23].In rural areas, particularly in Southeast Asia and Ecuador, mangroves have been rapidly converted to shrimp aquaculture ponds [58,59].Conversion to croplands other than rice, coconut and oil palm in Southeast Asia [26,58], and pastures have been less recorded in the literature but seem to be more prevalent in Africa and the Caribbean [55].Therefore, here we reported a novel transition or regime shift previously unnoticed.Mangrove loss associated with shoreline retreat has been recorded in many locations worldwide [60,61].

Mangrove Neighborhood Fragmentation
Anthropogenic edge proportion was significantly larger in the Peri-urban context than in other contexts (Table 2).It was nearly twice that in Rural and Military-protected contexts and >10-fold than that in the Wild context, and it was independent of mean patch area.Patch density was greatest in the Peri-urban context, 10-50 times greater than in other contexts.In the Peri-urban context, the anthropogenic edge was promoted by the construction of wooden houses and other infrastructure, while in the Rural context (to a lesser extent in the others), it was promoted by the establishment of pastures.As a consequence, urban edge length and patch density were the distinctive characteristics of the mangroves in the Peri-urban context, while pasture edge length distinguished those in the Rural context in the NMDS space (Figure 4).Military and Wild contexts exhibited opposite conditions for those characteristics.High total area and mean patch area characterized non Peri-Urban contexts.Annual loss rates of mangroves, however, were not associated with the variables describing neighborhood fragmentation across the study contexts.The length of anthropogenic edge was proportional to the size of the anthropogenic matrix or element surrounding the terrestrial vicinity of the mangrove neighborhoods.This pattern has been broadly documented for terrestrial forests [20,21].The type of the anthropogenic matrix may ease human access to mangroves, as observed in terrestrial forests for the case of hunters and loggers [29,30].In the Peri-urban context, human access is readily explained by the presence of inhabitants close to mangrove areas, while in the Rural context, pastures have favored human access.Further ecological implications of context will be explained in the following sections.A high edge effect is expected, understood as the ratio between the perimeter to core area, therefore easing human access to the mangrove interior.In terms of the classical definition of ecological edge effect, human access and other influences derived from either the presence of pasture or urban lands venture longer distances under unprotected contexts, or marginally under Military or Wild contexts.1 and 2) and ecological condition (Tables 3 and 4, in addition to transect-based ADI data pooled in Figure 5 but shown in Figures 6 and 7).
High fragmentation and small patch area were observed in Peri-urban and Rural contexts.Therefore, military protection and remoteness seem to reduce mangrove fragmentation, as observed in the literature for terrestrial forests [20,21].Fragmentation reinforces the ecological edge effects, because the core (center) of small patches can be easily accessed by people and a greater number of small patches could be impacted.Therefore, most of the extent of small patches is subjected to edge effects or cannot be considered as an ecological or functional interior.

ADI Relative to Context
The ADI scores were greater under the Peri-urban context, while they were lower in both Rural and Military-protected contexts, and the lowest in the Wild context (Figure 5).In particular, trampling and logging evidences were greater in Peri-urban and Rural contexts than in the Military and more significantly than in the Wild.This evidence of human access and tree logging were negatively associated with mangrove stump diameter of R. mangle and L. racemosa, the most logged species region-wide (Table 3).Such association was not observed in A. germinans.Human structures such as fences, drainage canals and small cement structures were more frequent in the Rural context.Solid waste did not show a pattern related with land context but with proximity to Turbo City.Although the ADI seemed to play a secondary role explaining the ordination of the study contexts when the fragmentation descriptors were included (Figure 4), it was significantly correlated with mean patch area (−) and mean patch density (+) (Figure 6).The ADI was greater at Peri-urban and Rural contexts, and therefore seemed to be influenced by anthropogenic edge length, in turn influenced by patch area.It confirms the hypothesis that, as observed in terrestrial forests, surrounding anthropogenic matrices or elements facilitate people access (mostly through land than water) to mangrove areas particularly for logging.Field observations suggest that loggers enter the mangroves from the land to extract firewood, while they access by boat to extract poles, R. mangle and L. racemosa being the target species.

Mangrove Vegetation Relative to ADI
The ADI was a significant predictor of mangrove vegetation structure across the study area as a result of anthrome context (Table 4, Figure 6).Transect-scale density and mean DBH of L. racemosa (white mangrove) showed strong correlations, respectively positive and negative, with ADI.No significant trend was observed for the IVI of both L. racemosa and R. mangle (red mangrove).The NMDS also showed that natural features such as high IVI of P. rhizophorae (piñuelo mangrove) and A. germinans (black mangrove) were distinctive of mangroves in Wild and Rural contexts, respectively (Figure 4).Finally, Peri-Urban mangroves were featured by high densities of L. racemosa, while mangroves across the Rural, Military and Wild gradient exhibited a progressive increase in mean DBH for all species (Figure 4).Anthropogenic disturbance, as a result of anthrome-dependent fragmentation, was a strong driver of mangrove structure in Turbo Municipality.L. racemosa may be employed as the indicator species of logging influence and overharvesting upon R. mangle, which tended to exhibit a reduced IVI in Peri-urban and Rural contexts.Poor ecological condition in mangrove vegetation has been observed in Peri-urban locations worldwide, being the finest examples recorded in Kenya [18].In a previous paper [36], we reported the region-wide negative correlation between the IVI of L. racemosa and R. mangle, and hypothesized that such a species transition was the result of overlogging on the latter.Here, we demonstrate that the degree of selective logging across an urban-to-wild gradient is the mechanism explaining such transition.
The poor ecological conditions observed in mangroves patches surrounded by Peri-urban and Rural contexts could be detrimental not only at a local scale but also for mangroves located near the military-protected area.For instance, Kairo et al. [62] found that mangrove volume and extent adjacent to Kaiunga MPA in Kenya were important for providing regeneration opportunities inside the reserve and sustainable exploitation off the boundaries.Therefore, mangroves within the military-protected area near Turbo City might become progressively ecologically isolated if overlogging and fragmentation continue in unprotected mangroves.
The presence of P. rhizophorae, a vulnerable species according to the International Union for the Conservation of Nature [31,63], in both Military and Wild contexts highlight the importance of law enforcement and geographical isolation, respectively, for preserving imperiled species.This finding is consistent with the high densities of Endangered Species Act-listed taxa and imperiled species in military lands in Hawaiian terrestrial forests [15].
Finally, our results provide the first indicators of mangrove vegetation features under the anthrome classification [3,4].Mangroves in Peri-urban and Rural contexts in our study area may be classified as used ecosystems dominated by L. racemosa, an otherwise secondary species elsewhere in the Urabá Gulf [34,35].Unfortunately, mangroves in the Wild context showed some evidence of anthropogenic influence, particularly at the landscape level (e.g., area loss and proximity to banana crops within <5 km inland and likely access by boat from Turbo City); therefore, it best fits the definition of a remote biome rather than a wild one.

Snail Presence Relative to ADI
The snail presence showed no correlation with the ADI, because it was highly variable across transects, even within the same context (Figure 7).Moreover, high percentages of presence were equally recorded in all contexts, despite the lowest values being recorded outside the Wild context.This result was inconsistent with our previous findings in the rural mangroves to the North of Turbo City where snail frequency linearly declined from the mangrove interior to the adjacent pastures, and it was lower in heavily logged mangrove interiors [34].Therefore, an edge effect upon snail frequency was inconsistent with the land use and tenure categories analyzed in the present study.A closer analysis suggested that snail frequency was greater in the Peri-urban mangroves because deforestation did not contribute to mangrove desiccation and, on the contrary, the houses are built on pilings to avoid flooding, and hence the intertidal habitat seemed little affected.Most of the human occupation on Peri-urban areas took place in fringing mangroves.Therefore, the present study indicates that edge effects are seemingly restricted to the interior-pasture ecotone rather than to the mangrove fringe.Therefore, an analysis of within-patch patterns should be the objective of future research.In addition, a greater effect on logging has been reported on the tree-climbing pulmonate gastropod Melampus coffeus in black mangrove inner zones adjacent to pastures in the vicinity of Turbo [31,32].Regardless of our results, we recommend testing for edge effects on ground-dwelling and tree-climbing snail populations and the ecosystem processes that they control in Caribbean and African mangroves, where gastropod species such as Melampus coffeus and Terebralia palustris, respectively, contribute large biomasses to the benthic consumer-level [64,65].In addition, we recommend testing for edge effects on ground-dwelling crabs, as we recently reported an edge effect on the population of a vulnerable mangrove blue crab (Gecarcinidae: Cardisoma guanhumi) in the same rural mangroves in the vicinity of Turbo City [37].Bioturbation by gecarcinid and grapsid crabs is an important ecosystem process well documented in many mangroves world-wide and it can be impacted by both natural and anthropogenic disturbances [66].

Conclusions and Perspectives
This study concluded that, although different mangrove change rates  were observed in different land-tenure contexts, they were inconsistent with the present-day patterns of mangrove neighborhood fragmentation.However, Peri-urban mangroves were distinctively fragmented as demonstrated by long edges and high density of patches, Rural mangroves were characterized by long edges surrounded by pastures, and both protected and wild mangroves were dominated by large patches.The main effect of fragmentation in Peri-urban and Rural mangroves was the high density of an otherwise secondary species, the white mangrove L. racemosa, promoted by logging upon red mangrove trees.No effect was found on a dominant benthic gastropod.These findings provide lessons for mangrove management as follows.Firstly, surveillance was important for preserving mangrove landscape integrity and, consequently, forest physiognomy in public lands.This is justified because human access to mangroves for logging was eased by the presence of pastures and built areas in the vicinity, and it was the main driver of a poor ecological condition (thin diameter trees) in unprotected land-tenure contexts in Turbo Municipality.Therefore, the main lesson from our study is the requirement of law enforcement (even through co-management) to prevent conversion of mangroves and surrounding freshwater wetlands to pastures and other land covers outside the limits of protected areas and conservation targets.This action would reduce access for illegal logging and extraction of non-woody products, particularly in the proximity to major human settlements.Secondly, this study demonstrated that isolation from major urban centers, by itself, provides a foundation for mangrove conservation.Mangroves in remote and wild locations should be a priority for legal conservation actions such as MPA declaration, particularly where vulnerable fauna and flora are found, such as the piñuelo mangrove (P.rhizophorae) in our case.Thirdly, and equally important, this study demonstrates that mangroves in rural (unprotected) settings may still contain valuable features that are worthy of conservation efforts.These features are the following: (a) the presence of large trees of the black mangrove (A.germinans), a species of localized distribution, indicator of the seasonally dry conditions in the Eastern coast of the Urabá Gulf, and scant in the wetter West Coast, (b) the dominance of the white mangrove (L.racemosa), offsetting the overlogging upon the red mangrove (R. mangle), therefore proving resiliency to the ecosystem, and (c) the habitat of the blue crab that sustains an important folk catchery.
At a global scale, it is urgent to conduct regional and national assessments of mangrove ecological condition.Since we observed that fragmentation (increased edge length and decreased patch area) is an important driver of human access and mangrove ecological condition, we hypothesize that cumulative small-scale edge effects might scale up to whole patches and neighborhoods, and, ultimately, to the entire regions, as proposed for terrestrial forests [67].We recommend a thorough field assessment of understory ecological conditions and targets for conservation in coastal-scapes, where remote sensors allow identifying sharp mangrove-to-pasture edges.It is an action of utmost importance in the Caribbean region where mangroves have been encroached upon by cattle ranching or urban expansion over various decades.
As a final remark, bearing in mind that the expansion of rangelands, villages and dense settlements are a major sign of the Anthropocene, this process and the resulting spatial mosaics offer a new paradigm for studying and preserving the biosphere, and therefore imposing new responsibilities, both conceptually and methodologically [10].A change in the mindset of managers and conservation biologists is needed to accept that urban and rural landscapes, and areas outside MPAs, in general, are new scenarios for the conservation, that may harbor fairly high species richness that provide valuable ecosystem services [5,9].As the new responsibility is to sustain biodiversity and people at the same time [7], conservation ecologists need to move from studying only pristine areas and MPAs (sensu [68]), in our case, to learning scientific lessons and finding conservation opportunities for mangroves in anthropogenic coastal-scapes.

Figure 1 .
Figure 1.Location of the Southeast coast of the Urabá Gulf (Colombian Caribbean), Turbo Municipality and various mangrove neighborhoods.The urban boundary of Turbo City is shown.The study was conducted in Punta Yarumal, Punta Las Vacas and Punta Coquito.

Figure 2 .
Figure 2. Location of sampling points on mangrove neighborhoods within four land use and tenure contexts in the Southeast coast of the Urabá Gulf: (a) Rural, (b) Peri-urban, (c) Military, and (d) Wild.
Codes: Annloss: annual loss rate; UrbanE: urban edge; PastureE: pasture edge; TotalE: total edge; AnthroE: anthropogenic edge; TotPatch: total number of patches; PatchD: patch density; MeanArea: mean patch area; TotArea: total neighborhood area; ADI: Cumulative Anthropogenic Disturbance Index; RmStumpD: Rhizophora mangle stump mean diameter; d: mean tree diameter for R. mangle, L. racemosa and A. germinans; D: mean density for each species; IVI: mean Importance Value Index for each species (including Pelliciera rhizophorae); Snails: mean percent presence for Neritina virginea; Stress: nearly zero when single values for each context were included, due to the reduced number of rows; 0.090 when the 31 individual transect values for ecological condition variables within each context were included.The ordination with the first data set is shown for clarity due to the overlap of many transects within the same context.

Figure 4 .
Figure 4. Non-Metric Multi-Dimensional Scaling of mangrove neighborhoods in different contexts in 2011 (Black letters).The ordination was based on 22 variables (Red letters) describing fragmentation (Tables1 and 2) and ecological condition (Tables3 and 4, in addition to transect-based ADI data pooled in Figure5but shown in Figures6 and 7).

Table 3 .
Number and diameter (mean and range) of stumps in three common mangrove species in mangrove neighborhoods in different contexts in 2011..:Not available.T: trampling; L: logging; W: waste; S: human structures.See the dispersion of ADI values across the individual transects in each context in Figures 6 and 7.

Figure 5 .
Figure 5. Cumulative scores of the Anthropogenic Disturbance Index (ADI) in mangrove neighborhoods in different contexts in 2011.

Figure 6 .
Figure 6.Pearson's correlation and linear regression parameters between mangrove fragmentation variables, ADI and vegetation structure parameters in mangrove neighborhoods in Turbo Municipality in 2011.Regressions of ADI relative to mean patch area and density, and of R. mangle IVI, and IVI, density and DBH of L. racemosa relative to ADI are shown.

Figure 7 .
Figure 7. Neritina virginea presence in mangrove neighborhoods in different contexts in Turbo Municipality

Table 1 .
[48]rove area dynamics relative to context between 1938 and 2009.The annual rate of change (r) was calculated according to[48].

Table 2 .
Patch-based metrics for mangrove neighborhoods within different contexts in the Turbo Municipality in 2009.

Table 4 .
Mangrove structure variables for neighborhoods in different contexts in 2011.Mean and standard deviation values are shown.Pelliciera rhizophorae.Note the presence of the latter species only in Military and Wild contexts.The IVI was expressed as the percent of the maximum potential number (300), based on the maximum values of relative density, dominance and frequency.See number of replicates (sampling points and transects) in text.
*: A deviation was not computed because a single tree was found.