Twenty-Eight Years of Plant Community Development and Dynamics in the Balize Mississippi River Delta, Louisiana, USA

: Deltaic landscapes go through cycles of birth, growth, decline, and death governed by intertwined geological, biological, and ecological processes. In this study, we tracked deltaic lobes in the Balize Mississippi River Delta, Louisiana, USA, over 28 years (1984–2012). Hydrologic and geomorphic patterns as well as sustained patterns of wetland plant richness, diversity, and biomass are described. Plant diversity and biomass were modeled by nMDS ordination. Taxa (53) were harvested and dried (116,706 g) from 965 (0.25 m 2 ) plots and divided into three groups: I. four foundation species, corresponding to 78.9% of the total harvest; II. nine pioneer species, corresponding to 13.6% of the total harvest; and III. all other taxa, corresponding to 7.5% of the total harvest (eight miscellaneous grasses, eight miscellaneous sedges, and twenty-four miscellaneous herbs). Autogenic/allogenic processes (sedimentation, subsidence, plant colonization, and succession events) affected composition and biomass. Eleven important species were identiﬁed. Taxon richness increased on mudﬂats during primary succession (ﬁfteen to twenty-ﬁve taxa per site), then declined to fewer than ﬁve taxa per site. The niche space theory explained patterns of community change, with a similar total biomass/yr (~500 g/m 2 /yr) at all study sites. Quantile regression analyses showed that the water quality and quantity of the Mississippi River inﬂuenced biomass, especially in springtime waters. Stochastic events (storms, herbivory, salt burn, and ﬂood pulses) impacted biomass. Long-term studies like this are required in a future of climate unknowns.


Introduction
Deltaic lands are intriguing ecological environments comprising only 0.57% of the world's land area [1].Such lands are wetland ecotone environments between terrestrial and fully aquatic environments [2].For such wet terrestrial ecosystems, primary plantsuccession begins with aerial biotic life on virgin mud substrates.Their vegetative life is forced, especially by abiotic filtering, as the ecosystem assembles (sensu lato [3,4]) into communities shaped by substrate feedback over an extended period of time.Delta environments are existentially water-driven [5] and especially dynamic because the life-supporting river's qualities are far from static [6].The details of how deltaic wetland plant communities are created, maintained, and even environmentally impacted by their river source are particularly relevant to their ecology [7][8][9].Deltas are highly threatened [10].Knowledge about the parameters central to their nature is important for their protection [8], and gaps in such knowledge represent a restorative impediment [7,9,10].Humanity has generally impacted deltas in two opposing ways: by upstream deforestation, adding erosional sediment, and sediment deprivation by damming, which drives land loss [11].However, a few of the largest deltas contribute to global net land gain [11].Significant increases in river discharge also caused by human activities have had impacts [12] that continue to this day.
Many river deltas exist worldwide, both freshwater and coastal/estuarine [10].Deltas form when sediments are dropped as a river's water slows upstream from its terminus; more so in its receiving basin.Deltas can be defined using a suite of major abiotic attributes, including size, age, discharge, channel sinuosity, substrate quality, and elevation, all of which are seasonally and spatially variable.An obvious geological feature of deltas is that the dropped sediment creates land that is, at times, suitable for plant colonization to eventually form plant communities.Deltas vary considerably worldwide, but all have cycles with births, lives, and deaths dependent upon the geology of the local environment [13].Delta plant communities support a food chain that the river translates into systems of high productivity.Beyond the overall ecological importance of deltas, humanity has ties to them.Deltas are buffers for coastal marine storms and transportation byways, and they create human food.Deltaic land has considerable worth to humans [14].
One of the largest deltaic areas of coastal land in the world, created by the Mississippi River (MR), is within the State of Louisiana, USA (Figure 1).Much of the theory and modeling of a general delta cycle from geological, biological, and human perspectives is based on an enormous synthesis of work on such lands [15,16].Large coastal deltaic landscapes consist of multiple large land lobes, their sublobes, and even interior splays (essentially mini deltas) and go through cycles created by the accumulating sediments [13].Individual lobes supporting wetlands are established, grow, are maintained, and retrograde to open water as the receiving pond's elevation gradients decrease due to a lack of sediment buildup and subsidence.The same process occurs at multiple temporal and spatial levels (Figure 1).With the reduction in gradient, the sediment supply channels to individual splays shift location through small (a few meters across) "crevassing" events to events of river scale (avulsion), moving towards a steeper, more energy-favorable gradient.The cycle repeats over wide-ranging time scales: many centuries for large deltas, when the main river radically changes course many times [17,18], or over considerably less time, spanning just decades, years, months, or portions thereof, for a single lobe or a splay of multiple lobes.
To fully describe the details of the wetland complexes associated with any deltaic region like the MR Delta coastal area requires a blend of basic knowledge of both geological and botanical processes.The geology of the region concerns the balance between periods of sediment delivery (increased elevation, i.e., accretion) [19] and periods of its absence (decrease in elevation), coupled with persistent shallow belowground to deeper platform subsidence [20].The durations between accretion and subsidence in shallow estuarine waters produce mud platforms for plant communities to take hold according to their particular elevation requirements [20].Below the threshold is a mudflat, whereas above it is a wetland.Wetland types are categorized by both duration and plant makeup, as dictated by the abiotic and biotic environment [20][21][22].
The establishment and maintenance of vegetation on new wetland substrate platforms are known to be driven by both allogenic and autogenic processes, which often act in concert, particularly during early succession and even over longer periods of considerable vegetation stability [23].Riverine wetlands are especially known for the dominance of allogenic processes because of the river's dynamic nature and the fact that patterns of seasonality force vegetation to respond, although the local vegetation (autogenic processes) often plays a minor but still important role [24].The mud (mineral and organic) platform substrates that support the emergent growing-season vegetation are defined as a biotic attribute [20].Clear allogenic processes are often visible in large rivers from "pulsing events" [25,26], as demonstrated by the river forcing a sudden vegetation regime shift, a "change-in-state" event from a large pulse of flood waters with an often-linked sedimentation event [27].Riverine wetlands can even return to an earlier vegetation stage after such an event, especially after a significant disturbance.sedimentation event [27].Riverine wetlands can even return to an earlier vegetation stage after such an event, especially after a significant disturbance.The entire Louisiana coast, the coastal plain, and its 3 geomorphological west-to-east divisions-the Chenier Plain [28]; the Delta Plain [13]; and, lastly, the easternmost Pearl Plain (advocated for in [29]).The largest is the Delta Plain, made up of seven geologically distinct deltaic regions formed sequentially (youngest to oldest: "A" through "G", with "B" essentially an extension of delta "C" (Plaquemines Delta)) over the past ~7500 yrs [30].The newest delta region (A) in the Delta Plain is the region around the terminus of the Atchafalaya River, which is composed of both the Atchafalaya and Wax Lake Deltas [31,32].The Balize Delta (B), the focus of this study, is shown in Panel II (modified from [33]).The Balize Delta began its growth ~550 yrs ago and is composed of four distinct, smaller-scale subdeltas identified sequentially by age ("a, b, c, and d").Panel III: Embedded within the four subdeltas of the extant Balize Delta are many smaller "inner deltaic splays", filling in shallow ponds, indicative of a senescing Balize Delta.Panel IV: Inner deltaic splays are "mini deltas", each composed of multiple lobes; each lobe likewise has its own birth, growth, and decline phases.The entire Louisiana coast, the coastal plain, and its 3 geomorphological west-to-east divisions-the Chenier Plain [28]; the Delta Plain [13]; and, lastly, the easternmost Pearl Plain (advocated for in [29]).The largest is the Delta Plain, made up of seven geologically distinct deltaic regions formed sequentially (youngest to oldest: "A" through "G", with "B" essentially an extension of delta "C" (Plaquemines Delta)) over the past ~7500 yrs [30].The newest delta region (A) in the Delta Plain is the region around the terminus of the Atchafalaya River, which is composed of both the Atchafalaya and Wax Lake Deltas [31,32].The Balize Delta (B), the focus of this study, is shown in Panel II (modified from [33]).The Balize Delta began its growth ~550 yrs ago and is composed of four distinct, smaller- The actual vegetation dynamics of splays or any deltaic complex, from wetland birth to death, have never been fully tracked in detail.However, the importance of sediment budgets (the overall elevation dynamics) for wetland growth, maintenance, and dynamics in young Louisiana deltaic landscapes has been widely studied for years [34].The growing Atchafalaya Delta [35], the adjacent Wax Lake deltas [32], and the much smaller inner Mississippi Balize Delta lobes/splays [36] all show similarities in wetland plant establishment and general community dynamics.
The overarching complex question addressed by this unprecedented long-term study is: How do plant communities assemble themselves from birth (colonization), early development, and establishment (primary succession), through a relatively stable existence, to their decline and disappearance (death)?This complex question was answered by a series of incremental questions: Are there common patterns of plant species richness among different sites?Do patterns exist in the annual study-wide plant community's productivity according to particular species or groups thereof?Do the river water attributes affect the deltaic wetland's quality and quantity?Can infrequent stochastic events have significant impacts on the phytosociology of a delta?Tied to these questions, significant contributions of specific plant species during the delta cycle were identified.

Study Region and Sites
This field study was designed close to 4 decades ago as a simple descriptive project focused on wetland establishment on new mud substrate platforms within the Balize Delta.Such substrates were particularly rare at the time in the far greater Louisiana coastal landscape experiencing wetland loss [37].Federal efforts to increase wetland areas had begun, and some had been completed [38].Within the Balize Delta, an area of numerous splays of different ages, sizes, and origins was identified in 1984 to monitor wetland plant colonization, community establishment, and community change over time (Figure 1) [36].Each splay developed down-current from crevasses along the area's many kilometers of shoreline.Lobes within these splays have a considerably attenuated elevation gradient down <1 m to the platform surface (Figure 2) [20].In addition, lobes decline in elevation perpendicular to the shore of the sediment supplying feeder pass and are typically assumed to be delineated at the emergent vegetation boundary, although the lobe's substrate continues to slope towards the bottom of the receiving pond, often over considerable distance [39].The entire elevation change across a lobe only in the vegetated areas is rarely more than 50 cm.Because splay lobes can be 1 km or more in length and usually not more than 100 m in width, the elevation gradients are imperceptible without biotic markers.Therefore, defining the precise width of a lobe is difficult because of its gradual transition to unvegetated flooded muck.The vegetated width of lobes grows far slower compared to the length axis because of the location of the more stable and coarser (sandier) sediment along the feeder pass edge.
The largest splay in all of the Balize Delta was one chosen for our study: the Brant Splay (BS).This area was documented by substantial NASA aerial photographic records as a shallow pond prior to 1978, the year its crevasse opened along the eastern shore of Brant Pass.MR water has flowed into this large (~10 km 2 ), shallow receiving pond.The cyclical vegetative nature of the area is demonstrated by the fact that in the early 20th century, the exact area supported freshwater wetlands described, drawn, and pictured [40] as vast areas of Typha spp.The exotic, Alternanthera philoxeroides was reported (authorship for all species collected in the study is provided in Table 1), and in the interspersed lower, more watery areas, Eichornia crassipes (Mart.)Solms was reported.The lands were all bordered by higher ridge/natural levees of Panicum repens L., Zizaniopsis miliacea, and Salix nigra.Interspersed within the wetlands were significant ponds-remnants of the original vast sediment-receiving areas.This historical wetland disappeared because of reduced sediment input coupled with the generally high subsidence rates [33,38].The White's Splay (WS) area began as a linear region of mudflat associated with a man-made (1982) crevasse along Octave Pass.The artificially created (not natural) mudflat was at grade with neighboring wetland areas northeast of the BS in the receiving pond, into which the dredging extended ~500 m away from the pass as a channel to accept sediment-laden water.The event engineered by the United States Fish and Wildlife Service was the first major attempt to build significant wetlands within Louisiana's Delta Plain by construction of an artificial crevasse [38].
In 1984, three sites were chosen for study.Two sites were established on 2-year-old (determined from NASA photographic records) rapidly elongating lobes of the BS: a north site (BSn) and a south site (BSs; Figure 1, panel IV).A third site was established ~150 m inland from Octave Pass on the east side of the southeasterly dredged channel.This  20] thanks to J. Lynch) with elevation change and spatial relationships depicting the spatial and height relationships of the wetland habitat on its platform range of <1 m from the shallow pond bottom (depth of ~50 cm) at the lobe's low end of the emergent mudflat (+17 cm) to the low marsh (+26 cm), the high marsh (+30 cm), and the most upland Salix nigra (SN) forest (>6 cm; mean of three laser level readings).Thus, the flood tolerance for the four predominant community-identifying species (SN > SD > SL > SP) extends through only ~62 cm elevation, regardless of the lobe's length.A lobe's length can range from 50 meters to kilometers, as in the case of the BSs lobe.Key: SN = Salix nigra; SD = Schoenoplectus deltarum-dominated marsh; SL = Sagittaria latifolia low marsh; SP = Sagittaria platyphylla lowest marsh; and SAV = submerged aquatic vegetation.The White's Splay (WS) area began as a linear region of mudflat associated with a man-made (1982) crevasse along Octave Pass.The artificially created (not natural) mudflat was at grade with neighboring wetland areas northeast of the BS in the receiving pond, into which the dredging extended ~500 m away from the pass as a channel to accept sediment-laden water.The event engineered by the United States Fish and Wildlife Service was the first major attempt to build significant wetlands within Louisiana's Delta Plain by construction of an artificial crevasse [38].
In 1984, three sites were chosen for study.Two sites were established on 2-yearold (determined from NASA photographic records) rapidly elongating lobes of the BS: a north site (BSn) and a south site (BSs; Figure 1, panel IV).A third site was established ~150 m inland from Octave Pass on the east side of the southeasterly dredged channel.This constructed or created WSc site was originally a designed test location for wetland development to monitor the two natural sites on BS.Years later, another lobe showed potential for growth on the same growing splay down from WSc but was not influenced by the dredged material from a decade before.In 1994, the area was designated as the new WSn site (Figure 1, Panel IV) and subsequently grew southwest, splitting off the then-most distal end of the WSc's southernmost, still growing lobe.Its sediment filled more of the northern portion of the same large receiving pond south of both Brant and Octave Passes.A few years into the growth of the WSn lobe, a declining phase began, caused by channel sediment backfilling enabled by high statured wetland plant invasion into the shallow feeder channel, considerably reducing the flow of water (see below).The elevated and most proximal mud substrates of the WSn lobe initially supported wetland communities similar, when of the same age, to those on the two natural BS lobes in previous years.
At BSs, BSn, and WSc in 1984 and at WSn in 1994, study site transect locations for plant collection on the newly emerged sites were marked by permanent stakes in the mid-summer near the shore at the most-distal growing ends of each site's lobe.The transect locations were down gradient from earlier deposited aerial land substrates, which supported more mature wetlands on the same lobe.The area of interface elevation between the previously deposited substrate of new mud arises when late-spring highly sediment-laden MR water recedes to expose the accommodating substrate for plant growth.The growth rates of the splays were determined in previous years from birth to 2000.With each year of substantial accretion (which forced elongation of the lobes), additional transects were added: two by 1986 at WSc and BSn and two by 1989 at BSs.All lobes had to have elongated ~50 m or more for a new transect.No further transects were deemed necessary after sufficient plots were established to describe the lobe plant community (see below; usually three transects with five plots each).
Later, in 1994, at the start of the WSn study, two transects were placed on substrates of two ages: a distal mudflat of just accreted land and an upslope, i.e., proximal large, area observed to have emerged the previous summer and colonized later that summer by a typical pioneer plant community similar to those found at the other sites [36].The WSn lobe grew no further after the establishment year.Coincidentally, the relative homogeneity of the developed wetland community meant that only two transects of harvest data were sufficient to be fully representative of the area's wetlands.
At the time of usual flowering (late August to early September), the peak plant standing crop (referred to as aboveground biomass (AGB)) along each transect was clipped in five 0.25 m 2 (50 cm × 50 cm) plots ~5 to 10 m apart.The distance between plots was gauged by lobe width coupled with the diversity of the immediate community so as to represent the plants of the general transect area.For a few transects in a few years, the widths of the wetlands across a lobe increased further into the downslope mud (muck) area at their far edge which supported sparse vegetation (Figure 1; Note 2).Throughout the study, the distant edge community was always of the same species as the most inland harvested transect plot.With confidence that the harvested plants along all transects represented the majority of the wetland vegetation of the complete lobe area each year (the specific plot locations were determined by D.W. for all years).
The harvested AGB of plants/plots were sorted to the lowest taxonomic group possible, i.e., species, in all years except the summer of 1984.During the year of site establishment that was not sampled, the mudflats at the first three sites were not recognized as first-summer wetlands because they were sparsely vegetated.The fortuitously captured 1984 and 1985 site photographs confirmed the 1984 areas as nascent months-old wetlands.Therefore, additional plant harvesting took place in 1985 from similar nascent areas at the three sites to capture the early wetland AGB total, although not at the taxon level, then hind-casted for the assumed 1984 total.
All plant material, including unidentifiable dead material, was dried to constant weight and entered into Excel spreadsheets.Voucher material was deposited in the herbarium of Tulane University, New Orleans, LA, USA and has since been moved to the collections at Louisiana State University, Baton Rouge, LA, USA.The transect-by-transect, plot-by-plot data were used to document and decipher the wetland community events over the 28-year study period.Events included plant colonization; primary succession; and other significant phytosociological processes, such as declines in coverage, wetland disappearances, and plausible causes.

Study-Wide Harvested Vegetation
Over the course of the study, a total of 53 species or taxonomic groups were harvested, of which 38 were assigned to species (Table 1).All taxa were divided into three major groups: foundation species (FS [41][42][43], four taxa), pioneer species (PIO, nine taxa), and the remaining other species and groups thereof (OTH, 40 sorted into three subgroups: 8 Msc.grasses, 8 Msc.sedges, and 24 Msc.herbs).The 40 other species are a collection of the leastabundant taxa study-wide.The FS are Salix nigra, Sagittaria latifolia, Sagittaria platyphylla, and Schoenoplectus deltarum.Salix nigra was the only fully woody species harvested (only as saplings) study-wide.
The group labeled PIO included Ammannia coccinea, Bacopa monnieri, Cyperus difformis, Cyperus erythrorhizos, Echinochloa walteri, Eclipta prostrata, Leptochloa panicoides, Panicum dichotomiflorum, and Sphenoclea zeylanica.The PIO species were chosen generally using a set of criteria including site abundance (AGB) and frequency (evenness, an aspect of diversity sensu [21]) in the harvested plots during the primary succession year (PSY) of each site following mudflat emergence [36].The criteria were (a) significant AGB at one site and present in at least three sites; or (b) significant AGB at two sites and present in another site over those PSYs.The total number of PSYs harvested study-wide was 10: 3 years each for the early established sites (WSc, BSs, and BSn) and a single year for WSn.Panicum dichotomiflorum was an exception because analyses showed it should be designated a PIO species (see below).

Analyses of Splay Lobe Phytosociological Dynamics
First, the major vegetation change at each of the four study sites (lobes) was detailed over the decades by scoring 965 plots according to their most dominant taxa using the AGB metric (sorted from all 53 into 11; one of the four FS, four of the nine PIO species of greatest study-wide AGB: E. walteri, C. difformis, C. erythrorhizos, L. panicoides, or one of the three Msc.groups, i.e., herb, grass, or sedge).The five least abundant PIO species were placed accordingly into their Msc.group, including the study-wide AGB of all species.Only in 23 of the 965 plots was one of the five least abundant PIO species the most dominant in the plot, whereas among 28 other plots, only eight of the other rare taxa within the Msc.groups dominated (see below).Major splay-wide community patterns (including subtle patterns) of species' presence or absence over time were captured by this method and proved important for the description of the dynamic nature of the MR Balize Delta wetlands.
Nonmetric multidimensional scaling (nMDS) statistical analysis for large datasets [44] was used to further examine vegetation changes.The Bray-Curtis similarity index was applied, with 50 iterations and 1000 permutations used to create a fully integrated ordination of all lobe plant communities over time using species AGB data from the 965 harvested plots.A trajectory ordination plot of the year-by-year vegetation change was then created for each lobe under study in the Balize Delta.
In addition to the nMDS procedure, a distance-based linear model (DistLM [44,45]) was run to determine which species ordinated as statistically significant (p > 0.05) over a complete cycle of the Balize Delta lobe.DistLM is a non-parametric permutation routine used to sort out relationships within a complex multivariate constellation (in this case, comprising the 53 harvested taxa).The marginal test output identifies the taxa of statistical importance to the recorded vegetation change over time (in this case, decades).
Species richness was determined by counting the total number of species in the plots harvested from each lobe.Trends over time were determined using linear regression.
Interlobe patterns common to lobe-wide wetland plant communities and their total aboveground biomass (TAGB) levels were determined by calculating the average AGB for each species on each lobe across all years.To reveal broad-scale similarities and differences in the vegetation communities over time, i.e., annually in each studied wetland, the TAGB values from the previous analysis were parsed by the four FS, the two taxa groups (PIO and OTH species), and the TAGB.Additionally, the yearly totals over common times were statistically compared using correlation analysis between sites to learn whether patterns existed irrespective of the vegetation quality.
Using available high-quality aerial NASA photographs, which covered most of the study period (1983, 1995, and 2000), splay-wide wetland increase was determined using Geographic Information System software (GIS; ArcView 8.3, ESRI-Environmental Systems Research Institute, Redlands, CA, USA).The area coverage of the two prominent wetland regional vegetation types that make up the splays (forest and herbaceous) were calculated to determine rates of growth per year [46].
Finally, quantile regression was employed to address the following question: Does the MR's water quality and quantity influence the prevalence of the dominant species on the splays under study in the Balize Delta?The regressions might explain some of the large-scale patterns in the observed annual fluctuations in a lobe's wetland AGB and community change.Quantile regression analyses were performed using SAS 9.2 [47].For the analyses, the annual AGB values of the most dominant species, Sc. deltarum and Sag.Latifolia, collected from 1989 to 2008 were applied as dependent variables for the following water variables of the MR: (1) spring discharge; (2) growing season nutrient (nitrogen) load; and (3) a specific MR temperature (when cold spring river waters transition to 15 • C and above), which was converted into a metric for growing season length using GDDs (growing degree days [48]).The variables were mined from three public governmental web sources.
(1) Most of the water temperature data were obtained from 8AM readings collected by the United States Army Corps of Engineers' MR gauge in New Orleans (Carrollton Gauge; River km 165.4).The relatively few missing gaps in these data were filled as per [49].
(2) The daily discharge data were obtained from the Corps of Engineers for Tarbert Landing (River km 493), based on which the average discharge for the period of March 1 through May 31 was calculated/yr to quantify the spring flood amount.(3) The nutrient load data were obtained from the Louisiana Department of Environmental Quality from three long-established sites transecting the MR at Pointe a la Hache, LA (River km 78.4).The obtained variables were nitrate + nitrite and TKN monthly concentrations.The average monthly concentrations (ppm) were multiplied by the average discharge (cfs) to obtain a monthly loading.The full summer loads were then calculated as the sum of the average monthly loads from April to August.
The second group of nine PIO species accounted for only 13.6% of the total harvest.Four of the five PIO graminoids with the highest harvested AGB in this group study-wide (E.walteri, C. difformis, C. erythrorhizos, and L. panicoides) made up 82.4% of the total PIO species.The 12 herbaceous species from the FS and PIO groups (excluding woody FS Sal.nigra) accounted for 90.1% of all harvested AGB (including Sal.nigra, the AGB equals 92.4% study-wide for all FS and PIO).PIO E. walteri, the third most abundant species in this study, had a marginally higher harvested AGB than FS Sag.platyphylla.
Of the last of the three groups, OTH species (eight grasses, eight sedges, and 26 herbs; 42 in total), all but two fully herbaceous species, except the semi-woody Baccharis halimifolia and Sesbania drummondii, were considerably rarer, accounting for only 7.5% of all plants harvested.The three most abundant taxa in OTH species, all Msc.herbs identified at the species level, were Vigna luteola, Polygonum punctatum, and T. domingensis.The combined AGB of the herbs equaled 42.3% of the OTH group's total but only 3.2% study-wide.A total of 15 taxa included in OTH species listed in Table 1 were not identified to the species level (three grasses, four sedges, and eight herbs = only 3.1% of all harvested; and a dicot group of 3 = <0.3%AGB).The combined AGB of two taxa in this 15-member group (Panicum spp.and Msc.Poaceae) made up 68.7%; thus, the 13 remaining taxa were exceptionally rare.Only 4 of the 40 taxa in the OTH group (Pluchea odorata, Po. punctatum, T. domingensis, and V. luteola) became important in DistLM (see below) because of their relatively large total AGB; two PIO species had a lower study-wide AGB than these four OTH species.

Plot Taxa Dominance over Time
The principal inter-and intra-lobe vegetation differences and variation over time were revealed by our analysis.The harvested AGB captured the vegetation diversity, even within the greater wetland area of each plot (Figure 3), because of the attention focused on the plot's precise field location along each transect, as previously described, and the relative homogeneity of the vegetation.The three herbaceous FS were the most prevalent on all the splay lobes under study and tended to follow a pattern of increasing dominance after the first year of study (see below).Although several plots are dominated by one species, biomass of other species is still present (Figure 4).A clear example occurred in 2005 at WSn, where four plots contained some SD biomass, but SL was dominant in 100% of the plots.
at the species level, were Vigna luteola, Polygonum punctatum, and T. domingensis.The combined AGB of the herbs equaled 42.3% of the OTH group's total but only 3.2% study-wide.A total of 15 taxa included in OTH species listed in Table 1 were not identified to the species level (three grasses, four sedges, and eight herbs = only 3.1% of all harvested; and a dicot group of 3 = <0.3%AGB).The combined AGB of two taxa in this 15-member group (Panicum spp.and Msc.Poaceae) made up 68.7%; thus, the 13 remaining taxa were exceptionally rare.Only 4 of the 40 taxa in the OTH group (Pluchea odorata, Po. punctatum, T. domingensis, and V. luteola) became important in DistLM (see below) because of their relatively large total AGB; two PIO species had a lower study-wide AGB than these four OTH species.

Plot Taxa Dominance over Time
The principal inter-and intra-lobe vegetation differences and variation over time were revealed by our analysis.The harvested AGB captured the vegetation diversity, even within the greater wetland area of each plot (Figure 3), because of the attention focused on the plot's precise field location along each transect, as previously described, and the relative homogeneity of the vegetation.The three herbaceous FS were the most prevalent on all the splay lobes under study and tended to follow a pattern of increasing dominance after the first year of study (see below).Although several plots are dominated by one species, biomass of other species is still present (Figure 4).A clear example occurred in 2005 at WSn, where four plots contained some SD biomass, but SL was dominant in 100% of the plots.Figure 3 shows the appearance and subsequent dominance of Sal.nigra in only the WSc plots.In the early years, Sal.nigra grew to seedlings only, then disappeared, returning the following 4 years to grow into a wetland of predominantly small trees.At the three initial sites (WSc, BSs, and BSn), our analysis showed general common colonization and primary succession vegetation during the PSYs based on plot presence during the 3 years after the mudflat's emergence.
All three sites had the same number of taxa (five) from the seven key PIO species but differed slightly depending on the community.All included C. erythrorhizos, L. panicoides, and Msc.grasses, but E. walteri, C. difformis, and multiple Msc.sedges were found at a single site only.At WSn, at the end of its third growing season (the first year of sampling), 50% of the lobe's harvested plots along the most distal transect were dominated by just two PIO species (C.erythrorhizos and L. panicoides), and like the other sites during their third PSY, the FS were dominant in the community.
Abundance is more descriptive of the early wetland communities than presence during the PSYs (Figure 4).The top five most abundant of the nine PIO species from the 119 harvested PSY plots were C. erythrorhizos (total = 357 g), L. panicoides (301 g), C. difformis (174 g), E. walteri (140 g), and Sphenoclea zeylanica (114 g) (in order of mean AGB/m 2 ).The total harvested AGB of these species over the PSY (9 combined years) amounted to more than a quarter (27.7%) of the AGB collected across all sites for the 50 combined study years (Table 1).The first-year substrates supported mostly PIO species (75.3% of the total AGB), with a large mix of very few OTH species (6.8%).More specifically, for the first year of the PSY at these sites, the five species averaged 67.8% of the total AGB on the then-nearly bare mudflats.The same year, the last four PIO species (Bacopa monnieri, Panicum dichotomiflorum, Ammannia coccinea, and Eclipta prostrata) made up just 7.5% of the mean site total AGB, and the first two were not harvested at every site.Our validation of early species colonization showed that two members of the PIO group (C.erythrorhizos and S. zeylanica) were most commonly observed mid-summer on newly emerged splay lobes well before their usual harvest in August.
Our analysis also showed that the four PIO species with the largest total studywide AGB, on occasion, reappeared post-PSY (PPSY) at the three initial sites (Figure 3): C. difformis, particularly in 1993 at BSs in a significant eruption and with a smaller eruption in 1988 at BSn; E. walteri similarly at BSs in 1995/2004; L. panicoides, with a large increase at BSs in 1992; and C. erythrorhizos at the WSc site in 1994.In 2003, Msc.herbs at both BSs and BSn showed considerable eruptions, as well as periodic lower levels during other years at the three sites.Msc.grasses dominated the PSYs in a few plots only at WSc and BSn and for a few PPSYs at BSs. Lastly, Msc.sedges were rare in terms of plot dominance and during isolated PPSYs, which were mostly at WSc but also at BSn (Figure 3).
The percentages of plots/site/study PPSYs that were dominated by species other than FS are listed as follows (from high to low): WSc = ~19.7%;BSs = ~17.4;BSn = ~11.7;and WSn = 0. Sagittaria platyphylla was particularly found along lobe edges (~35% of the edge areas), and Sag.latifolia was often observed in large monospecific stands across the wider lobe at its adapted elevation.The other 65% AGB of Sagittaria platyphylla was observed in areas of either new mudflat growth where seeds were plentiful, well removed from upland species' seeds, or areas of the same quality with reversion to the lower elevation threshold.The increase in the mudflat edge at the site also became a niche for Sag.latifolia on more firm edge areas of its adapted elevation space.Such edge morphodynamics are an attribute of deltas [50].

Plot Integration over Time
The nMDS analyses created trajectories of vegetation change across and within the wetland lobe areas, which confirmed the common patterns of change on Balize Delta splays.Analyses of variances (ANOSIM) for our large, study-wide efforts first revealed no major plant community differences among years or sites within the created species cloud (ANOSIM: p = 0.323).When the nMDS was run separately for each site, some significant differences were revealed (p = 0.029) in plant communities among sites at the same sampling site/time, especially among years (p = 0.004), as shown in the overlay of the four-site ordination model of trajectory plots presented in Figure 5, but in only two-dimensional space.The four plots along the single x-axis space are a composite of the three phases of vegetation change during the typical delta lobe cycle.Ordered from right to left along each third of the axis are (1) the start of colonization on nascent mudflats through the primary succession autogenic period; (2) a longer, much more dynamic community change period; and (3) a senescent phase culminating in wetland loss at only the site (WSn) with lands reverting to Mudflat.
portant five-plant taxonomic axis space (Figure 5) because of the known decreasing flood tolerance rankings for the four FS: SP > SL > SD > SN.As well as, the combined nine PIO species implied along the axis (from right to left) over the decades of vegetation change.Two other phytosociological implications demonstrated in the two dimensions of the figure can be described.First, Sal.nigra occurred usually as seedlings in plots associated with the PIO species during the first one to two years among the PSY community; if not, during PPSYs, Sal.nigra occurred on the higher substrates, usually as seedlings.Second, the trending community trajectory for both the WSc and BSs sites (the two most elevated sites) ended close to the same vegetation space (indicated by the arrows) because of the similar vegetation state-change, which resulted in the termination of the field sampling (coincidentally, ten years apart: 1998 and 2008).The two-dimensional space of the nMDS trajectory became definable by the important five-plant taxonomic axis space (Figure 5) because of the known decreasing flood tolerance rankings for the four FS: SP > SL > SD > SN.As well as, the combined nine PIO species implied along the axis (from right to left) over the decades of vegetation change.Two other phytosociological implications demonstrated in the two dimensions of the figure can be described.First, Sal.nigra occurred usually as seedlings in plots associated with the PIO species during the first one to two years among the PSY community; if not, during PPSYs, Sal.nigra occurred on the higher substrates, usually as seedlings.Second, the trending community trajectory for both the WSc and BSs sites (the two most elevated sites) ended close to the same vegetation space (indicated by the arrows) because of the similar vegetation state-change, which resulted in the termination of the field sampling (coincidentally, ten years apart: 1998 and 2008).
A DistLM marginal test identified the 11 top species jointly on the four lobes of significance based on their high AGB among the 53 total taxa (p > 0.05; three herbaceous FS, six PIO, and two OTH; Table 2) using the ordination's proportion rankings.The rankings define the structural change in the lobe vegetation during the decades of the study of the delta cycle.If only the face value of the AGB metric (abundance) was used in the analysis, taxa with high study-wide harvest values would with statistically significant importance commensurate with the Proportional metric; however, few did in the final derived constellation cloud.The best example was the third most abundant taxa collected and the most abundant PIO species: E. walteri (5566 g; Table 1), with a p-value below that of less abundant taxa.Incongruence was also observed in Sal.nigra, the fourth FS not listed among the 11, although it was the sixth most abundant taxon.Sal.nigra's p-value = 0.146, which is far below significance, likely because of the presence of its seedlings during PSYs and many PPSYs, not unlike the vast majority of the non-PIO herbs within the explanatory Water 2023, 15, 3481 14 of 29 model.Table 2 shows that only the top two species' AGB metric rankings (Sc.deltarum and Sag.latifolia) align with the proportion and p-value listings, the rankings of which fully match.1).d Presence defined as percentage of the 965 plots harvested in which the species was present.e Ranking based on the sum of the rank for abundance and presence.f Dominance defined as the percentage of the 965 plots harvested in which the species was dominant.g Ranking based on the sum of the rank for abundance, presence, and dominance.h Two species had equal values; the highest ranking was assigned to species with the highest TAGB.
The DistLM analysis is mostly reflective of the harvested data for each of the taxa, as well as their study-wide and site-wide abundance, but the AGB metric is apportioned with additional phytosociological data according to year, site, and plot through two other attributes: a taxon's plot-by-plot presence over wetland area (its ecological evenness) and plot dominance among the 965 total tabulated plots.The combination of these three metrics (abundance, presence, and dominance) summarizes the relative importance of the taxa in the ordination, particularly for those reaching statistical significance.
The three herbaceous FS were highly present across plots and evenly spread over the sampled lobes: Sc. deltarum, with a mean presence (among the four lobes) = 75.1%;Sag.latifolia = 34.9%;and Sag.platyphylla = 9.2%.Salix nigra had a mean presence of 13.5%, which is below that of two of the top-four-ranked PIO species (E.walteri = 14.0% and C. difformis = 22.8%;Table 1).Most of the other taxa that were ranked proportionally high in the top 20 (not 11) DistLM marginal test results (Po.punctatum, T. domingensis, Pl. odorata, Panicum spp., Poaceae, or Cyperaceae-all with a presence > ~5.0) were so only during short eruptions over decades and therefore did not exhibit significant community composition change.
The metric plot dominance (Table 2) is an attribute that represents the area commonality among the taxa-a concept concerning the vegetation in the vicinity of each plot (rather than simply the abundance) but still including each plot.Such a metric is especially informative because each plot was locally area-representative according to its selection criteria (see Methods).
The three herbaceous FS were understandably highly plots, i.e., the lobe wetland area was dominant: Sc. deltarum was found in 53.1% of the 965 plots in the four lobes (the largest proportion); Sag.latifolia = 18.7%; and Sag.platyphylla = 8.2% and nearly fully dominant at WSn in its last 6 yrs (52 of 60 harvested plots = 86.7%).Salix nigra was found in 1.5% and only dominant at WSc, below the mean dominance of two of the top-four-ranked PIO species (E.walteri = 2.4%; and C. difformis = 2.9%).Sc. deltarum and Sag.latifolia are clearly predominant species on the study lobes (Table 2).

Lobe-Wide Richness Patterns
At the initial three study sites (BSs, BSn, and WSc), the species richness increased considerably between the first and second years of the harvest period (Figure 6).The richness then generally plateaued for a few years near the years 12-16, except for two significant short-term shifts: a large one-year spike to 22 species in 1988 at WSc, and a 2-year dramatic decline at BSn in 1989/1990.A large initial increase in richness at the final site (WSn) in its third year of vegetation establishment (1994) would have been recorded had it been monitored for all of its PSYs.At WSn, no plateau in richness occurred as at the other sites, but a 2-year decrease to a plateau level of +five species through 2008 was observed, followed by a decline to a monospecific stand of Sag.platyphylla in its final year before converting to open water.
equal values; the highest ranking was assigned to species with the highest TAGB.
The three herbaceous FS were understandably highly plots, i.e., the lobe wetland area was dominant: Sc. deltarum was found in 53.1% of the 965 plots in the four lobes (the largest proportion); Sag.latifolia = 18.7%; and Sag.platyphylla = 8.2% and nearly fully dominant at WSn in its last 6 yrs (52 of 60 harvested plots = 86.7%).Salix nigra was found in 1.5% and only dominant at WSc, below the mean dominance of two of the top-four-ranked PIO species (E.walteri = 2.4%; and C. difformis = 2.9%).Sc. deltarum and Sag.latifolia are clearly predominant species on the study lobes (Table 2).

Lobe-Wide Richness Patterns
At the initial three study sites (BSs, BSn, and WSc), the species richness increased considerably between the first and second years of the harvest period (Figure 6).The richness then generally plateaued for a few years near the years 12-16, except for two significant short-term shifts: a large one-year spike to 22 species in 1988 at WSc, and a 2-year dramatic decline at BSn in 1989/1990.A large initial increase in richness at the final site (WSn) in its third year of vegetation establishment (1994) would have been recorded had it been monitored for all of its PSYs.At WSn, no plateau in richness occurred as at the other sites, but a 2-year decrease to a plateau level of +five species through 2008 was observed, followed by a decline to a monospecific stand of Sag.platyphylla in its final year before converting to open water.The 1994 richness at WSn was high (16 species) compared to the other sites during the initial year of harvest because WSn's first year of harvest was year 3 of the wetland's existence (see text), i.e., the year understood to be the last of the primary succession period [36].The 1994 richness at WSn was high (16 species) compared to the other sites during the initial year of harvest because WSn's first year of harvest was year 3 of the wetland's existence (see text), i.e., the year understood to be the last of the primary succession period [36].
Therefore, except for a few short-duration sharp increases in richness during the PSYs, the overall mean study-wide and across the four sites declined significantly (p < 0.001).However, BSn exhibited a precipitous decline in 1989, which resulted from the salt burn impact from tropical winds and an extensive herbivore eat-out, both with likely consequences lasting through collection time the following year.

Lobe Total Biomass Patterns
On average, a remarkably similar TAGB (~500 g/m 2 /study yrs) was harvested at the three initial sites, i.e., WSc, BSs, and BSn (Figure 7).The TAGB of the WSn site was significantly lower over the study period (~350 g/m 2 /study yrs, p < 0.01).Importantly, when the site totals were parsed by the four FS and two other groups (PIO and OTH; Table 1), inter/intra-site patterns were found, with mean relative elevations as ranked in Figure 7 (WSc BSs > BSn > WSn).These rankings were determined by coupling the metric of FS dominance/site represented by the area coverage/abundance of the four FS (mean TAGB/total study years) with fidelity to the actual elevation (the metric for flood tolerance; see Figure 2) on the lobes.The total abundance/site of Sal.nigra in Figure 7 tracked the relative elevations of the sites.Note that at WSc, the SN's 1995-1998 data were not included at Transect 1 because of the state change to a shaded woody community of Sal.nigra plus mature individuals of the semi woody shrub Se. drummondii; thus, harvesting after the first year of seedlings meant no future AGB tracking at WSc Transect 1.At WSc, Sal.nigra was the most abundant usually as first-year seedlings (in 72 of 182 plots) and the least abundant at WSn (found in only four of its 176 plots).

Lobe Total Biomass Patterns
On average, a remarkably similar TAGB (~500 g/m 2 /study yrs) was harvested at the three initial sites, i.e., WSc, BSs, and BSn (Figure 7).The TAGB of the WSn site was significantly lower over the study period (~350 g/m 2 /study yrs, p < 0.01).Importantly, when the site totals were parsed by the four FS and two other groups (PIO and OTH; Table 1), inter/intra-site patterns were found, with mean relative elevations as ranked in Figure 7 (WSc > BSs > BSn > WSn).These rankings were determined by coupling the metric of FS dominance/site represented by the area coverage/abundance of the four FS (mean TAGB/total study years) with fidelity to the actual elevation (the metric for flood tolerance; see Figure 2) on the lobes.The total abundance/site of Sal.nigra in Figure 7 tracked the relative elevations of the sites.Note that at WSc, the SN's 1995-1998 data were not included at Transect 1 because of the state change to a shaded woody community of Sal.nigra plus mature individuals of the semi woody shrub Se. drummondii; thus, harvesting after the first year of seedlings meant no future AGB tracking at WSc Transect 1.At WSc, Sal.nigra was the most abundant usually as first-year seedlings (in 72 of 182 plots) and the least abundant at WSn (found in only four of its 176 plots).FS Sc. deltarum's abundance generally mirrored the rankings of Sal.nigra, although the total collected AGB of the species was somewhat higher at BSs compared to that at WSc (Figure 7).Holding true to the elevation fidelity, both Sag.latifolia, the FS of middle substrate elevations, and Sag.platyphylla, the FS common to the lowest lobe elevations, followed the site rankings but conformingly opposite to that of Sal.nigra and Sc.deltarum.
Not only is this fidelity to elevation validated by the FS abundances; it was a wetland splay lobe community attribute that was fully mirrored when combining the total mean/yr AGB of the PIO species of each site with that of its OTH species.The WSc site accounted for 28.6% of the total AGB from these two groups, BSs accounted for 22.9%, BSn accounted for 21.3%, and WSn accounted for 4.8%.The total AGB of these two groups individually matched the site elevation rankings, with the exception of OTH species at BSn.Therefore, lobe elevation (the degree of flooding) is an important index of gross vegetation change.

Lobe Decadal Vegetation Patterns
Tracking the inter-and intra-site AGB yearly totals (Figure 4; Table 3) revealed seven patterns across the common years of harvest of all possible pairings of sites: 1.
A highly significant correlation among the three most elevated sites (WSc, BSs, and BSn) when the three common PSYs were included.The one exception was for the 15 common harvest years between WSc and BSn-a period of no significance, although the AGB did show a correlation during the PPSY (i.e., after 1987); 2.
The lowest mean elevation site (WSn) paired with the other three sites never showed significance across the common years of study (Table 3); 3.
The most striking pattern of the annual total AGB values occurred in 1996, when all lobe sites experienced a sharp one-year increase in community biomass, particularly for the three most elevated sites (Figure 4a); 4.
Large single-year increases were common across all sites (in 1986 at BSn, 1987 at WSc, 2005 at BSs, and 2006 at WSn), with matched increases in other years (BSs-BSn in 2000, BSs-BSn in 2004), and likely through 2005 if BSn had been harvested (the channel access was blocked by more than 0.5 km of floating and aquatic vegetation).The yearly trends in the total harvested AGB followed the abundances of Sc. deltarum and Sag.latifolia (Figure 4b-e); 5.
An intrasite pattern was statistically substantiated.Over the full study period at BSn, an inverse relationship between the presence of deltarum and Sag.latifolia occurred (CC = −0.534;p < 0.05).At BSn, Sc. deltarum and the combined members of OTH species showed positive relationships over the study years (CC = 0.551; p < 0.05).An inverse significance was observed among FS at the WSn site: between Sc. deltarum and Sag.latifolia (CC = −0.545;p < 0.05) and at the same site between Sag. latifolia and Sag.platyphylla (CC = −0.584;p < 0.05; Figure 4d).Hints of this pattern were visible at WSn between Sc. deltarum and Sag.latifolia (Figure 4e) but not among herbaceous FS at WSc.Because WSc was the most elevated site, its construction likely had impacts, although Sc. deltarum was common at that site; 6.
Several FS experienced eruptions in particular years, such as the 1996 increase in Sc. deltarum at all sites and two large increases of the same species at the three sites that were then under study: BSs, BSn, and WSn in 2000 and BSs and WSn in 2006.In 2004, at the same sites, Sag.latifolia exhibited substantial increases.The fourth and woody FS, Sal.Nigra, showed a massive peak in abundance at WSc during the fourth year of study (1987, its first PPSY) and primarily as seedlings during all years of the study.During the PPSYs, Sal.nigra seedlings were harvested in only 15 plots of the 763 other sites and grew into young shading trees the last 4 years at WSc, along with the semi woody Se. drummondii (Figure 4). Figure 4 does not adequately depict the extensive diversity of understory woody to semi woody shrubs and the occasional other tree taxa.

MRBD Splay Wetland Growth
Winter NASA infrared images (Figure 8) reveal the three unique habitats of the Brant Splay and White's Splay.The forest community indicated in red is dominated by Sal.nigra, and the white to blueish white colors correspond to bare winter mudflats that transform during the growing season into the herbaceous marsh lands harvested in the study and the generally aquatic water areas mostly dissected interiorly by channel feeder passes to open-water shallow ponds.The areas of the two terrestrial habitats minus any open-water areas among them were used to calculate the splay growth rates of the areas that were at times supportive of submergent vegetation [46].For 22 years, the two Brant Pass terrestrial habitats grew to 5.87 km 2 , with forested and herbaceous lands occupying 1.23 km 2 and 4.64 km 2 , respectively.Both land types grew more quickly over the first measured interval (1978)(1979)(1980)(1981)(1982)(1983), then slowed (1983-1995), followed by an increase in growth (1995)(1996)(1997)(1998)(1999)(2000).The herbaceous wetlands grew four times faster than the forested areas during each of the three intervals.

River Quality and Wetland Impacts
Quantile regression analyses showed that both the water quality and quantity of the MR influence the prevalence and abundance of the two most ubiquitous herbaceous species (Sc.deltarum and Sag.latifolia), with significant impacts on the total AGB production of the splay lobe wetlands (Table 4).The summer river temperature had no effect on AGB, although the growing season length had a negative effect on both the total and Sc.deltarum AGB (Table 4).Discharge from the MR was identified as the only environmental factor that significantly affected both TAGB and the AGB of each species.The best predictor of AGB was spring discharge.Although both the linear and quadratic relationships were significant based on the likelihood ratio, the parameters for the relationships were only significant for the linear relationship during the spring season.Both the total and Sc.Deltarum AGB generally declined with increasing spring discharge, whereas the Sag.latifolia AGB increased with increasing spring discharge.The total AGB and that of Sc. deltarum were also significantly affected by nutrient loading, with the peak AGB occurring in the middle of the observed range.

Discussion
In this detailed long-term study, we followed both abiotic and biotic landscape changes within a delta lobe cycle from open water to open water [51] shaped by the MR.Because such cycles occur on multiple scales, our documentation of these changes can serve as a model for the ecology and phytosociology of Louisiana's Delta Plain.The four splays studied in the Balize Delta were initially colonized by PIO species, which autogenically modified the area (see [36]).After their colonization, the three herbaceous FS grew to dominate substrate areas according to elevation.Only at WSc, the most elevated site, did FS Sal.nigra persist past the early germinated seedling stage.During its study years, WSc supported this woody species past the seedling stage, although along the shoreline, its numerous seedlings were grazed upon, likely by deer, nutria, and/or wild boar.Away from the shoreline, growth into young trees confirms that WSc is the most elevated site overall.
Schoenoplectus deltarum became dominant on the two most elevated lobes over most years of each lobe's multidecadal existence.Consistently, on the lowest two lobes, Sag.latifolia became dominant or subdominant over many years, despite annually alternating with SD on higher-elevation areas and with SP on the lower lobe areas.Such switches in dominance depended on the inherent accretion/subsidence balance of the yearly substrate growing season along the lobe's subtle elevation gradient.The last step before submergence was only observed at the lowest mean elevation site, i.e., WSn, where near pure stands of Sag.platyphylla appeared in the last years before the lobe disappeared.The highest species richness occurred during the earliest years of splays.Both the richness and community diversity decreased as the FS became dominant.At the WSc and WSn sites, a general decline in AGB was observed over time, which may have been related to the manmade crevasse being less efficient in delivering sediments and nutrients to the WSc site.At the WSn site, the decline in AGB was related to the silting-in of the channel that delivered sediments and nutrients.White and Visser concluded that the total AGB on Brant's Splay was negatively affected by river discharge [49].When considering the two most common FS, the AGB of Sc. deltarum declined with increased discharge, whereas that of Sag.latifolia increased.
Deltaic rivers create abiotic physical landscape changes through seasonal sediment delivery, sediment allocation patterns, and annual flow, all of which alter the physical environment, with considerable ecological impacts [34,52].Deltaic vegetation responds to the river's seasonal temperature and nutrient loads and even the biota itself [49,[53][54][55].Vegetation affects flow patterns and sediment deposition [55], and, in turn, elevation change affects the vegetation composition [34].In addition, abiotic ecological processes of surface sedimentation and deep landform subsidence [56] can be allogenic agents of natural biotic change in all active deltaic areas.Any driver that impacts soil quality can force vegetation change, especially if the vegetation ultimately significantly impacts elevation [34,57].The phytosociology of each active cycle is not a perfect taxonomic timeline match due to interconnected processes governed by season, climate, geomorphology, and ecological events, each with timeline differences.Every forced taxonomic change is assuredly set by these inherent natural acts governed by the totality of the hydrological conditions that inform the coastal Louisiana delta and other deltas worldwide.For example, Sal. nigra stands grew into shrub-sized, plot-dominant individuals (Figure 3) over several years after an exceptionally high-water late spring in 1995, which brought new sediment that elevated the distal lobe substrates.
The concordant abundance and presence of vegetation changes on each lobe by at least the larger groups (FS and subdominants) confirm that considerable parallel elevation changes occurred on each lobe over decades (Figure 8b and [20]).However, portions of a lobe that remained monospecific over years do not fully preclude some change in elevation as long as they are within the species' tolerance (ecological niche).Interlobe vegetation dynamism, especially solely as the result of environmental change, is a delta wetland norm.Such common inter-site patterns over both short periods up to years and even longer periods are forced by sedimentation coupled with subsidence, as well as the abiotic impacts of tropical storms on elevation [58,59], which result in similar processes (sediment deposition, erosion, and compaction) common to coastal wetlands [57,60].All of the most common taxa on mud substrates of the Atchafalaya Delta [32] were harvested in this study.However, the most common FS, Sc. deltarum, has never been found in the Atchafalaya Delta of the Wax Lake Delta [32,57,61,62].
The identification of multiple significant lobe changes in this study brought on after flood pulse events could have only occurred over a long study period.Regional pulse events of different scales (i.e., large events [63,64] and many local events) are a common occurrence in deltaic lands [25,63].Carle et al. [65] showed the dramatic change in vegetation cover of the Wax Lake Delta following the record flood of 2011 that resulted in 35% of the delta shifting to higher-elevation communities.The largest event of this sort in our study occurred in the mid-1990s, which birthed the WSn site, enabling its detailed 18-year study.Another dramatic vegetation change occurred at the remaining three sites: a state change (sensu [66,67]) to highly statured populations of either cane (Ph.australis) or a woody species.Ph. australis often appeared in and around the study lobes in patches, representing the most common species both in terms of biomass and coverage metrics, growing over hundreds of km 2 in the outer Balize Delta, particularly its eastern lobes, for many decades [68].
In the Balize Delta, the highest species richness occurred during the pioneer stages.Both species richness and community diversity decreased as the FS became dominant, generally matching the Atchafalaya Delta community trends [57] related to species numbers and distributions over a 6-year study period.A significant net upward elevation change at the sites with different taxa was also recorded, especially the dominant groups affected by extensive grazing.Our data indicate that lower lobe communities and elevations were far less changed after community establishment, even over nearly three decades of study, partially due to the much higher amount of deep subsidence in the Balize Delta compared to that in the Atchafalaya Delta [69].Our most commonly harvested species, Sc. deltarum, the superior competitor over a range of elevations, was revealed through its monospecific stands, which likely competitively decreased the richness on the large Balize Delta lobe areas.
The young Atchafalaya Delta colonizing wetlands were studied by Montz [70] using coverage classes (sensu [71]).His permanent quadrat survey collected 241 species, identifying a rich community during colonization.Montz sampled all elevations and labeled two habitats: marshes and ridges, with the ridges about twice as rich [70].Bevington et al. [62] found a lower level of species diversity compared to 1980s levels reported in the same area and blamed the lack of species diversity on the dominance ("coverage") of the introduced Colocasia esculenta, a species that is far less prevalent in the Balize Delta than in the Wax Lake Delta [72] and that occurs in areas adjacent to the Sal.nigra community's lower-down region of mucky substrates without flow.
In this study, the subtle causes of diversity change were impacted by the breadth of the FS monospecific dominance according to site (Figure 9).The plot dominance of FS at WSn after the first study year (1994, a PSY) was not found at the other three more elevated sites after their PSYs.Even later, in their first decade of study, taxa from OTH and PIO species were still present (Figures 3 and 4) and far more often found only near the shoreline, the most elevated and naturally disturbed mud areas of the lobes.They can be labeled as "subordinately opportunistic"; rarely abundant; and, at times, even less common than PIO species.The most elevated site throughout the later study years, WSc, supported far greater quantities of PIO and OTH species than the two medium-elevation sites.Collected examples of these "opportunistic" OTH taxa had fleeting appearances of more than a few species that were occasionally unexplained when large single-season eruptions of a few particular taxa, such as V. luteola, Po. punctatum, and several grasses, especially Panicum spp.(Figure 3), occurred, although all three only accounted for 4.5% of total harvested AGB.Persistent seed rain or propagule drift in the waters seem to be the most likely causes when a colonizable niche area opens for often short-term appearances of such opportunists during the lobe cycle.After the pioneer years, at times only for a year, the opportunists covered only small ephemeral bare substrate patches created by at least four common stochastic events that deltas can experience: 1.
Sediment pulses, which are instances of spotty deposition or erosion from small sediment pulse events, which create local patches of ±elevation changes, forcing vegetation changes [80].They fit into the category of regeneration niches through the re-establishment process after a disturbance [81] and contributed considerably to the increase in the lobe-wide alpha diversity reported in this study, furthering lobe plant community dynamism;

4.
Vegetation sediment impact: Some areas were fashioned in more subtle ways than others in this study at the landscape scale by autogenic processes, that is, within the lobe ecosystem itself by the taxa forcing changes in sedimentation [82][83][84].The clearest example of this did not appear in the sampled plots but in localized lobe areas where particular vegetation (especially tall, introduced haplotypes, mostly of Typha sp., but at times Ph.australis; D.W., pers.obs.) created pockets of substrate elevation differences, often resulting in areas available for secondary establishment and further impacting the wetland's alpha diversity.Models have shown that over time, overall vegetation composition (vegetation quality) enhances variability in deltaic mobility (the coevolution of vegetation and delta channel networks), that is, "feedback" [85], which adds to the elevation diversity.Therefore, an increased overall beta diversity corresponds to cumulative vegetation patchiness [82,86], which is a well-understood relationship [32,65,87].The effects of both sediment pulsing and the vegetation composition are fully integrated into deltaic ecology [55].
Scaling up from the increased local channeling effects, community-wide diversity is generally greater in larger deltas.Larger splays naturally support greater areas of the highest elevation, i.e., the areas of forest and mixed marsh with the highest beta diversity reported in this study.Channel diversions have been shown to impact species richness [88].Generally, with significantly increased scaling in ecosystems, an uptick in the immigration of additional species logically follows [21].White and Visser [49] analyzed data from the BSs and BSn sites and found that AGB increased over time, which seemed mostly to be a function of the reported 0.9 °C/decade increase in river water temperature.However, their study only included the years when FS dominated.By adding years in this study, the trend was obscured by the inclusion of the high AGB observed in the second year after colonization, when PIO species dominated and benefitted from local seed input.Carle and Sasser [89] showed that in the Wax Lake Delta, the vegetated area increased over time, as well as the average biomass of the vegetated areas (estimated from NDVI, i.e., satellite imagery of vegetation vigor).They also revealed that the impacts of hurricanes on NDVI were relatively short-lived [89].
At the WSn site, the reduction in AGB was related to the silting-in of the channel that delivered sediments and nutrients, which resulted in a decline in the more elevated substrates with a greater diversity.Carle and Sasser [90] showed that NDVI was positively correlated with the water level along the Louisiana coast.They speculated that the effect Models have shown that over time, overall vegetation composition (vegetation quality) enhances variability in deltaic mobility (the coevolution of vegetation and delta channel networks), that is, "feedback" [85], which adds to the elevation diversity.Therefore, an increased overall beta diversity corresponds to cumulative vegetation patchiness [82,86], which is a well-understood relationship [32,65,87].The effects of both sediment pulsing and the vegetation composition are fully integrated into deltaic ecology [55].
Scaling up from the increased local channeling effects, community-wide diversity is generally greater in larger deltas.Larger splays naturally support greater areas of the highest elevation, i.e., the areas of forest and mixed marsh with the highest beta diversity reported in this study.Channel diversions have been shown to impact species richness [88].Generally, with significantly increased scaling in ecosystems, an uptick in the immigration of additional species logically follows [21].White and Visser [49] analyzed data from the BSs and BSn sites and found that AGB increased over time, which seemed mostly to be a function of the reported 0.9 • C/decade increase in river water temperature.However, their study only included the years when FS dominated.By adding years in this study, the trend was obscured by the inclusion of the high AGB observed in the second year after colonization, when PIO species dominated and benefitted from local seed input.Carle and Sasser [89] showed that in the Wax Lake Delta, the vegetated area increased over time, as well as the average biomass of the vegetated areas (estimated from NDVI, i.e., satellite imagery of vegetation vigor).They also revealed that the impacts of hurricanes on NDVI were relatively short-lived [89].
At the WSn site, the reduction in AGB was related to the silting-in of the channel that delivered sediments and nutrients, which resulted in a decline in the more elevated substrates with a greater diversity.Carle and Sasser [90] showed that NDVI was positively correlated with the water level along the Louisiana coast.They speculated that the effect was due to tidal flushing and increased soil oxygenation.In this study, no significant effect of flooding on AGB was found, which may be due to the difference in species composition in the two deltas.Biomass in the Balize Delta is dominated by the contributions of Sc. deltarum, whereas in the Wax Lake Delta, biomass was dominated by Sag.latifolia in the 1980s [61] and by Nelumbo lutea Willdenow in the 2000s [50].Nelumbo lutea is a species usually found on inundated to mucky substrates and seen only occasionally, apart from in our Balize Delta study areas in moderately sized and short-lived stands.
Carle and Sasser [89] showed no significant effect of river discharge on NDVI in the Wax Lake Delta.River discharge is a proxy for the amount of sediment that is deposited on the splays, which can bury tubers from the previous season, possibly delaying the emergence of Sc. deltarum and thereby reducing its biomass, giving a competitive advantage to Sag. latifolia.Such a phenomenon may be most important at the highest-elevation sites near the channels where Sc. deltarum dominates and where sediment deposition is the highest.AGB does not necessarily reflect the total annual productivity [90,91], especially for taxa that have a high level of biomass turnover (e.g., Sag.latifolia), which is common in freshwater marsh species [92].
The histories of all deltaic areas have a scale of growth (Figure 8; for this study, >2 decades = ~0.25 km 2 /yr), in which each unique land area for particular wetland vegetation is maintained over periods of time, even within short geological time frames if provided with sediment.Such growth was understandably much slower than for the very young Wax Lake Delta (1-3 km 2 /yr) [93], although far closer to its much higher earlier expected growth rate [31] reported in [93].Cahoon [59] predicted that a several-kilometerslong SBP lobe would require ~1.4 cm/yr of sediment uniformly deposited over the entire lobe after vegetated marsh establishment for it to remain static, which is an improbable scenario both in terms of supply and distribution.Such a scenario would equate to complete FS dominance to balance subsidence from compaction and deeper geological processes.Therefore, implications related to the quality and quantity of subsidence rates over time and tropical storm impacts do increase deltaic dynamism [59], resulting in vegetation changes.
With the additional impacts of potential salinity intrusion for deltas near ocean coastal areas, their dynamism becomes more complex [75].Soil organic and mineral amendments of the original deposited sediment appear to be the norm below ground, increasing over time as a result of production by the surface plant life, especially during the PSYs, as shown in previous studies of these wetlands [20,36], although likely many years in the future [84,[94][95][96].Such "biogenic accretion" [97], i.e., autogenic impact, contributes to significant elevation and vegetation changes.The general "Fluvial Biogeomorphic Succession Theory" [98] uniquely details the feedbacks that forced the considerable changes observed in our study.The local-to-regional hydrological patterns within the Balize Delta create and support wetland plant composition over decades, as well as the diversity of the larger community.The plants themselves create change based on changes in water as well as where and when they occur [99].Elevation dynamics are a key driving factor for vegetation dynamics within deltaic regions created mostly by the region's hydrology [32,62] coupled with its sediment loads in water [1].When sedimentation is restricted by changes in channel morphology or low river flow and falls out of balance with subsidence (subsidence may dominate), the niche space of the existing vegetation communities/species naturally changes.Finally, when channels silt in or change course, a site, a lobe, or a splay revert to open water.
Niche space theory suggests that the diversity (taxa) and quantity (AGB) of harvested vegetation in a studied wetland convey important information about the deltaic environmental conditions through time.To discover the reasons for the large area dominated by only a few species (FS), even in an ecosystem with considerable species diversity over a delta cycle, a large-scale study such as that reported herein is required, confirming the geological complexity of deltas [100].The accommodation space for the taxa in all deltaic environments controlled by sedimentation is well established, with a clear impact of climate on sedimentation [101].Large similar wetland areas have been noted as markers of global climate change [102].In this study, data on the number of individuals of each species were not collected in our plots, making it impossible to calculate standard dominance indices (e.g., SIMPSON [103] or Berger-Parker [104]).We recommend that future long-term studies include such a measurement.Measurements of stem density, diameter, and height, which are important parameters for hydrologists that model flow through and over vegetated surfaces, are also recommended.
The stochastic nature of colonization and ultimate establishment by deltaic plants impact deltaic patterns, leading to in situ channel formations.Frans et al. [105] showed that in addition to the impacts of climate change on sediment, runoff levels from the replacement of natural cover species to agricultural crops in the upper Mississippi's drainage area have increased nutrient contents (particularly nitrogen) within the river water.Therefore, impacts on vegetation quantity and diversity are likely to follow.
With major global increases in sediment flux [106], river flow (discharge, not fully considering all other river attributes), and sea levels [107] as a result of climate change [12,108], the future of the Louisiana coastline, with its usual dynamism, will become more important to parse.The inherently short-term, subtle impacts of climate necessitate long-term studies [109].Scales of change over the decadal time periods reported in this study are illuminating.Hughes et al. [110] put it correctly when they stated: "long-term studies do indeed contribute disproportionately to ecology and environmental science".Deltaic ecosystems suffer from both natural and human perturbations.Both hinder a full ecological understanding of impactful stochastic events.Past human impacts on freshwater systems including river-dominated deltas indeed have memorable effects [111].
LA. Finally, Tulane University's Stephen Darwin must be acknowledged for his devotion to plant taxonomy; he made this detailed taxonomic study possible.
Water 2023, 15, x FOR PEER REVIEW

Figure 1 .
Figure 1.Louisiana's southern coastline: a landscape of nested deltaic patterns.Panel I:The entire Louisiana coast, the coastal plain, and its 3 geomorphological west-to-east divisions-the Chenier Plain[28]; the Delta Plain[13]; and, lastly, the easternmost Pearl Plain (advocated for in[29]).The largest is the Delta Plain, made up of seven geologically distinct deltaic regions formed sequentially (youngest to oldest: "A" through "G", with "B" essentially an extension of delta "C" (Plaquemines Delta)) over the past ~7500 yrs[30].The newest delta region (A) in the Delta Plain is the region around the terminus of the Atchafalaya River, which is composed of both the Atchafalaya and Wax Lake Deltas[31,32].The Balize Delta (B), the focus of this study, is shown in Panel II (modified from[33]).The Balize Delta began its growth ~550 yrs ago and is composed of four distinct, smaller-scale subdeltas identified sequentially by age ("a, b, c, and d").Panel III: Embedded within the four subdeltas of the extant Balize Delta are many smaller "inner deltaic splays", filling in shallow ponds, indicative of a senescing Balize Delta.Panel IV: Inner deltaic splays are "mini deltas", each composed of multiple lobes; each lobe likewise has its own birth, growth, and decline phases.The study splay and location of the study sites are marked with stars.Panel V: An individual lobe; the BSn site during winter in a later year of growth from Panel IV.Note 1: Panel IV clearly shows the shoreline crevasse through Octave Pass, which created the eventual WSc splay.Panel IV also labels four nascent "natural" splays: two active (AS) and two senescing (SS).Note 2: The NASA images of Panels IV and V clearly indicate an attenuated gradual slope down and away from the lobe shoreline feeder pass, as represented by the near-white color of the most elevated exposed mud, with the shallowest mud represented by the darker blue color of the receiving pond.Note 3: The border colors of Panels II through V correspond to the increasingly larger ovoid outlines around the green dot representing the study location.
Figure 1.Louisiana's southern coastline: a landscape of nested deltaic patterns.Panel I:The entire Louisiana coast, the coastal plain, and its 3 geomorphological west-to-east divisions-the Chenier Plain[28]; the Delta Plain[13]; and, lastly, the easternmost Pearl Plain (advocated for in[29]).The largest is the Delta Plain, made up of seven geologically distinct deltaic regions formed sequentially (youngest to oldest: "A" through "G", with "B" essentially an extension of delta "C" (Plaquemines Delta)) over the past ~7500 yrs[30].The newest delta region (A) in the Delta Plain is the region around the terminus of the Atchafalaya River, which is composed of both the Atchafalaya and Wax Lake Deltas[31,32].The Balize Delta (B), the focus of this study, is shown in Panel II (modified from[33]).The Balize Delta began its growth ~550 yrs ago and is composed of four distinct, smaller-scale subdeltas identified sequentially by age ("a, b, c, and d").Panel III: Embedded within the four subdeltas of the extant Balize Delta are many smaller "inner deltaic splays", filling in shallow ponds, indicative of a senescing Balize Delta.Panel IV: Inner deltaic splays are "mini deltas", each composed of multiple lobes; each lobe likewise has its own birth, growth, and decline phases.The study splay and location of the study sites are marked with stars.Panel V: An individual lobe; the BSn site during winter in a later year of growth from Panel IV.Note 1: Panel IV clearly shows the shoreline crevasse through Octave Pass, which created the eventual WSc splay.Panel IV also labels four nascent "natural" splays: two active (AS) and two senescing (SS).Note 2: The NASA images of Panels IV and V clearly indicate an attenuated gradual slope down and away from the lobe shoreline feeder pass, as represented by the near-white color of the most elevated exposed mud, with the shallowest mud represented by the darker blue color of the receiving pond.Note 3: The border colors of Panels II through V correspond to the increasingly larger ovoid outlines around the green dot representing the study location.

Figure 1 .
Figure 1.Louisiana's southern coastline: a landscape of nested deltaic patterns.Panel I:The entire Louisiana coast, the coastal plain, and its 3 geomorphological west-to-east divisions-the Chenier Plain[28]; the Delta Plain[13]; and, lastly, the easternmost Pearl Plain (advocated for in[29]).The largest is the Delta Plain, made up of seven geologically distinct deltaic regions formed sequentially (youngest to oldest: "A" through "G", with "B" essentially an extension of delta "C" (Plaquemines Delta)) over the past ~7500 yrs[30].The newest delta region (A) in the Delta Plain is the region around the terminus of the Atchafalaya River, which is composed of both the Atchafalaya and Wax Lake Deltas[31,32].The Balize Delta (B), the focus of this study, is shown in Panel II (modified from[33]).The Balize Delta began its growth ~550 yrs ago and is composed of four distinct, smallerscale subdeltas identified sequentially by age ("a, b, c, and d").Panel III: Embedded within the four subdeltas of the extant Balize Delta are many smaller "inner deltaic splays", filling in shallow ponds, indicative of a senescing Balize Delta.Panel IV: Inner deltaic splays are "mini deltas", each composed of multiple lobes; each lobe likewise has its own birth, growth, and decline phases.The study splay and location of the study sites are marked with stars.Panel V: An individual lobe; the BSn site during winter in a later year of growth from Panel IV.Note 1: Panel IV clearly shows the shoreline crevasse through Octave Pass, which created the eventual WSc splay.Panel IV also labels four nascent "natural" splays: two active (AS) and two senescing (SS).Note 2: The NASA images of Panels IV and V clearly indicate an attenuated gradual slope down and away from the lobe shoreline feeder pass, as represented by the near-white color of the most elevated exposed mud, with the shallowest mud represented by the darker blue color of the receiving pond.Note 3: The border colors of Panels II through V correspond to the increasingly larger ovoid outlines around the green dot representing the study location.
Figure 1.Louisiana's southern coastline: a landscape of nested deltaic patterns.Panel I:The entire Louisiana coast, the coastal plain, and its 3 geomorphological west-to-east divisions-the Chenier Plain[28]; the Delta Plain[13]; and, lastly, the easternmost Pearl Plain (advocated for in[29]).The largest is the Delta Plain, made up of seven geologically distinct deltaic regions formed sequentially (youngest to oldest: "A" through "G", with "B" essentially an extension of delta "C" (Plaquemines Delta)) over the past ~7500 yrs[30].The newest delta region (A) in the Delta Plain is the region around the terminus of the Atchafalaya River, which is composed of both the Atchafalaya and Wax Lake Deltas[31,32].The Balize Delta (B), the focus of this study, is shown in Panel II (modified from[33]).The Balize Delta began its growth ~550 yrs ago and is composed of four distinct, smallerscale subdeltas identified sequentially by age ("a, b, c, and d").Panel III: Embedded within the four subdeltas of the extant Balize Delta are many smaller "inner deltaic splays", filling in shallow ponds, indicative of a senescing Balize Delta.Panel IV: Inner deltaic splays are "mini deltas", each composed of multiple lobes; each lobe likewise has its own birth, growth, and decline phases.The study splay and location of the study sites are marked with stars.Panel V: An individual lobe; the BSn site during winter in a later year of growth from Panel IV.Note 1: Panel IV clearly shows the shoreline crevasse through Octave Pass, which created the eventual WSc splay.Panel IV also labels four nascent "natural" splays: two active (AS) and two senescing (SS).Note 2: The NASA images of Panels IV and V clearly indicate an attenuated gradual slope down and away from the lobe shoreline feeder pass, as represented by the near-white color of the most elevated exposed mud, with the shallowest mud represented by the darker blue color of the receiving pond.Note 3: The border colors of Panels II through V correspond to the increasingly larger ovoid outlines around the green dot representing the study location.
Voucher materials at LSU herbarium, Baton Rouge, LA, USA.Taxa are listed in order of total aboveground biomass harvested during the study (TAGB).Abbreviations: OTH = other species, PTB = percentage of total aboveground biomass harvested over the 28-year study provided >1% of 116,709 g total.* Species defined as common pioneers; such species represent 82.4% of the pioneer species' TAGB.** Species only identified at the family level represent 1.7% of TAGB.*** Species only identified at the genus level were sterile and found on only one transect, representing 0.2% of TAGB.

Figure 2 .
Figure 2. Illustration of a typical splay lobe (modified from[20]  thanks to J. Lynch) with elevation change and spatial relationships depicting the spatial and height relationships of the wetland habitat on its platform range of <1 m from the shallow pond bottom (depth of ~ 50 cm) at the lobe's low end of the emergent mudflat (+17 cm) to the low marsh (+26 cm), the high marsh (+30 cm), and the most upland Salix nigra (SN) forest (>6 cm; mean of three laser level readings).Thus, the flood tolerance for the four predominant community-identifying species (SN > SD > SL > SP) extends through only ~62 cm elevation, regardless of the lobe's length.A lobe's length can range from 50 meters to kilometers, as in the case of the BSs lobe.Key: SN = Salix nigra; SD = Schoenoplectus deltarum-dominated marsh; SL = Sagittaria latifolia low marsh; SP = Sagittaria platyphylla lowest marsh; and SAV = submerged aquatic vegetation.

Figure 2 .
Figure 2. Illustration of a typical splay lobe (modified from[20]  thanks to J. Lynch) with elevation change and spatial relationships depicting the spatial and height relationships of the wetland habitat on its platform range of <1 m from the shallow pond bottom (depth of ~50 cm) at the lobe's low end of the emergent mudflat (+17 cm) to the low marsh (+26 cm), the high marsh (+30 cm), and the most upland Salix nigra (SN) forest (>6 cm; mean of three laser level readings).Thus, the flood tolerance for the four predominant community-identifying species (SN > SD > SL > SP) extends through only ~62 cm elevation, regardless of the lobe's length.A lobe's length can range from 50 meters to kilometers, as in the case of the BSs lobe.Key: SN = Salix nigra; SD = Schoenoplectus deltarum-dominated marsh; SL = Sagittaria latifolia low marsh; SP = Sagittaria platyphylla lowest marsh; and SAV = submerged aquatic vegetation.
Figure 2. Illustration of a typical splay lobe (modified from[20]  thanks to J. Lynch) with elevation change and spatial relationships depicting the spatial and height relationships of the wetland habitat on its platform range of <1 m from the shallow pond bottom (depth of ~50 cm) at the lobe's low end of the emergent mudflat (+17 cm) to the low marsh (+26 cm), the high marsh (+30 cm), and the most upland Salix nigra (SN) forest (>6 cm; mean of three laser level readings).Thus, the flood tolerance for the four predominant community-identifying species (SN > SD > SL > SP) extends through only ~62 cm elevation, regardless of the lobe's length.A lobe's length can range from 50 meters to kilometers, as in the case of the BSs lobe.Key: SN = Salix nigra; SD = Schoenoplectus deltarum-dominated marsh; SL = Sagittaria latifolia low marsh; SP = Sagittaria platyphylla lowest marsh; and SAV = submerged aquatic vegetation.

Figure 3 .
Figure 3. Histograms of the major dominant vegetation (eight species plus three Msc.groups thereof) at the four lobes in terms of percentage of plots dominated (greatest AGB) per study year.The AGB values of the five least harvested PIO species of the nine total (only 2.4% of total harvested) were added to their corresponding Msc groups to simplify the trends.Sites are labeled a to d (high to low) according to relative composite elevation.

Figure 3 .
Figure 3. Histograms of the major dominant vegetation (eight species plus three Msc.groups thereof) at the four lobes in terms of percentage of plots dominated (greatest AGB) per study year.The AGB values of the five least harvested PIO species of the nine total (only 2.4% of total harvested) were added to their corresponding Msc groups to simplify the trends.Sites are labeled a to d (high to low) according to relative composite elevation.

Figure 5 .Figure 5 .
Figure 5. Nonmetric multidimensional scaling (nMDS) model of the aggregate trajectories for changes in plant community composition at the four lobe sites under study (WSc = White's Splay constructed; BSs = Brant Splay south; BSn = Brant Splay north; and WSn = White's Splay new) based on the harvested aboveground biomass (AGB) data of 13 species: 4 foundation species and 9 pioneers.A DistLM marginal test identified the 11 top species jointly on the four lobes of significance based on their high AGB among the 53 total taxa (p > 0.05; three herbaceous FS, six PIO, and two OTH; Table2) using the ordination's proportion rankings.The rankings define the structural change in the lobe vegetation during the decades of the study of the delta cycle.If only the face value of the AGB metric (abundance) was used in the analysis, taxa with high study-wide harvest values would with statistically significant importance

Figure 6 .
Figure 6.Number of species and species groups harvested in 965 plots at four splay lobe sites over the 28-year study.The solid mean trend line shows the statistically significant decline in richness after the first late summer of mudflat plant colonization (R 2 = 0.781; p-value = 9.82 × 10 -10 ).Abbreviations: WSc = White's Splay constructed; BSs = Brant Splay south; BSn = Brant Splay north; WSn = White's Splay new.Note 1: The significant drop in richness at BSn during 1989/1990 occurred because of a serious salt burn event just days before harvest.Note 2:The 1994 richness at WSn was high (16 species) compared to the other sites during the initial year of harvest because WSn's first year of harvest was year 3 of the wetland's existence (see text), i.e., the year understood to be the last of the primary succession period[36].

Figure 6 .
Figure 6.Number of species and species groups harvested in 965 plots at four splay lobe sites over the 28-year study.The solid mean trend line shows the statistically significant decline in richness after the first late summer of mudflat plant colonization (R 2 = 0.781; p-value = 9.82 × 10 −10 ).Abbreviations: WSc = White's Splay constructed; BSs = Brant Splay south; BSn = Brant Splay north; WSn = White's Splay new.Note 1: The significant drop in richness at BSn during 1989/1990 occurred because of a serious salt burn event just days before harvest.Note 2:The 1994 richness at WSn was high (16 species) compared to the other sites during the initial year of harvest because WSn's first year of harvest was year 3 of the wetland's existence (see text), i.e., the year understood to be the last of the primary succession period[36].

Water 2023 , 31 Figure 8 .
Figure 8.(a) Winter set of NASA infrared images of the two studied deltaic splays that show the period from birth after two natural levee breaches to maximum growth (modified from [46] in [20]).The splays began growing in a ~10 km 2 shallow pond (an average of < 1 m deep) east of the Mississippi River (see Figures 1 and 2 for location context).The color shades in the 1995 image abruptly changing from reddish black to blue/white identify the two clearly dominant splay wetland types: forests up-slope and the more distal winter surface bare mudflats of summer marshes, with the shallow aquatic open water being the most distal.(b) Five-year elevation change data of five sites along the South Brant Pass lobe (see 2000 image) [20].Relative elevations to pond bottom by site are

Figure 8 .
Figure 8.(a) Winter set of NASA infrared images of the two studied deltaic splays that show the period from birth after two natural levee breaches to maximum growth (modified from [46] in [20]).

Water 2023 , 31 Figure 9 .
Figure 9. Photos of the typical dominant wetland communities of a splay lobe within the Balize Delta, Louisiana, USA.All photos were taken by D.W. Note: the mixed highest marsh image depicts a large coverage that is generally a less extensive transition zone between forest and other marshes.* See [20].

Figure 9 .
Figure 9. Photos of the typical dominant wetland communities of a splay lobe within the Balize Delta, Louisiana, USA.All photos were taken by D.W. Note: the mixed highest marsh image depicts a large coverage that is generally a less extensive transition zone between forest and other marshes.* See [20].

Table 1 .
Species groups and total aboveground biomass for 53 vascular plant species and groups thereof harvested (1984-2012) on deltaic lobes within the Balize Delta of the Mississippi.

Table 2 .
Distance-based linear model (DistLM) marginal test results, relative importance components, and their rankings.
Notes:The table includes 11 species (of the 53 total harvested taxa) with statistical significance of p < 0.05.Two additional species had a study-wide AGB greater than 2 among the 11 listed species (Salix nigra, 2.3% (6th) and Typha domingensis, 0.90% (11th; see Table1), but neither showed statistical significance in the DistLM.Two multiple-species taxonomic groups (Panicum spp., 1.1%; and Poaceae, 1.1%) had a slightly higher AGB than 3 of the 11 taxa; the total AGB of these four taxa made up just 5.4% study-wide, with Salix nigra (2.3%) accounting for the largest proportion.Abbreviations: P = statistical probability (all > 0.05); R = ranking according to the values of the previous column; FS = foundation species, PIO = pioneer species; OTH = other species.a Species identified as significant by DistLM marginal test.b Proportion of the community composition explained by the DistLM marginal test.c Abundance defined as the percentage of the total aboveground biomass harvested study-wide (see Table

Table 3 .
Correlation among sites in total aboveground biomass during their common years of sampling after primary succession years.
Notes: Abbreviations: WSc = White's Splay created; BSs = Brant Splay south; BSn = Brant Splay north; WSn = White's Splay new; and PPSY = post primary succession year.a Significant values (p < 0.05) are indicated in bold.b Primary succession years are included.

Table 4 .
Quantile regression analysis results for the aboveground biomass (AGB) response to environmental change.Values are presented as the p-value of the likelihood ratio of the fit.Significant values (p < 0.05) are indicated in bold.Only parameters with at least one significant effect are included. Notes: