Enhanced Weathering and Erosion of a Cohesive Shore Platform Following the Experimental Removal of Mussels

The organisms inhabiting intertidal platforms can affect their weathering and erosion rates. Research on biotic influences on platform integrity has traditionally emphasized the role of bioeroders (i.e., organisms that scrap or bore into platforms via mechanical and chemical means). Yet, recent studies illustrate that covers of sessile organisms on the surfaces of intertidal platforms can have bioprotective effects by reducing the efficacy of physical weathering and erosion agents. Eroding cliffs fronted by cohesive shore platforms are a pervasive feature along the continental Argentinean coastline (37–52°S). In this study, we investigated how mussel (Brachidontes rodriguezii) cover mediates weathering and erosion of a cohesive, consolidated silt platform at Playa Copacabana (5 km north of Miramar, Buenos Aires Province; 38° 14′ S, 57° 46′ W). By means of mussel removal experiments, we found that mussel cover attenuates variations in platform surface temperatures, enhances moisture retention during low tide, reduces rates of salt crystallization within the pores of the platform material, and attenuates hydrodynamic forcing on the platform surface. Mussel removal also led to a 10% decrease in surface hardness and a 2-mm reduction in platform height after 5 months. Collectively, our findings indicate that mussel beds limit substrate breakdown via heating-cooling, wetting-drying, and salt crystallization and provide some of the first experimental field evidence for the direct impacts of biotic cover on platform erosion. As intertidal platforms protect the cliffs behind from the hydraulic impact of waves, which may be enhanced with future sea-level rise, we posit that the protection of platforms by mussels indirectly moderates coastline retreat, especially on soft cohesive shores.

The organisms inhabiting intertidal platforms can affect their weathering and erosion rates. Research on biotic influences on platform integrity has traditionally emphasized the role of bioeroders (i.e., organisms that scrap or bore into platforms via mechanical and chemical means). Yet, recent studies illustrate that covers of sessile organisms on the surfaces of intertidal platforms can have bioprotective effects by reducing the efficacy of physical weathering and erosion agents. Eroding cliffs fronted by cohesive shore platforms are a pervasive feature along the continental Argentinean coastline (37-52 • S). In this study, we investigated how mussel (Brachidontes rodriguezii) cover mediates weathering and erosion of a cohesive, consolidated silt platform at Playa Copacabana (5 km north of Miramar, Buenos Aires Province; 38 • 14 S, 57 • 46 W). By means of mussel removal experiments, we found that mussel cover attenuates variations in platform surface temperatures, enhances moisture retention during low tide, reduces rates of salt crystallization within the pores of the platform material, and attenuates hydrodynamic forcing on the platform surface. Mussel removal also led to a 10% decrease in surface hardness and a 2-mm reduction in platform height after 5 months. Collectively, our findings indicate that mussel beds limit substrate breakdown via heating-cooling, wetting-drying, and salt crystallization and provide some of the first experimental field evidence for the direct impacts of biotic cover on platform erosion. As intertidal platforms protect the cliffs behind from the hydraulic impact of waves, which may be enhanced with future sea-level rise, we posit that the protection of platforms by mussels indirectly moderates coastline retreat, especially on soft cohesive shores.

INTRODUCTION
Shore platforms are horizontal or gently sloping rock surfaces usually backed by cliffs that are widespread throughout the world's intertidal zones (Trenhaile, 1987;Sunamura, 1992). They have traditionally been viewed as being shaped by waves (e.g., hydraulic action as well as abrasion by mobile clasts), subaerial weathering (i.e., in situ breakdown of rock driven by climatic variations and cycles of tidal exposure), or some combination of the two (Trenhaile, 2002;Mottershead, 2013;Coombes, 2014;Moses, 2014). In comparison, biotic influences on shore platform dynamics remain relatively understudied, especially given the diversity of taxa and biogeomorphic mechanisms that affect them (see reviews in Naylor et al., 2012;Coombes, 2014). Indeed, research on shore platform biogeomorphology has placed focus on the role of bioeroders (i.e., organisms that scrape consolidated substrates or bore into them via mechanical and chemical means; see Healy, 1968;Trudgill, 1987;Spencer, 1988;Pinn et al., 2005;Bagur et al., 2014;Coombes, 2014;Dodge-Wan and Nagarajan, 2020), whereas the potential "bioprotective" roles (sensu Carter and Viles, 2005) of sessile organisms that cover the surfaces of shore platforms have been considered only relatively recently (e.g., Coombes et al., 2013Coombes et al., , 2017La Marca et al., 2014;Gowell et al., 2015;Naylor and Coombes, 2015).
The evidence available to date indicates that covers of sessile organisms, such as seaweeds or barnacles, can limit the impacts of physical weathering agents on shore platforms via a range of mechanisms. This includes moisture retention and buffering effects on temperature variations, which are thought to reduce the efficacy of substrate weathering via heatingcooling, wetting-drying, and salt crystallization (Stephenson and Kirk, 2000;Moura et al., 2012;Gowell et al., 2015). Through these mechanisms and their potential scale linkages that are thought to mediate wave-driven erosion (e.g., Naylor et al., 2012), bioprotection may reduce overall rates of platform lowering and retreat, particularly on shores and at tidal heights where biological covers are well developed. There is, however, no manipulative experimental field evidence for the direct impacts of biotic cover on platform weathering and erosion (but see Moura et al., 2012 for a comparative study). Evaluating such effects is especially relevant on relatively soft cohesive shore platforms (i.e., those formed by poorly consolidated sediments with a high silt and clay content), where erosion can proceed quickly relative to the time taken for the recovery of some sessile organisms (e.g., 5 years for mussels and erect macroalgae; Micheli et al., 2016) following mortality and/or dislodgement (Charman et al., 2007;Moses, 2014). Improving our understanding of biogeomorphological processes is also crucial to inform predictions of coastline responses to climate change, including shifts in the spatial arrangement of organisms and physical processes as a result of sea-level rise.
Eroding cliffs fronted by cohesive shore platforms are a pervasive feature along the Argentinean coastline (37-52 • S). At mid intertidal elevations, the cohesive platforms are dominated by bed-forming mussels (either Brachidontes rodriguezii, Perumytilus purpuratus, or both; see Arribas et al., 2013;Trovant et al., 2015). Mussels are known to physically modify environments (via physical ecosystem engineering; sensu Jones et al., 1994) in ways that affect the abiotic conditions and resources available to other organisms (Gutiérrez et al., 2003(Gutiérrez et al., , 2011. For instance, flow speeds, sediment transport, and temperature and moisture fluctuations are reduced in the interstices of mussel beds relative to adjacent bare rock (e.g., Carrington et al., 2008;Silliman et al., 2011;Bagur et al., 2016;Gutiérrez et al., 2019). While the ecological effects of physical ecosystem engineering by mussels are well understood (e.g., Borthagaray and Carranza, 2007;Silliman et al., 2011;Bagur et al., 2016;Gutiérrez et al., 2019), their potential to reduce weathering and erosion of the underlying rock remains largely unknown (see Baxter, 2015). Yet, mussel cover may protect the underlying rock from the mechanical impacts of waves and sediments as well as buffer variations in temperature and moisture, thus reducing the efficacy of weathering processes linked to rock expansion-contraction cycles and the crystallization of salts within rock pores.
Examining the effects of mussel cover on the weathering and erosion of cohesive shore platforms is particularly important in view of multiple, ongoing sources of mussel cover loss, including human harvesting and trampling as well as increased dislodgement/mortality due to the increased frequency, magnitude, and duration of storms and heatwaves (e.g., Tsuchiya, 1983;Erlandsson et al., 2006Erlandsson et al., , 2011Micheli et al., 2016;Mendez et al., 2017;Seuront et al., 2019). Here, we investigated mussel cover impacts on the weathering and erosion of intertidal cohesive platforms using removal experiments on the Argentinean coastline. Specifically, we evaluated the effects of mussel cover on water flow and abrasion, oscillations in rock temperature and moisture, and the accumulation of salts within the platform material. We predicted that (a) the hydraulic and abrasive impacts of sand-laden waves and currents, as well as the potential for thermal, wetting-drying, and salt-related weathering, would be higher on bare relative to the mussel-covered substrate, and (b) enhanced weathering and erosion following the removal of mussels would lead to substrate softening and platform lowering.

Study Site
This study was carried out on intertidal cohesive platforms at Copacabana (38 • 14 S, 57 • 46 W), 5 km northeast of Miramar, Buenos Aires Province, Argentina ( Figure 1A). These platforms are located at the base of an 8-12 m high cohesive cliff (known as Barranca Parodi in the stratigraphic literature; Figure 1B) developed in the Chapadmalal Formation (3.2-4.5 million years BP), which is primarily composed of reddish-brown, sandy siltstones originating from reworked loess deposits and sparse indurated caliche levels and nodules (Zárate and Fasano, 1989;Isla et al., 2015;Rico et al., 2020). The shore in this area faces southeast, which is the direction of the strongest swells (Fiore et al., 2009), and can be classified as "exposed" according to the MarLin wave exposure categories (i.e., an open coast facing away from prevailing winds but with a long fetch, and where strong Frontiers in Marine Science | www.frontiersin.org winds are frequent) 1 . Current velocities on the surface of the study platforms can be as high as 13 m s −1 (Gutiérrez et al., 2018). The tidal regime is semidiurnal and microtidal (0.80 m mean amplitude; Servicio de Hidrografía Naval, Argentina) 2 . Current coastal retreat rates in this area average 0.8 m yr −1 (Isla et al., 2018).
As is common for cohesive and soft rock platforms (e.g., Ledesma-Vázquez and Johnson, 1994;Gómez-Pujol et al., 2014;Moses, 2014), the studied platforms are cut by interconnected rills that run mostly perpendicular to the shoreline and are, thus, fragmented into elongated benches of varying size from a few square decimeters to a few square meters (see Figure 1C). The relatively small mussel, Brachidontes rodriguezii (up to 55 mm in length, most individuals <30 mm length), forms dense and primarily single-layered beds (200-800 ind. dm −2 ; Arribas et al., 2015;Gutiérrez et al., 2015) along the top and along the upper sides of the benches located within the mid intertidal zone ( Figure 1D). The height of these mid intertidal benches varies between 50 and 80 cm, though their top surfaces are relatively level with each other (i.e., height differences are mostly associated with variations in rill depth). Their lower sides and the bottoms of rills are periodically buried and abraded by sand and gravel (both sourced from nearby beaches; see Lamarchina et al., 2021) and, thus, lack biotic cover ( Figure 1D). Abrasion and undercutting of the benches is likely a significant mechanism of overall platform erosion independently of the biotic platform-lowering effects investigated here.

Field Experiment
Sixteen benches (0.4-0.8 m −2 top surface area) were selected on June 29, 2018, and mussels were removed from half while the remaining benches were left as mussel-covered controls. The mussel layer was pulled by hand or with the aid of a putty knife when strongly attached, but taking care not to scrape the bench surface. Treatments were assigned systematically to ensure adequate interspersion. Platform height, substrate hardness, flow/abrasion rates, near-surface temperatures, substrate moisture, and subsurface salt concentrations were subsequently measured over the course of the experiment (see section "Measurements of Response Variables"). The experiment was ended on November 21, 2018, following the incipient recruitment of barnacles, Balanus glandula, to the mussel-removal benches. While the size (<1 mm) and surface area covered by these barnacles (<5%) were small so as not to interfere with our measurements and sampling around this date, we predicted that barnacle coverage would become significant in the following weeks due to continuing recruitment and individual growth. This was confirmed in late December when the barnacle cover reached 20-30%.

Platform Height
Changes in the height of the mussel-removal and control benches were evaluated monthly for approximately 5 months until November 20, 2018. Variations in bench height were assessed as changes in surface level relative to benchmarks established within the platform. For this, we used a cost-effective, modified version of the micro-erosion meter (Stephenson and Finlayson, 2009), which we tailored to the specific needs of this research. Four nails (1.5 inches long) were driven into the substrate at each corner of a 15-cm 2 quadrat normal to the top of each bench, leaving ca. 1.5-cm section protruding from the surface (see Supplementary Figure 1A). At the time of each measurement, a square, removable PVC plate (sides = 15.5 cm) was mounted on top of the four nails (Supplementary Figure 1B) and the height of the substrate relative to the PVC plate (i.e., the benchmark level) was measured by inserting the depth measuring probe of a Vernier caliper through a tightly fitting hole made in the center of the plate (Supplementary Figure 1C). The PVC plate was removed after measurement so that only the nails remained in the field throughout the study period. To ensure that the PVC plate was mounted in the same position each time and relative height data were obtained at the same points on the bench surfaces, circular notches were made in each corner of the plate into which the nail heads could be fitted. The nail heads themselves could not be used as benchmarks (i.e., as erosion pins; e.g., Turowski and Cook, 2017) as surface cracking and pitting occurred within a ca. 0.5-cm radius around the nails. It was also not possible to use screw bolts instead of nails as the substrate easily disaggregated when attempting to secure them.

Substrate Hardness
The surface hardness of the benches was measured immediately after the final height measurements (on November 21, 2018) as this involved destructive sampling. Specifically, a substrate chip (∼300-400 cm 3 ) was extracted from each experimental bench with a hammer and chisel. To minimize any initial differences in moisture content, which may have affected the hardness measurements, the samples were sun-dried outdoors for a total of 10 h over the two subsequent days (in dry, sunny conditions) and stored indoors overnight to avoid moisture uptake from dew and mist. The samples were not oven-dried to constant weight as they readily disaggregate, which would prevent the subsequent hardness measurements.
A Shore A durometer was subsequently used to measure the hardness of the upper (i.e., exposed) face of the dry samples. The Shore A durometer is designed to measure resistance to indentation in flexible-to-hard rubber and semi-rigid plastics but has also been used to measure the surface hardness and, indirectly, cohesion of other low-strength consolidated materials including cement-stabilized clay, plasters, and limebased mortars (e.g., Wang et al., 2015;Santos et al., 2018Santos et al., , 2020. The A-type durometer has a small pin at the end (a 35 • truncated cone with a 0.79-mm diameter tip and a 1.40-mm base) that indents the test material under a fixed-force spring. Hardness values are measured in "Shore A Units" (ranging from 0 to 1 indicating low and high indentation strength, respectively), read from a digital LCD. Although several alternative methods for measuring surface hardness in situ and non-destructively are increasingly applied in geomorphological research, including on rock platforms (e.g., Knight and Burningham, 2020), we used destructive sampling and the durometer both for simplicity and to control for the potential impact of moisture content on substrate hardness and surface roughness, which are known to influence other measurement devices (e.g., Viles et al., 2011;Wilhelm et al., 2016;Desarnaud et al., 2019).

Flow/Abrasion Rates
The effects of mussel cover on water flow and abrasion at the surface of the benches were evaluated on August 27-29, 2018. For this, oven-dried (48 h at 60 • C) and pre-weighed plaster discs (5 cm in diameter and 1.5 cm in height, approximating the thickness of the mussel bed) were nailed to the musselremoval and control benches (one disc per experimental bench, n = 8; Supplementary Figure 2). These discs were prepared using three parts Plaster of Paris and four parts water and were molded using discoidal casts. The discs were installed at least 40 cm from the points where the height measurements were made to minimize any possible interference. On the musselcovered benches, the discs were attached after a small number of mussels were removed, so that they could be secured directly onto the substrate with their sides shielded by the surrounding mussels (see Supplementary Figure 2). After 48 h of exposure (i.e., four semidiurnal tides), the discs were removed, taken to the laboratory, oven-dried, and re-weighed. Plaster weight loss was used as a surrogate measure of relative flow rates/turbulence and abrasion around and within the mussel beds and cleared areas (e.g., Guichard and Bourget, 1998;Baxter, 2015).

Substrate Surface Temperature
Waterproof temperature loggers (Hobo Pendant MX2201) were installed on the mussel-removal and control benches (two benches per treatment, with a single logger per bench). As these loggers are small and flat (3 cm in diameter and 2 cm in thickness), they allowed us to measure air temperature at the surface of the substrate, and between the rock and the overlying mussel bed. Holes were drilled in the substrate into which the loggers fitted, so that their tops were flush with the surface (i.e., the holes were the same width and depth as the loggers). The loggers were then nailed in place, within their mounting hole. For the mussel-covered control surfaces, a square piece of the mussel layer (sides = 15 cm) was first carefully removed with the aid of a putty knife (Supplementary Figure 3A). A hole was then made into which a logger was fitted and fastened (Supplementary Figures 3B,C), and the square piece of mussel bed was then replaced in its original position (Supplementary Figure 3D). To our knowledge, this is the first time that temperature data have been obtained on a shore platform beneath an existing cover of sessile invertebrates, similar to previous studies on seaweed canopies (Coombes et al., 2013).
Data were recorded on an exceptionally hot late-spring day on November 21, 2018, when the sky was clear, the platform was air-exposed at noon (low tide slack occurred at 11:50 PM), and temperatures reached 31 • C as recorded by a nearby weather station (Mar del Plata Airport, 37 • 55 S, 57 • 34 W, ca. 40 km from the study site). On this day, temperature data were obtained during the low-tide air exposure period between 11:10 AM and 2:05 PM, when the platform began to be flooded again by the rising tide. Therefore, the data obtained reflect the influence of mussel cover on rock-surface temperatures on a typical hot summer's day in this area (mean maximum daily temperature in January = 25 • C; Giampietri and Piccolo, 2000).

Substrate Moisture Content
Two rock samples per experimental bench were obtained on November 21 in addition to those used for the hardness analysis (see section "Substrate Hardness") using a 2-cm diameter steel core driven to a depth of 1.5 cm. The day was clear and hot, and the platform was air-exposed at noon; the samples were taken between 1:00 PM and 2:00 PM, more than two and a half hours after the platform was first exposed to the air. All of the samples were sliced into three 0.5-cm depth sections with the aid of a putty knife. One sample per bench (i.e., one set of three 0.5-cm sections) was analyzed for salt content (see section "Subsurface Salt Content"), and moisture content was determined using the remaining samples, calculated as the difference between their wet and dry weights. All wet-weight measurements were made in the field immediately after the samples were collected using a portable scale (precision = 0.001 g). Their dry weights were subsequently determined in the laboratory after oven-drying (60 • C) to constant weight.

Subsurface Salt Content
The dried sample sections from each bench (at three depths of 0-0.5, 0.5-1.0, and 1.0-1.5 cm) were ground to dust with a glass mortar and pestle and subsequently analyzed for chloride and sulfate ions using ion chromatography (Dionex IC DX500; see Coombes et al., 2017).

Data Analyses
Repeated measures ANOVA (Neter et al., 1990) was used to evaluate the effect of Treatment (i.e., mussel removal; fixed factor) and Time (i.e., sampling date; random repeated measures variable) on platform height loss. Beta regression (i.e., a technique to model continuous random variables that assume values in the 0-1 interval; see Cribari-Neto and Zeileis, 2010) was used to evaluate the effect of Treatment on flow/abrasion rates (proportional weight losses from the plaster discs) and substrate hardness (Shore A units ranging between 0 and 1), as well as the effect of Treatment and Depth (fixed variable) on substrate moisture (proportional weight loss after oven drying) and salt content (salt weight/sample weight). All analyses were run in R (R Core Team, 2020) using the "lme4" (Bates et al., 2015) and "betareg" (Cribari-Neto and Zeileis, 2010) packages for the repeated measures ANOVA and beta regression, respectively.

RESULTS AND DISCUSSION
Our findings illustrate that both weathering regimes and wave-driven processes can be attenuated by mussel cover on intertidal cohesive shore platforms. Specifically, we found that (a) mussel cover buffers multiple subaerial weathering agents, i.e., temperature and moisture variations, and salt crystallization, and (b) mussel cover modulates the hydrodynamic/abrasive forces acting on the platform surface that contribute to platform erosion via their own mechanical impacts and/or by removing previously weathered material. We have also observed that (c) mussel removal from these platforms is followed by decreases in surface hardness and platform height. In the following sections, we discuss these mechanisms in detail and evaluate their roles in the context of local and temporal factors that may mediate their occurrence and efficacy both on our study platform and elsewhere (e.g., substrate type, wave climate, and seasonal influences on temperature and desiccation).

Effects of Mussels on Weathering Regimes
Mussel cover efficiently buffered high temperatures during air exposure; near-surface temperatures on the mussel-covered benches during a period of exposure on a hot spring day (11:10 AM to 2:05 PM) were lower than on the cleared bench surfaces (Figure 2). This is in line with previous studies showing the dampening effects of algae or sessile invertebrate covers on the temperature of intertidal consolidated substrates during daytime air exposure (Coombes et al., 2013(Coombes et al., , 2017Baxter, 2015;Gowell et al., 2015;Scrosati and Ellrich, 2018). The heat-attenuating effects of mussels were striking, with reductions ranging between 7 and 12 • C (Figure 2). This is attributed to substrate shading and, possibly, to increased evaporative cooling due to the water retained in the interstitial spaces of the mussel beds as well as the mussel bodies themselves (see Helmuth, 1998). Evaporative cooling mechanisms have also been proposed to explain lower substrate temperatures in the presence of seaweeds and barnacles (Gowell et al., 2015;Coombes et al., 2017).
As well as general reductions in temperature, mussel cover was found to reduce short-term, seconds-to-minute fluctuations in temperatures during air exposure, indicated by the more jagged temperature curves for the mussel-removed surfaces in Figure 2. By modulating these rapid temperature fluctuations (attributable to passing cloud and wind gusts) as well as the amplitude of daily thermal cycles (associated with tidal regimes FIGURE 2 | Air temperatures at the surface of mussel-covered and mussel-removal (exposed) benches during a 3-h air-exposure period at noon on an exceptionally hot late-spring day (November 21, 2018). Each curve plots the data registered every 15 s by temperature loggers installed in two mussel-covered and two exposed benches (one logger per bench). The loggers attached to the mussel-covered benches were placed underneath the mussel layer. and day-night cycles, e.g., Coombes, 2011;Pappalardo and D'Olivo, 2019), mussels are expected to reduce the frequency, magnitude, and rate of substrate expansion/contraction. While not likely sufficient to initiate breakdown via thermal "shock" effects (see Hall and Thorn, 2014), the buffering of such repeated, low-magnitude thermal changes is thought to limit the potential for thermal "fatigue" in a range of environments including on rocky coasts Coombes et al., 2013;Hall and Thorn, 2014;Molaro et al., 2020). From the data in Figure 2, it can also be inferred that sudden temperature decreases occur when air-exposed substrates heated by the sun are rapidly cooled by the incoming tidal waters (the typical water temperature in this area and time of the year is 16 • C) and that mussel layers can also modulate these changes. It remains to be evaluated whether buffering of such higher-magnitude and rapid thermal changes by mussels and other forms of biota alters the potential for thermal shock-related breakdown on coastal platforms (Coombes, 2011(Coombes, , 2014. The dampening effects of mussel cover on high near-substrate temperatures were associated with concurrent decreases in substrate desiccation and salt crystallization relative to cleared areas. The moisture content of the mussel-covered substrates was significantly higher than that of the cleared substrates after a few hours of exposure to hot air temperatures at low tide, and this difference was consistent across the three 5-mm depth ranges of the substrate samples (Table 1 and Figure 3). Asterisk indicates significant effects (p < 0.05).
FIGURE 3 | Substrate moisture content (water mass as a percentage of sample mass) of mussel-covered and mussel-removal (exposed) bench materials at three depth ranges. Substrate samples were collected on a clear and exceptionally hot late-spring day (November 21, 2018) 2.5-3 h after the tide had receded and the platform was air-exposed. The vertical lines above the bars indicate 1 SD. Asterisks indicate significant treatment effects at the corresponding depth level (p < 0.05). This indicates that mussel covers slow substrate desiccation rates and reduce net substrate desiccation, at least during hot periods of tidal exposure. This has implications for the efficacy of wetting and drying as a weathering mechanism across the intertidal and splash zones in association with the repeated expansion and contraction and microtopographic fluctuation of the substrate (e.g., Kanyaya and Trenhaile, 2005;Trenhaile, 2006;Porter and Trenhaile, 2007). Indeed, the modulation of wetting and drying cycles (and associated expansion and contraction) by biota, including microorganisms (e.g., Mayaud et al., 2014), may have significant weathering implications for a range of substrates, especially clay-bearing types including mudstones (e.g., Stephenson and Kirk, 2000;Baxter, 2015) and the sandy siltstones of our study platform composed of 10-20% coarse-to-medium silts and 60-70% fine silt and clay (Rico et al., 2020). The subsurface salt content of the platform material further indicates the moisture-modulating role of mussels. Specifically, we found higher concentrations of both chloride and sulfate ions in the upper 5-mm layer of the substrate from the mussel-removal benches relative to mussel-covered controls, while there was no difference in the deeper layers (5-15 mm) (Table 2 and Figure 4). This indicates that salt crystallization is reduced in the uppermost layer of the platform when covered by mussels-probably in association with their temperature-buffering and moisture-retaining functions. While such effects have been proposed based on laboratory simulations (see Gowell et al., 2015), this provides some of the first direct field evidence that biotic cover can reduce the efficacy of salt crystallization in the surface layers of intertidal platforms (also see Baxter, 2015 andCoombes et al., 2017). Salt weathering is thought to be important on rocky shores during hot weather, when desiccation during daytime tidal exposure causes the crystallization of dissolved salts in confined rock pores and associated destructive forces (Mottershead, 1989(Mottershead, , 2013Stephenson and Kirk, 2000;Porter and Trenhaile, 2007;Coombes, 2014). As well as reduced amounts of salt overall, lower crystallization frequencies and the less frequent hydration and thermal expansion/contraction of salt crystals within the upper layers of mussel-covered substrates due to the observed moisture/temperature-buffering effects may limit weathering relative to exposed substrates (e.g., Gowell et al., 2015;Coombes et al., 2017).

Effects of Mussels on Hydrodynamic and Abrasive Forces
Based on the significantly greater weight losses of the plaster discs placed on mussel-covered benches relative to those without mussels (χ 2 = 12.81, df = 1, p < 0.001; Figure 5), our results indicate that hydrodynamic and abrasive forces are lower within mussel beds. This conforms to the findings of Carrington et al. (2008), who observed that water flow in the interstitial spaces of mussel beds is low relative to above mussel beds and in mussel-free estuarine bottoms. This provides some initial support for mussel beds reducing the mechanical impacts of sedimentladen waves and tidal flows on platform rock-especially very soft rock as studied here-both by physically armoring the surface and moderating near-surface water flows. Similar findings were reported by Baxter (2015) for mussel-colonized and artificially cleared patches on an intertidal mudstone platform in Wales, United Kingdom. Further work is now needed to determine whether such effects translate into altered rates of erosion FIGURE 4 | Subsurface salt content of the mussel-covered and mussel-removal (exposed) bench substrates at three depth ranges. Substrate samples were collected on a clear and exceptionally hot late-spring day (November 21, 2018) 2.5-3 h after the tide had receded and the platform was air-exposed. The vertical lines above the bars indicate 1 SD. Asterisks indicate significant treatment effects at the corresponding depth level (p < 0.05).
FIGURE 5 | Flow/abrasion rates at the surfaces of mussel-covered and mussel-removal (exposed) benches based on the weight losses of plaster discs (5 cm diameter, 1.5 cm height). Discs attached to the mussel-covered benched had their sides shielded by the surrounding mussels. The vertical lines above the bars indicate 1 SD. Differences among treatments were significant (p < 0.001). and direct loss of material from the surfaces of colonized intertidal rock.

Effects of Mussels on Platform Height and Substrate Hardness
Mussel removal led to significant decreases in platform height ( Table 3). Height losses in the mussel-removal treatments averaged 2 mm after 6 months, while the height of musselcovered benches remained relatively constant over the same period, with minor variations attributable to measurement error as well as possible swelling and shrinking in response to wetting-drying effects (Figure 6). The observed decreases in platform height following mussel removal cannot directly be ascribed to the weathering processes discussed in section "Effects of Mussels on Weathering Regimes" as the height losses mostly occurred during the first 2 months of the experiment and under winter conditions (see Figure 6). Compared to other times of the year, wetting-drying and salt weathering are expected to be less efficient during winter, when desiccation is limited due to relatively low temperatures and irradiation. On the other hand, thermal fatigue due to short-term temperature variations (i.e., second-to-minute fluctuations attributable to winds and clouds; see section "Effects of Mussels on Weathering Regimes") might continue to operate in winter and may, therefore, have contributed to FIGURE 6 | Variations in the height of mussel-covered and mussel-removal (exposed) benches during the study period. The experiment began on June 27, 2018, and the differences among treatments were significant on all subsequent sampling dates (p < 0.05). The vertical lines above and below the data points indicate 1 SD.
the observed lowering of the mussel-removal benches. The likelihood that frost weathering occurred during these months is also low as minimum air temperatures fell below the freezing point of seawater (−2.5 • C) on only two occasions, when the platform was exposed to such temperatures for less than 2 h (see Table 4). This suggests that the duration of freezing periods and the number of freeze-thaw cycles were well below those required for the weathering of consolidated intertidal substrates (see Robinson and Jerwood, 1987;Dewez et al., 2015). Nevertheless, the potential impacts of biotic cover on the ice-driven breakdown of coastal rocks warrants further investigation, as freezing processes can be important in some situations (e.g., Trenhaile and Rudakas, 1981;Trenhaile and Mercan, 1984;Robinson and Jerwood, 1987;Dewez et al., 2015). Mussels may be especially suited to protect intertidal rocks from frost weathering as they have evolved local adaptations to survive freezing temperatures (e.g., down to −13 • C in the case of Arctic populations of Mytilus edulis; Thyrring et al., 2020). The hydraulic and abrasive impacts of sediment-laden waves are also unlikely to be the cause of the measured lowering of the mussel-removal benches during the first 2 months of the experiment. Although storm conditions conducive to erosion in the study area tend to be more frequent during the winter (i.e., the early phases of our experiment), they occur during all seasons (see Fiore et al., 2009). Within the specific timeframe of this experiment, we found no evidence of increased storm frequency or intensity during the earlier, winter months. Indeed, based on an analysis of positive storm surges (i.e., positive deviations of the observed tide relative to the astronomic tide that exceeded 30 cm and persisted for 6 h or more; see Fiore et al., 2009), storms and wave energy were evenly distributed across the experimental period (see Figure 7) and cannot, therefore, directly explain the temporal pattern of height loss on the musselremoval benches. 4 | Dates with air temperatures below the freezing point of seawater (i.e., −2.5 • C) and associated platform exposure times quantified as the number of hourly air-exposed temperature records <−2.5 • C.

Date
Minimum temperature ( • C) Timing of temperatures <−2.5 • C Air exposure period Platform exposure events at <−2.5 • C (h) 10-July-2018 −2.6 8:00 am 0:00 am-3:00 pm 1 25-July-2018 −3.2 7:00 am, 8:00 am 0:00 am-3:00 pm 2 Hourly temperature data from Mar del Plata Airport (37 • 55 S, 57 • 34 W, ca. 40 km from the study site) were obtained from https://www.smn.gob.ar/descarga-de-datos. The timing of the platform air-exposure period was estimated from hourly tidal measurements from the Mar del Plata Port (ca. 35 km from the study site) provided by the Servicio de Hidrografía Naval, Ministerio de Defensa, Argentina, and assuming a tidal elevation of 0.84 m as the threshold for platform flooding based on an empirical logistic model (see Supplementary Material). The long platform exposure periods on both dates are attributable to tidal forcing by offshore winds (see Fiore et al., 2009).
One possible explanation is that mussels weaken the rock they are attached to, for example by enhancing chemical weathering via the retention of water and organic matter (see Moura et al., 2012), but continue to keep these surfaces consolidated as long as they remain in place. Once detached or removed, lowering rates would be relatively rapid at first and then stabilize. This process is analogous to the rapid, catastrophic loss of rock and stone in the terrestrial environment following the removal of protective weathering crusts, which exposes the softened, highly erodible material underneath (e.g., Smith et al., 2010). While it is difficult to confirm such effects, it is possible that the impacts of mussels and other forms of biotic cover on intertidal rock weathering may involve a complex balance of bioprotective (when in place) and bioerosive (once detached/ removed) effects.
Another possible explanation is that the experimental removal of the mussels weakened the surfaces of the benches. Mussel detachment involves the failure of the collagenous byssal threads they produce to adhere to the substrate. Typically, byssus failure occurs either because the thread itself or its basal portion breaks (thread and root failure, respectively), the adhesive plaque at the base of the thread detaches from the substrate (plaque failure), or the substrate to which the adhesive plaque attaches is dislodged (substrate failure) (see Bell and Gosline, 1997). Individually, a substrate failure implies the removal of a small substrate chip from the platform surface (i.e., at scales of tens to hundreds of micrometers, based on adhesive plaque size). Yet, a high number of substrate failures per unit area can occur during the mass detachment of mussels. For instance, mussels 15-20 mm in length at our study area have approximately 80 byssal threads on average, and 30% of the failures observed during their detachment occur at the substrate level (Gutiérrez et al., 2018). Considering that mussels within this size range can occur at densities of 200 ind. dm −2 at this site (Soria, 2020), their mass detachment (either due to experimental removal or storm wave action) could cause approximately 4,800 substrate failures per dm −2 . Although an underestimate is based on a subset of mussel size classes occurring in the field, this calculation illustrates that these highly friable bench surfaces could be densely chipped and weakened during mussel detachment.
While most of the observed decrease in platform height occurred within the first 2 months after mussel removal, with a slower rate of lowering thereafter, it is not possible to rule out the potential for a subsequent period of rapid lowering. Unfortunately, barnacle recruitment prevented us from continuing the experiment over the summer (see section "Field Experiment") under higher temperature, irradiation, and desiccation conditions. However, our substrate hardness measurements, which were carried out at the end of the experiment in late spring, indicated significantly lower surface hardnesses in the mussel-removal patches compared to the mussel-covered controls (χ 2 = 12.36, df = 1, p < 0.001; Figure 8). This suggests that progressive surface weathering during the 6 months following mussel removal was linked to the development of weaker bench surfaces relative to those with mussels still attached.
FIGURE 7 | Positive deviations in the tide relative to tidal predictions during the experimental period. The black bars above the curve indicate storms, defined as periods with deviations >30 cm (gray dashed line) lasting for more than 6 h (see Fiore et al., 2009). Sampling dates are indicated as triangles on the x-axis. Tidal records were obtained from the tide gauge at Puerto Mar del Plata (38 • 02 08 S, 57 • 31 52 W, ∼ 35 km north of the study site). Tidal records and predictions were provided by the Servicio de Hidrografía Naval, Ministerio de Defensa, Argentina.
FIGURE 8 | Substrate hardness of the mussel-covered and mussel-removal (exposed) benches. Hardness was measured in Shore A units, which range from 0 (low) to (1) high. The vertical lines above the bars indicate 1 SD. Differences among treatments were significant (p < 0.001).

CONCLUSION AND PROSPECTS
This study shows that mussel cover has clear bioprotective potential by (1) modulating temperature, desiccation, and salt crystallization regimes; and (2) reducing the impacts of water flow and abrasion, which is in line with previous findings for seaweeds and barnacles (Coombes et al., 2013(Coombes et al., , 2017Gowell et al., 2015). Additionally, we illustrate for the first time that biotic cover can limit the erosion of shore platforms based on observed rates of surface lowering. Aside from the possible mechanisms of platform lowering following mussel removal (i.e., loss of protective mussel cover and/or exposure of relatively weaker/softer surfaces, weakening of the platform surface during mussel detachment, or a combination of these; see section "Effects of Mussels on Platform Height and Substrate Hardness"), our observations highlight the importance of persistent biotic covers for the stability of friable and soft platform materials. Indeed, in the case of semi-consolidated shore platforms, which can undergo weathering and erosion over relatively short timescales, this implies that the loss and recovery dynamics of ecological communities can have particularly important geomorphological implications. Studies involving other forms of biotic cover and different physical settings are now needed to explore these linkages.
In addition to their stabilizing effects on cohesive shore platforms, mussels may also have indirect stabilizing impacts on the cliffs they front. Shore platforms are a by-product of cliff retreat but also act as a natural defense against wave attack (Balson et al., 2006;Brew, 2007;Moses, 2014;Payo et al., 2015). Consequently, the protective role of platforms gradually decreases as they are progressively lowered by weathering and erosion (Trenhaile, 2002). It follows that, by stabilizing the studied platforms, mussels, which typically form beds at midintertidal elevations and persist all year round, help maintain the protective function of these platforms with respect to wave attack and cliff erosion.
Threats to the persistence and extent of mussel cover on intertidal platforms are varied and increasing in frequency and magnitude globally (e.g., human harvesting and trampling, storms, pollution, climate change, and heat waves; see Tsuchiya, 1983;Erlandsson et al., 2006Erlandsson et al., , 2011Micheli et al., 2016;Mendez et al., 2017;Seuront et al., 2019). In the light of our findings, these changes are expected to accelerate erosion rates on cohesive shore platforms along with the loss of potential protective functions for backing cliffs. Mussel bed conservation and enhancement measures could, therefore, contribute to coastline stability alongside the maintenance of associated biodiversity in the face of a range of anthropogenic threats including climate change and sea-level rise.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
JLG, MC, MP, and JAG contributed to the conception and design of the study. JAG, JLG, and SS conducted the field research. JAG and MC analyzed substrate samples in the lab. JAG and JLG analyzed the data. JLG, JAG, MC, and FI wrote the original draft. All authors contributed to manuscript revision, read, and approved the submitted version.

FUNDING
This study was supported by the National Geographic Society (Grant No. CP-120R-17, to JLG, MC and MP) and Agencia Nacional de Promoción Científica y Tecnológica (PICT-2019-03904 to all authors).