Productivity and Susceptibility Analysis (PSA) in Fisheries: A Comprehensive Framework for Risk Assessment in Data-Limited Contexts

Abstract

This article provides a comprehensive examination of Productivity and Susceptibility Analysis (PSA), a pivotal semi-quantitative risk assessment tool for evaluating the vulnerability of fish stocks and sensitive species in data-poor fisheries. Tailored for researchers, scientists, and development professionals, it explores PSA's foundational principles, detailing how life-history productivity and fishery-specific susceptibility attributes are scored and combined. The scope covers methodological applications across diverse global fisheries, from Peruvian groundfish to Mediterranean bycatch, and addresses critical challenges in data quality, scoring thresholds, and the integration of alternative knowledge systems like fishers' expertise. Finally, it evaluates the validation of PSA outcomes against independent stock status information and compares its utility with other assessment frameworks, offering insights into its role in guiding sustainable management and conservation priorities.

Understanding PSA Foundations: Core Principles, Components, and Applications in Data-Poor Fisheries

Productivity and Susceptibility Analysis (PSA) is a semi-quantitative risk assessment methodology designed to evaluate the vulnerability of fish stocks to overfishing, particularly in contexts where data are insufficient for formal stock assessment models [1]. It operates by measuring two core dimensions: biological productivity, which reflects the stock's capacity for growth and recovery, and susceptibility, which indicates its exposure to and impact from fishing pressures [2]. By integrating scores from these dimensions, PSA produces a composite vulnerability score. This framework forms an integral component of risk-based fishery management strategies, such as the Marine Stewardship Council's Risk-Based Framework, and is a critical tool for prioritizing management actions in data-poor scenarios common in small-scale and developing world fisheries [1] [3].

Core Components of the PSA Framework

The PSA framework decomposes the overarching concept of vulnerability into two independently scored indices, each defined by specific, measurable attributes.

• Productivity (P) Index: This index quantifies the innate biological capacity of a species or stock to withstand and recover from fishing mortality. It is derived from life history traits. Key attributes include growth rate, natural mortality, age at maturity, maximum age or size, fecundity, and breeding strategy [2]. Species with "fast" life histories (e.g., high growth, early maturity) typically receive lower risk scores as they are intrinsically more productive and resilient.
• Susceptibility (S) Index: This index measures the potential for a stock to be negatively impacted by a specific fishery. It evaluates the overlap and interaction between the stock and the fishing operation. Attributes include geographic overlap, seasonal overlap, the selectivity and mortality of fishing gear, the existence and effectiveness of management measures, and the fishery's economic and social value [2].

Each attribute within the Productivity (typically 10 attributes) and Susceptibility (typically 12 attributes) indices is assigned a risk score, usually on a scale from 1 (low risk) to 3 (high risk) [1]. The overall scores for P and S are calculated, often as the Euclidean distance from the origin (√(P² + S²)) or another aggregation method, to generate a final Vulnerability (V) score [3].

Table 1: Standardized PSA Scoring Attributes and Metrics

Index	Representative Attributes	Score Range	Data Sources
Productivity (P)	Growth rate (k), Natural mortality (M), Age at maturity, Maximum size, Fecundity, Trophic level [2]	1 (Low Risk) to 3 (High Risk)	Scientific literature, FishBase, local studies, expert knowledge [2]
Susceptibility (S)	Geographic overlap, Seasonal overlap, Gear selectivity/post-capture mortality, Management effectiveness, Habitat impact, Economic value [2]	1 (Low Risk) to 3 (High Risk)	Fisher interviews, logbooks, observer data, management plans

Application in Data-Poor Fisheries Research: A 2024 Case Study

A seminal application of PSA is demonstrated in a 2024 study assessing ten coastal groundfish species in Peru, a classic data-poor scenario [3]. The research adapted the PSA methodology to integrate susceptibility from multiple gear types into a single score per attribute. Given that official data quality was rated between "limited" and "no data," the study relied significantly on structured expert knowledge to inform scoring [3].

The results successfully stratified species by vulnerability, providing clear priorities for management. The analysis found that four species, including the broomtail grouper and Pacific goliath grouper, showed extremely high vulnerability (V > 2.2), highlighting them as urgent conservation priorities [3]. This application validated the utility of PSA thresholds for informing management even in extremely data-limited contexts.

Table 2: PSA Vulnerability Assessment of Peruvian Coastal Groundfish (2024 Study) [3]

Species	Vulnerability Score (V)	Risk Category
Broomtail grouper	2.57	Extremely High
Grape-eye seabass	2.50	Extremely High
Pacific goliath grouper	2.28	Extremely High
Galapagos sheephead wrasse	2.24	Extremely High
Black snook	2.13	High
Harlequin wrasse	2.05	High
Chino	2.01	High
Mulata	1.89	Medium
Pacific beakfish	1.89	Medium
Bumphead parrotfish	1.73	Low

Detailed Experimental Protocol for PSA

The following protocol provides a standardized methodology for conducting a PSA, synthesizing guidelines from established toolkits and applications [2] [4].

Phase 1: Pre-Assessment Scoping

• Define Assessment Unit: Identify the target stock(s) and the specific fishery or gear types to be evaluated.
• Assemble Team: Form a multidisciplinary team comprising biologists, fishery scientists, fishery managers, and where possible, experienced fishers.
• Gather Baseline Data: Collate all available information, including scientific papers, government reports, fishery statistics, and database entries (e.g., FishBase).

Phase 2: Data Collection and Attribute Scoring

• Productivity Scoring:
  • For each productivity attribute (e.g., growth rate, age at maturity), extract the best available quantitative or qualitative value.
  • Assign a score of 1, 2, or 3 based on pre-defined, attribute-specific thresholds. For example, a high natural mortality rate (M > 0.8) may score a 1 (low risk), while a very low rate (M < 0.2) may score a 3 (high risk).
  • Document the data source and assign a data quality score (e.g., 1=good, 2=limited, 3=poor/no data) for each attribute.
• Susceptibility Scoring:
  • For each susceptibility attribute, analyze fishery operations. For gear selectivity, score the likelihood of capture and mortality if captured.
  • For geographic/seasonal overlap, estimate the proportion of the stock's distribution and annual cycle that intersects with the fishery.
  • Score each attribute from 1 (low susceptibility) to 3 (high susceptibility) and document data quality.

Phase 3: Calculation and Analysis

• Calculate Indices: Compute the overall Productivity (P) and Susceptibility (S) scores. A common method is the Euclidean distance: ( P = \sqrt{\sum (pi^2)/n} ) and ( S = \sqrt{\sum (si^2)/m} ), where (pi) and (si) are individual attribute scores, and (n) and (m) are the number of attributes.
• Calculate Vulnerability: Compute the overall vulnerability score, typically ( V = \sqrt{P^2 + S^2} ).
• Plot and Interpret: Plot stocks on a two-dimensional PSA plot (P on the y-axis, S on the x-axis). Quadrants or diagonal risk lines can be drawn to categorize stocks as low, medium, high, or very high risk [2].
• Sensitivity Analysis: Test the sensitivity of scores to changes in uncertain attributes or to different weighting schemes.

Phase 4: Reporting and Management Integration

• Document Uncertainty: Clearly report data quality scores and the role of expert judgment.
• Prioritize Actions: Use the risk categorization to recommend management actions. High-vulnerability stocks may need immediate catch limits or spatial protections, while low-vulnerability stocks may be suitable for increased harvest [3].
• Identify Knowledge Gaps: Use the data quality assessment to prioritize future research and monitoring needs.

The Scientist's Toolkit: Essential Research Reagents and Materials

Table 3: Key Research Reagents and Tools for PSA Implementation

Item	Function in PSA	Example/Source
Standardized PSA Worksheet	Provides the structured framework for data entry, scoring, and calculation, ensuring consistency [4].	EDF Fishery Solutions Center tool (available in .xls format) [4].
Life History Parameter Database	Primary source for scoring productivity attributes when local data are absent [2].	FishBase (www.fishbase.org), SeaLifeBase.
Expert Elicitation Protocol	Structured method (e.g., Delphi technique) to formally gather and quantify expert knowledge for scoring in data-poor contexts [3].	Custom survey instruments, facilitated workshops.
Geographic Information System (GIS)	Analyzes spatial data to quantify geographic overlap between stock distributions and fishing effort, a key susceptibility attribute.	ArcGIS, QGIS, R packages (sf, raster).
Statistical Software	Performs sensitivity analyses, visualizes results (PSA plots), and manages data.	R, Python (Pandas, Matplotlib), Stata.
Data Quality Scoring Rubric	A predefined scale to consistently document the reliability of information used for each attribute score [3].	Typically a 3-point scale (Good, Limited, Poor/No Data).

Visualizing the PSA Framework and Workflow

PSA Methodological Workflow

PSA Scoring Logic and Aggregation

Productivity-Susceptibility Analysis (PSA) is a semi-quantitative risk assessment framework developed to evaluate the vulnerability of fish stocks to overfishing, particularly in data-limited circumstances [1]. It functions as a critical component within hierarchical risk-based frameworks, such as the Marine Stewardship Council's (MSC) Risk-Based Framework, to prioritize species for more resource-intensive, quantitative assessment [5]. The foundational premise of PSA is that a stock's overall vulnerability (V) is a function of two core properties: its Productivity (P), representing the intrinsic capacity for population growth and recovery based on life-history traits, and its Susceptibility (S), representing the intensity and nature of its interactions with a fishery [5] [1].

The method involves scoring a standardized set of attributes for productivity and susceptibility, typically on a scale from 1 (low risk) to 3 (high risk) [1]. The overall productivity score (P) is calculated as the arithmetic mean of its attribute scores, while the overall susceptibility score (S) is often calculated as the geometric mean, reflecting the multiplicative nature of risk factors [5]. The final vulnerability is derived from the Euclidean distance from the origin: V = √(P² + S²) [5]. This score places a stock within a risk matrix, categorizing it as having low, medium, or high vulnerability to inform management priorities [3].

This document provides detailed application notes and protocols for deconstructing and applying the core productivity and susceptibility attributes within a broader PSA research context, integrating contemporary critiques and methodological advances.

Deconstructing Productivity Attributes (Life-History)

Productivity attributes are proxies for a population's intrinsic rate of increase and resilience to depletion. They are derived from life-history parameters, which are themselves responsive to both fishing pressure and environmental drivers like temperature [6].

Core Attribute Definitions & Scoring Protocols

The following table details common productivity attributes, their biological basis, typical scoring criteria, and protocols for data collection or estimation.

Table 1: Core Productivity Attributes, Scoring Criteria, and Measurement Protocols

Attribute	Biological Basis & Role in Productivity	Typical Scoring Criteria (Risk: 1=Low, 3=High)	Measurement & Estimation Protocols
Average Age at Maturity (Amat)	Indicates generation time and speed of reproductive turnover. Earlier maturation generally supports higher productivity.	1 (Low): < 5 yrs [5].2 (Medium): 5-15 yrs [5].3 (High): > 15 yrs [5].	Protocol: Estimate from age-length keys using histological examination of gonads. Use representative sampling across seasons. Data-poor approach: Use published values from congeneric species or regional databases, applying a precautionary bias toward higher risk if uncertain.
Average Maximum Age (Amax)	Proxy for natural mortality rate (M). Shorter-lived species with higher M typically have higher potential population growth rates.	1 (Low): < 10 yrs.2 (Medium): 10-25 yrs.3 (High): > 25 yrs.	Protocol: Estimate from the oldest observed individual in age-validated samples (e.g., via otoliths). Caution: Heavily fished populations may exhibit truncated age distributions, confounding this metric.
Average Maximum Size (Lmax or L∞)	Correlated with fecundity, mortality, and growth rate. Smaller asymptotic size often associates with faster life histories.	1 (Low): < 30 cm.2 (Medium): 30-100 cm.3 (High): > 100 cm.	Protocol: Estimate von Bertalanffy growth parameter L∞ from length-frequency analysis or mark-recapture data.
Von Bertalanffy Growth Coefficient (K)	Describes the rate at which asymptotic size is approached. Faster growth (higher K) supports higher productivity.	1 (Low): > 0.3 yr⁻¹.2 (Medium): 0.1-0.3 yr⁻¹.3 (High): < 0.1 yr⁻¹.	Protocol: Estimate concurrently with L∞ from ageing data. Note: K and L∞ are strongly inversely correlated.
Fecundity (F)	Average number of eggs produced per female per spawning event. Higher fecundity generally increases reproductive potential.	1 (Low): > 20,000 [5].2 (Medium): 100-20,000 [5].3 (High): < 100 [5].	Protocol: Conduct gravimetric fecundity counts from mature ovaries prespawning. Adjust for batch spawning where applicable.
Spawning Strategy	Frequency and periodicity of reproduction affect the population's ability to compensate for losses.	1 (Low): Multiple batch spawner, prolonged season.3 (High): Total spawner, brief seasonal window.	Protocol: Determine via gonadosomatic index (GSI) tracking over an annual cycle and histological analysis.
Natural Mortality Rate (M)	The instantaneous rate of death from non-fishing causes. A primary determinant of population turnover rate.	1 (Low): M > 0.8 yr⁻¹.2 (Medium): M = 0.2-0.8 yr⁻¹.3 (High): M < 0.2 yr⁻¹.	Protocol: Estimate using empirical relationships (e.g., Pauly's formula linking M to K, L∞, and temperature), or from catch curve analysis in pre-exploitation data.

Advanced Application Note: Climate-Driven Trait Dynamics

Life-history traits are not static. A 2020 meta-analysis of 332 Indo-Pacific fish species demonstrated that increasing temperature systematically shifts traits toward "faster" life histories: growth coefficient (K) and natural mortality (M) increase, while asymptotic length (L_∞) and age at maturity (A_mat) decrease [7]. However, the net effect on population growth potential is not uniform and depends on a species' position on the fast-slow life-history continuum [7].

Table 2: Effect of +1°C Warming on Life-History Traits and Implications for PSA Scoring [7]

Trait	Average Direction of Change	Magnitude of Change (Mean)	Implication for PSA Productivity Scoring
Growth Coefficient (K)	Increase	+0.05 yr⁻¹ °C⁻¹	Tendency toward lower risk score (higher productivity).
Natural Mortality (M)	Increase	+0.05 yr⁻¹ °C⁻¹	Tendency toward lower risk score (higher productivity).
Asymptotic Length (L∞)	Decrease	-0.02 cm °C⁻¹	Tendency toward lower risk score (higher productivity).
Age at Maturity (A50)	Decrease	-0.04 yr °C⁻¹	Tendency toward lower risk score (higher productivity).
Net Population Growth (r)	Divergent	Varies by group	Fast-life-history species: Population growth rate tends to decrease. Slow-life-history species: Population growth rate tends to increase.

Protocol Note: For climate-resilient PSA, researchers should source life-history parameters from temperature-appropriate studies or apply known thermal coefficients [7]. For long-lived, slow-growth species (inherently high-risk), warming may marginally improve productivity scores, but this may be offset by other climate impacts (e.g., oxygen stress, habitat loss). The trait changes themselves can be used as a sensitivity analysis within the PSA framework.

Deconstructing Susceptibility Attributes (Fishery Interaction)

Susceptibility attributes measure the nature and extent of a population's interaction with fishing operations, determining the fraction of the population exposed to fishing mortality.

Core Attribute Definitions & Scoring Protocols

Susceptibility is often decomposed into components of availability, encounterability, and post-capture mortality [5]. The following table outlines key attributes and assessment methodologies.

Table 3: Core Susceptibility Attributes, Scoring Criteria, and Assessment Protocols

Attribute	Definition & Role in Susceptibility	Typical Scoring Criteria (Risk: 1=Low, 3=High)	Assessment Protocol
Spatial Overlap	Degree of geographic co-occurrence between the stock's distribution and the fishery's operational area.	1 (Low): Minimal overlap (<10% of range).3 (High): Complete overlap.	Protocol: Analyze fishery logbook data (set locations) against species distribution models from survey data or expert maps. Use GIS to calculate percentage area overlap.
Seasonal Overlap	Temporal coincidence between fishery activity and the stock's presence in the fishing grounds.	1 (Low): Fishery active outside peak abundance season.3 (High): Year-round operation during peak abundance.	Protocol: Compare monthly fishery effort curves with seasonal abundance indices from surveys or fishery-independent data.
Habitat Vulnerability	Susceptibility based on the habitat occupied (e.g., depth, substrate) relative to gear selectivity.	1 (Low): Inhabits untargeted habitat (e.g., deep reef, rough bottom).3 (High): Inhabits primary targeted habitat (e.g., flat sand).	Protocol: Based on known species habitat preferences and gear specifications (e.g., trawl footrope contact). Expert judgment is often required.
Gear Selectivity & Encounterability	Probability of capture given co-occurrence, based on size, behavior, and gear characteristics.	1 (Low): Size/behavior avoids gear (e.g., too small, escapes meshes).3 (High): Gear is highly size-selective for the species.	Protocol: Analyze size-frequency data in catch vs. population. Review gear studies (e.g., selectivity curves). For bycatch, use observer data on escapee numbers.
Post-Capture Survival (Discard Mortality)	Probability of survival if the individual is captured and discarded.	1 (Low): High survival (>75%) after release.3 (High): Very low survival (<25%).	Protocol: Conduct tag-and-release vitality studies, or hold animals in tanks post-capture. Condition factors (barotrauma, injury score) can be used as proxies.
Management Effectiveness	Capacity of existing regulations to control exploitation on the stock.	1 (Low): Effective, enforced catch/effort limits, spatial closures.3 (High): No relevant management measures.	Protocol: Qualitative assessment based on review of management plans, existence of monitoring control and surveillance (MCS), and compliance studies.

Application Note: Integrating Multiple Gears in Data-Poor Contexts

In complex small-scale fisheries, a single species may be targeted or captured by multiple gear types. A 2024 Peruvian groundfish PSA developed a protocol to integrate multi-gear susceptibility [3]:

• Score by Gear: Assess susceptibility attributes separately for each major gear type (e.g., gillnet, handline, trawl) interacting with the stock.
• Weight by Relative Impact: Assign a weight (wᵢ) to each gear type based on its relative estimated catch or effort for the species, where Σwᵢ = 1.
• Calculate Composite Score: For each susceptibility attribute, calculate a composite score: S_comp = Σ (wᵢ × Scoreᵢ). Use these composite attribute scores to calculate the overall geometric mean susceptibility [3].

Integrated PSA Workflow & Computational Protocol

The following diagram and protocol detail the step-by-step process for conducting a PSA, from data assembly to risk categorization.

Comprehensive PSA Protocol:

Phase 1: Data Assembly

• Compile species-specific data for all productivity and susceptibility attributes.
• Prioritize peer-reviewed literature, regional databases (e.g., FishBase), stock assessment reports, and fishery-independent survey data.
• For missing parameters, initiate structured expert elicitation: convene a panel of ≥3 independent experts, provide background data, use a modified Delphi process to reach consensus on point estimates and data quality scores [3].

Phase 2: Attribute Scoring

• For each attribute, assign a risk score (1, 2, or 3) based on predefined categorical thresholds (see Tables 1 & 3).
• Simultaneously, assign a Data Quality (DQ) score for each attribute (e.g., 1=High quality, empirical data; 2=Medium, inferred data; 3=Low, guess/analogy) [3].
• Document the rationale and source for every score.

Phase 3: Score Aggregation

• Calculate the overall Productivity score (P) as the arithmetic mean of the 7 productivity attribute scores.
• Calculate the overall Susceptibility score (S) as the geometric mean of the 4 susceptibility attribute scores: S = (s₁ * s₂ * s₃ * s₄)^{¹/₄}.
• Calculate the Vulnerability score (V) using the Euclidean distance formula: V = √(P² + S²) [5].

Phase 4: Risk Categorization & Uncertainty Analysis

• Categorize risk using standard thresholds: Low Risk (V < 2.64), Medium Risk (2.64 ≤ V ≤ 3.18), High Risk (V > 3.18) [5].
• Conduct sensitivity analysis: Recalculate V by varying scores for attributes with low DQ scores or by testing alternative weighting schemes.
• Conduct Monte Carlo simulation: Treat attribute scores as probability distributions (e.g., uniform across possible categories) to propagate uncertainty and generate probability distributions for P, S, and V. This quantitatively expresses classification certainty.

Critical Evaluation & Integration within a Research Thesis

The PSA is a prioritization tool, not a definitive stock assessment. Its assumptions and limitations must be explicitly addressed within a thesis context [5].

Key Strengths:

• Rapid Application: Allows assessment of many species with limited data [1].
• Transparent Structure: Makes rationale for risk categorization explicit via scores.
• Effective Prioritization: Successfully identifies species with high-risk life histories and high exposure, directing resources to where they are most needed [3].

Key Limitations & Critical Assumptions:

• Static Life History: Assumes productivity traits are fixed, ignoring density-dependent and climate-driven responses [7] [6].
• Additivity & Independence: Assumes attributes are independent and combine linearly/multiplicatively, which may not reflect biological reality [5].
• No Population Dynamics: Does not model actual stock size, recruitment, or depletion level. A "low vulnerability" species could be overfished, and a "high vulnerability" species could be abundant [5].
• Threshold Sensitivity: Arbitrary category thresholds can drastically alter scores [5].

Protocol for Contextualizing PSA within a Thesis:

• Comparative Analysis: Where possible, compare PSA outcomes for a stock with trends in independent abundance indices (e.g., CPUE, survey biomass) to test its predictive power [3].
• Link to Operating Models: As noted by [5], the data required for a PSA often approximates that needed to populate a simple population dynamics model. A robust thesis should treat PSA as a first step, followed by developing a simple stochastic population projection model for high-priority species using the same life-history parameters. This allows testing harvest strategies and quantifying risks in terms of biomass and probability of overfishing.
• Climate Integration: Incorporate climate scenarios by adjusting life-history input parameters according to thermal relationships (Table 2) and re-running the PSA and/or population projections [7].

The Researcher's Toolkit: Essential Materials & Reagents

Table 4: Key Research Reagent Solutions for PSA and Life-History Studies

Category	Item / Solution	Function & Application in Protocol
Field Sampling & Morphology	1. Otolith Sectioning Kit:Low-speed saw, mounting epoxy, grinding/polishing papers, microscope slides.	To prepare thin sections of otoliths (ear stones) for precise age determination, validating Amax and growth parameters.
	2. Gonadal Histology Suite:Neutral buffered formalin (10%), ethanol series (70%, 95%, 100%), paraffin embedding system, microtome, H&E stain.	To fix, section, and stain gonadal tissue for microscopic determination of sex, maturity stage, and spawning frequency.
	3. Digital Calipers & Balance	For accurate measurement of length (mm) and weight (g) for growth and fecundity analyses.
Laboratory Analysis	4. Image Analysis Software (e.g., ImageJ)	To analyze digital images of otolith sections (for annuli counts) and ovary sections (for oocyte counts and sizing).
	5. Fecundity Estimation Kit:Gilson's fluid (or formalin), Petri dishes, subsampling grid, fine forceps.	To chemically break down ovarian connective tissue and facilitate the counting of oocytes for gravimetric fecundity estimation.
Data Analysis & Modeling	6. Statistical Software (R, Python with pandas)	For data cleaning, analysis, growth curve fitting (e.g., using 'FSA' or 'fishmethods' packages in R), and statistical testing.
	7. Population Modeling Platform (e.g., R popbio, MSEtool)	To construct age/stage-structured matrix models or operating models for population projections and LTRE analysis beyond PSA [6].
	8. GIS Software (e.g., QGIS, ArcGIS)	To analyze spatial overlap between species distributions (from survey data or SDMs) and fishing effort layers from VMS/logbooks.
Expert Elicitation	9. Structured Elicitation Protocol Template	A standardized framework (e.g., based on the IDEA protocol) to systematically gather, weight, and aggregate expert judgments for missing data parameters [3].

Deconstructing productivity and susceptibility attributes is fundamental to applying the PSA framework. This document provides detailed protocols for attribute scoring, data estimation, and workflow execution, informed by current methodological applications and critiques. For thesis research, the greatest analytical value lies in moving beyond static PSA scoring. Researchers should employ PSA as a structured prioritization tool, then deepen their investigation by integrating climate-driven trait dynamics, propagating uncertainty, and linking results to dynamic population models. This approach addresses the core limitations of PSA while fulfilling its original purpose: focusing sophisticated analytical resources on the species and systems at greatest risk.

Productivity and Susceptibility Analysis (PSA) serves as a critical risk assessment framework in fisheries science, particularly where traditional, data-intensive stock assessments are unfeasible. By evaluating a stock's biological productivity and its susceptibility to fishing pressure, PSA provides a semi-quantitative measure of vulnerability to overfishing [1]. This methodology is indispensable for prioritizing management interventions in data-poor regions, which constitute a significant portion of the world's fisheries [3]. Framed within a broader thesis on advancing fisheries productivity research, this article details the application, protocols, and essential toolkit for implementing PSA, positioning it as a cornerstone for evidence-based management and sustainable aquatic food systems under FAO's Blue Transformation agenda [8].

The Global Context: Data Gaps and Management Imperatives

Effective fisheries management is the most powerful tool for conserving resources, yet its implementation is uneven globally [8]. The latest FAO assessment of 2,570 marine fish stocks reveals that while 64.5% are fished within biologically sustainable levels, 35.5% are classified as overfished [8]. This disparity underscores a stark management gap.

Table 1: Regional Variation in Fishery Sustainability (FAO 2025 Assessment) [8]

FAO Fishing Area	Region	% Stocks Fished Sustainably	Key Management Context
Areas 48, 58, 88	Antarctic	100%	Ecosystem-based management & strong international cooperation (CCAMLR).
Area 67	Northeast Pacific	92.7%	Long-term investment, robust & precautionary management frameworks.
Area 81	Southwest Pacific	85%	Effective regional management institutions.
Area 37	Mediterranean & Black Sea	35.1%	Early signs of recovery; fishing pressure dropped 30% since 2013.
Area 87	Southeast Pacific	46%	Limited institutional capacity, fragmented governance, major data gaps.
Area 34	Eastern Central Atlantic	47.4%	Central to food security; constrained by data and capacity gaps.

Sustainability rates far exceed the global average in regions with strong, science-based governance [8]. Conversely, areas with limited capacity and major data gaps—often in low and middle-income countries where fisheries are crucial for food security and livelihoods—face significant challenges [8]. These data-poor contexts, prevalent in small-scale and coastal fisheries, are where PSA is most vitally applied to inform precautionary management [3].

PSA Methodology: Framework and Calculation

PSA is a semi-quantitative risk assessment method that evaluates vulnerability based on two dimensions: Productivity (P), the biological capacity of a stock to withstand and recover from fishing pressure, and Susceptibility (S), the likelihood and intensity of its encounter with fisheries [1] [2].

Core Attributes and Scoring

Scores from 1 (low risk) to 3 (high risk) are assigned to a standardized set of attributes.

Table 2: Standard PSA Attributes and Scoring Criteria

Dimension	Typical Attributes	Score 1 (Low Risk/Vulnerability)	Score 2 (Medium Risk/Vulnerability)	Score 3 (High Risk/Vulnerability)
Productivity	Growth rate, Age at maturity, Fecundity, Natural mortality, Trophic level [2].	High growth rate, Early maturity, High fecundity.	Moderate life-history traits.	Slow growth, Late maturity, Low fecundity.
Susceptibility	Geographic overlap, Seasonal overlap, Selectivity/gear catchability, Management effectiveness, Post-capture mortality [2].	No/minimal overlap with fishery, Gear avoids species, Strong management.	Partial or seasonal overlap.	Complete overlap, Non-selective gear, No effective management.

Integration and Vulnerability Classification

The overall Productivity (P) and Susceptibility (S) scores are calculated, often as the arithmetic mean of their respective attributes. These two scores are then combined into a Vulnerability (V) score, typically using the Euclidean distance formula: V = √(P² + S²) [3]. The resulting score categorizes risk:

• Low Vulnerability: V < 1.8
• Moderate Vulnerability: 1.8 ≤ V ≤ 2.0
• High Vulnerability: V > 2.0 [2]

Application Notes: PSA in Data-Poor Scenarios

A 2024 study on Peruvian coastal groundfish exemplifies PSA's application in extremely data-poor contexts [3]. Facing a near-total absence of formal stock data, researchers adapted PSA by integrating expert knowledge from scientists and fishers to score attributes. The study assessed ten species, revealing a spectrum of vulnerability and providing the first scientific basis for prioritizing management.

Table 3: PSA Results for Selected Peruvian Coastal Groundfish [3]

Species (Common Name)	Calculated Vulnerability (V)	Risk Category	Management Implication
Broomtail grouper	2.57	High	Highest priority for immediate regulatory action (e.g., catch limits).
Pacific goliath grouper	2.28	High	High priority; large, slow-growing species highly vulnerable to depletion.
Black snook	2.13	High	Priority for data collection and precautionary catch restrictions.
Mulata	1.89	Moderate	Monitor catch trends; consider moderate management measures.
Bumphead parrotfish	1.73	Low	Lower management priority relative to others; sustainable harvest possible.

This case confirms that PSA effectively identifies vulnerable stocks even with "limited" to "no data" quality scores, and that established risk thresholds remain useful for guiding management priorities [3].

Experimental Protocols for PSA Implementation

Objective: To gather all necessary qualitative and quantitative data for scoring PSA attributes in the absence of robust fisheries-independent data.

Materials: PSA worksheet (paper or digital), literature databases (FishBase, SeaLifeBase), recording device for interviews, geographic information system (GIS) for overlap analysis.

Procedure:

• Desktop Review: Compile available biological parameters (e.g., maximum length, growth rate) for the target species from global databases and local scientific literature [2].
• Stakeholder Identification: Identify key informants: fishers, fishery officers, local biologists, and NGO staff with long-term area knowledge.
• Structured Interview: Conduct one-on-one or small-group interviews using the PSA attributes as a structured guide. For each attribute (e.g., "seasonal overlap"):
  • Present the scoring criteria (1-3).
  • Ask open-ended questions about the fishery and species behavior.
  • Collaboratively assign a score based on the informant's evidence.
• Data Triangulation: Compare and contrast scores from multiple experts. Resolve discrepancies through facilitated discussion to reach a consensus score for each attribute.
• Data Quality Recording: For each scored attribute, document the data quality (e.g., "Empirical Study," "Local Expert," "Inferred from Similar Species") [3].

Protocol 4.2: Scoring, Calculation, and Sensitivity Analysis

Objective: To calculate final Productivity, Susceptibility, and Vulnerability scores and test the robustness of conclusions.

Materials: Completed PSA data sheets, statistical software (e.g., R, Python) or spreadsheet software (e.g., Excel, Google Sheets).

Procedure:

• Input Data: Enter consensus scores for all Productivity (P) and Susceptibility (S) attributes into a spreadsheet.
• Calculate Dimension Scores: Compute the mean score for all P attributes (P̄) and all S attributes (S̄).
• Calculate Vulnerability: Apply the Euclidean distance formula: V = √(P̄² + S̄²) [3].
• Categorize Risk: Classify the stock as Low (V<1.8), Moderate (1.8≤V≤2.0), or High (V>2.0) vulnerability [2].
• Sensitivity Analysis:
  • Vary the scores for the most uncertain attributes within a plausible range (e.g., ±0.5).
  • Recalculate the Vulnerability score for each variation.
  • Determine if the risk category changes under different plausible scenarios. This defines the confidence in the primary result.
• Visualization: Plot the stock on a PSA plot (see Diagram 1) and generate a results table (see Table 3).

Visualization of the PSA Framework and Workflow

Diagram 1: PSA Implementation Workflow. This diagram outlines the process from data collection to management output.

Diagram 2: PSA Risk Plot and Interpretation. Stocks are plotted based on their P and S scores. Contour lines (conceptual) represent equal vulnerability levels, with thresholds at V=1.8 and V=2.0 separating risk categories [2].

The Scientist's Toolkit: Research Reagent Solutions

Table 4: Essential Toolkit for Conducting a PSA Study

Tool/Resource	Function/Description	Application in PSA
Standardized PSA Worksheets	Blank scoring sheets in multiple languages (e.g., English, Spanish) with predefined attributes [2].	Provides a consistent framework for data collection and scoring across all analysts and stocks.
Global Life History Databases	Online repositories (e.g., FishBase, SeaLifeBase) of species-specific biological parameters.	Primary source for scoring Productivity attributes (growth, maturity, fecundity) when local studies are absent [2].
Expert Elicitation Protocols	Structured interview guides and facilitation techniques for engaging fishers and local scientists.	Key methodology for scoring Susceptibility attributes and refining Productivity scores in data-poor contexts [3].
Geographic Information System (GIS)	Software for visualizing and analyzing spatial data (e.g., QGIS, ArcGIS).	Used to objectively assess geographic and seasonal overlap between fish stocks and fishing effort.
Sensitivity Analysis Scripts	Code (e.g., in R or Python) to automate the re-calculation of vulnerability scores under different scoring assumptions.	Tests the robustness of PSA conclusions and quantifies uncertainty, a critical step for credible results [3].
Data Quality Scoring Rubric	A defined scale (e.g., Empirical, Expert, Inferred, No Data) for rating the confidence in each attribute score.	Ensures transparency and allows managers to weight the confidence of the overall assessment [3].

Productivity and Susceptibility Analysis (PSA) is a semi-quantitative, risk-based framework used to assess the vulnerability of fish stocks to overfishing in data-limited contexts [1]. It functions as a critical component within broader fisheries productivity and susceptibility research, enabling the prioritization of management actions when traditional, data-intensive stock assessments are not feasible [3]. The method evaluates two primary dimensions: Productivity, representing the biological capacity of a stock to withstand and recover from fishing pressure (e.g., growth rate, age at maturity), and Susceptibility, representing the likelihood and intensity of its encounter with fishing operations (e.g., geographic overlap, gear selectivity) [2]. By integrating these scores, PSA generates a composite Vulnerability score, categorizing stocks from low to high risk [1]. This approach is particularly vital for assessing unregulated small-scale fisheries, bycatch species, and emerging fisheries where scientific data is sparse, often necessitating the integration of expert ecological knowledge and local stakeholder observations to fill critical information gaps [3] [9].

Application Notes: Key Use Cases and Findings

PSA is deployed across diverse fishery contexts to inform precautionary management. Its flexible, modular structure allows for adaptation to specific assessment needs, from single-species targets to complex multi-species assemblages.

• Assessing Data-Poor Target Stocks: In Peruvian coastal small-scale fisheries, a 2024 PSA of ten groundfish species revealed stark vulnerabilities. Expert knowledge was essential to complete assessments where data quality was rated between "limited" and "no data" [3]. Results prioritized four species for immediate management, including the critically vulnerable Pacific goliath grouper.
• Evaluating Bycatch and Ecosystem Vulnerability: PSAs are central to ecosystem-based management, evaluating the risk to non-target species. The methodology underpins risk assessments for protected, endangered, and threatened species (e.g., cetaceans, sea turtles) caught as bycatch. Surveys like the GoMMAPPS program collect distribution and abundance data that feed into susceptibility parameters for these species [10].
• Informing Niche and Baitfish Fisheries: For fisheries with highly specialized markets, such as the Marine Aquarium Trade (MAT), PSA is used to assess sustainability. A 2025 study applied PSA to the top 250 MAT species, finding significant data limitations for many, requiring generalizations at the genus or family level. The results provide stakeholders with a comparative sustainability ranking to guide harvest decisions [11].
• Multi-Scale and Participatory Assessments: PSA outcomes can vary with spatial scale and knowledge source. A 2022 study in Oregon integrated local fisher knowledge with scientific data for nearshore species. While results generally aligned with coastwide federal PSAs, local scores often indicated lower vulnerability, highlighting locally-specific stock conditions or fishing practices that broader assessments may miss [9].

Table 1: Comparative Vulnerability Scores from Recent PSA Case Studies

Case Study Context	Species / Stock	Productivity Score (1-3)	Susceptibility Score (1-3)	Vulnerability (V) Score	Risk Category
Peruvian Coastal Groundfish [3]	Broomtail Grouper	Not Specified	Not Specified	2.57	Extremely High
	Pacific Goliath Grouper	Not Specified	Not Specified	2.28	Extremely High
	Bumphead Parrotfish	Not Specified	Not Specified	1.73	Low
Oregon Nearshore Fishes [9]	Various (10 species)	Derived from life-history	Based on local expert knowledge	Generally lower than coastwide PSA	Low to Moderate
Marine Aquarium Trade (Top 250 spp) [11]	Species constituting 92.5% of trade	Modeled from life-history traits	Based on harvest factors & trade volume	Clustered via Gaussian mixture model	Spectrum from Low to High

Table 2: Standard PSA Scoring Thresholds and Risk Interpretation [1] [2]

Metric	Score Range	Qualitative Description	Management Implication
Attribute Scores	1	Low Risk / High Productivity / Low Susceptibility	Favorable trait requiring minimal management intervention.
	2	Medium Risk / Moderate Productivity or Susceptibility	Precautionary monitoring and management advised.
	3	High Risk / Low Productivity / High Susceptibility	Trait confers high vulnerability; warrants immediate management focus.
Final Vulnerability (V)	V < 1.8	Low Vulnerability Stock	Likely sustainable under current fishing pressure.
	1.8 ≤ V ≤ 2.0	Moderate Vulnerability Stock	Requires monitoring and potential precautionary catch limits.
	V > 2.0	High Vulnerability Stock	High risk of overfishing; needs stringent regulatory measures.

Experimental Protocols for PSA Execution

Conducting a robust PSA requires a structured, iterative protocol that rigorously documents data sources, scoring justifications, and uncertainty.

Protocol 1: Core PSA Workflow for a Target Stock

• Problem Formulation & Species Selection: Define the assessment's spatial scale (e.g., coastwide, local management area) [9] and select the stock or species complex for evaluation.
• Data Assembly & Gap Analysis: For each of the standard PSA attributes, compile available data from peer-reviewed literature, fisheries databases, and unpublished reports.
  • Productivity Attributes (e.g.,): Maximum age, size at maturity, fecundity, natural mortality rate, trophic level [2].
  • Susceptibility Attributes (e.g.,): Seasonal/geographic overlap with fishery, availability to gear, post-capture survival, management effectiveness [3] [2].
  • Data Quality Scoring: Rate the quality of information for each attribute (e.g., "Good," "Limited," "No Data") [3].
• Expert Elicitation Workshop: Convene a panel of experts (scientists, fishery managers, experienced fishers) to review compiled data, fill knowledge gaps using expert judgment, and independently score each attribute (1-3) [3] [9]. Use a modified Delphi approach to reach consensus on contentious scores.
• Vulnerability Calculation & Visualization: Calculate geometric mean scores for Productivity (P) and Susceptibility (S). Compute the composite Vulnerability score: V = √(P² + S²) [1]. Plot results on a bivariate PSA plot (P vs. S) and categorize risk using predefined thresholds (e.g., V > 2.0 = High) [2].
• Sensitivity & Uncertainty Analysis: Perform sensitivity tests by varying key scores within plausible ranges, especially for attributes with "Limited" data quality. Document how uncertainty propagates to the final risk categorization [11].
• Reporting & Management Integration: Document all scores, justifications, data sources, and expert panel composition. Translate findings into prioritized management recommendations (e.g., catch limits, spatial closures, gear restrictions for high-vulnerability stocks) [3].

Protocol 2: Integrated Visual Analysis for Bycatch Susceptibility

• Movement & Effort Data Fusion: Integrate geospatial datasets, including:
  • Species Movement Data: Satellite telemetry or acoustic tracking data for bycatch species (e.g., cetaceans, turtles) [12].
  • Fishing Effort Data: Vessel Monitoring System (VMS) or Automatic Identification System (AIS) data to map fishing activity in high resolution.
• Sequential Pattern & Rule Mining: Apply sequential pattern mining algorithms (e.g., VMSP) and rule-growth algorithms to the fused spatiotemporal dataset to identify frequent interaction pathways between fishing assets and marine animals [12].
• Interactive Visual Analytics: Use a multi-view visualization platform to allow analysts to interactively explore the mined patterns. This should include:
  • A geographic view displaying animal movement tracks and fishing effort heatmaps.
  • A timeline view showing temporal co-occurrence.
  • A ruleset browser displaying extracted behavioral rules (e.g., "if fishing set occurs in area X during period Y, species Z is observed with 75% confidence") [12].
• Susceptibility Parameterization: Use the identified patterns and rules to quantitatively inform PSA susceptibility attributes, such as "spatial overlap," "seasonal overlap," and "encounterability," moving beyond static, expert-based scores to dynamic, data-driven estimates [10] [12].

Figure 1: Core PSA methodology workflow from problem formulation to management integration, highlighting the iterative consensus-building process.

Figure 2: Multi-scale PSA framework showing integration of different knowledge sources and assessment scales for a comprehensive risk view [9].

The Scientist's Toolkit: Essential Research Reagent Solutions

Table 3: Key Reagents, Materials, and Tools for PSA Implementation

Tool/Reagent Category	Specific Item / Solution	Primary Function in PSA Protocol
Data Curation & Analysis	Standardized PSA Attribute Worksheets [2]	Provides structured template for consistent data compilation and scoring across all assessed stocks.
	Life-History Parameter Databases (e.g., FishBase, SeaLifeBase)	Primary source for scoring productivity attributes when local studies are absent [2].
	Geographic Information System (GIS) Software (e.g., QGIS, ArcGIS)	Essential for analyzing and visualizing spatial overlap between stock distributions and fishing effort to score susceptibility attributes.
Expert Elicitation	Delphi Method or Structured Expert Judgment Protocol	Formal process for eliciting, discussing, and reaching consensus on attribute scores from a diverse expert panel [3] [9].
	Data Quality Scoring Rubric (e.g., Good/Limited/No Data) [3]	Standardized metric to document and propagate uncertainty associated with each attribute score.
Modeling & Visualization	PSA Calculation Software (e.g., R packages psa, custom scripts)	Automates calculation of Productivity, Susceptibility, and composite Vulnerability scores from attribute matrices.
	Visual Analytics Platform for Spatiotemporal Data [12]	Enables interactive exploration of movement and fishery interaction data for advanced bycatch susceptibility analysis.
	Colorblind-Friendly Plotting Palette (e.g., Tableau's "Color Blind" palette) [13]	Ensures accessibility of published PSA graphs and charts for all audiences, adhering to WCAG guidelines [14] [15].

PSA Methodology in Action: Step-by-Step Application and Case Studies from Global Fisheries

Defining the assessment scope is the foundational step in fisheries science, determining the feasibility, relevance, and ultimate utility of the research. This process involves the systematic selection and bounding of the species, stocks, and fisheries to be evaluated, particularly within the context of a Productivity and Susceptibility Analysis (PSA). PSA is a semi-quantitative risk assessment method used to evaluate the vulnerability of a stock to overfishing when traditional, data-intensive stock assessments are not available [1]. The scope must be carefully circumscribed to balance scientific rigor with practical constraints, ensuring that management and conservation priorities are effectively addressed [16] [17].

Core Principles and Criteria for Selection

The selection process is guided by a combination of ecological, socio-economic, and governance criteria. A holistic system view is essential, considering not only the target species but also the associated habitats, fishing fleets, dependent communities, and regulatory frameworks [17]. The primary goal is to identify units where assessment will yield the greatest impact for sustainable management and conservation.

Key selection criteria include:

• Management Priority: Stocks subject to existing fishery management plans or legal mandates (e.g., under the Magnuson-Stevens Act) [16].
• Perceived Risk: Species exhibiting signs of decline, those with high economic value driving fishing pressure, or those with life history characteristics (e.g., low productivity) that make them inherently vulnerable [3].
• Data Availability: The presence of historical catch data, biological information, or abundance indices that can inform an assessment, while also recognizing that PSA is specifically designed for data-poor situations [3] [1].
• Stakeholder Concern: Fisheries of high social, economic, or cultural importance to communities and other stakeholder groups [17].

Comparative Analysis of Assessment Methodologies

The choice of assessment methodology is directly influenced by the defined scope and the data available for the selected stock. The following table compares the core approaches.

Table 1: Comparison of Fisheries Assessment Methodologies

Methodology	Description	Data Requirements	Best Suited For	Outputs
Full Stock Assessment	A quantitative analysis using mathematical models to estimate stock biomass and fishing mortality [16].	High-quality, time-series data on catch, abundance, and biology (age, growth, reproduction) [16].	Data-rich, high-value, or heavily exploited stocks in developed management systems.	Estimates of stock size, overfishing status, sustainable catch limits (quotas).
Productivity & Susceptibility Analysis (PSA)	A semi-quantitative, risk-based framework scoring a stock's inherent productivity and its susceptibility to fisheries [1].	Life history parameters (productivity) and fishery interaction details (susceptibility). Can incorporate expert judgment where data is limited [3].	Data-poor stocks, multispecies assessments, preliminary risk screening, and small-scale or emerging fisheries [3].	A vulnerability score (1-3) and ranking, identifying high-risk stocks for priority management.
Population Trend Analysis	Monitoring changes in relative abundance or catch-per-unit-effort (CPUE) over time.	Fishery-dependent catch data or fishery-independent survey indices.	Stocks with consistent monitoring data but insufficient information for full modeling.	Indicators of population increase, decrease, or stability.
Ecosystem-Based Assessment	Evaluates the stock within the context of its ecological community and habitat.	Data on species interactions, habitat dependencies, and environmental drivers.	Stocks where bycatch, trophic interactions, or habitat impacts are a major management concern.	Integrated advice considering species interactions and broader ecosystem health.

Application Notes: Implementing a PSA Framework

PSA is particularly valuable for the initial scoping and prioritization phase. Its implementation involves two parallel tracks of evaluation [1].

Productivity Attributes reflect the biological capacity of the stock to recover from depletion. Common attributes include:

• Average Age at Maturity: Earlier maturity generally indicates higher productivity.
• Fecundity: Number of offspring produced.
• Average Maximum Size: Smaller maximum size often correlates with faster growth and higher turnover.
• Natural Mortality Rate: Higher natural mortality suggests a stock is adapted to higher population turnover.

Susceptibility Attributes measure how exposed and vulnerable the stock is to the fishery. These include:

• Spatial Overlap: The degree to the stock's range overlaps with fishing grounds.
• Seasonal Overlap: Concurrent timing of fishing activity and stock aggregation (e.g., spawning).
• Gear Selectivity & Encounterability: How effectively the fishing gear catches the species and the stock's ability to avoid it.
• Post-Capture Mortality: Survival rate after escape or release from gear.

Table 2: PSA Attribute Scoring Framework (Simplified Example)

Attribute Category	Specific Attribute	Scoring Scale (1-3)	Data Quality Indicator
Productivity	Average Age at Maturity	1=Low (<2 yrs), 2=Medium (2-5 yrs), 3=High (>5 yrs)	Empirical, Published, Expert Judgment, Unknown
Productivity	Fecundity	1=High (>1M eggs), 2=Medium, 3=Low (<10,000 eggs)	Empirical, Published, Expert Judgment, Unknown
Susceptibility	Spatial Overlap	1=Low (<25% overlap), 2=Medium, 3=High (>75% overlap)	Empirical, Published, Expert Judgment, Unknown
Susceptibility	Gear Encounterability	1=Low (easily avoids gear), 2=Medium, 3=High (cannot avoid)	Empirical, Published, Expert Judgment, Unknown

Each attribute is scored from 1 (low risk/vulnerability) to 3 (high risk/vulnerability). The scores are then aggregated (often as a Euclidean distance from the origin) to produce an overall Vulnerability Score [3] [1].

Detailed Experimental Protocol: Conducting a PSA

Protocol Title: Semi-Quantitative Risk Assessment for Data-Poor Stocks Using Productivity and Susceptibility Analysis (PSA)

Objective: To determine the relative vulnerability to overfishing for one or more fish stocks where data are insufficient for a full quantitative stock assessment.

Materials & Pre-Assessment Preparation:

• Stock Identification: Clearly define the assessment unit (e.g., genetic stock, management stock, geographic population) [17].
• Fishery System Mapping: Document the relevant fisheries, gear types, seasons, and key spatial areas [17].
• Literature Review: Compile all available biological, ecological, and fisheries data for the target stock and similar species.
• Expert Identification: Assemble a panel of experts (scientists, fishery managers, experienced fishers) to provide informed judgment for missing data.

Procedure:

• Attribute Selection: Adopt or adapt a standardized list of productivity (P) and susceptibility (S) attributes (e.g., 10 P and 12 S attributes as per common frameworks) [1].
• Data Compilation: For each attribute, gather all available empirical data. Document the source and reliability.
• Scoring Workshop: Convene the expert panel.
  • Present compiled data for each attribute.
  • For attributes with robust data, guide the panel to assign a consensus score (1-3).
  • For attributes with poor or no data, facilitate a structured expert elicitation to assign a plausible score and a Data Quality rank (e.g., Empirical, Published, Expert, Unknown).
• Vulnerability Calculation:
  • Calculate the average score for all Productivity attributes (P̄).
  • Calculate the average score for all Susceptibility attributes (S̄).
  • Compute the Vulnerability Score (V) for each stock: V = √(P̄² + S̄²).
  • Plot results on a bivariate plot (P̄ on Y-axis, S̄ on X-axis). Quadrants of the plot represent different risk categories.
• Sensitivity & Uncertainty Analysis:
  • Test the effect of weighting different attributes.
  • Examine how scores change if expert-derived values are replaced with default or proxy values.
• Reporting: Document scores, data quality, expert rationale, and final vulnerability rankings. Contrast results with any independent information on stock status for validation [3].

Case Study: Application in a Data-Poor Context

A 2024 study applied an adapted PSA to ten coastal groundfish species in Peru, a classic data-poor environment [3]. The research streamlined the process to integrate multiple gears into a single susceptibility score per attribute.

Table 3: PSA Results for Selected Peruvian Coastal Groundfish [3]

Species (Common Name)	Productivity Score (P̄)	Susceptibility Score (S̄)	Vulnerability (V)	Risk Category
Broomtail Grouper	2.4	1.9	2.57	Extremely High
Grape-eye Seabass	2.3	2.0	2.50	Extremely High
Pacific Goliath Grouper	2.5	1.6	2.28	Extremely High
Black Snook	2.1	1.9	2.13	High
Mulata	1.8	1.9	1.89	Medium
Bumphead Parrotfish	1.7	1.6	1.73	Low

Key Findings & Validation: The assessment found data quality to be between "limited" and "no data," underscoring the reliance on expert knowledge [3]. The high-vulnerability species identified (e.g., groupers) matched general ecological understanding of slow-growing, late-maturing reef fish being at greater risk. This consistency with external knowledge validated the PSA thresholds and confirmed the method's utility for prioritization in management-deficient regions [3].

Visualizing the Process: Workflows and Logic

Scope Definition for Fisheries Assessment

PSA Scoring and Vulnerability Calculation Workflow

Stock Assessment Prioritization Decision Logic

The Scientist's Toolkit: Essential Reagents & Materials

Table 4: Key Research Reagent Solutions for Fisheries Assessment

Item/Category	Function/Description	Application in Scope Definition & PSA
Fishery-Dependent Data	Logbooks, landing receipts, and observer records providing catch (landings and discards) by species, area, and gear.	Fundamental for defining fishery extent, effort distribution, and scoring susceptibility attributes like spatial/seasonal overlap [16].
Fishery-Independent Survey Data	Scientifically designed surveys (trawl, acoustic, aerial) to estimate abundance, distribution, and biological parameters independent of fishing.	Provides critical biomass indices and biological data (size, age, maturity) for productivity scoring and trend analysis [16].
Biological Sample Kits	Tools for collecting otoliths (age determination), gonads (maturity staging), tissue samples (genetics, stable isotopes), and stomachs (diet analysis).	Generates core data for productivity attributes (growth, maturity, fecundity) and understanding stock structure [18].
Expert Elicitation Protocols	Structured frameworks (e.g., Delphi method, questionnaires) to formally gather and quantify expert judgment.	Essential for scoring PSA attributes in data-poor contexts and validating/scoring socio-economic factors [3].
Geographic Information System (GIS)	Software for spatial analysis and visualization.	Used to map fishing grounds, stock distributions, habitat features, and human uses to quantify spatial overlap and define system boundaries [17].
Stock Assessment & Statistical Software	Specialized platforms (e.g., R packages stockassessment, FishLife, JABBA) for data analysis and modeling.	Used for data exploration, running PSA calculations, and conducting more complex assessments if data allows [16] [3].
Stakeholder Engagement Materials	Facilitation guides, multi-lingual surveys, and participatory mapping tools.	Critical for inclusive system specification, understanding socio-economic dimensions, and validating fishery operational details [17].

Within the framework of Productivity and Susceptibility Analysis (PSA) for fisheries, the phase of data sourcing and scoring is a critical methodological bridge. It transforms raw, often disparate, information into structured, scorable parameters that assess a stock's vulnerability. This process systematically moves from gathering existing knowledge through comprehensive literature reviews to generating new, consensus-based data via formal expert elicitation. The integrity of the entire PSA hinges on the rigor, transparency, and repeatability of this step, ensuring that vulnerability scores are grounded in the best available evidence and expert judgment, thereby supporting sustainable fisheries management and certification processes [19].

Application Notes

Application Note: Conducting a Structured Literature Review for PSA Parameters

Purpose: To systematically identify, collate, and synthesize quantitative and qualitative data for predefined PSA productivity and susceptibility attributes (e.g., intrinsic rate of increase, area occupied, management effectiveness) from published and gray literature.

Protocol Overview: The review follows a modified systematic review process tailored for rapid yet comprehensive evidence gathering. It involves protocol registration, systematic searching across academic and institutional databases, and structured data extraction into a standardized template. The focus is on relevance to the specific fishery stock and geographic region, with data quality assessed based on source credibility and methodological rigor.

Key Considerations:

• Source Hierarchy: Prioritize peer-reviewed literature, followed by technical reports from recognized bodies (e.g., FAO, RFMOs), and finally other gray literature [19].
• Data Gaps Identification: Clearly document parameters with insufficient or conflicting data; these are prime candidates for resolution via expert elicitation.
• Version Control: Maintain a detailed log of the search strategy, inclusion/exclusion decisions, and extracted data with protocol versioning to ensure auditability [20].

Application Note: Developing and Populating a PSA Data Scoring Rubric

Purpose: To translate literature-derived data into consistent, ordinal scores (e.g., Low, Medium, High) for each PSA attribute using a predefined, objective rubric.

Protocol Overview: Scoring is performed by multiple independent reviewers. The rubric defines clear thresholds (e.g., metric values, qualitative descriptions) for each score level per attribute. Reviewers extract relevant data points from the synthesized literature review output and assign scores. Discrepancies are flagged, and a consensus meeting is held to resolve differences, documenting the rationale for final scores.

Key Considerations:

• Rubric Calibration: Conduct initial calibration exercises with the review team using sample data to ensure consistent interpretation of scoring thresholds.
• Handling Uncertainty: The rubric must include guidance for scoring when data is uncertain, incomplete, or conflicting, potentially using a "precautionary" default or flagging for expert input.
• Data Structuring: Extracted data and scores must be structured in a "tidy" format, where each row represents a unique stock-attribute observation, facilitating analysis [21].

Purpose: To obtain calibrated, quantitative judgments for PSA parameters where literature data is absent, highly uncertain, or contradictory.

Protocol Overview: This formal protocol, inspired by clinical trial standards, ensures the elicitation is unbiased and reproducible [20]. It involves selecting a diverse panel of experts, training them on the elicitation method and PSA context, conducting individual and/or group elicitation sessions (e.g., using the Delphi method or modified nominal group technique), and finally aggregating the judgments into consensus distributions or values.

Key Considerations:

• Expert Selection: Define clear expertise criteria, seek balanced representation (e.g., academia, industry, government), and manage conflicts of interest [20].
• Elicitation Method: Choose a method (e.g., 4-step elicitation, Delphi) that minimizes cognitive biases and group dominance. The SPIRIT framework emphasizes pre-defining the objectives, methods, and analysis plan for such "trials" [20].
• Documentation and Transparency: Record all materials, responses, and aggregation procedures. The protocol, like a trial protocol, should be accessible, detailing how the public or stakeholders were involved in design or review [20].

Detailed Experimental Protocols

Protocol: Systematic Data Acquisition and Structuring

Objective: To gather and organize all relevant data for a PSA from available literature into an analyzable format.

• Define the Data Schema: Create a master table with columns for: Stock ID, PSA Attribute, Data Value, Data Type (continuous, ordinal, categorical), Units, Source Citation, Source Type (journal, report, etc.), and Relevance Score.
• Execute Search Strategy:
  • Develop search strings using keywords related to the target species/fishery and each PSA attribute.
  • Search academic databases (e.g., Web of Science, Scopus), institutional repositories (e.g., FAO, IUCN), and fishery management body websites.
  • Screen titles/abstracts against inclusion criteria (e.g., geographic region, species, date range).
• Extract and Structure Data:
  • For each included document, extract relevant data points into the predefined schema.
  • Convert all data to consistent units.
  • For qualitative information (e.g., "management is considered moderately effective"), note the verbatim description and source context.
• Perform Data Cleaning and Validation:
  • Identify and correct outliers or entry errors by cross-referencing sources.
  • Classify the granularity of each data row (e.g., individual specimen data, stock-level summary) [21].
  • Generate summary statistics and distributions (e.g., using histograms) for each continuous attribute to visualize data range and central tendency [22].

Objective: To derive robust, consensus estimates for specific, data-poor PSA parameters through structured expert judgment.

• Protocol Finalization and Registration:
  • Finalize the elicitation protocol detailing background, objectives, elicitation method, expert panel composition, and statistical analysis plan. Reference the SPIRIT 2025 checklist for item completeness [20].
  • Register the protocol on an open-science platform or repository to ensure transparency [20].
• Expert Preparation and Training:
  • Recruit experts based on pre-defined criteria. Provide them with a background dossier containing the PSA framework, literature review summary for the target parameters, and definitions.
  • Conduct a training session on the elicitation format (e.g., estimating 5th, 50th, and 95th percentiles) and common cognitive biases.
• Elicitation Session:
  • Individual Elicitation (Round 1): Experts provide their judgments independently and anonymously, along with rationale.
  • Controlled Feedback (Round 2): Provide anonymized summary statistics (e.g., aggregate distribution) and rationales from Round 1. Experts may revise their estimates.
  • Consensus Discussion (Optional Round 3): For unresolved disparities, convene a structured, facilitated discussion focusing on reasons for differences, not debate.
• Data Aggregation and Documentation:
  • Aggregate final expert estimates using a pre-specified method (e.g., equal-weighted linear pooling).
  • Document the entire process, including all communications, raw elicitations, and final aggregated values, adhering to open science data sharing principles [20].

Data Presentation and Scoring Tables

Table 1: Categorization and Sources of Quantitative Data for PSA

Data Category	Description	Common Sources	Example PSA Attributes	Typical Format
Biological/Life History	Intrinsic species traits affecting population growth.	Peer-reviewed journals, stock assessment reports.	Age at maturity, fecundity, intrinsic rate of increase (r).	Continuous metrics (e.g., r = 0.28 yr⁻¹).
Spatial-Distribution	Metrics related to the geographic range and occupancy of the stock.	Species distribution models, fishery logbooks, telemetry studies.	Area occupied, overlap with fishery, depth range.	Area (km²), presence/absence grids.
Fishery-Dependent	Data collected directly from fishing activities.	Catch/landing statistics, observer programs, vessel monitoring.	Catch per unit effort (CPUE), selectivity, discard rate.	Time series, ratios.
Management & Compliance	Metrics reflecting the regulatory environment and its implementation.	Fishery management plans, compliance review reports, expert surveys [19].	Management strategy evaluation (MSE) use, enforcement capacity.	Ordinal scores (e.g., Low/Med/High), binary.

Table 2: Scoring Rubric for a Sample PSA Susceptibility Attribute ("Management Effectiveness")

Score (1-3)	Qualitative Description	Quantitative Indicators (if available)	Literature/Evidence Triggers
1 (Low Susceptibility)	Robust, adaptive management system with high compliance.	Formal MSE in place [19]; Compliance rate >90%; Recovery plans for depleted stocks.	Reports citing "effective enforcement" and "stakeholder compliance."
2 (Medium Susceptibility)	Basic management measures exist but are static or compliance is variable.	Total Allowable Catch (TAC) set but not based on formal harvest strategy; Compliance rate 70-90%.	Literature notes "management exists but lacks adaptability" or "partial compliance."
3 (High Susceptibility)	Inadequate or absent management, with low enforcement.	No catch controls; No monitoring; Compliance rate <70%.	Reports state "open access" or "widespread illegal fishing."

Table 3: SPIRIT 2025 Checklist Items Adapted for Expert Elicitation Protocol Design [20]

SPIRIT Section	Item No.	Original Item Description	Adaptation for Expert Elicitation Protocol
Administrative	1a	Title stating trial design, population, interventions.	Title stating elicitation method, expert panel, target parameters.
Open Science	4, 5, 6	Trial registration; Protocol access; Data sharing.	Register protocol; Make protocol & final data openly accessible.
Introduction	9a, 10	Scientific background; Specific objectives.	PSA background & rationale; Specific parameters for elicitation.
Methods	11	Patient/public involvement.	Stakeholder involvement in protocol review or question formulation.
Methods	12, 18	Trial design; Outcomes.	Elicitation design (e.g., Delphi); Target estimands (e.g., median, range).
Methods	22, 26	Statistical methods; Data management.	Methods for aggregating expert judgments; Data handling plan.

Visualizations and Workflow Diagrams

PSA Data Sourcing and Scoring Workflow

Methodology for Translating Data to PSA Scores

The Scientist's Toolkit: Essential Research Reagent Solutions

Table 4: Key Tools and Materials for PSA Data Sourcing and Scoring

Tool/Material Category	Specific Item/Software	Function in PSA Data Workflow
Literature Management	Reference Managers (e.g., Zotero, EndNote)	Centralized storage, deduplication, and citation of sourced literature; enabling team sharing.
Systematic Review	Covidence, Rayyan	Facilitating the blind screening of titles/abstracts and full texts during the structured literature review phase.
Data Structuring & Analysis	R, Python (pandas), Tableau Prep [21]	Cleaning, transforming, and structuring raw extracted data into a tidy format for analysis; generating summary statistics and visual distributions [22].
Expert Elicitation Platforms	Elicit.org, DelphiManager, Custom Google Forms	Conducting anonymous online elicitation rounds, aggregating responses, and providing controlled feedback to experts.
Consensus Facilitation	Video Conferencing (e.g., Zoom, Teams), Virtual Whiteboards (e.g., Miro)	Hosting structured discussion rounds for expert panels, especially in remote settings, to resolve scoring discrepancies.
Data Visualization & Accessibility	Graphviz (DOT language), Color Contrast Checkers [23]	Creating clear workflow and methodology diagrams; ensuring all visual elements meet WCAG 2.1 AA contrast requirements (≥4.5:1 for text, ≥3:1 for graphics) [24] [25] for inclusive science communication.
Protocol & Documentation	SPIRIT 2025 Checklist [20], Electronic Lab Notebooks	Ensuring expert elicitation protocols are comprehensive and transparent; maintaining an auditable record of all decisions and data transformations.

Productivity-Susceptibility Analysis (PSA) provides a semi-quantitative, risk-based framework for assessing the vulnerability of fish stocks to overfishing, particularly in data-limited scenarios [3]. The core analytical step translates biological productivity and fishery susceptibility attributes into a single, comparable vulnerability score. This score facilitates the prioritization of management actions by ranking stocks according to their relative risk [3] [11].

The fundamental formula for calculating the vulnerability (V) of a stock in a PSA is a Euclidean distance measure from the origin in a two-dimensional space defined by productivity and susceptibility:

V = √(P² + S²)

Where:

• P = Productivity Score. A composite score derived from life history traits that determine a population's intrinsic capacity to recover from depletion (e.g., fecundity, age at maturity, growth rate).
• S = Susceptibility Score. A composite score representing the stock's exposure to and likelihood of capture by a fishery, based on attributes like availability, encounterability, and selectivity of gears [3].

This geometric calculation treats productivity and susceptibility as orthogonal axes. A higher score in either dimension increases the overall vulnerability distance. The result is a continuous score, typically bounded between 1 (lowest vulnerability) and 3 (highest vulnerability), which is then interpreted using established thresholds to categorize risk (e.g., low, medium, high, very high) [3].

Advanced Computational Integration For enhanced predictive power, PSA frameworks can be integrated with data-driven mathematical models. A seminal approach involves using longitudinal data (e.g., time-series catch or biomarker data) to calibrate model parameters through profile-likelihood analysis, which determines identifiable parameters and their confidence intervals from available data [26]. This calibration allows for the generation of virtual populations or "clones," enabling the prediction of outcomes like treatment failure in medical contexts or stock collapse in fisheries, and the analysis of progression-free survival (PFS) within heterogeneous cohorts [26]. Furthermore, Principal Component Analysis (PCA) can be employed to reduce dimensionality and reveal the dominant biological and fishery factors (e.g., net expression rates of key biomarkers, growth rates) that distinguish populations within a cohort, thereby informing which parameters are most critical for accurate vulnerability assessment [26].

Experimental Protocols for PSA Data Collection and Analysis

Protocol for Data Collection in Data-Poor Fisheries

This protocol outlines a standardized method for gathering and preparing data for a PSA in contexts where formal stock assessments are not feasible [3].

Objective: To systematically collect and codify information on productivity and susceptibility attributes for target species. Materials: Species checklist, data extraction forms, access to scientific literature (e.g., FishBase, peer-reviewed journals), expert elicitation questionnaire. Procedure:

• Species Selection & Scoping: Define the list of species or stocks to be assessed. For broad-scale analyses, this may include the top 250 species by trade volume to capture >90% of fishery activity [11].
• Productivity Attribute Scoring:
  • For each species, search primary and gray literature for quantitative estimates of key life-history parameters: maximum age, von Bertalanffy growth coefficient (K), natural mortality (M), age at maturity, and fecundity.
  • Where species-specific data are absent, use the phylogenetic extrapolation method: fill data gaps with values from congeneric or confamilial species, explicitly documenting the source and taxonomic level of extrapolation [11].
  • Score each parameter using a predefined scoring system (e.g., 1-3, with 1 being low productivity/high resilience). Apply a data quality score (e.g., "good," "limited," "no data") to each attribute to quantify uncertainty [3].
• Susceptibility Attribute Scoring:
  • Map all fisheries and gears interacting with the stock. For each gear type, assess three sub-attributes:
    • Availability: Spatial/temporal overlap between stock and fishery.
    • Encounterability: Probability of capture given co-occurrence.
    • Selectivity: Size/age selection characteristics of the gear.
  • Score each sub-attribute (e.g., 1-3, with 3 being high susceptibility). Integrate scores from multiple gears into a single composite score per sub-attribute, often using a maximum score approach (worst-case) or a weighted average based on catch proportion [3].
• Expert Elicitation: Convene a panel of fisheries biologists and local fishing experts. Present compiled data and scores for review. Use a structured, anonymous voting or discussion process to reconcile scores where literature is ambiguous or absent. Record final agreed-upon scores and rationales.

Protocol for Computational Vulnerability Analysis and Clustering

This protocol details the steps for calculating vulnerability scores and analyzing patterns across multiple species using statistical clustering [26] [11].

Objective: To compute composite vulnerability scores and identify groups of species with similar risk profiles. Software Requirements: Statistical computing software (R, Python), PSA calculation script, clustering library (e.g., mclust for R). Procedure:

• Data Preparation: Compile finalized productivity (P) and susceptibility (S) attribute scores into a matrix format (species × attributes). Impute any remaining minor gaps using median values from the respective attribute column.
• Composite Score Calculation:
  • Calculate the mean score for all productivity attributes to derive the final Productivity Score (P) for each species.
  • Calculate the mean score for all susceptibility attributes to derive the final Susceptibility Score (S) for each species.
• Vulnerability Score Calculation: For each species, compute the vulnerability score using the Euclidean distance formula: V = √(P² + S²).
• Threshold Classification: Categorize species into vulnerability classes based on pre-defined thresholds (e.g., Low: V < 2.0; Medium: 2.0 ≤ V < 2.2; High: 2.2 ≤ V < 2.5; Very High: V ≥ 2.5) [3].
• Gaussian Mixture Model Clustering:
  • To visualize comparative sustainability and identify natural groupings not constrained by linear thresholds, apply a Gaussian Mixture Model (GMM) clustering algorithm to the bi-variate data (P, S) [11].
  • Standardize P and S scores. Determine the optimal number of clusters using the Bayesian Information Criterion (BIC).
  • Fit the GMM and assign each species a cluster membership probability.
  • Visualize results as an ellipsoid cluster plot overlaid on the PSA matrix.
• Uncertainty Propagation: Incorporate data quality scores into the final presentation. Species assessed with "limited" or "no data" scores should have their vulnerability scores visually flagged (e.g., with hatched patterns or confidence intervals) on the PSA plot to reflect higher uncertainty [3] [11].

Data Presentation and Results

Table 1: Productivity and Susceptibility Attributes with Scoring Criteria This table defines the common attributes used in fishery PSA and their scoring logic, where a lower score (1) indicates higher resilience/lower risk.

Dimension	Attribute	Score = 1 (Low Risk/High Resilience)	Score = 2 (Moderate)	Score = 3 (High Risk/Low Resilience)
Productivity	Maximum Age	Shorter-lived (< 5 yrs)	Medium (5-15 yrs)	Longer-lived (> 15 yrs)
	Growth Rate (K)	Fast (> 0.3)	Medium (0.1-0.3)	Slow (< 0.1)
	Natural Mortality (M)	High (> 0.2)	Medium (0.1-0.2)	Low (< 0.1)
	Age at Maturity	Early (< 2 yrs)	Medium (2-5 yrs)	Late (> 5 yrs)
	Fecundity	High (> 10^6 eggs)	Medium (10^4 - 10^6)	Low (< 10^4 eggs)
Susceptibility	Availability	No spatial/temporal overlap	Seasonal/partial overlap	Full spatial/temporal overlap
	Encounterability	Avoids gear (e.g., habitat)	Neutral behavior	Attracted to/does not avoid gear
	Selectivity	Gear selects non-target sizes	Catches some target sizes	Gear perfectly targets vulnerable sizes
	Post-Capture Mortality	Survival after release likely	Mixed survival	Total mortality (e.g., bottom trawl)

Table 2: Illustrative PSA Results for Peruvian Coastal Groundfish Adapted from a study assessing data-poor stocks, this table shows final composite scores and vulnerability classifications [3].

Species	Productivity Score (P)	Susceptibility Score (S)	Vulnerability (V)	Classification	Data Quality
Broomtail Grouper	2.4	1.0	2.57	Very High	Limited
Pacific Goliath Grouper	2.2	0.8	2.28	Very High	Limited-No Data
Black Snook	1.8	1.2	2.13	High	Limited
Mulata	1.6	1.0	1.89	Medium	Limited
Bumphead Parrotfish	1.4	0.9	1.73	Low	Limited

Visualizing the PSA Workflow and Matrix

PSA Vulnerability Calculation Workflow

Plotting Results on the PSA Matrix

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Reagents and Materials for PSA Research

Item	Function/Description	Application in PSA Protocol
Life-History Database Access	Subscription or access to curated biological trait databases (e.g., FishBase, SeaLifeBase).	Provides foundational data for scoring productivity attributes; critical for data-poor species via phylogenetic extrapolation [11].
Expert Elicitation Framework	Structured protocol (Delphi method, anonymous voting software) for formalized expert judgment.	Standardizes the reconciliation of ambiguous data and scoring of attributes where empirical information is lacking [3].
Statistical Software Suite	Platforms with PSA, clustering, and visualization packages (R with fishmethods/mclust, Python with sklearn/MDAnalysis).	Performs core vulnerability calculations, Gaussian Mixture Model clustering, and generates the PSA matrix plots [26] [11] [27].
Uncertainty Quantification Module	Scripts to implement data quality scoring and propagate uncertainty through calculations.	Assigns and visualizes confidence levels to final vulnerability scores, ensuring transparent communication of assessment reliability [3] [11].
Path Similarity Analysis (PSA) Toolkit	Computational framework (e.g., MDAnalysis.psa) for quantifying geometric similarity between trajectories.	While developed for molecular dynamics, its metrics (Hausdorff/Fréchet distances) and clustering approaches offer analogies for comparing stock recovery or depletion pathways across different management scenarios [27].

Coastal groundfish stocks in Peru represent a critical component of both marine biodiversity and socio-economic stability for local communities. Species such as the Peruvian grunt (Anisotremus scapularis) are highly valued in national markets but face significant threats from overexploitation and illegal, unreported, and unregulated (IUU) fishing [28]. The assessment and management of these species are hindered by a classic "data-poor" scenario, characterized by irregular landings, limited formal stock assessments, and high dependency on artisanal fishing sectors [28]. This context creates an urgent need for innovative assessment frameworks that can function effectively with limited data.

This case study is framed within the broader thesis on Productivity and Susceptibility Analysis (PSA) for fisheries research. PSA provides a risk-based framework ideal for data-poor situations, as it systematically evaluates a species' inherent productivity (biological capacity to withstand fishing pressure) against its susceptibility to various fishery-related threats. The application of PSA to Peruvian groundfish, supplemented by emerging technologies and cooperative management, offers a pathway to sustainable management where traditional, data-intensive stock assessment models are not feasible.

Peruvian Grunt (Anisotremus scapularis) Fishery and Biology

Table 1: Key biological and fishery data for the Peruvian Grunt.

Parameter	Value / Description	Source / Notes
Common Names	Chita, Sargo, Roncador, Corcovado	[28]
Geographic Range	Ecuador (Manta) to Chile (Antofagasta & Cocos Island)	[28]
Habitat	Rocky areas, intertidal pools (juveniles); sandy bottoms (adults); depth down to 25m	[28]
Maximum Size	80 cm length, 7.5 kg weight	[28]
Diet	Juveniles: Omnivorous (macroalgae, copepods). Adults: Euryphagous (algae, amphipods, mollusks)	[28]
Annual Landings (2017-2021)	119 - 144 Gross Metric Tons (GMT)	Highly irregular, among smallest Peruvian artisanal fisheries [28]
Fishing Gear	Artisanal (string, seine, curtain, spinel)	[28]
Economic Value	High market value, appreciated for organoleptic quality	[28]
Primary Threat	Overexploitation and IUU fishing	Status not fully evaluated [28]

Contextual Fishery and Environmental Data

Table 2: Regional fishery and environmental context.

Parameter	Value / Description	Relevance to Groundfish Assessment
Anchovy Biomass Survey Method	Hydroacoustics (echo sounders on research vessels)	Non-invasive method provides model for assessing data-poor groundfish stocks [29].
Key Ecosystem	Humboldt Current Large Marine Ecosystem (HCLME)	Highly productive system supporting biodiversity and fisheries [29].
Threat to Marine Reserves	Illegal fishing, use of explosives, hunting, sanitation issues	Highlights enforcement challenges that impact all fishery management, including for groundfish [30].
Major Climate Stressor	El Niño and La Niña phenomena	Alters species distribution and reduces catch, affecting stock assessments and community livelihoods [31].
Contribution of Fisheries to Peruvian GDP	~1.3%	Underscores economic importance of sustainable management [28].

Detailed Methodological Protocols

Protocol for Hydroacoustic Biomass Estimation (Adapted for Groundfish)

This protocol is adapted from the binational standardized methodology developed for anchovy assessment in the HCLME [32], demonstrating a transferable approach for data-poor groundfish stocks.

Objective: To estimate the biomass and spatial distribution of coastal groundfish stocks using non-invasive hydroacoustic methods.

Pre-Survey Preparation:

• Vessel and Equipment Calibration: Utilize a dedicated research vessel equipped with calibrated scientific echo sounders (e.g., 38, 120, 200 kHz). Calibration must be performed pre- and post-survey using standard sphere targets [32].
• Survey Design: Develop a systematic parallel transect grid covering the target area (e.g., rocky shore and shelf habitats). Transect spacing is determined based on the expected school patchiness of the target groundfish species.

Data Collection:

• Acoustic Data Acquisition: Continuously collect raw acoustic backscatter data (SV, volume backscattering strength) along all transects. Record precise GPS position, time, and environmental data (temperature, salinity) concurrently [29] [32].
• Biological Sampling (Groundtruthing): Conduct periodic mid-water or bottom trawls (as appropriate for demersal species) at predetermined stations aligned with acoustic transects. All catches are identified to species level. For target groundfish, collect biological samples including length, weight, sex, and maturity stage. Otoliths are extracted for age determination [32].
• Species Discrimination: Apply advanced acoustic processing to differentiate species based on frequency response (multi-frequency analysis), school morphology, and depth distribution. Match acoustic signatures to species composition confirmed by biological trawls [32].

Data Processing and Analysis:

• Echo Integration: Process acoustic data to integrate backscatter over defined depth layers and geographic intervals (Elementary Distance Sampling Units, EDSUs).
• Target Strength (TS) Application: Apply a species-specific TS-to-length relationship (e.g., TS = m log10(L) + b) to convert acoustic backscatter to fish density. For Peruvian grunt, a proxy relationship from a similar species may be used initially.
• Biomass Estimation: Calculate biomass (B) for each EDSU: B = (Mean Density) × (Survey Area Volume). Sum across all EDSUs for total area biomass estimate, with associated confidence intervals derived from geostatistical analysis.

Standardization: This protocol emphasizes alignment with neighboring institutions (e.g., Chile's IFOP) on key parameters: cruise timing, equipment specs, frequencies, and data processing algorithms to ensure comparable and shared stock assessments [32].

Experimental Protocol for Aquaculture-Based Research on Peruvian Grunt

Objective: To establish captive breeding and grow-out protocols for A. scapularis to support stock enhancement, provide alternative livelihoods, and reduce fishing pressure on wild stocks [28].

Broodstock Management:

• Collection & Acclimation: Wild adults are captured using hook-and-line or traps from rocky shore habitats. They are acclimated in large, flow-through seawater tanks with simulated natural photoperiod and shelter.
• Spawning Induction: Administer hormonal implants (e.g., LHRHa) to ripe females and males. Place induced broodstock in spawning tanks with gentle currents.
• Egg Collection & Incubation: Collect buoyant eggs from tank overflow. Incubate eggs in conical tanks with gentle aeration and UV-sterilized seawater at 18-20°C.

Larval Rearing:

• Initial Feeding: Upon yolk-sac absorption (3-4 days post-hatch), feed larvae with rotifers (Brachionus plicatilis) enriched with omega-3 fatty acids for 15-20 days.
• Weaning: Co-feed Artemia nauplii from day 10, gradually transitioning to inert microdiets by day 30.
• Environmental Conditions: Maintain water quality: DO > 5 mg/L, pH 7.8-8.2, ammonia < 0.1 mg/L.

Grow-Out in Recirculating Aquaculture Systems (RAS):

• System Setup: Utilize RAS with biofiltration, foam fractionation, and UV sterilization. Implement a 12-hour light:12-hour dark photoperiod.
• Stocking & Feeding: Stock juveniles at 5-10 kg/m³. Feed with high-protein formulated pellets (45-50% protein) to satiation twice daily.
• Monitoring: Record daily feed intake, weekly growth (weight, length), and mortality. System performance is tracked via a bioeconomic model integrating growth rates, feed conversion ratios, and market prices to assess economic viability [28].

The Scientist's Toolkit: Key Research Reagent Solutions

Table 3: Essential materials and reagents for groundfish assessment and related research.

Item Category	Specific Item / Technology	Function in Research
Field Survey Equipment	Scientific Echo Sounder (multi-frequency)	Transmits sound waves and measures returning echoes to detect and quantify fish biomass [29] [32].
Field Survey Equipment	Calibrated Trawl Nets	Captures biological samples for species identification, size distribution, and age structure to "groundtruth" acoustic data [32].
Laboratory Analysis	Otoliths (fish ear bones)	Used for aging fish via annuli counts, providing critical data on growth and mortality for stock assessment models.
Laboratory Analysis	Hormone Implants (e.g., LHRHa)	Induces spawning in captive broodstock for aquaculture and stock enhancement trials [28].
Live Feed Culture	Rotifers (Brachionus sp.) & Microalgae	Essential first food for marine fish larvae in hatchery production [28].
Aquaculture System	Recirculating Aquaculture System (RAS)	Technology for intensive, land-based fish farming with controlled water quality, used for grow-out studies and bioeconomic modeling [28].
Data Analysis	Bioeconomic Model	Integrates biological growth data with costs and market prices to assess the economic feasibility of aquaculture or management interventions [28].

Visualizing Analytical Frameworks and Workflows

PSA Framework for Data-Poor Groundfish Assessment

Diagram 1: PSA framework for data-poor fisheries.

Integrated Groundfish Assessment Workflow

Diagram 2: Integrated assessment workflow for data-poor stocks.

The Bay of Bengal represents a critical marine ecosystem where intense fishing pressure intersects with significant populations of elasmobranchs (sharks, rays, and skates). In Bangladesh, these species are primarily captured as bycatch in fisheries targeting high-value teleosts, such as the Indian threadfin (Leptomelanosoma indicum) using large-mesh gillnets known as Lakkha nets [33] [34]. This case study applies Productivity and Susceptibility Analysis (PSA), a cornerstone semi-quantitative method for ecological risk assessment in data-limited fisheries research, to evaluate the vulnerability of elasmobranch bycatch in this region [34] [35].

The broader ecosystem context underscores the urgency of this assessment. Recent surveys indicate a 78.6% decline in small pelagic fish stocks in the Bay of Bengal over seven years, from 158,100 tonnes in 2018 to 33,811 tonnes in 2025, primarily due to overfishing [36]. This precipitous drop signals severe ecosystem stress that directly impacts predators like elasmobranchs. Furthermore, stock assessments for aggregated elasmobranch catches in Bangladesh indicate overexploitation, with current biomass below sustainable levels (B/BMSY < 1.0) and fishing pressure exceeding sustainable limits (F/FMSY > 1.0) [37]. Within this stressed system, the PSA framework provides a critical tool for prioritizing management by identifying species most vulnerable to extirpation due to their life-history traits (productivity) and interaction with fisheries (susceptibility) [33] [35].

Application Notes: PSA Framework and Regional Adaptation

Core Principles and Rationale for Application

PSA is designed for fisheries where species-specific catch, abundance, and life-history data are sparse—a condition that defines the elasmobranch fishery in the Bay of Bengal [34] [37]. The method evaluates vulnerability (V) as a function of two components:

• Productivity (P): The capacity of a stock to withstand fishing pressure, based on life-history attributes like age at maturity, fecundity, and natural mortality.
• Susceptibility (S): The likelihood and intensity of a stock's encounter with fishing gear, based on factors like spatial overlap, seasonality, and gear selectivity [34] [35].

The advantage of PSA in this context is its ability to use expert elicitation, local fisher knowledge, and patchy biological data to derive relative risk scores, thereby informing conservation priorities without requiring intensive stock assessment data [34].

Bay of Bengal-Specific Context and Challenges

Applying PSA in this region involves navigating specific challenges:

• Data Poverty: For elasmobranch bycatch, 76% of species received a "poor" data quality score for productivity attributes, necessitating heavy reliance on regional analogues and expert judgment [33].
• Fishing Dynamics: The primary gear of concern is the Lakkha net, a large-mesh (300-400 mm) drift gillnet deployed for 5-6 hours in open sea [34]. Its non-selectivity makes it a major source of megafauna bycatch.
• Socioeconomic Drivers: Elasmobranchs have local market value and are part of the international fin trade, but are not primary targets. Conservation measures must therefore balance ecological risk with the livelihoods of artisanal fishers [38].
• Multiple Stressors: Beyond bycatch, elasmobranchs in the region face threats from habitat degradation and bioaccumulation of metals (e.g., arsenic, mercury), which may compound population vulnerabilities [39].

Species Composition and Catch Analysis

The foundational PSA study in the Bangladesh region of the Bay of Bengal identified 40 elasmobranch species as bycatch in the Lakkha net fishery [33] [34].

Table 1: Elasmobranch Bycatch Species Composition in Lakkha Net Fishery

Taxonomic Group	Number of Species	Number of Families	Key Orders/Notes
Sharks	13	3	Primarily pelagic species
Rays	27	6	Majority (60%) belong to Myliobatiformes (eagle rays, stingrays)
Total Elasmobranchs	40	9	All high-vulnerability species are IUCN threatened

Productivity and Susceptibility Scoring Results

Scores were assigned on a scale (commonly 1-3), with lower values indicating higher productivity (resilience) or lower susceptibility [34].

Table 2: PSA Scoring and Vulnerability Outcomes for Bay of Bengal Elasmobranchs

Metric	Score Range	Mean/Median	Key Findings
Productivity (P)	1.27 – 2.73	Not specified	Scores reflect generally low productivity (k-selected traits).
Susceptibility (S)	1.50 – 2.63	Not specified	Target fish (Indian threadfin) scored highest, followed by pelagic sharks/eagle rays.
Vulnerability (V)	Derived from P & S	Not specified	31.7% (n=13) species = High Vulnerability 43.9% (n=18) species = Moderate Vulnerability 24.4% (n=10) species = Low Vulnerability

Supporting Stock Assessment Data

Broader stock-level analysis corroborates the high risk to the elasmobranch population [37].

Table 3: Stock Assessment Parameters for Aggregated Elasmobranchs in Bay of Bengal, Bangladesh

Parameter	CMSY Model Estimate	BSM Model Estimate	Management Implication
Intrinsic Growth Rate (r)	0.282 year⁻¹	Derived from model	Low resilience, typical for elasmobranchs.
Carrying Capacity (k)	119,000 mt	134,000 mt	Potential unexploited biomass.
Maximum Sustainable Yield (MSY)	8,420 mt	5,110 mt	Sustainable catch limit.
2022 Catch	7,017 mt	7,017 mt	Exceeds BSM MSY, indicating overfishing.
Biomass Status (B/BMSY)	< 1.0	< 1.0	Stock is depleted.
Fishing Pressure (F/FMSY)	> 1.0	> 1.0	Fishing mortality is unsustainable.

Experimental Protocols

Protocol 1: Field Sampling and Bycatch Identification for PSA

Objective: To systematically collect, identify, and record elasmobranch bycatch from the Lakkha gillnet fishery to create a species-specific dataset for PSA [34].

• Study Site & Gear Selection:
  • Focus on major landing sites in Chattogram and Cox's Bazar districts, Bangladesh [34].
  • Sample catches from the Lakkha net (large-mesh drift gillnet, 300-400 mm mesh size, 5-6 hour soak time) [34].
• Data Collection:
  • Conduct random sampling at landing sites over a minimum two-year period to account for seasonal variation [34].
  • For each elasmobranch specimen, record: species, total length, disc width (rays), sex, and maturity state.
  • Document fishing effort data (net length, soak time, location) from fishers.
• Species Identification:
  • Use a combination of morphological keys, regional field guides, and photographic documentation.
  • Collect tissue samples (fin clips) from ambiguous specimens for future genetic barcoding to resolve taxonomic uncertainties [34].
• Data Management:
  • Compile data into a relational database linking species records with corresponding fishing effort and location data.

Protocol 2: Scoring Productivity and Susceptibility Attributes

Objective: To assign standardized scores to productivity and susceptibility attributes for each identified elasmobranch species [34] [35].

• Attribute Selection:
  • Productivity Attributes (7): Maximum age, size at maturity, fecundity, reproductive strategy, natural mortality rate, trophic level, and geographic range [34].
  • Susceptibility Attributes (5): Spatial overlap with fishery, seasonal overlap, encounterability, selectivity, and post-capture mortality [34] [35].
• Scoring Criteria Definition:
  • Establish a 3-point scale (e.g., 1=Low, 2=Medium, 3=High) for each attribute. For productivity, a low score (1) indicates high resilience (e.g., early maturation, high fecundity). For susceptibility, a low score (1) indicates low exposure to the fishery [35].
  • Define clear, evidence-based thresholds for each score based on literature, global databases (FishBase), and regional studies.
• Expert Elicitation Workshop:
  • Convene a panel of local fisheries scientists, biologists, and experienced fishers.
  • For each species and attribute, present available data. Where data is missing (a common occurrence), guide the panel to assign scores based on expert consensus and analogy with similar species [33].
• Score Calculation:
  • Calculate the geometric mean of scores for all productivity attributes to derive the final Productivity (P) score.
  • Calculate the geometric mean of scores for all susceptibility attributes to derive the final Susceptibility (S) score.
  • Compute the overall Vulnerability (V) score using the formula: ( V = \sqrt{(P^2 + S^2)} ) [34].

Protocol 3: Data Quality Assessment and Sensitivity Analysis

Objective: To evaluate the uncertainty in the PSA and test the robustness of vulnerability rankings [33].

• Data Quality Scoring:
  • Assign a data quality rank (e.g., High, Medium, Low, Poor) to each attribute score for every species, based on the source and reliability of the information used.
  • Summarize to identify major data gaps (e.g., "76% of species had poor data quality for productivity attributes") [33].
• Sensitivity Analysis:
  • Systematically vary the susceptibility scores (e.g., ± 0.5 points) to simulate changes in fishing practices or effort.
  • Recalculate the vulnerability score and classification (High/Moderate/Low) for each species under the altered scenarios.
  • Identify "threshold species" whose vulnerability classification changes with minor score adjustments, highlighting priorities for targeted research and precautionary management [33].

Visualizations

PSA Methodology Workflow for Elasmobranch Bycatch

PSA Workflow for Bycatch Risk Assessment

Determinants of Gillnet Susceptibility for Elasmobranchs

Determinants of Gillnet Bycatch Susceptibility

The Scientist's Toolkit: Research Reagent Solutions

Table 4: Essential Materials and Tools for Elasmobranch Bycatch PSA Fieldwork

Item Category	Specific Item/Reagent	Function in Research
Field Sampling & Morphology	Dial Calipers (300+ mm), Measuring Boards, Spring Scales	Accurately records morphometric data (length, disc width) and weight, essential for life-history parameter estimation [34].
Species Identification	Regional Taxonomic Keys (e.g., FAO Species Catalogues), Digital Camera, Tissue Sampling Kit (scalpel, forceps, vials, 95% ethanol)	Enables morphological identification and preserves tissue for genetic barcoding to resolve taxonomic uncertainties common in elasmobranch bycatch [34].
Data Recording	Waterproof Field Notebooks, GPS Unit, Barcoded Sample Tags	Ensures accurate, spatially-explicit metadata collection and maintains chain of custody for biological samples.
Biological Sampling	Dissection Kit (scalpels, scissors), Histology Cassettes, 10% Neutral Buffered Formalin	Allows for onboard collection of gonads, vertebrae, and other structures for age, growth, and reproductive analysis, filling critical data gaps [33].
Interview & Survey	Structured Questionnaire (translated to local language), Digital Audio Recorder	Captures local ecological knowledge (LEK) from fishers on species occurrence, seasonality, and fishing practices, a crucial data source in PSA [34].
Data Analysis	PSA Scoring Software (e.g., R package rPSA), Genetic Analysis Software (e.g., Geneious, BOLD Systems)	Facilitates standardized vulnerability score calculation and performs genetic analysis for species identification and population structure studies.

This protocol details the application of a Productivity-Susceptibility Analysis (PSA) framework to assess the vulnerability of sensitive marine megafauna to bycatch in Mediterranean multi-gear fisheries. The PSA is a semi-quantitative Ecological Risk Assessment (ERA) tool, designed for data-limited situations, that evaluates a species' or species group's inherent capacity to withstand fishing pressure (Productivity) against its likelihood of encountering and being captured by fishing operations (Susceptibility) [40]. The outcome is a Vulnerability score, which prioritizes species groups and fisheries for targeted research and management action [40]. This case study focuses on five key megafauna groups: marine mammals (MM), seabirds (SB), pelagic elasmobranchs (PE), demersal elasmobranchs (DE), and sea turtles (ST), across major Mediterranean and Black Sea fishing gears [40].

Materials and Methodology

Defined Study Parameters

Table 1: Defined Study Parameters for the PSA Framework [40]

Parameter Category	Specific Elements	Operational Definition & Notes
Study Area	Mediterranean Sea, Black Sea	Subdivided into FAO-GFCM sub-regions: Adriatic Sea, Central Med., Eastern Med., Western Med., Black Sea.
Target Species Groups	Marine Mammals (MM), Seabirds (SB), Pelagic Elasmobranchs (PE), Demersal Elasmobranchs (DE), Sea Turtles (ST)	78 relevant species from the FAO list of vulnerable species were categorized into these 5 functional groups for management-relevant analysis.
Assessed Fishing Gears	Bottom Trawl (OTB), Pelagic Trawl (TM), Drifting Longline (LLD), Set Longline (LLS), Set Nets (GEN), Purse Seine (PS)	Gears were selected based on documented potential for bycatch interaction. Other gears (e.g., dredges) were excluded due to minimal perceived risk.
Core PSA Components	Productivity (P): Biological capacity for recovery. Susceptibility (S): Exposure and capture probability. Vulnerability (V): Combined risk score (√(P² + S²)).	A semi-quantitative scoring system (1-3) is applied to attributes within P and S.

Core Experimental Protocol: PSA Scoring and Calculation

Step 1: Productivity (P) Scoring. For each species group, assign a score (1=Low, 2=Medium, 3=High) to each productivity attribute based on averaged life-history traits of representative species within the group [40]. The overall Productivity score is the average of the attribute scores.

• Key Attributes: Average age at maturity, average maximum size, average fecundity, natural mortality rate.

Step 2: Susceptibility (S) Scoring. For each combination of species group and fishing gear, assign a score (1=Low, 2=Medium, 3=High) to each susceptibility attribute based on literature review and expert judgment [40].

• Key Attributes: Spatial Overlap: Co-occurrence of species habitat and fishing effort. Temporal Overlap: Co-occurrence in time. Morphological/Behavioral Selectivity: Probability of capture given encounter (e.g., size selectivity, ability to avoid gear). Post-Release Mortality: Survival rate after escape or release.

Step 3: Vulnerability (V) Calculation. Calculate the two-dimensional Euclidean distance for each species group-gear combination [40]: V = √(P² + S²) The resulting score (range: √2 to √18) is classified as:

• Low Vulnerability: V ≤ 2.5
• Medium Vulnerability: 2.5 < V < 3.5
• High Vulnerability: V ≥ 3.5

Step 4: Data Integration and Hotspot Analysis. Spatially integrate PSA results with fishing effort data (e.g., from Vessel Monitoring Systems - VMS) and species distribution maps. Areas with high effort and high-vulnerability combinations are identified as bycatch risk hotspots [41]. This can be supplemented with Local Ecological Knowledge (LEK) from fisher questionnaires to ground-truth and refine hotspot maps [41].

Diagram 1: PSA Methodology and Risk Assessment Workflow (87 characters)

Supporting Field Protocol: Onboard Observer Data Collection

For empirical validation and calibration of susceptibility scores, systematic onboard monitoring is essential [42].

Procedure:

• Trip Selection: Implement a stratified random sampling design to cover different seasons, sub-regions, and vessel sizes. Target a minimum coverage of 0.5-1% of total fishing effort, with ideal coverages of 2-7% for robust estimation [40].
• Data Recording per Haul/Set:
  • Fishing Effort: Gear specifications, set/haul start and end times, coordinates, depth, speed.
  • Catch Examination: All megafauna bycatch is identified to species level.
  • Biological Data: Length, weight, sex determined where possible.
  • Fate Assessment: Record condition (alive/dead, injured) and release method [42].
• Data Management: Store data in a relational database (e.g., MySQL, PostgreSQL) with linked tables for trips, hauls, catch, and biological data to ensure integrity and facilitate analysis [42].

Results and Data Synthesis

PSA Vulnerability Outcomes by Region and Gear

Table 2: Synthesis of PSA Vulnerability Outcomes for Mediterranean Sub-Regions [40]

Species Group	High-Risk Gears & Regions	Medium-Risk Gears & Regions	Key Management Insights
Sea Turtles (ST)	LLD, OTB, GEN, LLS across most regions.	PS in specific coastal areas.	High risk pervasive across net and longline gears. Requires broad mitigation.
Pelagic Elasmobranchs (PE)	LLD (all regions), OTB (Adriatic, Central Med).	TM, GEN.	Drifting longlines are a critical threat to pelagic sharks.
Demersal Elasmobranchs (DE)	OTB, GEN, LLS (all regions).	TM.	Bottom-contact gears (trawls, nets) pose greatest threat to bottom-dwelling rays and sharks.
Marine Mammals (MM)	GEN (Black Sea - Harbour porpoise).	OTB, LLS, LLD across most regions.	Set nets are an acute threat to low-productivity populations. Depredation on LLD can lead to bycatch [43].
Seabirds (SB)	Generally low risk across gears.	Potential low-medium risk from LLD in specific areas.	Lower priority group for bycatch mitigation in Mediterranean.

Quantitative Bycatch and Discard Data from Case Studies

Table 3: Empirical Bycatch and Discard Metrics from Regional Case Studies

Metric	Value	Context / Fishery	Source
Average Discard Rate (Weight)	70.7%	Deep-sea bottom trawl fishery, Antalya Bay (Eastern Med)	[44]
Bycatch Rate (Cetaceans)	0.071 individuals/1000 hooks	Pelagic longline fishery, Cyprus (Eastern Med)	[43]
Reported Depredation Rate	Up to 100% (Mean: 17%) per trip	Pelagic longline (albacore tuna), Cyprus	[43]
Estimated Annual Economic Loss	~259,272 EUR	Pelagic longline fleet, Cyprus, due to cetacean depredation	[43]
Monitoring Coverage	3-7% of annual effort	Midwater pair trawl, Adriatic Sea (14-year program)	[42]

The Scientist's Toolkit: Research Reagent Solutions

Table 4: Essential Materials, Technologies, and Analytical Tools for PSA Implementation

Tool / Solution Category	Specific Item / Platform	Function in PSA and Bycatch Research
Field Data Collection	Standardized Observer Forms (Electronic/Paper)	Ensures consistent recording of fishing effort, bycatch counts, and biological data onboard vessels [42].
	Remote Electronic Monitoring (REM) Systems	Video/Audio systems to validate logbook data and observer reports, covering more fishing effort [40].
Spatial Analysis	Vessel Monitoring System (VMS) Data	Provides high-resolution fishing effort location data for spatial overlap analysis in Susceptibility scoring [45].
	Species Distribution Models (SDMs)	Predicts species habitat using environmental variables to map spatial overlap with fisheries.
Data Management & Analysis	Relational Database (e.g., MySQL, PostgreSQL)	Essential for storing, managing, and querying complex, linked fishery and bycatch datasets [42].
	R/Python Statistical Packages (dplyr, ggplot2, sf)	Used for data cleaning, PSA score calculation, visualization, and spatial analysis.
Alternative Data Sources	Local Ecological Knowledge (LEK) Questionnaires	Structured interviews with fishers to gather data on bycatch hotspots and seasonal patterns where observer coverage is low [41] [43].
	Fisher Logbooks (Standardized)	Fisher-reported data on catch and bycatch, useful for large-scale trend analysis when cross-validated [43].

Diagram 2: PSA Risk Logic and Conservation Impact (81 characters)

Advanced Protocol: Integrated Hotspot Mapping

This protocol combines PSA outputs with spatial data to guide area-based management [41].

Procedure:

• Data Layer Preparation:
  • Risk Layer: Rasterize PSA results (V scores) for each species group-gear combination.
  • Effort Layer: Aggregate VMS or logbook data into a spatial grid (e.g., 10km x 10km) showing fishing days or gear sets per unit area.
  • Habitat Layer: Incorporate known critical habitats (rookeries, foraging grounds) from tracking studies or models [41].
• Spatial Overlay Analysis: Use Geographic Information System (GIS) software to perform a weighted overlay of the layers. For example: Composite Risk Score = (V_Score * Weight_A) + (Fishing_Effort * Weight_B) + (Habitat_Value * Weight_C)
• Hotspot Validation: Cross-reference generated hotspot maps with Local Ecological Knowledge (LEK) from fisher surveys to confirm identified areas and identify drivers [41].
• Output: Generate maps identifying "Priority Intervention Zones" where high vulnerability, high effort, and critical habitat converge.

This protocol provides a standardized, semi-quantitative framework for assessing bycatch risk in complex, multi-gear fisheries. The PSA results prioritize mitigation efforts for high-vulnerability combinations, such as:

• Technical Measures: Implementing turtle excluder devices (TEDs) in trawls, using circle hooks and bird-scaring lines in longlines.
• Spatial-Temporal Closures: Enacting fishing restrictions in identified hotspot areas, particularly during key life-history periods (e.g., turtle nesting seasons).
• Monitoring Enhancement: Directing limited onboard observer resources to fisheries and regions flagged as high-risk by the PSA.

The integrated, multi-method approach—combining PSA, spatial analysis, and LEK—strengthens the scientific basis for evidence-based management by regional fisheries bodies, supporting the conservation of sensitive megafauna in the Mediterranean Sea.

Navigating PSA Challenges: Optimizing Data Quality, Scoring, and Knowledge Integration

The global fight against extreme poverty faces a critical juncture. After decades of progress, the rapid decline in the number of people living in extreme poverty—from 2.3 billion in 1990 to an estimated 808 million in 2025—is projected to stall and potentially reverse due to economic stagnation in the world's poorest regions [46]. This stagnation is acutely felt in fisheries-dependent communities, which are often characterized by high levels of poverty and a critical lack of data for sustainable management [47]. Within this context, Productivity and Susceptibility Analysis (PSA) emerges as an essential, semi-quantitative risk assessment method for evaluating fishery stocks when traditional, data-intensive stock assessments are impossible [1].

This article provides application notes and detailed protocols for implementing PSA and complementary strategies in environments of extreme data poverty. It is framed within a broader thesis positing that effective resource management is a fundamental pillar for alleviating poverty in dependent communities. By providing structured methodologies for decision-making in the absence of robust data, these protocols aim to bridge the critical gap between information scarcity and the urgent need for sustainable, productivity-enhancing fisheries management.

Quantitative Landscape: Poverty and Data Gaps

Effective strategy formulation requires an understanding of the scale of both economic and data poverty. The following tables summarize the current quantitative landscape.

Table 1: Global Extreme Poverty Trends and Projections (1990-2040) [46] [47]

Year	People in Extreme Poverty (Millions)	Global Population (%)	Key Notes
1990	2,310	N/A	Baseline using previous poverty line.
2022	~1,500 (under $3/day line)	N/A	Based on revised $3/day 2021 PPP line [47].
2025	808 - 831 [46] [47]	9.9% [47]	Projection; over 75% live in Sub-Saharan Africa or fragile states.
2030	793 (projected)	N/A	Projected decline based on current trends [46].
Post-2030	Projected increase [46]	N/A	Expected due to stagnation in poorest economies.
Primary Driver	Economic growth (or lack thereof) in poorest countries [46].
SDG 1.1 Outlook	Eradication by 2030 is "highly unlikely" [47].

Table 2: Characterization of Fishery Data Poverty Tiers

Data Poverty Tier	Typical Data Available	Management Context	Recommended Assessment Approach
Tier 1: Data-Limited	Catch time-series (potentially unreliable), limited life history parameters.	Small-scale, artisanal, or developing fisheries.	Core PSA [1], Catch-only models, traditional ecological knowledge (TEK) integration.
Tier 2: Data-Poor	Partial catch data, sporadic effort data, rough biological parameters from related species.	Many tropical, multi-species, or remote fisheries.	Enhanced PSA with uncertainty scoring, hybrid models using analogous data.
Tier 3: Extreme Data Poverty ("No Data")	No reliable time-series data. Anecdotal catch, undefined fishing grounds, unknown species composition.	Pre-development fisheries, unregulated informal sectors, post-conflict regions.	Rapid PSA based on TEK and expert elicitation, participatory resource mapping, pilot data collection protocols.

Core Protocol: Productivity and Susceptibility Analysis (PSA) in Data-Poor Contexts

PSA is a semi-quantitative risk assessment that evaluates a stock's vulnerability by scoring its inherent productivity (biological capacity to recover) and its susceptibility to the fishery [1]. The following protocol is adapted for contexts where data is absent or highly uncertain.

Phase 1: Preparatory Framework and Attribute Selection

• Objective: Define the fishery system and select/scorable attributes with available information.
• Procedure:
  • System Definition: Clearly define the stock (species complex or single species) and the fishery (gear types, sectors).
  • Attribute Customization: Use the standard PSA attributes (10 productivity, 12 susceptibility) [1]. For each attribute, define what constitutes Low (1), Medium (2), and High (3) risk for your specific context. Example: For "Maximum Age," define thresholds based on regional knowledge (e.g., >15 years = Low risk, <5 years = High risk).
  • Data Quality Matrix: Create a companion score (e.g., High/Medium/Low) for the quality of information underlying each attribute score. This explicitly documents uncertainty.

• Objective: Populate the PSA matrix where conventional data is absent.
• Procedure:
  • Expert and Local Knowledge Elicitation: Conduct structured interviews with fishers, traders, and processors. Use ordinal scales (e.g., low/medium/high) to translate qualitative knowledge into scores for attributes like "Spatial Distribution," "Aggregation," and "Seasonal Migration."
  • Use of Analogous Data: For biological productivity attributes (e.g., "Natural Mortality," "Fecundity"), use values from closely related species or populations in similar ecosystems. Document the source and justify the analogy.
  • Participatory Mapping: For susceptibility attributes related to "Overlap" and "Habitat," use participatory GIS (PGIS) methods to map fishing grounds and habitats based on community knowledge.
  • Default Scoring and Flagging: If no information can be sourced for an attribute, assign a precautionary default score (typically higher risk) and flag it for priority investigation.

Phase 3: Calculation, Visualization, and Risk Prioritization

• Objective: Synthesize scores to identify highest-risk stocks and inform data collection.
• Procedure:
  • Calculate the average score for all Productivity (P) and Susceptibility (S) attributes.
  • Plot results on a bivariate risk plot (P vs. S). The overall vulnerability score can be calculated as the Euclidean distance from the origin: Risk = √(P² + S²).
  • Prioritization: Stocks in the high-susceptibility, low-productivity quadrant are highest priority for management intervention and targeted data collection.
  • Sensitivity Analysis: Re-run scoring using alternative, reasonable assumptions to test the robustness of the risk ranking.

Supplemental Protocol: Structured Data Reconstruction and Standardization

When a minimal foothold of data exists, a systematic approach to its organization and enhancement is critical. This protocol is based on standardized fishery data frameworks [48].

Implementing a Standardized Two-File Data System

• Objective: Create a future-proof, clear, and machine-readable data record.
• Procedure:
  • Create Data Input File: Populate a standardized spreadsheet (CSV or Excel) with all available quantitative data [48]. The structure includes:
    • Metadata: Species name, region, last historical year of data.
    • Biological Parameters: Estimates for natural mortality (M), growth (Linf, K), maturity, even if from analogues.
    • Time-Series: Any catch, effort, or relative abundance index data. Gaps are left blank.
    • Composition Data: Any length or age frequency samples.
  • Create Data Documentation File: A corresponding text (.md) file that provides narrative for every data entry in the Input File [48]. For each parameter or time-series, document:
    • Source: "Fisher interviews, 2023," "Analogous species (X) from region Y."
    • Method of Derivation: "Averaged from three elder fishers," "Scaled from species X using a length-weight ratio."
    • Uncertainty Statement: "Considered highly uncertain due to seasonal variability."

Gap-Filling via Iterative Stock Reduction Analysis (iSRA)

• Objective: Reconstruct a plausible historical catch time-series from fragmented data.
• Procedure:
  • Anchor with Known Points: Use any reliable recent catch estimate or effort index as an anchor.
  • Define Priors: Use life history parameters (from PSA) to define prior distributions for intrinsic population growth rate (r) and carrying capacity (K).
  • Run iSRA Model: Use a tool (e.g., iSRA in R DLMtool) to estimate a trajectory of historical catches and biomass that is consistent with the anchor point and life history priors.
  • Sensitivity Testing: Run multiple scenarios with different anchors and priors to produce an envelope of plausible histories. This output directly informs the uncertainty scored in the PSA.

The Scientist's Toolkit: Essential Reagents for Data-Poor Research

Table 3: Research Reagent Solutions for Data-Poor Fisheries Science

Item / Solution	Function in Addressing Data Poverty	Key Considerations & Protocols
Standardized Data Template [48]	Provides a consistent framework for organizing all available quantitative fishery data (catch, biology, effort) and linking it directly to narrative documentation.	Use the DLMtool/MSEtool DataInit() function to generate templates. Ensures data is machine-readable for stock assessment and Management Strategy Evaluation (MSE).
PSA Attribute Scoring Rubric	Translates qualitative, anecdotal, or analogue information into semi-quantitative scores for risk assessment. Enables comparison across species/fisheries.	Must be customized locally for each fishery system. Requires participatory development with stakeholders to define risk thresholds.
Expert Elicitation Protocol	Structured method to formally gather and quantify the knowledge of fishers, processors, and other local experts.	Uses calibrated questionnaires, ordinal scoring sheets, and facilitated workshops. Essential for scoring susceptibility attributes like aggregation and behavior.
Analogous Life History Database	Curated database of biological parameters (M, Linf, K, maturity) for related species/regions. Used to fill critical gaps when target species data is absent.	Must document sources and include uncertainty metrics. Should be regionally specific to improve relevance.
Participatory Mapping Kit	Enables the spatial documentation of fishing grounds, habitat types, and seasonal patterns using local knowledge.	Includes printed base maps, GPS devices, and markers. Outputs shapefiles for GIS analysis of spatial overlap (a key PSA susceptibility attribute).
Rapid Catch Monitoring App	Simple digital tool for recording catch at landing sites. Establishes a baseline time-series where none exists.	Must be low-cost, offline-capable, and designed for easy use by enumerators or fishers. Data feeds directly into the Standardized Data Template.

Visualization and Communication Guidelines

Communicating findings from data-poor contexts requires exceptional clarity to build trust and inform action [49] [50].

• For PSA Results: Use a bivariate scatter plot (Productivity vs. Susceptibility) with points colored/shaped by overall risk score. This immediately highlights priority stocks [50].
• For Data Quality: Incorporate data quality scores directly into visualizations (e.g., using point transparency or halo size) to visually convey uncertainty [49].
• Color Palette Application:
  • Use a sequential palette (e.g., light to dark blue) for ordered data like risk scores [49].
  • Use a diverging palette (e.g., red-white-green) cautiously, only for data with a meaningful midpoint (e.g., deviation from a target) [49].
  • Critical - Contrast Rule: Ensure a minimum contrast ratio of 4.5:1 for standard text and 3:1 for large text or graphical elements against their background [51] [52]. Use specified palette colors like #202124 (dark gray) on #F1F3F4 (light gray) or #FFFFFF (white) on #4285F4 (blue).
• Narrative Integration: Combine charts with brief, direct text annotations that tell the story of the data—highlighting key gaps, assumptions, and management implications [50].

Productivity and Susceptibility Analysis (PSA) serves as a semi-quantitative risk assessment tool within Ecosystem-Based Fisheries Management (EBFM). Designed for data-limited scenarios, it rapidly evaluates the relative vulnerability of target and non-target species to fishing pressures by scoring biological (productivity) and fishery-interaction (susceptibility) attributes [5]. The core output is a vulnerability score, typically calculated as the Euclidean distance from the origin on a plot of averaged productivity (P) and susceptibility (S) scores, which are then categorized into Low, Medium, or High risk [5]. Its primary utility lies in prioritizing management actions and directing research resources toward species at highest perceived risk [53].

However, the PSA methodology is inherently subjective. Critical challenges include the selection and weighting of attributes, the definition of scoring thresholds for each attribute, and the aggregation of scores into final risk categories. Studies have demonstrated that subjective weightings can significantly alter model outcomes and increase uncertainty in risk classification [53]. Furthermore, a foundational quantitative evaluation has questioned the underlying assumptions of PSA, indicating that its predictive performance can be poor and that the information required is often comparable to that needed for more robust, quantitative stock assessment models [5]. This document provides application notes and detailed protocols for optimizing PSA attribute scoring to manage these sources of subjectivity and establish robust, defensible thresholds, thereby improving the reliability of this rapid assessment tool within a comprehensive ecological risk assessment strategy [54].

Protocol for Optimized Attribute Selection and Redundancy Analysis

A critical step in reducing subjectivity is to refine the attribute list to minimize redundancy without losing biological or fishery significance.

Objective

To identify and eliminate redundant biological and susceptibility attributes from a standard PSA matrix to create a simplified, non-redundant attribute set that maintains the analytical power of the assessment [53].

Materials & Required Data

• Initial Attribute Matrix: A comprehensive list of potential productivity (e.g., maximum age, growth rate, fecundity, age at maturity) and susceptibility (e.g., spatial overlap, post-capture mortality, gear selectivity) attributes.
• Species Data Set: Life history and fishery interaction parameter estimates for a broad spectrum of species (e.g., 50-100 species) representative of the fishery system.
• Statistical Software: Capable of performing multivariate analysis (e.g., R, PRIMER, SPSS).

Step-by-Step Methodology

• Data Compilation and Scoring: For each species in the training dataset, score all initial attributes using predefined, traditional thresholds.
• Correlation Analysis: Calculate a pairwise correlation matrix (e.g., Pearson or Spearman correlation) for all attributes.
• Redundancy Identification: Flag attribute pairs with correlation coefficients (|r|) > 0.7. This high correlation suggests the attributes convey similar information.
• Biological Justification Review: For each flagged pair, convene an expert panel to determine if the correlation is biologically meaningful or coincidental. Retain the attribute deemed more reliable to estimate or more directly relevant to the core concept (productivity or susceptibility).
• Multivariate Validation: Perform a Principal Component Analysis (PCA) on the attribute scores. Attributes that cluster tightly together in the PCA biplot confirm redundancy. The goal is to achieve a set of attributes that load on different principal components, explaining unique variance.
• Final Set Definition: Establish the final, reduced attribute set. Document the justification for removing each redundant attribute.

Expected Outcome & Application

The IATTC's application of this protocol for tuna fisheries demonstrated that the number of biological attributes could be reduced effectively [53]. This optimization decreases the data burden, simplifies the assessment, and mitigates the risk of double-counting certain life-history traits, leading to a more robust and streamlined model.

Table 1: Example of Attribute Redundancy Analysis for Productivity Traits

Attribute Pair	Correlation (r)	Recommended Action	Rationale
Maximum Age vs. Natural Mortality (M)	-0.92	Retain Natural Mortality (M); remove Maximum Age.	M is a more direct and widely used input for population models. The inverse relationship is well-established.
Fecundity vs. Size at Maturity	0.65	Retain both.	Correlation is moderate. Expert review confirms they represent distinct reproductive strategy elements.
von Bertalanffy K vs. Age at Maturity	-0.81	Retain Age at Maturity; remove growth coefficient K.	Age at maturity is a more stable and directly measurable parameter for data-limited species.

Protocol for Establishing Standardized Scoring Thresholds

Subjectivity in scoring arises from arbitrary breakpoints between risk categories (e.g., Low, Medium, High). This protocol establishes data-driven thresholds.

Objective

To derive objective, data-informed numerical thresholds for scoring continuous life-history attributes (e.g., fecundity, age at maturity) into discrete risk categories.

Materials & Required Data

• Global Life History Database: e.g., FishBase, PANTHERA.
• Regional Species List: The list of species subject to PSA.
• Statistical Software: For percentile analysis and distribution fitting.

Step-by-Step Methodology

• Data Assembly: Extract parameters for the relevant attributes from the global database for all species within the taxonomic class/order relevant to your region (e.g., all Actinopterygii).
• Distribution Analysis: Plot the frequency distribution (histogram) for each attribute (e.g., log-transformed fecundity). Analyze the distribution to identify natural breaks, modes, or percentiles.
• Percentile-Based Thresholding: Define scoring thresholds based on the distribution of the global dataset, not the assessed species list. A proposed framework is:
  • Low Risk (Score=1): Attribute value > 66th percentile of global data (high productivity/low susceptibility).
  • Medium Risk (Score=2): Attribute value between 33rd and 66th percentile.
  • High Risk (Score=3): Attribute value < 33rd percentile of global data (low productivity/high susceptibility).
• Regional Calibration (Optional): Compare the distribution of your regional species list against the global percentiles. If a significant bias exists, document it, but prefer global percentiles to maintain consistency and avoid introducing local bias.
• Threshold Documentation: Explicitly publish the final numerical thresholds for every scored attribute. For example: "Fecundity: Low Risk (>10,000 eggs), Medium Risk (100–10,000 eggs), High Risk (<100 eggs)" [5].

Expected Outcome & Application

This replaces expert-opinion-based breakpoints with a reproducible, transparent standard. It acknowledges that a "low productivity" species in one region is biologically "low productivity" everywhere, allowing for comparable risk outcomes across studies and regions [5].

Table 2: Data-Driven Scoring Thresholds Based on Global Percentiles

Attribute	High Risk (3) - Low Productivity	Medium Risk (2)	Low Risk (1) - High Productivity	Basis
Age at Maturity (years)	> 10	3 - 10	< 3	Global percentiles for marine fishes [5].
Fecundity (# of eggs)	< 100	100 - 10,000	> 10,000	Common literature values, aligned with percentile logic [5].
Maximum Size (cm)	> 200	50 - 200	< 50	Adapted from size-based productivity classifications.

Workflow for an Optimized PSA Implementation

The following diagram illustrates the integrated workflow for conducting an optimized PSA, incorporating the protocols for attribute selection and threshold standardization.

Diagram: Workflow for an Optimized PSA Implementation

The Researcher's Toolkit for PSA Studies

Implementing a robust PSA requires specific analytical tools and conceptual frameworks. The following toolkit is essential for researchers.

Table 3: Essential Research Toolkit for PSA Implementation

Tool / Resource	Category	Primary Function in PSA	Key Consideration
Global Life History Databases (FishBase, SeaLifeBase)	Data Source	Provide parameter estimates (e.g., growth, maturity, fecundity) for scoring attributes, especially for data-limited species.	Quality is variable; use as priors or for comparative analysis, not absolute truths.
Multivariate Statistical Packages (R, PRIMER, SPSS)	Analysis Software	Execute correlation analysis, PCA, and other methods for attribute redundancy analysis and validation.	Essential for implementing the attribute optimization protocol.
Hierarchical Risk Assessment Framework (e.g., ERAEF) [5]	Conceptual Model	Positions PSA as a Tier 2 screening tool to prioritize species for deeper, Tier 3 quantitative assessment.	Clarifies PSA's role as a screening tool, not a definitive stock assessment.
Reference Libraries on Data-Limited Methods [54]	Guidance	Provide context on alternative and complementary assessment techniques (e.g., length-based, catch-only models).	Allows researchers to situate PSA within the broader toolbox and understand its limitations.
Expert Elicitation Protocols	Methodology	Structured process to gather informed judgments on attribute weighting or scoring when data are utterly absent.	Must be structured and documented to minimize individual bias. Use sparingly.

Optimizing attribute scoring through redundancy removal and threshold standardization significantly enhances the transparency, reproducibility, and objectivity of PSA. These protocols directly address the core criticisms of subjectivity and arbitrary assumptions [53] [5].

It is imperative to recognize that PSA provides a relative risk ranking, not an absolute measure of stock status or sustainability. The IATTC's experience concluded that PSA can require substantial data while offering only relative risk outputs, leading to a shift towards more quantitative approaches like EASI-Fish for ecological risk assessment [53]. Therefore, an optimized PSA is most powerful as a prioritization and screening tool within a tiered risk assessment strategy [5] [54]. Species flagged as "high risk" should be candidates for more data collection and rigorous, quantitative assessment techniques to inform definitive management actions.

Within the framework of Productivity and Susceptibility Analysis (PSA) for fisheries research, a critical challenge is the assessment of data-poor stocks and the incorporation of spatially explicit, local-scale information into management. Traditional PSA relies on biological productivity attributes and fishery susceptibility characteristics to estimate a stock's vulnerability to overfishing [9]. However, coastwide assessments often fail to capture fine-scale ecological and fishery dynamics, leading to potential mismatches between broad-scale management and local stock status [9]. This application note posits that the integration of fishers' expertise with conventional scientific data creates a more robust, nuanced, and actionable PSA. This integrated approach enhances the accuracy of vulnerability scores, informs management at appropriate spatial scales, and builds essential trust and collaboration between scientific and fishing communities [55]. The protocols detailed herein provide a methodological pathway for this integration, directly contributing to the broader thesis that participatory, multi-knowledge frameworks increase the effectiveness and legitimacy of fisheries assessments.

Quantitative Data Synthesis from Integrated Assessments

The value of integrating knowledge sources is quantitatively demonstrated by differences in vulnerability outcomes across spatial scales and knowledge types. The following tables summarize key findings from recent studies.

Table 1: Comparative PSA Vulnerability Scores: Coastwide vs. Local Integrated Assessment [9]

Species (Example)	Coastwide PSA Vulnerability Score	Local Integrated PSA Vulnerability Score	Knowledge Inputs for Local Assessment
Nearshore Rockfish A	High	Medium-High	Scientific life-history data + Fisher interviews on local abundance & fishing pressure
Nearshore Rockfish B	Medium	Medium	Published growth parameters + Survey data on local catch composition
General Trend	Higher vulnerability	Moderated vulnerability	Local fishers' knowledge often reflected healthier local stock conditions or different susceptibility factors, lowering final scores.

Table 2: PSA Vulnerability of Data-Poor Peruvian Coastal Groundfish (Expert-Driven Assessment) [3]

Species	Vulnerability Score	Vulnerability Category	Primary Knowledge Source
Broomtail Grouper	2.57	Extremely High	Expert biological & fishery knowledge
Grape-eye Seabass	2.50	Extremely High	Expert biological & fishery knowledge
Pacific Goliath Grouper	2.28	Extremely High	Expert biological & fishery knowledge
Black Snook	2.13	High	Expert biological & fishery knowledge
Bumphead Parrotfish	1.73	Low	Expert biological & fishery knowledge
Overall Data Quality	"Limited" to "None"	-	Highlights critical reliance on expert knowledge in data-poor contexts.

Experimental Protocols for Knowledge Integration and Biological Analysis

This protocol outlines a systematic process for gathering and quantifying fishers' expertise for integration into Susceptibility attributes of a PSA.

Objective: To collect spatially resolved, fishery-specific knowledge on species susceptibility (e.g., encounterability, selectivity, post-capture mortality) and qualitative data on productivity (e.g., observed growth rates, habitat preferences) to complement scientific data.

Materials:

• Pre-defined PSA scoring rubric (1-3 scale for low-high susceptibility/productivity).
• Survey instruments (online and physical) with visual aids (species photos, maps).
• Recording equipment (with consent) for interviews.
• GIS software for spatial data mapping.

Procedure:

• Participant Identification & Recruitment: Partner with fishing cooperatives, harbormasters, and community leaders to identify knowledgeable fishers. Use stratified sampling to ensure representation across gear types, vessel sizes, and fishing regions [9].
• Structured Workshops & Surveys: Conduct facilitated group meetings and one-on-one interviews.
  • Present the PSA framework and scoring system.
  • For target species, guide participants through scoring susceptibility attributes (e.g., "On a scale of 1-3, how likely are you to encounter this species in your fishing grounds?").
  • Use participatory mapping to document fishing grounds, species hotspots, and habitat changes [55].
• Data Processing & Quantification:
  • Transcribe and code qualitative responses.
  • Calculate median or modal scores from fishers' ratings for each PSA attribute.
  • Georeference spatial data to create layers of fishing effort and species distribution.
• Integration with Scientific Data:
  • Convergence: Where scientific data (e.g., survey CPUE) and fisher scores align, use a weighted average to produce a final score.
  • Divergence: Where they differ, document the discrepancy. Investigate through follow-up or scientific analysis. This divergence often reveals locally important, science-blind processes [55].
• Validation & Feedback: Present integrated results back to participants for review and validation, closing the feedback loop and ensuring accurate interpretation.

Protocol B: Machine Vision Analysis of Fish Behavior for Stressor Response

This protocol details a quantitative, non-contact method to assess fish behavior as a biomarker of sub-lethal stress (e.g., from climate change or pollution), which can inform "Productivity" metrics related to resilience [56].

Objective: To quantify changes in locomotion, spatial distribution, and anxiety behaviors in fish exposed to environmental stressors (e.g., ammonia, altered pH) using a lightweight deep learning model.

Materials:

• Experimental aquaria/tanks with controlled environmental systems.
• High-resolution video cameras.
• Computing workstation with GPU.
• Software: Python, PyTorch/TensorFlow, DeepMOT tracker or similar.
• LAMP-SAG-YOLOv8 model architecture (pruned YOLOv8 for efficient fish detection) [56].

Procedure:

• System Setup & Acclimation:
  • Set up tanks with control (e.g., pH 8.1) and treatment conditions (e.g., pH 7.8, simulating ocean acidification) [57].
  • Acclimatize test fish (e.g., black sea bream, bass) for two weeks.
  • Position cameras orthogonally to record the entire tank volume.
• Video Acquisition & Model Inference:
  • Record video segments (e.g., 10-minute intervals) across multiple days.
  • Process video frames using the pre-trained LAMP-SAG-YOLOv8 model to detect and locate individual fish.
  • Feed detection outputs into the DeepMOT tracker to maintain individual identities across frames, generating trajectory data [56].
• Behavioral Metric Extraction:
  • Swimming Distance & Velocity: Calculate from trajectory coordinates.
  • Spatial Distribution: Analyze heatmaps of position frequency (e.g., time spent near surface vs. bottom).
  • Anxiety Behavior: Quantify erratic movement patterns (e.g., rapid turning, burst swimming).
• Statistical Analysis:
  • Compare behavioral metrics between control and treatment groups using multivariate statistics (e.g., PERMANOVA).
  • Correlate behavioral shifts with physiological data (e.g., growth, feed efficiency) [57].

Protocol C: Geometric Morphometrics for Population Structure Analysis

This protocol provides a workflow for analyzing shape variation in fish, which can indicate population sub-structure, local adaptations, or responses to environmental gradients—critical information for defining assessment units in PSA [58].

Objective: To quantify and compare body shape morphology among fish samples from different geographic regions using landmark-based geometric morphometrics.

Materials:

• Standardized digital photographs of fish specimens (lateral view).
• Software: tpsDig2, MorphoJ, R with geomorph and Momocs packages [58].
• Calibration scale.

Procedure:

• Image Preparation:
  • Photograph fish on a neutral background with a fixed scale and orientation (lateral view, left side preferred).
  • Use ImageJ or an AI background remover to standardize images.
• Landmarking:
  • Digitize Type I (anatomical) and Type II (geometric) landmarks on each image using tpsDig2 (e.g., tip of snout, insertion of fins, points of maximum curvature) [58].
  • For outline data, place semi-landmarks along the body contour.
• Shape Analysis:
  • Perform Generalized Procrustes Analysis (GPA) in MorphoJ to align landmarks, removing effects of size, position, and rotation.
  • Conduct Principal Component Analysis (PCA) on Procrustes coordinates to visualize major axes of shape variation.
  • Perform Canonical Variate Analysis (CVA) or Discriminant Function Analysis (DFA) to test for significant shape differences between pre-defined groups (e.g., fishing regions).
• Visualization & Interpretation:
  • Generate deformation grids (thin-plate splines) to visualize shape changes associated with PCA or CVA axes.
  • Interpret shape differences in an ecological context (e.g., body depth related to habitat type).

Visualizations of Workflows and Relationships

Integrated Knowledge PSA Workflow (71 characters)

AI Fish Behavior Analysis Pipeline (37 characters)

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Materials for Integrated Fisheries Assessment Research

Item/Category	Function/Application	Example/Note
Participatory Elicitation
Structured Survey Platforms	Quantifies fisher and expert knowledge into scalable scores for PSA attributes.	Online forms (e.g., Qualtrics), physical scorecards with visual aids [9] [55].
Participatory GIS (PGIS) Tools	Geospatially references local ecological knowledge, mapping fishing grounds, hotspots, and habitat changes.	Open-source QGIS with plugin for community mapping; printed nautical charts for manual annotation [55].
Biological Analysis
Environmental Control System	Maintains precise water chemistry (pH, temperature, ammonia) for stressor exposure experiments.	Computer-controlled CO2 bubbling systems for acidification studies [57].
Machine Vision Setup	Enables non-contact, high-throughput behavioral phenotyping for stress response.	RGB or IR cameras, GPU workstation, LAMP-SAG-YOLOv8 software model [56].
Geometric Morphometrics Software	Quantifies shape variation for population structure analysis.	tpsDig2 (landmarking), MorphoJ & R geomorph (statistical shape analysis) [58].
Sample Analysis
Histology & Microscopy Reagents	Assesses sub-lethal physiological impacts of stressors (e.g., ocean acidification) on fish health.	Fixatives (formalin), stains (H&E), SEM/TEM preparation kits for gill/liver tissue [57].

Weighting Attributes and Conducting Sensitivity Analysis to Identify Key Risk Drivers

Productivity and Susceptibility Analysis (PSA) serves as a cornerstone risk assessment framework in fisheries science, designed to evaluate the relative vulnerability of species or stocks to overexploitation. The core objective of PSA is to systematically integrate information on a stock's biological productivity—its capacity to recover from depletion—with its exposure to anthropogenic pressures, primarily fishing. The analytical rigor and management utility of a PSA are fundamentally dependent on two interrelated processes: the statistically sound weighting of input attributes and a thorough sensitivity analysis to identify which factors most strongly govern the output risk scores. Mismanagement of these processes can lead to biased assessments, misdirected conservation efforts, and ultimately, unsustainable fisheries. This document provides detailed application notes and protocols for these critical steps, framed within a broader thesis on enhancing the objectivity and reliability of PSA outcomes for research and management.

Theoretical Foundations and Error Categorization

A robust PSA must account for various sources of uncertainty and error. Thorson et al. (2023) provide a critical framework by categorizing errors in integrated assessments into four distinct types [59]:

• Sampling Bias: Systematic error in data collection (e.g., non-representative survey coverage).
• Sampling Imprecision: Random error due to finite sample sizes.
• Assessment Model Bias: Error from incorrect structural model assumptions (e.g., wrong natural mortality rate).
• Assessment Model Imprecision: Error from estimating too many parameters, often addressed with random effects.

A key diagnostic derived from this framework is the Percent Excess Variance (PEV). The PEV quantifies the net effect of unmodeled processes (bias) by comparing the variance expected from sampling alone (Input Sample Size, ESS) to the variance observed in model residuals (Effective Sample Size, ESS) [59]: PEV = 100% × (1 - Input SS / Effective SS)

A high PEV indicates that data conflicts with the model, signaling potential misspecification that should be investigated rather than simply resolved by down-weighting data [59].

Table 1: Categorization and Diagnostics for Error Sources in Stock Assessment (adapted from [59])

Error Category	Definition	Common Source in PSA	Proposed Diagnostic
Sampling Bias	Systematic deviation between observed and true value.	Non-representative spatial coverage of surveys; gear selectivity.	Analysis of residual patterns across strata.
Sampling Imprecision	Random error due to finite sampling effort.	Low sample sizes for life-history parameters (e.g., maturity, fecundity).	Calculation of coefficient of variation (CV) for input data.
Assessment Model Bias	Error due to incorrect model structure or assumptions.	Use of inappropriate productivity or susceptibility indices.	Percent Excess Variance (PEV) > 20-30% [59].
Assessment Model Imprecision	Error from estimating poorly informed parameters.	Including highly correlated or uncertain attributes without regularization.	Wide confidence intervals for key output parameters.

Application Note 1: Protocols for Weighting Attributes in PSA

Data Pre-processing and Input Sample Size Estimation

Objective: To assign appropriate statistical leverage to each data type (e.g., abundance index, age-composition data) within the PSA model. Background: Composition data (e.g., proportion-at-age) are often overweighted in stock assessments because their true variance, which includes both sampling and process error, is underestimated [59].

Protocol: Estimating Input Sample Size for Composition Data

• Data Expansion: Using design-based estimators, expand survey catch data to generate abundance indices and expanded numbers-at-age (n_(a,y)) [59].
• Variance Estimation: For composition data derived from surveys, calculate the Input Sample Size (Input SS). This represents the amount of information, adjusted for survey design, and serves as a prior for the likely sampling variance [59]. For length- or age-stratified sampling, this can be derived from the design's effective sample size.
• Likelihood Specification: Fit the model using a likelihood (e.g., Dirichlet-multinomial) that incorporates the Input SS as a data-weighting parameter. The model will estimate an Effective Sample Size (Effective SS) from the residuals.
• Diagnostic Calculation: Compute the Percent Excess Variance (PEV) as shown in Section 2. A high PEV for a fleet or data component suggests a conflict requiring model investigation, not just re-weighting [59].
• Iterative Refinement: If PEV is high, investigate sources of model bias (e.g., time-varying selectivity, misspecified natural mortality) and refine the model structure before finalizing weights.

Objective: To formalize the integration of expert judgment for attributes lacking robust quantitative data. Background: In data-limited PSA contexts, many attributes (e.g., habitat specialization, management effectiveness) rely on semi-quantitative scores. Weighting these attributes subjectively is common but must be structured to minimize bias.

Protocol: Structured Expert Elicitation for Attribute Weights

• Selection: Convene a panel of 5-10 independent experts covering biology, fisheries, and ecology.
• Calibration: Train experts on the PSA framework and provide them with a consistent set of background data and definitions for each attribute.
• Elicitation: Using a modified Delphi process:
  • Round 1: Experts privately assign initial weight (e.g., 1-10) and uncertainty range to each attribute.
  • Round 2: Provide anonymized summaries (median, range) of Round 1 responses. Experts discuss rationales and revise their weights.
  • Round 3: Conduct a final, private weighting. The median or mean of final weights is adopted.
• Documentation: Record all weights, rationales, and measures of expert agreement (e.g., inter-quartile range). This audit trail is critical for transparency and review [60].

Application Note 2: Protocols for Sensitivity and Global Sensitivity Analysis (GSA)

Local Sensitivity Analysis (One-at-a-Time)

Objective: To understand the localized effect of changing a single input attribute on the PSA risk score. Protocol:

• Define a base case using the best available point estimates for all n attributes.
• For each attribute i, define a plausible range of variation (e.g., ±25% of its base value, or its 95% confidence interval).
• For attribute i, set its value to the lower and upper bounds of its range, while holding all other n-1 attributes at their base value.
• Re-run the PSA model and record the change in the output risk score. Calculate the normalized sensitivity index, S_i: S_i = (ΔOutput / Output_base) / (ΔInput_i / Input_i_base)
• Rank attributes by the absolute magnitude of S_i. This identifies parameters where precise estimation is most critical.

Limitation: This method does not account for interactions between attributes and may miss effects visible only when multiple factors change simultaneously.

Global Sensitivity Analysis (GSA) Using Novel Probability-Weighted Moments (PWM)

Objective: To apportion output uncertainty to the uncertainty in all input attributes, including interactions, across their entire joint distribution. Background: Variance-based Sobol' indices are common but can be insufficient if output distributions are non-Gaussian. A novel, robust GSA measure uses Probability Weighted Moments (PWMs), which are less sensitive to outliers and more reliable with small samples than classical moments [61] [62].

Protocol: PWM-Based Global Sensitivity Analysis [61]

• Define Input Distributions: For each of the n PSA input attributes (X_i), specify a probability distribution (e.g., uniform, triangular, normal) representing knowledge uncertainty.
• Generate Sample Matrix: Create an N × 2n sample matrix using a Quasi-Monte Carlo method (e.g., Sobol' sequence). The matrix is composed of two independent sampling matrices, A and B.
• Compute PWM of Output: Run the PSA model for all N samples in matrix A to generate output vector Y_A. Estimate the k-th order PWM (β_Y^k) for the unconditional output distribution [61].
• Compute Conditional PWM: For each input X_i, create a hybrid matrix where X_i is taken from matrix B and all other variables are from matrix A. Run the model to get output Y_(B~i). Estimate the conditional PWM, β_(X~i)^k (Y|X_i), averaged over all values of X_i [61].
• Calculate Sensitivity Index: The main effect PWM-based sensitivity index for attribute X_i and moment k is: ω_i^k = 1 - [E_(X_i)(β_(X~i)^k (Y|X_i))] / [β_Y^k (Y)] This measures the fractional reduction in the k-th output PWM variance expected from learning the true value of X_i [61].
• Interpretation: High ω_i^k values indicate attributes that are key drivers of uncertainty in the output's distribution, as characterized by its PWMs. This method is particularly effective for identifying drivers of extreme (tail) behavior in risk scores.

Table 2: Comparison of Sensitivity Analysis Methods for PSA

Method	Key Principle	Strengths	Weaknesses	Best Use Case in PSA
Local (One-at-a-Time)	Vary one input at a time around a baseline.	Simple, intuitive, computationally cheap.	Misses interactions; depends on baseline choice.	Screening analysis; prioritizing data collection for specific stocks.
Variance-Based (Sobol')	Decompose output variance into contributions from inputs.	Model-free, accounts for interactions, provides total effect indices.	Computationally expensive; assumes variance is sufficient output descriptor.	Comprehensive analysis for well-sampled stocks with moderate uncertainty.
PWM-Based [61]	Decompose uncertainty in Probability Weighted Moments of output.	Robust to outliers and small samples; informative for tail risks.	Novel, less established in ecological literature; requires specialized code.	Analysis for data-limited stocks or when assessing risk of extreme outcomes.
Interval & Reliability Analysis [63]	Treat inputs as intervals (p-boxes) to find bounds on failure probability.	Handles deep epistemic (knowledge-related) uncertainty.	Computationally very intensive; results are intervals, not point estimates.	When inputs are highly uncertain (e.g., new fishery, climate change impacts).

Diagram 1: Workflow for PWM-Based Global Sensitivity Analysis (GSA) [61]

Integrated PSA Workflow for Identifying Key Risk Drivers

This protocol integrates the weighting and sensitivity analysis steps into a coherent PSA workflow, using a Bayesian state-space model as the integrating engine, similar to frameworks used in data-limited stock assessments like CMSY/BSM [64] and WHAM [59].

Protocol: Integrated Risk Driver Analysis

• Problem Formulation & Data Assembly:
  • Define the stock(s) and fishery system.
  • Assemble data: catch time series (essential), any abundance index, and life-history parameters (e.g., L_∞, K, M) [64] [65].
  • For qualitative attributes, conduct structured expert elicitation (Sec. 3.2).

• Base Model Construction:

  • Implement a Bayesian state-space model (e.g., a Schaefer production model) [64].
  • Weighting: Use formal methods to assign likelihood weights. For composition data, use the Dirichlet-multinomial with Input SS [59]. For indices, use assumed CVs.
  • Prior Specification: Assign informative priors based on life-history invariants or meta-analysis for parameters like r (intrinsic growth rate) and K (carrying capacity) [64].

• Model Fitting & Diagnostic Checking:

  • Fit the model using MCMC sampling.
  • Check convergence (Gelman-Rubin statistic).
  • Calculate PEV diagnostics for all data components. Investigate and rectify sources of high PEV (e.g., add time-varying parameters) [59].

• Sensitivity Analysis:

  • Using the fitted model as the base, perform Local Sensitivity Analysis on key posterior estimates (e.g., B/B_MSY, F/F_MSY) to input data and priors.
  • Perform a Global Sensitivity Analysis (preferably using the PWM-based method for robustness) [61] treating major uncertain inputs (e.g., M, priors on r) as random variables with defined distributions.

• Synthesis and Reporting of Key Risk Drivers:

  • Primary Drivers: Attributes with high global sensitivity indices (ω_i^k) and high inherent uncertainty. These are the core risk drivers (e.g., natural mortality M).
  • Secondary Drivers: Attributes with moderate sensitivity or high sensitivity but low uncertainty.
  • Management Levers: Attributes with high sensitivity that are directly controllable (e.g., fishing mortality F). The analysis quantifies the risk reduction achievable by managing this lever.

Diagram 2: Integrated PSA Workflow for Risk Driver Identification

Case Study Application: Seabream Stock in Oman

A study on the seabream fisheries of Oman applied the CMSY and Bayesian Schaefer Model (BSM) methods [64], which align with the integrated workflow described.

• Context: Data-limited stock, with primary catch time series from 1988-2021.
• Method/Weighting: The BSM is a Bayesian state-space model that uses priors for r and K. Likelihood weighting is applied to the catch and abundance index data.
• Sensitivity Analysis Implied: The Bayesian framework inherently produces posterior distributions for stock status benchmarks (e.g., B/B_MSY, F/F_MSY). The width of these posteriors reflects combined input uncertainty.
• Key Risk Driver Identified: The results indicated a 53% probability that the stock was overfished and undergoing overfishing (B/B_MSY = 0.96, F/F_MSY = 1.25) [64]. Fishing pressure (F) was identified as the primary, controllable risk driver.
• Management Implication: The study concluded that reducing fishing mortality is the essential management action to lower risk, directly informed by the quantitative analysis [64].

Table 3: Key Research Reagent Solutions for PSA

Tool Category	Specific Tool / Resource	Function in PSA	Source / Example
Assessment Software	Woods Hole Assessment Model (WHAM)	Flexible, state-space assessment model for integrating diverse data with random effects; ideal for implementing PEV diagnostics.	[59]
Assessment Software	CMSY & BSM (JABBA, JABBA-Select)	Bayesian state-space production models designed for data-limited settings; core engines for integrated PSA.	[64]
Sensitivity Analysis	sensobol R package	Computes variance-based Sobol' indices for GSA.	(Common statistical tool)
Sensitivity Analysis	Custom PWM-GSA Code	Implements the novel PWM-based sensitivity indices for robust GSA.	Based on algorithms in [61]
Uncertainty Analysis	rrisk / pbox R packages	Facilitates modeling with probability boxes (p-boxes) for interval reliability analysis under deep uncertainty.	Related to [63]
Expert Elicitation	Sheffield Elicitation Framework (SHELF)	Structured protocol and R tools for conducting and aggregating expert judgments.	(Widely used framework)
Data & Priors	FishBase / SeaLifeBase	Global database for fish life-history parameters (e.g., L_∞, K, M) essential for priors and productivity scoring.	[65]
Data & Priors	RAM Legacy Stock Assessment Database	Source of meta-analytic relationships for priors on r, K, and other population dynamics parameters.	(Widely used database)

Validating PSA Outcomes: Comparing Results, Frameworks, and Management Utility

Benchmarking PSA Vulnerability Scores Against Independent Stock Status Indicators

This document provides detailed application notes and protocols for the systematic benchmarking of Productivity and Susceptibility Analysis (PSA) vulnerability scores against independent indicators of stock status. Developed within the context of a broader thesis on advancing PSA methodologies for fisheries research, this protocol addresses the critical need to ground semi-quantitative risk assessments in empirical, observational data. The PSA is a semi-quantitative risk assessment tool used to estimate a stock's vulnerability to overfishing by scoring its biological productivity and susceptibility to fishing pressure [1]. While invaluable in data-poor scenarios, its predictive accuracy relative to actual stock condition must be validated [3]. This protocol outlines a standardized framework for pairing PSA outputs with independent status metrics—such as Catch Per Unit Effort (CPUE) trends, length-frequency analyses, and fishery-independent survey data—to calibrate risk thresholds, improve management prioritization, and bolster confidence in using PSA for preliminary assessments and decisions [66].

Productivity and Susceptibility Analysis (PSA) has emerged as a cornerstone methodology for ecological risk assessment in data-limited fisheries worldwide [67]. It functions by scoring a series of attributes related to a species' life history productivity (e.g., growth rate, age at maturity) and its susceptibility to a given fishery (e.g., geographic overlap, management controls) [2]. These scores are synthesized into a composite vulnerability score, which is then categorized (e.g., low, moderate, high) to inform management priorities [2].

However, a PSA score represents an inherent risk of overfishing, not a direct measure of current stock depletion [66]. A stock with high vulnerability may not yet be overfished if fishing pressure is low, while a low-vulnerability stock may be depleted under intense exploitation. This distinction is critical for effective management. Consequently, benchmarking PSA-derived vulnerability against independent indicators of actual stock status is an essential step in the assessment cycle. This process validates the PSA model, refines its risk categorization thresholds, and transitions the assessment from a theoretical risk framework to a more concrete evidence-based evaluation of fishery impacts [3] [67].

This protocol is designed for researchers and fishery scientists operating in resource-limited environments where full, quantitative stock assessments are not feasible. It aligns with the hierarchical approach of frameworks like the Ecological Risk Assessment for the Effects of Fishing (ERAEF), where PSA serves as a mid-tier, semi-quantitative analysis [67].

Protocol I: Conducting a Productivity and Susceptibility Analysis

Objectives and Prerequisites

The objective of this protocol is to derive a standardized vulnerability score for a fish stock or species. Prerequisites include: 1) Defining the specific stock and fishery to be assessed; 2) Assembling a multidisciplinary team with expertise in fishery biology, gear technology, and local fishery operations; 3) Gathering all available biological, ecological, and fishery operational data.

Step-by-Step Methodology

Step 1: Data Collection and Attribute Scoring Collect data for a standardized set of attributes. Each attribute is scored on a scale of 1 (Low Risk/Vulnerability), 2 (Medium Risk), or 3 (High Risk).

• Productivity Attributes (Biological): Score the following life-history parameters [2]:

  • Average age at maturity
  • Average maximum age / longevity
  • Fecundity (average number of offspring)
  • Reproductive strategy (e.g., broadcast spawning vs. live-bearing)
  • Growth rate (e.g., von Bertalanffy k parameter)
  • Natural mortality rate (M)
  • Trophic level
  • Habitat specificity
  • Population size/stock structure

• Susceptibility Attributes (Fishery-specific): Score the following based on the interaction between the stock and the fishery [2]:

  • Geographic overlap between stock and fishery
  • Seasonal overlap
  • Selectivity and encounterability of the fishing gear
  • Post-capture mortality (including discard mortality)
  • Management effectiveness (e.g., existence and enforcement of regulations)
  • Economic incentive/value driving fishing pressure
  • Existence of bycatch reduction devices or other technical measures

Step 2: Data Quality Assessment Concurrently with scoring, assign a Data Quality index (e.g., "Good," "Limited," "Poor," "No Data") for each attribute. This documents uncertainty and highlights critical data gaps [3].

Step 3: Calculation of Composite Scores Calculate the arithmetic mean for all productivity attributes (P) and all susceptibility attributes (S). The overall Vulnerability Score (V) is calculated as the Euclidean distance from the origin: V = sqrt(P² + S²) [3] [2].

Step 4: Categorization Categorize the stock based on its final V score. A common classification scheme is [2]:

• Low Vulnerability: V < 1.8
• Moderate Vulnerability: 1.8 ≤ V ≤ 2.0
• High Vulnerability: V > 2.0

Outputs and Documentation

Primary outputs include: the completed PSA scoring sheet, the calculated P, S, and V scores, the vulnerability category, and a detailed record of data quality for each attribute. Results should be visualized on a two-dimensional plot with P on the y-axis and S on the x-axis, with isoclines representing V scores [2].

Experimental PSA Results: Peruvian Coastal Groundfish Case Study

A 2024 PSA of ten data-poor coastal groundfish species in Peru yielded the following vulnerability scores and categorizations [3]:

Table 1: PSA Vulnerability Scores for Peruvian Coastal Groundfish (2024)

Species	Productivity Score (P)	Susceptibility Score (S)	Vulnerability Score (V)	PSA Category
Broomtail Grouper	-	-	2.57	Extremely High
Grape-eye Seabass	-	-	2.50	Extremely High
Pacific Goliath Grouper	-	-	2.28	Extremely High
Galapagos Sheephead Wrasse	-	-	2.24	Extremely High
Black Snook	-	-	2.13	High
Harlequin Wrasse	-	-	2.05	High
Chino	-	-	2.01	High
Mulata	-	-	1.89	Medium
Pacific Beakfish	-	-	1.89	Medium
Bumphead Parrotfish	-	-	1.73	Low

Note: The study noted that data quality across all stocks ranged from "limited" to "no data," necessitating heavy reliance on expert judgment [3].

Protocol II: Generating Independent Stock Status Indicators

Rationale for Independent Validation

PSA scores indicate potential risk. Independent status indicators provide evidence of actual impact. Using multiple methods with different data streams increases confidence in conclusions, which is especially critical if results are to inform management actions [66]. No single indicator is perfect; a weight-of-evidence approach is recommended.

Selection of Appropriate Methods

Choose methods based on data availability. The FISHE framework recommends using several corroborating methods where possible [66].

Table 2: Independent Stock Status Assessment Methods

Method	Core Data Input	Output Indicator	Key Assumptions & Caveats
CPUE Trend Analysis	Catch and effort data over time (>5 years ideal) [66].	Standardized CPUE as an index of relative abundance.	Assumes catchability is constant. Can be misleading if fisher behavior, gear efficiency, or fish distribution changes [66].
Length-Frequency Analysis	Biological sampling: lengths of fish in the catch.	Indicators of size structure (e.g., % mature, % mega-spawners, mean length).	Assumes the catch is a representative sample of the population. Affected by gear selectivity and discarding practices [66].
Fishery-Independent Surveys	Scientific surveys (e.g., trawl, acoustic, visual).	Absolute or relative biomass/density estimates.	Provides the most direct measure of stock status. Often resource-intensive and expensive [66].
Spatial Comparison	Density data from inside vs. outside Marine Protected Areas (MPAs).	Depletion ratio (Densityopen / DensityMPA).	Assumes habitats are comparable and the MPA is well-enforced and old enough for recovery [66].

Step-by-Step Protocol for Length-Frequency Analysis

As a widely applicable, relatively low-cost method, length analysis is detailed here.

Objective: To estimate exploitation status and detect growth overfishing by analyzing the size structure of the catch.

Materials: Calibrated measuring boards, data recording sheets or electronic devices, species identification guides.

Procedure:

• Stratified Sampling: Develop a sampling plan to collect length measurements randomly across different vessels, gears, seasons, and areas to avoid bias.
• Data Collection: For each specimen, record: species, total length (or standard length, consistent for the species), sex (if possible), and date/location of capture.
• Data Analysis: a. Construct a length-frequency histogram for each major time period (e.g., by year). b. Calculate indicators: * Percentage of mature individuals: Proportion of the catch larger than the known length at maturity (L_m). * Percentage of mega-spawners: Proportion of the catch larger than a target length (e.g., L_opt or a large percentile). * Mean length in the catch: Track changes over time.
• Interpretation: A declining trend in mean length, a low percentage of mature fish, or a vanishing right-hand side (larger sizes) of the length distribution are strong indicators of overfishing [66].

Caveats: This method is less suitable for species with indeterminate growth or complex life histories [66].

Protocol III: Integrated Benchmarking of PSA Against Status Indicators

Benchmarking Workflow

This protocol integrates the outputs from Protocol I (PSA) and Protocol II (Status Indicators) to test the correlation between predicted vulnerability and observed stock status.

Step 1: Tabulate Paired Data For each stock assessed, create a master table with the calculated PSA Vulnerability Score (V) and the results from each independent status indicator (e.g., CPUE trend direction, depletion ratio, % mature in catch).

Step 2: Perform Comparative Analysis

• Visual Analysis: Plot each stock on a 2D graph: for example, PSA V score on the x-axis and a status indicator (e.g., biomass depletion) on the y-axis. Look for a positive correlation where higher V is associated with greater depletion.
• Threshold Testing: Test the standard PSA risk categories (e.g., Low: V<1.8, High: V>2.0) against status categories (e.g., "Healthy," "Overfished"). Calculate the rate of true positives, false positives, true negatives, and false negatives to assess the predictive power of the PSA thresholds.

Step 3: Calibration and Refinement Based on the analysis, the PSA risk categorization thresholds may be adjusted for the local fishery context. For instance, if most stocks with V > 1.9 show clear signs of depletion, the "High Vulnerability" threshold could be lowered accordingly.

Step 4: Management Prioritization Synthesis Generate a final, integrated prioritization matrix. Stocks with high PSA vulnerability AND poor status indicators constitute the highest priority for urgent management intervention. Stocks with high PSA vulnerability but currently healthy status require precautionary management to prevent future decline.

Integrated Assessment Workflow Diagram

Diagram Title: PSA Benchmarking and Management Prioritization Workflow

Case Study: Benchmarking Application

The Peruvian groundfish study [3] performed an initial benchmark by contrasting its PSA vulnerability scores (Table 1) with the limited available information on stock statuses. The authors concluded that the scores were "mostly consistent with contrasted information" and that standard PSA thresholds (e.g., V > 2.0 for high vulnerability) remained useful for informing management in that data-poor context. This process validated the PSA's identification of high-priority species like the Broomtail Grouper (V=2.57).

The Scientist's Toolkit: Essential Research Reagents & Materials

Table 3: Essential Research Toolkit for PSA Benchmarking Studies

Item / Category	Function / Purpose	Key Considerations
Standardized PSA Scoring Sheets	Provides a consistent framework for scoring productivity and susceptibility attributes across multiple assessors and stocks [2].	Should include columns for scores, data quality indices, and notes on evidence sources.
Expert Elicitation Protocols	Formalizes the process of gathering qualitative data and expert judgment to fill critical data gaps in scoring [3].	Reduces bias through structured interviews and Delphi techniques.
Biological Sampling Kit	For collecting length-frequency and other biological data. Includes calipers/measuring boards, digital scales, dissection tools, and tissue sample storage (e.g., vials, ethanol) [66].	Calibration of tools is essential. Requires standardized species ID guides.
Fishery Logbook/Datasheets	Structured forms for fishers or observers to record catch, effort, location, and gear specifics for CPUE calculation [66].	Design must be simple, clear, and logistically feasible for the fishery context.
Data Analysis Software	For statistical analysis and visualization (e.g., R, Python with packages like fishboot, TropFishR; or general tools like SPSS, PRISM).	Expertise in standardizing CPUE (e.g., using GLMs) and analyzing length-frequency distributions is required.
Scientific Survey Equipment	For generating independent biomass estimates (e.g., calibrated trawls, underwater video systems, acoustic equipment) [66].	Highly resource-intensive. Often involves specialized vessels and personnel.
Reference Databases	Online life history databases (e.g., FishBase, SeaLifeBase) and local scientific literature to inform productivity attribute scores [2].	Local studies should be prioritized over global averages to account for regional variation.

Data Visualization for Accessible Communication

Effective communication of benchmarking results is crucial. Adhere to the following principles derived from data visualization and accessibility best practices [68] [69]:

• Contrast for Clarity: Ensure a minimum 3:1 contrast ratio between adjacent data elements (e.g., bars in a chart) and a 4.5:1 ratio between text and its background [69]. The specified Google palette (Blues, Reds, Yellows, Greens, Grays, White) provides a strong basis for this.
• Do Not Rely on Color Alone: Use patterns (hatching), shapes, or direct data labels in addition to color to convey information [69]. This makes visualizations accessible to color-blind individuals and robust in black-and-white print.
• Recommended Chart Types:
  • Scatter Plot: The primary tool for benchmarking, plotting PSA V score against a status metric (e.g., % biomass depletion).
  • Bar Chart: To compare final vulnerability scores or status indicators across multiple species (as in Table 1).
  • Heat Map/Highlight Table: To visually summarize the multi-method, multi-species assessment matrix, using color intensity to show level of risk or concern [68].
• Provide Alternative Data Formats: Always accompany key charts with accessible data tables or descriptive text to ensure the information is available to all users [69].

This document provides a detailed methodological comparison between Productivity and Susceptibility Analysis (PSA) and traditional stock assessment models within a fisheries research framework. It outlines their respective applications, detailed protocols, and essential toolkits to guide researchers in model selection and implementation.

Productivity and Susceptibility Analysis (PSA) is a semi-quantitative, risk-based framework designed for the rapid vulnerability assessment of data-poor species [1]. It evaluates a stock's intrinsic capacity to recover (Productivity) and its exposure to a fishery (Susceptibility) [70]. In contrast, traditional stock assessment models are quantitative analytical tools that use time-series data to estimate population size, fishing mortality, and sustainable catch limits for data-moderate to data-rich stocks [71].

Table 1: Fundamental Comparison of PSA and Traditional Stock Assessment Models

Aspect	Productivity and Susceptibility Analysis (PSA)	Traditional Stock Assessment Models (e.g., Statistical Catch-at-Age)
Primary Objective	Rapid, relative risk ranking to prioritize species for management or further assessment [72] [1].	Quantitative estimation of absolute stock status (biomass, mortality) to determine sustainable harvest levels [73].
Data Requirements	Low. Relies on life history parameters (e.g., max age, fecundity) and general fishery interaction attributes. Tolerates qualitative or missing data [9] [3].	High. Requires long-term, high-quality data: catch, abundance indices, age/length compositions, and fishery selectivity [71] [73].
Methodological Core	Expert-informed scoring (1-3) of pre-defined attributes. Calculates a vulnerability score (V = √(P² + S²)) [70].	Statistical fitting of a population dynamics model to observed data to estimate key biological and fishery parameters.
Key Output	Relative vulnerability score and risk category (Low, Medium, High) [70].	Absolute metrics: stock biomass relative to targets/limits (B/BMSY), fishing mortality rate (F/FMSY), and overfishing status [73].
Temporal Dynamics	Static "snapshot" based on general life history and fishery characteristics.	Explicitly models population changes over time in response to fishing and recruitment.
Validation Evidence	Can be precautionary. One study showed a 50% misclassification rate vs. quantitative assessments [72].	Considered the "best scientific information available" when data is sufficient, validated through peer review [73].
Ideal Application Context	Data-poor stocks, bycatch species, initial screening of many species, and small-scale or developing fisheries [9] [3].	Assessed target species in managed fisheries where investment in data collection is feasible and required for precise management [71].

Detailed Experimental Protocols

Protocol for Conducting a Productivity and Susceptibility Analysis (PSA)

This protocol is adapted from the Hobday et al. framework and contemporary applications for data-poor fisheries [9] [3] [70].

Step 1: Definition of Scope and Formation of Expert Panel

• Define the spatial scale (e.g., regional, local), the fishery/gears assessed, and the list of species [9].
• Assemble a panel including fisheries biologists, ecologists, and local fishing experts. Local knowledge is critical for accurately scoring susceptibility attributes in data-poor contexts [9].

Step 2: Data Collation and Quality Grading

• For each species, compile available information for standard PSA attributes. Productivity attributes often include maximum age/size, growth rate (k), natural mortality (M), fecundity, and age at maturity. Susceptibility attributes include spatial/temporal overlap with the fishery, catchability, post-capture mortality, and management effectiveness [70].
• Assign a Data Quality score (e.g., 1=High, 2=Medium, 3=Low/No Data) for each attribute based on source reliability [3].

Step 3: Attribute Scoring and Weighting

• Score each attribute on an ordinal risk scale: 1 (Low Risk/Vulnerability), 2 (Medium Risk), 3 (High Risk) using predefined breakpoints [70]. For missing data, a precautionary score of 3 is often assigned, or expert judgment is used [3].
• (Optional) Apply weighting to attributes if certain factors are deemed more influential for a specific fishery. Note that subjective weighting can increase outcome uncertainty [53].

Step 4: Calculation of Vulnerability

• Calculate the mean Productivity (P) score (arithmetic mean of productivity attributes).
• Calculate the mean Susceptibility (S) score (geometric mean of susceptibility attributes is common).
• Compute the overall Vulnerability Score: ( V = \sqrt{P^2 + S^2} ) [70].
• Classify species into risk categories based on V-score thresholds (e.g., Low: V < 2.5, Medium: 2.5 ≤ V ≤ 3.0, High: V > 3.0) [3].

Step 5: Sensitivity Analysis and Reporting

• Conduct sensitivity analysis by varying scores for data-poor attributes to test the robustness of the risk ranking.
• Document all scores, data sources, quality grades, and expert rationale transparently. The results should prioritize high-vulnerability species for management action or further data collection [53].

Diagram 1: Workflow for a standard PSA.

Protocol for a Traditional Benchmark Stock Assessment

This protocol outlines key phases, as demonstrated in the 2025 American Lobster Benchmark Stock Assessment [73].

Phase 1: Pre-Assessment Preparation and Data Compilation

• Define Stock Structure: Establish assessment units based on biology and fishery (e.g., Gulf of Maine/Georges Bank vs. Southern New England lobster stocks) [73].
• Compile Core Data Time Series: Gather and standardize critical datasets:
  • Fishery-Dependent: Total catch/landings (weight), discards, effort data.
  • Fishery-Independent: Scientific survey indices of abundance by age/size.
  • Biological Samples: Age and length compositions from catches and surveys.
  • Auxiliary Data: Environmental covariates (e.g., water temperature) [73].

Phase 2: Model Selection, Configuration, and Fitting

• Select a state-space statistical catch-at-age model (e.g., Stock Synthesis). The model mathematically represents population processes: recruitment, growth, natural mortality, and fishing mortality [71].
• Configure the Model: Define initial conditions, selectivity curves for each fishery/survey, and stock-recruitment relationships.
• Fit the Model: Use statistical maximum likelihood or Bayesian methods to estimate parameters (e.g., annual recruitment, fishing mortality, stock biomass) that best explain the observed catch, survey, and composition data.

Phase 3: Estimation of Status and Reference Points

• Calculate stock status: Estimate current (e.g., 2021-2023 average) spawning stock biomass and fishing mortality rate [73].
• Compare estimates to biologically based reference points (e.g., BMSY, FMSY) to determine status: Is the stock depleted? Is overfishing occurring? [73].

Phase 4: Retrospective Analysis, Projection, and Peer Review

• Perform retrospective analysis to check for systematic trends in model estimates, which may indicate model misspecification.
• Conduct management strategy evaluations (MSEs) to project future stock outcomes under different harvest policies [73].
• Submit the entire assessment report and analysis code for formal external scientific peer review before use in management [73].

Diagram 2: Phases of a traditional benchmark stock assessment.

The Scientist's Toolkit: Essential Research Reagents and Materials

Table 2: Essential Toolkit for Comparative Stock Assessment Research

Tool/Reagent Category	Specific Item or Software	Primary Function in Assessment
PSA-Specific Frameworks	Modified Hobday PSA worksheets; IATTC PSA protocol [53] [70].	Provides standardized attribute lists, scoring breakpoints, and calculation formulas.
Expert Elicitation Tools	Structured survey questionnaires (e.g., online Delphi surveys); Workshop facilitation guides [9].	Systematically captures qualitative expert knowledge for scoring PSA attributes, especially in data-poor settings.
Statistical Stock Assessment Packages	Stock Synthesis (SS3), A Stock Assessment Model (ASAM), Extended Survival Analysis (XSA) [71].	Provides pre-built, peer-reviewed population dynamics models for rigorous quantitative assessment.
Programming & Statistical Environments	R (with packages r4ss, ggplot2), AD Model Builder, Template Model Builder (TMB), Python (with pyStockAssessment).	Core platforms for data processing, model fitting, statistical analysis, and generating diagnostic plots.
Data Management & Version Control	Git/GitHub/GitLab, SQL databases, RShiny/Dash applications.	Ensures reproducibility, tracks model code changes, manages complex datasets, and creates interactive data summaries.
High-Performance Computing (HPC)	Access to HPC clusters or cloud computing (AWS, Google Cloud).	Enables complex model runs, extensive sensitivity analyses, and management strategy evaluation simulations.
Reference Material & Metadata	FishBase/SeaLifeBase; Previous assessment reports; Data quality standards (ICES, NOAA).	Sources for life history priors; context for model configuration; ensures data integrity and comparability.

Application Notes: Contextual Selection and Integration

Selecting PSA: Choose PSA for scoping studies, assessing bycatch in multispecies fisheries, or evaluating stocks in small-scale or developing fisheries where quantitative data is extremely limited [3]. Its strength lies in rapid prioritization, not precise status determination. Researchers must acknowledge its tendency toward precautionary misclassification [72] and the sensitivity of results to expert judgment and attribute weighting [53].

Selecting Traditional Models: Employ traditional assessments for key commercial stocks where reliable catch, survey, and biological data exist, and where management requires precise estimates of sustainable catch [73]. This approach is resource-intensive but provides the foundation for legally defensible harvest control rules.

Integrative Approach: A hierarchical or tiered approach is scientifically defensible. Use PSA to screen a large number of species and identify those at potentially high risk. These high-priority species can then be targeted for enhanced monitoring and data collection, with the ultimate goal of graduating them to a quantitative assessment framework when resources allow [72] [70]. This aligns with the core thesis of using PSA as a productivity tool within research programs to strategically focus limited assessment resources.

Productivity and Susceptibility Analysis (PSA) is a semi-quantitative risk assessment method used to evaluate the vulnerability of fishery species or stocks when conventional, data-intensive stock assessments are not available [1]. It operates by scoring a standardized set of life history attributes (productivity) and fishery interaction attributes (susceptibility) to determine a composite risk score [70]. This methodology is a cornerstone of precautionary management, particularly within the Marine Stewardship Council's (MSC) Risk-Based Framework (RBF), where it functions as a critical tool for assessing data-limited fisheries [74]. Furthermore, its principles are increasingly integrated into Ecosystem-Based Fisheries Management (EBFM), providing a mechanism to assess cumulative impacts and inform holistic management strategies that consider entire socio-ecological systems rather than single species in isolation [75] [76] [77].

PSA within the MSC Risk-Based Framework (RBF)

The MSC Risk-Based Framework is designed to maintain program accessibility for data-limited fisheries, including small-scale and developing world fisheries, where the full, quantitative requirements of the default MSC Fisheries Standard present a barrier [74].

Role and Application of PSA in the RBF

Within the RBF, PSA is one of four complementary assessment tools. Its specific role is to assess stock recovery likelihood and the interaction probability of non-target species with fishing gear [74]. The RBF is applied selectively: a certification body may use it for specific performance indicators where insufficient data exists for default assessment, while using the standard method for all other indicators [74]. The decision to employ the RBF must be announced at the assessment's outset, allowing for stakeholder comment [74]. It is important to note that the RBF does not cover the assessment of fishery management systems (MSC Principle 3), which are evaluated using the standard method regardless of data limitations [74].

The RBF Toolkit and Workflow

The framework employs a suite of tools, each addressing different ecological components [74]:

• Consequence Analysis (CA): Uses available data to assess trends in target stocks.
• Productivity Susceptibility Analysis (PSA): Assesses vulnerability of target and non-target species.
• Consequence Spatial Analysis (CSA): Identifies potential impacts of fishing on habitats.
• Scale Intensity Consequence Analysis (SICA): Aims to identify the fishery's impacts on the wider ecosystem.

The following workflow diagram illustrates how these tools integrate within the MSC assessment process for a data-limited fishery.

Diagram 1: MSC Risk-Based Framework assessment workflow.

Core PSA Methodology and Protocol

The foundational PSA protocol involves scoring a series of attributes related to a species' inherent capacity to recover (productivity) and its exposure to a fishery (susceptibility).

Attribute Scoring and Calculation

The standard PSA evaluates a species against a fixed set of attributes [1] [70]. Each attribute is assigned a risk score from 1 (low risk) to 3 (high risk).

Table 1: Core PSA Attributes and Scoring Criteria

Category	Attribute	Score 1 (Low Risk)	Score 2 (Medium Risk)	Score 3 (High Risk)
Productivity	Age at Maturity	< 5 years	5 - 15 years	> 15 years
	Fecundity	> 20,000 eggs/year	100 - 20,000 eggs/year	< 100 eggs/year
	Natural Mortality (M)	High	Medium	Low
	Reproductive Strategy	Broadcast spawner, low parental investment	Intermediate	Live bearer, high parental investment
Susceptibility	Spatial Overlap	Low overlap with fishery	Moderate overlap	High/Complete overlap
	Catchability	Low vulnerability to gear	Moderate vulnerability	High vulnerability
	Post-Capture Mortality	High survival if released	Medium survival	Low survival (always dies)
	Management Effectiveness	Effective management & enforcement	Moderate management	Ineffective or no management

The overall Productivity (P) score is the arithmetic mean of all productivity attribute scores. The overall Susceptibility (S) score is typically the geometric mean of all susceptibility attribute scores, reflecting the multiplicative nature of risk factors [70]. The final Vulnerability (V) score is calculated as the Euclidean distance from the origin: V = √(P² + S²) [70]. Species are then classified as Low, Medium, or High risk based on thresholds applied to the V score.

Extended PSA (ePSA) for Targeted Stocks

An extended PSA (ePSA) has been developed for application to targeted stocks, incorporating additional and re-weighted attributes [70]. It includes up to 10 productivity and 12 susceptibility attributes, the latter grouped into catchability and management sub-categories. A key feature of the ePSA is a weighting system (default weight of 2 for all attributes, adjustable based on expert judgment) and a data quality scoring system to account for information uncertainty [70].

Table 2: Comparison of Standard PSA (sPSA) and Extended PSA (ePSA)

Feature	Standard PSA (sPSA)	Extended PSA (ePSA)
Primary Use	Bycatch species assessment [70]	Targeted stock assessment [70]
Number of Attributes	7 Productivity, 4 Susceptibility [70]	10 Productivity, 12 Susceptibility [70]
Scoring for Missing Data	Precautionary: scored as high risk (3) [70]	Attribute omitted from analysis [70]
Weighting	Attributes equally weighted [70]	Adjustable weighting system [70]
Data Quality Metric	Not formally included [70]	Formal data quality score accompanies result [70]

PSA in Ecosystem-Based Fisheries Management (EBFM)

EBFM is a holistic management approach that considers interactions within an entire ecosystem, including human communities, rather than managing single species in isolation [75] [77]. Its goal is to maintain ecosystems in a healthy, productive, and resilient condition [75].

PSA as a Bridge to EBFM

PSA supports EBFM by providing a rapid, risk-based tool to evaluate a wide range of species impacted by fishing, which is essential for meeting legal mandates to minimize bycatch and protect ecosystem structure [70]. It aligns with the EBFM guiding principles of addressing cumulative impacts and considering trade-offs among ecosystem components [77]. PSA outputs can feed directly into Integrated Ecosystem Assessments (IEAs), which are structured scientific processes supporting EBFM by scoping issues, defining indicators, evaluating risks, and monitoring outcomes [75] [77].

Protocol for Cumulative Impact Assessment Using PSA

A critical advancement is the modification of PSA to assess the cumulative risk posed by multiple, co-occurring fisheries to a shared species or ecosystem. This addresses a key EBFM objective [76].

Experimental Protocol: Cumulative Risk Assessment for Multiple Fisheries

Objective: To quantify the aggregated vulnerability of a marine species to all fisheries operating within its range.

Materials: Life history data for the species; spatial, effort, and gear data for all relevant fisheries; PSA scoring sheets.

Procedure:

• Define the System: Identify the geographic region and all fisheries (f=1...n) operating within it that interact with the target species.
• Conduct Individual PSAs: Perform a separate, conventional PSA for the target species against each individual fishery (f). This yields a susceptibility score (S_f) for each fishery.
• Calculate Aggregate Susceptibility (AS): For the target species, calculate the Aggregate Susceptibility as the Euclidean mean of the susceptibility scores from all fisheries: AS = √(Σ(Sf)² / n) [76]. *Rationale:* The Euclidean mean ensures that the addition of any new fishery (with *Sf* > 0) increases the overall AS, properly reflecting cumulative pressure [76].
• Determine Cumulative Vulnerability: Combine the species' Productivity score (P) (which is constant across fisheries) with the Aggregate Susceptibility (AS) to calculate the Cumulative Vulnerability score: V_cum = √(P² + AS²).
• Compare and Interpret: Compare the cumulative vulnerability ranking (V_cum) with the rankings from single-fishery PSAs. Species for which risk classification shifts significantly (e.g., from Medium to High) only under the cumulative assessment are priority candidates for cross-fishery management attention [76].

The following diagram illustrates this cumulative assessment logic.

Diagram 2: Logic of cumulative PSA for multiple fisheries.

Application Note: A study applying this method to 81 species in Baja California fisheries found the proportion of species classified as "high risk" increased from 22.2% (under single-fishery analysis) to 38.3% when cumulative effects were considered, demonstrating the critical importance of the aggregated approach [76].

Conducting a robust PSA requires specific data inputs and analytical resources. The following table details key "research reagent solutions" for implementing the protocols described.

Table 3: Research Reagent Solutions for PSA Implementation

Item / Resource	Function / Purpose	Application Notes
Life History Parameter Database (e.g., FishBase, SeaLifeBase)	Provides standardized data for scoring productivity attributes (age at maturity, fecundity, growth parameters).	Essential for initial scoring. Data quality from these sources should be noted for ePSA data scoring [70].
Fishery Observer Data & Logbooks	Provides critical information on spatial overlap, catch composition, gear interactions, and post-capture mortality for susceptibility scoring.	Key for estimating spatial overlap and catchability attributes.
Geographic Information System (GIS)	Analyzes spatial overlap between species distribution (from surveys or models) and fishery effort layers.	Used to quantify the "Spatial Overlap" susceptibility attribute.
Expert Elicitation Framework (e.g., structured workshop, Delphi method)	Elicits quantitative or semi-quantitative judgments to score attributes when empirical data are lacking.	A formal, documented process is crucial for transparency and reproducibility within the RBF and EBFM contexts [74].
PSA Calculation Software / Script (e.g., R package psa, custom spreadsheet)	Automates the calculation of Productivity (P), Susceptibility (S), and Vulnerability (V) scores from attribute inputs.	Should be configured to allow for both arithmetic (P) and geometric (S) means, and Euclidean distance (V) calculation [70].
Colorblind-Friendly Visualization Palette	Ensures that risk matrices and output diagrams are accessible to all audiences.	Use palettes like Okabe-Ito or Tableau's built-in scheme. Avoid red/green/brown/orange combinations for critical distinctions [13] [78].

Productivity and Susceptibility Analysis (PSA) is a semi-quantitative risk assessment framework designed to evaluate the vulnerability of fish stocks to overfishing, particularly in data-limited or data-poor contexts [1]. The method is founded on assessing two core components: a stock's intrinsic productivity (its biological capacity for growth and recovery) and its external susceptibility (its exposure and sensitivity to a fishery) [1] [5]. By scoring a series of predefined attributes for productivity and susceptibility, managers can derive a composite vulnerability score, which serves to rank species by their relative risk and inform subsequent conservation priorities and regulatory actions [3] [5].

This framework has become a cornerstone of the MSC Risk-Based Framework (RBF) and similar hierarchical ecological risk assessment systems [1]. Its primary utility lies in its ability to rapidly screen a large number of species—including target, bycatch, and data-poor stocks—to identify those most in need of detailed assessment or immediate management intervention [3] [5]. The following application notes and protocols detail the implementation of PSA, its role in shaping management outcomes, and its integration within a broader scientific thesis on sustainable fisheries.

Core Methodological Framework and Application Notes

Foundational Protocol: The PSA Scoring System

The standard PSA protocol involves scoring a series of biological and fishery-specific attributes on a scale from 1 (low risk) to 3 (high risk) [1] [5].

• Productivity Attributes (Biological): These reflect the stock's intrinsic capacity for population growth and recovery. Common attributes include:

  • Average Age at Maturity: Species maturing later (>15 years) receive higher risk scores (3) [5].
  • Fecundity: Species with lower reproductive output (<100 eggs/year) score higher risk (3) [5].
  • Average Maximum Size, Growth Rate, Natural Mortality, and Trophic Level.

• Susceptibility Attributes (Fishery-specific): These measure the stock's exposure and sensitivity to the fishery.

  • Availability: Overlap of stock distribution with fishing grounds.
  • Encounterability: Behavior (e.g., position in water column) making encounter with gear likely.
  • Selectivity: Match between gear and stock size.
  • Post-Capture Mortality: Survival rate after capture [1] [5].

Calculation of Vulnerability:

• The overall Productivity score (P) is the arithmetic mean of all productivity attribute scores.
• The overall Susceptibility score (S) is the geometric mean of all susceptibility attribute scores [5].
• The final Vulnerability score (V) is calculated as the Euclidean distance from the origin: V = √(P² + S²) [5]. This yields a score ranging from approximately 1.41 (lowest risk) to 4.24 (highest risk). Thresholds are then applied to categorize risk as Low, Medium, or High [5].

Application Note: PSA in Extremely Data-Poor Contexts

In scenarios with severe data limitations, the PSA framework can be adapted to integrate expert knowledge and handle multiple gear types. A 2024 study of Peruvian coastal groundfish demonstrated this adaptive approach [3].

Protocol for Data-Poor Adaptation:

• Problem Scoping: Identify target species and all interacting fisheries/gears.
• Data Inventory and Gap Analysis: Systematically review available biological, catch, and effort data for each attribute. Document data quality (e.g., "limited," "no data") [3].
• Structured Expert Elicitation: Convene a panel of experts (biologists, fishery managers, fishers) to fill data gaps. Use a standardized questionnaire to score each PSA attribute, documenting the rationale and confidence for each score [3].
• Susceptibility Aggregation: For stocks impacted by multiple gears, calculate a susceptibility score per attribute for each gear, then integrate into a single composite score for management purposes [3].
• Validation and Threshold Calibration: Contrast PSA results with any independent, albeit limited, indicators of stock status (e.g., CPUE trends, anecdotal reports from fishers) to evaluate the accuracy of risk categorization thresholds [3].

Table 1: PSA Vulnerability Assessment of Peruvian Coastal Groundfish (Adapted from [3])

Species (Common Name)	Vulnerability Score (V)	Risk Category	Primary Management Implication
Broomtail grouper	2.57	Extremely High	Priority for immediate regulatory action (e.g., catch limits, spatial closures).
Grape-eye seabass	2.50	Extremely High	Priority for immediate regulatory action.
Pacific goliath grouper	2.28	Extremely High	Priority for immediate regulatory action and conservation listing.
Galapagos sheephead wrasse	2.24	Extremely High	Priority for immediate regulatory action.
Black snook	2.13	High	Priority for data collection and precautionary management.
Harlequin wrasse	2.05	High	Priority for data collection and precautionary management.
Chino	2.01	High	Priority for data collection and precautionary management.
Mulata	1.89	Medium	Requires monitoring; candidate for inclusion in future management plans.
Pacific beakfish	1.89	Medium	Requires monitoring.
Bumphead parrotfish	1.73	Low	Lower management priority; maintain status quo with observation.

Workflow Visualization: From Data to Management Action

The following diagram outlines the logical workflow for conducting a PSA and translating its results into regulatory and conservation outcomes.

PSA to Management Decision Workflow

Linking PSA Outcomes to Regulatory and Conservation Priorities

Informing Hierarchical Management Frameworks

PSA is explicitly designed as a tiered risk assessment tool. Its results directly dictate the level of management resource allocation:

• Low-Vulnerability Stocks: Maintain under general monitoring; no urgent regulatory action required [3].
• Medium-Vulnerability Stocks: Flag for enhanced monitoring and potential inclusion in developing fishery management plans (FMPs) [3].
• High to Extremely High-Vulnerability Stocks: Trigger immediate Priority Actions:
  • Development of Targeted Regulations: Implementation of species-specific catch limits (quotas), size limits, or gear restrictions [3].
  • Spatial Management: Establishment of time-area closures to protect critical habitats (spawning aggregations, nursery grounds) or reduce fishing pressure [3] [79].
  • Catalyzing Formal Stock Assessment: PSA outcomes provide the justification for investing in more complex, data-intensive quantitative stock assessments [5].

Setting Scalable Biodiversity Conservation Targets

PSA-derived vulnerability rankings can be integrated with systematic conservation planning to maximize ecological and climate benefits. Research on scalable priorities in Asia demonstrates a model for using species risk data to inform protected area network design [79].

Protocol for Integrating PSA into Conservation Planning:

• Spatial Overlay: Map the distributions of high and very high-vulnerability species identified by PSA.
• Identification of Synergy Zones: Use spatial analysis (e.g., in GIS software) to identify areas where high-vulnerability species hotspots overlap with other critical ecological values, such as high carbon sequestration ecosystems (e.g., mangrove forests, seagrass beds) [79].
• Target Setting: Establish explicit, measurable conservation targets. For example: "Protect 30% of the habitat for 'Extremely High' vulnerability species within the next 5 years" or "Designate marine protected areas (MPAs) in the top 20% of synergy zones identified between vulnerability and carbon storage." [79].
• Nested Scaling: Define complementary targets at regional, national, and local scales to ensure ecological representativeness and practical implementability [79].

Table 2: Translating PSA Risk Categories into Conservation & Regulatory Actions

PSA Risk Category	Immediate Regulatory Response	Medium-Term Conservation Priority	Research & Monitoring Mandate
Extremely High / High	Implement immediate catch or effort controls; consider emergency listing on protected species ordinances.	Designate as a priority species for spatial protection (MPAs, habitat reserves); integrate into regional biodiversity strategy.	Mandate a full stock assessment; establish a dedicated fishery-independent monitoring program.
Medium	Include in broader fishery management plans with precautionary reference points; implement data reporting requirements.	Include as secondary priority in MPA network design; target for habitat restoration projects.	Fund life history studies; improve catch and effort data collection.
Low	General compliance with existing, non-specific fishery regulations.	Consider in broad ecosystem-based management strategies.	Continue routine surveillance within fishery-dependent data streams.

Critical Evaluation and Limitations of PSA

While PSA is a valuable screening tool, its assumptions and predictions require careful consideration. A 2018 quantitative evaluation mapped PSA logic to age-structured population models and found that its qualitative assumptions could be inappropriate, leading to poor performance under a range of exploitation scenarios [5]. Key limitations include:

• Oversimplification of Dynamics: The framework does not explicitly account for population dynamics, stock recruitment relationships, or density-dependent processes [5].
• Arbitrary Scoring Thresholds: The breakpoints between low/medium/high scores for attributes (e.g., age at maturity) are often generalized and may not be appropriate for all ecosystems or taxonomic groups [5].
• Static Nature: A basic PSA presents a snapshot of risk but does not project future risk under alternative management scenarios.

Therefore, PSA is best used as a prioritization and communication tool, not as a definitive assessment of stock status. It efficiently identifies which stocks require more sophisticated, dynamic simulation modeling (e.g., Management Strategy Evaluation - MSE) for robust management advice [5].

The Scientist's Toolkit: Essential Research Reagents & Materials

Conducting a rigorous PSA and translating it into effective policy requires a suite of conceptual and material tools.

Table 3: Research Toolkit for PSA and Management Evaluation

Tool / Reagent Category	Specific Item or Protocol	Function & Application
Data Collection & Synthesis	Standardized PSA Attribute Spreadsheet	A template for compiling biological parameters (max size, growth rate) and fishery metrics for all scored species [3] [1].
	Historical Catch & Effort Database	Time series data for analyzing trends and calibrating susceptibility scores related to fishery overlap and intensity [3].
Expert Elicitation	Structured Interview/Delphi Protocol	A formal questionnaire and process for systematically gathering and reconciling expert judgments to score attributes where empirical data are lacking [3].
	Confidence Scoring Rubric	A scale (e.g., 1-5) for experts to indicate their certainty in each provided score, allowing for uncertainty analysis [3].
Spatial Analysis	Geographic Information System (GIS) Software	To map species distributions, fishing effort, and habitat data to assess spatial susceptibility and design spatial management [79].
	Species Distribution Models (SDMs)	To predict range and habitat suitability for data-poor species, informing availability and encounterability scores.
Decision Support & Visualization	Risk Plot (P vs. S Scatterplot)	A scatterplot graphing Productivity (P) against Susceptibility (S) to visually cluster species by risk category and identify key drivers [5].
	Conservation Prioritization Software (e.g., Marxan, Zonation)	To identify optimal networks of protected areas that meet biodiversity (high-vulnerability species) and carbon storage targets efficiently [79].

Integrated Protocol for Evaluating Management Outcomes

This final protocol outlines a comprehensive approach for using PSA not just for initial risk ranking, but for the ongoing evaluation of management effectiveness.

Title: Post-Implementation Monitoring and Evaluation of PSA-Informed Management. Objective: To assess whether management actions triggered by PSA risk categorizations are leading to improved stock status and reduced conservation risk.

Procedure:

• Establish a Monitoring Baseline: At the time of PSA completion and management implementation (T₀), document:
  • The PSA vulnerability score and risk category.
  • The specific management measures enacted (e.g., quota level, MPA boundaries).
  • Available indicators of stock status (e.g., mean size in catch, CPUE, spawning biomass index if available).

• Define Evaluation Metrics and Timelines:

  • Short-term (1-3 years): Monitor compliance with new regulations (e.g., catch compliance in logbooks, VMS data for closures).
  • Medium-term (3-5 years): Track changes in fishery-dependent indicators (e.g., trends in standardized CPUE for the stock).
  • Long-term (5+ years): Re-evaluate key PSA productivity attributes through targeted biological sampling (e.g., has the average age/size in the catch increased?).

• Conduct a Follow-up PSA (at 5-year interval): Repeat the PSA scoring process using updated data.

  • Hypothesis: Effective management should reduce Susceptibility scores (e.g., through improved selectivity or reduced post-capture mortality) and may eventually lead to improvements in Productivity proxies (e.g., increased average size).
  • Success Criterion: A reduction in the overall Vulnerability score (V) and/or a downgrade in the risk category (e.g., from High to Medium).

• Adaptive Management Feedback Loop: Integrate the evaluation results into the management cycle.

  • If vulnerability remains high or increases, investigate causes (e.g., non-compliance, inadequate measures, environmental factors) and strengthen management measures.
  • If vulnerability decreases, consider the potential for cautiously relaxing measures or reallocating monitoring resources to higher-priority stocks.

The integration of this evaluation protocol ensures that PSA is embedded within a dynamic, evidence-based management framework, directly linking scientific assessment to regulatory action and ultimately, to measurable conservation outcomes.

Conclusion

Productivity and Susceptibility Analysis stands as an indispensable, flexible tool for proactive risk assessment in the world's numerous data-limited fisheries. By synthesizing biological productivity and fishery exposure, PSA effectively prioritizes vulnerable species, from Peruvian groundfish to Mediterranean elasmobranchs, for urgent management and research. The methodology's strength is enhanced by integrating fishers' knowledge and transparently addressing data quality issues. While not a replacement for full quantitative assessments, validated PSA outcomes provide a critical scientific basis for implementing precautionary measures, guiding monitoring efforts, and advancing toward ecosystem-based fisheries management. Future directions should focus on standardizing attribute thresholds, expanding its integration with climate vulnerability assessments, and strengthening its formal role in regional fishery management organization decision-making processes to safeguard marine biodiversity and fishery sustainability.

References

• 1. Productivity and Susceptibility Analysis [PSA] [https://www.iss-foundation.org/glossary/productivity-and-susceptibility-analysis/]
• 2. Productivity and Susceptibility Analysis - FISHE [https://fishe.edf.org/node/95]
• 3. Productivity susceptibility analysis in extremely data-poor ... [https://www.sciencedirect.com/science/article/pii/S0165783624001085]
• 4. Productivity and Susceptibility Analysis PSA | EDF Fishery [https://fisherysolutionscenter.edf.org/tools/productivity-and-susceptibility-analysis-psa]
• 5. A quantitative evaluation of a qualitative risk assessment... | PLOS One [https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0198298]
• 6. May life history traits be used as indices of population ... [https://www.sciencedirect.com/science/article/abs/pii/S1385110100000411]
• 7. Life histories determine divergent population trends for ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC7428017/]
• 8. FAO releases the most detailed global assessment of ... [https://www.fao.org/newsroom/detail/fao-releases-the-most-detailed-global-assessment-of-marine-fish-stocks-to-date/en]
• 9. Integration of scientific and local expertise to develop risk ... [https://www.sciencedirect.com/science/article/pii/S0165783621002812]
• 10. Cetacean visual observations using line-transect survey ... [https://www.fisheries.noaa.gov/inport/item/65919]
• 11. "Addressing the Productivity-Susceptibility Analysis in Data ... [https://digitalcommons.odu.edu/biology_fac_pubs/643/]
• 12. Visual analysis of oceanic data for marine ecosystems [https://www.sciencedirect.com/science/article/pii/S1574954124003042]
• 13. 5 Tips on Designing Colorblind-Friendly Visualizations [https://www.tableau.com/blog/examining-data-viz-rules-dont-use-red-green-together]
• 14. Accessible Color Palette Generator | WCAG Compliant [https://venngage.com/tools/accessible-color-palette-generator]
• 15. Contrast Checker [https://webaim.org/resources/contrastchecker/]
• 16. Population Assessments and Fish Stocks [https://www.fisheries.noaa.gov/topic/population-assessments/fish-stocks]
• 17. Step 1: Specify the Fishery System [https://climateresilientfisheries.net/step-1-specify-the-fishery-system/]
• 18. Guidelines for the Use of Fishes in Research [https://fisheries.org/policy-media/science-guidelines/guidelines-for-the-use-of-fishes-in-research/]
• 19. The MSC Fisheries Standard Review policy development ... [https://www.sciencedirect.com/science/article/pii/S0308597X23003883]
• 20. SPIRIT 2025 statement: Updated guideline for protocols of ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12037212/]
• 21. Structure Data for Analysis [https://help.tableau.com/current/pro/desktop/en-us/data_structure_for_analysis.htm]
• 22. 11 Summarising quantitative data | Scientific Research and ... [https://bookdown.org/pkaldunn/SRM-Textbook/SummariseQuantData.html]
• 23. Understanding WCAG 2 Contrast and Color Requirements [https://webaim.org/articles/contrast/]
• 24. Text has enhanced contrast | ACT Rule | WAI [https://www.w3.org/WAI/standards-guidelines/act/rules/09o5cg/]
• 25. Understanding Success Criterion 1.4.11: Non-text Contrast [https://w3c.github.io/wcag21/understanding/non-text-contrast.html]
• 26. Evaluating PSA dynamics for predicting androgen ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC12130521/]
• 27. Examples and data for performing path similarity analysis ... [https://github.com/Becksteinlab/PSAnalysisTutorial]
• 28. Economic Viability of the Production of Peruvian Grunt ... [https://www.mdpi.com/2076-2615/15/1/48]
• 29. Chile and Peru team up to measure anchovy biomass in ... [https://bioenergytimes.com/chile-and-peru-team-up-to-measure-anchovy-biomass-in-shared-waters/]
• 30. Study finds wildlife and industrial fishing can coexist in ... [https://thefishsite.com/articles/study-finds-wildlife-and-industrial-fishing-can-coexist-in-perus-national-reserve]
• 31. Sustainable Strategies for Artisanal Fishing in the Pacific - RGSA [https://rgsa.openaccesspublications.org/rgsa/article/view/10310]
• 32. Chilean and Peruvian researchers met to define the ... [http://www.fis-net.com/fis/worldnews/worldnews.asp?monthyear=8-2024&day=15&id=131413&l=e&country=142&special=&ndb=1&df=0]
• 33. Assessing the vulnerability of Elasmobranch species in ... [https://www.sciencedirect.com/science/article/pii/S2405844024133622]
• 34. Assessing the vulnerability of Elasmobranch species in ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC11408832/]
• 35. Productivity and Susceptibility Analysis (PSA) of ... [https://www.preprints.org/manuscript/202509.2048/v1]
• 36. Fish stocks in the Bay drop nearly 80% in 7 years [https://www.thedailystar.net/news/bangladesh/news/fish-stocks-the-bay-drop-nearly-80-7-years-4047516]
• 37. Stock Assessment of Marine Elasmobranchs (Sharks and ... [https://www.mdpi.com/2077-1312/13/6/1126]
• 38. Overfishing leads to decline in Bangladesh marine fish ... [https://news.mongabay.com/2024/01/overfishing-leads-to-decline-in-bangladesh-marine-fish-stocks-diversity/]
• 39. The global issue of metal contamination in sharks, rays and ... [https://www.sciencedirect.com/science/article/pii/S0147651324014349]
• 40. Assessing the vulnerability of sensitive species in ... [https://www.frontiersin.org/journals/marine-science/articles/10.3389/fmars.2024.1411033/full]
• 41. Combining methods for detection of bycatch hotspot areas ... [https://www.sciencedirect.com/science/article/pii/S0141113625003563]
• 42. A 14-year time series of marine megafauna bycatch in the ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC8844044/]
• 43. A Case Study From the Cyprus' Pelagic Longline Fishery [https://www.frontiersin.org/journals/marine-science/articles/10.3389/fmars.2022.868464/full]
• 44. Sustainable Fisheries and Non-Target Species Management [https://www.mdpi.com/2071-1050/17/11/5040]
• 45. Large-scale distribution of a deep-sea megafauna community along ... [https://scientiamarina.revistas.csic.es/index.php/scientiamarina/article/view/1794/0?articlesBySameAuthorPage=3]
• 46. The end of progress against extreme poverty? [https://ourworldindata.org/end-progress-extreme-poverty]
• 47. Extreme poverty persists, affecting one in ten people worldwide. [https://unstats.un.org/sdgs/report/2025/goal-01/]
• 48. Fishery Data Guide [https://dlmtool.github.io/DLMtool/Data/Data.html]
• 49. Designing Effective Data Visualizations [https://guides.library.jhu.edu/datavisualization/design]
• 50. Data Visualization Tips and Best Practices [https://www.tableau.com/visualization/data-visualization-best-practices]
• 51. Text has minimum contrast[proposed] | ACT Rule | WAI [https://www.w3.org/WAI/standards-guidelines/act/rules/afw4f7/proposed/]
• 52. Color contrast - Accessibility - MDN Web Docs - Mozilla [https://developer.mozilla.org/en-US/docs/Web/Accessibility/Guides/Understanding_WCAG/Perceivable/Color_contrast]
• 53. Develop and update Productivity-Susceptibility Analyses (PSAs) of tuna fisheries in the EPO | IATTC [https://www.iattc.org/en-us/Research/Project/Detail/L-2-a]
• 54. Appendix I. Stock Assessment and Data-limited Techniques [https://wildlife.ca.gov/Conservation/Marine/MLMA/Master-Plan/Stock-Assessment]
• 55. Integrating scientific data, local knowledge, and expert ... [https://ecologyandsociety.org/vol30/iss4/art14/]
• 56. A new quantitative analysis method for the fish behavior ... [https://www.sciencedirect.com/science/article/pii/S004484862500078X?dgcid=rss_sd_all]
• 57. Experimental Studies on the Impact of the Projected Ocean ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC8614255/]
• 58. Quantitative Analysis of Fish Morphology Through ... [https://en.bio-protocol.org/en/bpdetail?id=5087&type=0]
• 59. Data weighting: An iterative process linking surveys ... [https://www.sciencedirect.com/science/article/abs/pii/S0165783623001558]
• 60. Risk identification and perception in the fisheries sector [https://www.sciencedirect.com/science/article/abs/pii/S0308597X10000989]
• 61. A Novel Global Sensitivity Measure Based on Probability ... [https://www.mdpi.com/2073-8994/13/1/90]
• 62. Robust and non asymptotic estimation of probability ... [https://arxiv.org/abs/2306.10806]
• 63. Interval reliability sensitivity analysis using Monte Carlo ... [https://www.sciencedirect.com/science/article/abs/pii/S1568494623003344]
• 64. Assessment of Seabream Fisheries Stock of Oman Using ... [https://www.mdpi.com/2071-1050/15/22/15692]
• 65. Yield-per-recruit Analyses [https://fishbase.se/MANUAL/fishbaseyieldperrecruit_analyses.htm]
• 66. Step 5: Initial Fishery Assessment | FISHE [https://fishe.edf.org/framework/step-5-initial-fishery-assessment]
• 67. Ecological risk assessment of marine resources caught as ... [https://www.frontiersin.org/journals/marine-science/articles/10.3389/fmars.2025.1490894/full]
• 68. 17 Important Data Visualization Techniques - HBS Online [https://online.hbs.edu/blog/post/data-visualization-techniques]
• 69. Data Visualizations, Charts, and Graphs - Digital Accessibility [https://accessibility.huit.harvard.edu/data-viz-charts-graphs]
• 70. A quantitative evaluation of a qualitative risk assessment ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC5983465/]
• 71. Key drivers of model choice by fisheries scientists and their ... [https://www.sciencedirect.com/science/article/pii/S0308597X24005839]
• 72. Ecological risk assessments for the effects of fishing [https://www.sciencedirect.com/science/article/abs/pii/S0165783616302272]
• 73. American Lobster Benchmark Stock Assessment Finds ... [https://asmfc.org/news/press-releases/american-lobster-benchmark-stock-assessment-finds-gom-gbk-stock-not-depleted-but-experiencing-overfishing-sne-stock-significantly-depleted-but-not-experiencing-overfishing/]
• 74. Our approach to data-limited fisheries [https://www.msc.org/en-us/for-business/fisheries/developing-world-and-small-scale-fisheries/our-approach-to-data-limited-fisheries]
• 75. Understanding Ecosystem-Based Fisheries Management [https://www.fisheries.noaa.gov/insight/understanding-ecosystem-based-fisheries-management]
• 76. A risk-based framework for assessing the cumulative ... [https://www.sciencedirect.com/science/article/abs/pii/S0006320714002274]
• 77. Ecosystem-based Fisheries Management [https://noaa-edab.github.io/presentations/20200611_PEP-EBFM_Large.html]
• 78. Color blind friendly palettes for data visualizations with ... [https://thenode.biologists.com/data-visualization-with-flying-colors/research/]
• 79. Regional scalable priorities for national biodiversity and ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC8388611/]