Reviews and syntheses : 210 Pb-derived sediment and carbon accumulation rates in vegetated coastal ecosystems – setting the record straight

Vegetated coastal ecosystems, including tidal marshes, mangroves and seagrass meadows, are being increasingly assessed in terms of their potential for carbon dioxide sequestration worldwide. However, there is a paucity of studies that have effectively estimated the accumulation rates of sediment organic carbon (Corg), also termed blue carbon, beyond the mere quantification of Corg stocks. Here, we discuss the use of the 210Pb dating technique to determine the rate of Corg accumulation in these habitats. We review the most widely used 210Pb dating models to assess their limitations in these ecosystems, often composed of heterogeneous sediments with varying inputs of organic material, that are disturbed by natural and anthropogenic processes resulting in sediment mixing and changes in sedimentation rates or erosion. Through a range of simulations, we consider the most relevant processes that impact the 210Pb records in vegetated coastal ecosystems and evaluate how anomalies in 210Pb specific activity profiles affect sediment and Corg accumulation rates. Our results show that the discrepancy in sediment and derived Corg accumulation rates between anomalous and ideal 210Pb profiles is within 20 % if the process causing such anomalies is well understood. While these discrepancies might be acceptable for the determination of mean sediment and Corg accumulation rates over the last century, they may not always provide a reliable geochronology or historical reconstruction. Reliable estimates of Corg accumulation rates might be difficult at sites with slow sedimentation, intense mixing and/or that are affected by multiple sedimentary processes. Additional tracers or geochemical, ecological or historical data need to be used to validate the 210Pbderived results. The framework provided in this study can be instrumental in reducing the uncertainties associated with estimates of Corg accumulation rates in vegetated coastal sediments. Published by Copernicus Publications on behalf of the European Geosciences Union. 6792 A. Arias-Ortiz et al.: Reviews and syntheses: 210Pb-derived sediment and carbon accumulation rates

Abstract. Vegetated coastal ecosystems, including tidal marshes, mangroves and seagrass meadows, are being increasingly assessed in terms of their potential for carbon dioxide sequestration worldwide. However, there is a paucity of studies that have effectively estimated the accumulation rates of sediment organic carbon (C org ), also termed blue carbon, beyond the mere quantification of C org stocks. Here, we discuss the use of the 210 Pb dating technique to determine the rate of C org accumulation in these habitats. We review the most widely used 210 Pb dating models to assess their limitations in these ecosystems, often composed of heterogeneous sediments with varying inputs of organic material, that are disturbed by natural and anthropogenic processes resulting in sediment mixing and changes in sedimentation rates or erosion. Through a range of simulations, we consider the most relevant processes that impact the 210 Pb records in vegetated coastal ecosystems and evaluate how anomalies in 210 Pb specific activity profiles affect sediment and C org accumulation rates. Our results show that the discrepancy in sediment and derived C org accumulation rates between anomalous and ideal 210 Pb profiles is within 20 % if the process causing such anomalies is well understood. While these discrepancies might be acceptable for the determination of mean sediment and C org accumulation rates over the last century, they may not always provide a reliable geochronology or historical reconstruction. Reliable estimates of C org accumulation rates might be difficult at sites with slow sedimentation, intense mixing and/or that are affected by multiple sedimentary processes. Additional tracers or geochemical, ecological or historical data need to be used to validate the 210 Pbderived results. The framework provided in this study can be instrumental in reducing the uncertainties associated with estimates of C org accumulation rates in vegetated coastal sediments.

Introduction
Recognition of the globally significant role of vegetated coastal habitats, including tidal marsh, mangrove and seagrass, as sinks of carbon dioxide (CO 2 ) (Duarte et al., 2013) has generated interest in the potential to mitigate CO 2 emissions through the management of these ecosystems, an approach termed the "blue carbon" strategy (Duarte et al., 2013;Mcleod et al., 2011;Nellemann et al., 2009). However, efforts to include vegetated coastal ecosystems in existing carbon mitigation strategies have met with an important limitation: there is a paucity of estimates of C org sequestration rates, particularly in seagrass habitats (Johannessen andMacdonald, 2016, 2018;Macreadie et al., 2018).
Two interrelated measurements of importance to the rate of C org sequestration are the sediment C org content and the sediment accumulation rate. To date, most research has focused on the first term, which provides information about the C org stock sequestered in sediments (Howard et al., 2014;Pendleton et al., 2012). However, C org stocks alone cannot be used to fully assess the C org storage capacity or to establish comparisons among sites. Measurements of C org accumulation rates (CARs) address the question of how much C org is sequestered in a specified time period and quantify the ongoing sink capacity. In general, CAR is obtained in one of two ways: by measuring the concentration of C org in sediments and ascribing dates to either the entire profile of interest or to specific intervals; or by estimating sediment accumulation rates. Determination of mean CAR is partially dependent on the timescale of interest and the dating methods used. 210 Pb has been shown to be an ideal tracer for dating aquatic sediments deposited during the last ca. 100 years, providing a time frame compatible with management actions (Marland et al., 2001) and enabling the determination of CAR and its changes with time due to natural or human impacts. Due to the relatively long integration period (decades to a century), mean 210 Pb-derived CAR estimates are not affected by interannual variability, hence allowing for the assessment of shifts in C org accumulation from the "baseline" condition (such as the C org that naturally cycles through an ecosystem; Howard et al., 2017). Although several review papers have elaborated the applications of 210 Pb as a tracer in lacustrine and marine environments (Appleby, 2001;Baskaran et al., 2014;Du et al., 2012;Kirchner and Ehlers, 1998;Mabit et al., 2014;Sanchez-Cabeza and Ruiz-Fernández, 2012;Smith, 2001), little attention has been paid to the potential limitations of the 210 Pb dating method in vegetated coastal sediments. Experience shows that vegetated coastal environments often prove to be more challenging than lake or marine sediments (Saderne et al., 2018).
Vegetated coastal ecosystems may act as closed systems in which sediment accumulation is mainly associated with the buildup of autochthonous organic and inorganic material (McKee, 2011). In this situation, 210 Pb is deposited primarily from atmospheric fallout at steady state, with no postdeposi-tional mobility except for physical or biological mixing of the sediments (e.g., Alongi et al., 2004;Cochran et al., 1998;Marbà et al., 2015). In some cases, however, the process responsible for incorporating 210 Pb into the sediments might be more complex. Vegetated coastal ecosystems may receive both autochthonous and allochthonous sediments from the upstream catchment, coastal erosion, from the offshore zone during storm events (Turner et al., 2007) or in response to land use change (Mabit et al., 2014;Ruiz-Fernández and Hillaire-Marcel, 2009). Their sediments might be reworked through the action of fauna (bioturbation), tides, currents and waves as well as through boat anchoring, dredging or fishing activities (e.g., Mazarrasa et al., 2017;Sanders et al., 2014;Serrano et al., 2016a;Smoak et al., 2013). Effects associated with climate change, such as sea level rise and extreme climatic events, may also have an impact on rates of production and decomposition of organic matter (OM) and on sediment and C org accumulation (Alongi et al., 2008;Mudd et al., 2010). In such instances, sediment redistribution processes and complex accretion dynamics may violate some of the assumptions of the 210 Pb dating models, producing erratic 210 Pb specific activity profiles that are difficult to interpret.
Sediments of vegetated coastal ecosystems are known to be heterogeneous, consisting of coarse-grained deposits or bedrock covered by fine-grained sediments that settled as vegetation established Olff et al., 1997). The percentage of living (e.g., roots) and recently formed organic material is greatest in the upper 10 cm and may be affected by varying inputs of detrital sediment within vegetated coastal ecosystems and by its relative rate of decomposition. While tidal marsh and mangrove sediments have relatively high organic matter content (on average 25 %) (Breithaupt et al., 2012;Cochran et al., 1998), mineral deposits account for the majority (> 85 %) of the accumulated substrate in seagrass sediments (Koch, 2001; (Table 1). 210 Pb has a strong affinity for fine sediments (Chanton et al., 1983;Cundy and Croudace, 1995;He and Walling, 1996a) and organic matter (Wan et al., 2005), and thus any changes in these parameters due to sediment redistribution processes or to natural heterogeneity may also result in unique types of 210 Pb specific activity profiles in sediment cores of vegetated coastal ecosystems, adding complexity to the determination of sediment model age and sedimentation rates.
Here, we present how the processes of mixing and changes in the sedimentation rate, erosion, grain size heterogeneity and OM decay impact the depth distribution of 210 Pb in vegetated coastal sediments. We compare sediment and C org accumulation rates derived from disturbed and ideal 210 Pb profiles and assess the discrepancy between them. Through a set of simulations based on examples from the literature and using various 210 Pb dating models, we assess the limitations that apply to the determination of last-century C org accumulation rates in such ecosystems. Finally, we provide guidance Table 1. Common values of the main parameters for vegetated coastal sediments (seagrass, mangrove and tidal marshes): average dry bulk density (DBD), range of sedimentation rates (mean values are provided in brackets) and organic matter (OM) content, median organic carbon (C org ) contents, and decay rate of buried C org (from aboveground biomass to refractory sediment C org ).  ; mangrove (Donato et al., 2011) andtidal marsh (Craft, 2007;Hatton et al., 1983). b Seagrass (Duarte et al., 2013), mangroves (Breithaupt et al., 2012) and tidal marsh (Ouyang and Lee, 2014). c Seagrass (Koch, 2001); mangrove (Breithaupt et al., 2012); tidal marsh (Cochran et al., 1998;Ember et al., 1987). d Seagrass ; mangrove (Breithaupt et al., 2012); tidal marsh (Chmura et al., 2003). e Seagrass, mangrove and tidal marsh (Lovelock et al., 2017b).
on complementary analyses to accompany the 210 Pb dating technique that can improve sediment and derived C org accumulation rate estimates.

210 Pb dating models
210 Pb is a naturally occurring radionuclide with a half-life of 22.3 years. Particles that sink and accumulate in the bottom of aquatic systems scavenge the 210 Pb present in the water column due to the decay of 222 Rn in (1) the atmosphere and ulterior dry and wet deposition and (2) in the water. This is known as "excess 210 Pb" ( 210 Pb xs ) and is added to the "supported 210 Pb" ( 210 Pb sup ), which is continuously produced by the in situ decay of 226 Ra in bottom sediments. The accumulation of sediments over time ideally generates a decreasing distribution of 210 Pb specific activity as a function of depth (or cumulative mass in g cm −2 ) governed by the decay of the 210 Pb xs , as illustrated in Fig. 1. The 210 Pb dating models are based on the interpretation of the 210 Pb xs rate of decline in a sediment profile. Under ideal conditions, 210 Pb can accurately date sediments back to about seven halflives, i.e., about 150 years (the "dating horizon"), whereby the measurement uncertainty becomes too large to detect any 210 Pb xs . However, chronologies reaching back that far might be rarely achievable in vegetated coastal sediments as these contain relatively low concentrations of 210 Pb. The basis of the distribution of 210 Pb xs in sediments can then be described as (Koide et al., 1972) where ρ is sediment bulk density (g cm −3 ), C is the specific activity of 210 Pb xs (Bq g −1 ), z is depth below the sediment-water interface (cm), D b is the sediment mixing rate (cm 2 yr −1 ), r is the sedimentation rate (cm yr −1 ), λ is the 210 Pb decay constant (yr −1 ) and t is time (years). Commonly, depth (z) is represented as cumulative mass (m) to correct for compaction. Cumulative mass (g cm −2 ) results from the multiplication of z and ρ, and sedimentation rates are expressed as mass accumulation rates (MAR) in g cm −2 yr −1 , which can be described as MAR = ρ(v + q), where v and q are the accretion and compaction velocities, respectively (Abril, 2003) where k m is an effective mixing coefficient (g 2 cm −4 yr −1 ). The 210 Pb technique was first applied by Koide et al. (1972) to date marine sediments. Since then, a family of dating models has been used to interpret the 210 Pb xs depth distribution in marine and freshwater sediment cores, increasing in variety and complexity and involving a large diversity of postdepositional redistribution processes (Table 2). However, there are three models that are most widely used and described here: the constant flux-constant sedimentation (CF-CS) model (Krishnaswamy et al., 1971), the constant rate of supply (CRS) model  and the constant initial concentration (CIC) model (Robbins, 1978). Although these three models each have specific assumptions, they share the following: (1) the deposition of 210 Pb xs is at steady state, (2) there is no postdepositional mobility of 210 Pb and (3) the deposition of 210 Pb xs is ideal; i.e., new radioactive inputs are deposited above the previously existing material. In the simplest of the cases, the CF-CS model assumes constant 210 Pb xs depositional flux ( ) and MAR. In this case the 210 Pb specific activity at the surface (C 0 : Bq kg −1 ) is constant and decreases exponentially with cumulative mass. The depth of burial m is related to the elapsed time since burial through the rate of sedimentation (MAR) ( Table 2). If there is mixing at the surface of the core, the mean MAR can be calculated from the 210 Pb xs specific activity profile below the surface mixed layer (SML). If the specific activity of 210 Pb xs declines in sections, showing two or more exponentially decaying segments, then a mean MAR can be derived for each segment (Goldberg et al., 1977). In this way the model is to some degree able to cope with temporal variations in mass accumulation rates.
Variations in the accumulation rate may occur in response to natural processes or anthropogenic influences. Under these 6794 A.    conditions, the CRS or CIC models could be suitable. The CRS model assumes a constant flux of 210 Pb ( ) to the sediments over time. The initial specific activity is variable and inversely related to MAR (higher MAR leads to lower 210 Pb xs specific activity and vice versa). The dating is based on the comparison of 210 Pb xs inventories (I m ; Bq m −2 ) below a given depth (integration of 210 Pb xs specific activity as a function of the cumulative mass) with the overall 210 Pb xs inventory in the sediment core (I ) ( Table 2). The accurate determination of the 210 Pb xs inventories is of critical importance and required for the application of the CRS model (Appleby, 2001).
The CIC model could be a better choice at locations where sediment focusing is a major factor, where event deposit layers are present, or if significant hydrologic changes have occurred or there are hiatuses in the sediment record caused by erosion events (Appleby, 2008). The CIC model assumes that the initial 210 Pb xs specific activity at the sediment-water interface is constant with time irrespective of the sedimentation rate so that the 210 Pb xs flux covaries with MAR. This model permits an estimation of the age (t) at any depth at which 210 Pb has been measured if the initial specific activity is known (Table 2). However, the CIC model requires a monotonic decrease in 210 Pb xs specific activity down-core for age reversals to be avoided, which is rare in most vegetated coastal sediments. In that event, the calculation of mean accumulation rates alone using the CF-CS model would be a more reasonable approach, as it might be too ambitious to calculate a detailed stepwise chronology based on an often limited number of data points decreasing monotonically.
While the CIC and CF-CS models have typically been used in the marine environment, the CRS model is the most preferred in lake sediments and it is becoming widely applied in estuarine environments and vegetated coastal ecosystems (Andersen, 2017;Breithaupt et al., 2014). Some of the reasons could be that it suffers less from problems associated with non-monotonic features in the 210 Pb record and is relatively insensitive to mixing (Appleby et al., 1983;Appleby and Oldfield, 1992;Oldfield et al., 1978). The selection and use of a particular model should be based on the nature of the 210 Pb xs specific activity and MAR. For further details on the main aspects relevant to the application of 210 Pb dating models in lake or estuarine environments we recommend two detailed and comprehensive book chapters by Appleby (2001) and Andersen (2017). Here, we focus specifically on analysis of the 210 Pb dating of sediments in vegetated coastal ecosystems.

Methods
We performed a literature review of studies on sediment accumulation in vegetated coastal ecosystems in the Web of Science ™ (http://webofknowledge.com, last access: 23 August 2018) with the following keywords: mangrove sediment; salt marsh, saltmarsh or tidal marsh sediment; seagrass sediment AND (Boolean operator) 210 Pb, Pb-210 or lead-210. The search produced 86, 223 and 27 results, respectively, all of them using one or more of the three models described above, probably due to their simplicity, with some exceptions such as Klubi et al. (2017) that additionally uses the TERESA model (Table 2). From the literature review, we identified the most common sedimentary processes that result in anomalous types of 210 Pb xs specific activity profiles with depth. These could be summarized in five main processes: mixing, increasing sedimentation, erosion, changes in sediment grain size and decay of organic matter (OM). Then, we simulated the target processes on initial undisturbed seagrass, mangrove and tidal marsh sediments to determine the potential discrepancies in MAR between ideal and anomalous 210 Pb xs profiles and analyze the limitations of the 210 Pb dating technique to derive CAR in these ecosys- Table 2. Summary of the main 210 Pb-based models for sediment dating (adapted from Mabit et al., 2014).

C r
; I ref = local fallout 210 Pb inventory; C r = initial 210 Pb xs in catchment-derived sediment He and Walling (1996b) IP-CRS: initial penetration, constant rate of supply (2), initial mobility of 210 Pb xs downward; two compartments 0 to z k and z k to ∞ z m z m−1 r m C m dz; f i = 1 See reference for constants Olid et al. (2016) TERESA: time estimates from random entries of sediments and activities (1), 210 Pb xs fluxes are governed by horizontal inputs, correlation with MAR Botwe et al. (2017) (1) Non-postdepositional redistribution; (2) constant 210 Pb xs fluxes at the SWI; (3) constant MAR. C m : 210 Pb xs activity concentration in sediments at mass depth m. I: total inventory of 210 Pb xs . I m : excess inventory accumulated below depth m; k m : effective mixing coefficient (Dρ 2 ); m a : mass thickness of top sediment zone.
: Flux of 210 Pb xs onto the sediment; g: fraction of 210 Pb xs flux distributed within a certain mass depth. F i : additional supply of 210 Pb xs to layer i.
tems. Throughout this work, we refer to "deviation" in MAR or CAR as the difference between the value computed and the value estimated from the ideal 210 Pb xs profile.

Numerical simulations
All simulations started from an ideal 210 Pb xs profile, complying with all assumptions, that was then manipulated to reflect the potential effect of each process. The ideal 210 Pb xs profile was modeled considering the following: (1) a constant flux of 210 Pb xs ( ) of 120 Bq m −2 yr −1 , i.e., the average global atmospheric flux reported by Preiss et al. (1996); (2) a MAR of 0.2 g cm −2 yr −1 and dry bulk density (DBD) of 1.03 g cm −3 to represent seagrass sediments; and (3) a MAR of 0.3 g cm −2 yr −1 and DBD of 0.4 g cm −3 to represent mangrove-tidal marsh sediments based on the typical values representative of these ecosystems (Duarte et al. 2013) (Table 1). Simulated surface activity per unit area of 210 Pb xs (I 0 ; in Bq m −2 ) in ideal profiles was estimated through Eq. (3). Then Eq. (4) was applied to estimate 210 Pb xs activities per unit area along the ideal profile (Table S1 in the Supplement).
Activities of 210 Pb xs per unit area (I m ) were then converted to activities per unit mass, C m in Bq kg −1 , by dividing I m Table 3. Summary description of the numerical simulations conducted to test for the effects of sedimentary processes on 210 Pb xs specific activity profiles in seagrass and mangrove-tidal marsh sediments. MAR and CAR results derived from simulated profiles were compared with MAR and CAR estimates derived from the ideal 210 Pb xs profiles reported here; k s is the decay rate of the refractory sediment organic matter (OM) under anoxic conditions and k ox is that in oxic conditions; k lb is the decay constant of the labile OM derived from seagrass and mangrove-tidal marsh ecosystems (0.01 and 0.03 yr −1 , respectively).

Influencing
Scenario Description MAR ideal profile CAR ideal profile factor (g cm −2 yr −1 ) (g C org m −2 yr −1 ) Seagrass by the cumulative mass (m) and converting grams to kilograms at each layer. Ideal profiles were then altered to simulate the following processes or scenarios: mixing (surface and deep mixing), increasing sedimentation (by 20 %, 50 %, 100 % and 200 %), erosion (recent and past), changes in sediment grain size (coarse and heterogeneous) and OM decay (under anoxic and oxic conditions and with labile OM contribution in sediments containing 16.5 % and 65 % OM) (Table 3; refer to Appendix A for a detailed description of the methodology used to conduct each simulation). The CF-CS and CRS dating models were applied to the simulated 210 Pb xs profiles to determine the average MAR for the last century ( Table 2). The CIC model was excluded from the simulations presented in this study because in erratic 210 Pb xs profiles (1) the CRS model would lead to more reasonable approaches when the flux of 210 Pb xs is constant and (2) when that is not the case (e.g., simulations of erosion or heterogeneous grain size), the determination of mean accumulation rates alone by the CF-CS model would be a more reasonable approach. The models were applied in accordance with the simulated process. For instance, MAR was determined below the surface mixed layer in mixing simulations using the CF-CS model and piecewise in those with a change in average MAR (Appendix A). However, the models were also applied considering that (1) 210 Pb xs profiles of mixing simulations were generated by increasing MAR and vice versa, and (2) erosion was not a factor in simulated scenarios H to J. This was done to test the potential disagreement in MAR and CAR relative to the ideal profile if the incorrect process was assumed and dating models were ap-plied. Once the dating model was established, the C org accumulation rate (CAR) was estimated through Eq. (5) assuming average sediment C org contents of 2.5 % in seagrass and 8 % in mangrove-tidal marsh in both ideal and simulated sediment profiles. Under ideal conditions, CARs were 50 and 240 g C org m −2 yr −1 in seagrass and mangrove-tidal marsh sediments, respectively. While this overall model structure was used in all simulated scenarios, MAR and CAR rates under ideal conditions varied from those reported above in increasing MAR and OM decay simulations to represent real increases in accumulation, changes in OM content and associated losses of sediment mass with depth (Table 3).
where (% C orgi · m i ) is the mass per unit area of C org at layer i (g C org m −2 ), m t is the cumulative mass over the period (t) (g m −2 ) and MAR t is the mass accumulation rate of the period of interest (n − t) (g m −2 yr −1 ). When CAR is examined over the last 100 years, m t is the cumulative mass down to the 210 Pb xs horizon (i.e., the depth at which 210 Pb xs activity approaches zero) and MAR t is the mean mass accumulation rate.
3 Results and discussion

Types of 210 Pb xs specific activity profiles
Seven distinct types of 210 Pb xs specific activity profiles can be identified in vegetated coastal sediments based on examples from the literature (Fig. 2). Type I is produced by constant sediment accumulation in steady-state conditions (i.e., "ideal" profiles). The other six types of 210 Pb xs profiles summarize the most common disturbances encountered in vegetated coastal sediments that are related to the presence of mixing (physical or bioturbation), increasing MAR, erosion or alteration by intrinsic features of sediments such as heterogeneous grain size distribution and decay of OM. Type II (dark blue profile) illustrates a moderate decrease in the slope of 210 Pb xs specific activities in the upper part of the sediment core, which is often attributed to higher MAR (Cearreta et al., 2002;Haslett et al., 2003;Swales and Bentley, 2015), but can also be related to a mixing process .
Type III, showing constant 210 Pb xs specific activities along the upper part of the core overlaying an exponential decaying trend, is usually interpreted as the outcome of mixing as a result of bioturbation or sediment resuspension, redeposition, and reworking (Jankowska et al., 2016;Sanders et al., 2010a;Serrano et al., 2016a;Sharma et al., 1987;Smoak and Patchineelam, 1999). This profile type has also been related to rapid accumulation of homogeneous sediment or to recent increases in MAR (Walsh and Nittrouer, 2004).
Type IV profiles show a reverse 210 Pb xs pattern at the surface and have been attributed to a variety of factors. Similar to type III, this profile type can be caused by mixing processes in vegetated coastal ecosystems (Sanders et al., 2010b;Serrano et al., 2016a;Yeager et al., 2012) or by the deposition of allochthonous older material (Johannessen and Macdonald, 2018). It could also be produced by an acceleration of the sedimentation rate, as interpreted by Greiner et al. (2013), Smoak et al. (2013) and Bellucci et al. (2007) in seagrass, mangrove and tidal marsh, respectively, or by the decay of OM, as modeled by Chen and Twilley (1999) and Mudd et al. (2009) and observed by Church et al. (1981) in tidal marsh sediments containing > 30 % OM in top layers. Additionally, type IV profiles could also be explained by nonideal deposition (i.e., a fraction of the new 210 Pb xs input onto the sediment is not retained at the surface but penetrates to deeper layers), a process reported in peatlands and in sediments with very high porosities (> 90 %) at the sedimentwater interface (Abril and Gharbi, 2012;Olid et al., 2016).
Type V profiles show scattered 210 Pb xs specific activities, which might reflect the periodic occurrence of processes that can cause type III or IV profiles and are often interpreted as evidence of repetitive reworking in the overall mixed sediment column (Alongi et al., 2001;Serrano et al., 2016a;Smoak and Patchineelam, 1999). However, this profile form has also been explained by the deposition of 210 Pb xs outpacing its decay (λ = 0.03111 yr −1 ) (Alongi et al., 2005) or by a heterogeneous grain size sediment distribution with depth (Chanton et al., 1983;Kirchner and Ehlers, 1998;Sanders et al., 2010a), which could indicate varying 210 Pb xs fluxes due to flood events, major land use changes or changes in vegetation cover (Appleby, 2001;Marbà et al., 2015).
Types VI and VII represent low 210 Pb xs activities with depth, apparently showing low, negligible modern net accumulation of sediments. Such profiles are usually related to an abundance of coarse sediments or to erosion processes, as shown in tidal marsh sediments (Ravens et al., 2009) and bare sediments that were previously vegetated with seagrass in Greiner et al. (2013), Marbà et al. (2015) and Serrano et al. (2016a).
These examples identified from the literature reveal that various sedimentary processes might produce similar types of 210 Pb xs specific activity profiles. Any particular 210 Pb xs profile can accommodate a range of mathematical modeling approaches (see below), which leads to the development of differing chronologies and MAR estimates. Hence, the identification of the process driving accumulation and causing variation in the 210 Pb xs record aids in the determination of the C org accumulation rates.

Simulated sediment and C org accumulation rates (MAR and CAR)
We ran simulations for sedimentary processes (mixing, enhanced sedimentation, erosion) and heterogeneous sediment  composition with depth (grain size distribution and OM decay). Results of the modeled 210 Pb xs profiles are summarized in Figs. 3 and 4 and Tables S1-S7. We estimated mean 100year MAR and CAR for the simulated profiles by applying the CF-CS and CRS models, and results were compared with those from their respective ideal non-disturbed 210 Pb xs profiles. The estimated deviations in accumulation rates from those expected under ideal conditions are shown in Fig. 5 for seagrass and mangrove-tidal marsh ecosystems. These deviations are driven by variations in MAR caused by anomalies in 210 Pb profiles as the C org fraction considered to be the same in both ideal and disturbed sediment profiles.

Mixing
Simulations of surface mixing (A and B in Fig. 3a) yielded 210 Pb xs profiles similar to types II and III (Fig. 2), while deep mixing (scenario C) led to stepwise 210 Pb xs profile forms similar to type V. Calculated MAR and CAR devi-ated between 4 % and 80 % from the ideal value in seagrass sediments, while such deviation was negligible (≤ 7 %) in mangrove-tidal marsh sediments due to the smaller proportion (5 %-10 %) of the 210 Pb xs profile affected by mixing ( Fig. 5a and c). In both cases, higher deviations from the expected rates were associated with deep mixing and with the use of the CF-CS model since this model interprets any divergence from the ideal exponential decrease in the 210 Pb xs activity with depth to reflect random variation. In contrast, the CRS model is based on the 210 Pb xs inventory (I ) that is unaffected by vertical mixing. Profiles of mixing in sediments could be equally explained by an increase in the sedimentation rate in recent years. If the incorrect process was assumed and dating models were applied accordingly, mean MAR and CAR were largely overestimated in seagrass sediments, by 20 %, 30 % and 95 %, using the CF-CS model in surface (scenario A, B) and deep mixing simulations, respectively (Fig. 5b). In mangrove-tidal marsh sediments, overestimation in mean MAR and CAR was substantial (30 %) when deep mixing was considered to be caused by an increase in MAR (Fig. 5d). A process mismatch between mixing and increased sedimentation in re- cent years did not cause large deviations (between 2 % and 5 %) in MAR and CAR derived by the CRS model. The CRS model outputs are similar if mixing or changes in accumulation rates are present, albeit ages within the mixed layer cannot be reported if mixing occurs.

Increasing sedimentation rates
Simulated increases in MAR from 20 % to 200 % (scenarios D to G, Fig. 3b) resulted in similar profile forms to those simulated with surface mixing. Increases in MAR were modeled over the last 30 years, a period over which more than a 100 % increase (2-fold) was needed to produce a reversal of 210 Pb xs specific activities with depth under this simulation (type IV profiles). The influence of a change in MAR was better captured with the CRS model. The CF-CS model, in contrast, failed to account for rapid and large increases in MAR. Deviations from the ideal value ranged from 0 % to 15 % in scenarios D and E (20 % to 50 % increase in MAR) and were up to 60 % for a 100% increase in MAR (scenario F). Calculated MAR in scenario G (200% increase in MAR) was underestimated by a 30%, as piecewise dating is not applicable in profiles with constant or reversed activities of 210 Pb xs with depth. In such situations, additional tracers or time markers are required to estimate MAR and CAR in the layer of constant 210 Pb xs activities (see Sect. 4.2). Deviations from the expected value ranged from 0 % to 4 % when using the CRS model ( Fig. 5a and c). Results were similar for both ecosystem types. If the recent increase in MAR was interpreted as mixing, the mean MAR and CAR were underestimated by between 10 % and 30 % in both habitat types using the CF-CS model ( Fig. 5b and d). In contrast, the deviation from the ideal value was ≤ 5 % if the CRS model was applied.

Erosion
We ran three simulations (H, I and J) to represent recent (H) and past erosion events (I and J) (Fig. 3c). Simulations of erosion yielded lower 210 Pb xs specific activities than those of the ideal reference profile (type VII, Fig. 2), and 210 Pb xs dating horizons were found at shallower depths in these simulations (Fig. 3c). Consequently, 210 Pb xs inventories (I ) in eroded profiles were lower than expected (reference ideal profile I ref : 3900 Bq m −2 ). Inventories of simulated seagrass sediments  Figure 5. Ratio of average 100-year C org accumulation rates (CAR) between simulated and ideal 210 Pb profiles produced by various sedimentary processes in seagrass (a, b) and mangrove-tidal marsh habitats (c, b). (a, c) The correct process is assumed and models are applied. (b, d) An incorrect process is assumed and models are applied accordingly. Error bars represent the result of error propagation. Uncertainties for mean MAR were derived from SE of the regression and SE of the mean using the CF-CS and CRS models, respectively. Ratios of simulated to ideal sedimentation rates (MARs) are equal to those for CAR, determined from multiplying MAR by the fraction of C org in sediments (Eq. 5), which was considered constant between ideal and simulated profiles. In simulations of increasing sedimentation and organic matter decay, new MAR and CAR were estimated for ideal 210 Pb profiles to represent real changes in accumulation, organic matter decay and associated changes in sediment mass with depth. had a deficit of 2400 Bq m −2 (60 %), 1250 Bq m −2 (30 %) and 600 Bq m −2 (15 %) in erosion scenarios H, I and J, respectively, while these deficits were of 900 Bq m −2 (22 %), 700 Bq m −2 (19 %) and 600 Bq m −2 (15 %) in mangrovetidal marsh sediments. Because seagrass ecosystems have lower sedimentation rates, a greater proportion of the 210 Pb xs inventory was found in the top 10 cm of the sediment column and thus missing because of erosion. Simulations of past erosion events, which can be identified deeper in the profile, produced breaks in the slope of the 210 Pb xs profiles similar to those of type II (Fig. 2). Simulated erosion scenarios did not result in a large impact in the CF-CS-derived MAR and CAR estimates under the conditions of this simulation (Fig. 5). Derived MAR and CAR were underestimated by 7 % and 2 % in seagrass and mangrove habitats, respectively. The CRS model cannot be applied to eroded 210 Pb xs profiles unless the missing inventory is known and the total (I ) and depth-specific (I m ) 210 Pb xs inventories can be corrected. Assuming erosion was not a factor, the application of the CRS model to our simulated profiles underestimated MAR and CAR by up to 25% in seagrass and by 10 % in mangrovetidal marsh sediments ( Fig. 5b and d). Therefore, we caution against the use of the CRS model in profiles that show 210 Pb xs inventories that are lower than those expected via atmospheric 210 Pb xs deposition or relative to nearby undisturbed sites (Fig. 3c). The magnitude of erosion is better estimated by the deficit in inventories of 210 Pb xs , rather than by sedimentation rates. The comparison between sediment 6802 A. Arias-Ortiz et al.: Reviews and syntheses: 210 Pb-derived sediment and carbon accumulation rates records can provide information about the degree of erosion (see Sect. 4.3). In our simulations, the C org stocks over the last 100 years were 20 % and 5 % lower in seagrass and in mangrove-tidal marsh sediments, respectively, compared to the corresponding ideal non-eroded profile. Part of this is likely related to the fact that the concentration of C org is not changed, which in reality may actually change since fine sediments, in which C org is more efficiently adsorbed, are more easily eroded and OM is remineralized when exposed to oxic conditions during resuspension (Burdige, 2007;Lovelock et al., 2017a) (see simulations in Sect. 3.2.4 and 3.2.5). Consequently, losses of sediment C org could be significantly larger, as shown in some recent studies (Macreadie et al., 2013(Macreadie et al., , 2015Marbà et al., 2015;Serrano et al., 2016a).

Sediment grain size distribution
Coarse sediments are often unsuitable for 210 Pb dating as they may lead to very low 210 Pb xs specific activities. We simulated 210 Pb xs profiles in a sediment consisting of 70 % coarse sand (scenario K, Fig. 4a). This led to diluted 210 Pb xs specific activities with depth similar to those produced by erosion processes. In contrast to erosion simulations, coarse but homogeneous grain size distribution with depth did not have any impact on MAR and CAR estimated by the CF-CS model, since the dilution effect did not cause any anomaly in the slope of the 210 Pb xs profile. However, the CRS model underestimated the sedimentation rate by 10 % in both habitats (Fig. 5). The reduction of 210 Pb xs specific activities may cause the limits of detection of 210 Pb xs (0.35 Bq kg −1 in our simulations) to be reached at shallower depths than in the ideal profile. In this simulation, the limits of detection were 4 and 6 cm shallower in seagrass and in mangrove-tidal marsh sediments, respectively (Table S5a and b in the Supplement). This was conducive to the overestimation of the sediment age at bottom layers by the CRS model and underestimated mean MAR due to the omission of a higher fraction of the integrated 210 Pb xs activity per unit area from I m and I at depths greater than those at which the limit of detection was reached (MacKenzie et al., 2011). This effect is known as the "old date error" of the CRS model and can be corrected as described in Binford (1990) and Appleby (2001). Because we have assumed the same C org content in ideal as in simulated profiles, CAR estimates vary similarly to MAR. However, C org content would likely covary with grain size, and we therefore expect lower C org content with coarser sediments (Dahl et al., 2016;Sanders et al., 2012).
Simulations of varying grain size distribution with depth (scenarios L, M and N) led to stepwise 210 Pb xs profile forms (Fig. 4b). A sharp increase in 210 Pb xs specific activities in surface layers can be produced by the presence of finer sediments in which 210 Pb is preferentially associated (scenario L). As a result, sedimentation rates were 2 % to 20 % lower than those estimated for the ideal profile in both habitat types using the CF-CS and the CRS models, respectively (Fig. 5). The CRS model assumes that 210 Pb xs specific activities are inversely related to the sedimentation rate, and thus higher 210 Pb xs activities led to lower accumulation rates.
When coarser sediments dominated at the surface layers (scenario M), the simulated profiles obtained were similar to those with mixing and accelerated sediment accumulation in recent years (types II, III and IV). The dilution of 210 Pb specific activities caused by the deposition of coarse sediments in surface layers was interpreted by the CRS model as an increase in the sedimentation rate; however, this effect was compensated for in part by the old date error. With coarser sediments at surface layers, the CF-CS model applied piecewise overestimated average MAR and CAR by only 1 % in both habitat types, while the CRS model resulted in a 5 % overestimation (Fig. 5). If changes in grain size are considered throughout the entire 210 Pb xs profile (scenario N), the deviations in accumulation rates were up to 10 % using both models in both habitat types. Indeed, the deposition of coarse sediments may indicate exceptional increases in sedimentation in the case of storm surge deposits or pulsed sediment deliveries. In these cases, the CF-CS model could be applied if the event layer can be identified and can be subtracted to produce a corrected depth profile from which to determine the CF-CS-derived ages and mean mass accumulation rate. However, the presence of coarse sediments is often related to a reduction in the deposition of fine particles or to the transport and erosion of these in high-energy environments, leading to a variation in the 210 Pb xs flux onto the sediment surface that is considered constant through time by the two dating models. Where heterogeneous sediment layers are present, some corrections, such as the normalization of 210 Pb xs specific activities, are required before the application of any of the 210 Pb dating models to obtain more accurate estimates of MAR and CAR (see Sect. 4.4).

Organic matter decay
Two different scenarios with low and high sediment organic matter (OM) content (16.5 % and 65 %, respectively) were modeled in relation to OM decay. In both scenarios simulated MAR and CAR were overestimated relative to those derived from ideal profiles. Variation in OM decay (from a starting level of 16.5 %) only slightly affected 210 Pb xs specific activity profiles (Fig. 4c), causing a small overestimation of MAR and CAR of between 2 % and 5 % in both habitats and by both models under any of the rates of decay considered in this simulation (0.00005, 0.0005 and 0.01 − 0.03 d −1 ) (Fig. 5a and c). OM decay in very rich organic sediments (65 % OM) caused increased 210 Pb xs specific activities at surface (scenarios R and T) and subsurface sediments in which the decay of OM is greater, leading to a reversal of the 210 Pb xs profile in simulated scenario S. Derived CARs were 20 %-30 % higher as estimated by the CF-CS model and 10 %-20 % using the CRS model in both habitat types (Fig. 5a and c).
Mass accumulation in vegetated coastal ecosystems is the result of the balance between material accretion (detritus and sediment) from autochthonous and allochthonous sources, decomposition, and erosion (e.g., Mateo et al., 1997). Assuming there is no erosion, the estimates of MAR and CAR by means of 210 Pb are the net result of mass accumulation with time and hence integrate both the burial and decomposition of organic matter over a centennial timescale. Because mean CAR rates are based on the C org presently available and not the amount originally deposited, their determination will be dependent on the timescale over which they are calculated.

General remarks
Among the various ecosystems considered here, average last 100-year MAR and CAR derived from the CF-CS and the CRS models were less vulnerable to anomalies in mangrovetidal marsh compared to seagrass sediments. Higher sedimentation rates lead to deeper 210 Pb xs dating horizons and thus the fraction of the 210 Pb profile affected by anomalies was lower in mangrove-tidal marsh than in seagrass sediments. Anomalies caused by deep mixing or 2-to 3-fold acceleration in sedimentation had larger effects on the CF-CS-derived accumulation rates, while alterations caused by heterogeneous grain size composition or underestimation of 210 Pb xs inventories primarily affected the CRS model results (Fig. 5). Care must be taken in these cases since the discrepancy in mean MAR relative to the ideal value could range between 20 % and 80 %. Our simulations showed that the decay of OM results in an overestimation of the accumulation rate, which was most severe in very rich organic sediments regardless of the model used (> 50 % OM). However, this effect could reasonably be ignored in most cases since vegetated coastal ecosystems rarely contain OM concentrations > 25 % (Table 1), for which the deviation in computed MAR was below 10 %. Overall, simulations showed that the variability in MAR and hence CAR due to sedimentary processes and differences in sediment composition was moderately low when appropriate dating models were applied and interpreted. Deviations in the determination of MAR and CAR, generally within 20 %, confirmed that the 210 Pb dating technique is secure (Fig. 5). However, failure to account for the correct process affecting 210 Pb specific activity profiles could lead to deviations in mean MAR and CAR exceeding 20 % (Fig. 5c, d).
MAR and CAR were most overestimated, from 20 % to 95 % in simulations with low accumulation rates, when acceleration was interpreted in mixed 210 Pb xs profiles and the CF-CS model was applied piecewise. Deep mixing confounded with an increase in MAR generated the largest overestimation of mean CAR in both habitat types. In contrast, if mixing was assumed in 210 Pb xs profiles showing a recent increase in MAR, mean accumulation rates were underestimated by up to 30 % using the CF-CS model below the surface mixed layer. Indeed, the CRS model was less sensitive to anomalies in 210 Pb xs profiles; however, its application requires accurate determination of the 210 Pb xs inventory at each depth (I m ) and in the entire record (I ), which can be problematic, for instance, when all samples along a sediment core have not been analyzed or when sediment erosion has occurred at the core location. When the total 210 Pb xs inventory is underestimated, be it through erosion, poor detection limits or insufficient core length, this generates erroneous dates and underestimation of average MAR and CAR. Underestimation of accumulation rates will depend largely on the proportion of the missing fraction of the 210 Pb xs inventory from I m and I . In our simulations, MAR and CAR were underestimated by 10 % to 25 %. While uncertainties within 20 % might be acceptable for the determination of mean MAR and CAR over a centennial timescale, they may not allow for the determination of a detailed geochronology, historical reconstruction, or ascertained rates of change and fluxes at specific times. In that event additional tracers or geochemical, ecological and historical data need to be used to validate the 210 Pb-derived results and reduce uncertainties caused by anomalies in 210 Pb xs profiles in vegetated coastal sediments.

Approaches and guidelines
Retrieving reliable CAR depends on the correct determination of MAR and the diagnosis of the intervening sedimentary process. However, similarities in simulation outcomes and variations associated with anomalies in 210 Pb xs profiles point to the need for additional sources of evidence to discriminate between alternative processes and constrain 210 Pbderived estimates. 137 Cs or other independent radioactive tracers can be used to corroborate 210 Pb geochronologies. However, in its absence, geochemical information combined with knowledge of events related to land use and/or environmental changes (e.g., by means of aerial photographic evidence; Swales et al., 2015) can also be used as a tool to validate 210 Pb geochronologies and interpret 210 Pb xs profiles. In Fig. 2 we have summarized the steps to characterize 210 Pb profiles and the sedimentary processes most likely involved and suggest several techniques to complement the 210 Pb dating method to obtain reliable MAR and CAR.
Prior to analysis, researchers can have control over some factors, such as coring, sampling or sample handling, that can create artifacts in 210 Pb profiles and therefore contribute to dating error. Guidelines for core sampling for the analysis of 210 Pb and other radionuclides have been described in detail, for example in Brenner and Kenney (2013) and in the technical report IAEA- TECDOC-1360TECDOC- (2003. Some knowledge of the expected sedimentation rate is useful to decide how to section a sediment core for 210 Pb measurements, as is knowledge of the length that a core must have to reach the depth of the 210 Pb xs horizon. Low sedimentation rates 6804 A. Arias-Ortiz et al.: Reviews and syntheses: 210 Pb-derived sediment and carbon accumulation rates (∼ 1-2 mm yr −1 ) and/or coarse sediments may imply that the 210 Pb-datable part of sediment cores is limited to the very top centimeters. In such situations, fine sectioning intervals (0.5 cm) would be required. Longer cores (of about 100 cm) should be collected if high sedimentation rates are expected (several mm yr −1 ) so that the entire 210 Pb xs inventory is captured and the CRS model can be applied. These can be sliced at thicker intervals without compromising the temporal resolution of the 210 Pb record. If the order of magnitude of the sedimentation rate is not known a priori, it is best to choose fine sampling intervals (e.g, at 0.5 cm along the upper 20 cm, at 1 cm from 20 to 50 cm, and at 2 cm below 50 cm) to ensure sufficient resolution.
After collection, a visual description (e.g., color, sediment texture, presence of roots, organisms or layers) of the sediments and measurement of parameters such as water content, OM and grain size are relatively low-cost actions that provide information to interpret 210 Pb distribution and the pattern of accumulation. Indeed, the type of sediment (e.g., fine vs. coarse, rich in carbonates, homogeneous or with organic debris embedded) is a factor that should be considered (IAEA-TECDOC-1360. Coarse particles or coarsegrained carbonates for which 210 Pb xs is less preferentially adsorbed (Wan et al., 2005) may hinder the detection of any 210 Pb xs in vegetated coastal sediments. In such situations, the analysis of 210 Pb in the smaller sediment fraction (i.e., < 63 µm or < 125 µm) is recommended to concentrate 210 Pb and reduce the dilution effect caused by coarse fractions. However, the application of the CRS model would then be limited to those cases in which the 210 Pb xs is contained entirely in the sieved sediment fraction. Sieving combined with 210 Pb dating has been applied to mangrove sediments from arid regions where 210 Pb xs flux is low (Almahasheer et al., 2017) and to carbonate-rich seagrass sediments in Florida Bay (Holmes et al., 2001). Similarly, large organic material such as roots and leaves should be removed from the sediment samples prior to 210 Pb analyses as these may contribute to the dilution of the 210 Pb xs specific activity.
The analytical methods for 210 Pb measurements can also be chosen depending upon the amount of sample available and its expected specific activity. While the indirect determination of 210 Pb through alpha spectrometry of its granddaughter 210 Po requires a small amount of sample (150-300 mg) and will provide a significantly better limit of detection (< 1 Bq kg −1 ), the direct determination of 210 Pb by gamma spectrometry can simultaneously provide data for supported 210 Pb ( 226 Ra) and relevant radionuclides, such as 137 Cs, 228 Th, 7 Be and 40 K, to validate the 210 Pb geochronologies. For a detailed description of the analytical methods and their advantages and disadvantages, see Corbett and Walsh (2015) and Goldstein and Stirling (2003).

Artificial radionuclides
Independent validation of the chronology is essential to ensure a high level of confidence in the results (Smith, 2001). Varved sediments used to validate chronologies in lakes do not occur in vegetated coastal sedimentary sequences, and thus transient signals such 137 Cs or 239+240 Pu become the most commonly used option to validate 210 Pb chronologies (Lynch et al., 1989;Sanders et al., 2010a). 137 Cs and 239+240 Pu were released to the environment through the testing of high-yield thermonuclear weapons in the 1950s to 1960s and can be used as chronometers in sediments either by assuming that the peak in activity corresponds to the fallout peak in 1963 or 1965 in the Northern and Southern Hemisphere, respectively, and/or that the depth of its first detection corresponds to the onset of fallout in the mid-1950s (Ribeiro Guevara and Arribére, 2002;Stupar et al., 2014). In addition, 137 Cs can also display a peak of elevated activity in sediment cores from Europe, corresponding to the emissions caused by the Chernobyl accident in 1986, which can also help to validate 210 Pb chronologies (Callaway et al., 1996).
However, the use of 137 Cs might have some limitations in vegetated coastal sediments. Two-thirds of the 137 Cs activity released due to the tests in the atmosphere decayed after 6 decades, rendering the identification of peaks and their correspondence to the mid-1950s and 1960s depths more difficult to determine. In addition, the detection of 137 Cs is more difficult in sediment cores from habitats located in the Southern Hemisphere and near the Equator. The low 137 Cs bombtest fallout and Chernobyl inputs in these regions (Kelley et al., 1999;Ruiz-Fernández and Hillaire-Marcel, 2009), the greater solubility of 137 Cs in seawater (K d = 10 2 to 10 3 , Bruland, 1983), and the presence of sands and carbonates, particularly in seagrass sediments (Koch, 2001), are conditions that do not favor the adsorption of 137 Cs (He and Walling, 1996a) and may lead to its mobility (Davis et al., 1984). This effect could be intensified in the intertidal zone, which is not permanently submerged due to periodic changes in the water table. High contents of OM can also affect the distribution of 137 Cs in sediments as it is preferentially accumulated in leaf litter and may be absorbed by living roots (Olid et al., 2008;Staunton et al., 2002). In addition, decomposition of the organic phase in organic-rich sediments may cause the mobility of this radionuclide (Davis et al., 1984). These factors together may compromise the use of 137 Cs to validate 210 Pb geochronologies in vegetated coastal ecosystems. In contrast, Pu isotopes ( 239 Pu half-life = 24 100 years and 240 Pu half-life = 6500 years), although they are also dependent on the distribution of bomb-test fallout, would appear to offer several advantages over 137 Cs in these environments, since 239+240 Pu is relatively immobile under both freshwater and saltwater conditions (Crusius and Anderson, 1995). For instance, Sanders et al. (2016) determined sedimentation rates and 239+240 Pu penetration depths to study nutrient and CAR in intertidal mangrove mudflats of Moreton Bay, Australia. Nevertheless, because of the limitations in validating older 210 Pb dates near the base of the core and the low inventories of bomb-test fallout in sediments of Southern Hemisphere latitudes, alternative tracers might need to be used.

Geochemical information on sediments
Besides the irregular shape of 210 Pb xs profiles, the absence of a secondary radioactive tracer to validate 210 Pb results can make interpretation even more complicated. However, geochemical information in sediments can provide the potential for an additional temporal frame and can also help to explain sedimentary processes that could be misinterpreted (e.g., mixing, increasing MAR, higher primary productivity or reduction of sediment supply). Analyses of additional proxies (pollen, diatom, nutrient concentrations, stable isotopes or trace metal records; López-Merino et al., 2017) that are based on well-described historical events at the study sites (e.g., pollution, crops and land clearance) could be used in the absence of secondary radioactive tracers to corroborate 210 Pbderived dates and accumulation rates. For instance, stable Pb isotopes or total Pb concentrations in sediments are related to the history of the use of leaded gasoline in the area and can be used to identify age marks corresponding to peaks in its use or changes in lead sourcing. An example can be found in seagrass sediment cores from Florida Bay, USA (Holmes et al., 2001), and in Gehrels et al. (2005) that combines a marsh elevation reconstruction with a precise chronology derived from pollen analysis, stable isotopes ( 206 Pb, 207 Pb), 210 Pb and artificial radionuclides ( 137 Cs, 241 Am). Additionally, profiles of trace and heavy metals and of carbon δ 13 C and nitrogen δ 15 N isotopic composition of OM provide information about environmental changes for which historical information may be well known, i.e., human settlement, onset of tourism industry, temporal evolution of cropland areas or histories of variation in plant communities (Garcia-Orellana et al., 2011;Mazarrasa et al., 2017;Ruiz-Fernández and Hillaire-Marcel, 2009;Serrano et al., 2016b).

Mixing or rapid sedimentation
The methods described above for the general validation of 210 Pb models can also serve to discriminate between mixing or increasing MAR in recent years. 137 Cs and 239+240 Pu can also be used as tracers of bioturbation (Crusius et al., 2004) or acceleration of MAR during the past 50 years (Appleby, 1998;Cearreta et al., 2002;Lynch et al., 1989;Sharma et al., 1987). For instance, a demonstration of acceleration versus fast mixing could be supported when it is possible to find the distinct 137 Cs or 239+240 Pu peaks in the same zone where 210 Pb xs activities are constant (Appleby, 2001). Changes along the profiles of geochemical elements consistent with shifts in 210 Pb xs specific activities can often be associated with changes in sedimentation or erosion processes. For instance, instantaneous depositional event layers can be identified in the sedimentary record as isolated minima of 210 Pb xs specific activities (Jaeger and Nittrouer, 2006;Smoak et al., 2013), but also as variations in grain size composition, OM, water content or dry bulk density (Smoak et al., 2013;Walsh and Nittrouer, 2004) (Sect. B1). Changes in sediment mineralogy can be discerned trough X-ray radiographs, X-ray fluorescence and CAT scans (described below), but also through other radionuclides, like 226 Ra and 40 K, the profiles of which can be measured together with those of 210 Pb through gamma spectrometry. In particular, 40 K is also part of the mineral matrix and is often used as a surrogate for the lithogenic sediment fraction (Garcia-Orellana et al., 2006;Peterson, 2009;Xu et al., 2015).

Geophysical analyses
Prior to core sectioning and subsampling, nondestructive geophysical analyses such as X-ray radiographs, X-ray fluorescence (XRF), CAT scans (computerized axial tomography) or magnetic susceptibility can be conducted to identify changes in the composition of sediments with depth and in MAR or provide evidence of mixing. For instance, using X-ray radiographs many features and physical sedimentary structures may be visible (Sun et al., 2017) and, if preserved, could support the interpretation of a rapid increase in sedimentation (Walsh and Nittrouer, 2004). Pulsed sediment deliveries or erosion could be identified by discontinuous physical stratification, and sediment mixing by the presence of active burrows or the absence of sedimentary stratification (Chanton et al., 1983).

Short-lived radionuclides ( 234 Th, 228 Th, 7 Be)
Radionuclides such as 234 Th, 7 Be and 228 Th with properties such as particle reactivity and relatively short half-lives (24.1 days, 53.3 days and 1.9 years, respectively) are suitable to quantify sedimentation processes at scales from several months ( 234 Th and 7 Be) to a decade ( 228 Th). Excess 234 Th and 7 Be are sensitive indicators of mixing in the zone of constant, scattered or reversed 210 Pb xs specific activity profiles and are the most widely used (types II, III and IV, Fig. 2) (Cochran and Masqué, 2005;Sommerfield and Nittrouer, 1999). Demonstrating the presence of any of these short-lived radionuclides can give confidence that there is little material missing from the top of the sediment record and no recent erosion, which is essential for the application of the CRS model. An example is documented by Smoak and Patchineelam (1999) for a 210 Pb xs profile affected by bioturbation in a mangrove ecosystem in Brazil (Sect. B2).
Recent increases in MAR could be estimated in vegetated coastal ecosystems from the slope of the best-fit line of the plot of excess 228 Th against cumulative mass, as Alongi et al. (2005) showed in a mangrove ecosystem in Jiulongjiang 6806 A. Arias-Ortiz et al.: Reviews and syntheses: 210 Pb-derived sediment and carbon accumulation rates Estuary, China (Sect. B3). However, the use of 228 Th to derive recent increases in sedimentation is restricted to habitats with high accumulation rates (i.e., > 4 mm yr −1 ) with the last 10 years found in the upper centimeters. Its application is also limited due to the often very low excess activities in coastal waters and the constraint that sediments must contain a significant lithogenic-to-detrital fraction. In general, the use of short-lived radionuclides might be indicated to assess the magnitude of mixing or recent erosion in vegetated coastal sediments.
Mixing, either due to bioturbation or hydrodynamic energy, is the most common process affecting vegetated coastal sediment records. Although the presence of vegetation and anoxic sediments tends to reduce the depth of sediment mixing (Duarte et al., 2013), the mixed layer can extend to depths of 10-15 cm in marine sediments (Boudreau, 1994). If surface mixing occurs, valid estimates of sedimentation rates (within 5 % variability as shown in Sect. 3.2.1) can still be obtained using the dating models described above; however, this is only possible in sediments for which 210 Pb xs is buried below the mixed layer prior to decay; i.e., the residence time of sediments in the mixed layer must be shorter than the effective dating timescale (∼ 100 years) (Crusius et al., 2004). In the example from Smoak and Patchineelam (1999) (Sect. B2) in which mixing extends to a depth of 11 cm, the sedimentation rate had to be higher than 1.1 mm yr −1 in order for 210 Pb to be a useful chronometer (residence time in the mixed layer is 110 mm to 1.1 mm yr −1 = 100 years, which is within the effective dating timescale of 210 Pb).

Maximum penetration depth of 210 Pb xs
A chronology cannot be estimated if mixing affects the whole or the vast majority of the sediment record, such as in the deep mixing simulation in seagrass sediments in this work. However, information such as the total historical inventory of elements, like nutrients accumulated at a site, and the maximum conservative sedimentation rate can still be estimated. The penetration-depth method (Goodbred and Kuehl, 1998;Jaeger et al., 2009) uses the maximum penetration depth of 210 Pb xs (depth of disappearance) as a marker horizon for sediments that are ∼ 100 years old. Low surface 210 Pb xs specific activities can greatly restrict the age of the 210 Pb dating horizon; therefore, this is an issue that should be considered when establishing the age of the 210 Pb xs horizon. For surficial specific activities less than ∼ 100 Bq kg −1 this could be as few as three to four 210 Pb half-lives, i.e., 65-90 years. By locating the dating horizon independently of the alteration of sedimentary processes and of assumptions of the CF-CS or CRS models, an upper estimate of the average MAR can be derived. Note that by using this method, the rates of change or fluxes cannot be estimated and these types of 210 Pb xs profiles may be of little use in establishing chrono-stratigraphies since they are unlikely to have good records of other environmental parameters.

Erosion: 210 Pb xs inventories (I )
Assessing the extent of erosion requires the comparison of the 210 Pb xs inventories between reference, i.e., undisturbed locations (I ref ), and eroded sites (I ). Because 210 Pb xs is particle reactive, once deposited in sediments, its subsequent lateral redistribution is primarily controlled by resuspension and transport processes, and thus a deficit in 210 Pb xs inventories relative to undisturbed sediments may indicate loss or mobilization of sediment particles. This approach has been used in terrestrial soils (Martz and Jong, 1991;Walling et al., 2003) and more recently to assess erosion of seagrass sediments (Greiner et al., 2013;Marbà et al., 2015;Serrano et al., 2016a) (Sect. B4). Because the 210 Pb xs inventories at a reference undisturbed location may be spatially variable, we recommend the use of a reference inventory value based on several cores (i.e., mean ± 2 SE). The consistency of the resulting reference inventory value can then be assessed by comparing it with the 210 Pb xs inventory measured in a nearby undisturbed soil characterized by minimal slope or with that expected from the local atmospheric flux of 210 Pb xs ( ). See Preiss et al. (1996) for global and regional ranges of atmospheric fluxes. The expected inventory (I ref ) can then be derived as I ref = /λ, where λ is the decay constant of 210 Pb (0.03111 yr −1 ).

Normalization of 210 Pb xs specific activity
Dating models assume rapid and nondiscriminatory removal of radionuclides from the water column regardless of major changes in grain size or OM content along a sediment record. Radionuclide adsorption onto sediments is strongly governed by the binding capacity of the settling particles (Cremers et al., 1988;Loring, 1991), and thus its scavenging is increased by fine-grained texture (He and Walling, 1996a) and OM particles (Yeager and Santschi, 2003). Variations in the influx of these particles into vegetated coastal sediments may proportionally affect the influx of particle-bound 210 Pb xs (as long as it is still available), thus violating the assumption of constant flux of the CRS and CF-CS models and leading to subsections and irregularities in 210 Pb xs profiles. Constant or reversed patterns in 210 Pb xs activity profiles, which could be easily mistaken for reworked deposition, could be caused, for instance, by vertical fluctuations of grain size due to seasonal variations in sediment discharge or reoccurring tidal currents. Sediment studies often attempt to minimize these effects by normalizing radionuclide specific activities to granulometric or geochemical parameters that reduce the influence of preferential adsorption by fine sediments and OM (Álvarez-Iglesias et al., 2007;Loring, 1991;Wan et al., 2005), allowing researchers to obtain 210 Pb xs profiles showing an exponential decreasing trend with depth (Kirchner and Ehlers, 1998;Sun et al., 2017). Radiometric applications in coastal sediments have traditionally opted for grain size normalizers such as the < 4, < 63 µm fraction or Al content (Álvarez-Iglesias et al., 2007;Sanders et al., 2010b;Sun et al., 2017;Walsh and Nittrouer, 2004), while in dynamic, sandy-rich coastal systems where the mud fraction is small, normalization by OM content has been shown to also be effective (Van Eaton et al., 2010). Equation (6) can be used to normalize 210 Pb xs specific activities ( 210 Pb xs-NORM in Bq kg −1 ) by grain size fractions, OM content or other geochemical parameters that control the variation of the input of 210 Pb xs . 210 Pb xs-NORM = 210 Pb xs-MEAS (NP AVG /NP m ), where 210 Pb xs-MEAS is the measured specific activity of the bulk sample at depth m, and (NP AVG /NP m ) is the ratio between the core average normalizing parameter and its content at depth m. For instance, multiplication by this ratio corrects measured 210 Pb activities for variations in OM with respect to an average core value.

226 Ra specific activity profiles
210 Pb xs specific activity is determined by subtracting supported 210 Pb, assuming it is in equilibrium with 226 Ra, from total 210 Pb specific activity. This is straightforward when gamma spectrometry is employed since the total 210 Pb and 210 Pb sup (i.e., 226 Ra) can be quantified simultaneously. On some occasions, particularly when 210 Pb is determined by alpha spectrometry, 226 Ra is not measured, and supported 210 Pb is most often determined from the region of constant and low 210 Pb specific activities at depth or, alternatively, from a number of determinations of 226 Ra via gamma spectrometry or liquid scintillation counting (LSC) along the core. This method assumes that 226 Ra and 210 Pb sup are constant throughout the sediment core (Binford, 1990). However, this might not be always the case, especially in heterogeneous profiles consisting of a variety of sediment types (Aalto and Nittrouer, 2012;Armentano and Woodwell, 1975;Boyd and Sommerfield, 2016) or in records containing episodes of rapid sedimentation (Chanton et al., 1983). In addition, equilibrium of 210 Pb sup with 226 Ra might be compromised in surface sediments for which 222 Rn is deficient (Appleby, 2001). Although variations in 226 Ra specific activities with depth are small in most cases, accurate determination of 226 Ra might be crucial in sediments with low total 210 Pb (e.g., due to the presence of coarse sediments), for which slight variations in the 210 Pb sup fraction may result in significant errors in the estimation of 210 Pb xs and derived accumulation rates (Diemer et al., 2011). Therefore, we recommend measuring 226 Ra specific activity profiles or, at least, using depth-specific 226 Ra values at several depths along a sediment profile to estimate 210 Pb xs .

Conclusions
The 210 Pb dating technique provides crucial information for the study of carbon sequestration in vegetated coastal sediments and can also provide accurate geochronologies for the reconstruction of environmental processes. However, 210 Pbbased geochronologies may be difficult to conduct in mangrove, tidal marsh and seagrass ecosystems where unaltered sedimentary records are rare. Shallow vegetated coastal sediments are often affected by a number of processes such as mixing and bioturbation, accelerated sedimentation, or erosion and might be composed of heterogeneous sediments. These factors may lead to anomalies in the 210 Pb xs specific activity profiles and thus produce erroneous geochronologies and biased mean lastcentury MAR and CAR. Discrepancies in mean MAR and CAR between irregular and ideal 210 Pb profiles simulated in this study are within 20 % if the intervening sedimentary process is well diagnosed. Otherwise, these deviations may range between 20 % and 100 %, with higher errors associated with the application of (1) the CF-CS model in records showing intense mixing or large increases in MAR and (2) the CRS model in incomplete 210 Pb xs sediment profiles. Additional tracers or geochemical, ecological or historical data can be used to identify the process causing anomalies in 210 Pb xs profiles and reduce uncertainties in derived accumulation rates. Using the procedures in Sect. 4, researchers have been able to obtain reliable chronologies and CAR in vegetated coastal sediments. Special caution, however, should be applied for sites at which sediments might be altered by multiple processes (leading to profile types V or VI shown in this study) and when other chronological tools or time markers are not available (e.g., 137 Cs). Sites that have slow accumulation rates and/or intense mixing may not be datable and derived CAR estimates may be largely overestimated. Mistakes would include assigning discrete ages in mixed sediments or extrapolating an age-depth model for a core that should be considered undatable to depths down the core or to nearby sites. While attention should be paid to the limitations of 210 Pb-derived results in vegetated coastal ecosystems, the guidelines provided here should help interpret complex 210 Pb profiles obtained from vegetated coastal sediments and develop a strategy to strengthen the evaluation of MAR and CAR.
Data availability. Data to replicate the dating models and formulas used to conduct the simulations in this article can be found in the Supplement, Tables S1-S7.

6808
A. Arias-Ortiz et al.: Reviews and syntheses: 210 Pb-derived sediment and carbon accumulation rates Appendix A: Simulation methods

A1 Mixing
To simulate surface mixing (scenarios A and B), we estimated the accumulated 210 Pb xs activity per unit area over the top 5 cm of the ideal 210 Pb xs profile (I 5 cm : 2126 Bq m −2 in seagrass and 723 Bq m −2 in mangrove-tidal marsh sediments) (Table S1a and b in the Supplement). We split this inventory within the upper 5 cm using a random function, the outputs of which fell within the standard deviation (±SD) of the mean of the 210 Pb xs activities in the upper 5 cm (±107 Bq m −2 in seagrass; ±9 Bq m −2 in mangrove-tidal marsh sediments). To simulate deep mixing (scenario C), we followed the same methodology but we randomly split the 210 Pb xs inventory within the upper 15 cm, which is a depth reported as deep mixing in seagrass (Serrano et al., 2016a), mangroves and tidal marshes (Nittrouer et al., 1979;Smoak and Patchineelam, 1999) and is characteristic for marine sediments globally (Boudreau, 1994). We ran the simulation several times until we obtained three scenarios (A, B, C) of mixing encompassing a range of surface mixed layers (SMLs) ( Table S2a and b). Mixing A (k m : ∞ g 2 cm −4 yr −1 ) consisted of constant 210 Pb xs specific activities with depth in surface layers; mixing B (k m : 20-23 g 2 cm −4 yr −1 ) was characterized by a decrease in the slope of 210 Pb xs in top layers; and mixing C represented deep mixing from the sediment surface down to 15 cm (k m : 6-25 g 2 cm −4 yr −1 ). 210 Pb xs activities per unit area (A) were converted to 210 Pb xs specific activities (C) in Bq kg −1 , which we averaged every two layers to represent smooth transitions. Sedimentation and derived CAR were estimated from the modeled profiles using the CF-CS and the CRS models. The CF-CS model was applied below the depth of the visually apparent SML (3 cm) in scenarios A and B to avoid overestimation of MAR. The CF-CS model was applied to the entire profile in deep mixing scenario C in seagrass sediments and below the apparent mixed layer (13 cm) in mangrove sediments. Deep mixing affected 10 % of the entire 210 Pb xs profile of mangrove-tidal marsh sediments and 45 % of seagrass sediments. To account for the deviations in mean MAR and CAR associated with a process mismatch (i.e., as if considering the actual mixing to be caused by an increase in MAR), we applied the CF-CS model piecewise (scenarios B and C) and to the entire profile (scenario A). In the case of the CRS model, ages were determined at each layer and average centennial MAR was estimated by dividing the mass of sediment accumulated (g cm −2 ) down to 100-year depth by its age (i.e., 100 years) in all cases.

A2 Increasing sedimentation
We simulated an enhancement of MAR that could result, for instance, from increased sediment runoff due to coastal development by increasing the basal MAR (0.2 and 0.3 g cm −2 yr −1 in seagrass and mangrove-tidal marsh, respectively) by different magnitudes (20 %, 50 %, 100 % and 200 %). Increases in MAR were simulated over the top 6 cm and 23 cm of the idealized 210 Pb xs specific activity profiles, which represent the last 30 years of accumulation in seagrass and mangrove-tidal marsh sediments, respectively. Last-century mass accumulation rates expected for ideal profiles were estimated by dividing the accumulated mass down to a 100-year depth (derived from gradual increases in MAR) by its age (Table S3a and b). 210 Pb xs specific activities (C m ) as a result of increased MAR were estimated through Eq. (A1) for each layer. Simulations of increasing MAR generated four profiles per habitat type (scenarios D, E, F and G) (Fig. 3b). Average MAR and CAR were estimated from the modeled profiles using the CF-CS and CRS models. The CF-CS model was applied piecewise in scenarios D, E and F and below the layer of constant 210 Pb xs in scenario G.
where λ is the decay constant of 210 Pb (0.0311 yr −1 ) and I m is the 210 Pb xs inventory accumulated at layer m. The value 10 allows for the unit conversion to Bq kg −1 . We also estimated mean MAR and CAR assuming that the process causing scenarios D, E, F and G was mixing. For this, we applied the CF-CS model below the surface mixed layer (6 and 23 cm in seagrass and mangrove-tidal marsh sediments, respectively). The shift in the slope of the 210 Pb xs profile in scenario D in seagrass sediments was minimal, and hence we applied the CF-CS model to the entire 210 Pb xs profile, as this would likely be the method applied by most researchers in a real case. The CRS model was run similarly if mixing or changes in accumulation rates are expected, ages were determined at each layer and average centennial MAR was estimated by dividing the mass of sediment accumulated (g cm −2 ) down to 100-year depth by its age (i.e., 100 years). If mixing is expected, ages within the mixed layer cannot be reported.

A3 Erosion
Erosion in vegetated coastal sediments can occur due to highenergy events (Short and Wyllie-Echeverria, 1996), vegetation loss and subsequent destabilization of sediments  or mechanical disturbances (e.g., Serrano et al., 2016c). We ran three simulations to represent recent (H) and past erosion events (I and J) (Fig. 3c). We started with an ideal 210 Pb xs profile with a total initial 210 Pb xs inventory of 3900 Bq m −2 . To simulate erosion, we removed the 210 Pb xs inventory accumulated in the top 0-5 cm (H), middle 5-10 cm (I) and 10-15 cm sections (J) in sediments from both habitat types (mangrove-tidal marsh and seagrass). The resulting 210 Pb xs activity per unit area (Bq m −2 ) was converted to 210 Pb xs specific activities (Bq kg −1 ) by dividing by the corresponding mass depth (g cm −2 ) at each section af-ter correcting the latter for the loss of sediment layers (Table S4a and b). 210 Pb specific activities were averaged every two layers to simulate smooth transitions rather than a sharp discontinuity after an erosion event. We estimated the resulting average MAR and CAR using the CF-CS model (applied piecewise in erosion scenarios I and J). The CRS model should not be applied in simulated erosion scenarios since the overall core inventories (I ) are incomplete. However, we ran the CRS model to test the errors associated with its application in eroded sediments assuming that erosion is not a factor.

A4 Changes in sediment grain size
We simulated various 210 Pb xs profiles with changes in sediment grain size distribution using the approach described by He and Walling (1996), in which the specific surface area of particles exerts a primary control on the 210 Pb xs adsorbed: where C is 210 Pb xs specific activity (mBq g −1 ), S sp is the specific surface area of the sediment particles (m 2 g −1 ) and µ is a constant scaling factor depending upon the initial 210 Pb xs activity per unit area (mBq m −2 ). The 210 Pb xs specific activity in bulk sediments can also be represented by Eq. (A2), replacing S sp by the mean specific surface area S mean (m 2 g −1 ) of the bulk sample. In this work, we estimated µ at each layer of an ideal 210 Pb xs profile in seagrass and mangrove-tidal marsh sediments if S sp throughout the core is 0.07 m 2 g −1 , corresponding to a mean particle size of 63 µm. The surface area can be estimated as (Jury and Horton, 2004) where ρ is the density of the sediment particles and r is the mean radius of sediment particles, which are considered spherical. We estimated the weighted-mean specific surface area of a very coarse sediment composed of 70 % coarse sand (500-1000 µm), 20 % medium sand (250-500 µm) and 10 % silt (4-63 µm) (S mean = 0.0153 m 2 g −1 ) through Eq. (A3) (size scale: Wentworth, 1922). The bulk density (ρ) of sediment fractions was considered: 1.03 g cm −3 for silt, 1.6 g cm −3 for medium sand and 1.8 g cm −3 for coarse sand. First, we simulated 210 Pb xs profiles as a function of the specific surface area by applying Eq. (A2) to an ideal 210 Pb xs specific activity profile (scenario K) ( Table S5a and b). Second, we simulated a shift to sandy and clayey sediments in surface layers, which could result after the restoration or loss of vegetated coastal ecosystems. The percentages of sands and clay along the core were changed using a random function (from 60 ± 20 % in surface to 15 ± 5 % in bottom layers; scenarios L and M) ( Table S6a and b). The shift was simulated at the same age depth (30 years before collection) in all scenarios and habitat types. Finally, we simulated a heterogeneous grain size distribution along the entire sediment profile intercalating sand and clay layers randomly with depth (scenario N) (Table S6a and b). The mass depth term was corrected in each case for changes in grain size, which lead to variations in DBD with depth. The bulk density (ρ) of sediment fractions was considered: 0.4 g cm −3 for clays and 1.6 g cm −3 for medium sands. In addition, the value of µ was readjusted at each sediment depth of the ideal profile to represent non-monotonic variations in cumulative dry mass. 210 Pb xs specific activity profiles were estimated as a function of the specific surface area that was estimated at each layer according to the various proportions of clay and sand. The average MAR was estimated using the CF-CS and CRS models. The CF-CS model was applied piecewise in simulated scenarios L and M.
A5 Organic matter decay 210 Pb xs in vegetated coastal sediments is deposited in association with mineral particles but also with organic particulates (Krishnaswamy et al., 1971;Yeager and Santschi, 2003). Once buried, sediment organic matter (OM) content usually decays with sediment depth and aging due to remineralization of labile fractions, leading to an enrichment of 210 Pb xs specific activities. We simulated the resultant 210 Pb xs profiles derived from this process in two sediments with different OM contents (16.5 % and 65 %). The first value (16.5 % OM) is within the usual range of tidal marsh-mangrove and in the high range for seagrass sediments   (Table 1). The second value (65 % OM) represents an extreme scenario based on existing studies in seagrass and mangrove ecosystems (Callaway et al., 1997;Serrano et al., 2012). The simulations were run under three OM decay constants assuming the following: (1) the whole pool of OM is refractory under anoxic conditions, decaying at a rate of 0.00005 d −1 in seagrass and in mangrove-tidal marsh sediments (Lovelock et al., 2017); (2) 50 % of the refractory pool is exposed to oxic conditions, decaying at a rate of 0.0005 d −1 in mangrove-tidal marsh sediments; and (3) 50 % of the OM pool is labile and decaying fast, although exposed to anoxic conditions, at 0.01 and 0.03 d −1 in seagrass and mangrove-tidal marsh sediments, respectively (Lovelock et al., 2017).
The 210 Pb enrichment factor (η) can be determined for a given time after deposition as where χ s is the mineral fraction of sediments, χ org is the organic fraction of sediments at time 0, k org is the decay constant of the OM in sediments and t is time and can be estimated as m/MAR. As time (t) increases the exponential term tends to zero, and hence the OM stored in the sediment reaches a constant value at which it is no longer decomposed. We assume that the remineralized OM leaves the sediment 6810 A. Arias-Ortiz et al.: Reviews and syntheses: 210 Pb-derived sediment and carbon accumulation rates as CO 2 , but in fact a fraction (f ) would transform to mineral matter as χ s (t) = χ s(0) + f · χ org(0) · 1 − e −k org ·t . In our simulations f = 0 was assumed. Then, the 210 Pb xs specific activity of a sample of age t with initial specific activity C 0 is and the total mass accumulated with depth (M) above a layer of age t is M = MAR · χ S · t + MAR · χ org · e −k org ·t · t.
MAR was estimated using the CF-CS and CRS models. The CF-CS model was applied below the 210 Pb xs reversed profile in scenario S. CAR was estimated through Eq. (5). Organic matter (%OM) in mangrove-tidal marsh sediments was transformed to %C org using Eq. (A7) . In seagrass sediments we applied the relationship reported by Fourqurean et al. (2012) (Eq. A8) ( For this simulation new MAR and CAR were estimated as derived from ideal 210 Pb profiles to represent changes in organic matter content due to decay and associated losses of sediment mass with depth. This resulted in lower ideal MAR in seagrass and mangrove-tidal marsh sediments (seagrass: 0.17 and 0.07 g cm −2 yr −1 ; mangrove-tidal marsh: 0.25 and 0.10 g cm −2 yr −1 in OM decay simulations starting at 16.5 % and 65 % OM, respectively) (Tables 3 and S7).

Appendix B B1 Case study of a sedimentation event
Hurricanes and cyclones can lead to the sudden delivery of large amounts of sediments and nutrients to mangroves and tidal marshes, which in turn can result in enhanced production (Castañeda-Moya et al., 2010;Lovelock et al., 2011). Smoak et al. (2013) obtained a 210 Pb xs specific activity profile consistent with a large pulse of sediment delivered to fringing mangroves in the Everglades, Florida (Fig. B1). The specific activity of 210 Pb xs was vastly different (several times lower) in sediments accumulated during the event. The sediment accumulation rate estimated by the CRS model for the upper part of the sediment record was 6 times that of background levels, resulting in a doubled accretion rate due to the high bulk density of the delivered sediments (Castañeda-Moya et al., 2010). C org concentrations in the abruptly accumulated sediments were lower (5 %) than those of the sediments beneath the event layer (20 %-25 %). In fact, event deposits could consist of coarse sediments (for instance, sand and shell sediment layers deposited during storm events characteristic of offshore environments; Swindles et al., 2018), but also of fine sediments that could present lower 210 Pb xs specific activity compared to surrounding layers (e.g., siltation events due to clearing of the catchment area; Cambridge et al., 2002;Serrano et al., 2016d). Indeed, if the initial 210 Pb xs specific activity (C 0 ) is known, the CIC model could be useful to constrain dating when it is difficult to precisely define the thickness of such deposits. Otherwise the CF-CS model could be applied if the event layer is identified (e.g., using XRF, 226 Ra or granulometry) and can be subtracted to produce a corrected depth profile from which to determine derived CF-CS ages and mean mass accumulation rates.

B2 Case study of mixing
An example of bioturbation processes is documented by Smoak and Patchineelam (1999) in which they showed a mixed 210 Pb xs profile down to 11 cm of depth in a mangrove ecosystem in Brazil evidenced from the 210 Pb, 234 Th and 7 Be specific activity profiles (Fig. B2). The 210 Pb xs activities decreased exponentially below the surface mixed layer, resulting in an estimated accumulation rate of 1.8 mm yr −1 . In the upper layers the 210 Pb xs followed a complex pattern, with alternate relative maxima and minima, which could be representative of varying conditions of fluxes and sediment accumulation rates, the presence of coarse sediments, or physical or biological mixing. However, 7 Be penetrated down to 4 cm of depth and excess 234 Th was detected only in the surface layer. Sediments that are buried for a period of more than 6 months will have undetectable 7 Be, and hence its presence at 4 cm of depth indicated that the activity of benthic communities had remobilized it downwards to a much greater degree than sedimentation.  Figure B2. 210 Pb xs specific activity profile affected by bioturbation. Short-lived 7 Be and excess 234 Th specific activity profiles are indicators of mixing in the zone of constant 210 Pb xs (0-5 cm) (Adapted from Smoak and Patchineelam, 1999).

6812
A. Arias-Ortiz et al.: Reviews and syntheses: 210 Pb-derived sediment and carbon accumulation rates B3 Case study of rapid sedimentation rates Alongi et al. (2005) studied the rates of sediment accumulation at three mangrove forests spanning the intertidal zone along the south coastline of the heavily urbanized Jiulongjiang Estuary (China). Mass accumulation rates (MARs) were rapid and one of the 210 Pb xs specific activity profiles showed scattered concentrations with depth. This could be related to either a very high MAR during the last decades or an intense mixing down-core. However, the excess 228 Th specific activity profile, determined from the difference between the total 228 Th and 228 Ra activities in the sediment, showed a clearly decaying trend down to 15 cm (Fig. B3). The exponential decay curve fitted to the excess 228 Th profile yielded an accumulation rate of 10 cm yr −1 , which was consistent with the 210 Pb specific activity profile. Therefore, the evidence provided by excess 228 Th indicated that a very high MAR was the most plausible process responsible for the sediment record.  Figure B3. Vertical specific activity profiles of excess 210 Pb and 228 Th in core 3564 from Alongi et al. (2005), produced by a rapid mass accumulation rate.

B4 Case study of erosion
Incomplete inventories of 210 Pb xs indicative of erosion can be illustrated by the measured 210 Pb specific activity profiles in sediments from Oyster Harbor (Albany, Western Australia), some of which have been devoid of seagrass vegetation since the 1980s due to eutrophication . The measured 210 Pb xs specific activities in unvegetated sediments were relatively low, and the horizon of 210 Pb xs was detected at a shallower sediment depth than in neighboring sediments, in which seagrass meadows persisted (Fig. B4). The inventory of 210 Pb xs in the unvegetated sediment exhibited a deficit of 722 Bq m −2 compared to that in the vegetated site. This deficit results not only from the lack of accumulation of 210 Pb xs while sediments were unvegetated (30 years; atmospheric flux of 25 Bq m −2 yr −1 ), but also from the subsequent sediment erosion. These results, combined with C org analyses, showed that unvegetated sediments had an average deficit in accumulated C org stocks of 2.3 kg C org m −2 compared to vegetated sediments over the last ca. 100 years. This deficit has been produced due to seagrass loss in the 1980s but is equivalent to a loss of approximately 90 years of C org accumulation.  Figure B4. Comparison of 210 Pb specific activity profiles and inventories of 210 Pb xs and organic carbon (C org ) between a vegetated and unvegetated site. The grey area indicates supported 210 Pb specific activity (Adapted from Marbà et al., 2015).