A 1900 Year Sediment Record Suggests Recent Establishment of Black Mangrove ( Avicennia germinans ) Stands within a Salt Marsh in St. Augustine, Florida, USA

: Salt marshes and mangroves are currently being affected by rising temperatures. Mangroves thrive below − 29 ◦ N latitude in Florida, USA, and have a low tolerance for extreme cold events, whereas salt marshes dominate further north. One potential effect of climate change is a reduction in the frequency of extreme cold events, which may lead to mangrove expansion into salt marsh systems. Our research identiﬁed sediment proxy indicators of salt marsh and mangrove environments. These indicators were applied to soil cores from intertidal wetlands near the current northern limit of mangrove presence on the east coast of Florida, to determine if mangrove expansion into salt marsh environments has precedence in the deeper past. Our ﬁndings suggest that mangrove and salt marsh sediments can be distinguished using a combination of stable carbon isotope ratios of sedimentary organic matter and macroscopic plant fragments, and our results showed that a mangrove stand that we cored established only recently. This result is consistent with other work in the southeastern United States that suggests that mangroves established at the current boreal limit only recently after the end of the Little Ice Age, and that the current mangrove expansion may be fueled by anthropogenic climate change. Department of Earth and Stable Isotope Laboratory. Powdered core samples were measured for total organic carbon, nitrogen, and carbon isotope composition ( δ 13 C) using a Costech Elemental Analyzer coupled to a Thermo Delta V Plus Mass Spectrometer (Thermo Scientiﬁc, Waltham, MA, USA) at the University of North Carolina Wilmington’s Center for Marine Science Isotope Ratio Mass Spectrometry Core Facility. Carbonates were removed from both surface and core samples through HCl pretreatment prior to analysis. All δ 13 C values were reported relative to the V-PDB standard. Multiple samples from all three coring locations were sent to International Chemical Analysis for radiocarbon dating of the bulk peat. Indicated sample depths represent the mid-point of the sediment sections sent for analysis. Radiocarbon analysis was performed by following standard procedures for organic sediments, and the calibrations were calculated using CALIB Rev. 7.1.0 (Stuiver, Reimer, and Reimer, Belfast, UK) and the IntCal13 calibration database [32]. These standard procedures included examination of the samples and manual removal of contaminants like rootlets during the physical pretreatment. Addi-tionally, samples underwent Acid–Alkali–Acid pretreatment, which included acid (HCl) treatment to remove acid soluble compounds and secondary carbonates, as well as a base (NaOH) treatment to remove humic acids, and a ﬁnal acid treatment (HCl) to eliminate atmospheric CO 2 . Linear interpolation was used to assign ages between the dated horizons, and the age /depth models were based on the mid-point of the sediment section used for radiocarbon analysis and the median, 2 σ age in cal BP (present = 1950). Carbon sequestration rates were determined for different core segments by multiplying the sedimentation rate of that core segment (determined through the radiocarbon dating) by the average carbon content for that segment and the average bulk density for that core segment. This was done only for core sections where the paleoenvironmental interpretation salt marsh or mangrove) was robust.


Introduction
Current climate change, caused largely by anthropogenic increases in greenhouse gas concentration, is affecting Earth's temperatures, weather patterns, and ecosystems. Increasing levels of greenhouse gasses in the atmosphere have led to increases in global temperatures and decreases in extreme cold events [1]. These climate changes have resulted in impacts on different ecosystems and the services they provide.
Coastal saline wetlands in Florida, USA, contain a diverse mosaic of woody mangroves and herbaceous saltmarsh vegetative communities. On a global basis, mangrove communities tend to dominate in lower latitudes, generally below −29 • N, while salt marshes dominate higher latitude settings [2]. Physiological tolerance to cold is a primary control on the position of the ecotone between mangroves and saltmarshes, with cold-intolerant mangrove communities largely being restricted to warmer, tropical, and sub-tropical latitudes [3,4]. At a local level, the mangrove/saltmarsh ecotone position may also result from elevation changes, hurricanes, increased predation, or rising sea levels [5][6][7][8][9]. Cavanaugh et al. [10] found that a temperature-related ecological threshold of −4 • C provides a control on mangrove distribution and favors stable-state dominance by herbaceous saltmarsh ecosystems in northern Florida. Osland et al. [11] found that while leaf damage occurs close to −4 • C for Avicennia germinans (black mangrove, the northern-most species in Florida), mortality occurs near −7 • C. However, Cavanaugh et al. [10] also found, through analysis of satellite imagery, substantial increases in mangrove area and the northward

Materials and Methods
Surface soil samples were collected in May of 2017 at Ponce Preserve Park in the town of Ponce Inlet and at the Marine Discovery Center in New Smyrna Beach, Florida (see Figure 1). Both of these locations occur within the historic transition zone between salt marsh and mangrove environments and currently contain both ecosystem types. We took surface samples from Spartina alterniflora (Salt Marsh Cordgrass) dominated salt marshes at Ponce Inlet Park, where the marsh vegetation was extensive and mostly unmixed with mangroves, and mangrove surface samples from both Ponce Inlet Park and at the Marine Discovery Center, where there were extensive, well-developed mangrove systems dominated by black mangrove (see Supplementary). Rhizophora mangle (red mangrove) and Laguncularia racemose (white mangrove) were also present in minor abundance near the black mangrove surface sampling locations [30].
Three soil cores were removed from black mangrove stands in St. Augustine, Florida on 23 March 2018, using a half-spoon, 1 m length, and 4 cm diameter Eijkelkamp push corer. Boat reconnaissance along the Intracoastal Waterway north of St. Augustine revealed that our core sampling locations were at the northern limit of sizable mangrove-dominated stands in the area, although individual black mangroves and small clumps of black mangroves were found north of our coring locations, and other researchers have reported extant black mangrove stands as far north as Fort George Inlet, near Jacksonville, FL [16]. Two of the soil cores were taken at Nease Park near St. Augustine, Florida (29 • 55 35 N, 81 • 17 51 W) within a stand of black mangroves that was bounded on one side by Spartina alterniflora or Batis maritima (saltwort)-dominated salt marsh and on the other side by upland (Figure 1). Another core was taken within a black mangrove stand surrounded by Spartina alterniflora-dominated salt marsh near the Vilano Boat Ramp, St. Augustine (29 • 54 43 N, 81 • 18 22 W; Figure 1). Surface soil samples from the Marine Discovery Center were collected from random points within the mangrove stands and were collected along transects at Ponce Preserve, through salt marsh and mangrove separately. A total of 30 surface samples were collected from mangrove stands, and an additional 30 samples were collected from the salt marsh by either sub-sampling from the top of shallow soil cores using a modified 1.5 cm diameter plastic syringe with volume indicators, or by using the plastic syringe to sample the surface sediments directly. In this way sampled sediment volume was recorded for later bulk density measurements. Surface samples were described in the field for Munsell color, texture (using the texture-by-feel method, [31]), and macroscopic sediment characteristics, such as the presence of woody versus grass fragments. Samples from each location were returned to the lab for further, immediate processing. Three soil cores were removed from black mangrove stands in St. Augustine, Florida on 23 March 2018, using a half-spoon, 1 m length, and 4 cm diameter Eijkelkamp push corer. Boat reconnaissance along the Intracoastal Waterway north of St. Augustine revealed that our core sampling locations were at the northern limit of sizable mangrovedominated stands in the area, although individual black mangroves and small clumps of black mangroves were found north of our coring locations, and other researchers have reported extant black mangrove stands as far north as Fort George Inlet, near Jacksonville, FL [16]. Two of the soil cores were taken at Nease Park near St. Augustine, Florida (29°55'35" N, 81°17'51" W) within a stand of black mangroves that was bounded on one side by Spartina alterniflora or Batis maritima (saltwort)-dominated salt marsh and on the other side by upland (Figure 1). Another core was taken within a black mangrove stand surrounded by Spartina alterniflora-dominated salt marsh near the Vilano Boat Ramp, St. Augustine (29°54'43" N, 81°18'22" W; Figure 1).
Surface soil samples from the Marine Discovery Center were collected from random points within the mangrove stands and were collected along transects at Ponce Preserve, through salt marsh and mangrove separately. A total of 30 surface samples were collected from mangrove stands, and an additional 30 samples were collected from the salt marsh by either sub-sampling from the top of shallow soil cores using a modified 1.5 cm diameter The core samples were transported to the Environmental Science lab at Stetson University, where they were processed further. Modified plastic syringes were used to sub-sample the cores at 3.5 to 5 cm intervals based on total core length, with shorter cores having a tighter sub-sampling interval. Sub-sample volumes were recorded, and then samples were transferred to aluminum dishes and weighed. Samples were then dried to a constant weight at 65 • C. Dry weight was recorded, and bulk density was then calculated by dividing sample weight by volume. Samples were then powdered using a mortar and pestle. The samples were split, with half of the sediment fraction sent for bulk C, N, and δ 13 C analysis and the other half retained for loss on ignition measurements. Surface samples were processed and split in the same way. Samples were measured for loss on ignition through combustion at 550 • C for 2 h.
Samples sent for bulk C, N, δ 13 C, and also for radiocarbon analysis were processed and analyzed in a similar way as described in Tanner et al. [27] and as summarized again here. Powdered core samples were measured for total organic carbon, nitrogen, and carbon isotope composition (δ 13 C) using a Costech Elemental Analyzer (Costech Analytical Technologies Inc., Valencia, CA, USA) coupled to a Thermo-Finnigan Delta+ XL Mass Spectrometer (Thermo Finnigan LLC, San Jose, CA, USA) at the University of Tennessee's Department of Earth and Planetary Sciences Stable Isotope Laboratory. Powdered core samples were measured for total organic carbon, nitrogen, and carbon isotope composition (δ 13 C) using a Costech Elemental Analyzer coupled to a Thermo Delta V Plus Mass Spectrometer (Thermo Scientific, Waltham, MA, USA) at the University of North Carolina Wilmington's Center for Marine Science Isotope Ratio Mass Spectrometry Core Facility. Carbonates were removed from both surface and core samples through HCl pretreatment prior to analysis. All δ 13 C values were reported relative to the V-PDB standard.
Multiple samples from all three coring locations were sent to International Chemical Analysis for radiocarbon dating of the bulk peat. Indicated sample depths represent the mid-point of the sediment sections sent for analysis. Radiocarbon analysis was performed by following standard procedures for organic sediments, and the calibrations were calculated using CALIB Rev. 7.1.0 (Stuiver, Reimer, and Reimer, Belfast, UK) and the IntCal13 calibration database [32]. These standard procedures included examination of the samples and manual removal of contaminants like rootlets during the physical pretreatment. Additionally, samples underwent Acid-Alkali-Acid pretreatment, which included acid (HCl) treatment to remove acid soluble compounds and secondary carbonates, as well as a base (NaOH) treatment to remove humic acids, and a final acid treatment (HCl) to eliminate atmospheric CO 2 . Linear interpolation was used to assign ages between the dated horizons, and the age /depth models were based on the mid-point of the sediment section used for radiocarbon analysis and the median, 2σ age in cal BP (present = 1950).
Carbon sequestration rates were determined for different core segments by multiplying the sedimentation rate of that core segment (determined through the radiocarbon dating) by the average carbon content for that segment and the average bulk density for that core segment. This was done only for core sections where the paleoenvironmental interpretation (either salt marsh or mangrove) was robust.
Microsoft Excel v. 2013 software with the data analysis tool pack (Microsoft Corporation, Redmond, WA, USA) was used to test for significant differences in independent variables (water content, loss on ignition, bulk density, carbon content, C/N, and δ 13 C) between mangrove and salt marsh groups for the sampled surface sediments. The same software was also used to test for correlation between the measured variables using the coefficient of determination (r 2 ). We performed multiple t-tests (homoscedastic, two-tailed distribution) to determine significant differences in means of the measured variables between the two groups (marsh and mangrove) at a significance level of α = 0.01. A Bonferroni Correction was applied to the alpha value as a conservative method to counteract the problem of multiple comparisons by dividing that value (0.01) by the number of variables tested (6). Therefore, a significance of α = 0.01 was achieved by assessing each comparison at α = 0.0017.

Results
Thirty surface samples were recovered from Spartina alterniflora (with minor Batis maritima and lesser abundance of other species) dominated salt marsh and 30 surface samples were also recovered from mangrove (black mangrove dominated with minor red mangrove and Batis maritima and lesser abundance of other species) at Ponce Inlet, FL. Texture-by-feel grain size analysis revealed clay-rich sediments in both salt marsh and mangrove environments. Macroscopic description of surface sediments showed that the deposits can be visually distinguished by the abundant presence of grass fragments in salt marsh sediments and the presence of woody material in mangrove sediments ( Figure 2). maritima and lesser abundance of other species) dominated salt marsh and 30 surface samples were also recovered from mangrove (black mangrove dominated with minor red mangrove and Batis maritima and lesser abundance of other species) at Ponce Inlet, FL. Texture-by-feel grain size analysis revealed clay-rich sediments in both salt marsh and mangrove environments. Macroscopic description of surface sediments showed that the deposits can be visually distinguished by the abundant presence of grass fragments in salt marsh sediments and the presence of woody material in mangrove sediments ( Figure 2).

Figure 2.
The top photograph (a) shows a soil core from a salt marsh site with abundant grass fragments. The bottom photograph (b) shows a soil core from a mangrove site with fibrous roots, sticks, and woody fragments. Both of these images include the sediment surface (right side of each).
The surface samples from the mangrove sites in Volusia County, Florida, contained primarily silty clay (30% of samples), silty loam (27%), sandy clay (17%), and clay (13%) and consisted of fibrous roots, sticks and woody fragments ( Figure 2). The surface samples from the salt marsh sites contained primarily silty clay (47% of samples), clay (33%), and silty loam (10%). The soil cores from St. Augustine, Florida, had similar texture and plant fragment content to the surface samples ( Figure 3). In all three soil cores, woody fragments were noted in the upper sections of the cores and only grass fragments were found down the core. The surface samples from the mangrove sites in Volusia County, Florida, contained primarily silty clay (30% of samples), silty loam (27%), sandy clay (17%), and clay (13%) and consisted of fibrous roots, sticks and woody fragments ( Figure 2). The surface samples from the salt marsh sites contained primarily silty clay (47% of samples), clay (33%), and silty loam (10%). The soil cores from St. Augustine, Florida, had similar texture and plant fragment content to the surface samples ( Figure 3). In all three soil cores, woody fragments were noted in the upper sections of the cores and only grass fragments were found down the core. Table 1 shows values for measured surface sediment physical and chemical characteristics, including water content, loss on ignition, bulk density, percent carbon, C/N, and δ 13 C of organic matter. Compared to mangrove surface sediments, salt marsh surface sediments had a significantly higher (p < 0.01) water content (average 79.3% to 58.6%), loss on ignition (average 42.2% to 24.6%), carbon percentage (average 20.2% to 10.5%) and δ 13 C value (average −19.3‰ to −23.91‰). Conversely, mangrove sediments showed a significantly higher (p < 0.01) bulk density (average 0.55 g/cm 3 to 0.23 g/cm 3 ). C/N showed no significant difference between salt marsh and mangrove surface sediments (p = 0.14). Loss on ignition values and percent carbon values were highly correlated (r 2 = 0.67), as would be expected (Table 2). Therefore, the percent carbon values are presented in Figure 3 without the loss on ignition values in order to avoid redundancy and to save space. Those measurements were also correlated with soil water content and bulk density (Table 2), as would be expected. Soils with a higher carbon content had a higher loss on ignition, higher water content, and lower bulk density. Black mangrove surface sediment samples showed a range of δ 13      The soil cores from Nease Park were collected to refusal at 76 cm depth (core 1) and 70 cm depth (core 2). The soil core from the Vilano Boat Launch was collected to refusal at 205 cm depth. Radiocarbon ages for all cores are presented in Table 3 along with supporting information. The upper sample from Nease Park, core 1, and the sole sample from core 2 both returned modern ages. The deeper 66.5 cm sample depth from core 2 returned a median age of 737 cal BP. We did not attempt to date a deeper sample from core 2 because of a lack of datable organic material. Given that only one sample returned an age across both Nease Park cores, an age/depth model was not attempted, and all data are presented vs. depth for that site. Four radiocarbon ages were obtained for the Vilano Boat Launch core due to its greater length. All radiocarbon ages for the Vilano core were in stratigraphic order, and the median ages ranged from 604 cal BP at 52.5 cm depth to 1828 cal BP at 192 cm (Table 3). An age/depth model for that core was constructed using linear interpolation (Figure 4), and the age depth model was applied to the proxy data from the core (Figure 3). For the linear interpolation, the surface sample (1 cm) was assumed to represent the time of core collection (−78 cal BP). Table 3. Radiocarbon data are presented for the three cores. Sample depths are in cm below the ground surface and represent the mid-point of the sediment section sent for dating. The calibration procedures are described in the methods section. The soil cores from Nease Park and the Vilano Boat Launch all showed a general decrease in organic carbon content and water content with depth, while bulk density generally increased with depth (Figure 3), as would be expected with temporal decomposition and compaction of organic matter. C/N varied between 12.5 and 21.8 among the three soil cores, showing no clear trends. The δ 13 C values were almost all in the range of overlap between mangrove and salt marsh sediments for both Nease Park cores. For the Vilano Boat Launch core, δ 13 C values were lower just at the surface (within the range of overlapping mangrove and salt marsh sediment values) and were higher with depth (in the range of salt marsh sediment values), in line with the trend of woody and grassy fragments recorded for that core.

Conventional
The data from the well-dated Vilano Boat Launch core allowed for the calculation of carbon sequestration rates for both the mangrove-dominated upper sediments and the remainder of the core, dominated by salt marsh (Table 4). Data were divided into sections representing interpreted mangrove vs. salt marsh deposition (based on Figure 3) and then by differences in sedimentation rate between bracketed radiocarbon ages (see Table 3). Carbon sequestration rates were generally similar between the upper mangrove sediments (average = 20.01 g C m −2 yr −1 ) and the underlying, older salt marsh sediments (total average for all depth intervals = 21.53 g C m −2 yr −1 ), with a relatively persistent sequestration rate over long time periods.  The soil cores from Nease Park and the Vilano Boat Launch all showed a general decrease in organic carbon content and water content with depth, while bulk density generally increased with depth (Figure 3), as would be expected with temporal decomposition and compaction of organic matter. C/N varied between 12.5 and 21.8 among the three soil cores, showing no clear trends. The δ 13 C values were almost all in the range of overlap between mangrove and salt marsh sediments for both Nease Park cores. For the Vilano Boat Launch core, δ 13 C values were lower just at the surface (within the range of overlapping mangrove and salt marsh sediment values) and were higher with depth (in the range of salt marsh sediment values), in line with the trend of woody and grassy fragments recorded for that core.
The data from the well-dated Vilano Boat Launch core allowed for the calculation of carbon sequestration rates for both the mangrove-dominated upper sediments and the remainder of the core, dominated by salt marsh (Table 4). Data were divided into sections representing interpreted mangrove vs. salt marsh deposition (based on Figure 3) and then by differences in sedimentation rate between bracketed radiocarbon ages (see Table 3). Carbon sequestration rates were generally similar between the upper mangrove sediments (average = 20.01 g C m −2 yr −1 ) and the underlying, older salt marsh sediments (total average for all depth intervals = 21.53 g C m −2 yr −1 ), with a relatively persistent sequestration rate over long time periods. Figure 4. The radiocarbon age/depth curve is plotted for the Vilano Boat Launch core, where multiple radiocarbon dates were obtained. The data points represent the mid-point of the sample interval sent for radiocarbon age determination and the calibrated, 2σ median age in cal BP. The solid line connecting the points represents the linear interpolation between dated horizons. The dashed lines indicate the calibrated 2σ age ranges for the radiocarbon ages. The surface sample is assigned an age of −78 cal yr BP given that the samples were collected in 2018 and present equals 1950 by convention. Table 4. Carbon sequestration data for the Vilano Boat Launch core are presented for depth intervals based on interpreted wetland environment (mangrove vs. salt marsh) and by differences in sedimentation rate between bracketed radiocarbon ages. The presented carbon percentages, sedimentation rates, and bulk densities represent averages for those sample depth intervals.

Discussion
Surface sediment samples collected from salt marsh and mangrove environments showed significant differences (p < 0.01) in water content, loss on ignition, bulk density, carbon percentage, and δ 13 C (Table 1). In addition, salt marsh surface sediments contained abundant grassy fragments, whereas mangrove surface sediments contained woody material ( Figure 2). Therefore, water content, loss on ignition, bulk density, carbon percentage, and δ 13 C all have potential as proxy indicators to distinguish between salt marsh and mangrove sediments in core deposits, and several of these indicators have been used in other locations in this regard [17,18,33]. Our data suggest that C/N is not a metric that is useful to distinguish between salt marsh and mangrove sediments for our study site (Table 1; p = 0.14), although others working in the region have found that mangrove deposits have a lower C/N, possibly due to additional algal inputs related to increased tidal inundation [22]. Recorded sediment textures were similar between salt marsh and mangrove surface sample sites, suggesting that soil texture cannot be effectively applied to core samples to distinguish salt marsh and mangrove environments for our study sites.
Of the measured proxy indicators that were potentially useful to distinguish salt marsh from mangrove sediments (water content, loss on ignition, bulk density, carbon percentage, δ 13 C, and plant fragments), four logically change with depth. Loss on ignition and carbon percentages relate to soil bulk density and tend to decline with depth due to decomposition and soil compaction [34]. Our data showed that water content and bulk density are strongly correlated ( Table 2) with an inverse relationship (data from Table 1). Declines of organic carbon and bulk density with depth and increases in water content with depth were seen in our core samples (Figure 3). Given that these four potential proxies were strongly related to depth down the soil profile as opposed to marsh vs. mangrove sediment input, they have limited utility in allowing discrimination between salt marsh and mangrove sediments in core deposits. Therefore, the best proxy indicators that we employed for distinguishing mangrove and salt marsh deposits included a combination of δ 13 C values of sedimentary organic matter and plant fragment information.
Radiocarbon age determinations were obtained for bulk sediments for the Nease Park and Vilano Boat Launch cores. Plant macrofossils represent an ideal material used for radiocarbon dating because their carbon source is known and does not represent an admixture of multiple carbon sources [35]. However, we were not able to identify aboveground mangrove plant macrofossils that would be suitable for radiocarbon dating in the sediments. Only woody material was found in the sediments, which could have represented either above-ground timber or below-ground roots that are not suitable for radiocarbon age control [36]. Salt marsh grass fragments identified in the sediment deposits were fragile and we were not able to recover them intact. Bulk sediment is commonly used for radiocarbon chronologies in mangrove sediments [33,36]. Studies have shown that organic matter in both mangrove and salt marsh sediments is generally reflective of autochthonous sources [22,[36][37][38][39]. In mangrove environments, the primary carbon source to the sediments is often represented by the mangrove plants themselves, followed by macro-and micro-algae, and finally transport and deposition from upstream materials and materials from the adjacent coastal zone [37]. A consistent age-depth model for the Vilano Boat Launch core with in-sequence radiocarbon ages ( Figure 4) and reasonable calculated sedimentation rates (Table 4) suggested that the bulk radiocarbon ages for that core were robust and reflective of true ages of deposition for that core.
The two soil cores from Nease Park were relatively short (each <80 cm long), and only one non-modern radiocarbon age was returned, representing the deeper sediment sample from core 2 ( Table 3). The lack of non-modern radiocarbon ages confounded the development of age/depth models for core 2. However, plant fragment data indicated the presence of a mangrove environment underlain by salt marsh deposits for both cores (Figure 3). Carbon isotope data for these two cores were ambiguous, however, as δ 13 C values for both cores were mostly in the field of overlapping salt marsh and mangrove values, albeit towards the mangrove side of the overlapping range for both cores, and especially for core 2. The paleoenvironmental interpretation was therefore difficult for both of the Nease Park cores, given that plant fragment information suggested a mangrove to salt marsh transition for each core with depth, whereas δ 13 C values were more suggestive of a persistent mangrove environment (or at the very least, a mixed salt marsh and mangrove environment) for both cores.
The interpretation of the Vilano Boat Launch core was facilitated by proxy indicators that were in good agreement ( Figure 3). Also, all four radiocarbon dates obtained for the core were in stratigraphic order, and the age/depth model obtained through linear interpolation suggested a relatively constant sedimentation rate for the past −1900 cal BP at the core location (Figure 4), without any obvious depositional hiatuses. The core was collected within a black mangrove stand that was surrounded by Spartina alterniflora-dominated salt marsh. The upper part of the core reflected deposition in a mangrove environment, with woody fragments and lower δ 13 C values observed in sediments representing the past 100+ years (Figure 3). Only grass fragments were found deeper in the core, and the shift to higher δ 13 C values for the remainder of the core, which were all in the range of C4 plant-dominated Spartina alterniflora salt marsh [24], suggested insignificant input from mangroves to the sediment pool for the remainder of the core. It was thus apparent that this core sampling location transitioned to mangrove only recently, after almost 1800 years of salt marsh deposition. It should be noted that coarse sediments (sand-dominated) were observed towards the middle and bottom parts of the core (Figure 3), indicating higherenergy deposition. This was perhaps due to closer proximity of a migrating tidal channel to the core site during times of coarse sediment deposition, or perhaps deposition from a high-energy event like a tropical cyclone, which would both increase the coarse sediment fraction and the deposition rate [40]. If the coarse sediment fraction was driven by a rapid event, such as a tropical, a high sedimentation rate for that depth interval would be expected. However, the sedimentation rate for the Vilano Boat Launch core was not highest where the sand layers were present and was instead higher in between the sand layers ( Figure 3, Table 4).
The results from the well-dated Vilano Boat Launch core suggested a persistent salt marsh environment from approximately 1900 cal BP to approximately 50 cal BP, when significant mangrove sedimentation was then indicated to the present. A nearby study based on historical records and a climate-related "mangrove suitability" index suggested that mangrove stands were present in the region at least as far back as 1860 AD [16]. Our results do not contradict this record given the uncertainties inherent in radiocarbon dating and in the development of age/depth models. Our results were also consistent with research from coastal Louisiana, where researchers suggested the development of mangrove stands at the current boreal limit only after the Little Ice Age, ending between 1850 and 1890 AD [17].
Our data allowed for the calculation of carbon sequestration rates for the Vilano Boat Launch core (Table 4). These measured rates were similar to those reported from a nearby site for mangrove sediments and were higher than those reported for nearby salt marsh sediments [22]. They were lower by about one order of magnitude than average global rates of carbon storage for both salt marsh (218 ± 24 g C m −2 yr −1 ) and mangrove (226 ± 31 g C m −2 yr −1 ) systems, although there is considerable variability in carbon storage measurements across the globe [41]. Our carbon sequestration rate results, although limited, suggested that ecosystem services related to the efficiency of carbon burial are similar between salt marsh and mangrove systems in the area of our study over longer time periods (decadal to century scale). The higher percentage of carbon measured for the mangrove-dominated sediment was offset by the higher bulk density measured for the salt marsh sediments. Vaughn et al. [22] suggested that recalcitrant lignin associated with woody mangrove sediments allows for the potential for higher long-term carbon sequestration rates given slower potential decomposition. Our data suggest that salt marsh sediments can also store carbon efficiently over longer timescales in the region of our study.

Conclusions
Our results from surficial sediments suggest that mangrove sediments in northern Florida can be distinguished from C4 plant-dominated salt marsh sediments using a combination of stable carbon isotope ratios of sedimentary organic matter combined with macroscopic plant fragment material (woody vs. grassy). These indicators were applied to three soil cores recovered from near the modern black mangrove boreal limit in northern Florida, from within stands of black mangrove. Paleoenvironmental interpretation for two of the soil cores (from Nease Park) was ambiguous, and the short sediment record preserved at that site confounded the development of an age/depth model. The paleoenvironmental proxy record developed from a core from the other site (Vilano Boat Launch) suggested C4 plant-dominated salt marsh deposition from the bottom of the core to around 50 cal BP, where there was a transition to mangrove-dominated sediments that persisted to the present time. This finding was consistent with work by others that suggests mangrove presence at its current boreal limit in the southeastern United States only within the past several hundred years and is also consistent with the hypothesis that the current observed mangrove expansion in the region is due to anthropogenic climate change.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/quat5010002/s1, Figure S1: photograph taken at Ponce Preserve Park showing mangrove and salt marsh, Figure S2: photograph of core collection at Nease Park, Figure S3: photograph of Vilano Boat Launch core showing sampling scheme.
Author Contributions: B.T. and J.E. conceptualized the study and were also involved with the study methodology, formal analysis, investigation, writing-original draft preparation, writing-review and editing, project administration, and funding acquisition. J.C. and C.S. conducted formal analysis, investigation, writing-original draft preparation, and writing-review and editing. All authors have read and agreed to the published version of the manuscript.