Oyster Reefs in Northern Gulf of Mexico Estuaries Harbor Diverse Fish and Decapod Crustacean Assemblages: A Meta-Synthesis

Oyster reefs provide habitat for numerous fish and decapod crustacean species that mediate ecosystem functioning and support vibrant fisheries. Recent focus on the restoration of eastern oyster (Crassostrea virginica) reefs stems from this role as a critical ecosystem engineer. Within the shallow estuaries of the northern Gulf of Mexico (nGoM), the eastern oyster is the dominant reef building organism. This study synthesizes data on fish and decapod crustacean occupancy of oyster reefs across nGoM with the goal of providing management and restoration benchmarks, something that is currently lacking for the region. Relevant data from 23 studies were identified, representing data from all five U.S. nGoM states over the last 28 years. Cumulatively, these studies documented over 120,000 individuals from 115 fish and 41 decapod crustacean species. Densities as high as 2800 ind m-2 were reported, with individual reef assemblages composed of as many as 52 species. Small, cryptic organisms that occupy interstitial spaces within the reefs, and sampled using trays, were found at an average density of 647 and 20 ind m-2 for decapod crustaceans (range: 18 - 2843) and fishes (range: 2 - 56), respectively. Both groups of organisms were comprised, on average, of 8 species (range: 3 - 15). Larger-bodied fishes captured adjacent to the reef using gill nets were found at an average density of 6 ind m-2 (range: 1- 68) which came from 23 species (range: 9 - 46). Decapod crustaceans sampled with gill nets had a much lower average density, < 1 ind m-2, and only contained 2 species. On average, seines captured the greatest number of fish species (n = 33), which were made up of both facultative residents and transients. These data provide general gear-specific benchmarks, based on values currently found in the region, to assist managers in assessing nekton occupancy of oyster reefs, and assessing trends or changes in status of oyster reef associated nekton support. More explicit reef descriptions (e.g., rugosity, height, area, adjacent habitat) would allow for more precise benchmarks as these factors are important in determining nekton assemblages, and sampling efficiency.


INTRODUCTION
Understanding the impacts of habitat change on natural resources remains a key component for informing restoration and management policy (Barbier et al., 2011;Bennett et al., 2015). Fisheries policy in particular, through the 1996 amendment to the U.S. Magnuson-Stevens Act (passed 1976), introduced a mandate on defining and protecting Essential Fish Habitat, which include "all waters and substrate necessary to fish for spawning, breeding, feeding, or growth to maturity" (Magnuson Stevens Fishery Management Conservation Act, 1996, 16 U.S.C. section 1801-1804. This placed high importance on specieshabitat associations (Magnuson Stevens Fishery Management Conservation Act, 1996: 50 CFR, sections 600.805-930). Since then, data on habitat-specific fisheries and associated species have been used to justify, inform, and guide policy and management activities (Rondinini and Chiozza, 2010;Vasconcelos et al., 2015;NRC, 2017). Thus, documenting species richness and assemblages within a habitat type is a first crucial step to defining habitat support and ensuring collection of proper baseline data. These baseline data can be used to support modeling, evaluate restoration outcomes (Dufrêne and Legendre, 1997;NRC, 2017), and, as highlighted recently, to evaluate resource injury (i.e., under the U.S. Oil Pollution Act, 1990 1 , section 1006I(1) for coastal resources).
Reefs built by the eastern oyster, Crassostrea virginica, have been recognized as Essential Fish Habitat (Coen and Grizzle, 2007). Historically valued for their economic impact as a direct food commodity, oysters also create reefs, providing valuable three-dimensional habitat within coastal environments. In recent decades, significant efforts to conserve, and restore oyster reefs have been justified based on the valuable ecosystem services they provide, including water quality improvements, shoreline protection, and habitat creation for commercially and recreationally important fisheries (Coen and Grizzle, 2007;Grabowski et al., 2012). Importantly, for fisheries and restoration policy, oyster reefs are recognized as key biogenic habitat for a diverse assemblage of fishes and decapod crustaceans (e.g., Mobius, 1877;Frey, 1946;Wells, 1961;Coen et al., 1999;Grizzle, 2007, 2016). The reported functional decline in oyster reefs (Beck et al., 2011;Zu Ermgassen et al., 2012) is likely to have broad consequences for habitat provision, and thus biodiversity and fisheries production (Peterson et al., 2003;Humphries and La Peyre, 2015;Zu Ermgassen et al., 2015), but a lack of established benchmarks hinders our ability to assess impacts.
Recent guidelines for monitoring oyster reefs suggest setting explicit goals for assessing habitat support (Baggett et al., 2015;NRC, 2017). Goals for restoration projects may focus simply on provision of habitat (i.e., La Peyre et al., 2014a;Baggett et al., 2015), or, they may be expanded to include local enhancement of ecosystem services (i.e., Coen and Luckenbach, 2000). Central to either of these goals, however, is an understanding of what the nekton assemblages look like on the desired habitat in terms of expected assemblages, abundance, and biomass of species.
As such, many oyster reef restoration planning documents call for the use of project-specific reference sites (Coen et al., 2004;SER Society for Ecological Restoration, 2004;Baggett et al., 2015;NRC, 2017); however, such sites are often not available, or monitoring efforts fail to collect these data. An alternative is to establish desirable conditions, or benchmarks to measure changes in resource status (Ehrenfeld, 2000;Kentula, 2000).
Benchmarks provide a tool to assess the status of natural, managed, or restored ecosystems (Angermeier and Karr, 1994;McClanahan et al., 2019). Ideally, benchmarks are set using baseline data derived from natural or pristine systems, which are defined by functional and evolutionary limits of the ecosystem (Pickett et al., 1992). For many resources and regions, effects of altered landscapes (e.g., river management, climate change), and a lack of historical data confound efforts to set benchmarks representing pristine conditions (Toledo et al., 2011). Lacking these historical data, useful benchmarks can still be established for current ecosystem status against which future assessments of natural resource condition or management effects can be measured. Across the northern Gulf of Mexico (nGoM), significant coastal restoration, river management, and climate change impacts have affected local hydrology, water quality, and landscape configurations (Pendleton et al., 2010;Montagna et al., 2011). These changes alter oyster population dynamics and reef characteristics which ultimately impacts the oyster reef 's provision of ecosystem services, including nekton habitat. Benchmarks representing reef-resident and transient assemblages based on the current range or means of values from recent assessments would provide useful tools to assess changes in status of natural reefs, as well as managed and restored reefs.
The nGoM supports extensive natural oyster reefs (Kilgen and Dugas, 1989), and recently has experienced significant restoration activities related to oyster reefs, with a goal of providing fisheries habitat (La Peyre et al., 2014a). Reports from the Atlantic coast of the United States document occupancy of eastern oyster reefs by 50 to over 300 species, with the wide range of reported species depending on whether the surveys include fish, decapod crustaceans, molluscs or other organisms, such as protozoans, sponges and flatworms (e.g., Frey, 1946;Wells, 1961;Coen and Grizzle, 2007). Similar data do not exist for the nGoM, despite the exceptional fisheries support and habitat characteristics of the region. For example, the nGoM is characterized by a microtidal regime (<0.5 m) where water level changes are driven mostly by winds, precipitation, and runoff (Kim and Park, 2012). As such, depth and exposure of Gulf coast oyster reefs can vary and be difficult to predict (Kilgen and Dugas, 1989). Importantly, water depth and substrate or reef exposure influence reef formation, which impacts oyster population dynamics (Ridge et al., 2015). Oysters in the nGoM have higher growth rates and reduced mortality rates compared to other regions, despite experiencing lower salinity, and higher temperatures than other regions (Lowe et al., 2017). These differences in reef formation, water quality, and growth characteristics likely influence the associated nekton communities, suggesting a need for region specific benchmarks.
This study aims to quantify fish and decapod crustacean assemblages associated with nGoM oyster reefs. The goal is to provide benchmarks using available data collected on nekton occupancy across a broad spectrum of oyster reefs. These benchmarks will be useful as a guideline to assess policy and restoration outcomes in the nGoM, as well as inform future injury assessments (Geist and Hawkins, 2016). These needs have been identified in previous studies and planning documents and represent a critical information gap for assessing restoration and policy goals (Baggett et al., 2015;Coen and Humphries, 2017;NRC, 2017).

Peer Reviewed Study Selection
A literature search for fish and decapod crustacean sampling on oyster reefs across the nGoM was conducted by searching Google Scholar using the following search terms: "oyster reef " AND (Texas OR Louisiana OR Mississippi OR Alabama OR Florida) AND (macrofauna OR nekton), for all papers dated through April 2019. Only peer-reviewed, regional nGoM shallow water (<3 m) studies that provided oyster reef location, sample size, sample methods (i.e., gear), species identities, and density information were included in the final database. Also, only studies targeting oyster reef fish and/or decapod crustacean communities (and not individual species focused studies) were included for the targeted evaluation of community composition, and not specific species.

Nekton Variables
Nekton variables were collected by gear type, and standardized across studies. To address the issue of unequal sampling effort, all nekton density data were standardized to individuals m −2 for all gear types, except gill net, which was standardized to individuals hr −1 of net soak time. Sampling effort was calculated by multiplying the number of samples taken by either the total area sampled, or total hours fished.

Reef and Water Quality Characteristics
Variables describing reef characteristics and water quality, where available, were added to the database. Data availability was examined for tidal height, reef elevation, exposure, and water quality for selected study reefs and locations. Water quality data for salinity, temperature ( • C), and dissolved oxygen (mg L −1 ) were the only three variables consistently available. We recorded reef type (natural vs. constructed) as there were few data sets for natural reefs, and using data from only natural reefs resulted in a very limited database. For the purposes of this study, we define natural reefs as reefs that had not, to the best of our knowledge, been anthropogenically enhanced, restored, or created in any way. We define "constructed" reefs as those that were anthropogenically enhanced, restored, or created using shell or other materials such as rock or cement with the goal of reef restoration and not oyster production for human consumption. Reefs were classified as intertidal or subtidal based on information from the original studies. Any reef not explicitly identified as intertidal was classified as subtidal, although we acknowledge that some subtidal reefs in this region are occasionally exposed from weather events (i.e., La .

Final Database
The final database included density or catch per unit effort (CPUE) of fish and decapod crustaceans, number of fish and decapod crustacean species, list of species collected, sample size, sampling gear, reef type (natural, constructed), tidal status (intertidal vs. subtidal), and a few discrete environmental parameters (water temperature, salinity, dissolved oxygen).
Data are presented and discussed using statistical means and distributions by gear type. We use the means for fish and decapod crustacean density and species number separately, by dominant gear type quantifying "on-reef " and "near reef " values to suggest potential benchmarks based on the range and means of values currently found across the nGoM (R 3.5.1; R Core Team, 2018).

RESULTS
From a total of 1,072 search results, 23 articles were selected for this analysis based on their overall goals and data availability (Figure 1; Table 1). The majority of studies were excluded because they focused on a region other than the Gulf of Mexico, the study was not related to oyster reefs, or data were not reported for all species. Within the Gulf of Mexico, there were many studies examining various aspects of fisheries, nekton production, and trophic changes on oyster reefs, but they failed to meet our criteria for inclusion in the final dataset. A number of studies focused exclusively on indicator species, economically important species, a subset of the reef community, or failed to report all catch data. While valuable, such studies were not included. The studies selected for analyses covered all five Gulf Coast states and used multiple gear types. Texas had the most data (9 studies), followed by Louisiana (6 studies), Florida (3 studies), Alabama (2 studies), and Mississippi (1 study), and 2 studies covered multiple states (Figure 2). Across the 23 studies, seven different gear types were used with most studies using more than one gear type on a single oyster reef. Trays (TR) were the most frequently used gear type (n = 12). Gill nets (GN) were the next most used gear type (n = 8), although the lengths, mesh sizes and soak times varied. Drop-samplers (DS), throw traps (TT), and lift nets (LN) were grouped in one category (DS/TT/LN; n = 8). Seines (SN; n = 4) and epibenthic sleds (ES; n = 3) represented the remainder of sampling gear included in these analyses. We classified TR, DS/TT/LN and ES as sampling "on-reef, " and GN and SN as sampling "near-reef " assemblages.
Reefs sampled across the selected studies included a mix of reef types, tidal status, and location. Studies sampled either natural reefs (n = 12) or constructed reefs (n = 9), with two studies comparing natural and constructed. Constructed reef habitats (n = 9) were built at different times  using different base materials that included shell, limestone, and concrete base materials. There was a mix of intertidal and subtidal reefs (intertidal = 5; subtidal = 18). Due to lack of site-specific water level data, we were unable to determine exposure levels of reefs or exposure time. Intertidal reefs were all located nearshore (identified as < 25 m from shore), along with a few of the subtidal reefs (n = 3); the remainder of reefs were located offshore (>25 m from shore).
Water quality parameters commonly reported included salinity, temperature, and dissolved oxygen. Oyster reefs were FIGURE 2 | Map of study sites used to examine the diversity and abundance of fish and decapod crustacean communities on oyster reefs across the northern coast of the Gulf of Mexico (Texas, Louisiana, Mississippi, Alabama, and Florida Gulf coast). See Table 1 for explanation of study codes.
sampled through a broad range of salinities (0-39.5) from Florida to Texas with most studies reporting ranges of salinities >20, reflecting the highly variable nature of the estuarine environments where many oysters are found in the nGoM ( Table 1). Temperature reflected the annual range in the region, from 2.0 • to 40.6 • C, encompassing winter and summer samples, while dissolved oxygen conditions ranged from 0.4 to 17.3 mg L −1 . No other environmental parameters were reported consistently across studies.
A total of 115 fish species, based on over 32,000 individual fish, were reported on or near oyster reefs ( Table 2,  Supplementary Table 1). Of the 115 reported fish species, 22 species (∼19%) were collected in eight or more of the studies (Table 3), and 40 species were only reported in one study (Supplementary Table 1). The number of fish species collected ranged from 4 to 46 per study ( Table 2; Figure 3;  Supplementary Table 1). Species assemblages varied among gear types with TR and DS/TT/LN capturing more known reef residents, such as gobies and blennies, which included 15 species found only in these gear types. Gill nets captured larger and more mobile transient fish species, such as sheepshead (Archosargus probatocephalus) and red drum (Sciaenops ocellatus), while seines captured a mix of both the transient fish species and some of the smaller more ubiquitous resident species (Table 3;  Supplementary Table 1). Seine and gill net data included 51 fish species not captured with any of the other methods. Of these, 23 species were captured exclusively by gill net, including many of the shark species, and 15 species were captured exclusively by seine. Combined, the SN and GN gear types uniquely captured 44% of the fish species (Supplementary Table 1). Densities reported were gear-dependent, with TR and TT/DS/LN having mean densities of 20.2 ± 4.3 ind m −2 (mean ± SE) and 15.6 ± 2.8 ind m −2 , respectively (Figure 4). GN, capturing larger bodied fish, reported mean densities of 6.7 ± 8.4 ind m −2 . SN and ES reported low densities (< 1 ind m −2 ).
A total of 56 decapod crustacean species, based on over 90,000 individuals were identified ( Table 2, Supplementary Table 2). The number of decapod crustacean species ranged up to 26 per study. Of the 56 reported crustacean species, 11 species (snapping shrimp, blue crabs, depressed mud crabs, brown shrimp, white shrimp, Gulf stone crabs, Palaemonetes pugio and P. vulgaris grass shrimp, Panopeid mud crabs, green porcelain crabs, Harris mud crabs) were listed in eight or more of the studies ( Table 4) while 27 species were only listed in one study (Supplementary Table 2). Species composition differed between gear types, with TR, TT/DS/LN, and ES capturing the most species. TR consistently captured highest densities (mean ± SE: 647.9 ± 245.9; Figure 4), while GN and SN captured few decapod crustaceans (Figure 4, Supplementary Table 2).

DISCUSSION
Modern day natural and constructed reefs across the nGoM support diverse and dense assemblages of fishes and decapod crustaceans. Nekton density was as high as 2,800 ind m −2 and some reefs had >50 species of fish and decapod crustaceans. Using the available published data, gear-dependent means for density and species richness of dominant "on-reef " sampling (tray) and "near reef " sampling (gill net) serve as useful benchmarks. These benchmarks provide quantitative measures that can be used to assess changes in the status of existing and constructed oyster reefs. Management of coastal habitat and fisheries requires data on species-habitat associations to delineate Essential Fish Habitat, better implement restoration, perform injury assessments, and set policy goals. These suggested benchmarks, based on current reef status, may be used to assess  trends in oyster reef habitat support of fisheries and ecosystem functioning for this region.
Relating specific assemblages or densities to reef characteristics or across locations was problematic due to a lack of gear standardization, established gear conversion factors for oyster reefs, and reef habitat metrics. For example, some studies only used gill nets and seines to sample whereas others used drop samplers; given differences in catch efficiencies based on species' identity and size, and variable efficiency of different gear types across reefs, results from different gear types are not directly comparable (Zimmerman et al., 1984). As more data are collected, methodologies are standardized, or conversion factors developed, these fish and decapod crustacean benchmarks will become more targeted to specific locations, reef complexity or reef type.

Reef Characteristics and Complexity
Reef characteristics, such as adjacent habitats, connectivity, habitat redundancy, and water quality all affect fish and decapod crustacean assemblages, and are thus vital metrics for interpreting benchmark data. In general, data synthesized for this study lacked consistent reporting of reef location, in terms of adjacent habitats, connectivity or habitat redundancy, or complexity even though these factors influence nekton assemblages. For example, in Texas, Nevins et al. (2014) hypothesized that low faunal densities observed on natural oyster reefs may be a result of adjacent habitats, as well as difficulties sampling the complex habitat. Tolley et al. (2006) highlighted the impacts of salinity and freshwater flow on reef assemblages identifying flow rates and salinity as key location characteristics influencing reef communities. For constructed reefs, Gregalis et al. (2009) demonstrated that reef height affected fish abundance, while resident species abundance and transient fish assemblages varied by reef location. In some cases, reefs have been proposed as potentially redundant habitat due to their location adjacent to or near other high-quality habitat (i.e., Geraldi et al., 2009;La Peyre et al., 2014a;Heck et al., 2017). Reporting distance to adjacent habitats, reef exposure, and water quality/level at sampling would be useful for studies to better understand drivers of fish and decapod crustacean occupancy of oyster reefs.
Habitat complexity may also affect nekton assemblages through direct and indirect effects on trophic cascades, predation and habitat availability (e.g., Grabowski et al., 2008;Humphries et al., 2011a). Quantifying complexity, particularly across various constructed reef materials remains challenging. In particular, while habitat complexity may be important, some studies have suggested a potential complexity threshold due to failure of finding increasing nekton numbers with increased complexity (e.g., Humphries et al., 2011b;George et al., 2015). For example, George et al. (2015) found no difference in nekton assemblages or prey mortality in experimental studies comparing five possible substrate materials. Similar densities of resident nekton despite increasing oyster density led the authors of one study to hypothesize that there might be a low threshold for habitat complexity (Beck and La Peyre, 2015). These studies highlight the difficulty of relating nekton density and assemblage composition to habitat complexity; some of the differences, however, might also relate to faulty comparisons across studies as a result of sampling gear issues.

Impact of Sampling Gear
Sampling gear are selective for specific size ranges or species, and are not equally effective across, or within, complex habitats (Rozas and Minello, 1997). For oyster reefs in particular, the natural and constructed reefs include a wide range of reef complexity, heterogeneity, reef sizes, reef history, and reef locations. All of these factors influence the assemblages found on or around oyster reefs (Grabowski et al., 2005;Luckenbach et al., 2005;Geraldi et al., 2009;Nevins et al., 2014;Beck and La Peyre, 2015), as well as the effectiveness of gear types. While gear-dependent benchmarks for fish and decapod crustacean density and species number provide general region-wide values, the effects of the interaction of gear type with reef complexity, type, and location remain unknown.
Our study highlighted that different gear types capture different assemblages and numbers of nekton. The highest densities of fish and decapod crustaceans were consistently quantified using methods sampling "on-reef " despite reporting the lowest sampling effort. Specifically, trays, and/or enclosure samplers (lift nets, throw traps, drop-samplers), report sampling "on-reef, " where they capture species generally occupying interstitial spaces within the reef structure. Some variance between these gear types likely reflects differences in location as they cannot be used across all reef types (i.e., throw traps and drop-samplers require water depth < 1. 5 m), or require small patches of reefs where the sampler can fully enclose the reef (i.e., Stunz et al., 2010;Humphries et al., 2011a). Substrate trays sampled similar resident faunal assemblages as other "on-reef " gear, but can be used at greater depths and on larger reef patches. Trays are often criticized for allowing organisms to escape during retrieval in deeper waters, although modifications including nets that can be drawn closed have been suggested (Beck and La Peyre, 2015). Overall, resident decapod crustaceans (e.g., panopeid crabs) and resident fishes (e.g., gobies, blennies) were ubiquitous on the nGoM oyster reefs. In coral reef systems, these cryptobenthic fauna have been shown to provide as much as 70% of the energy consumed in the ecosystem (Brandl et al., 2019). However, many of these species were captured only with tray sampling, suggesting that without these sampling approaches, an important part of the oyster reef community would not be captured.
In contrast, lower numbers were generally captured using gear sampling adjacent to the reef; specifically, by seines and gill nets. Due to the nature of the reef, these techniques limit sampling to near the reef. These approaches however captured the larger commercially and recreationally important fish species of interest, such as red drum (Sciaenops ocellatus), and Spanish mackerel (Scomberomorus maculatus). For both gill net and seine sampling, variation in mesh size, time of day, water flow, and distance covered all influence catch rates (e.g., Vandergroot et al., 2011;Hubert et al., 2012). While techniques exist for standardizing data for some of these differences (i.e., mesh size; Shoup and Ryswky, 2016), other details may be more difficult to reconcile (i.e., interaction with hydrology, reef characteristics), and are often not fully reported.
Gear conversion factors have been developed for comparing different sampling gear, using a variety of statistical techniques (i.e., Pelletier, 1998;Gibson-Reinemer et al., 2016). This might be a useful technique to standardize reported densities and species' numbers for oyster reefs, if effects of habitat location and complexity can be reconciled. For example, even when similar gears can be used, their catch efficiencies and deployment techniques can vary dramatically across conditions. Drop samplers were found to have over 90% catch efficiency at small Texas oyster reefs (Zimmerman et al., 1984;Stunz et al., 2010), but this rate is highly correlated to water clarity and wave action (Rozas and Minello, 1997), two variables that are not frequently reported. Improved reporting of full environmental conditions, standardization of gear type for specific habitats, and, as data permit, exploring possible gear conversion factors could help in developing more targeted reef benchmarks than what we present here.

Reef Restoration, Monitoring, and Benchmarks
If we consider the current state of oyster reefs to be the desired reference condition for assessing changes in status and setting goals and policy for restoration (Coen and Luckenbach, 2000;NRC, 2017), this synthesis provides a set of benchmarks based on means and ranges of current values found across the region. These benchmarks help provide an abundance, density and composition of fish and decapod crustaceans expected to occupy a reef. We suggest these benchmarks as gear-dependent goals, using trays and gill nets, the two most versatile gear types used to sample "on-reefs" (tray), and "near-reefs" (gill nets). To effectively develop and use this approach, habitat characteristics that may be driving the differences between reef assemblages, densities, and species richness need to be better documented. Specifically, habitat characteristics, including reef location, complexity, water quality, and reef exposure critically influence occupancy of reefs by nekton, and benchmarks should be developed to reflect these variations. These same habitat characteristics, along with sampling conditions (i.e., winds, currents, waves, tidal conditions) may also impact sampling gear efficiency. Suggested benchmarks should be treated as a general guide for this region, with adjustments made based on knowledge of local reef habitat characteristics and conditions during sampling. Ideally, future data collection will provide for more targeted benchmarks for this region.
Over the last decade, the oyster reef restoration community has developed and increasingly follows detailed guidance for selection of restoration sites (Coen et al., 2004, and further discussed in Coen and Humphries, 2017) and monitoring reef restoration (i.e., Baggett et al., 2015). However, the difference between these original criteria for guiding site selection, and assessing occupancy by fish and decapod crustaceans, is that the characteristics of the species assemblage using a reef may be dependent on variables that do not necessarily limit oyster reef development (e.g., structural complexity). As a result, better quantification of occupancy of oyster reefs by fish and decapod crustaceans would significantly help in predicting effects of reef complexity, or location (Gilby et al., 2018).
With enormous investments targeted for habitat restoration in the nGoM, and continued emphasis on habitat-species linkages through Essential Fish Habitat policy, quantitative and standardized baseline data to establish benchmarks are increasingly important for managers and policy-makers (Baggett et al., 2015;NRC, 2017;Blomberg et al., 2018a). In the nGoM, over 6 billion USD has been designated for restoration of ecosystems (Environmental Law Institute, 2016) as a result of injury settlements. To ensure effective use of these funds, and to support existing policies related to fisheries and habitat management, generation of standardized monitoring metrics and clear benchmarks to help assess restoration and policy outcomes remains critical.

DATA AVAILABILITY STATEMENT
Datasets were taken from the published literature identified in Table 1.