Plant Community Establishment in a Coastal Marsh Restored Using Sediment Additions

A goal of wetland restoration is the establishment of resilient plant communities that persist under a variety of environmental conditions. We investigated the role of intraspecific and interspecific variation on plant community establishment in a brackish marsh that had been restored by sediment addition. Plant growth, sediment accretion, and surface elevation change in planted, not-planted, and nearby reference sites (treatments) were compared. Four perennial macrophytes were planted: Bolboschoenus robustus, Distichlis spicata, Phragmites australis, and Schoenoplectus californicus. There was 100% survival of the planted species, and all exhibited rapid vegetative spread. Intraspecific variation in stem height and cover was identified, and interspecific comparisons also indicated differences in species cover. Treatment comparisons revealed that final total cover at not-planted sites was equivalent to that at reference sites, and was highest at planted sites where P. australis became dominant. Species richness was initially highest at the reference sites, but final richness was equivalent among treatments. Soil surface elevation was greater at planted compared to not-planted and reference sites. Because of the rapid cover and increased surface elevation generated by planted species, the resiliency of restored coastal marshes may be enhanced by plantings in areas where natural colonization is slow and subsidence is high.


Introduction
Effective wetland restoration requires that hydrologic regimes and soils support the growth and persistence of wetlandadapted plant communities. While hydrology and soils are recognized as major influences on wetland structure and function, other factors that affect the outcome of wetland restoration efforts include disturbance regimes, topography, nutrient availability, and seed and propagule supply (Zedler 2000). Deposition of sediments dredged from open water or riverine habitats is a primary method used to create or restore wetlands. Dredging sediments within wetlands became common in the USA in the 1930's, when both bucket and low-pressure hydraulic dredges were used to excavate canals for oil and gas exploration in coastal Louisiana (Cahoon and Cowan 1988). The sediments were stacked along the canal sides, contributing to high wetland loss through altered hydrologic patterns and increased saltwater intrusion (Bass and Turner 1997). The first projects using dredged sediments to create or restore marshes were based in North Carolina in the 1960s (Craft 2016, p 289). Marsh restoration using sediment addition consists of distributing sediments within degraded marshes (i.e., marshes characterized by high plant mortality); these projects have the goal of increasing soil surface elevation to a level that supports vegetation growth (Turner and Streever 2002).
The results of adding sediments, or sediment subsidy, in marsh restoration have been studied extensively (Elsey-Quirk et al. 2019). Sediment slurries, defined as sediments with a high fluid content of about 85% liquids and 15% solids Electronic supplementary material The online version of this article (https://doi.org/10.1007/s13157-019-01217-z) contains supplementary material, which is available to authorized users. (Mendelssohn and Kuhn 2003), are often used. Wetlands that have been restored using sediment addition include degraded salt marshes (Slocum et al. 2005;Croft et al. 2006;VanZomeren et al. 2018), marshes affected by large-scale vegetation dieback following drought (Schrift et al. 2008;Stagg and Mendelssohn 2010), tidal wetlands impacted by urban development (Haltiner et al. 1997), and deteriorating brackish marshes (La Peyre et al. 2009). The influence of surface elevation on soil physical and chemical properties (e.g., bulk density, reductionoxidation potential, nutrient concentration) in sediment-restored marshes has been noted (Croft et al. 2006;Stagg and Mendelssohn 2011;Kongchum et al. 2017;VanZomeren et al. 2018). Soil characteristics in turn influence the establishment and growth of restored marsh plant communities.
Transplanting rooted seedlings and seeding have been used to accelerate plant establishment in sediment-restored marshes. Early studies focused on use of Spartina alterniflora and Spartina patens in salt and brackish habitats, respectively (Seneca et al. 1976; Barko et al. 1977). Detailed information on seeding techniques, plant propagation, and plant spacing intervals have been developed for a variety of species (e.g., Thunhorst 1993). Planting at 1-m intervals was estimated to require 50 h of labor per ha (Broome et al. 1988), at a cost of up to $8500 per ha (Utomo et al. 2010). Natural plant colonization and succession in restored or created wetlands have also been studied. For example, S. alterniflora colonized and became dominant at a created salt marsh in Louisiana, where it was also dominant in natural marshes . However, other created marshes at the site, which had comparatively higher elevation, were colonized by the high marsh species Spartina patens and Distichlis spicata and shrubs . Following sediment addition in a degraded Louisiana salt marsh, Schrift et al. (2008) found that plant cover at low elevation sites was equivalent to that of healthy reference marshes at the end of the second growing season, while cover was lower at higher-elevation restored sites. In general, research has indicated that, if soil surface elevations are appropriate, natural colonization by wetland plants will occur following sediment addition, but the time frame and pace will vary depending on the availability of seeds and propagules and site characteristics.
Alterations in environmental factors in coastal areas resulting from climate change, including increased salinity and depth of tidal flooding associated with rising sea levels, sediment erosion, and altered precipitation and temperature regimes, can act as stressors affecting the establishment and persistence of plant communities. Variation in stress tolerance can exist within and between plant species, affecting overall community structure and function. Genetic diversity in clonal wetland species has been documented (e.g., Silander 1984;Brix 1999;Seliskar and Gallagher 2000;Ivey and Richards 2001;Travis et al. 2002;Richards et al. 2004;Sweetman et al. 2013;Hughes and Lotterhos 2014;Kettenring 2016). As a reflection of this diversity, intraspecific variation in response to environmental stressors exists (e.g., Lessmann et al. 1997;Hester et al. 1998;Smith and Proffitt 1999;Howard and Rafferty 2006;Howard 2010). Genetic variation among keystone marsh species has been noted to affect ecosystem properties (Seliskar et al. 2002). Zerebecki et al. (2017) found that intraspecific variation in the salt marsh foundation species Spartina alterniflora affected plant-plant and plant-herbivore interactions, and Tumas et al. (2018) noted the importance of genotypic diversity in preserving the adaptive potential of species following wetland restoration. Additional consequences of high genetic diversity in wetlands include the promotion of species coexistence (Vellend 2006), resiliency to disturbance (Brewer and Bertness 1996;Hughes and Stachowicz 2011), increased productivity and nutrient retention (Zedler 2005;Reynolds et al. 2012), increased seed production (Sweetman et al. 2013), and altered species interactions (Proffitt et al. 2005). Field studies that examined the performance of different genotypes within species following restoration used transplanted S. alterniflora Hughes 2014) and Zostera marina, (Hughes and Stachowicz 2011), and seeds of Z. marina (Reynolds et al. 2012). To our knowledge, however, there are no published studies that examined the effect of multiple individuals within a suite of species on plant community composition and structure at wetland restoration sites.
We investigated the influence of intra-and interspecific variation on plant community establishment and structure in a brackish marsh that was restored using sediment addition. The overall goal was to provide a base for the development of a vigorous plant community capable of sustaining marsh elevation, thereby preventing further marsh degradation and conversion to open water habitats. Plant community establishment, species composition, soil accretion and soil surface elevation change were documented. Marsh treatments included sediment-addition sites with vegetation plantings, sedimentaddition sites without plantings, and nearby reference sites that were not affected by sediment addition. Our study addressed the following questions: (1) does intraspecific variation within planted species result in growth differences among individuals of the species? (2) does interspecific variation affect plant community composition and structure on planted sediment addition sites? (3) do species cover and composition vary among marsh treatments? and (4) do soil properties, accretion rates, and surface elevation change vary among treatments?

Study Area
The study was conducted in a brackish marsh along Bayou Dupont in the Barataria Basin of southeast Louisiana, USA (Fig. 1). The nearest tide station, located about 35 km south at Grand Isle, Louisiana, experiences a diurnal micro-tidal range of 0.32 m and a relative sea-level rise rate of 9.08 ± 0.42 mm yr. −1 (NOAA 2019). Prior to restoration the marsh was in a degraded state, as evidenced by numerous shallow (< 0.75 m deep) interior open-water areas. The dominant emergent plant species was Spartina patens; other commonlyoccurring species included Distichlis spicata, Lythrum lineare, Ipomoea sagittata, Kosteletzkya pentacarpos, Solidago sempervirens, and Seutera angustifolia. The restoration goal was to raise soil surface elevation to a level capable of supporting the growth of emergent vegetation. A small, mobile, low-pressure hydraulic dredge (pipe diameter of 30.5 cm) was used in the summer of 2000 to pump sediment from the adjacent bayou into open water areas of the marsh; sediment was added to an elevation about 15 cm above the remaining marsh surface to compensate for dewatering and compaction. An area of about 65 ha received sediment during the dredge operation.

Plant Material
The following commonly-occurring macrophytes native to coastal Louisiana marshes were included in the study: Bolboschoenus robustus (sturdy bulrush), D. spicata (saltgrass), Phragmites australis (common reed), and Schoenoplectus californicus (California bulrush). All are perennial, clonal species that regenerate through both vegetative propagation and sexually through seed production. They occur across a range of salinities in Louisiana marshes; B. robustus and D. spicata are found primarily in brackish (5-20 psu) to salt (20-30 psu) marshes while P. australis and S. californicus occur primarily in fresh to intermediate marshes (0.5-5 psu); however, P. australis has been found in brackish marshes up to 20 psu (Penfound and Hathaway 1938;Chabreck 1972). To obtain a variety of genetic individuals, five to ten individual plants of each species were collected across the approximately 340 km-long Louisiana coastline along the Gulf of Mexico. Several culms connected by a single rhizome were retained to represent a single individual plant, and just one individual per species was collected at a location. Genetic analyses were not conducted; we assumed that genetic variation was present due to the distance between collection locations. The plants were transported to a greenhouse facility in Lafayette, Louisiana, where they were separated into ramets and vegetatively propagated through several generations under common garden conditions (freshwater, flooded to sediment surface). The original plants collected were assigned identification numbers; the propagated ramets  that identification number (individual number) and were used in a series of experiments. Additional details on field collection procedures are provided in Howard and Rafferty (2006). For this study, experimental plants were conditioned to 13 psu salinity before transport to the field site.

Experimental Design
Three locations along Bayou Dupont were selected in May 2001, approximately 9 months after sediment addition. Four sites were established at each location; two sites were randomly assigned as planted treatments (hereafter P) and two as not-planted treatments (hereafter NP). All sites on the sediment deposits were void of vegetation at the experiment's initiation. A healthy-appearing marsh site that was not affected by sediment addition and supported high vegetation cover was selected near each location as a reference treatment (hereafter R). The study therefore included six replicates of the P and NP and three replicates of R, for a total of 15 sites.
Each site (P, NP, and R treatments) was protected by an exclosure to prevent herbivory by the introduced rodent Myocastor coypus (nutria), a species that causes extensive vegetation damage in Louisiana wetlands (Taylor and Grace 1995;Evers et al. 1998;Ford and Grace 1998). The exclosures were constructed of vinyl-coated wire about 1.5 m in height that was sunk 30 cm into the substrate to prevent M. coypus from burrowing beneath. The dimensions of the exclosures (hereafter called units) were 9.5 m × 7.75 m. Because of the unconsolidated substrate, boardwalks were placed within each unit to prevent damage during data collection. The boardwalks divided the unit length into four columns, and seven rows in each column were designated; this provided a grid system with 28 cells (see Howard et al. 2008 for details). In the planted treatment each species was randomly assigned to one of the four columns, and within each column different individual plants within the species were randomly assigned to rows 1-6; row 7 was reserved for destructive sampling of soil accretion and elevation (see below). A 1-m 2 area at the center of each cell was identified, and two ramets of an individual plant of the assigned species were planted 30-cm apart. Six different individuals of B. robustus, P. australis, and S. californicus were planted; the planting of an individual of Paspalum sp. mistakenly identified at collection as D. spicata resulted in just five individuals of this species in the experiment.
Rod surface elevation tables (RSETs) were installed at all treatments to measure relative soil surface elevation, or vertical land motion (VLMw) of the wetland surface over time (Cahoon et al. 2002;Cahoon 2015). The RSET device consists of a series of connected steel rods driven into the sediment to refusal, serving as a benchmark; see Cahoon (2015) for details on measurement techniques. Marker horizons (MH), consisting of a clay (feldspar) layer spread on the soil surface, were placed near each RSET to allow measurement of sediment accretion (Cahoon and Turner 1989). The RSET-MH stations were placed in row 7 cells in one of the two planted units at each location. Each designated cell 7 was planted as described above with two ramets of a randomlyselected individual of the species assigned to that column. This resulted in 12 RSET-MH installations in the planted treatment (four installation at each of three locations). A single RSET-MH station was installed in a randomly-determined row-7 cell of one of the not-planted treatments at each location, as well as at each reference treatment. These six additional installations resulted in a total of 18 RSET-MH stations.

Data Collection
Plant Growth Data were collected within a 1-m 2 quadrat positioned at the cell centers at approximately monthly intervals through two growing seasons (May 2001to November 2002 and in April 2003, for a total of 19 measurement times. All cells were sampled in P (n = 28); in NP and R, two cells in each column were randomly selected as permanent sampling locations (n = 8). Cover class by species, total cover, planted species maximum height, and planted species distance of lateral spread were measured. The cover class of bare substrate was also recorded. Cover class was estimated using a modified Braun-Blanquet system as follows: (P) <1%, (1) 1-10%, (2) 11-25%, (3) 26-50%, (4) 51-75%, and (5) 76-100%. Because vegetation often occurred in layers, it was possible for total cover in a quadrat to exceed 100%. The number of stems and stem height of all species were measured within a 15 cm × 15 cm (0.0225 m 2 ) sub-quadrat located at the center of the 1-m 2 quadrat in September 2001and in April 2002; only live stems with at least 50% green tissue were included. At the conclusion of the experiment a soil core (10.2 cm diameter × 35 cm deep) was collected from one randomly-selected cell in each column of the 15 units. The cores were placed on ice until transported to the laboratory, washed to remove soil, and dried at 80°C to a constant mass for determination of total belowground biomass.
Salinity and Soil Characteristics Surface and soil interstitial water salinity were measured during each site visit. Interstitial water was collected at 15 cm depth at two to three randomly selected locations within each unit using a syringe attached to thin, rigid plastic tubing inserted into the soil. A soil core (10.2 cm diameter × 30 cm deep) was collected outside of and adjacent to each of the 15 units in July 2001 and March 2003. The cores were divided into depth categories (0-10 cm, 10-20 cm, and 20-30 cm) and stored on ice for transportation to the laboratory. Wet weight of each core section was determined, and samples were placed in a forced-air oven at 80°C and dried to a constant mass for determination of bulk density. The dried soil was ground in a commercial mill, and subsamples were used for determination of percent organic matter using the loss on ignition method (modified from Ben-Dor and Banin 1989) and percent sand, silt, and clay (Sheldrick and Wang 1993).

Soil Surface Elevation and Sedimentation
The sites were allowed to recover from installation disturbance for four weeks prior to baseline RSET measurements, which were obtained in June 2001. Four compass positions were selected for VLMw measurement, and at each position nine elevation measurements were obtained (see Cahoon 2015 for details). At the time of initial VLMw measurements two marker horizons were placed in each cell in the vicinity of the RSET benchmark. Soil accretion was first measured during the subsequent VLMw measurement, in October 2001. A hollow copper rod was inserted into the soil surface and liquid nitrogen gas was pumped into the rod, creating a layer of frozen soil attached to the rod. The depth of soil above the marker horizon was measured to the nearest millimeter. Elevation and accretion measurements were obtained using these techniques at approximately three month intervals over two years, with final measurements in April 2003; a total of seven postbaseline data collections were made.

Statistical Analyses
Plant Growth All cover class variables were converted to midpoint values for analyses. Intraspecific variation was examined for the first growing season only; because the plants grew rapidly, individual plants could not be visually distinguished after that time (i.e., individuals of a species grew into adjoining cells). To investigate intraspecific variation, species were analyzed separately using repeated measures analysis of variance (ANOVA) with individual plant identification as the independent variable. Cover of individual plants and their maximum height were dependent variables using the first four measurement times (May to September). Maximum lateral spread as a response was examined using data collected two times in the first growing season, in August and September, also using repeated measures ANOVA. Data quantifying individual plant stem number and stem mean height were examined for September of the first growing season using ANOVA.
Individuals within a species were pooled for analyses of interspecific variation. To ensure sufficient degrees of freedom to conduct analyses of cover variables, a subset of six of the 19 total measurement times was used. The subset consisted of early, mid, and late growing-season dates throughout the study (July and September 2001;April, July, and September 2002;and April 2003). The first date was about eight weeks following planting. Dependent variables in the repeated measures ANOVA were planted species percent cover and relative percent cover (i.e., the percent of total cover in a quadrat that was the planted species), and planted species maximum height.
Maximum lateral spread of the planted species was analyzed using four measurement times in the first growing season, from early August to late October. All four measurement times for stem characteristics were used in repeated measures ANOVA with planted species stem number and mean height as dependent variables.
To examine treatment differences (i.e., P, NP, R), the 1-m 2 quadrats within the units were designated as subsamples; cells containing RSET-MH stations were not included in these analyses because of destructive sampling within them. Subsample number was therefore 24 in P and eight in NP and R. The six measurements times described above were included in the repeated ANOVA analyses for treatment comparisons. Response variables examined were total cover and species richness, and percent cover of P. australis, Spartina alterniflora, Spartina patens, and bare ground. Stem number and mean height were also examined over the four measurement times by using repeated measures ANOVA. Similarity of vegetation between treatments was examined by using the Sørensen coefficient (equivalent to the Steinhaus coefficient), as described in Kent (2012). Both presence/absence and cover data (mean by species per treatment) were used. Analyses were conducted using two collection dates: the peak of the second growing season after planting (August 2002) and early in the third growing season (April 2003). We also investigated similarity between treatments at these two dates using analysis of similarity (ANOSIM) and multiresponse permutation procedures (MRPP). Treatment differences in belowground biomass were analyzed by using one-way ANOVA with the four samples within a unit treated as subsamples.
Salinity and Soil Characteristics Surface and interstitial salinity treatment comparisons were made using two-way ANOVA with treatment and time as independent variables. Marsh treatment comparisons of bulk density, percent organic matter, and sand, silt and clay content were conducted separately by soil depth class and time using ANOVA. In addition, two-way ANOVA with treatment and time as independent variables was used to examine percent sand, silt, and clay responses with soil depths pooled.
Soil Surface Elevation and Sedimentation Invalid VLMw measurements (e.g., pin touched unmovable debris) were discarded. The mean of pin measurements for each compass position was determined, and elevation for the RSET station was defined as the mean of the four compass positions. Sediment accretion in the vicinity of each RSET was calculated as the mean of the two marker horizons. The four RSET-MH stations in planted treatments were included as subsamples in these analyses. Data were first examined within each treatment to determine if significant change occurred over time (i.e., if conditions varied from baseline). VLMw data were expressed as cumulative change from baseline conditions and were examined using linear regression analysis. Accretion data were analyzed in the same manner as VLMw data. Linear rates of change among treatments were compared using regression analysis.
All repeated measures ANOVA (PROC MIXED), ANOVA (PROC GLM), and regression analysis (PROC REG), were conducted using SAS software (ver. 9.3; SAS Institute, Inc., Cary, North Carolina). Response variables that did not meet required assumptions for normality and homogeneity were transformed. The significance level for analyses was α ≤ 0.05; for multiple comparisons the significance level was adjusted using the Bonferroni correction. ANOSIM and MRPP analyses were conducted using R software, vegan package (Oksanen et al. 2019).

Plant Growth
Intraspecific Variation All individuals of the four species survived the first growing season following planting. There were no differences among B. robustus individuals and few differences among individuals of the other species (Table 1). One D. spicata individual (plant 5) had higher cover (59.46 ± 5.13%; mean + SE) throughout the season compared to one other (plant 6; 46.60 ± 4.62%) regardless of measurement time; cover of this individual, however, was equivalent to that of the remaining three. The time/plant interaction found for D. spicata maximum height was expressed by greater height in plant 5 compared to two other individuals in July and three others in August. The time/plant interactions for P. australis and S. californicus maximum height were expressed as increasing values for some individuals throughout the summer, while in others the values did not differ between July and August. The maximum distance of lateral spread increased over time for all species, but there were no differences among individual plants within a species (Table 1). Analysis of stem number and mean stem height data collected in September of the first growing season also revealed no significant differences between individuals in any of the species (results not shown).
Interspecific Variation Species differences in growth responses were identified for all variables examined (Table 2). In the first growing season lateral spread of B. robustus and P. australis exceeded that of S. californicus and D. spicata (Fig. 2), while spread of S. californicus was less than that of all other species. The species with highest percent cover varied over time. For example, in July of the first growing season S. californicus cover (21.44 ± 4.72%) was equivalent to D. spicata (47.17 ± 3.27%), and lower than that of B. robustus (57.03 ± 6.91%) and P. australis (58.14 ± 5.42%). In July of the second growing season, however, P. australis cover (88.00 ± 0%) exceeded that of the other species (B. robustus, 3.69 ± 2.00%; D. spicata, 20.92 ± 4.88%; S. californicus, 21.44 ± 4.72%), and this relationship persisted throughout the remainder of the study. The decrease in abundance for species other than P. australis was also illustrated by the significant species/time interaction for relative percent cover (Fig. 3a). Maximum stem height of P. australis and S. californicus were similar early in the experiment, but height of the former species exceeded all other species by the end of the first growing season (Fig. 3b)

Treatment Comparison
We identified 33 plant species that occurred within the quadrats over the 23-month study (Table S1). Colonization at NP was slow during the first growing season, as indicated by low total cover values (Fig. 4a). By spring of the second growing season, however, there was rapid establishment of plants at NP, and total cover increased rapidly. Total cover by the conclusion of the study was significantly greater at P (110.38 ± 1.82%) compared to both R (64.92 ± 4.10%) and NP (56.40 ± 4.03%); cover at the latter two was equivalent (Fig. 4a). Species richness at NP also remained low during the first growing season and, along with richness at P, was significantly lower than richness at R early in the study (Fig. 4b). However, final richness was equivalent among the treatments (NP = 2.48 ± 0.33; P = 3.12 ± 0.39; R = 4.08 ± 0.48).
Common plant species varied with treatment. While an equal number of ramets of three of the species (B. robustus, S. californicus, P. australis; two fewer ramets of D. spicata were planted) established in P, P. australis became dominant by the beginning of the second growing season (Fig. 5). Aboveground runners of P. australis originating in P allowed this species to establish in some of the NP units; P. australis did not occur in R. Spartina alterniflora was dominant at NP, where its cover at the end of the second growing season was approximately 60%. Cover of S. alterniflora was low in R and negligible in P (Fig. 5). Spartina patens occurred only in R, where it displayed an overall pattern of decreasing cover over time. Cover of S. patens in R was significantly lower at the conclusion of the experiment (12.52 ± 9.32%) compared to July of the first growing season (40.56 ± 25.27%); the dominant cover species in R at the end of the experiment was D. spicata (40.50 ± 19.22%) (  Table 3). These differences were no longer present at the conclusion of the study (R, 13.23 ± 6.97%; P, 5.82 ± 3.67%; NP, 17.88 ± 5.13%).
Total stem number was significantly lower in NP (0.27 ± 0.21) than in P (7.45 ± 0.41) in September of the first growing season; at all other measurement times, however, stem number was equivalent among the treatments. Mean stem height was also lower in NP (0.84 ± 0.68 cm) than in P (59.95 ± 6.75 cm) in September of the first growing season. In April of the second growing season mean height was lower in NP (8.04 ± 3.72 cm) than in P (90.79 ± 12.01 cm) and R (51.55 ± 3.33 cm). In September of the second growing season there were no significant differences among treatments, but by April of the third growing season stem height in P (152.19 ± 15.77 cm) exceeded that in R (45.18 ± 0.58 cm) and NP (57.09 ± 6.09 cm).
The Sørensen index identified low similarity between all treatment comparisons at both times examined (Table 4). Similarity based on presence/absence data tended to be higher than that based on cover data. ANOMSIM and MRPP analyses indicated the treatments were significantly different (P < 0.0001, all comparisons). Three species -D. spicata, S. alterniflora, and Symphyotrichum tenuifoliumwere found in all treatments at the conclusion of the study (Table S2). Belowground biomass differed by treatment (P < 0.0001);

Salinity and Soil Characteristics
Surface water salinity ranged between 1.2 and 9.3 psu during the study, and no treatment differences were found (P = 0.3850). Variation in surface salinity occurred over time (P < 0.0001). Interstitial water was difficult to collect at the sedimentaddition sites, resulting in an insufficient number of data points for analysis. Values for interstitial water salinity ranged from 1.2 to 13.0 psu. Soil bulk density varied among treatments in the upper 10-cm soil layer at both measurement times (P = 0.0037, July 2001; P < 0.0001, March 2003); this property was higher at NP and P than at R (Table 5a). Percent organic matter also varied in the upper 10 cm soil layer at both times (P = 0.0115, July 2001; P = 0.0003, March 2003), and was higher in R compared to NP and P (Table 5b). No treatment differences in soil sand, silt, or clay content were identified.

Sedimentation and Elevation Change
The depth to resistance for steel rod placement during RSET installation varied between 14.6 and 18.3 m. Significant accretion was measured at all treatments, while only P displayed significant VLMw (i.e., an increase in elevation) ( Table 6). Positive subsurface change, indicating subsidence, was documented at all treatments (Table 6). Treatment comparisons identified no differences in accretion (P = 0.2071) or subsidence (P = 0.3919). However, a significant difference was found for VLMw (P = 0.0186); elevation change at P was greater than that at R. Elevation change at NP was intermediate compared to P and R, and did not differ from either (Table 6).

Discussion
Sediment addition to the degraded marsh at Bayou Dupont raised soil surface elevations in interior ponded areas to levels that supported the growth of emergent marsh macrophytes, and the plant communities that established persisted throughout our 23-month study. The manipulation of plant succession through plantings revealed interesting insights concerning the utility of such plantings in coastal wetland restoration. We documented intraspecific growth variation in three of the planted species -D. spicata, P. australis, and S. californicus during the first growing season following planting. These findings are consistent with results from greenhouse (Vasquez et al. 2005;Howard and Rafferty 2006;Howard 2010), and field (Seliskar and Gallagher 2000) studies. The relative magnitude of intraspecific variation found, however, was low. The difference in cover value between two D. spicata individuals had little practical significance; cover of those two overlapped with the other individuals, resulting in a lack of a plant with consistently greater growth. The intraspecific differences in P. australis and S. californicus responses also had minor biological significance; again, an individual plant that displayed consistently greater growth was not identified. Therefore, the data provide little support for a conclusion that intraspecific variation influenced plant community structure. It should be noted, however, that environmental conditions at the study area apparently imposed little stress, as evidenced by the lack of mortality and rapid spread of the planted species. Surface water salinity at the study sites was less than 10 psu, which is within the normal range for these intermediate to brackish marsh species, and there was no evidence of extended flooding (e.g., no water marks on unit support posts). Greenhouse studies that identified higher degrees of intraspecific growth variation in these species imposed salinity up to  (Howard and Rafferty 2006;Howard 2010). We recognize that strength of our conclusions regarding the role of intraspecific variation in structuring plant communities is limited by the lack of genetic analyses confirming the presence of different genets of the plants studied (but, see below for an exception concerning P. australis).
Our assumption that the individual plants collected represented separate genets within a species may not have reflected actual conditions. We believe, however, that the assumption of intraspecific variation among the plants was reasonable based on studies that have indicated high levels of genetic diversity in clonal wetland plants (e.g., Richards et al. 2004). Although intraspecific variation identified among the planted species was of little consequence during initial community establishment, competitive variation among P. australis individuals was noted during the course of the study. After this experiment was initiated, information was published describing the presence of an invasive Eurasian lineage of P. australis in North America (Saltonstall 2002). Genetic analyses of specimens from our original plant collection revealed that, of the six P. australis individuals planted, three were the Eurasian haplotype M and three were haplotype I (Howard et al. 2008). The later lineage is also found in South America and Asia, and genetic analyses by Lambertini et al. (2012) indicated it is related to the African species Phragmites mauritianus. While we planted the two lineages at about the same density, genetic analyses of plant tissues collected in the middle of the second growing season revealed that haplotype M plants had become dominant and represented between 67% and 100% of the plants sampled (Howard et al. 2008). The  Fig. 3 Planted species (a) relative percent of cover (i.e., percent of total cover that is the planted species) and (b) maximum stem height (cm). Bars indicate the standard error of the mean values superior competitive performance of the M haplotype identified conflicts with findings from a common garden study conducted in freshwater conditions with individually potted plants (Mozdzer et al. 2016). In that study, haplotype I had greater total, aboveground, and belowground biomass and maximum stem height compared to haplotype M, but lower stem density (Mozdzer et al. 2016). The difference between that study and ours may be related to the exclusion of competitive interactions in the greenhouse study. Several studies have documented greater growth in haplotype M compared to native lineages of P. australis (e.g., Vasquez et al. 2005;League et al. 2006;Price et al. 2014).
Regarding our second question related to the effects of interspecific variation, we found a major influence on plant community composition and structure at P. Two species, B. robustus and P. australis, spread most rapidly during the initial few months of the study. A contraction in the spread of B. robustus and D. spicata at the end of the first growing season (Fig. 2) may have indicated the onset of winter senescence. In September, five months after planting, it became apparent that P. australis was displacing the other planted species. Although not verified experimentally, the mechanisms for the superior competitive ability of P. australis were likely its greater mean height (1.65 m), leading to light limitation on growth of the shorter species, and prolific belowground production, as documented in other studies (League et al. 2006;Moore et al. 2012;Mozdzer and Megonigal 2012). The competitive interactions between plant species at the restoration site may have been affected by salinity and nutrient dynamics. The surface water salinity measured, which was consistently lower than 10 psu, likely provided conditions more favorable to the growth of Not Planted Planted Reference Fig. 4 Mean (a) total cover (%) and (b) species richness over time by treatment (not planted, planted, reference). Bars indicate the standard error of the means P. australis, a species adapted to low to moderate salinity, than to B. robustus and D. spicata, which are tolerant of moderate to high salinity. The restoration site is located in the Mississippi River Deltaic Plain, where nutrient loading through agricultural activities and urbanization have altered water quality over the past two centuries (Turner and Rabalais 2003). We did not measure soil and water nutrient levels, but these factors can have a strong influence on interactions between plant species (Zedler and Kercher 2004). The invasive P. australis lineage was noted to have greater floret and inflorescence production compared to native lineages under elevated nutrient conditions (Kettenring et al. 2011), and eutrophication played a central role in the expansion of invasive P. australis along the northeastern Atlantic coast of the United States (Holdredge et al. 2010). Nutrient enrichment, therefore, may have been an additional factor promoting P. australis growth at our sites. Aerial imagery of the restoration site showed that P. australis was present in large patches eight years after planting, and that percent cover increased up to 22 times over six years (Howard and Turluck 2013). The fact that we intentionally planted P. australis at the restoration site may be viewed as controversial by restoration ecologists based on the invasive nature of the non-native lineages. The ecosystem impacts  Bold font indicates the significant main factors or highest-order interactions examined for effect on the response variable of P. australis are often described as negative, and include reduced species richness and altered nutrient cycling (reviewed in Meyerson et al. 2000). Native P. australis is considered to be threatened along the Atlantic coast of the United States because it is being displaced by the invasive lineage (Meadows and Saltonstall 2007). Other studies, however, have identified positive ecosystem effects, including high rates of carbon sequestration and soil stabilization (reviewed in Kiviat 2013). In Louisiana P. australis is considered to be a native species. An early vegetation survey of southeastern Louisiana wetlands noted that Phragmites communis (syn. P. australis) was a typical species in the Bcane zone^, often as a co-dominant in association with Spartina cynosuroides (Penfound and Hathaway 1938). By the late 1960's, P. communis was described as a dominant species in both fresh and brackish marshes of the Mississippi River Delta (MRD) (Chabreck and Palmisano 1973), and analyses of survey data from 1968 documented the occurrence of P. australis across the Barataria Basin ).
As noted previously, when the experiment was initiated in 2001 information on the existence of invasive, non-native lineages of P. australis had not yet been published. Because our goal was to encourage the development of a plant community that would be resilient under stressful conditions along the Louisiana coast (i.e., high subsidence and rapid sea-level rise), P. australis was considered a reasonable candidate for inclusion in the study based on its characteristics of readily trapping sediments and contributing to elevation gain (Rooth and Stevenson 2000). In addition to having high competitive ability, P. australis has been described as very resilient and capable of maintaining vigor in stressed ecosystems (Hershner and Havens 2008), attributes that may be considered desirable in coastal areas threatened by rising sea levels. Lambertini et al. (2012) confirmed that invasive P. australis lineages are now dominant in the MRD. A recent dieback of P. australis in the MRD has generated widespread concern over wetland loss because dead stands have often converted to open water (Knight et al. 2018). This highlights a situation where, although invasive, P. australis is providing critical ecosystem services in the form of storm protection and wildlife habitat. Marsh treatment comparisons revealed significant differences in plant communities composition and structure. Colonization at NP was slow during the first growing season, with mean cover of just 5.80 ± 2.52% in early October, 14 months following sediment addition. This low establishment rate was likely the result of a depauperate seed bank, which has been noted for other restored wetlands (Morzaria-Luna and Zedler 2007). Also, it is unlikely that rhizome fragments of clonal marsh species were present in the sediments added to the degraded marsh because those sediments were dredged from the bottom of the adjacent canal. Rapid colonization (within six months) in a New Jersey salt marsh restored with sediment addition was attributed to rhizome recruitment (VanZomeren et al. 2018). At our site, S. alterniflora was the primary colonizing species and initially established during the first growing season through seed germination. Active seed dispersal mechanisms have been noted for S. alterniflora in Louisiana (Utomo et al. 2009), with wind as the primary means for local dispersal (Elsey-Quirk et al. 2009). Salt marsh restoration using sediment addition has been noted to increase sexual reproductive capacity in S. alterniflora at the patch scale (Jones et al. 2019). Other plants that established at NP during the first growing season were Elodea canadensis, Eleocharis sp., and Paspalum sp. The only planted species found at R was D. spicata. We did not conduct a vegetation survey in the vicinity of the restored sites, so it is not known if the other planted species were initially present as possible nearby seed sources. However, three of the planted species -B. robustus, D. spicata, and P. australisestablished at NP during the second growing season, likely through both seed germination and vegetative propagation.  Rapid growth of S. alterniflora and other colonizing species through the second growing season resulted in equivalent total cover at NP and the other treatments at the conclusion of experiment. A similar 2-year time frame for cover equivalency was found at Louisiana salt marsh restoration areas (Mendelssohn and Kuhn 2003;Schrift et al. 2008). Sediment additions that raised elevation 2-11 cm above MSL were identified as supporting marsh stability, with additions above the 11 cm threshold leading to decreased stability (Stagg and Mendelssohn 2011). Elevation at the Bayou Dupont site was initially raised 15 cm above the surrounding intact marsh; plant establishment and growth documented during our study indicates that surface elevation was decreased to a suitable range, probably through soil compaction over time.
Despite similar final total cover values between marsh treatments, Sørensen similarity and ANOSIM analyses revealed that species composition varied among the treatments. Dominant plant species at the experiment's conclusion were S. alterniflora, P. australis, and D. spicata in NP, P, and R, respectively. The declining trend in S. alterniflora cover at the end of the second growing season at NP (Fig. 5) was likely due to its replacement by species with greater competitive ability at the moderate salinity levels that characterized the study site. The initial dominant species at R was S. patens, and the reason for the shift to D. spicata is unknown. Both are typical high marsh species and have similar patterns of aboveground and belowground biomass allocation; however, D. spicata has superior colonizing ability compared to S. patens (Brewer et al. 1998a). Distichlis spicata responded positively to wrack disturbance while S. patens responded negatively (Brewer et al. 1998b). The species shift we documented may therefore be related to disturbance (e.g., hurricane effects). Belowground biomass was lower at not-planted compared to planted and reference sites. This is consistent with a study that found belowground biomass at a restored salt marsh remained lower than that at reference sites seven years following sediment addition (Tong et al. 2013). Similarly, seven years were required before belowground biomass in restored brackish-marsh ponds approached that of reference marshes (La Peyre et al. 2009).
Properties of the sediment used for wetland restoration have an important influence on restoration outcome (Mendelssohn and Kuhn 2003;Kongchum et al. 2017;VanZomeren et al. 2018). We found differences in soil bulk density and organic matter between marsh treatments, but those differences occurred only in the top 10-cm layer. The reference marsh had lower bulk density and higher organic matter content compared to sediment-addition marshes. This was expected because plants contribute organic material to the soil over time, reducing soil bulk density. Surface bulk density in R was higher than the 0.08-0.11 g cm −3 range identified in a degrading Louisiana brackish marsh (Nyman et al. 1994). In salt marshes, bulk density below 0.20 has been identified as a threshold below which growth of S. alterniflora is not supported (DeLaune et al. 1990). The increase in bulk density at our site resulting from sediment addition likely improved conditions for plant growth.
Vertical accretion and subsurface change were significant at all treatments. We did not document treatment variation in accretion or subsurface change, but soil surface elevation differences were found. Because surface elevation did not change at R and NP, accretion compensated for the significant subsidence that occurred. In contrast, elevation increased at planted sites. As previously noted, P. australis, a species that can generate high root and rhizome biomass, became dominant at the planted sites. It can be concluded that growth of P. australis belowground structures drove root zone expansion at planted sites. The contribution of P. australis to soil processes was also noted in Chesapeake Bay in Maryland, USA, where the large amount of surface litter in a P. australis stand acted to trap organic and mineral matter; the vertical accretion rate was found to be 3-4 mm yr. −1 greater in the P. australis stand than in adjacent stands of Typha spp. and Panicum virgatum (Rooth et al. 2003). If P. australis continues to spread aggressively at the Bayou Dupont restoration site, it is expected that this species will have a significant impact on soil surface elevation dynamics over time.

Conclusions
We identified a strong influence of interspecific variation on plant community composition and structure in the restored marsh at Bayou Dupont in southeast Louisiana. The planted sites became dominated by P. australis, likely a result of the unintended planting of an invasive lineage. The presence of invasive P. australis, which is considered undesirable by many ecologists, makes forming conclusions concerning the efficacy of planting problematic. We do not know if the positive results of rapidly-increasing cover and increased surface elevation would have occurred in the absence of invasive P. australis. Cover at not-planted sites was equivalent to that at reference sites after two years, so the labor-intensive and costly practice of transplanting rooted plants may be considered unnecessary. In addition, dissimilarity of plant communities between the treatments may be a concern. The primary species that naturally colonized the sediment-addition sites, S. alterniflora, had relatively low cover at the reference sites. It is probable that competitive interactions will lead to species replacement over time and increase the similarity between not-planted and reference sites; this trajectory, however, is contingent on marsh elevation gain exceeding the rate of sea-level rise. In contrast, the dominance of invasive P. australis at planted sites precludes the possibility that these sites will transition over time to more closely resemble the species assemblage at reference sites. Finally, geomorphic setting may be an important consideration when determining if planting is advantageous. Plantings that rapidly increase cover and increase elevation may serve a critical role in wetland restoration if natural colonization is slow under conditions of high subsidence and rising sea levels, such as occur in coastal Louisiana.