Density dependence in total egg production per spawner for marine fish

A paradigm of fisheries science holds that spawning stock biomass (SSB) is directly proportional to total egg production (TEP) of fish stocks. This “SSB–TEP proportionality” paradigm has been a basic premise underlying the spawner–recruitment models for fisheries management and numerous studies on recruitment mechanisms of fish. Studies on maternal effects on reproductive potential of a stock have progressed during the last few decades, leading to doubt concerning the paradigm. Nonetheless, a direct test of the paradigm at multidecadal scales has been difficult because of data limitations in the stock assessment systems worldwide. Here, we tested the paradigm for marine fish based on a novel combination of two independent 38‐year time series: fishery‐dependent stock assessment data and fishery‐independent egg survey data. Through this approach, we show that the SSB–TEP proportionality is distorted by density dependence in total egg production per spawner individual (TEPPS) or spawner unit weight (TEPPSW) at a multidecadal scale. The TEPPS/TEPPSW exponentially declined with biomass and thus was density‐dependent for Japanese sardine, a small pelagic species exhibiting a high level of population fluctuation, in the western North Pacific. By contrast, the TEPPS/TEPPSW was sardine‐density‐dependent for Japanese anchovy, another small pelagic species exhibiting a moderate level of population fluctuation well‐known for being out of phase with sardine. Our analysis revealed intraspecific (sardine) and interspecific (anchovy) density dependence in TEPPS/TEPPSW, which was previously unaccounted for in spawner–recruitment relationships. Such density‐dependent effects at the time of spawning should be considered in fisheries management and studies on recruitment mechanisms.     

Influence of larval transport and temperature on recruitment dynamics of North Sea cod (Gadus morhua) across spatial scales of observation

The survival of fish eggs and larvae, and therefore recruitment success, can be critically affected by transport in ocean currents. Combining a model of early‐life stage dispersal with statistical stock–recruitment models, we investigated the role of larval transport for recruitment variability across spatial scales for the population complex of North Sea cod (Gadus morhua). By using a coupled physical–biological model, we estimated the egg and larval transport over a 44‐year period. The oceanographic component of the model, capable of capturing the interannual variability of temperature and ocean current patterns, was coupled to the biological component, an individual‐based model (IBM) that simulated the cod eggs and larvae development and mortality. This study proposes a novel method to account for larval transport and success in stock–recruitment models: weighting the spawning stock biomass by retention rate and, in the case of multiple populations, their connectivity. Our method provides an estimate of the stock biomass contributing to recruitment and the effect of larval transport on recruitment variability. Our results indicate an effect, albeit small, in some populations at the local level. Including transport anomaly as an environmental covariate in traditional stock–recruitment models in turn captures recruitment variability at larger scales. Our study aims to quantify the role of larval transport for recruitment across spatial scales, and disentangle the roles of temperature and larval transport on effective connectivity between populations, thus informing about the potential impacts of climate change on the cod population structure in the North Sea.     

Selection on the timing of migration and breeding: A neglected aspect of fishing‐induced evolution and trait change

Fishing can drive changes in important phenotypic traits through plastic and evolutionary pathways. Size‐selective harvest is a primary driver of such trait change, has received much attention in the literature and is now commonly considered in fisheries management. The potential for selection on behavioural traits has received less study, but mounting evidence suggests that aggression, foraging behaviour and linked traits can also be affected by fishing. An important phenomenon that has received much less attention is selection on reproductive phenology (i.e., the timing of breeding). The potential for this type of “temporal selection” is widespread because there is often substantial variability in reproductive phenology within fish populations, and fisheries management strategies or fishermen's behaviours can cause fishing effort to vary greatly over time. For example, seasonal closures may expose only early or late breeding individuals to harvest as observed in a range of marine and freshwater fisheries. Such selection may induce evolutionary responses in phenological traits, but can also have demographic impacts such as shortened breeding seasons and reduced phenotypic diversity. These changes can in turn influence productivity, reduce the efficacy of management, exacerbate ongoing climate‐driven changes in phenology and reduce resilience to environmental change. In this essay, we describe how fisheries management can cause temporal variability in harvest, and describe the types of selection on temporal traits that can result. We then summarize the likely biological consequences of temporally selective fishing on populations and population complexes and conclude by identifying areas for future research.     

Larval dispersal and population connectivity among a network of marine reserves

A major challenge in marine ecology is describing patterns of larval dispersal and population connectivity, as well as their underlying processes. We used a biophysical model to simulate dispersal of eastern oyster, Crassostrea virginica, larvae and connectivity among a network of 10 no‐take reserves in a shallow, wind‐driven estuary to assess the relative importance of spawning location, spawning date, larval behavior, larval mortality, and adult reproductive output to predicted dispersal and connectivity patterns. The location (i.e., natal reserve) and date of spawning relative to physical processes, particularly frequency of wind reversals, were the dominant drivers of dispersal and connectivity patterns. To a lesser extent, larval behavior (i.e., 3D vertical advection and ontogenetic depth regulation) and mortality modified dispersal and connectivity, whereas spatiotemporal variability in adult reproductive output was of minimal importance. Over a 21‐day larval duration, mean dispersal distance of passive surface particles ranged from 5 to 40 km. Reserves were too small (1 km²) relative to mean dispersal distances to promote extensive local retention (median 2%) and spaced too far apart (typically ~50 km) to promote extensive inter‐reserve connectivity (median 2%). Limited connectivity and local retention may preclude the network from being self‐sustainable, thereby limiting its long‐term conservation and management benefits. In reserve systems characterized by limited connectivity, management efforts should focus on increasing connectivity by increasing the number or size of reserves to realize the benefits of improved adult demographics within reserves.     

Do albacore exert top‐down pressure on northern anchovy? Estimating anchovy mortality as a result of predation by juvenile north pacific albacore in the California current system

Quantifying the mortality of marine fishes is important for understanding spawner–recruit relationships, predicting year‐class strength, and improving fishery stock assessment models. There is increasing evidence that pelagic predators can exert a top‐down influence on prey, especially during critical early life‐history stages. The objective of this study was to quantify predation by North Pacific albacore on Northern anchovy in the California current system (CCS). I estimated the abundance of juvenile albacore in the CCS from 1966–2005 using stock assessment models and spatially explicit catch‐per‐unit‐effort time series. Anchovy abundance (1966–93), both recruitment and total biomass, was obtained from a stock assessment model. Annual rates of anchovy consumption by albacore were calculated using diet studies of albacore in the CCS, an age‐structured bioenergetics model, and regional estimates of albacore abundance. The range of estimates was large: albacore may remove from less than 1% to over 17% of anchovy pre‐recruitment biomass annually. Relationships between predation and recruitment biomass were consistent with expectations from top‐down effects, but further study is required. This is the first attempt to quantify a specific source of mortality on anchovy recruits and to demonstrate potential top‐down effects of predation on anchovy.     

Where do you think you’re going? Accounting for ontogenetic and climate‐induced movement in spatially stratified integrated population assessment models

Understanding spatial population structure and biocomplexity is critical for determining a species’ resilience to environmental and anthropogenic perturbations. However, integrated population models (IPMs) used to develop management advice for harvested populations have been slow to incorporate spatial dynamics. Therefore, limited research has been devoted to understanding the reliability of movement parameter estimation in spatial population models, especially for spatially dynamic marine fish populations. We implemented a spatial simulation–estimation framework that emulated a generic marine fish metapopulation to explore the impact of ontogenetic movement and climate‐induced distributional shifts between two populations. The robustness of spatially stratified IPMs was explored across a range of movement parametrizations, including ignoring connectivity or estimating movement with various levels of complexity. Ignoring connectivity was detrimental to accurate estimation of population‐specific biomass, while implementing spatial IPMs with intermediate levels of complexity (e.g. estimating movement in two‐year and two‐age blocks) performed best when no a priori information about underlying movement was available. One‐way distributional shifts mimicking climate‐induced poleward migrations presented the greatest estimation difficulties, but the incorporation of auxiliary information on connectivity (e.g. tag‐recapture data) reduced bias. The continued development of spatially stratified modelling approaches should allow harvested resources to be better utilized without increased risk. Additionally, expanded collection and incorporation of unique spatially explicit data will enhance the robustness of IPMs in the future.     

Is spawning stock biomass a robust proxy for reproductive potential?

Reference points based on fishing mortality (F) and spawning stock biomass (SSB) are a requirement of many fisheries management frameworks. SSB is assumed to be a proxy for stock reproductive potential (SRP). Limit reference points based on SSB are used to indicate the level of biomass below which productivity is affected. SSB fails to account for changes in fecundity, egg viability and sex ratio, and it has been argued that total egg production (TEP) provides a better reflection of SRP. We explore how accounting for TEP impacts limit reference points and evidence for a relationship between stock and recruit. Time series of SSB and TEP are compared for three North Sea stocks: cod (Gadus morhua), herring (Clupea harengus) and plaice (Pleuronectes platessa). Dynamics based on TEP are different from those based on SSB for cod and plaice, but the stock–recruit relationships were not ‘improved’ using TEP. Shifts in productivity (spawner per recruit) occur in all three time series and SSB underestimated uncertainty. Yet again, it was shown that assumptions of stationarity about fish population productivity are incorrect. We argue that the use of TEP does improve the realism in our understanding of stock dynamics, and demographically, more complex management strategy evaluation is required to develop management procedures that are robust to uncertainty and integrate F and the demographic health of a stock. Empirical feedback control systems based on fisheries independent indices including surveys of eggs, larvae, recruits, juveniles or spawning adults should be evaluated and compared to traditional approaches.     

Oceanographic drivers of sablefish recruitment in the California Current

Oceanographic processes and ecological interactions can strongly influence recruitment success in marine fishes. Here, we develop an environmental index of sablefish recruitment with the goal of elucidating recruitment‐environment relationships and informing stock assessment. We start with a conceptual life‐history model for sablefish Anoplopoma fimbria on the US west coast to generate stage‐ and spatio‐temporally‐specific hypotheses regarding the oceanographic and biological variables likely influencing sablefish recruitment. Our model includes seven stages from pre‐spawn female condition through benthic recruitment (age‐0 fish) for the northern portion of the west coast U.S. sablefish stock (40°N–50°N). We then fit linear models and use model comparison to select predictors. We use residuals from the stock‐recruitment relationship in the 2015 sablefish assessment as the dependent variable (thus removing the effect of spawning stock biomass). Predictor variables were drawn primarily from ROMS model outputs for the California Current Ecosystem. We also include indices of prey and predator abundance and freshwater input. Five variables explained 57% of the variation in recruitment not accounted for by the stock‐recruitment relationship in the sablefish assessment. Recruitment deviations were positively correlated with (i) colder conditions during the spawner preconditioning period, (ii) warmer water temperatures during the egg stage, (iii) stronger cross‐shelf transport to near‐shore nursery habitats during the egg stage, (iv) stronger long‐shore transport to the north during the yolk‐sac stage, and (v) cold surface water temperatures during the larval stage. This result suggests that multiple mechanisms likely affect sablefish recruitment at different points in their life history.     

Interannual variation in pelagic juvenile rockfish (Sebastes spp.) abundance – going with the flow

We report results from 28 yr of a midwater trawl survey of pelagic juvenile rockfish (Sebastes spp.) conducted off the central California coast. The fishery‐independent survey is designed to provide pre‐recruit indices of abundance for use in groundfish stock assessments. Standardized catch rate time series for 10 species were developed from delta‐generalized linear models that include main effects for year, station, and calendar date. Results show that interannual fluctuations of all 10 species are strongly coherent but highly variable, demonstrating both high‐ and low‐frequency components. A similarly coherent result is observed in the size composition of fish, with large fish associated with elevated catch rates. In contrast, spatial and seasonal patterns of abundance show greater species‐specific differences. A comparison of the shared common trend in pelagic juvenile rockfish abundance, derived from principal components analysis, with recruitments from five rockfish stock assessments shows that the time series are significantly correlated. An examination of oceanographic factors associated with year‐to‐year variability indicates that a signature of upwelled water at the time of the survey is only weakly related to abundance. Likewise, basin‐scale indices (the Multivariate El Niño‐Southern Oscillation Index, the Pacific Decadal Oscillation, the North Pacific Gyre Oscillation, and the Northern Oscillation Index) are poorly correlated with abundance. In contrast, sea level anomalies in the months preceding the survey are well correlated with reproductive success. In particular, equatorward anomalies in the alongshore flow field following the spawning season are associated with elevated survival and poleward anomalies with poor survival.     

Trends and management implications of human‐influenced life‐history changes in marine ectotherms

Evidence is accumulating that many marine ectotherms are undergoing rapid changes in their life‐history characteristics. These changes have been variously attributed to fisheries‐induced evolution, inhibited adult growth rate due to oxygen limitation at higher temperatures, and plastic responses to density dependence or changes in ocean productivity. Here, we review the diverse underlying mechanisms by which plastic and evolutionary responses to climate change and fisheries are likely to produce similar life‐history trends in harvested marine ectotherms, leading to faster life‐histories with earlier maturation and smaller adult size‐at‐age. While mechanistically understanding these growth and maturation changes may be difficult, it is becoming clear that changing life‐histories will lead to modified population dynamics, productivity and natural mortality of the affected species. We discuss how the observed and expected life‐history changes could affect the assumptions and uncertainty within single and multispecies models currently used in marine ecosystem management, highlighting that models which allow for dynamic life‐history traits often report significantly different estimates of stock biomass. Given that both climate‐ and harvest‐induced life‐history changes are likely to intensify and possibly amplify each other, there is an urgent need to adequately assess the implications of faster life‐histories for marine ecosystem management. This is especially true for data‐poor stocks, where growth and maturation are not regularly assessed. Targeted monitoring can be used to inform responsive management, but for improved sustainability outcomes, a precautionary approach to management that is robust to life‐history trends is advised.     

Historical spatial baselines in conservation and management of marine resources

Increased knowledge on the spatial distribution of marine resources is crucial for the implementation of a true ecosystem approach to management and the conservation of marine organisms. For exploited fish species characterized by aggregation behaviour during spawning time, the identification and tracking of spawning areas is essential for a correct assessment of their productivity and population abundance. To elucidate this concept, we reconstructed the spatio‐temporal distribution of adult plaice (Pleuronectes platessa, Pleuronectidae) during spawning time along the 20^th century. Historical data reveal that not only the abundance but also the former population richness was much higher than previously estimated and has declined because of protracted over‐exploitation during the last 30 years. We conclude that forecast of stock recovery to former levels of abundance neglecting spatial reorganizations might be over‐optimistic and shaded by a lost memory of the past population richness. These results reinforce the importance of managing exploited marine resources at a greater spatial resolution than has been carried out in the history of fishery management.     

Compensatory density dependence in fish populations: importance, controversy, understanding and prognosis

Density‐dependent processes such as growth, survival, reproduction and movement are compensatory if their rates change in response to variation in population density (or numbers) such that they result in a slowed population growth rate at high densities and promote a numerical increase of the population at low densities. Compensatory density dependence is important to fisheries management because it operates to offset the losses of individuals. While the concept of compensation is straightforward, it remains one of the most controversial issues in population dynamics. The difficulties arise when going from general concepts to specific populations. Compensation is usually quantified using some combination of spawner–recruit analysis, long‐term field monitoring or manipulative studies, and computer modelling. Problems arise because there are limitations to each of these approaches, and these limitations generally originate from the high uncertainty associated with field measurements. We offer a hierarchical approach to predicting and understanding compensation that ranges from the very general, using basic life‐history theory, to the highly site‐specific, using detailed population models. We analyse a spawner–recruit database to test the predictions about compensation and compensatory reserve that derive from a three‐endpoint life‐history framework designed for fish. We then summarise field examples of density dependence in specific processes. Selected long‐term field monitoring studies, manipulative studies and computer modelling examples are then highlighted that illustrate how density‐dependent processes led to compensatory responses at the population level. Some theoretical and empirical advances that offer hope for progress in the future on the compensation issue are discussed. We advocate an approach to compensation that involves process‐level understanding of the underlying mechanisms, life‐history theory, careful analysis of field data, and matrix and individual‐based modelling. There will always be debate if the quantification of compensation does not include some degree of understanding of the underlying mechanisms.     

A metapopulation perspective for salmon and other anadromous fish

Salmonids are an important component of biodiversity, culture and economy in several regions, particularly the North Pacific Rim. Given this importance, they have been intensively studied for about a century, and the pioneering scientists recognized the critical link between population structure and conservation. Spatial structure is indeed of prime importance for salmon conservation and management. At first glance, the essence of the metapopulation concept, i.e. a population of populations, widely used on other organisms like butterflies, seems to be particularly relevant to salmon, and more generally to anadromous fish. Nevertheless, the concept is rarely used, and barely tested. Here, we present a metapopulation perspective for anadromous fish, assessing in terms of processes rather than of patterns the set of necessary conditions for metapopulation dynamics to exist. Salmon, and particularly sockeye salmon in Alaska, are used as an illustrative case study. A review of life history traits indicates that the three basic conditions are likely to be fulfilled by anadromous salmon: (i) the spawning habitat is discrete and populations are spatially separated by unsuitable habitat; (ii) some asynchrony is present in the dynamics of more or less distant populations and (iii) dispersal links populations because some salmon stray from their natal population. The implications of some peculiarities of salmon life history traits, unusual in classical metapopulations, are also discussed. Deeper understanding of the population structure of anadromous fish will be advanced by future studies on specific topics: (i) criteria must be defined for the delineation of suitable habitats that are based on features of the biotope and not on the presence of fish; (ii) the collection of long‐term data and the development of improved methods to determine age structure are essential for correctly estimating levels of asynchrony between populations and (iii) several key aspects of dispersal are still poorly understood and need to be examined in detail: the spatial and temporal scales of dispersal movements, the origin and destination populations instead of simple straying rates, and the relative reproductive success of immigrants and residents.     

An investigation of the biological basis of return variability for sockeye salmon (Oncorhynchus nerka) from Great Central and Sproat lakes,Vancouver Island

We tested whether variations in stock characteristics (spawner and smolt abundance) and biotic conditions (prey variability, predation, competition) during the early marine period explained variations in the return of sockeye salmon (Oncorhynchus nerka) to Great Central and Sproat lakes, adjacent lakes on the west coast of Vancouver Island. There are two freshwater age groups in each lake; fish spend 1 or 2 yrs in freshwater after hatching. We tested the influences of stock and biotic factors on the return of each of the two age groups from each of the two lakes. Results of regression analyses showed that prey biomass variability best explained the variation in return for all lake‐age groups. Euphausiid (Thysanoessa spinifera) and cladoceran (Evadne) prey biomass variability explained between 0.75 and 0.95 (adjusted R²) of the variation in return. There appear to be instances of a mismatch between the seasonality of prey productivity and the apparent critical period of feeding for juvenile sockeye.     

Population reproductive potential: Evaluation of long-term stock productivity

In order to avoid recruitment overfishing, fish stocks must have sufficient reproductive ability. The spawning stock biomass (SSB), which ignores the value of immature fish, is widely used as an index of stock sustainability. From the perspective of sustainability, immediate reproduction, as well as future spawning, must be considered. We developed an index of long-term stock productivity, called the population reproductive potential (PRP). PRP is defined as the expected total reproductive value of the standing stock. We used PRP to assess the western Atlantic bluefin tuna (WBT) stock. The trends in SSB, numbers (N), biomass and PRP of WBT are inconsistent when compared to each other, due to fluctuation in age composition. We evaluated the long-term productivity of WBT by computer simulation and compared the result with trends in the abundance indices. The result of the computer simulation was highly consistent with the trend in the PRP. Short-term trends in SSB and N often do not reflect long-term stock trends, because they are highly sensitive to age-composition dynamics. The PRP is useful for evaluating stock trends, especially when the age composition is unstable.     

Is Guam a regional source,destination, or stepping‐stone for larvae of three fisheries species?

The pelagic larval duration (PLD) period of fish can influence dispersal, recruitment, and population connectivity, thereby potentially informing best strategies for fisheries management. Computer models were used to simulate the dispersal of larvae of three species, representing a range of PLDs, from the Pacific island of Guam and neighboring islands for a 9‐year period (2004–2012) to gain insight into the best management strategies for these species. The species included two springtime spawners with shorter and longer PLDs, scribbled rabbitfish (Siganus spinus; 33 days) and yellowfin goatfish (Mulloidichthys flavolineatus; ~90 days maximum), and a fall spawner with a similarly long PLD, bluespine unicornfish (Naso unicornis; ~94 days maximum). An ocean circulation model coupled with a particle dispersal model provided simulated numbers of larvae settling at each island in relation to the island where they were spawned. Graph analysis was used to examine generational connections between islands. For S. spinus, self‐seeding was the dominant means of replenishment at Guam. Local management actions to maintain adequate spawning stock may be a primary control on long‐term sustainability for that fishery. In contrast, N. unicornis and M. flavolineatus populations at Guam were reliant on outside sources for 92%–98% of larval supply. For them, identifying and negotiating the preservation of upstream spawning potential in the Marshall Islands and Federated States of Micronesia will be needed. Guam played a relatively minor role in generational connectivity across the region. Shortest paths spanning the region often did not pass through Guam, or there were equally short paths through other islands.     

The slow rise of technology: Computer vision techniques in fish population connectivity

1. Technological advancements in data collection and analysis are producing a new generation of ecological data. Among these, computer vision (CV) has received increased attention for its robust capabilities for rapidly processing large volumes of digital imagery.
2. In marine ecosystems, the study of fish connectivity provides fundamental information for assessing fisheries stocks, designing and implementing protected areas and understanding the impact of habitat loss. While the field of fish connectivity has benefited from technological advancements, the extent to which novel techniques, such as CV, have been utilized has not been assessed. To inform future directions and developments, this study reviewed the current use of CV in fish connectivity research, quantified how the implementation of such technology in fish connectivity research compared with other areas of marine research and described how this field could benefit from CV.
3. The review found that the use of remote camera systems in fish connectivity research is increasing, but the implementation of automated analysis of digital imagery has been slow. Successful implementation and expansion of CV frameworks in aquaculture and coral reef ecology suggest that CV techniques could greatly benefit fish connectivity research.
4. A case study of potential use of CV in fish connectivity research, scaling up optimal foraging models to predict marine population connectivity, highlights how beneficial it could be.
5. The capacity for CV techniques to be adopted alongside traditional approaches, the unparalleled speed, accuracy and reliability of these approaches and the benefits of being able to study ecosystems along multiple spatial–temporal scales, all make CV a valuable tool for assessing connectivity. Ultimately, these technologies can assist data-driven decisions that directly influence the health and productivity of marine ecosystems.

     

Variations in reproductive potential between nearshore and offshore spawning contingents of hogfish in the eastern Gulf of Mexico

Hogfish, Lachnolaimus maximus (Walbaum), have multiple traits that confound measuring reproductive potential: they are protogynous, relatively long‐lived harem‐forming fish that spawn daily for months. Additionally, recent evidence demonstrates that size, age and timing of sex change vary on a spatial scale within the study area (West Florida shelf, USA). This study investigates the effect of this spatial variation on hogfish reproductive potential by evaluating spawning seasonality, spawning frequency and batch fecundity using an indeterminate egg production model. Offshore females were larger than nearshore females, and batch fecundity was related in a log‐linear manner to female size. Gonad histology demonstrated a more protracted reproductive period for females offshore (8 months) than nearshore (4 months). Spatial variations in size coincide with ontogeny because hogfish move offshore with growth; however, even after accounting for fish size, offshore females spawned more. In areas where male removal rates are elevated, spawning harems are disrupted; thus, greater fishing effort nearshore may further reduce the reproductive potential of these females. These nearshore and offshore spawning components of the population are not genetically distinct, but instead represent two contingent spawning strategies that likely enhance total population stability and resilience of this stock in the eastern Gulf of Mexico.     

Spatial abundance patterns and recruitment of a virus‐affected commercial mollusc fishery

Infectious pathogens figure prominently among those factors threatening marine wildlife. Mass mortality events caused by pathogens can fundamentally alter the structure of wild fish stocks and depress recruitment rates and yield. In the most severe instances, this can precipitate stock collapses resulting in dramatic economic losses to once valuable commercial fisheries. An outbreak of a herpes‐like virus among commercially fished abalone populations in the south‐west fishery of Victoria, Australia, during 2006–2007, has been associated with high mortality rates among all cohorts. Long‐term records from fishery‐independent surveys of blacklip abalone Haliotis rubra (Leach) enabled abundance from pre‐ and post‐viral periods to be analysed to estimate stock density and biomass. The spatial distribution of abundance in relation to physical habitat variables derived from high‐resolution bathymetric LiDAR data was investigated. Significant differences were observed in both measures between pre‐ and post‐viral periods. Although there was some limited evidence of gradual stock improvement in recent years, disease‐affected reefs have remained below productivity rates prior to the disease outbreak suggesting a reduction in larval availability or settlement success. This was corroborated by trends in sublegal sized blacklip abalone abundance that has yet to show substantial recovery post‐disease. Abundance data were modelled as a function of habitat variables using a generalised additive model (GAM) and indicated that high abundance was associated with complex reef structures of coastal waters (<15 m). This study highlights the importance of long‐term surveys to understand abalone recovery following mass mortality and the links between stock abundance and seafloor variability.