Effects of Small-Scale Armoring and Residential Development on the Salt Marsh-Upland Ecotone

Small-scale armoring placed near the marsh-upland interface to protect single-family homes is widespread but understudied. Using a nested, spatially blocked sampling design on the coast of Georgia, USA, we compared the biota and environmental characteristics of 60 marshes adjacent to either a bulkhead, a residential backyard with no armoring, or an intact forest. We found that marshes adjacent to bulkheads were at lower tidal elevations and had features typical of lower elevation marsh habitats: high coverage of the marsh grass Spartina alterniflora, high density of crab burrows, and muddy sediments. Marshes adjacent to unarmored residential sites had higher soil water content and lower porewater salinities than the armored or forested sites, suggesting that there may be increased freshwater input to the marsh at these sites. Deposition of Spartina wrack on the marsh-upland ecotone was negatively related to elevation at armored sites and positively related at unarmored residential and forested sites. Armored and unarmored residential sites had reduced densities of the high marsh crab Armases cinereum, a species that moves readily across the ecotone at forested sites, using both upland and high marsh habitats. Distance from the upland to the nearest creek was longest at forested sites. The effects observed here were subtle, perhaps because of the small-scale, scattered nature of development. Continued installation of bulkheads in the southeast could lead to greater impacts such as those reported in more densely armored areas like the northeastern USA. Moreover, bulkheads provide a barrier to inland marsh migration in the face of sea level rise. Retaining some forest vegetation at the marsh-upland interface and discouraging armoring except in cases of demonstrated need could minimize these impacts.


Introduction
Humans have been living near, and protecting themselves from, the ocean for millennia (Doody 2004;Popkin 2015). Although it was historically assumed that coastal areas would accrete land on the seaward side of shoreline armoring or seawalls, creating more upland (Doody 2004), the modern understanding of the land-sea border suggests the opposite. Instead, hard structures steepen and shorten intertidal habitats, leading to a loss of area in a phenomenon known as coastal squeeze (Pethick 2001;Dugan et al. 2011). At present, 14% of the tidal shoreline within the continental US is armored, and armoring is expected to increase in the next century (Gittman et al. 2015).

Communicated by Carolyn A. Currin
Electronic supplementary material The online version of this article (doi:10.1007/s12237-017-0300-8) contains supplementary material, which is available to authorized users.
Armoring can be less effective at protecting the shoreline from erosion than natural defenses. For example, marshes provide greater erosion protection against the effects of a Category 1 storm than bulkheads (Gittman et al. 2014). In addition, the decrease in structural complexity associated with armoring often supports fewer species than found in natural shorelines (Chapman 2003;Gittman et al. 2015), and the introduction of novel substrate can facilitate species invasions (Landschoff et al. 2013). Many studies of coastal development and armoring have focused on the effects of extreme examples: development near dense human populations, or large seawalls in high-energy environments (Silliman and Bertness 2004;Long et al. 2011). Previous work also often examined armoring placed at low tidal elevation, or associated with new upland created by filling large areas (Long et al. 2011;Balouskus and Targett 2012;Lowe and Peterson 2014;Lowe and Peterson 2015). While these are important studies, they confound the effects of the armoring per se with the effects of intertidal habitat loss (but see Bozek and Burdick 2005).
Few studies have examined how the more commonplace, low-density residential development and more modest types of shoreline armoring affect coastal marsh habitats (Walters et al. 2010;Bozek and Burdick 2005). Small-scale forms of armoring, such as bulkheads, are commonly used to protect single-family homes that are adjacent to salt marshes in the southeastern US. Historically, these structures were installed to fill and reclaim land (Doody 2004). Although rules vary from state to state, filling of marshes is now generally prohibited. However, homeowners still often place bulkheads at the marsh-upland ecotone in order to guard against erosion, sea level rise, and flooding (Scyphers et al. 2014). These types of structures are typically about 1 m tall and located at or just above the high-tide line. The few studies of modest armoring in low-density residential developments have typically found few or very subtle effects rather than large effects, with effects concentrated in the high marsh (Bozek and Burdick 2005;Walters et al. 2010). These studies, however, may have suffered from limited replication (four to five pairs of developed and control sites) and only considered one type of development.
In this study, we aimed to evaluate the effect of placing a hard substrate at the upland-marsh ecotone by studying bulkheads that were placed above the high-tide mark adjacent to salt marshes. It is likely that these bulkheads have less of an effect than those placed lower in the tidal profile or in higher energy environments (Dugan et al. 2017). However, we hypothesize that they still alter the flow of fresh water and associated nutrients from the upland to the marsh, and impede the movement of animals in both directions. Severing sediment supply from the upland may result in lower elevations next to bulkheads, which could alter plant and invertebrate communities. For example, armoring can sequester sediments previously supplied by an eroding upland, leading to sediment starvation of environments seaward of the structure (Nordstrom et al. 2009;Nordstrom and Jackson 2013). It is possible, however, that not all the effects of armoring are negative: armoring may protect the upper marsh by limiting runoff from upland development.
Even if homeowners do not install bulkheads, residential development alone may have impacts on adjacent marsh ecosystems (McClelland et al. 1997;Bertness et al. 2002;Fitch et al. 2009). In the northeastern US, upland development has been linked to eutrophication and changes in the plant communities in the upper marsh (Bertness et al. 2002;Bozek and Burdick 2005;Fitch et al. 2009). However, development intensity is substantially lower in the southeast and the marshes are larger, so it is unclear whether results from the northeast can be extrapolated to the US east coast as a whole. The small amount of work that has been conducted in the southeastern US suggests that, in fact, these marshes show subtler and different responses to development than do those in New England (Walters et al. 2010). The southeastern US coast is projected to have the highest rate of human population growth from 2010 to 2020 in the coastal US (Crossett et al. 2005;Bamford 2013), so it is important both to understand impacts of current coastal development and predict the effects of more intense, future development in this area.
To separate the effects of low-intensity residential development and shoreline armoring, we compared salt marshes adjacent to upland that was (1) armored and developed (Barmoredŝ ites); (2) unarmored and developed (Bunarmored^sites); and (3) unarmored and forested (Bforested^sites). We hypothesized that upland modifications at either type of developed site would alter the extent and composition of the high marsh community, with marshes adjacent to bulkheads exhibiting the greatest effects because the upland was both developed and armored. Focusing on the upper salt marsh, we evaluated how site type affected the following: (1) physical and environmental characteristics; (2) biological characteristics; (3) the relationship between physical and biological characteristics; and (4) the use of terrestrial habitats by an organism that moves routinely between the upland and marsh.

Field Survey Design and Site Selection Methods
We surveyed high marsh characteristics at 20 blocked stations ( Fig. 1) along the Georgia coastline. We used GIS data on the locations of armored shorelines (Alexander 2010) and land use to select an armored, unarmored, and forested site at each station (20 stations × 3 site types = 60 sampling sites). We limited our armored sites to locations where the bulkhead was placed at the marsh-upland boundary, adjacent to a single-family home. None of the armored sites displayed obvious evidence of buildout (i.e., armoring that directly covers and replaces marsh) based on an evaluation of current and historical (1942, 1972, and 2009/2010) aerial photography. Bulkheads were between 0.7 and 1.77 m in height (average = 0.87 m), measured from the marsh surface on the seaward side to the top of the bulkhead. All 20 bulkheads were made of wood; seven also included vinyl siding as part of their structure. Armored and unarmored sites all had single-family homes with lawns adjacent to the marsh, with the main difference being the presence or absence of a bulkhead. Forested sites were selected at locations where no development existed adjacent to the marsh and forest vegetation was prevalent in the upland and extended to the marsh-upland boundary. Within each station, the defining land use characteristic of each of our three site types (i.e., bulkhead, lawn, or forest) had at least 20 m of frontage, and sites within a station were separated by at least 20 m. The marsh-upland boundary was used as the Bzero^line for each site, which was delineated either by the location of the bulkhead or the edge of the lawn or the forest. Two transects were run perpendicularly from the zero line into the upper marsh, with sampling points established at 2, 4, and 8 m from the marsh edge boundary (Fig. 1).
Question 1: How does site type affect the physical and environmental characteristics of the upper marsh?
We characterized upper marsh geomorphology (elevation) and stratigraphy, soil characteristics (grain size, organic matter content), and porewater (salinity, nutrients) at each site. We used a Trimble real time kinematic (RTK; model R6 and R8) GPS with a virtual reference network to measure the elevation and location (latitude and longitude) of each sampling point. Elevation was referenced to the North American Vertical Datum of 1988. We estimated the distance from the upland to the nearest creek (henceforth, Bupland-creek distance^) at each site using GIS, by measuring the shortest distance from the upland or structure to the first substantial creek.
Porewater was collected from surficial soils at each sampling point using Rhizon Core Solution Samplers with a 10cm hydrophilic porous polymer tube (Rhizosphere Research Products). Samples were frozen at − 80°C prior to conducting analyses for porewater nutrient content and salinity. Ammonium concentrations were determined using the phenol-hypochlorite method (Koroleff 1983) with a Shimadzu UV-1601 spectrophotometer. Nitrate + nitrite (reported as nitrate) and phosphate concentrations were measured on an Alpkem RFA-300 autoanalyzer. We used EPAapproved methods to analyze nitrate (4500-NO 3 − automated cadmium reduction method) and phosphate (4500-P automated ascorbic acid reduction method) (Rice et al. 2012). All nutrient samples were analyzed in triplicate. We measured salinity in collected porewater with a handheld refractometer (Vee Gee STX-3).
To evaluate sediment organic matter and water content, we collected a 10-cm sediment core from each sampling point. Samples were stored at ambient temperature and brought to the lab for processing. Sediment water content was measured by drying samples for 3 days at 60°C. Sediment organic content was determined by weight loss after combustion at 440°C overnight. We collected sediment samples from the marsh surface (0-2 cm) in each quadrat for grain size analysis. Sediment samples were wet-sieved using standard protocols through a 63-μm (4-phi) sieve (Alexander et al. 1986). The coarse fraction (> 63 μm) was then dried and sieved through stacked sieves starting at − 1 phi (2 mm) to separate gravel (larger than 2 mm) from sand (2 mm-63 μm) at 0.25-phi intervals. The percentage of mud (< 63 μm) was quantified by drying an aliquot of the total mud fraction captured during wet sieving. If sufficient quantities of mud existed (> 10% by weight), the silt and clay grain size distributions were determined with a Micromeritics Sedigraph 5100. If the sample . At each station we selected three sites varying in site type categorization (forested, unarmored, armored), with sampling points at three distances (black bars, 2, 4, and 8 m) from the marsh-upland boundary (starred), along two transects (one transect shown in illustration; middle panel). Illustrations at the right provide an example of each site type from within a single station of the survey contained < 10% mud, an additional aliquot was taken to quantify the percent silt and clay in the sample. Sediment statistics (e.g., mean size and sorting; standard deviation) were derived from these data using the method of moments (Griffiths 1967). At each sampling site, we collected a 30cm core 2 m from the marsh-upland boundary for stratigraphic analysis; we described cores using the grain size nomenclature of Folk (1980) to illustrate broad-scale, down-core trends in sediment character.
Statistical Analysis All statistical analyses were done using R 3.1.3 (R Development Core Team 2015). To evaluate the relationship of site type with physical and environmental characteristics, we fit linear mixed-effects models. First we evaluated the relationship between site type, elevation, and upland-creek distance with linear mixed-effects models, using site type as the fixed effect and station and sampling point along the transect nested in station as the random effect (packages lme4 and lmerTest; Bates et al. 2014;Kuznetsova et al. 2015). Because elevation and upland-creek distance are known to affect many other environmental variables in a salt marsh (Hladik and Alber 2014), we analyzed each environmental response with the linear mixed-effects model just described, but with the addition of elevation and upland-creek distance as fixed effects. All models were run initially with an interaction between elevation and site type, which was retained if the interaction term was significant or if including it qualitatively changed model results. To compare all site types, we conducted Tukey's post hoc analysis (package multcomp; Hothorn et al. 2008). Data were evaluated for model assumptions and some variables were log or logit transformed to satisfy assumptions.
Question 2: How does site type affect the biological community of the upper marsh?
We quantified flora (vegetation composition) and fauna (snail, bivalve, and crab abundance) in the upper marsh ecotone. We measured vegetative cover in two ways; first we took an overhead photograph of every sampling point to produce an estimate of total vegetative cover. Next, we measured percent cover of each plant species at each sampling point using a 0.5 × 0.5 m quadrat subdivided into 100 cells. Wrack and bare mud (no vegetation or wrack) were included as categories. We counted eastern melampus snails (Melampus bidentatus), and marsh periwinkles (Littoraria irrorata; Melampus and Littoraria, respectively, henceforth), and any other snails visible in a 0.25 × 0.25 m quadrat in each sampling point. If snails were rare, we used a larger quadrat size, up to 1 × 1 m. Crabs are highly mobile and affected by human presence, so we counted crab burrows (> 0.5 cm diameter) in a 0.25 × 0.25 m quadrat as a proxy for their density (Mouton and Felder 1995). We identified as many crabs as possible upon arriving at each sampling point. We counted the ribbed mussel Geukensia demisa in a 1 × 1 m quadrat.
Bulkheads can provide novel hard substrate, which can lead to new (and sometimes invasive) species recruiting to areas with armoring. At the bulkhead sites, we counted species living on the bulkhead by searching a 5-m length (usually between the two transects), and identifying benthic invertebrates such as barnacles, mobile crustaceans such as crabs, macroalgae, and any other species using the bulkheads as habitat.
Statistical Analysis To compare the biological communities among site types, we fit Bray-Curtis similarity matrices for abundance and percent cover data for all quadrats and visualized these matrices with Multidimensional Scaling (MDS). The advantage of MDS analysis is that it collapses the complex variability of species composition and abundance to major modes of variability among site types (Clarke and Warwick 2006;Siddon et al. 2011). To lower the influence of highly abundant species on the analyses, data were squareroot transformed prior to analysis. We applied multivariate ANOVA using the Bray-Curtis similarity matrices (package adonis) to test for the effect of site type on marsh biological communities, with site type as well as sampling point along the transect nested in station (R, package vegan; Oksanen et al. n.d.). We conducted a similarity of percentages (SIMPER) analysis to evaluate the contribution of each species to the Bray-Curtis similarity matrices to determine which species contributed most to differences among biological communities by site type.
Question 3: How does site type affect the relationship between the physical and biological characteristics in the upper marsh?
Statistical Analysis To evaluate the relationship between the environmental variables and biological community composition, we performed additional statistical analyses informed by the output from Question 2. MDS collapses the variance in the biological community to three axes, so we used a multivariate ANOVA (package car; Fox and Weisberg 2011) with the three axes of the MDS output as the response variables. Elevation, percent sand, soil water content, soil organic matter, porewater salinity, and porewater concentrations of ammonium, nitrate, and phosphate were used as environmental predictor variables.
To further evaluate the response of individual species, we used the outputs from the MDS and the multivariate ANOVA as a guide and correlated the distribution of the top species that contributed to the differences between site types and the top environmental variables that were associated with the biological community. We fit mixed-effects linear or generalized linear models to the top four species (or organism proxy in the case of crabs burrows) contributing to differences in the biological community based on output from the SIMPER analysis in Question 2. These were percent cover of Spartina alterniflora and Juncus romerianus (hereafter Spartina and Juncus, respectively), crab burrow counts, and Littoraria abundance. Because we used a gridded quadrat method, our percent cover data derives from cumulative presence/absence data, and so we used a linear model and logit transformed Spartina and Juncus percent cover. We used a Poisson model for the crab burrow counts and a linear model for log-transformed Littoraria abundance. We used environmental variables that were correlated with the biological community as predictor variables. Our study was conducted along the entire Georgia coast, and many of the variables included in our survey likely vary across this biogeographic range. To account for our sampling design and tease apart the relationships between the environmental variables and biological variables across the geographic range, we used station and sampling point along the transect (2, 4, or 8 m) as a random effect in each individual model (packages lme4 and lmerTest; Bates et al. 2014;Kuznetsova et al. 2015). To compare all site types, we conducted Tukey's post hoc analysis (package multcomp; Hothorn et al. 2008).
Question 4: Does site type affect the movement of organisms between the marsh and the upland?
To estimate the effect of a physical barrier on the use of upland habitats by marsh species, we documented use of the upland by the squareback crab Armases cinereum (Armases, henceforth). Because human presence changes crab behavior, one member of the research team walked the upland-marsh border immediately upon arrival at the site. Armases crabs were counted in 1 × 4 m quadrats in the upper marsh and in the upland (− 0.5 m and + 0.5 m from the Bzero^line or bulkhead; n = 3 per site). To verify the robustness of the quadrat counts, we employed a pitfall trap survey (Supporting Information).

Statistical Analysis
We used a generalized linear mixedeffects model to examine the number of Armases counted in the upland (Poisson distributed), with the number of Armases in the marsh as a covariate and station included as a random variable to account for geographic variation in the data (package lme4; Bates et al. 2014).

Field Survey
Spartina alterniflora (Spartina) and J. romerianus (Juncus) dominated the plant community at our sampling sites. However, Spartina cynosuroides, Spartina bakeri, Salicornia depressa, Borrichia frutescens, Iva frutescens, and Schoenoplectus sp. were all observed in at least one plot. The dominant invertebrates in marsh plots were L. irrorata (Littoraria) and crabs. Crab species associated with crab burrows included Armases cinereum, Uca pugnax, Uca minax, Uca pugilator, Sesarma reticulatum, and Eurytium limosum. There was no evidence of sessile marine invertebrates fouling any of the bulkheads. Armases were found on bulkheads at 11 locations, Littoraria were found on five bulkheads, and Uca pugnax were found on two bulkheads. In addition, we found Anolis carolinensis lizards on one bulkhead, Eumeces faciatus or Eumeces inexpectatus (five-lined skink or southeastern five-lined skink) on three bulkheads, and an unidentified snake on one bulkhead.
Question 1: How does site type affect the physical and environmental characteristics of the upper marsh?
Site types varied in terms of elevation and upland-creek distance, both of which were correlated with other environmental variables. Elevation was significantly lower at armored than unarmored or forested sites (Table 2, Fig. 2). Higher elevations were associated with lower salinity, soil organic matter, ammonium and phosphate concentration, and higher percent sand, soil water content, and nitrate concentration (Table 2; Figs. 3 and 4). Upland-creek distance was longer at forested than at armored or unarmored sites (Table 2; Fig. 2). Longer upland-creek distance was associated with higher ammonium and phosphate concentrations, more bare (unvegetated) space, and lower salinity (Table 2; Figs. 3 and 4).
After accounting for the effects of elevation and uplandcreek distance, several significant effects of site type remained. First, salinity was lower at the unarmored sites than at the forested sites (Table 2, Fig. 3A). Second, soil water content was higher at unarmored sites than at armored and forested sites (Table 2, Fig. 3C, d). Third, the fraction bare of vegetation was higher at unarmored sites than at armored sites ( Table 2, Fig. 4I, J). In one case, the interaction of Spartina wrack cover and elevation was significant (Table 2); wrack cover at armored sites was greater at lower elevations, whereas at unarmored sites it was greater at high elevations (Table 2, Fig. 4K).
Question 2: How does site type affect the biological community of the upper marsh?
Four members of the biological community accounted for 70% of the dissimilarity between marsh communities adjacent to the different site types: Spartina, Juncus, Littoraria, and crabs (as indexed by their burrows; Table 3). Although the Table 1 The mean, minimum, and maximum values from the survey for each of the variables by site type. It should be noted that the mean presented here may be misleading and may conflict with the statistical results, as it is the simple arithmetic mean calculated without regard to geographic location (station) or sampling point along the transect (i.e., 2, 4, or 8 m from the upland). For summarized data accounting for sampling point along the transect, see   (Table 4).
Question 3: How does site type affect the relationship between the physical and biological characteristics in the upper marsh?
Seven environmental variables were correlated with the structure of the biological community: elevation, uplandcreek distance, porewater salinity, soil water content, porewater concentrations of nitrate and phosphate, and wrack cover (multivariate ANOVA, Table 5).
Spartina had greater coverage at armored versus unarmored sites (generalized mixed modeling, Table 6, Fig. 5). It also had greater coverage at sites with lower elevations, higher salinity and soil water content, and lower wrack cover and porewater phosphate concentrations (Table 6, Fig. 5). Spartina coverage was lower at stations with longer uplandcreek distances. Juncus coverage was not affected by site type, but there was greater coverage at sites with lower soil water content and wrack cover ( Table 6, Fig. 5). There were more crab burrows at forested sites, followed by armored sites and then unarmored sites. More crab burrows were also found at stations with shorter upland-creek distance and sites with lower elevation, porewater phosphate, and wrack coverage, and L higher salinity, soil water content, and porewater nitrate concentrations (Table 6, Fig. 5). Littoraria densities were not affected by site type, but there were more Littoraria at sites with longer upland-creek distances and with lower nitrate concentration and wrack coverage (Table 6, Fig. 5).
Question 4: How does site type affect the use of terrestrial habitats by organisms that routinely move between the upland and marsh?
The number of Armases found in the upland varied by site type, with the highest densities at forested sites, then unarmored sites, and the lowest counts at armored sites (Table 7, Fig. 6). Additionally, the number of Armases found in the marsh was positively correlated to the number found in the upland (Table 7, Fig. 6).

Discussion
We found that shoreline armoring and development affected the environmental and biological structure of the upper marsh, but the effects were subtle rather than dramatic. Marshes adjacent to bulkhead armoring had lower elevations than those adjacent to unarmored or forested sites (Fig. 2), and the biological and physical characteristics of these marshes were consistent with a lower elevation. In particular, armored sites had greater Spartina coverage and crab burrow abundance than unarmored sites (Fig. 5). In addition, silty clay-a soil characteristic of tidal creek sediments-was only found adjacent to armored sites (Supporting Information).
The lower elevation associated with armored sites may be the result of increased erosion at bulkheaded sites after installation. Bulkheads and riprap without adjacent marsh habitats have also been shown to have lower elevation (measured as water depth; Kornis et al. 2017). However, it is also possible that the lower elevation reflects that the bulkheads themselves were built directly in marshland. Although we do not have preconstruction elevations to directly address this, installing a bulkhead directly in the marsh is against permit regulations that have been in place since 1970. Moreover, many of the bulkheads in this study were paired with neighboring homes in the same subdivision, and there is no a priori reason to expect that the observed differences in elevation occurred before the bulkheads were in place. Regardless of the history, local supply of sand from the upland can be cut off by the presence of the bulkheads, which would affect the ability of the upper marsh to accrete vertically (Edwards and Frey 1977;Frey and Basan 1985). The abrupt jump in elevation from high marsh to upland created by the bulkhead could remove transitional habitat for plant communities. We did not sample plants directly at the marsh-upland boundary, but Bozek and Burdick (2005) found that plant species richness at the marsh-upland boundary in Fig. 4 Mixed-effects model fit of percent soil organic matter (SOM, A, B), porewater ammonium (NH4), nitrate (NO3) and phosphate (PO4) concentrations in milligram/gram (C, H), the fraction of bare soil (Fraction Bare, I, J), and the probability of wrack cover (Pr(Wrack), K, L) as a function of elevation (left column) and upland-creek distance (right column). Armored sites are denoted by black solid lines, unarmored sites by gray dashed lines, and forested sites by light gray solid lines. In contrast to Fig. 3, partial residuals are not shown because these models are binomial fits and log-transformed response variables. Internal tic marks denote the distribution of observations with respect to elevation (left column) or distance (right column) Table 4 Test of the effect of site type (armored (A), unarmored (U), and forested (F)) and sampling point along the transect (2, 4, or 8 m from the upland) on the salt marsh biological communities as measured by multivariate ANOVA (package adonis)  Table 3 The cumulative contributions of the top four biological variables, Spartina alterniflora and J. romerianus (% cover), and crab burrows and L. irrorata (number m −2 ), that together account for over 70% of the dissimilarity between marshes adjacent to the site types (as quantified using a SIMPER analysis)  Table 5 Relationship between salt marsh environmental variables and biological community structure (represented by MDS output), indicating which variables are significantly correlated with differences in the biological community as measured by multivariate ANOVA. Significant variables were included in the Question 3, New Hampshire was reduced by 50% in areas with rock bulkheads, and it is likely that a similar effect occurred at our sites. The most dramatic difference we found among site types was the effect on Armases, which is the one species we measured that readily moves across the ecotone. Armases was most abundant in the upland associated with the forested sites, with intermediate densities in unarmored sites and the fewest in armored sites ( Fig. 6 and Table 7). These results were supported by a pitfall trap study that showed that Armases moved further into the upland at forested sites than at the other site types (Supporting Information). Armases likely prefers heavily wooded areas because it experiences less desiccation in shaded, cooler habitats, which may explain the decrease in Armases found at the unarmored and armored sites. Although field observations suggested that Armases is not able to climb some bulkhead materials, such as vinyl siding, only 7 of our 20 bulkheads were constructed with vinyl siding and Armases were regularly found on and inside wooden bulkheads. There are several other species that move across the upland-high marsh ecotone, including butterflies, grasshoppers, birds, and raccoons; future studies should evaluate how these species are affected by upland development (with and without armoring).
Contrary to our hypothesis that forested, unarmored, and armored sites would present a response gradient, we found that unarmored development, i.e., without a protective bulkhead, often had different, opposing effects on the upper marsh community than armored development. For example, wrack deposition adjacent to armoring was low and increased at lower elevations farther from the bulkhead, whereas at unarmored sites wrack deposition was enhanced with increasing elevation (Table 2, Fig. 4K). Wrack acts as a disturbance in salt marshes, as it can smother the underlying vegetation if it is present for a long enough period of time (Bertness and Ellison 1987;Valiela and Rietsma 1995;Li and Pennings 2016;Li and Pennings 2017), and in fact higher wrack cover was associated with lower Spartina and Juncus coverage and a decreased density of crab burrows and Littoraria (Table 6, Fig.  5). However, it is also a subsidy in that it provides a source of food and shelter for invertebrates such as isopods and amphipods, which represent food for higher trophic levels (e.g., spiders; Zimmer et al. 2002;Buck et al. 2003). The decrease in wrack associated with bulkheads observed here is similar to patterns that have been found adjacent to armored structures on open-coast beaches, where armored beaches poorly retained wrack deposits, subsequently leading to lower bird populations (Dugan et al. 2008;Sobocinski et al. 2010).
Marshes adjacent to unarmored sites in this study were characterized by higher soil water content than either armored or forested sites and had lower salinity than forested sites ( Table 2, Fig. 3A, C). Taken together, this suggests that overland and groundwater input of freshwater to the marsh could be increased by development in the absence of a bulkhead. Increased freshwater input has been associated with unarmored development in South Carolina as well as the Table 6 The effects of site type (armored (A), unarmored (U), and forested (F)) and the top environmental variables correlated with change in the biological community structure (from Table 5) on the top four species driving the changes in the biological communities (from Table 3) as tested through mixedeffects modeling. Variance is reported for the random variable of station, sampling point along the transect (2, 4, or 8 m from the upland) nested in station, and βcoefficients are reported for the fixed variables. Significant coefficients are indicated by * (i.e., p < 0.05) northeastern US (Silliman and Bertness 2004;Walters et al. 2010). Armoring could ameliorate this effect by blocking water flow paths from the upland to the marsh or by decreasing runoff from the upland by reducing the land-surface gradient adjacent to the marsh. If so, it could be possible to alleviate the effect of unarmored development on the upper marsh by reducing freshwater input in other ways, for example through tying into a sanitary sewer system, collecting precipitation with rain barrels, and limiting the watering of lawns. Forested sites were characterized by increased porewater salinity, higher density of crab burrows, and longer upland-creek distances compared to other site types (Tables 2 and 4, Figs. 2, 3, and 5). The longer upland-creek distances observed at forested sites may reflect a predilection for development to occur in areas with shorter distances to the water, making it easier to install docks. Although the blocked design of this study should have controlled for large-scale geographic variability in marshupland border characteristics, we could not control for preexisting differences among site types within a local area. This highlights the limitations of this type of field survey. If possible, we encourage coastal managers to require studies before development is initiated, enabling a before-after control-impact design that would better isolate the effects of development on marsh communities (Stewart-Oaten et al. 1986).
The effects of armoring and development on the upper marsh ecosystems in this study were characterized by subtler changes than those previously reported in the northeastern US or in open-coast systems (Wahl et al. 1997;Bertness et al. 2002;Dugan et al. 2011). There are several possible reasons for this difference. First, marsh ecosystems are relatively low energy compared to the open-coast systems where other armoring research has been conducted. Effects of shoreline armoring across soft-sediment environments appears to vary by energy and armoring type, and our work supports the hypothesis that the effects of armoring will be subtler in low energy systems (Dugan et al. 2017;Bozek and Burdick 2005). Second, the southeastern US coastline is relatively undeveloped (Crossett et al. 2005;Gittman et al. 2015). Thus, the southeast may not yet exhibit the cumulative effects of development that are present in highly developed coastal regions like New England (Walters et al. 2010). Third, we standardized the selection of bulkheads in our study, sampling only at structures that had not obviously reclaimed upland from the marsh. This excluded some of the structures with the greatest potential impacts, such as cases where the filled and armored area completely replaced the upper marsh, thereby removing this ecotone entirely. We eliminated these cases during site selection because such placement is no longer permitted in Georgia and because we were interested in studying the effects of bulkheads as barriers per se, without associated habitat destruction.
The coastal southeastern US is expected to see intense development in the future (Crossett et al. 2005). The impacts we measured could therefore become more widespread and increase in magnitude as more bulkheads are built. Our results suggest that the bulkheads block the supply of sand from upland areas, potentially resulting in vertical loss of elevation in the upper marsh (Edwards and Frey 1977). Moreover, the presence of a physical barrier limits the ability of marshes to migrate horizontally onto the upland, resulting in what is known as Bcoastal squeeze^ (Doody 2004). In addition to limiting armoring, regulatory agencies can work to minimize the effects of residential development that we observed. For example, forest vegetation could be retained along the high marsh ecotone to provide a buffer that could serve to minimize freshwater runoff into the marsh and provide room for potential upward marsh migration. Forested vegetation could also potentially provide habitat for organisms such as Armases that routinely move between the marsh and the upland. Because the southeastern US presently has lower population densities and a limited extent of armoring in coastal marsh environments, proactive policy in this part of the country could help prevent the stronger effects of armoring and development that have been seen elsewhere. Table 7 Generalized linear mixed-effects model evaluating Armases movement into the upland as a function of the number found in the marsh, and site type (A = armored, U = Unarmored, F = forested). Variance is reported for the random variable of Station, and β-coefficients are reported for the fixed variables. Significant coefficients are indicated by * (i.e., p < 0.05)

Variable Armases in marsh U/A F/A U/F Station Distribution
Armases 0.018* 1.44* 0.65* −0.8* 1.88 Poisson  Fig. 6 Generalized linear mixed-effects model of the total number of Armases cinereum expected to be found in the upland as a function of the number of Armases found in the marsh and the site type, with station as a random variable. Armored sites are denoted by black solid line, unarmored sites by gray dashed line, and forested sites by light gray solid line