Assessment of Exploitation Intensity of Commercial Species and Associated Benthic Communities, in Chilean Marine Management Areas of North Patagonia

The Aysén region of Chile (North Patagonia), has had limited studies on the effectiveness of management and exploitation areas of benthic resources, and performance relative to open access areas in this region has never been evaluated. We evaluated seven management areas (MAs) and five open access areas (OAAs) between 43.9°S and 45.2°S for exploitation intensity of three commercial species (Concholepas concholepas, Loxechinus albus, and Ameghinomya antiqua) together with characterization of the benthic community. Indicators based on size, density and weight were used to evaluate exploitation intensity of commercial species. Associated benthic communities were evaluated considering density, species composition, and community structure. We found a high species richness and a community structure with low variability between MAs and OAAs. Low densities and small sizes classes of C. concholepas in most of the areas indicated high exploitation intensity in both MAs and OAAs. In this context, a permanent ban to harvest C. concholepas within OAAs may need to be reevaluated since with no enforcement and monitoring, the exploitation status of this species remains unclear in these areas. L. albus in most areas were absent in the harvestable sizes, which could be indicating high exploitation intensity in both regimes. High densities and small sizes of C. concholepas and L. albus in some MAs, indicated a potential recruitment zones which bears further investigation. A. antiqua, showed better conditions than other commercial species evaluated, with no significant differences in densities and size-based indicators when comparing OAAs and MAs. Benthic communities were dominated numerically by the Echinoidea class in both MAs and OAAs, with L. albus, Arbacia dufresnii and Pseudechinus magellanicus being the dominant species. High densities of sea urchins co-occurring with low coverage of macroalgae found in MAs-Gala could indicate that a sea urchin barren was dominant during the study period. On the other hand, high densities of Cosmasterias lurida, a predatory sea star, in conjunction with low densities of C. concholepas in most of the studied areas suggested that a shift in predator roles is occurring. No differences were estimated in terms of fisheries indicators and benthic community structure across the two management regimes, suggesting the poor performance of MAs in Aysén region. Our data also support the need to improve monitoring of MAs especially with respect to associated benthic community incorporating a broader spatial scale.

The Aysén region of Chile (North Patagonia), has had limited studies on the effectiveness of management and exploitation areas of benthic resources, and performance relative to open access areas in this region has never been evaluated. We evaluated seven management areas (MAs) and five open access areas (OAAs) between 43.9 • S and 45.2 • S for exploitation intensity of three commercial species (Concholepas concholepas, Loxechinus albus, and Ameghinomya antiqua) together with characterization of the benthic community. Indicators based on size, density and weight were used to evaluate exploitation intensity of commercial species. Associated benthic communities were evaluated considering density, species composition, and community structure. We found a high species richness and a community structure with low variability between MAs and OAAs. Low densities and small sizes classes of C. concholepas in most of the areas indicated high exploitation intensity in both MAs and OAAs. In this context, a permanent ban to harvest C. concholepas within OAAs may need to be reevaluated since with no enforcement and monitoring, the exploitation status of this species remains unclear in these areas. L. albus in most areas were absent in the harvestable sizes, which could be indicating high exploitation intensity in both regimes. High densities and small sizes of C. concholepas and L. albus in some MAs, indicated a potential recruitment zones which bears further investigation. A. antiqua, showed better conditions than other commercial species evaluated, with no significant differences in densities and size-based indicators when comparing OAAs and MAs. Benthic communities were dominated numerically by the Echinoidea class in both MAs and OAAs, with L. albus, Arbacia dufresnii and Pseudechinus magellanicus being the dominant species. High densities of sea urchins co-occurring with low coverage of macroalgae found in MAs-Gala could indicate that a sea urchin barren was dominant during the study period. On the other hand, high densities of Cosmasterias lurida, a predatory sea star, in conjunction with low densities of C. concholepas in most of the studied areas suggested that a shift in predator roles is occurring. No differences were INTRODUCTION Management and exploitation areas of benthic resources (MAs) were established along the Chilean coast to contribute to the ecological sustainability of benthic resources; their establishment was driven by the biological collapse in the 80s of the muricid mollusk Concholepas concholepas, a previously profitable benthic fishery (Castilla, 1994;Stotz, 1997;Gonzalez et al., 2006). This administrative system gives territorial user rights (TURFs) to artisanal fisheries organizations that are legally constituted for the co-management (with the state) of benthic resources within a limited geographic area (Zuñiga et al., 2008;Gelcich et al., 2010). After two decades of MAs functioning (Stotz, 1997;Gonzalez et al., 2006;Gelcich et al., 2016), the ecological effects of this management approach have been somewhat diverse. In Central Chile, higher densities, larger sizes, and lower mortality rates of commercial species (C. concholepas, Loxechinus albus, and Fissurella spp.) were found in MAs compared to Open-Access Areas (OAAs, i.e., no access restriction) (Castilla and Fernandez, 1998;Gelcich et al., 2010;Defeo et al., 2014;Andreu-Cazenave et al., 2017). Positive effects were also reported in the associated macroinvertebrate communities where higher species richness, biomass and densities were observed within MAs compared to OAAs (Gelcich et al., 2008, 2012, Biggs et al., 2016. In contrast, performances of MAs in northern Chile were classified as poorly sustainable, with few stocks showing stability in densities, sizes, and catches over time (Arias and Stotz, 2020). Moreover, to maintain productivity, studies showed that at least half of the catch was extracted from areas outside MAs (Gonzalez et al., 2006). The performance of MAs has been heterogeneous and appears to betoo complex to generalize for a country such as Chile, with an extensive coast that spans ∼38 • of latitude, where a variety of environmental, social, geographical, and economical factors combine to produce different outcomes. In this context local research is needed to integrate site-specific particularity into the global analyses of benthic management.
The Aysén region of Chile (43 • 38 -49 • 16 S) harbors many small coves (fishers' villages) throughout islands, fjords and channels in which the subsistence primarily relies on artisanal fisheries. Their fishing grounds are situated within a marine ecosystem considered not only pristine compared to northern regions (Godoy et al., 2010;Navarrete et al., 2010), but also characterized by a high species richness (Fernandez et al., 2000;Häussermann and Försterra, 2009;Försterra et al., 2016;Betti et al., 2017;Bertolino et al., 2020). The presence of diverse and unique habitats not only support species richness but also host high abundances and biomass of important commercial species (Flores et al., 2020;Pardo et al., 2020). These highly productive ecosystems have been the target of illegal fishing, and regional fisheries authorities have needed to co-operate with local fishers to address the constant influx of fishers from neighboring regions. In this context, the introduction of MAs provided the potential for exclusive privileges in accessing benthic resources within a limited geographic area, and the interest and cooperation of local fishers has therefore increased in recent years. Currently in Aysén, there are 75 established MAs (data obtained from the Chile's Undersecretariat of Fisheries and aquaculture, Subpesca, 2020a), with 56% (43 total MAs) having resource management plans approved and 18 of which have extracted the estimated catch quota over the last years (2016-2019). Prior to this (2010)(2011)(2012)(2013)(2014)(2015), only 7 MAs contributed to landings, showing that there is an increased interest from fishers to harvest under this management system. Moreover, there are now 210 MAs pending approval in Aysén. The remoteness and intricate geography of Aysén has limited studies on the effectiveness of MAs (Moreno and Revenga, 2014;Romero et al., 2019), and no studies have yet evaluated their performance relative to OAAs in this region. On the other hand, follow up studies that considered an annual stock assessment in the management area provide data on catch quotas for commercial species but neither explicitly evaluated the intensity of exploitation, nor assessed the ecological effects on the associated benthic community. Even though, several studies have indicated the potential degrading impact of harvests upon habitats and biodiversity (Perez-Matus et al., 2017;De Juan et al., 2018;Contreras et al., 2019), and therefore on resilience of communities and populations in continuing to provide an ecosystem service. The responses of communities to harvest can be highly variable and complex since many factors are involved, on one hand there is the magnitude, timing, and areal extent of fishing (often can be well documented), on the other, the productivity of the ecosystems and the biological, and ecological characteristics of harvestable species (often with a poor level of understanding) which determine different effects on ecosystem processes and functions (Levin et al., 2009;De Juan et al., 2015). The effects of harvest on the benthic community structure can be both direct and indirect (Jennings and Kaiser, 1998;Pinnegar et al., 2000). Decrease in the abundances and changes in size structures of commercial species had been described as direct effects while some examples of indirect effects are benthic habitat degradation and changes on trophic interactions (Pinnegar et al., 2000). The monitoring of these cascade-type effects, such as a decrease in macroalgae biomass due to grazing effects by strict herbivores (e.g., L. albus) (Dayton, 1985;Buschmann et al., 2004;Wright et al., 2005;Contreras et al., 2019) or the detection of changes in the roles of top predators (e.g., C. concholepas) and their concomitant effect on benthic communities structure (Gaymer and Himmelman, 2008;Navarrete et al., 2010), are crucial to increase the knowledge related to responses of the ecosystem associated to impact of fishing.
The main objective of the present study was to evaluate MAs in the Aysén region by comparing their performance with OAAs in relation to: (i) exploitation intensity of three commercial species (Concholepas concholepas, Loxechinus albus, and Ameghinomya antiqua) and (ii) characterization of the overall benthic community structure associated with these commercial species. These results are discussed in the context of management of fisheries and benthic communities in Patagonian coastal waters.

Study Area
The study was carried out in Northern Patagonia, in an area located between 43.9-45.2 • S and 73.5-73.7 • W (Figure 1), situated between southern Corcovado Gulf (Guaitecas Islands) and southern Moraleda Channel (Meninea constriction, 45.2 • S-73.6 • W). The Moraleda Channel is 90 miles long, has an approximate average depth of 250 m, and is the main channel that separates the archipelago to the west from the continental coast to the east.
In the study area, seven MAs and five OAAs were sampled; two MAs (GUA1 and GUA2) localized northwest of Las Guaitecas archipelago near Melinka Cove (∼43.9 • S) with more oceanic influence (i.e., SAAW) and with greater exposure to wave action compared to the other studied areas (Sievers and Silva, 2008). The main substrates (>80%) of these areas were hard bottom (i.e., boulder and rock slabs) in GUA1 and mixed bottom (hard and soft bottom) in GUA2 (i.e., sand, shells cover, boulder, and rock, Table 1). The other areas were localized in the so-called "inland sea" in the central-south zone of the Moraleda Channel where there is a greater influence of fresh waters (main water masses are represented by EW and MSAAW). In the central study zone, two MAs (GALA1 and GALA2) were located outside the Jacaf Channel near Puerto Gala Cove (∼44.2 • S). Jacaf is a narrow and deep channel (>400 m) that connects with the Puyuhuapi Fjord, from which it receives freshwater influence (Schneider et al., 2014). The main substrate in these areas was hard bottom (rock slabs and rock) in Gala2 and soft bottom (sand) in Gala1. In the southern zone (Huichas), three MAs (HUIC1, HUIC2, and HUIC3) and five OAAs (OAAs1 to OAAs5) were located mainly in small, protected bays along the Moraleda Channel, close to Puerto Aguirre Cove (∼45.2 • S). The main substrate in these MAs (HUIC1, 2, 3) was hard bottom (rock, rock slabs, and boulder). In most of OAAs (OAA3, 4, and 5), substrate was dominated by mixed (sand, gravel, rock slabs, rock, and boulder), being soft bottom (sand) the main substrate only in OAA2, and hard bottom in OAA1. In MAs, substrate type was obtained from baseline studies that considered characterization of the total surface of each area ( Table 1). In the case of OAAs area, substrate type was characterized and quantified in each area. MAs were grouped considering their closeness to their coves since each of them has its socio-economic particularities that could influence their management, in this sense three zones were established: (1) Melinka (MAs-Guaitecas), (2) Gala (MAs-Gala), and (3) Aguirre (MAs-Huichas).
Open access areas were selected according to historical fishing grounds declared by local fishers of commercial species evaluated in the present study (C. concholepas, L. albus, and A. antiqua). In these areas, fishers who have registered commercial resources in the National Fishery Service can harvest them with the only restriction of following national regulations (e.g., minimum legal size of extraction, reproductive closure). Although only OAAs were studied around Huichas due to logistic problems, it is expected that these OAAs would allow for comparison among management regimes since these areas have a similar suitable substrate for the commercial species evaluated in MAs and its associated benthic communities.

Sampling
Commercial species were selected in MAs according to baseline studies (Subpesca, 2016) that described the species that constitutes the object of the management and exploitation plan according to their direct quantification (Subpesca, 1995). In this context, three species were considered; the red sea urchin Loxechinus albus (Molina, 1782) that principally inhabits kelp beds, rocky and mixed substrates (rock and sand), the gastropod mollusk Concholepas concholepas (Bruguiere, 1789) that is present in mixed/rocky habitats, and the bivalve mollusk Ameghinomya antiqua (King, 1832) which is mainly an inhabitant of soft substrates. All three species have a wide distribution along the Chilean coast (18 • -56 • S). These same species were evaluated in OAAs. Sampling was carried out during summer 2016 (9th of January till 14th of February) in most of the areas, except for three OAAs that were evaluated at the beginning of autumn 2016 (9th till 12th of April, Table 2).
In each area (MAs and OAAs), density of commercial species and mega-invertebrates were sampled by scuba diving surveys conducted from fishing boats. Five quadrats (0.25 m 2 ) were randomly placed within each transect (10 × 2 m) which were arranged parallel to the coast, at depths of up to 20 m (Supplementary Table S1). The total number of  transects per area was established according to their surface area, fluctuating between 25 and 42 for the selected commercial species (C. concholepas, L. albus, and A. antiqua), and between 15 and 29 for analysis of benthic community structure ( Table 2). On soft bottom benthic habitat, only the commercial species were evaluated with no analysis undertaken for the macroinfauna community. A total of 238 community transects were sampled, in which 19,818 individuals were identified to their lowest possible taxonomic level (Forcelli, 2000;Häussermann and Försterra, 2009); with the exception of a few taxa that could only be identified to higher taxonomic levels (e.g., Demospongiae, Actinaria, and Holothuroidea). Biological sampling (measured and weighed) of primary target resources (C. concholepas, L. albus and A. antiqua) was performed in all areas, except in the ones with low densities and scarce presence of individuals along transects (<20% of the evaluated transects). Between 269 and 526 individuals per each primary target resource per area were collected by scuba diving, immediately after they were measured and weighed on shore ( Table 2). For size of individuals test diameter was taken for L. albus and shell length for C. concholepas (peristome length: maximum length from the siphonal notch to the posterior edge of the shell) and A. antiqua (valvar length along the maximum axis). Size was measured using a Vernier caliper with a precision of 0.1 mm, and body wet weight (including the shell in mollusks) was measured using an electronic balance with a precision of 0.1 g.

Commercial Species
Performances of MAs relative to OAAs were evaluated using three indicators as a proxy for exploitation intensity of commercial species (C. concholepas, L. albus, and A. antiqua). Fishing regulations have established a legal minimum size of extraction (MLS) and fishers clearly target larger individuals; therefore, the absence of the largest harvestable individuals can be assumed as a proxy for exploitation intensity of a given species (Blanchard et al., 2005;Miethe et al., 2016). In this context, two fisheries indicators were chosen based on size. The first indicator considered the average size of the largest 5% of the sample (Lmax5) and was selected as a measure that is less affected by environmental effects and recruitment variability (Miethe et al., 2016). The second indicator relates the number of individuals in a specific size class to the total number of individuals collected [proportional stock density (PSD)]. The limit sizes adopted were the MLS established by Subpesca by the time the studied was carried out (C. concholepas = 100 mm; L. albus = 70 mm, and A. antiqua = 55 mm), and PSD therefore estimates the percentage of individuals whose size exceeds these MLS. A third indicator was not based on size, but instead considered average density of commercial species, estimated from the average number of individuals counted within each quadrat, divided by surface area (0.25 m 2 ). Average density of individuals for each area (ind.m −2 ) was computed from replicate transects. As a measure of the nutritional condition of each species, a fourth indicator was estimated, the Medium Condition Factor (MCF; Arana, 2006). This indicator relates the average weight of the individuals in a specific size range in one area to the average weight of all individuals in that specific size range in all the areas. To compare the weight of a similar number of individuals between areas, size range was chosen considering the class intervals that had the highest but even number of individuals between areas. Values exceeding one, indicate that individuals in one area have a better nutritional condition (i.e., suggests adequate food supply) than in other areas. This indicator had been used previously to assess gonadal development in L. albus over a time series (Arana, 2006). Since our study examined only spatial variability, sampling was carried out trying to study the areas during the same period to find individuals in similar reproductive conditions. According to reproductive cycle of L. albus, C. concholepas, and A. antiqua, during our sampling, post-spawning individuals in a reproductive rest state were collected and therefore it would be reasonable to assume that the MCF indicators are comparable reflecting relative nutritional condition. In the case of L. albus reproductive cycle in the northern neighboring region (41 • 45 ) and northern part of Aysén Region (Melinka, 43 • 53 ) indicated that a higher proportion of mature individuals can be found between spring and beginning of summer (Bay-Schmith et al., 1981;Arias et al., 1995, Molinet et al., 2016 which coincides with the reproductive seasonal closure established between 15th of October till 15th of January (D. Ext. N • 439, Subpesca, 2000). Studies of reproductive cycle of C. concholepas (41.8 • S; Manriquez et al., 2009) indicated maximum gonadal development in December. Seasonal closure of this species is established between the 1st of September and 31th of January (D. Ext. 697, Subpesca, 2011). A. antiqua, has no seasonal closure; maximum reproductive months had been observed between July and August in northern part of Aysén Region (43.8 • S), with spawning period occurring between spring and early summer (Canales et al., 2019).
Since data was not normally distributed (Kolmogorov-Smirnov, p < 0.01) and data presented no homogeneity of variance (Levene's test p < 0.05), Kruskal-Wallis non-parametric test (Statistica software 7.0) was used to independently compare the different indicators between groups.

Benthic Communities
Benthic communities were characterized, considering the four sampling sites described previously, relative to: (i) average density (number of individuals in each quadrat divided by its surface area of 0.25 m 2 with data standardized to ind.m −2 ) and (ii) dominance (% D, percentage of individuals of one species relative to the total number of individuals considering all species). Since density data was not normally distributed (Kolmogorov-Smirnov, p < 0.01), the Kruskal-Wallis non-parametric statistical test (Statistica software 7.0) was used to independently compare between the four groups. Additional measures of community structure were calculated, including total number of species (S), Shannon-Wiener diversity index (H ) and Evenness index (J ) (Pielou, 1977). The rarefaction method (ESn) was used to compare taxa richness in samples of unequal size for benthic macrofauna (Hurlbert, 1971). These measures were calculated using the software package PAST V4.03 (Hammer et al., 2001).
To determine which taxa contributed to the dissimilarity of the faunal assemblages between groups, an analysis of similarity was conducted (SIMPER, Clarke, 1993). Non-Metric Multidimensional Scaling (nMDS) and an analysis of similarity test (one-way ANOSIM) were used to evaluate similarities between groups, and for these analyses, species densities were transformed to the fourth root and the Bray-Curtis distance index was applied to estimate the degree of similarity between groups/species (Clarke, 1993). Transects were treated as replicates to provide an average density for each area. Differences between groups were evaluated using a one-way permutational multivariate variance analysis (PERMANOVA; Anderson, 2001). Since all OAAs were localized in Huichas, a more detailed statistical analyze was performed to compare them with MAs geographically closed. For this, nMDs and one-way ANOSIM analyses were performed using transects as replicas. SIMPER, nMDS, and ANOSIM analyses were completed using the software package PRIMER Version 6 (Clarke and Gorley, 2006).

Commercial Species
Concholepas concholepas was evaluated in all MAs studied since it is a primary target resource in all management plans. Suitable substrate in OAAs, indicated that four of the five areas studied should have had this resource, but it was present only in two of them. Mean densities varied between 0.4 and 3.8 ind.m −2 (Figure 2A) with significant differences (Kruskal-Wallis test, p < 0.05) in MAs (2.3 ± 5.6 ind.m −2 ) compared to OAAs (0.4 ± 1.3 ind.m −2 ). These differences were due to higher densities in MAs-Gala (2.6 ind.m −2 ) and MAs-Huichas (2.6 ind.m −2 ) compared to OAAs (Kruskal-Wallis test, p < 0.05). Within MAs-Huichas, the absence of this resource in MA-HUIC2 contrast with higher densities estimated in the other nearby areas (HUIC3 and HUIC1). In MAs-Gala high density was represented by small sized individuals (average 61 ± 11 mm) that resulted in Lmax5 of 84 (±15 mm) in comparison with MAs-Huichas where high densities were dominated by larger individuals (Lmax5 = 111 ± 11 mm) ( Figure 3A). Densities between both groups of MAs (Huichas and Gala) were not significantly different (Kruskal-Wallis test, p > 0.05) but when comparing Lmax5, significant difference were estimated (Kruskal-Wallis test, p < 0.05). Due to low densities in most of the areas, fisheries indicators (Lmax and PSD) could only be estimated in one OAAs preventing its statistical comparison with MAs. The percentage of individuals exceeding MLS (PSD) was 45% in MA-HUIC3, whereas in all other areas, PSD was <8%. MA-Gala1 represented an extreme example where none of the individuals sampled exceeded 100 mm (average size was 58 ± 9 mm). MCF, calculated over the size range of 70-89 mm, showed that individuals in MAs-Gala were also in poor nutritional condition compared to other areas (Table 3).
Proportional stock density which represented the percentage of individuals over the MLS (≥70 mm) showed high variability between areas (MAs; average 17 ± 15% and OAAs; average 10 ± 14%). In MAs, this indicator fluctuated between 0.3 and 55.7% (Table 2), with maximum values estimated in MA-HUIC3 and a minimum in MA-Gala2 (0.3%). In OAAs, PSD varied between 4.2 and 21.6%, with no significant differences between OAAs and MAs (Kruskal-Wallis test, p > 0.05). The Medium condition factor (MCF) considered size range between 40 and 59 mm showed that sea urchins in MAs-Gala had the lowest nutritional condition (MCF < 1, Table 2). No significant differences in MCF were found between MAs and OAAs (Kruskal-Wallis test, p > 0.05).

Benthic Communities
A total of 102 taxa were recorded in the 12 areas studied, 84 taxa/species in MAs and 71 taxa/species in OAAs. The most significant Phyla in terms of species richness were Mollusca (MAs n = 30; OAAs n = 24), Echinodermata (MAs n = 21; OAAs n = 16) and Arthropoda (MAs n = 12; OAAs n = 10). Of the total taxa registered in the survey, 21 corresponded to species of commercial interest, and most of these were mollusks ( Figure 4A and Supplementary Table S2).
The Echinoidea class represented mainly by Loxechinus albus, Arbacia dufresnii, and Pseudechinus magellanicus tended to dominate numerically (50-64% of total fauna density) in both MAs and OAAs. The bivalves Ameghinomya antiqua (11% of total fauna density) and Aulacomya atra (11% of total fauna density) also represented important numerical contributions to the overall density in MAs-Guaitecas and in OAAs, respectively ( Figure 4A and Supplementary Table S2).
Average densities of benthic species (considering all the community) showed significant differences in MAs-Gala (Kruskal-Wallis test, p < 0.05) in comparison with other areas ( Table 4). The Shannon-Wiener diversity index ranged between 2.09 and 2.33, indicating low variability between areas. Evenness index (J ) was ∼0.6 in all areas, indicating that the community was dominated by few species; no significant differences were estimated in either of the above indexes between the four groups ( Table 4; Kruskal-Wallis test, p > 0.05). Rarefaction curves indicate equal values of ES (100) = 22 both in MAs and OAAs (Figure 4B), although a trend of higher species richness was observed in MAs, with ES (2000) = 63 compared to OAAs with ES (2000) = 58. The shape of the curve for total-fauna, which considers all 102 taxa/species observed in MAs and OAAs, indicates that species richness was well represented by sampling effort in the studied areas ( Figure 4B).
Analysis of data using nMDS highlighted three main groups: (i) MAs-Gala, (ii) MAs-Guaitecas, and (iii) another cluster grouping all OAAs with MAs-Huichas ( Figure 4C). The one-way ANOSIM showed separation between groups (Global R = 0.32, p < 0.05) but only significant differences were found between MAs-Gala and OAAs (PERMANOVA test, p < 0.05). Results associated to the evaluation of the similarity between areas geographically closed (MAs-Huichas and OAAs) did not showed defined groups, associated to management regimes, in relation to the community composition and their abundances (one-way ANOSIM, Global R = 0.33, p < 0.05).
SIMPER similarity analysis was used to determine the species that contributed to the similarity (or dissimilarity) of the faunal assemblages within and between areas; high similarities within the areas (47-65%) were obtained, with group MAs-Gala having the highest percentage (65%). Dissimilarity between MAs and OAAs, showed MAs-Huichas as having  lower dissimilarity (49%) compared to OAAs and MAs-Gala as having greater dissimilarity (63%) compared to OAAs. The main species that contributed to these dissimilarities (MAs-Gala and OAAs) were Arbacia dufresnii, Loxechinus albus, Concholepas concholepas, Pentactella leonina, Tegula atra, Cosmasterias lurida, Argobuccinum pustulosum, Pseudechinus magellanicus, Ameghinomya antiqua, Actinaria unidentified, Aulacomya atra, and Metacarcinus edwardsii ( Figure 4D). However, only the first five of the above species showed significantly higher densities (Kruskal-Wallis test, p < 0.05) in MAs-Gala than in OAAs. In contrast to the other species, the sea star Cosmasterias lurida, showed significantly lower densities in MAs-Gala compared to MAs-Huichas and OAAs (Kruskal-Wallis test, p < 0.05).

Exploitation Intensity of Commercial Species
Concholepas concholepas is a species for which extraction has been banned from OAAs since 2002 (San Martin et al., 2010). Low densities in OAAs hinder the statistical comparison of fisheries indicators for C. concholepas between this regime and MAs. In most MAs, low densities combined with low values of the fisheries indicators could be indicative of a high exploitation intensity. To contribute to the discussion of these results, total allowable catch (TAC) for this resource and the others considered in this study, was estimated using the same methodology for MAs and OAAs (details of the methodology and quotas are given in Supplementary Table S3). In this context, only two areas obtained a quota for C. concholepas and only in one (MA-HUIC3), was high enough to be worth the extraction effort, suggesting that the exploitation intensity for this resource was high in most of the areas, despite their management regimes. Historical TACs for this resource in the studied MAs (obtained from Subpesca, 2020b) sustained this hypothesis since they had been decreasing since their first baseline study, except in HUIC3. Considering all MAs in Aysén, landings of C. concholepas have been decreasing since 2010, from 65,447 kg during the period 2010-2015, to 9,641 kg during the period 2016-2019, and in 2020 (up to October) only 1,180 kg has been harvested (data obtained from Servicio Nacional de Pesca, Sernapesca, 2020). OAAs were selected according to historical fishing grounds of the primary target resources, but no quotas were obtained suggesting illegal fishing activities. Data for illegal fishing shows that in the last 5 years, a total of 3,000 kg of C. concholepas has been confiscated (Sernapesca, 2020); although this is likely to be an underestimation of the true illegal harvest. In central Chile, a study of illegal extraction of this species concluded that fishers harvest illegally not only in OAAs, but also within their local MAs and within MAs of other organizations (Oyanedel et al., 2017). The same study estimated that official landings account for only 14-30% of the total C. concholepas catch. The massive extraction of this species within OAAs is not only of concern regarding total biomass, but also in the context of MLS; estimates suggest that up to 48% of this illegal catch is composed of undersized individuals, compared to only 8% in MAs (Fernandez et al., 2020). Our data therefore support the need to examine the effectiveness of the permanent ban on harvest of C. concholepas within OAAs (Bandin and Quiñones, 2014;Andreu-Cazenave et al., 2017). In contrast to the C. concholepas fishery, L. albus can be extracted from OAAs, but is managed through an annual TAC. Higher densities estimated in MAs compared to OAAs were mainly driven by small sizes individuals in MAs-Gala. Low values of the fisheries indicators in both MAs and OAAs, could be indicating that harvest had been intense in most of the areas, except in some MAs in Huichas. Estimated TACs supports these findings, since on average, low quotas were estimated for MAs and OAAs (11 ton and ≤3 ton in MAs and OAAs, respectively). In Aysén, increased exploitation of this resource in MAs had been registered since 2010 (Sernapesca, 2020). Indeed, from 2016, this species has contributed to 99.5% of landings within MAs. It appears that MAs present an opportunity for fishers to extract and commercialize this resource at a higher price in MAs (USD 0.7 ± 0.2 in OAAs and USD 1.5 ± 0.1 in MAs; A. Lafon, personal communication, December 12, 2020), at a time when the regional TAC had been already harvested in OAAs. MAs in Gala presented high densities of C. concholepas and L. albus but of small sizes individuals which could be indicating that this zone could be important in terms of recruitment. For both resources, the average sizes were in the lower limit of the first size at maturity; estimated between 70 and 90 mm for C. concholepas (41.8 • S; Manriquez et al., 2009) and between 41 and 44 mm for L. albus (Arias et al., 1995). A regular sampling (at least annually) that includes habitat characteristics (e.g., food supply, water exchange) will allow to determine if this situation is permanent and the area has become important in terms of recruitment or small sizes are the result of the intense exploitation of harvestable individuals.
Ameghinomya antiqua is a species restricted to habitats with sandy-muddy bottoms. This clam has only been targeted as a primary resource in four MAs (Huichas and Guaitecas) and was present in all OAAs. Unlike other commercial species evaluated in the present study, there are no restrictions for fishers to extract this species in OAAs, since no TACs had been implemented. Our results showed no significant differences in densities between management regimes. In both MAs and OAAs, size-frequency distribution was skewed toward larger sizes showing that exploitation intensity was low both in OAAs and MAs. Since suitable substrate for this species was generally restricted within these areas, estimated quotas were not high (average TACs: MAs = 26 ton; OAAs = 13 ton). MAs-Guaitecas were the only areas with historical TACs. In both MAs, these first studies estimated higher quotas (∼ 113 ton) that the ones estimated in our study. After that, stable quotas of around 30-50 ton (Subpesca, 2020b) could be indicating and supporting our findings of low exploitation intensity.

Benthic Community Structure
The benthic community was represented by a high number of species (102 taxa/species), dominated by Echinoidea (Arbacia dufresnii, Pseudechinus magellanicus, and Loxechinus albus). Diversity and evenness indexes did not show differences between areas with different management regimes. Only one group of  MAs (Gala) showed significant differences with OAAs regarding the densities of commercial species (L. albus and C. concholepas) and other secondary species (Arbacia dufresnii, Pentactella leonina, and Tegula atra).
High densities of L. albus (10 ± 15 ind.m −2 ) were comparable with those reported in northern Patagonia (11 ind.m −2 , Molinet et al., 2016;29 ind.m −2 , Contreras et al., 2019). L. albus is a strict herbivore Gonzalez et al., 2008) and therefore could have a controlling effect on macroalgae biomass and in surrounding communities, since macroalgae are three dimensional structures that favor settlement and serve as refuges for juveniles' stages (Dayton, 1985;Buschmann et al., 2004;Wright et al., 2005;Contreras et al., 2019). In MAs-Gala, high densities of sea urchins were observed, co-incident with low macroalgae abundance (<5%, Supplementary Table S4). These types of habitats dominated mainly by sea urchins and devoid of macroalgae, resemble the sea urchin barrens (Fillbee-Dexter and Schebling, 2014) that have been described in regions that support kelp beds as in Patagonia (Dayton, 1985;Fernandez et al., 2000). Sea urchin barrens can shift to a stage where kelp beds can recuperate or could maintain a stable state associated to longterm persistence of the parameters that drive and sustain high abundances of sea urchins (e.g., urchin grazing rate, kelp growth rate, recruitment rates, per capita predation rates) (Dayton, 1985;Fillbee-Dexter and Schebling, 2014). Threshold densities of sea urchins could determine these shifts: in fjords of the northern hemisphere, urchin densities of 45-75 ind.m −2 and 10-16 ind.m −2 , have been reported for barrens, and kelp beds, respectively (Fillbee-Dexter and Schebling, 2014). In northern and central Chile, barrens have been shown to alternate with kelps patches (Vasquez and Buschmann, 1997) and persistence of barrens has been associated with unregulated extraction of macroalgae (Vasquez, 2008). In the fjords and channels of Patagonia (45 • -54 • S) barrens dominated by L. albus have been reported (with densities up to 100 ind.m −2 ) (Dayton, 1985). In contrast, low densities of L. albus (<2 ind.m −2 ) have been reported in dense kelp forest in parts of the Magellan region (Friedlander et al., 2018). In MAs-Gala, high densities of L. albus (22 ± 14 ind.m −2 , with a maximum of 78 ind.m −2 ) and A. dufresnii (20 ± 14 ind.m −2 with a maximum of 51 ind.m −2 ) coincided with relatively low coverage of macroalgae which could be indicating a sea urchin barren in these areas. This contention is supported by juveniles of L. albus (39 ± 9 mm) and C. concholepas (61 ± 11 mm) dominating the community in these areas. High settlement in barrens has been reported, probably due to the lower structural complexity, which could limit available habitats for predators (Fillbee-Dexter and Schebling, 2014). Further research is needed to evaluate the persistence of this sea urchin barren state in this zone.
An unexpected finding of the present study was related to Cosmasterias lurida, a predatory sea star that inhabits the subtidal environment. This species was present in high frequency of occurrence (64% of transects) and higher densities (2 ± 3 ind.m −2 ) in the sampling sites compared to other studies in Patagonia Garrido, 2012). Recognized top predator such as C. concholepas (Navarrete et al., 2010) occurred infrequently and in low densities in our study, indicating that C. lurida would have no competitors being the highest-order predator. The intense harvesting of C. concholepas in the study area could be contributing to shifts on roles of predators that could have cascading effects on the structure of benthic communities (Fernandez et al., 2000;Pinnegar et al., 2000;Gaymer and Himmelman, 2008;Perez-Matus et al., 2017), and further evaluation of its effects is now a high priority.

Implications for Management of Fisheries and Benthic Environments in Patagonia
Our results indicate that MAs in the Aysén region are not fulfilling the objectives of sustainability regarding exploitation of benthic resources. Few of the MAs that were studied showed significantly lower exploitation intensities than OAAs, with most showing similar exploitation intensity. In both MAs and OAAs, benthic communities showed signs of disturbance characterized by the absence of the largest harvestable individuals (L. albus and C. concholepas). High densities of echinoderms (sea urchins and sea stars) were also indicative of changes in the community structure (e.g., trophic roles, competition, and dominances).
Currently, information provided by follow up studies is deficient in the sense that it does not permit evaluation of changes in benthic community structure associated to harvest. However, efforts have been made by the national fisheries institute (Spanish acronyms, IFOP) to evaluate MAs performance mainly in the northern regions. These results indicate that the community associated to C. concholepas presented two groups (north-central zone 28 • -33 • S and south zone 39 • -41 • S), with latitudinal differences in their abundances and composition of their functional groups. However, community indicators (e.g., density, diversity of functional groups, diversity, uniformity, and dominance indexes) resulted highly variable spatially and temporally (Ariz et al., , 2018Romero et al., 2019). Due to the high heterogeneity of these indicators, the use of a regional scale was recommended (Romero et al., 2019). A regional spatial scale will allow through environmental, biological and socioeconomic monitoring to identify more defined spatial patterns, and apply local management measures and policies (e.g., minimum extraction sizes). This monitoring should include integral spatial management since all commercial species considered in MAs management plans have long dispersive larval stages hence settlement, recruitment and growth are processes that probably greatly exceed the boundaries of MAs (Romero et al., 2019). In short, the productivity of one area is dependent on processes and management in other OAAs and MAs (Molinet et al., 2010;Arias and Stotz, 2020). It should be highlighted that our results are the first attempt to evaluate MAs' performance in comparison with OAAs regime on a regional scale with a community approach in Aysén (43.9 • -45.2 • S).
Illegal fishing has been recognized as a major impediment to marine conservation, and although MAs were designed to provide a disincentive, our data combined with regional statistics, highlight deficiencies in the existing system. The surveillance of fishers within MAs in the Aysén region has been almost impossible to implement because of long travel distances from fishing villages to their MAs (distances between coves and MAs in this study varied between ∼9 and ∼32 km and up to 180 km when considering all regional MAs). It should be highlighted that most fishers own small boats which resulted in 1 or 2 h sailing from their coves to their MAs. In this context, the only area (HUIC3) showing improved fisheries indicators for commercial species was one with year-round surveillance provided by the proximity of local fishers. In this case, fishermen asked for a land concession next to the area to diminish cost associated to surveillance. This surveillance system worked well because of social cohesion within their union, which highlighted the social capital as a key factor in the successful functioning of MAs (Marin and Gelcich, 2012). MAs-Gala were the only other management areas close to coves (∼9 km apart) but due to low social capital no vigilance has been implemented. Social capital is especially important in Aysén region because fishing coves are not only remote from MAs but also from markets. Clearly, the establishments of collaborative networks should be considered as a priority in the strengthening of this benthic management regime.

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
MH: conceptualization, methodology, writing -reviewing and editing, supervision, and visualization. PO: conceptualization, methodology, and writing -reviewing and editing. Both authors contributed to the article and approved the submitted version.

FUNDING
This study was financially supported by the Subpesca through grant CUI 2015-7-FAP "Asistencia para el desarrollo de Estudio de Situación Base de Área y Seguimiento en Áreas de Manejo de la Región de Aysén, Etapa 1". Regional Program ANID R17A10002 and R20F0002 partial funding.