Maintaining food production while sustaining productive ecosystems is among the central challenges of our time, yet, it has been for millennia. Ancient clam gardens, intertidal rock-walled terraces constructed by humans during the late Holocene, are thought to have improved the growing conditions for clams. We tested this hypothesis by comparing the beach slope, intertidal height, and biomass and density of bivalves at replicate clam garden and non-walled clam beaches in British Columbia, Canada. We also quantified the variation in growth and survival rates of littleneck clams (Leukoma staminea) we experimentally transplanted across these two beach types. We found that clam gardens had significantly shallower slopes than non-walled beaches and greater densities of L. staminea and Saxidomus giganteus, particularly at smaller size classes. Overall, clam gardens contained 4 times as many butter clams and over twice as many littleneck clams relative to non-walled beaches. As predicted, this relationship varied as a function of intertidal height, whereby clam density and biomass tended to be greater in clam gardens compared to non-walled beaches at relatively higher intertidal heights. Transplanted juvenile L. staminea grew 1.7 times faster and smaller size classes were more likely to survive in clam gardens than non-walled beaches, specifically at the top and bottom of beaches. Consequently, we provide strong evidence that ancient clam gardens likely increased clam productivity by altering the slope of soft-sediment beaches, expanding optimal intertidal clam habitat, thereby enhancing growing conditions for clams. These results reveal how ancient shellfish aquaculture practices may have supported food security strategies in the past and provide insight into tools for the conservation, management, and governance of intertidal seascapes today.

Funding: This work was funded by a National Geographic Research and Exploration Grant (8636-09) to Dana Lepofsky, a Natural Sciences and Engineering Research Council of Canada Discovery Grant (385921-2009) to Anne K. Salomon, a Social Sciences and Humanities Research Council Insight Grant (2011-0833) to Dana Lepofsky, and a David and Lucile Packard Foundation Grant (2008-32497) to Kirsten Rowell. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

In this study, we quantified the productivity of ancient clam gardens on Quadra Island, BC with surveys and an in situ transplant experiment. We measured how bivalve communities and beach morphology differs between clam gardens and non-walled beaches, and what environmental factors contribute to these differences. Specifically, we ask, do clam gardens have higher clam densities, biomass, and growth rates compared to non-walled beaches? And if so, what physical characteristics best explain these differences?

Researchers have recently turned their attention to ancient clam gardens along the Northwest coast highlighting another example of ancient marine enhancement and management. These human-engineered intertidal terraces, thought to have been constructed in the late Holocene, have been recorded from Alaska through British Columbia (BC) and into Washington State [28] . Made by building rock walls in the low intertidal of soft-sediment clam beaches, clam gardens are thought to stabilize sediments at a specific tidal height, presumably to enhance shellfish productivity [28] , [29] , [30] , [17] . Clam gardens walls exist at the mouths and along the edges of embayments, parallel to coastlines, and vary in shape, length, width, and intertidal height ( Fig 1A-D ). Although the age of these ancient features is currently unresolved, the immense shell middens associated with clam garden walls suggest the significance of shellfish as a staple food source for Northwest Coast First Nations for at least 5000 years [31] . The combination of the widespread occurrence of clam gardens on the northwest coast, their associated shell middens, and traditional knowledge of clam garden tending passed down in song, story, and practice [32] , [16] , underscores the importance of these features and suggests that they were constructed to increase clam yields. Knowing how and the extent to which clam gardens boost clam yields may offer insights into contemporary investments in food security for coastal communities.

Several ancient environmental engineering and resource management strategies have been documented among coastal indigenous peoples. Enforced harvest size restrictions of forest resources in Fiji [18] ; construction and tending of root gardens [16] ; transplantation and cultivation of berries to increase yields [19] ; and prescribed burns to clear land and magnify plant production [20] along the Northwest Coast provide examples of intentional management of terrestrial coastal resources. Globally, marine examples include stone ponds for fish aquaculture in Polynesia [21] ; complex wooden fish weirs in Brittany [22] ; and temporal and age-specific harvest restrictions to conserve reef fish in Oceania [15] . Along the northwest coast of North America, marine examples include intertidal stone fish traps at the mouths of salmon bearing streams and intertidal wooden fish weirs [23] , [24] in conjunction with size selective fishing practices to enhance salmon productivity [25] and the reduction of predatory sea otters to increase shellfish abundance [7] . These records of direct environmental manipulations, tending, and stewardship practices suggest complex systems of resource management which increased local food security [5] , [26] , [27] .

Humans have been altering, exploiting, and managing marine and terrestrial ecosystems for millennia [5] , [6] , [7] , [8] . Throughout history, human hunting and fishing in coastal ecosystems has caused declines in key species [9] , reduced prey size [10] , [11] , triggered trophic cascades [12] , [13] , and facilitated ecosystem regime shifts [11] . In other cases, the archaeological record indicates long term sustained yields, with little indication of resource depression [14] . Recent archaeological evidence and oral historical knowledge suggests that First Peoples around the world actively managed and enhanced nearshore ecosystems to maintain and increase productivity [15] , [16] , [17] , [8] .

Sustaining global food production presents one of the greatest environmental and humanitarian challenges of the 21st century. Given current global population and consumption trajectories, the world's food production must double by 2040 [1] , [2] and its footprint must shrink substantially to reduce the degradation of land, water, biodiversity, and climate. Consequently, society will need to develop clever ways to meet demands on terrestrial and marine resources and spaces efficiently, while maintaining ecosystem productivity and resilience. Fortunately, evidence from the past often offers solutions to contemporary quandaries [3] , [4] . Here, we provide empirical evidence of an ancient form of mariculture that magnified shellfish production in a limited space, providing practical insights into sustainable marine management techniques which may inform local food security strategies of today.

All field work was conducted with the following permits: a Contaminated Shellfish Collection for Scientific Purpose Permit from Department of Fisheries and Oceans Canada, a Park Use Permit from BC Parks, an Animal Care Protocol Exemption from Canadian Council on Animal Care for working with invertebrates, and permission from Laichwiltach Treaty Society. The individual pictured in this manuscript ( Figure 1C and Striking Image) has given written informed consent (as outlined in PLOS consent form) to publish their likeness. Our study did not involve endangered or protected species. All necessary permits were obtained for the described study, which complied with all relevant regulations.

We tested for differences in survival and growth rates of L. staminea transplanted in clam gardens (n = 6) and non-walled clam beaches (n = 5) in Waiatt Bay, using generalized GLMMs where beach type was a fixed effect and site was a random effect. To test for differences in L. staminea growth and survival across tidal stations within clam gardens and non-walled beaches, we constructed the same GLMMs as above, with beach type, bay and beach type*tidal station as fixed effects. Growth rates models were fit with REML, a Gaussian error distribution, and identity link function using the lme function in the nlme package [34] . Survivorship models were fit with a Laplace Approximation, a binomial error distribution, and logit link function using the lmer function and lme4 package [35] . Survivorship and growth were compared among different juvenile clam size classes using the same models described above on three size classes of transplanted L. staminea binned by 5mm increments. Bin size was determined by the size of the transplanted clams and number of clams in each bin that would yield sample sizes large enough for sufficient statistical power.

To assess if and how clam garden engineering altered intertidal height and optimal growing conditions for clams, we modeled the relationship between intertidal height and a) density and biomass of surveyed L. staminea and b) survivorship and growth of transplanted L. staminea in clam gardens and non-walled beaches in both bays, by fitting Gaussian models (Eq.1) to each metric of clam productivity (y) as a function of intertidal height where: (Eq.1)α (curve height) describes the magnitude of clam productivity (biomass, density, or growth rate), µ (curve mean) is the intertidal height at which productivity is greatest, and σ (curve width) describes the standard deviation in clam productivity. We then compared fitted model parameters across clam gardens and non-walled beaches in both bays based on our predictions ( Fig 3 ).

To test for differences in clam density and biomass between clam gardens (n = 11) and non-walled clam beaches (n = 10), we constructed general linear mixed effects models (GLMMs) where beach type was treated as a fixed effect and site was treated as a random effect. These models were constructed for the three most dominant species (L. staminea, S. giganteus, and Macoma spp.) and total clams. To test for differences in L. staminea density among different size classes, we ran the same models described above on five size classes of clam binned into categories based on 12 mm increments. Bin size was determined by the size range of the surveyed clams and number of clams in each bin that would yield sample sizes large enough for sufficient statistical power. Differences in clam density and biomass between clam gardens (n = 11) and non-walled clam beaches (n = 10) as a function of tidal station in both Kanish and Waiatt Bay were assessed using the same GLMMs as described above, with the additional fixed effects of bay and the interaction of beach type*tidal station. Beach type, bay, and type*tidal station were specifically chosen as treatments to be tested, beach type to detect a clam garden effect, bay to detect an effect of oceanographic context, and type*tidal station to detect our predicted across-beach effect of tidal station in clam gardens and non-walled beaches ( Fig 3b ). Clam biomass models were fit with a REML, a Gaussian error distribution, and identity link function using the lme function in the nlme package [34] . Clam density models were fit with Laplace Approximation, a Poisson error distribution, and log link function using the lmer function and lme4 package [35] . All GLM and GLMM modelling was conducted in R [36] .

To test for an effect of beach type (i.e., clam garden [n = 11] vs. non-walled beach [n = 10]) and an effect of bay (Kanish vs. Waiatt Bay) on beach slope, we used a general linear model (GLM). We used the same strategy to examine differences in heights of clam garden walls and mean heights of garden terraces between bays. Models of Slope, Wall Top Height, and Terrace Height were fit with a Restricted Maximum Likelihood (REML), a Gaussian error distribution, and identity link function using the lme function in the nlme package [34] .

We selected clam garden sites based on the presence of a complete rock wall spanning an embayment. Non-walled control beaches were beaches that lacked a rock wall but encompassed clam habitat, specifically intertidal soft sediment and the full array of clam species observed in the area. When comparing non-walled beaches and clam gardens, we assume that contemporary non-walled beaches are representative of the pre-construction state of clam garden sites.

By constructing a rock-walled terrace in the intertidal, we predicted that clam gardens expand optimal clam habitat by altering the slope of the beach and thus increasing clam habitat at the intertidal height at which clams grow and survive best ( Fig 3a ). In response to increased habitat and enhanced conditions we predicted that clam density, biomass, growth, and survival would increase at the first and last tidal stations ( Fig 3b ), i.e., the extreme high and low of intertidal heights ( Fig 3c ), in clam gardens compared to non-walled beaches. We also predicted that Gaussian models of optimal clam productivity as a function of intertidal height would peak at the same intertidal height in both clam gardens and non-walled beaches, but that the magnitude of productivity would be greater (i.e. due to increased water retention, differences in sediment composition, or both) and the variance smaller (i.e. due to the reduction in beach slope) in clam gardens relative to non-walled beaches ( Fig 3c ).

We chose L. staminea for our experiment because they are a biologically and culturally significant species in the area, and they possess distinct annuli, which have been verified by the annual temperature-driven δ 18 O signature observed between annuli (K Rowell, pers. comm.). For the same reasons L. staminea will also be used in a companion archeological study comparing ancient clam growth rates through deep time, both before and after clam garden construction.

To test if clam gardens increase the growth rates of L. staminea, we conducted a transplant experiment across six clam gardens and five non-walled clam beaches during the growing season from May to October 2011 [33] . Clams 11–34 mm in length were collected from Waiatt and Kanish Bay, labelled with two uniquely numbered vinyl tags, measured to the nearest 0.1 mm and weighed to the nearest 1.0 gram. Fifteen individuals, representing the range of sizes collected, were placed inside a 34 cm×24 cm Vexar mesh bag. Five bags of L. staminea were evenly spaced at five tidal stations along a single vertical transect, perpendicular to shore, from the top of clam habitat to ∼0 m intertidal height at non-walled beaches and at the edge of rock wall in clam gardens ( Fig 3a ). Each bag was buried approximately 10 cm below the surface, based on our natural history observations for clams of this size and recommendations from a contemporary aquaculture facility. Each transplant bag was then secured with a flagged and labelled rebar stake. Transplanted clams were left in situ for 160 days. Upon retrieving the transplanted clams, max lengths, widths, and weights were recorded. We noted all losses, mortalities, and evidence of predation.

A) We surveyed clams across a vertical transect (black line) from the top of clam habitat (H) to ∼1.0 m in non-walled beaches (left) and to the edge of the rock wall (W) in clam gardens (right). We then transplanted clams in mesh bags at 5 evenly stratified tidal stations (blue and green colored squares) across a vertical transect (purple line) from the top of clam habitat (H) to ∼0 m in 5 Non-Walled Beaches (left) and to the edge of the rock wall (W) in 6 clam gardens (right). B) Hypothesis 1: predicted clam productivity as a function of tidal station. Tidal station 1 = top of clam habitat, tidal station 15 = top of clam garden wall in clam gardens or ∼1.0 m tidal height in non-walled beaches. C) Hypothesis 2: predicted clam productivity as a function of intertidal height.

We characterized beach slope at 11 clam gardens and 10 non-walled beaches to quantify the physical differences in intertidal clam habitat between clam gardens and non-walled clam beaches. At each site, we established a vertical transect, perpendicular to shore, with 15 randomly stratified stations. The tidal height of each station was quantified using a total station or laser level in meters above Canada chart datum, lowest low water large tide (LLWLT). Transects began at the highest intertidal height at which clams were found in test pits. In clam gardens, the bottom tidal height of each transect was anchored by the landward edge of the human-made rock wall. At non-walled clam beaches, the bottom tidal height of each transect was anchored at ∼1.0 m above chart datum, the mean tidal height of the landward edge of the clam garden rock-walls.

We located 5 non-walled clam beaches in each bay that we deemed appropriate controls for our study. We then chose at least 5 clam gardens in each bay as treatments for comparison. The total number of non-walled beaches (n = 10) and clam gardens (n = 11) surveyed for physical and biological comparison was constrained by the number of low tide days available during the spring and summer field season of 2011.

We conducted our research on northern Quadra Island in British Columbia (BC), Canada, where an exceptionally high density of clam gardens have been documented [30] in Kanish (n = 45 clam gardens) and Waiatt Bays (n = 49 clam gardens) ( Fig 2 ). Quadra Island has an abundance of archaeological sites found throughout the landscape, with shell middens representing both permanent settlements and short term camps. Today, the northern part of the island falls within the traditional territories of the Northern Coast Salish and the Southern Kwakwaka’wakw (now Laich-kwil-tach) First Nations. Some of the descendants of the ancient settlements live in nearby Indian Reserves and town centers. Presently, Kanish and Waiatt Bays are only sparsely settled and bordered by second growth forests and active wood lots. The bays are popular recreation areas and anchorages, encompassing two provincial parks and an active scallop farm. Clam digging, once a mainstay of the dense human population in the bays, is now only conducted recreationally and sporadically.

Transplanted L. staminea grew significantly faster in clam gardens than non-walled beaches, and this effect varied as a function of tidal station ( Fig 8 , S3 , Table 4 , S3 ). Clams grew proportionally faster at tidal station extremes (the first and last tidal station) in clam gardens compared to non-walled beaches. In line with our expectations, the overall magnitude of growth rates as a function of tidal height was higher in clam gardens than non-walled beaches (Fig7e, S4 , Table 3 ). Size appears to be a major predictor of survivorship - small size classes of L. staminea (11−16 and 17−22 mm) were more likely to survive in clam gardens than non-walled beaches, although clam garden habitat did not appear to effect survivorship when all size classes were pooled. This size dependent effect varied as a function of tidal station ( Table 4 ).

As predicted ( Fig 3c ), the magnitude of L. staminea productivity (a) in terms of density, biomass, and growth, was higher in clam gardens than non-walled beaches, and the standard deviation (σ) was lower ( Fig 7a-d , Table 3 ). Contrary to expectations, the intertidal height at which little neck clams reach their maximum density, biomass, and growth (µ) was consistently higher in clam gardens than non-walled beaches in Kanish Bay ( Fig 7 , S4 , Table 3 ), suggesting that optimal habitat shifted ∼0.5 m higher up the beach in clam gardens in Kanish. In Waiatt, maximum L. staminea productivity (µ) did not differ substantially between beach types, and biomass did not conform to a Gaussian relationship within non-walled beaches.

By examining clam density and biomass as a function of tidal station and beach type in both bays, we found that L. staminea and S. giganteus densities and biomass were significantly greater in clam gardens than non-walled beaches and as predicted ( Fig 3b ), this relationship varied as a function of tidal station ( Fig 6 ). This is reflected by the significant interaction terms in Table 2 . Specifically, clam densities and biomass tended to be higher at the first 6–7 tidal stations. The effect of tidal station position was highly significant for total clam densities and densities of L. staminea (all sizes, 14–26 mm, 27–39 mm) ( Fig 6 , Table 2 ), S. giganteus, and V. philippinarum ( Table 2 ). The density and biomass of the invasive V. philippinarum varied significantly as a function of tidal station and was significantly higher in Kanish Bay ( Table 2 ).

We observed greater variation in the intertidal height of clam gardens in Waiatt Bay. There, clam garden terraces were located between 0.53–1.45 m above chart datum, with four of the six clam gardens having mean terrace heights between 0.78–1.16 m, and two outliers having mean terrace heights of 0.53 m and 1.45 m respectively ( Fig 4 , S1 ).

On Quadra Island, BC, clam garden terrace heights varied between bays, but their slopes were consistently shallower than unaltered beaches ( Fig 4 , S2 a, Table S1 , F (1,19) = 6.914, p = 0.017). On average, mean intertidal heights of clam garden terraces in Waiatt Bay were significantly lower than those in Kanish Bay ( Fig 4 , S1 , S2 a, Table S1 , F (1,9) = 15.848, p = 0.003). In Waiatt Bay, clam garden terraces were located on average at 0.97 m (+/– 0.31 SE) above chart datum, while the tops of rock wall features averaged 0.68 m (+/– 0.36 SE) in intertidal height ( Fig 4 , S1 ). In contrast, Kanish clam garden terraces were located on average at 1.57 m (+/–0.21 SE) above chart datum and the tops of rock wall features averaged 1.3 m (+/–0.19 SE) in intertidal height ( Table S2 ). Non-walled beach slopes and mean intertidal heights did not differ between bays ( Fig 4 , S2 a,b).

Discussion

Strong evidence from both our surveys and experimental transplants, suggests that ancient clam gardens increased clam productivity. By altering the slope of soft-sediment beaches (Fig 4), these human-made, intertidal terraces expanded the optimal intertidal habitat and enhanced growing conditions for clams. Specifically, we detected significantly greater densities of S. giganteus and L. staminea in clam gardens compared to non-walled beaches, particularly among smaller size classes of pre-reproductive clams (Fig 5). Overall, clam gardens contained 4 times as many butter clams (S. giganteus) and over twice as many little neck clams (L. staminea) relative to non-walled beaches (Fig S2). As predicted, the magnitude of this relationship varied as a function of intertidal height, whereby clam density and biomass was enhanced in clam gardens compared to non-walled beaches at the top and bottom of the beach, the areas where clam gardens extend optimal clam habitat (Fig 6,8). The pattern of increased clam productivity by clam gardens appears to be driven by the modification of intertidal height (Fig 7) and is supported by our experimental results indicating higher L. staminea growth rates within clam gardens (Fig 7e,8). Even though clam gardens on Quadra Island have not been actively tended or managed for decades, we detected significant signals of enhanced shellfish production across these engineered beaches simply due to their modified slopes. Furthermore, elevated clam densities, biomass, and growth rates at equivalent intertidal heights in clam gardens compared to non-walled beaches suggests that additional mechanisms, in addition to tidal height modification, appear to be magnifying secondary production in these clam terraces.

Mechanisms Enhancing Productivity in Clam Gardens In addition to altering beach slope and thereby extending the area of intertidal habitat at the optimum tidal height for clam survival and growth, clam gardens terraces may have enhanced clam productivity in multiple ways. For example, over the diurnal tidal exchange, we observed increased water retention over clam gardens relative to natural sloping non-walled beaches. Water retention may increase the opportunity and success of larval clam recruitment and survivorship in a clam garden. In fact, a Hul’qumi’num First Nation reported that clam gardens were a way to “trap the seeds and keep them here” [27]. Low water velocities create optimal conditions for larval settlement and recruitment [37]. Water retained over shallow sloping clam gardens could increase in temperature in the late spring and summer months thereby increasing phytoplankton growth rates [38] and fuelling secondary production[39], [40]. Furthermore, increased water temperatures in temperate intertidal systems are known to enhance bivalve growth rates [41], [42] and trigger bivalve spawning events [33]. By intentionally modifying substrate, reducing density dependence, and excluding competitors and predators in clam garden terraces, indigenous people may have enhanced clam productivity in clam gardens even further. On Quadra Island, we observed that clam garden substrate tended to be higher in gravel and shell hash content compared to non-walled beach substrate, which tended to have more fine sediments. Similar observations have been reported elsewhere [23], [28] and L. staminea, are commonly found on natural beaches in course sand or fine gravel mixed with mud, stones, or shells [43], [44]. According to a Heiltsuk First Nation knowledge holder; "We put gravel in the garden to increase the number of clams" (D Wilson, pers. comm.). Increasing the gravel content in clam gardens may have created larger interstitial spaces in the substrate which is likely important for porewater flow and reduced fine silts and clays that are known to smother newly settled L. staminea larvae [43]. The act of aerating beach sediments by rolling rocks, or “turning over beaches” is also commonly reported [32], [5], [27] and aims to reduce anoxic conditions that can reduce productivity. In addition, Hul’qumi’num First Nation knowledge holders report returning crushed and whole clam shells to clam gardens as a management practice [27] and adult bivalve shell has been shown to offer an important settling cue for shellfish like oysters and clams [45], [46]. In fact, it has been demonstrated that gravel and shell increase clam settlement and survival [47]. Finally, Turner [5] reports that clam gardens were “thinned,” reducing densities of large adult clams via harvest, giving smaller clams the space and resources to grow, and thus increasing overall yields [48]. We also hypothesize that predators such as sea stars (Pycnopodia helianthoides, Pisaster brevispinus.), large crabs (Metacarcinus magister) and mammalian coastal predators (river otters, sea otters), may have been intentionally excluded from these gardens to decrease both direct predator mortality and negative non-lethal predator effects on clam productivity [49], [50]. In addition to these ecological factors, several key social factors, including systems of tenure and control, may have equally enhanced and maintained the productivity of clam gardens. Indigenous peoples of the northwest coast had territorial governance systems and complex protocols that delineated access rights to the land and sea [51], [5], [52]. Clam gardens, like seaweed picking areas, root gardens, and fish traps, would have been embedded in these traditional systems of marine governance and tenure. Among the Heiltsuk First Nation, families owned and tended productive clam beaches (D Wilson, pers. comm.). Building and maintaining clam gardens were intentional acts, clearly showing cultural investment. Territorial access rights, via family-based proprietorship, established a governance system over common pool fisheries resources that likely conferred resilience to societies on the Northwest coast for millennia [52], [8]. Similarly, empirical evidence from contemporary fisheries management highlights the importance of designating access rights to enhance resource sustainability [53], [54].

Assumptions Documenting ancient resource management within contemporary landscapes presents a challenge for various reasons, one being the identification of adequate controls. Comparisons of areas with and without archaeological features are complicated by the uncertainty in why suitable land and seascapes were not modified. In this study, we assumed naturally shallow sloping non-walled clam beaches are appropriate clam garden controls; however a variety of alternative and non-mutually exclusive hypotheses could explain why some beaches were left unmodified. For example, non-walled beaches could have been of poor quality habitat, used for other purposes, owned by other title holders or may have been at an inappropriate tidal height due to past sea level change. To further examine the potential productivity of clam gardens and what constitutes a true control, several experiments and surveys could be performed. We recommend sampling clam garden beach sediment cores to quantify pre-modification beach characteristics. Indeed, building experimental clam gardens or deconstructing them, while quantifying ecological and physical responses before and after the modification at both treatment and control sites would be the ultimate way to test the effects of clam gardens on clam productivity. The experimental modification of substrate (i.e. adding shell hash), elimination of clam predators (i.e. Pycnopodia helianthoides and Pisastser brevispinus), and reduction of density dependence (i.e. reducing the density of juvenile clams) could also help tease apart other detailed mechanisms driving increased productivity in clam gardens.

Enhancing Food Security Confers Resilience to Social Ecological Systems An increased appreciation of the coupling between ecosystems and human well-being has triggered a paradigm shift in the applied ecological sciences towards a focus on understanding the dynamics of coupled social-ecological systems (SES), linked systems of people and nature [55], [56], [23]. In many marine systems, current management approaches have demonstrably failed to halt or reverse fisheries declines [57], in part due to the inadequate recognition of the strong links between social and ecological processes [58], [59]. Ancient clam gardens and their governance by coastal communities are an example of an adaptive strategy that likely enhanced regional food security and thus conferred resilience to these coupled human-coastal ocean ecosystems. Our observations on the variation in this ancient form of mariculture also highlight key aspects of resilient social-ecological-systems. The general uniformity in tidal height of clam garden walls on Quadra -- which likely date to different time periods and were owned and managed by different social groups -- reflects knowledge that was shared inter-generationally and across communities. As for variation, we suggest that the two clam gardens that fell outside the optimal tidal height range could be representations of engineering errors and learning, or perhaps, were built during times of differing sea levels. Alternatively, these features may have been built to target other clam species or may have had purposes other than shellfish harvest. Food security is not only a contemporary issue. It has motivated ingenuity and development of civilizations throughout time. Investigations of how ancient clam gardens work will provide information on possible solutions to local food security and economic resiliency of coastal communities. Based on our clam surveys on Quadra Island, densities of L. staminea and S. giganteus within clam gardens are elevated on average by 151% and 300%, respectively, within clam gardens (Table S2). Clam garden biomass of L. staminea and S. giganteus were elevated on average by 68% and 253%, respectively (Table S2). Estimates made from our experiments indicate that clam gardens within optimal clam habitat can enhance growth rates of L. staminea on average by 89%, meaning that clams reach harvestable size at a faster rate (Table S2). If we had chosen to transplant S. giganteus, we predict that we would have detected higher growth rates in clam gardens compared to non-walled beaches, but isolated to the lower half of clam gardens, within S. giganteus optimal habitat. The archaeological record is clear; abundant shellfish have supported large populations of people on the Northwest Coast through history [60], [8]. This new evidence helps emphasize the value of incorporating traditional management techniques into future strategies towards sustainable solutions, contributing to local food security efforts globally.