Key Life History Attributes and Removal Efforts of Invasive Lionfish (Pterois volitans) in the Flower Garden Banks National Marine Sanctuary, Northwestern Gulf of Mexico

Indo-Pacific lionfish (Pterois volitans and P. miles) are the first marine teleost to have successfully invaded and become established in the Western Atlantic Ocean of the United States, Gulf of Mexico (GoM), and Caribbean Sea. Pterois volitans were first reported in the Flower Garden Banks National Marine Sanctuary (FGBNMS), a protected coral reef system in the northwestern GoM, in 2011. Little is understood about the life history characteristics of lionfish in this ecosystem. This study assessed population characteristics (size, density, age and growth) and removal efforts of lionfish (n = 1,665) at two coral reef sites within FGBNMS for 2015, 2016, and 2018. The annual increment formation in sagittal otoliths was examined to assess the age and growth of lionfish collected in 2018 (n = 100). Lionfish ranged in size from 75 to 444 mm total length (TL) and 4–1,153 g in total weight (TW). Six hundred and ten fish were randomly dissected for sex determination (females = 256, males = 354), females ranged in size from 137 to 348 mm TL and 21–586 g, while males ranged from 118 to 444 mm TL and 18–1,153 g. Interannual variation in mean lionfish density ranged from 26.7 individuals per hectare (ind ha–1) in 2016 to 81.1 ind ha–1 in 2018, while removal effort significantly increased ranging from 1.92 to 5.42 kg diver h–1. Lionfish age ranged from 0 to 10 years, with a mean age of 3.9 years. The observed values of the asymptotic maximum total length (L∞) and Brody’s growth coefficient (K) were 345 mm and 0.30 for females and 415 mm and 0.18 for males. Results suggest lionfish from FGBNMS exhibit markedly lower mean densities, a lower L∞ and growth rate, but attain older ages than lionfish in the Caribbean Sea, Western Atlantic Ocean, and other ecoregions in the northern GoM. This study describes the first key life history parameters and removal efforts for lionfish in a protected, healthy coral reef system in the northwestern GoM that may provide insight into environmental population controls (e.g., ecological resilience). Metrics from this study could be integrated into mechanistic ecological models to determine if FGBNMS is in fact exhibiting natural resilience to the lionfish invasion.


INTRODUCTION
Indo-Pacific lionfish (Pterois volitans and P. miles, herein referred to as lionfish) are non-indigenous marine fish that were introduced in the mid-1980s along the Western Atlantic Coast of the United States (U.S.), most likely through release from aquaria (Semmens et al., 2004;Hunter et al., 2021). Their distribution has since expanded into the Caribbean Sea, throughout the Gulf of Mexico (GoM), and into the Southwest Atlantic Ocean along South America (Aguilar-Perera and Tuz-Sulub, 2010;Schofield, 2010;Ferreira et al., 2015;Campbell et al., 2022). More recently, lionfish have invaded the Mediterranean Sea (Azzurro et al., 2017;Bariche et al., 2017) via the Suez Canal. Lionfish appear to be thriving in the warm temperate, and subtropical, waters of the Western Atlantic Ocean, Caribbean Sea, and GoM, as their densities are an order of magnitude higher than observed in their native range in the Indo-Pacific (Grubich et al., 2009;Kulbicki et al., 2012;Elise et al., 2015). Several biological characteristics of lionfish have been well documented to explain their rapid expansion and sustained population levels in invaded ranges. Life history characteristics facilitating the population expansion of lionfish are rapid growth, early sexual maturity, and frequent spawning capabilities (Betancur-R et al., 2011;Cote et al., 2013a;Gardner et al., 2015;Fogg et al., 2017). Lionfish exhibit physiological tolerances to a breadth of habitat and environmental conditions (Kimball et al., 2004;Biggs and Olden, 2011;Nuttall et al., 2014;Jud et al., 2015) and are opportunistic generalist feeders (Peake et al., 2018). Finally, lionfish appear to be resistance to disease and parasites (Simmons, 2014;Sellers et al., 2015;Fogg et al., 2016;Montoya-Mendoza et al., 2017;Aguilar-Perera et al., 2018;Harris et al., 2020) and lack natural predators (Maljkovic et al., 2008;Ulman et al., 2021), rendering their populations virtually unchecked. Lionfish have the potential to alter marine communities by consuming herbivorous fishes responsible for controlling algal production (Hughes, 1994;Albins and Hixon, 2008;Barbour et al., 2011), competing for resources with native meso-predators (Raymond et al., 2015;Marshak et al., 2018), and reducing local reef-fish populations (Arias-Gonzalez et al., 2011;Green et al., 2012a;Albins, 2013;Cote et al., 2013b).
It is important to gather regional baseline demographic data to inform management of this invasive species. Population dynamics and age-specific life history parameters are fundamental to understanding factors affecting recruitment success and are essential inputs in population models to evaluate efficacy of fishery management policies (Campana, 1999;Mendoza, 2006;Pope et al., 2010;Kindsvater et al., 2016). These data become particularly important when non-native species establish in a new region, as this can inform management and assess the effects they have on native species and the environment they invaded. Changes in size, age structure, and sex ratio of fishes are common indicators of overfishing (Berkeley et al., 2004;Cope, 2013), but in the case of invasive species, may be used to determine the success of management strategies used to reduce their populations (e.g., targeted removals; Frazer et al., 2012) or whether natural control may be occurring (Harris et al., 2020).
The northwestern Gulf of Mexico (generally recognized as the area from the Texas-Mexico border to the state line between Louisiana and Mississippi; New World Encyclopedia, 2021) is marked by well-established coral communities and banks that encompass biologically and ecologically important features along the outer continental shelf (Schmahl et al., 2008). The Flower Garden Banks National Marine Sanctuary (FGBNMS) represents the northernmost coral reef system in the Gulf of Mexico (GoM) and was designated a national marine sanctuary by the Office of National Marine Sanctuaries, within the National Oceanic and Atmospheric Administration (NOAA), in 1992 (Federal Register, 1991). Originally comprised of three banks (East Flower Garden Bank (EFGB), West Flower Garden Bank (WFGB), and Stetson Bank), this marine sanctuary was expanded in 2021 to encompass 17 bank systems located 110-304 km offshore from Galveston, Texas (Federal Register, 2021). EFGB and WFGB are populated by a rugous, complex reef structure and are arguably some of the healthiest coral reefs in the Caribbean region, with 55-60% living scleractinian coral cover on the reef caps (16-40 m) that has remained stable or increased significantly over the last 30 years (Johnston et al., 2021). Lionfish were first observed at EFGB and WFGB in 2011 (Johnston et al., 2016b), though recent records estimate their invasion may have begun as early as 2008 . This study reports life history characteristics and removal efforts of lionfish collected from EFGB and WFGB across 3 years.
Lionfish density, weight-length relationship, age, and growth have been suggested to differ among invaded regions likely attributed to invasion chronology, environmental influences (e.g., resource availability and accessibility, and climate), and removal efforts (Barbour et al., 2011;Farquhar, 2016;Fogg, 2017). The aim of this study was to document key life history attributes of lionfish on a protected coral reef system (i.e., FGBNMS) in the northwestern GoM, and to determine if differences were evident for age at length and age structure composition among invaded ecosystems. In addition, catch per unit effort (CPUE) was estimated for each year as another mechanism to describe trends in the lionfish population.

Study Area
EFGB and WFGB are prominent geological features located near the edge of the continental shelf and were created by the uplift of underlying salt domes (Hickerson and Schmahl, 2005; Figure 1). The sanctuary boundaries of EFGB encompass approximately 72 km 2 , while WFGB comprises approximately 96.3 km 2 (Figure 1). Well-developed scleractinian corals (e.g., Orbicella franksi, Psuedodiploria strigosa, and Porites astreoides) comprise 55 and 60% of the benthic cover at EFGB and WFGB, respectively, in water depths ranging from 16 to 40 m. On reef slopes, coral cover increases to as high as 94% at both banks, where the dominant species transition to O. franksi and Montastraea cavernosa (Johnston et al., 2018(Johnston et al., , 2021. Mooring buoys are present at EFGB (n = 7) and WFGB (n = 5), anchored to the seafloor with permanent stainless-steel U-bolts.
In 2015, 2016, and 2018, research at EFGB and WFGB was conducted at buoyed sites that ranged in depth from 18.6 to 19.5 m and 23.8-24.4 m, respectively (Supplementary Table 1), to survey lionfish densities, collect samples for age-growth and size analyses, and calculate removal effort (Figure 1). Four sampling efforts occurred across multiple days over the following research cruises: 31 August-3 September 2015, 29 August-1    2016, 26-28 June 2018, and 27-30 August 2018 Each research expedition constituted one cruise. Research cruises visited the same buoys when possible. If buoys were missing, similar alternate sites were visited if time allowed, but some cruises were cut short due to poor weather and/or sea conditions. Percent coral cover, coral composition, depth range, and habitat complexity are comparable for all sites. However, it is important to note buoy 1 at EFGB (i.e., EFGB1), WFGB1, and WFGB3 are located at the edge of the reef cap which transitions into deeper reef area (>30.1 m) marked by coral cover as high as 94% (Johnston et al., 2021).

Lionfish Surveys and Collections
Lionfish surveys and collections were conducted in 16-30 m water depth by two separate dive teams: scientific divers who were experts in reef fish field identification and techniques completed surveys to count lionfish, while removal divers collected lionfish on a separate, subsequent dive. Surveys were conducted using SCUBA during daylight hours between 0800 and 1600 h, with underwater visibility ranging from 80 to 100 m. Starting from the buoy U-bolt, scientific divers laid transect tapes 100 m along the reef bottom in each cardinal direction for lionfish collections by removal divers (Supplementary Figure 1). On each transect, a lionfish survey following techniques described in Green et al. (2012c) was conducted within a 50 m × 20 m area, covering a total of 4,000 m 2 per buoy. Scientific divers waited 5 min after a transect line was laid, then began a survey at the 75 m mark on the transect tape and swam toward the U-bolt in a sinusoidal pattern hovering above the transect tape, surveying under ledges and over coral heads to count lionfish. Lionfish surveys ended at the 25 m mark on the transect tape. Each lionfish survey was considered one sample.
Following the survey, lionfish were collected using pole spears and retained within a hard containment unit on two dives separated by a 2-h surface interval. One survey dive and two removal dives were conducted at each buoyed site visited. In total, 12 samples across three buoys were conducted at EFGB and WFGB in 2015 (n = 24 surveys). In 2016, 12 samples were conducted at three buoys at WFGB and 8 samples from two buoys at EFGB (n = 20 surveys). In 2018, eight samples were collected from two buoys at each bank in June (n = 16 surveys), while 12 samples were conducted at three buoys at each bank in August (n = 24 surveys). Fewer lionfish surveys and removals were completed in 2016 and June 2018 due to inclement weather.

Density
Visual survey data from transects were entered following each dive and all data underwent assurance for quality control (QA/QC) by a second party, to ensure consistency between data recorded in the field and data entered into the database.
Survey area was computed by multiplying the transect length and transect width, for a total area of 1,000 m 2 surveyed per sample. Lionfish counts recorded in each cardinal direction (i.e., four samples per buoy; Supplementary Figure 1) were divided by total area surveyed to determine density per sample (ind. m −2 ). Mean density for each buoy was calculated by averaging lionfish density across the four samples and extrapolated to density per 1-ha (ind. ha −1 ). Mean density per bank was calculated by averaging lionfish density across all samples at the bank in a sampling period and extrapolated to ind. ha −1 .
A Welch's two-sample t-test was used to contrast mean density between banks to identify if significant differences existed.

Catch per Unit Effort
Catch per unit effort (CPUE) for removal divers was standarized by restricting teams to two divers with a 45 min total bottom time for each dive. The total catch for each removal team was computed by summing the TW of all lionfish removed (kg) during one cruise. CPUE (kg diver h −1 ) was calculated by: where TW kg is the total weight captured, n is the number of divers (n = 2), and t is the total dive time (number of dives * bottom time). An ANOVA was used to determine if there were significant differences in CPUE between cruises and banks. If differences were detected, a Tukey-post hoc analysis was conducted to determine homogenious subsets.

Weight-Length Relationship
At the end of each day, lionfish total length (TL, mm) was measured on the vessel and fish were retained in labeled bags and later transported to the NOAA Galveston Lab to measure total weight (TW, g) and other metrics. A random subset of lionfish (37% of total collected) were thawed and dissected to determine sex, following methods found in Green et al. (2012b). Not all lionfish could be dissected because of time and personnel limitations. In some instances, lionfish could not be retained for additional examination. All lionfish collected in 2015 (n = 304) and 2016 (n = 376) were processed in the NOAA Galveston Lab. In 2018, lionfish exceeded the available freezer capacity on the vessel, thus some had to be discarded after TL measurements were collected. In total, 985 lionfish were collected in 2018 and 684 had to be discarded before arrival at the lab. For discarded lionfish, TW was calculated using an allometric weight-length-relationship (WLR) developed from other lionfish collected at FGBNMS. TL-TW values for lionfish measured in the NOAA Galveston Lab were log transformed to derive parameters for the WLR calcuations. An Analysis of Covariance (ANCOVA) was used to determine if these weight-length relationships differed by sex, cruise, and bank with TL as a continuous variable. A Tukey-post hoc analysis tested for the differences among categorical groups (e.g., cruise) if significant differences were detected in the ANCOVA. If differences were detected, WLR was calculated by categorical group, rather than being pooled.
Lionfish weight-length-relationship was calculated by: where W = whole body wet weight in grams (g), L = length in centimeters, and a and b are parameters (Switzer et al., 2015). A linear regression of log transformed weight and length values failed tests for heteroskedasticity (i.e., error terms were heteroskedastic). Therefore, coefficients a and b were estimated by a linear regression with robust standard erros of logarithms from the TL and TW measurements of lionfish retained each year: logW = loga + blogL, with parameters defined above, whereby a is the intercept and b is the regression coefficient (slope).

Age and Growth Parameters
Sagittal otoliths were removed from a random subset sample of lionfish in June (n = 51) and August (n = 71) 2018 from EFGB (n = 65) and WFGB (n = 57) across five size bins (0-100 mm, 101-200 mm, 201-300 mm, 301-400 mm, ≥ 401 mm), following methods described in Green et al. (2012b). Otoliths were rinsed with fresh water, dried, and stored in labeled vials. Single otoliths were placed dorsal-ventrally in the center of 22 x 22 x 20 mm embedding molds, filled with West System 105 resin and 206 hardener mixture and left to dry for 24 h (h). When possible, left otoliths were sectioned along the transverse plane using a Buehler IsoMet saw with diamond wafering blade to expose the primordium and growth increments (i.e., annuli), mounted onto glass slides using Crystal Bond, allowed to dry for 24 h, then ground and polished to the vertical mid-sagittal plane with a graded series of Buehler silicon carbide paper (600, 800, and 1,200 grit) and Buehler micropolish alumina 0.3 µm. If the left otolith was not available (e.g., damaged), the right otolith was used for age analysis. Annuli were viewed under a Leica compound microscope at 45x magnification, and sections were enumerated using the Leica LAZ EZ image analysis software. Otoliths were examined independently by two readers to determine age. Upon completion, readers compared ages and if a consensus could not be made, the otolith was discarded from further analysis. Daily increments were not enumerated in this study.
Age and length data were used to model growth by estimating parameters for the von Bertalanffy (VB) growth equation (von Bertalanffy, 1957): where L t is the length at age (t), L ∞ is the asymptotic maximum total length, K is the Brody growth coefficient or rate of growth toward L ∞ , and t 0 is the theoretical age at which a fish would be 0 mm in length. The scaling factor, t 0 , is a modeling artifact used to adjust the model for the initial size of the fish by defining the age at which the organism would be of zero length if its growth were constant over-time (Beverton, 1954;Schnute and Fournier, 1980). An ANCOVA was used to compare TL between sex and bank while controlling for age to identify if significant differences were evident. If no significant differences were found, all data were pooled. If significant differences existed, a Tukey post-hoc test was used to determine homogeneous subsets. Additionally, a Welch's two-sample t-test was used to compare mean age by sex and bank to identify if significant differences existed. Estimates from VB growth parameters from this study were compared to previous studies and other regions. All mean densities are reported in individuals per hectare (ind. ha −1 ) ± standard error (SE).

Densities
Mean lionfish densities (individual ha −1 ) were computed for each buoy at each bank by cruise (Supplementary Table 1 Table 1). The Welch's t-test showed significant differences in mean density between banks, in which greater mean lionfish densities were observed at WFGB [t (184.43) = -3.31), p < 0.001].
Mean lionfish densities at FGBNMS appear to be markedly lower than in other invaded regions of the southern GoM, Caribbean Sea, and Western Atlantic Ocean, but slighlty higher than other natural reefs in the northeastern GoM ( Table 2).

Catch per Unit Effort
A total of 1,665 lionfish were collected from FGBNMS during the four sampling efforts.

Age and Growth Parameters
Otolith pairs were removed from lionfish ranging in TL 125-444 mm and TW 18-1,153 g. Lionfish were sexed, of which 84 were male, 35 were female, and 3 were undetermined. Of the 122 otoliths collected for aging, 12 were discarded due to poor cuts or over-polished sections. Two readers separately aged otoliths, in which 51% agreement was reached. Following this, readers collaboratively aged the otolith sections and agreed on 91% of the sample. Lionfish ages ranged from 0 to 10 years, with a mean of 3.9 years. Female lionfish ranged in age from 1 to 9 (mean 3.6 ± 2.3), while males ranged from 0 to 10 years old (yo; mean 3.9 ± 2.1; Table 5). ANCOVA analysis indicated significant differences in TL between males and females [ANCOVA: F (1, 97) = 10.57, p = 0.002], but not by bank (p = 0.894; Supplementary Table 7) when controlling for age. As a result of this difference, otolithbased age data were pooled by bank to generate sex-specific von Bertalanffy (VB) growth equations (Figure 3). Age-0 lionfish were removed from the sample (n = 2) as to not influence the VB equations by misrepresenting the growth curve from lack of daily increment estimations. Male lionfish had lower Brody growth coefficients (K), higher asymptotic maximum length (L ∞ ), and reached a greater length-at-age than females (Table 5). Mean size-at-age data shows males were larger than females for every age class, except age-2 lionfish (Figure 4). For example, age-1 male lionfish were on average 266 mm, but as large as 313 mm; age-1 females were on average 213 mm and as large as 273 mm (Figure 4).

DISCUSSION
Our study is the first to provide population demographic and CPUE data for lionfish in FGBNMS, a nationally significant marine system marked by some of the healthiest coral reefs in the U.S. Lionfish were observed at both banks on all four research cruises. The results of our study revealed significant differences  Only 610 lionfish were sexed as they were randomly selected during processing (approximately one of every three lionfish), thus there is a discrepancy between total collected ( Table 3) and sexed values.
in density, size, and age of lionfish at FGBNMS, as well as CPUE, between sampling periods.

Main Findings
Lionfish densities increased fourfold at EFGB from 2015 to August 2018, but remained nearly constant at WFGB in 2015, 2016, and June 2018, followed by a twofold increase from June 2018 to August 2018. This was surprising given the similarities in habitat structure, prey availability, and nearly absent fishingmortality for both reefs. Survey effort was consistent across the banks for all years apart from 2016, in which, four additional samples were collected at WFGB than at EFGB. This, however, would not account for the difference in density trends recorded between these two banks. Lionfish larvae are likely distributed to EFGB and WFGB via the Loop Current from southern GoM reefs in the Yucután Peninsula (Kitchens et al., 2017); however, retention and distribution of lionfish larvae between the two reefs have not been investigated. High local retention of larvae from dominant coral species (e.g., P. strigosa) has been recorded for EFGB and WFGB (Davies et al., 2017;Garavelli et al., 2018), largerly influenced by eddies known to spin off from the Loop Current (Limer et al., 2020). Limer et al. (2020) noted betweenbank exchange of coral larvae, with an eastward transport, whereby WFGB contributed coral larve to EFGB. Although this is highly influenced by variation in eddy movement, it could lead to a disproportionate mid-range larval dispersion, in which EFGB receives more larval input than WFGB. This is not to suggest that a direct comparison between inter-bank connectedness, retention, and transport of coral larvae can be made for lionfish larval exchange; however, it is worth noting that WFGB could act as an additional source for lionfish larvae at EFGB. This exchange needs to be confirmed using mesocale current modeling of lionfish larvae, but would likely provide insight as to why lionfish densities continually increase at EFGB while they remain nearly constant at WFGB. Fishing mortality is negligible at both banks and the fish community and benthic habitat composition are similar, which suggests there are other factors contributing to this disparity in densities.
EFGB1, WFGB1, and WFGB3 exhibited some of the highest mean densities documented across the sampling period. All three buoys are located near the edge of the reef cap which slopes to deeper water rather rapidly (Figure 1). These slopes and deeper areas of the reef (>30.1 m) are marked by more expansive coverage of scleractinian coral, as high as 94%, offering a more continuous and complex habitat than the bank crest. Lionfish have been noted to occur in higher densities on reef ecosystems lying 30-60 m, with higher abundances associated with greater structural complexity (Lesser and Slattery, 2011;Claydon et al., 2012;Nuttall et al., 2014;Andradi-Brown, 2019). On Honduran reefs, lionfish in depths > 30 m exhibited behavior consistent with diver wariness that was observed in the shallow reef areas (<30 m), suggesting potential exposure to culling. This does not explicitly suggest lionfish travel between the shallow reef cap to deeper areas; however, is notable given the lack of culling efforts that occur below 30 m (Andradi-Brown et al., 2017). Because of the high continuous coral coverage between the shallow (16-24 m) and deep reef (30-60 m) at the three buoys, and that lionfish have been recorded to travel up to 2 km (Dahl and Patterson, 2020), it is not unrealistic to assume that lionfish could be transiting along the reef slope at these locations. It is plausible this transition to deep reef plays a role in the higher densities observed at these three buoyed areas, as has been previously recorded in Little Cayman (Bejarano et al., 2014).
Overall, mean CPUE significantly increased with each sampling period, but did not differ between banks. In particular, mean CPUE nearly doubled between June 2018 and August 2018 cruises. It is not apparant why CPUE increased so drastically between these two periods. Removal divers were thoroughly vetted for these cruises, with preference given to those with removal experience and those who had been to FGBNMS before. We do not suspect this difference in CPUE can be attributed to differences in diver ability. CPUE and lionfish density increased overtime concurrently, suggesting the lionfish population is expanding at EFGB and WFGB. We would expect that an increase in effort would reduce lionfish densities if fishing mortality was  having a measurable impact on the population; however, this was not evident in our study.
We note, in general, lionfish shifted toward larger size classes at EFGB and WFGB during our study period. This is contradictory to previous studies in the northeastern GoM, in which lionfish mean size appeared to shift to smaller size classes (Dahl et al., 2019). Targeted fishing efforts can have a pronounced effect on the size structure of marine fish species, with a general shift toward smaller size classes with increasing fishing pressure (e.g., Bianchi, 2000;Tu et al., 2018). Lionfish are predominantly removed through spearfishing as they have low vulnerability to conventional hook-and-line techniques (Morris and Whitfield, 2009). Spearfishing is a restricted activity in FGBNMS, thus, fishing-mortality is nearly negligible at EFGB and WFGB and likely does not effect lionfish size class variation. Lionfish at EFGB and WFGB showed sexually dimorphic size variation consistent with previous research (Morris and Whitfield, 2009;Morris et al., 2011;Edwards et al., 2014;Dahl et al., 2019;Eddy et al., 2019;Fogg et al., 2019), where males were significanlty larger than females and comprised a greater percentage of the larger size classes across all sampling periods. Evidence of adult mediated density-dependent growth has been reported for lionfish on artificial reefs in the northeastern GoM (Dahl et al., 2019); however, we do not find evidence to support this phenomenon at FGBNMS. FGBNMS is absent of the ecological (e.g., warm temperatures, low resource availability) and biological (e.g., high density, disease) constraints described to influence densitydependent traits in adult lionfish Dahl et al. (2019). If lionfish densities at EFGB and WFGB increase to carrying capacity, density-dependence may affect the population structure in the future; however, this would need to be monitored to determine if this occurred.
Results of this study revealed significant differences in weightlength relationships between sex and cruise, but not by bank. Males were heavier per length than females for all sampling periods, consistent with results from other regions in the northern GoM (Dahl and Patterson, 2014;Fogg et al., 2019). Mean size was significantly different between cruises, generally increasing with each year. Similarly, Dahl and Patterson (2014) reported significant inter-annual increases in mean size for lionfish in the northeastern GoM on artificial and natural reefs. Lionfish weight-length relationships did not differ by bank; EFGB and WFGB are comparable in habitat and food resources which may allow for similar body conditions at the two reefs. Lionfish populations likely have not reached their peak at FGBNMS, as density and mean size increased throughout the study period.
Mean age of lionfish were significantly different between banks, but not by sex. Lionfish were older at WFGB than EFGB by approximately 1-year and WFGB was marked by more, older individuals. Males constituted more of the lionfish aged at WFGB than EFGB (78% vs. 54%) and likely contributes to the observed differences. Female lionfish tend to grow faster than males, but attain a smaller maximum size and length at maturity throughout their invaded range (Edwards et al., 2014;Dahl et al., 2019;Eddy et al., 2019;Fogg et al., 2019), which was also true for this study. A disproportion in larval exchange could contribute to the difference in lionfish age between the two banks, if WFGB in fact acts as a larval source for EFGB. If this holds true, we would expect to see more, older individuals at WFGB and a predominatly younger population at EFGB. Lionfish from FGBNMS were notably comprised of more, older individuals which may be attributed, in part, to fewer collections of smaller lionfish. Small lionfish are harder to detect in complex natural habitats (Harris et al., 2019); EFGB and WFGB are notably complex reef habitats that could provide refuge for smaller lionfish, rendering their detection and removal challenging. Lionfish at FGBNMS exhibited one of the slowest growth rates than in any of the other invaded regions. Daily increment analyses were not conducted in this study, and thus, removes the effect of juvenile growth characteristics in the VB growth curve estimates, which may lead to an underestimation of these values (Pusack et al., 2016;Dahl et al., 2019;Eddy et al., 2019;Fogg et al., 2019). Additionally, the sample size for the age analysis in this study was limited and may further skew the VB growth curve estimates. Future studies should prioritize collection of smaller lionfish to include juveniles and increase the sample size in age-based analyses to better quantify these parameters for FGBNMS.
Lionfish age structure at EFGB and WFGB most closely matched the age structure of lionfish collected in Bermuda , which is interesting given the difference in invasion chronology. However, we note FGBNMS and Bermuda have similar habitat types and are both northern latitude coral reefs with seasonal temperature fluctuations, which may contribute to the observed similarities. The first reported observation of lionfish in Bermuda occurred in 2000, while the age and size analysis were on lionfish collected between 2012 and 2016 . Here, we report age structure for lionfish from FGBNMS in which observations were first reported in 2011 and analysis conducted in 2018. Results from this study suggest lionfish invaded EFGB and WFGB 3 years prior to initials observations , but also shows lionfish can successfully attain older ages in their invaded region when fishing-induced mortality and predation are negligible. Lionfish grow considerably slower at EFGB and WFGB than Bermuda for pooled data and males, exhibiting two-fold lower Brody growth coefficients. However, male lionfish from FGBNMS attain similar asymptotic maximum total length as males in Bermuda (415 mm vs. 418 mm, respectively). We report a larger L ∞ for females at FGBNMS as compared to Bermuda (345 mm vs. 332 mm; Eddy et al., 2019), though this is likely a result of under-sampling for the age analysis. It is not obvious from our analysis whether FGBNMS has a greater proportion of larger, older lionfish than other invaded regions, or if this is an artifact of a small sample size. A more robust sample size would afford for more confidence in interpreting growth parameter values for lionfish from FGBNMS.
It should also be noted that it is not clear whether this differentiation between regional age and size structure is due to differences in species, in that, both species of lionfish (P. volitans, P. miles) have been confirmed in the Caribbean Sea and Western Atlantic Ocean, while only one species (P. volitans) has been confirmed in the GoM (Betancur-R et al., 2011;. There are conflicting reports on whether the distribution and expansion of the two species of lionfish have been homogeneous throughout their invasion of these regions (Perez-Portela et al., 2018) or if a genetic founder effect occurred (Betancur-R et al., 2011;Toledo-Hernández, 2014;. Nevertheless, genetic analyses may afford the ability to determine if differences reported in age/size and growth demographics are species-specific, and thus should be included in future research. Lionfish life history parameters characterized in this study have illustrated a successful invasion at both EFGB and WFGB in the northwestern GoM. When compared with lionfish population structure in other invaded reefs in the GoM, Western Atlantic Ocean, and Caribbean Sea, striking differences are evident. Seven years after the invasion, lionfish densities at EFGB and WFGB remain twofold lower than in Mexico and Belize (Hackerott et al., 2013;Sosa-Cordero et al., 2013), two to sevenfold lower than in the Bahamas and North Carolina (Morris and Whitfield, 2009;Darling et al., 2011;Green et al., 2012b), twofold lower than in the southern Caribbean Sea, and 10-fold lower than on artificial reefs in the northeastern GoM (Dahl and Patterson, 2014;Dahl et al., 2019). The lionfish invasion chronology was similar among these regions, in terms of time between first reported observation and date of survey, which suggests a slower colonization and growth rate of lionfish at EFGB and WFGB. Lionfish in higher latitudes (i.e., Bermuda, northern GoM) have shown seasonal spawning patterns, with peak spawning activity in the summer and fall and a subsequent decline or absence in the winter months Eddy et al., 2019). This seasonality in spawning could result in seasonal fluctuations in density. However, the limited surveys, predominantly conducted the same time of year, may yield skewed density estimates at EFGB and WFGB. More surveys that span a greater temporal scale would be necessary to determine if density exhibits a seasonal pattern at FGBNMS, particularly one that correlates with spawning seasonality.
Lionfish were smaller on average at FGBNMS (mean TL 281.9 mm ± 53.2) than those collected in Bermuda (mean TL 339 mm ± 1.6; Eddy et al., 2019) and North Carolina (mean TL 305 mm; Whitfield et al., 2007). The differences in the size structure between FGBNMS and the Western Atlantic Ocean may be attributed to different environmental factors (e.g., temperature) or regional differences in prey availability, sampling techniques (i.e. fewer smaller, younger lionfish sampled in the Atlantic), or that marine teleosts are typically larger in temperate regions (Luckhurst et al., 1992;Robertson et al., 2005).
The size classes of lionfish at EFGB and WFGB showed a continual shift toward having more, larger fish, particularly with respect to males. This differs from reports in the northeastern GoM, which showed lionfish shifting to smaller size classes (Dahl et al., 2019). Density-dependent growth processes may have acted as an important population-level regulator of lionfish on reefs in the northeastern GoM (Dahl et al., 2019); however, it is not apparent that this is occurring with lionfish populations at FGBNMS. Ecological conditions (e.g., rugous, complex reef structure, high species richness and prey biomass) are much different at EFGB and WFGB than on artificial reefs (Rooker et al., 1997;Dance et al., 2011;Patterson et al., 2014;Dahl et al., 2016Dahl et al., , 2019Johnston et al., 2021), for example, and we suspect that lionfish are not as susceptible to competition for food and habitat resources that would lead to density-dependent regulation of adults in this system. This suggests there are likely other mechanisms controlling lionfish populations at FGBNMS.

Invasion Impacts and Implications for Resilience and Management
Lionfish are subject to negligible fishing-induced mortality, lack observations of predation by native predators, do not exhibit traits of population-level regulation (i.e., density-dependent growth), and lacked signs of disease at FGBNMS. Lionfish are larger in the invaded range as compared to their native range in the Indo-Pacific (Pusack et al., 2016), which may lessen predation by native predators. Additionally, native predators may be wary of consuming lionfish because of the venomous dorsal spines that are uncharacteristic of other native prey species. Densitydependent population control is not evident at FGBNMS, as lionfish exhibited lower densities and were characterized by larger, older fish. Fewer or shorter recruitment events may, in part, be responsible for lower lionfish densities at EFGB and WFGB affected by seasonal spawning. Lionfish showed seasonal spawning in northern latitudes where ambient water temperature fluctuations followed a seasonal pattern Eddy et al., 2019). EFGB and WFGB also reveal seasonal variation in water temperature, with a 10 • C difference between maximum and minimum temperatures in the summer and winter months (Johnston et al., 2021). Although it cannot be deduced from data presented in this study if seasonal spawning or intermittent recruitment are leading to lower densities, given other similarites between lionfish in FGBNMS and Bermuda and that FGBNMS is a remote reef system, it could be assumed these are individually or collectively contributing factors. Future research would need to investigate this phenomenon to determine if recruitment or spawning is reduced at these banks. An ulcerative disease led to considerable declines in northern GoM lionfish populations in 2018 (Harris et al., 2020); however, evidence of this disease was not observed during surveys/removals of lionfish at FGBNMS during this study period. Although this may have impacted lionfish populations at EFGB and WFGB at some point, it was not a contributor to observations for this study. The population characterisitcs at these reefs may also indicate there are other environmental or ecological controls contributing to the observed differences between FGBNMS and other invaded regions.
Our study may support the idea that lionfish have been incorporated functionally into the existing ecosystem processes at EFGB and WFGB (e.g., lower densities, slower growth, older invidivuals) because few alterations, if any, have been observed to the extant ecosystem dynamics of these reefs (Johnston et al., 2016a(Johnston et al., , 2021. This could be attributed in part, to ecological resilience, in that the overall reef health improves the capacity to absorb the lionfish invasion with very little impact. Ecological resilience is the capacity of a system to withstand change while maintaining processess and structures (Holling, 1973), and offers a theoretical frame for understanding ecosystem-level controls for biological invasions. For example, EFGB and WFGB possess slightly different reef fish assemblages than other reef systems in the Caribbean Sea and GoM, maintaining stable resident populations dominated by pomacentrids and labrids (Wetmore et al., 2020;Johnston et al., 2021), with a limited presence of important predators such as lutjanids and haemulids (Rooker et al., 1997;Johnston et al., 2021). Lionfish are a known generalist meso-predator with a diet consisting primarily of small-bodied teleosts and crustaceans (Morris and Akins, 2009;Cote et al., 2013b;Eddy et al., 2016). Thus far, the presence of lionfish at EFGB and WFGB has not correlated with any measurable decrease or negative effect to biomass of prey fishes (Johnston et al., 2016a), as has been reported to have rapidly occurred in other invaded regions (e.g., Albins and Hixon, 2008;Green et al., 2012a;Benkwitt, 2014;Ballew et al., 2016), even with a diet that is broadly similar in composition and quantity (Johnston, unpublished data). This may suggest lionfish provide the same functional role as other important predators at EFGB and WFGB without reducing the ecological funtion of prey species that could otherwise lead to ecosystem-level impacts (Albins and Hixon, 2013). Despite the addition and increasing densities of lionfish at EFGB and WFGB, the fundamental ecological functions provided by reef fish seems to be conserved across the ecosystem.
EFGB and WFGB have shown evidence of ecological resilience to perturbations for several decades, being one of the healthiest coral reefs in the GoM and Caribbean regions (Gil-Agudelo et al., 2020;Johnston et al., 2021). In the early 1990's, a phase shift in coral reefs was described referring to a drastic change from a coral to algal dominance, particularly with coral reefs throughout the Caribbean Sea (Done, 1992;Hughes, 1994;Mumby et al., 2013). However, this has not occurred at EFGB and WFGB, rather it has remained a coral dominated reef over the last three decades even with significant increases in macroalgae cover, again suggesting resilience to disturbance (Schmahl et al., 2008;Johnston et al., 2021). Additionally, EFGB and WFGB appear to be resilient, in part, in their apparent recovery, with very minimal mortality, from major bleaching events , albeit the significant warming water temperature trend documented over the last 30 s (Johnston et al., 2021). Given the reefs' ability to respond favorably following extreme and longwithstanding disturbances, it may be the same has been true for the lionfish invasion. It is not clear whether the lionfish invasion will have negative impacts in the future, but thus far, their presence has not contributed to apparent degredation of the reef systems or ecosystem function (Johnston et al., 2016a(Johnston et al., , 2021. The presence and/or level of ecological resilience the lionfish population characteristcs confer cannot be deduced from our study; however, it may prove itself overtime with longer temporal observations. Other regions have revealed the lionfish invasion did not lead to signifcant negative impacts to prey fish communities, and though a definitive explanation was unclear, the health of reef systems and ecological parameters were suggested as potential contributors (Valdivia et al., 2014;Elise et al., 2015;Hackerott et al., 2017;Goodbody-Gringley et al., 2019). It cannot be ignored that ecological resilience of EFGB and WFGB may be a mechanism for controlling lionfish at these reefs, whereby managing the populations lie with the environment itself rather than through human intervention alone. However, it must also be understood that more time, information, and data are needed to evaluate whether or not ecological resilience is evident at FGBNMS.
Our present study of life history characteristics and removal effort of lionfish at EFGB and WFGB may suggest biological (e.g., seasonal recruitment) and/or ecological (e.g., resilience) mechanisms are contributing to controlling the invasion at this protected reef. This has the potential for very different management of lionfish at these reefs, because to date, most strategies have been reactionary and focused on control through mitigation and removal. Here, we argue, that gaining an understanding of the lionfish invasion, seasonality in recruitment and spawning, and population characteristics as it relates to ecological resilience will allow for more deliberate engagement with resilience-based approaches to governing (Cook et al., 2010;Foxcroft and McGeoch, 2011;Garmestani and Benson, 2013) their presence in reef ecosystems. In this sense, we argue lionfish management may be more successful if governance takes a dynamic approach that considers ecological systems jointly to achieve a strategy that accommodates adaptation. As such, we recommend resource managers employ resilience-based management practices that "use knowledge of current and future drivers influencing ecosystem function to prioritize, implement, and adapt management actions that sustain ecosystems and human wellbeing" (McLeod et al., 2019) for established invasive species. The data provided in this study are an example of empirical-based metrics that can be used to identify desirable characteristics that support resilience in a reef system, which can later be integrated into mechanistic ecological models to understand the processes of resilience and/or predict future coral reef trajectories under multiple stressors. We suggest FGBNMS is uniquely positioned to evaluate ecological resilience and the level of control it can provide after an invasive marine fish species establishes in the ecosystem. Furthermore, we propose resource managers at FGBNMS, and other invaded areas with similar characteristics (e.g., remote offshore location, healthy system, few anthropogenic disturbances), redirect their approach to management of lionfish (and other invasive species) to focus on understanding ecosystem dynamics as opposed to controlling a single species (e.g., Angeler et al., 2015). Additionally, these efforts would also be useful in protected marine areas where consistent, physical removal may not be an effective or practical mitigation strategy. By gaining an understanding of biological invasions through ecosystem function, governance approaches can be resiliencebased and coordinate management at scales relevant to ecosystem functions, rather than anthropocentric scales driven primarily by political and jurisdictional boundaries (Chaffin et al., 2016). It will be important that future studies prioritize ecosystem-based analyses that tailor to the understanding of ecosystem-scale functions that may be affected by the presence of lionfish. This is particularly critical for locations, like FGBNMS, that have minimal to no influence of fishery-induced mortality and will require other mechanisms of populationlevel control. With ecosystem-based analyses, resource managers and scientists may be able to better understand whether lionfish are having a negative impact to an invaded area or if in fact, the environment can transition into another functioning regime, whereby ecosystem resilience is the primary control mechanism.

CONCLUSION
The extent to which invasive lionfish may impact the unique offshore ecosystem of FGBNMS is relatively unknown, though no evidence suggests, to date, their presence has had any measurable impact on biodiversity, community structure, or ecosystem function (Johnston et al., 2016a(Johnston et al., , 2018(Johnston et al., , 2021. This study describes key life history attributes and removal effort of lionfish on a protected coral reef system in the northwestern GoM. Despite minimal human removal of lionfish at FGBNMS, lionfish populations have shown slower population expansion as compared to other invaded regions of the Western Atlantic Ocean, Caribbean Sea, and northeastern GoM. In particular, this study demonstrated lionfish at FGBNMS exhibit lower densities with significantly increasing CPUE, indicating population expansion. Additionally, a temporal shift toward larger size classes suggests an expanding population. This study showed longevity of lionfish with more, older individuals and slower growth rates comparable to those aged in other temperate regions (i.e., North Carolina and Bermuda). EFGB showed an increasing lionfish density, while WFGB remained relatively stable over the study period, likely attributed to a slight unidirectional larval dispersal between the two banks. EFGB had predominately younger lionfish as compared to WFGB, further suggesting this unidirectional transfer of larvae. Given the absence of fishing mortality (because spearfishing is a restricted activity at FGBNMS), expanding population of lionfish over a decade (introduced in 2008, study occurred in 2015, 2016, and 2018), lack of observed predation from natural predators, absence of disease, and that lionfish did not exhibit traits of densitydependent growth (Dahl et al., 2019), our study suggests there may be an environmental or ecological control for the lionfish populations at FGBNMS. Lionfish densities may be lower due to seasonal spawning or intermittent recruitment which has impacted the expansion of lionfish populations in northern latitudes (i.e., Bermuda and northern GoM), mostly attributed to lower ambient water temperatures, of which are experienced seasonally at FGBNMS. In contrast, we suggest EFGB and WFGB may present ecological resilience to invasive lionfish, that affords these reef systems the ability to absorb impacts on the established populations. This phenomenon would have very different implications for management at this protected reef, where consistent, physical removal is not practical. Ecosystem-based analyses would afford a better understanding of whether lionfish are having a negative impact at FGBNMS and require directed intervention, or if the environment can act as the primary control mechanism for this reef system.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

ETHICS STATEMENT
The animal study was reviewed and approved by the Texas A&M University's Animal Welfare Assurance Program Institutional Animal Care and Use Committee.

AUTHOR CONTRIBUTIONS
RB, AF, MJ, JR, and GJ contributed to the conceptualization of this study. RB and GJ designed the experiments. RB and AF analyzed the data. JR contributed the reagents, materials, and analysis equipment. RB wrote the initial draft of the manuscript. All authors contributed to the article, revised the manuscript, and approved the submitted version.

FUNDING
Funding for sampling efforts were funded, in part, by Texas State Aquarium, Ripley's Aquarium, the Georgia Aquarium, Texas Lionfish Control Unit, Reef Environmental Education Foundation, and the M/V Fling.

ACKNOWLEDGMENTS
We would like to thank the Texas State Aquarium, Ripley's Aquarium, Texas Lionfish Control Unit, Reef Environmental Education Foundation, and the M/V Fling for contributing funds for these research expeditions. We also greatly appreciate all of the volunteer divers that participated, as well as the crew members of the M/V Fling. We appreciate all of the volunteers that helped with lionfish dissections and measurements, and we appreciate graduate and undergraduate students in JR.'s lab. Lionfish collections were a permitted activity in FGBNMS, under permit numbers FGBNMS-2015-001 (2015), FGBNMS-2016-006 (2016), and FGBNMS-2018-002 (2018. We appreciate FGBNMS granting the permits for these expeditions. No Animal Use Protocol (AUP) was required for this study, as determined by Texas A&M University's Animal Welfare Assurance Program Institutional Animal Care and Use Committee.