Comparing Life Cycle Assessment ( LCA ) of Salmonid Aquaculture Production Systems : Status and Perspectives

Aquaculture is the fastest growing food sector worldwide, mostly driven by a steadily increasing protein demand. In response to growing ecological concerns, life cycle assessment (LCA) emerged as a key environmental tool to measure the impacts of various production systems, including aquaculture. In this review, we focused on farmed salmonids to perform an in-depth analysis, investigating methodologies and comparing results of LCA studies of this finfish family in relation to species and production technologies. Identifying the environmental strengths and weaknesses of salmonid production technologies is central to ensure that industrial actors and policymakers make informed choices to take the production of this important marine livestock to a more sustainable path. Three critical aspects of salmonid LCAs were studied based on 24 articles and reports: (1) Methodological application, (2) construction of inventories, and (3) comparison of production technologies across studies. Our first assessment provides an overview and compares important methodological choices. The second analysis maps the main foreground and background data sources, as well as the state of process inclusion and exclusion. In the third section, a first attempt to compare life cycle impact assessment (LCIA) and feed conversion ratio (FCR) data across production technologies was conducted using a single factor statistical protocol. Overall, findings suggested a lack of methodological completeness and reporting in the literature and demonstrated that inventories suffered from incomplete description and partial disclosure. Our attempt to compare LCA results across studies was challenging due to confounding factors and poor data availability, but useful as a first step in highlighting the importance of production technology for salmonids. In groups where the data was robust enough for statistical comparison, both differences and mean equalities were identified, allowing ranking of technology clusters based on their average scores. We statistically demonstrated that sea-based systems outperform land-based technology in terms of energy demand and that sea-based systems have a generally higher FCR than land-based ones. Cross-study analytics also strongly suggest that open systems generate on average more eutrophying emissions than closed designs. We further discuss how to overcome bottlenecks currently hampering such LCA meta-analysis. Arguments are made in favor of further developing cross-study LCA analysis, particularly by increasing the number of salmonid LCA available (to improve sample sizes) and by reforming in-depth LCA practices to enable full reproducibility and greater access to inventory data. Sustainability 2019, 11, 2517; doi:10.3390/su11092517 www.mdpi.com/journal/sustainability Sustainability 2019, 11, 2517 2 of 27


Introduction
Population growth, dietary shifts and resource challenges in capture fisheries and agriculture are driving the development of aquaculture production systems worldwide [1].By 2050, the global population is expected to reach 9.8 billion [2], and robust economic growth in developing nations is expected to lead to a dietary shift towards higher consumption of meat and dairy [3].These commodities are known to carry a high environmental burden [4].While one can hardly argue that social development is negative, population growth and higher affluence are expected to generate additional pressure on fragile ecosystems already severely threatened by a variety of adverse human impacts [5].Current intensive and linear production systems of food, feed, fiber, and bioenergy are largely exceeding the earth's capacity to regenerate the consumed biomass, exposing our societies to serious environmental and food security risks in the long run [6].Global food production cannot continue to depend on agriculture and capture fisheries to the same degree in the future.Fisheries have been stagnating since the 1990s, constrained by the natural production limits of wild populations [7].Meanwhile, agricultural systems are facing various challenges including decreasing soil fertility and arable land availability [8,9].
Aquaculture is one of the most efficient food production systems the world will rely upon to meet tomorrows demand for animal-sourced foods and reduce current environmental impacts associated with the supply of animal protein [1,10].Since the 1980s, one of the largest developments in aquaculture has been driven by carnivorous salmonids, especially Atlantic salmon (Salmo salar), rainbow trout (Oncorhynchus mykiss), and coho salmon (Oncorhynchus kisutch).Salmonid aquaculture is central to the global seafood landscape since it alone accounts for 18.1% of the value and 7.4% of the biomass of the world trade in fish and fish products [1].This production is especially important in Northern Europe (Norway, Ireland, Scotland, Faroe Islands, Iceland) and North and South America (Chile andCanada), in countries with abundant access to cold/temperate marine waters and extensive coastlines.Salmonids are primarily consumed in the United States, Europe, and East Asia.While consumption growth slows down in these traditional markets (but remain positive), consumption is now increasing the most in emerging markets such as Southeast Asia and Latin America [11].Salmonid species are particularly appreciated for their taste, reddish-orange colored fillets, high-quality proteins, and their marine omega-3 fatty acids content [12,13].Farmed salmon are also considered a sustainable alternative to terrestrial meat producing species, including beef cattle, hogs, and broiler chickens [14,15].Salmonids are efficient livestock with low energy expenditure due to their cold-blooded homeostasis physiology and the gravitational constraints opposed by buoyancy under water [16,17].Yet, despite these favorable characteristics, salmonid aquaculture faces a variety of environmental challenges.While aquaculture production addresses some of the issues of current food systems by using efficient marine livestock and (partly) moving production at sea, it does not address the underlying challenges of intensive and linear systems (inherited from agricultural practices).Future production increases will hinge on the industry's ability to solve challenges posed by parasites, diseases, nutrient emissions, and critical feed resource scarcity [18][19][20], among others.For instance, in Norway, the aquaculture industry has ambitions to increase salmonid production fivefold by 2050; however, it is increasingly clear that such a target cannot be achieved without identifying and addressing the underlying environmental issues of existing open sea-based systems [21].
Life cycle assessment (LCA) is one of the environmental accounting tools that can be used to provide the critical information necessary to improve the sustainability of such aquaculture systems.During the last 20 years, the emergence of LCA as a standardized method to measure the environmental impacts of products and services enabled researchers to assess complex seafood systems and provide sustainability guidelines to the industry and policymakers [22][23][24].Since salmonids are one of the most high-value types of aquaculture species produced worldwide [1], and the interest to increase production volumes is high, this finfish type is the most studied species group in the aquaculture LCA literature.In fact, it is the only group for which LCA studies exist for various species, countries, feed recipes, and production technologies.LCA has been employed to compare products with similar functions [25], measure differences between feed formulations [26,27], and to evaluate the environmental efficiency of trout [28][29][30], and salmon [31][32][33].Within this species group, the LCA methodology has also been employed repeatedly to compare the environmental impact of existing and emerging production technologies [34][35][36].
Following the accumulation of aquaculture LCA studies, several reviews recently analyzed this section of the literature, focusing on the application of the LCA methodology [37][38][39] and comparing production systems [40][41][42][43].In the latter ones, attempts to compare results focused primarily on the comparison of the system's intensities (intensive, semi-intensive, and extensive) and culture types (polyculture vs. monoculture).Overall, the authors strictly avoided cross-study single factor statistical comparisons, because LCA data are known to be subject to a wide range of intricate influential variables.The scoping of this review has a sensitive balance, since a broad scope increases sample sizes and the influence of confounding factors, while a narrow one reduces discrepancies and the number of data points available.While previous reviews in this field focused on the whole aquaculture sector, we believe farmed salmonids to be a suitable species group to do an in-depth, comprehensive literature review.Such a scope allows us to study how different production systems for similar species with similar requirements compare, and if such LCA comparisons are meaningful.We consequently analyzed three key features of current salmonid LCAs.Firstly, we provide an overview and compare methodological choices.The objective is to evaluate common methodological practices, identify shortcomings, and suggest recommendations to practitioners.Secondly, we map and compare the inventory and modeling protocol of each study.We strive to identify the strengths and weaknesses of models and datasets used in publications.Thirdly, we perform a cross-study analysis on the life cycle impact assessment (LCIA) and FCR results through the prism of production technologies.Within this third analysis, we are making a first attempt to perform a cross-study statistical comparison of LCIA and FCR scores.This statistical protocol voluntarily ignores confounding factors to specifically investigate if LCA results comparisons can be performed on relatively homogenous samples (similar salmonid species) despite the influence of other variables.With this, our objectives are to assess if overarching statistical differences between groups can be identified, if discrepancies between studies are simply too great to validate such comparisons, and if LCA methodological developments are required to enable more direct comparability.

Study Selection
A systematic analysis of 24 LCA studies (20 peer-reviewed articles and four gray literature research) based on salmonid aquaculture systems was performed in this review (Table 1).The systems assessed produced Atlantic salmon (Salmo salar), rainbow trout (Oncorhynchus mykiss), Arctic char (Salvelinus alpinus), Chinook salmon (Oncorhynchus tshawytscha), brook trout (Salvelinus fontinalis), and brown trout (Salmo trutta fario).All of these species are classified under the Salmonidae family.Article identification was performed using the online search engines of Science Direct, Web of Science, Google Scholar, and Google.The keywords used were a combination of "LCA", "life cycle assessment", "aquaculture", "farmed salmon", "trout", "rainbow trout", and "marine products."The timespan ranged from 2000 to 2018.We identified 355 studies of interest.Primary refinement focused on titles and abstracts, with "environmental assessment" and "aquaculture system" as the selection criteria to isolate 65 publications.Secondary refinement concentrated on methods and models to select research focusing on salmonid aquaculture systems and containing original LCA case studies, resulting in the final selection of 24 assessments (Section 2.1.1,"article selection", supplementary data).A study was deemed original when significant first-hand data collection of at least one of the main foreground processes (hatchery and grow-out) took place.Five studies qualifying for selection were excluded.The inventories of Buchspies et al. [44], Ytrestøyl et al. [45], and Hall et al. [46] revealed an absence of original data.Winther et al. [47] was excluded in favor of Ziegler et al. [48] since both are presenting the same data/results and the latter is a peer-reviewed publication.Samuel-Fitwi et al. [49] was removed from the list in favor of Samuel-Fitwi et al. [50] since both studies were built with the same dataset and the latter assessed different production technologies (Section 2.1.2,"study exclusion", supplementary data).

Analytical Procedure
This work was articulated around three main analytical sections: (1) The specific choices related to LCA methodology, (2) the construction of life cycle inventories, and (3) a cross-study technological comparison based on compiled LCIA and FCR results (Section 2.2.1, "scope of the review", supplementary data).All analyses were performed in Excel.The first section presented the main methodological choices made by LCA practitioners.This encompassed mapping and comparing selections of functional units (FU), system boundaries, characterization methods, characterization factors, multi-functionality, as well as checking for the inclusion of interpretation tools such as contribution, sensitivity, and uncertainty analyses.The second section measured and compared the data quality of inventories by reviewing the main foreground and background sources used by authors.We also estimated study completeness by mapping process inclusion and exclusion, including an evaluation of the LCA practitioner reporting practices.In this analysis, the literature was divided into two groups: Trout and salmon systems.This separation was performed to highlight significant inventory differences between the trout and salmon LCAs, and its application was strictly limited to this section.Overall, 14 independent processes were identified to fully cover the life cycle of salmonids (Table 2).
In the third section, LCIA and key inventory data in terms of FCR results were compared across studies according to their production technology.Although FCR scores are LCI data and as such differ significantly in nature from LCIA results, we found it meaningful to add this parameter to the analysis.The FCR indicates the quantity of feed use per unit of fish produced, and although it has been criticized lately in cross-species efficiency comparison [61], we estimated this parameter to be representative of the overall efficiency of salmonid production systems.It is also a parameter directly linked to the aquafeed process, which is by far the largest contributing component to the life cycle emissions of salmonids [51, 54,57].FCR has recently been identified as a useful indicator to follow up environmental performance in farmed aquaculture within a species and using the same feed [62].While it was not possible to standardize feed compositions across studies due to a lack of inventory disclosure, our comparison dealt with species of the same family, with similar nutrient requirements to enhance comparability of FCRs.For the LCIA data, the comparison was performed for the four most commonly used impact categories in the reviewed studies: Global warming potential (GWP), acidification potential (AP), eutrophication potential (EP), and cumulative energy demand (CED).Out of the 24 LCA studies selected, a total of 67 scenarios were identified.These LCA scenarios have variable focus, among others investigating impacts of different feed composition, electricity mix, or cage structure.Scenarios were selected based on their production technologies.Scenarios outside the scope of this review were excluded.When scenarios used similar production technology and only the cofounding factor differed, results were averaged.After selection, 35 scenarios remained for GWP, 30 for AP, 29 for EP, and 24 for CED.Of those scenarios, 33 disclosed their FCR scores, which were subsequently used in this analysis (Section 2.2.2, "LCIA scenario selection", supplementary data).While the comparability of LCIA results is uncertain due to the wide range of diverging parameters existing between studies (e.g., year, country, production site, feed composition, characterization method, assumption, etc.), we increased consistency by only comparing cradle-to-gate systems producing 1 t live-weight salmonids.Most studies already expressed their results in this format.Six studies used different system boundary and FU combinations [35,36,48,52,55,58].We performed a small system boundary adjustment of Ayer et al. [52] by adding hatchery impacts using contribution analysis results from Pelletier et al. [31].Cradle-to-farm gate impacts of 1 t FU produced were available in Liu et al. [35] and Ziegler et al. [48] and part of the secondary results.Newton and Little [55], Silvenius et al. [58], and Grönroos et al. [36] adjusted impacts were obtained through personal communication with the authors.Grönroos et al. [36] adjusted results differ significantly from the published results due to an inventory error discovered by the authors (Section 2.2.3, "adjustment calculations", supplementary data).
All scenarios were separated into four production technology clusters.Salmonid aquaculture production technologies revolve around two major characteristics: Sea-based vs. land-based (clusters A-B) and open vs. closed (clusters C-D) systems.This led us to the selection of four-production technology clusters: (1) Closed sea-based systems, (2) open land-based systems, (3) open sea-based systems, and (4) closed land-based systems.Cluster (1) contains experimental pilots of either a solid wall aquaculture system [33,56] or marine floating bag [34,36].These systems provide a controlled production environment at sea and a form of waste collection.Cluster (2) exclusively gathers flow-through aquaculture system (FTAS) technology.Such raceways usually divert flowing water (often from rivers using gravity) through successive concrete flumes before discharging back to the water body [63].While this is mostly a freshwater system used for rainbow trout production, a few studies also investigated the production of salmon and char using marine and brackish waters [27,34].Cluster (3) is the most common form of salmon production, usually taking place in net pens of various sizes and forms (mostly circular).In open cages, only nets separate the biomass from the environment, which means that salmonid metabolites and extra feed pellets are directly emitted to the surrounding waters [64].Cluster (4) is mostly recirculating aquaculture system (RAS) technology.This land-based design uses circular tanks to produce trout, char, and salmon, primarily in freshwater conditions.Water movement and oxygenation is insured by machinery and is recycled using drum and bio-filters [65].

Statistical Analysis
A statistical analysis of each production clusters for LCIA and FCR scores was performed in SPSS v24 (IBM Corp, Armonk, New York, United States).In this review, we deliberately opted for a single parametric statistical protocol to compare the means of the different technological clusters.This is a first attempt to make a statistical LCA cross-study comparison without involving multifactorial analysis to investigate if a simpler, non-discriminant, comparative protocol can be applied to identify overarching trends in LCA results.We are not trying to compare the groups independent of confounding factors, but to investigate if current LCA data can be statistically compared despite them.The equivalent non-parametric test to ANOVA (Kruskal-Wallis) was disqualified due to heterogeneous distributions in the data.The LCIA and FCR data was analyzed following a four-group comparison (clusters 1 to 4), and two subsequent two-group comparison testing data of clusters A-B and C-D (Figure 1).Throughout the analysis, the significance level α (referring to the probability of making a type one error) was set to 5%.An identical test protocol was applied to all data.Each group was tested for normality with Shapiro-Wilk and homogeneity of variances using a Levene's test.For Shapiro-Wilk, p-values were tested against the null hypothesis, assuming a normal distribution within clusters, while for Levene's test, the null hypothesis was supporting similar group variances.If homogeneity of variance was achieved, a one-way ANOVA t-test was used to compare the means of clusters.Where variances were found to be statistically different, a Welch t-test was performed to compare the homogeneity of the means.If the Welch-test indicated a statistically significant difference between some of the means, the Games-Howell post hoc test was performed to identify which of the cluster's means were different from one another.Results of both the ANOVA and Welch t-test were rejected if the group involved displayed were skewed.For these three statistical tests, the null hypotheses assume means equality.

Methodology Analysis
The location of the production systems reflects the current dominance of salmonid aquaculture by a handful of industrialized countries; Norway, the United Kingdom, Canada, Chile, and Australia dominate salmon production, while Finland and France concentrate on trout species.Dealing with multi-functionality remains controversial and heterogeneous, with practitioners using a mix of mass, economic, and energy-based allocation throughout the literature (Table 3).System expansion, recommended by ISO 14040 [66], is rarely employed.All studies follow the attributional LCA framework, except for Samuel-Fitwi et al. [50], employing the consequential approach.Out of the 24 studies reviewed, 20 are using different versions of SimaPro (PRé Sustainability, Amersfoort, the Netherlands) for modeling and calculations.CML-IA remains the most common characterization method (applied 14 times), followed by the ReCiPe midpoint framework (applied five times).All studies provide a process-based contribution analysis, and most of them include a sensitivity analysis to test for critical data sources and/or methodological choices.Only four studies [30,33,52,54] included an uncertainty analysis.The range of midpoint impact categories selected is widespread across studies.GWP, EP, AP, and CED are by far the most represented categories with application rates of respectively 100%, 87%, 78%, and 61%.Characterization factors used to calculate these midpoint scores evolve over time, but are mostly (in frequency) based on the work of IPCC [67], Heijungs et al. [68], Huijbregts [69], Frischknecht [70], and Frischknecht et al. [71].LCA practitioners often add biotic resource use (43%), water dependence (43%), and land occupation (39%) to measure additional impacts significant for salmonid aquaculture production.Results show strong discrepancies between the characterization methods, the characterization factors and their application by authors.For example, in CML-IA, the EP characterization factors (based on Heijungs et al. [68]) refer to phosphorus (kg PO4-e) without discerning marine and freshwater nutrient limitation, while ReCiPe expresses the same impacts based on nitrogen (kg N-e) and phosphorus (kg P-e) depending on the emissions medium [72].In different studies, water dependence is also named water depletion [30], water use [29,59], freshwater footprint [32], and consumptive water use [55].Land occupation is inconsistent in both naming and units between studies.

Methodology Analysis
The location of the production systems reflects the current dominance of salmonid aquaculture by a handful of industrialized countries; Norway, the United Kingdom, Canada, Chile, and Australia dominate salmon production, while Finland and France concentrate on trout species.Dealing with multi-functionality remains controversial and heterogeneous, with practitioners using a mix of mass, economic, and energy-based allocation throughout the literature (Table 3).System expansion, recommended by ISO 14040 [66], is rarely employed.All studies follow the attributional LCA framework, except for Samuel-Fitwi et al. [50], employing the consequential approach.Out of the 24 studies reviewed, 20 are using different versions of SimaPro (PRé Sustainability, Amersfoort, the Netherlands) for modeling and calculations.CML-IA remains the most common characterization method (applied 14 times), followed by the ReCiPe midpoint framework (applied five times).All studies provide a process-based contribution analysis, and most of them include a sensitivity analysis to test for critical data sources and/or methodological choices.Only four studies [30,33,52,54] included an uncertainty analysis.
The range of midpoint impact categories selected is widespread across studies.GWP, EP, AP, and CED are by far the most represented categories with application rates of respectively 100%, 87%, 78%, and 61%.Characterization factors used to calculate these midpoint scores evolve over time, but are mostly (in frequency) based on the work of IPCC [67], Heijungs et al. [68], Huijbregts [69], Frischknecht [70], and Frischknecht et al. [71].LCA practitioners often add biotic resource use (43%), water dependence (43%), and land occupation (39%) to measure additional impacts significant for salmonid aquaculture production.Results show strong discrepancies between the characterization methods, the characterization factors and their application by authors.For example, in CML-IA, the EP characterization factors (based on Heijungs et al. [68]) refer to phosphorus (kg PO4-e) without discerning marine and freshwater nutrient limitation, while ReCiPe expresses the same impacts based on nitrogen (kg N-e) and phosphorus (kg P-e) depending on the emissions medium [72].In different studies, water dependence is also named water depletion [30], water use [29,59], freshwater footprint [32], and consumptive water use [55].Land occupation is inconsistent in both naming and units between studies.It is alternatively referred to as agricultural land occupation [30], surface use [53], land use [54,55], and land competition [50,60], and is expressed in square meters per year [28][29][30]50,54,55,60] or square-meters [53] alternatively.Finally, biotic resource use is also called net primary production in about one-fifth of the studies [28,29,51,53,59] (see supplementary data, Section 3.1.1).
A vast majority of trout studies exclude egg production (broodstock reproduction), effluent treatment, processing, and distribution.Out of 10 studies reviewed, the inclusion frequency of these processes was 30%, 40%, 20%, and 10% (Figure 2).Reasons for these exclusions vary from the absence of reliable technology and legislation (e.g., effluent treatment) to the selection of narrower system boundaries (e.g., processing, distribution), and the estimated negligible impacts of a given process (e.g., egg production).None of the trout LCAs included the infrastructure end-of-life (EOL) process because of its versatility and estimated overall low environmental contribution (cut-off criterion).The inclusion of equipment and infrastructure and chemotherapeutants are relatively high, with a use frequency of 70%.Equipment and infrastructure construction of trout systems is intensive since this species is mostly farmed with FTAS and RAS land-based technology.In the literature, the use of chemicals and veterinary products is often reported under the overarching term "chemotherapeutant".However, in practice, this process has an uncertain connotation.In some studies, it refers to chemicals, disinfectants, limestone, and oxygen used in the hatchery and grow-out processes, while in others it also includes veterinary treatments such as antibiotics.Overall, energy carrier, transport, feed production, hatchery, fish production, and effluent are the most common processes (Figure 2).Practitioners expect these processes to have the highest contribution, which makes them less prone to cut-offs, system boundaries, and FU exclusions (see supplementary data, Sections 3.2.1.1 and 3.2.1.2).
[52] to specifically compare two types of infrastructure through their life cycle (nylon and copper net-pens).Effluent treatments are systematically excluded from open net-pen cages studies but are described by Ayer and Tyedmers [34] and McGrath et al. [33] in the modeling of closed sea-based systems equipped with filters (see supplementary data, section 3.2.2.1 and section 3.2.2.2).In the life cycle inventories of salmon LCA models, energy carrier, transport, feed production, hatchery, fish production, and effluent processes have a very high inclusion rate (Figure 3).Only hatchery and effluent are under the 100% inclusion with 86% and 79%, respectively.These lower rates can be explained by few process exclusions due to narrower system boundaries (e.g., hatchery) and/or to impact categories selection (e.g., effluent).As such, only carbon footprint studies [32,35,48] excluded effluent since it is estimated not to affect GWP.Chemotherapeutant, equipment, infrastructure, egg production, effluent treatment, infrastructure EOL, processing, and distribution have lower inclusion rates ranging from 7% to 50%.Here again, processing and distribution directly depend on the study's scope and its system boundaries.Infrastructure EOL is included in Ayer et al. [52] to specifically compare two types of infrastructure through their life cycle (nylon and copper net-pens).Effluent treatments are systematically excluded from open net-pen cages studies but are described by Ayer and Tyedmers [34] and McGrath et al. [33] in the modeling of closed sea-based systems equipped with filters (see supplementary data, Sections 3.2.2.1 and 3.2.2.2).

Overview of LCIA and FCR Results
GWP impacts vary widely across the reviewed literature with a delta between the maximum (13,622 kg CO2-e) [53] and the minimum (1157 kg CO2-e) [54] reaching 12,465 kg CO2-e (Figure 4).Within clusters 1 to 4, the situation is also heterogeneous with delta scores of respectively 1585; 4252; 7170; and 11,799 kg CO2-e.Overall, closed sea-based systems perform best with an average score of 2404 kg CO2-e per ton LW salmonids produced at the farm gate.Open land-based systems rank second with an average impact of 2613 kg CO2-e while open sea-based and closed land-based systems arrive in 3 rd and 4 th position with 2933 and as much as 6414 kg CO2-e for the same production output.FCR values seem to fluctuate independently of climate change impacts (Figure 4).
AP impacts are also heterogeneous.Grönroos et al. [36] reported as little as 8.8 kg SO2-e per ton LW trout while Ayer and Tyedmers [34] calculated a score as high as 63.4 kg SO2-e for a similar

Overview of LCIA and FCR Results
GWP impacts vary widely across the reviewed literature with a delta between the maximum (13,622 kg CO 2 -e) [53] and the minimum (1157 kg CO 2 -e) [54] reaching 12,465 kg CO 2 -e (Figure 4).Within clusters 1 to 4, the situation is also heterogeneous with delta scores of respectively 1585; 4252; 7170; and 11,799 kg CO 2 -e.Overall, closed sea-based systems perform best with an average score of 2404 kg CO 2 -e per ton LW salmonids produced at the farm gate.Open land-based systems rank second with an average impact of 2613 kg CO 2 -e while open sea-based and closed land-based systems arrive in 3rd and 4th position with 2933 and as much as 6414 kg CO 2 -e for the same production output.FCR values seem to fluctuate independently of climate change impacts (Figure 4).
AP impacts are also heterogeneous.Grönroos et al. [36] reported as little as 8.8 kg SO 2 -e per ton LW trout while Ayer and Tyedmers [34] calculated a score as high as 63.4 kg SO 2 -e for a similar output (Figure 5).Intra-cluster heterogeneity is comparable to that of climate change impacts.Deltas of 10.1 kg SO 2 -e, 23.3 kg SO 2 -e, 33.6 kg SO 2 -e, and 51.5 kg SO 2 -e are observed between minimums and maximums of clusters 1, 2, 3, and 4 respectively.Closed land-based systems perform worst with 26.7 kg SO 2 -e, while closed sea-based systems outperform all other clusters by generating 15.1 kg SO 2 -e per FU on average.Open land-based and sea-based systems reach the 2nd and 3rd positions with average scores of 16.3 kg SO 2 -e and 18.7 kg SO 2 -e.Graphic representation of the FCR seems to fluctuate independently of AP impacts scores (Figure 5).For instance, Samuel-Fitwi et al. [50] calculated emissions of 40.7 kg SO2-e for an FCR of 0.86 while Dekamin et al. [54] 6).Intra-cluster deltas remain high with differences ranging from 30.3 kg PO 4 -e, 51.2 kg PO 4 -e, and 62.6 kg PO 4 -e for systems 4, 2, and 3 respectively.Closed sea-based systems are the exception, displaying a low delta of 8.6 PO 4 -e for this impact category.Any direct correlation between FCR and EP impacts remain widely uncertain with studies displaying both low impacts and high FCR [34,54] as well as low FCR and high impacts [55,57].6).Intra-cluster deltas remain high with differences ranging from 30.3 kg PO4-e, 51.2 kg PO4-e, and 62.6 kg PO4-e for systems 4, 2, and 3 respectively.Closed sea-based systems are the exception, displaying a low delta of 8.6 PO4-e for this impact category.Any direct correlation between FCR and EP impacts remain widely uncertain with studies displaying both low impacts and high FCR [34,54] as well as low FCR and high impacts [55,57].Variability of CED scores is substantial throughout the reviewed literature.The delta for this category represents 206,100 MJ-e, which corresponds to a difference of factor 9.23 between the maximum of 233,000 MJ-e [34] and the minimum of 26,900 MJ-e [34] (Figure 7).Intra   Variability of CED scores is substantial throughout the reviewed literature.The delta for this category represents 206,100 MJ-e, which corresponds to a difference of factor 9.23 between the maximum of 233,000 MJ-e [34] and the minimum of 26,900 MJ-e [34] (Figure 7).Intra Variability of CED scores is substantial throughout the reviewed literature.The delta for this category represents 206,100 MJ-e, which corresponds to a difference of factor 9.23 between the maximum of 233,000 MJ-e [34] and the minimum of 26,900 MJ-e [34] (Figure 7).Intra-cluster impact deltas range from 172,569 to 29,764 MJ-e with closed, land-based systems reporting the most variable scores (indicating improvement potential).Closed sea-based systems displayed the most consistent results.Salmonid life cycle CED results show that land-based systems have the highest energy requirements with cluster averages equal to 133,220 MJ-e (closed) and 75,943 MJ-e (open).Sea-based systems, on the other hand, perform best with average CEDs of 37,913 MJ-e (open) and 54,620 MJ-e (closed) per ton LW salmonids produced at the farm gate (Figure 7).Again, a high FCR does not appear to correlate with high CED directly.Studies reporting low FCR and high CED [27,60] and some displaying the inverse relationship [30,59] [27,60] and some displaying the inverse relationship [30,59] are commonly found.4).Since FCR averages variations between impact categories are minimal and only imputable to sample size, the sub sequential statistical four-groups and two-groups analysis of FCR averages were performed independently of impact categories variations, on the single largest sample size available (also corresponding to the FCR GWP group).4).Since FCR averages variations between impact categories are minimal and only imputable to sample size, the sub sequential statistical four-groups and two-groups analysis of FCR averages were performed independently of impact categories variations, on the single largest sample size available (also corresponding to the FCR GWP group).Out of the 16 LCIA groups (Figure 1), cluster 2-3 for GWP and 2-4 for AP failed normality tests.The remaining clusters displayed normal distributions.Levene's results showed that all groups violated the null hypothesis, pointing to heterogeneous variances between groups.GWP and AP Welch t-test means comparison were rejected because of the skewed distributions of some of their clusters.This means that 50% of the data is not comparable according to this statistical protocol due to non-normal distributions.The means comparisons were validated for EP and CED.The test revealed statistically significant differences of the means between the clusters of EP (p = 0.004) but no differences were found between the clusters' means of CED (p = 0.066).Post hoc Games-Howell test identified the statistical difference within EP to be between clusters 1-2, 3-4, and 2-4.It means that overall, only significant differences were identified for EP between closed

(b) Two-Groups Analysis
Two two-group analyses were performed on both LCIA and FCR data.One was testing the data of sea-based (cluster A) vs. land-based (cluster B) studies while the other was testing open (cluster C) vs. closed (cluster D) systems.Each two-group LCIA analysis consists of eight clusters while each FCR one consists of two (Figure 1).Cluster A-B for GWP, B for AP, and A for EP failed the normality test.Homogeneity of variances was observed in AP and EP while the null hypothesis was rejected for GWP and CED.GWP, AP, and EP ANOVA and Welch means comparisons were rejected due to skewness in at least one of their clusters.This means that 75% of the data could not be tested and interpreted.The Welch t-test mean comparison for CED was accepted, showing significant differences between the means of closed, sea-based and open, land-based clusters for this impact category (p = 0.01).All clusters failed normality, except cluster D for EP in the open vs. closed data comparison.Levene's test demonstrated that clusters of GWP, EP, and CED had heterogenous variances while AP clusters displayed a homogeneous one.Since seven out of the eight clusters were skewed, results of both ANOVA and Welch t-tests were rejected.In this case, 100% of the data was rejected and could not be compared using this statistical protocol (see supplementary data, Sections 3.3.2.3 and 3.3.2.4).The FCR data of clusters A-B confirm the null hypothesis of Shapiro-Wilk, supporting that both samples are distributed normally.The same results could be observed for the FCR data of clusters C and D. Also, both A-B (sea-based and land-based) groups and C-D (open and closed) systems displayed heterogeneous variances.In the A-B comparison, the null hypothesis of Welch was rejected (p = 0.021), demonstrating a significant difference between the means of sea-based and land-based clusters.Opposite results were found for the C-D comparison (p = 0.654); open and closed systems FCR means were not statistically different from each other (Section 3.3.2.5, "FCR two-group stat.A-B", and Section 3.3.2.6, "FCR two-group stat.C-D", supplementary data).

Methodologies
The methodological variations identified in this LCA review are understandable considering the difficulty of developing a "one size fits all" LCA methodology for production systems that differ in time, space, and structure.Overall, the results showed relatively consistent FU and system boundaries but displayed variations in the way LCA practitioners deal with multi-functionality and LCIA.These variations were not necessarily negative since methods selections were often case-specific, determined by the scope of the study, and could also be driven by methodological developments.Keeping this flexibility enables practitioners to adapt the method to the production system instead of the opposite.However, method customization should not be achieved at the expense of clarity.For instance, for impact category selection, it is essential to cover the main environmental impacts of the system while also limiting the number of categories to keep the results comprehensible for readers.This issue is discussed in Steinmann et al. [80], who found that within current impact categories, seven indicators were estimated to cover 92% of the variance of product rankings from the Ecoinvent database (based on 976 products).It means that by selecting GWP, land use, ozone depletion, AP, EP, marine, and terrestrial ecotoxicity, the practitioner covers nearly the entire spectrum of impacts at the midpoint level.The authors also suggest using RACER (relevant, accepted, credible, easy, robust) criteria to refine the indicator selection, as well as encourage the use of regionalized impact categories [80].Based on Steinmann's perspective and the categories selection rates in salmonid LCA, we recommend practitioners assessing such aquaculture systems to concentrate primarily on: (1) GWP, (2) EP, (3) AP, (4) CED, (5) marine ecotoxicity, (6) terrestrial ecotoxicity, and (7) land-use.Biotic resource use and water use are also of interest, depending on the scope of the study.However, while these established impact categories are adapted to measure impacts of product and services, they are not necessarily well equipped for production technology comparison.The development of aquaculture specific impact categories is long wished for in the literature [40,50,54,58] and suggestions to include local ecological impacts have already been made years ago [81].In this perspective, accounting for the impacts of aquaculture on biodiversity and water quality suggested by Ford et al. [81] would significantly improve the environmental criteria by which we compare the different types of salmonid aquaculture production systems.
There is an overall lack of innovative methodological integration in LCA of salmonids production systems (Table 3).For instance, it is unclear why several recent publications [27,54,57] continue to apply characterization factors from CML 2001 [68,82,83] when new (and more comprehensive) indicators are available in ReCiPe [72,84,85].Lack of comprehensive and standardized methodological reporting is another shortcoming of the current literature.For example, since no rules exist for including and reporting contribution, sensitivity, and uncertainty analyses, the quantity and quality of information provided by LCA practitioners vary widely.While contribution analysis is commonly disclosed in LCA studies (although not always in numerical terms), the inclusion of uncertainty analysis remains rare.In this review, only one out of five LCIA results were calculated with an uncertainty range (Table 3).Systematic uncertainty assessments would benefit the comparisons conducted in this paper.The lack of overall methodological description is a challenge for replicability, transparency, and review.For instance, gathering the detailed characterization factors used for each impact category was complicated by the overall lack of data description (both in articles and supplementary data).We recommend that LCA practitioners pursue a flexible application of the LCA methodology to best suit each study's scope.It should entail considering the latest methodological LCA developments (when relevant) and going further than the LCA guidelines provided by the ISO framework [86] using, for instance, the procedure developed by the European Research Commission [87].However, this adaptive application of LCA necessitates detailed and systematic reporting to be fully exploitable (reproducible and comparable) by the research community.It is therefore critical that researchers document their LCA methodological choices extensively, especially in the supplementary data.We recommend practitioners to systematically disclose and argue for major modeling decisions (e.g., characterization factors) and assumptions (e.g., cut-off criteria) performed under their assessment.

Inventories
In both trout and salmon LCA studies, process completeness is directly influenced by the FU and system boundaries selection (e.g., hatchery, processing, distribution), by the choice of impact categories (e.g., effluent), by the type of technology studied (e.g., effluent treatment), and by the scope of the research (e.g., infrastructure EOL).The inclusion of equipment and infrastructure is lower in salmon LCI, which could be explained by the "light-weight" of open sea-based cages compared to RAS and FTAS systems, more frequently used in trout production.The inclusion of chemotherapeutic agents appears more frequently for trout than salmon (55% and 23%, respectively).This could be due to the importance of water quality treatments (input of oxygen, limestone, and other chemicals) in FTAS and RAS production systems.Often the exclusion of processes such as chemotherapeutants, equipment, infrastructure, and egg production is based on the assumption that those processes have an overall low contribution [35,54,55].Since process-based LCA remains resource intensive to conduct, practitioners often choose to perform such cut-offs, especially if data availability is low and/or quality poor.Interpreting the importance of such processes (when accounted for) is also challenging due to their variable definitions across studies.For instance, several authors include the production of salmonid eggs in their inventory [27,28,33,50], but the specific nature of this process remains uncertain because detailed descriptions are lacking.Considering the growing biological challenges faced by open sea-based salmon production [88,89], accounting for the environmental impacts of parasites and diseases treatments is becoming a necessity.So far only three of the 13 salmon LCAs reviewed accounted for a form of chemotherapeutant use (antibiotics and water quality treatments), and only for closed land-based production phases (hatcheries using FTAS or RAS systems) [32,34,57].The frequent salmon lice sea-based treatments using biological (cleaner fish), chemical (baths and feed), and mechanical (primarily warm-water baths) processes remain unaccounted for in current LCAs.In Norway, use of cleaner fish requires aquaculture and fisheries activities dedicated to produce and capture Lumpfish (Cyclopterus lumpus) and various species of Wrasses (Labridae).The value-chains involved are complex and require significant use of marine ingredient rich aquafeed to produce the cleaner fish and sustain them throughout the salmon production cycle.Overall, the use of chemicals is following a descending trend, mainly because of sea-lice resistances.Yet, the use of hydrogen peroxide and other substances like azamethiphos and cypermethrin in baths remain commonly used, generating unquantified marine ecotoxicity impacts.In addition, the use of both mechanical and chemical delousing treatments has proven deleterious effects on farmed salmon mortality [90].It is essential to determine how much direct and indirect impacts are linked to these various types of delousing treatments, especially when we know that in recent years, mortality represented between 6-9% of the total Norwegian production in terms of biomass [91].
Overall, inventories also suffer from limitations similar to those from the methodologies: Incomplete description and partial disclosure.While data confidentiality can explain why portions of the inventory cannot be fully disclosed, we still think process description and data reporting can be largely improved.Good research practices could consist to systematically disclose foreground systems, including the sources, time, and place of data collection.More descriptions of the production systems are also required.It would increase transparency and enable more subsequent analyses based on this data.Pelletier et al. [31] and White [59] are two good examples of how inventory data can efficiently be presented.When confidentiality is a requirement from companies, authors should work closely with the data provider to share a maximum of data while still safeguarding sensitive information.Aggregating inventory inputs can be one way to do so.
Another issue identified in this review concerns data representativity.Although most studies report collecting first-hand data for the most crucial processes (feed and grow-out), a significant fraction of the datasets are built with recycled data from previous studies.For example, the feed dataset developed in Pelletier and Tyedmers [92] and Pelletier et al. [31] was reportedly used in six subsequent salmonid LCAs, primarily because of its completeness (high representativity), availability (supplementary data), and perceived quality (highly cited paper).A few years ago, Henriksson et al. [93] suggested a protocol for the horizontal averaging of foreground and background data, accounting for data uncertainty generated by inaccurate measurement and means variability.This protocol is one of the tools LCA practitioners can use to circumvent cross-referencing of obsolete data.There is also a lack of open data on salmonid feed ingredients and the production processes of aquafeed, especially since it is estimated central for salmonid LCA to use data representative of the feed used.Substantial inclusion of an analog dataset tends to artificially increase the level of diversity and completeness (more case studies tested with similar data) of a given research area and reduces variability between studies.The construction of national inventories and open LCA databases is crucial to address these issues, and it has already been recommended in an aquaculture LCA review published so far [37,39].We, therefore, recommend LCA practitioners to support the development of open inventories such as AGRYBALYSE [94], the Seafood LCI database [95], or the Thai LCI database [96].We also suggest researchers to participate more actively in the GFLI (Global Feed LCA Institute) initiative and collaborate with the FEFAC (European Feed Manufacturers' Federation) to access to the latest data.
It can also be highlighted that the extensive use of case studies in salmonid LCAs (and the broader field of LCA) are not always well equipped to fully capture the environmental impacts of large production value-chains such as aquaculture.In salmonid LCA, data collection is limited in time (a production cycle) and often limited in space (one production site).Exceptions consisting of Chen et al. [28], Papatryphon et al. [24], Pelletier et al. [31], Ziegler et al. [48], and White [59] who covered significant proportions of the national productions they analyzed (global production in the case of Pelletier et al. [31]).Process-based LCA is well tailored to compare specific production parameters, sites, technology, or products, but is not optimized to achieve process completeness or to model whole sectors of an economy.Using environmental input-output models could address these issues by providing process completeness and resolution to national aquaculture sectors.Research to develop successful hybrid input-output/LCA models have been taking place for several years [97,98] and the recent development of EXIOBASE v3 [99] and its incorporations in a hybridized framework is promising [100,101].In practice, this could be done by replacing conventional LCA background processes by data from the multi-regional supply and use input-output databases.Not only could such frameworks account for side processes usually left out of process-based LCA (e.g., office activities, bank transactions, labor), but it could also allow practitioners to integrate processes with low inclusion levels such as chemotherapeutants, equipment, infrastructure, egg production, as well as the whole upper value-chain (processing and distribution).Unsurprisingly, this method will also bring its load of challenges.Among others, LCA practitioners will likely have more difficulties to disaggregate processes, work at local or regional scales, and choose restrictive system boundaries.

Outlying Results
Large intra-cluster differences observed throughout the four impact categories are due to the variable nature of production systems studied and of LCA methodology applied.Although we review similar species using standardized system boundaries and FU, LCAs are performed on specific systems, using different feed regimes, varying in time and space, and calculated with adaptive LCA methodology.This inevitably generates significant discrepancies.Yet, if we analyze the outlying results for each impact category, individual studies can be identified, and their extreme scores can be linked to methodological choices, specific inventory, or other particularities.Salmon impacts from Ayer and Tyedmers [34], White [59], and Parker [57] display consistently high outstanding impacts while d'Orbcastel et al. [53], Grönroos et al. [36], Nyhus [56], and Dekamin et al. [54] report much lower impacts for trout production.Samuel-Fitwi et al. [50] generated both fringe high and low scores depending on the impact category considered.In Ayer and Tyedmers [34] the inventory shows that both the FTAS and RAS scenario report elevated electricity consumption (13,400 kWh and 22,600 kWh, respectively).This high demand coupled to the Canadian electricity mix (61% hydro, 18% coal, 13% nuclear, 4% oil, and 4% natural gas) leads to unusually high impacts for GWP and AP (Figure 4; Figure 5).Parker [57] and White [59] studied the environmental impacts of Tasmanian salmon, commonly fed with high inclusion of animal by-products from fisheries, poultry, as well as mammals like swine and beef cattle.In these studies, the higher impacts from the feed composition are accentuated by high FCR levels (1.3 and 1.46).These factors combined explain the remarkably high GWP, AP, EP, and CED scores of these two studies.The borderline high results obtained by Samuel-Fitwi et al. [50] in GWP and AP are directly linked to the high electricity demand of the system (19,622 kWh) and the high percentage of fossil fuels in the electricity mix of the region.In fact, by replacing this mix by wind power, the RAS GWP impacts are cut by a factor of 10 [50].The reasons for the outlaying low scores for EP and CED remain uncertain.The authors report nearly zero nutrient emissions thanks to the mechanical and biological filtration systems, yet such filters are not a unique component of this study.The significant differences between attributional and consequential LCA frameworks could explain some of these differences.Rainbow trout LCA studies appear to generate most of the extreme low values across both technology clusters and impact categories (Figures 4-7).d'Orbcastel et al. [53] reported consistently low impacts, especially for the hypothetical RAS scenario displaying a low FCR and high fish density throughout production.The hypothetical nature of this production system and the reuse of fish slurry as biofertilizer could explain these outlier results.While the corrected data provided by the author [102] somewhat leveled the differences between Grönroos et al. [36] and other studies, this assessment still reports some of the lowest scores for the closed sea-based and open sea-based clusters AP.The significant difference with the other studies remains uncertain, but we suspect the lower AP scores to be linked to outdated data (1987 to 2000), to the feed composition, which was largely based on fish meal at this time, and potentially the use of Finnish characterization factors.Similarly to d'Orbcastel et al. [53], Dekamin et al. [54] report low GWP and AP scores for the FTAS scenario and a low EP for the RAS system.In the latter, water outlet is treated by drum and bio-filters as well as ion exchange, which could explain nutrient emissions as low as 8.5 PO4-e.Based on the contribution analysis, the low GWP and AP are due to an advantageous FCR and low electricity requirement [54].

Cross-Study Statistical Comparison
Although a significant portion of the LCIA statistical analysis was rejected, the current data state allowed us to perform some statistical comparisons of clusters means on both LCIA and FCR data.In the LCIA four-group analysis we statistically demonstrated that the clusters mean were similar for CED, which means that there was no statistical evidence that CED of salmonids production systems was varying in pair with the production technologies.Yet, by increasing the sample size of the comparison in the two-group LCIA A-B analysis we demonstrated that sea-based and land-based production technologies correlated with CED impacts.This means that on average, sea-based outperformed land-based systems for this impact category (41,768 vs. 96,771 MJ-e).Such results have been repeatedly observed in the literature, where land-based systems generally demonstrated to have higher energy consumption due to the land-based water circulation requirements [34,36,51,60].There is also evidence that clusters 1-  6).While there were strong indications that EP impacts correlate with the open and closed characteristics of productions systems, this was not strictly demonstrated since EP impacts of clusters 1-3 have equal means and the results of the two-group LCIA C-D analysis were rejected due to skewness.These results were also coherent with the scientific literature since open systems have no waste collection, and therefore emit more nutrients into the environment than their closed counterparts [42].
All the FCR data qualified for statistical analysis.This suggests that FCR data is less sensitive to some confounding factors since it is an inventory parameter and not a multivariable computed LCIA result.In the four-group FCR analysis, we demonstrated that cluster 2-3 have statistically different means while all the other cluster comparisons did not.It means that on average, open land-based outperform open sea-based systems (1.124 vs. 1.256) for the FCR.The technological characteristic determinant for this difference is identified in the two-group FCR A-B and C-D analyses, once the sample size has been expanded.The A-B comparison shows that the sea-based or land-based traits of the systems were responsible for this difference while the C-D analysis proved that the open or closed characteristic did not influence the results.Overall, it means that land-based systems have lower FCR than sea-based systems on average (1.124 vs. 1.257).We hypothesize that this significant difference is due to the more controlled environment unique to land-based facilities.At sea, salmonids are exposed to fluctuating seasonal temperatures, parasites, and latent or active disease outbreaks.These dynamic factors invariably lead to less predictable feeding, and mortality discrepancies, which directly affects the FCR of sea-systems negatively.

Data Quality for Meta-Analysis
One of the clear outcomes of the statistical analysis performed on LCIA and FCR data, is the bottleneck generated by the poor meta-data quality available from LCAs of salmonids.While it certainly exists other statistical protocols one could have used to compare the means of different clusters in various groups, the statistical tests we applied were selected for their robustness, simplicity and recognized comparative properties across scientific fields [103][104][105].We observed that 50% of the LCIA four-group, 75% of the LCIA two-group A-B, and 100% of the LCIA two-group C-D analyses were rejected due to the skewness of the data.Small sample sizes, data skewness, variability of distributions, and the multitude of confounding factors existing between studies represent a major limitation for such statistical meta-analysis.Increasing the scope of the comparison may improve sample sizes and distributions but it will inevitably reinforce influences of confounding factors.To improve the potential of cross-study results comparison in LCA, we identify two main levers.The first one consists of significantly increasing the quantities of LCA studies available per commodities to increase sample sizes and statistical power required for more robust statistical comparisons.Today, the inertia of LCA studies are such that for certain commodities, only a few case-studies are published sporadically through time.Process-based LCA analyses are expensive to conduct and often require several months of data collection, often largely depending on the voluntary participation of companies.This bottleneck could be partially addressed using hybridized IO-LCA methodology (see 4.2."Inventories") to model more efficiently background systems.It could also be tackled by means of direct access to data streams from production systems, reducing drastically time and resources applied to collect data manually from companies.The quantity of data collected by companies is increasing exponentially, and the adoption of industry 4.0 (automated manufacturing using connected sensors) could drastically increase the speed and ease of access to data [106].Such a technological revolution could also allow practitioners to perform LCA on digital twins of production systems to enable various forms of prospective environmental assessments.The second lever should provide the tools to enable full reproducibility.This implies re-thinking in-depth the current LCA methodology.Reaching full reproducibility would require major changes in LCA practices.Today, LCA reproducibility is drastically limited by data ownership, poor data availability, limited disclosure practices (see 4.1."Methodologies" and 4.2."Inventories"), and restrictive software licenses.To alleviate confounding factors limiting cross-study comparisons without using a strictly standardized LCA framework for all systems, practitioners need access to the full inventory in the LCA software.With direct access to the LCA models, it would be much more feasible for researchers to verify, reproduce, and modify LCAs for comparative purposes.This off-course raises issues since desegregated foreground data is usually subject to confidentiality agreements with companies, background data is subject to access to LCA databases (e.g., Ecoinvent), and most popular LCA software (e.g., SimaPro, Gabi) are licensed.A way around these challenges could consist to provide incentives for open data and open software (e.g., OpenLCA) and favorize the publications of LCAs in open access journals.Monetary data is freely open in most countries and represent one development lead.Open LCA databases are valuable initiative requiring more development (see 4.2."Inventories").Finally, a third alternative to democratize open data could be achieved by national manufacturing reporting schemes managed by governmental entities.For instance, in Norway, the aquaculture industry is required to report some production data (e.g., biomass produced, mortality, feed use, etc.) to the Directorate of Fisheries regularly.Such data is compiled and made available for stakeholders.With the emergence of industry 4.0, efficient data collection, and data transfer, one could imagine that authorities could collect detailed production data from companies to dynamically monitor the environmental impacts of the commodities they are producing.Emerging technologies such as distributed ledger technology could also be an integral part of such systems to certify that the data provided by the manufacturer has not been modified or tampered with [107,108].

Conclusions
Throughout our analyses, the idea that LCA practitioners should strive to achieve systematic reporting rather than looking for a single and standard way to apply LCA on salmonid aquaculture systems was strengthened.Using balanced levels of methodological flexibility gives researchers the freedom they need to adapt the methodology to their production system and specific goals, thereby allowing methodological innovation without compromising minimum requirements of comparability.Systematic reporting of methodological choices and inventory data is an essential comparability requirement currently lacking in the reviewed literature.We, therefore, urge LCA researchers to disclose their methodology and their foreground system inventory in their supplementary data.Such practices will facilitate future attempts to conduct reviews, but more importantly, it would allow the scientific community to access, analyze, adapt, compare, and reuse data for other purposes.Journal editors are in a position to make sure authors comply by enforcing full data and methods transparency as part of their publication requirements.
The lack of process and methodology completeness is another pitfall we identified across the literature.While sensitivity and uncertainty analyses are not required to perform an ISO-compliant LCA, these are essential research tools to measure uncertainty surrounding data, methodological choices or assumptions, and to identify the critical environmental hot spots of a system.Still, often authors omitted sensitivity and uncertainty analyses and excluded processes on uncertain basis.Although we cannot expect LCA practitioners to reach 100% completeness using process-based LCA, there is room for improvements.For instance, we expect current disease and parasite treatments performed on salmonid open sea-based systems, activities that to date have been excluded, to give significant contributions to total impacts.
Our first attempt to do a cross-study analysis based on salmonid LCA demonstrates two main results.Firstly, although the data quality remains weak, it proves possible to perform single factor comparisons across studies despite confounding factors.We were able to compare some of the data statistically and, based on differences measured, rank specific clusters according to their average EP, CED, and FCR scores.Secondly, we demonstrated that such meta-analyses are limited by the current state of LCA methods and data.To improve salmonid LCA cross-study comparison, researchers should thrive on increasing the quantity and reproducibility of LCA data and on making sure that LCA methodological developments lead towards more data availability for stakeholders.We believe that investigating the meta-analysis potential of LCA is critical for methodological progress, and more importantly, to drive sustainable development of products and services.Influencing the sustainability of production sectors through research is a challenging task for which isolated, highly specific environmental assessments are not necessarily well suited for.Industrial actors and policymakers need robust research-based information, built on overarching trends, to be able to make sustainable strategic choices for the future.

Figure 1 .
Figure 1.Structure of the life cycle impact assessment (LCIA) and feed conversion ratio (FCR) statistical analysis displaying the numbers of data points per clusters.

Figure 1 .
Figure 1.Structure of the life cycle impact assessment (LCIA) and feed conversion ratio (FCR) statistical analysis displaying the numbers of data points per clusters.

Figure 2 .
Figure 2. Process inclusion and exclusion in trout LCA studies (data collection dates are placed in brackets).

Figure 2 .
Figure 2. Process inclusion and exclusion in trout LCA studies (data collection dates are placed in brackets).

Figure 3 .
Figure 3. Process inclusion and exclusion in salmon LCA studies.

Figure 3 .
Figure 3. Process inclusion and exclusion in salmon LCA studies.

Figure 7 .
Figure 7. Salmonids cumulative energy demand (CED) impacts (MJ-e) and FCR based on production technology clusters.The average FCR values of clusters 1, 2, 3, and 4 are 1.262, 1.124, 1.256, and 1.125, respectively for GWP.While FCR fluctuates slightly across impact categories due to differences in sample sizes, results consistently show that open and closed land-based clusters have the lowest average FCR with scores ranging between 1.124-1.131and 1.113-1.131and that closed sea-based and land-based technologies oscillate in the upper range within intervals of 1.225-1.265and 1.249-1.291(Table4).Since FCR averages variations between impact categories are minimal and only imputable to sample size, the sub sequential statistical four-groups and two-groups analysis of FCR averages were performed independently of impact categories variations, on the single largest sample size available (also corresponding to the FCR GWP group).

Figure 7 .
Figure 7. Salmonids cumulative energy demand (CED) impacts (MJ-e) and FCR based on production technology clusters.The average FCR values of clusters 1, 2, 3, and 4 are 1.262, 1.124, 1.256, and 1.125, respectively for GWP.While FCR fluctuates slightly across impact categories due to differences in sample sizes, results consistently show that open and closed land-based clusters have the lowest average FCR with scores ranging between 1.124-1.131and 1.113-1.131and that closed sea-based and land-based technologies oscillate in the upper range within intervals of 1.225-1.265and 1.249-1.291(Table4).Since FCR averages variations between impact categories are minimal and only imputable to sample size, the sub sequential statistical four-groups and two-groups analysis of FCR averages were performed independently of impact categories variations, on the single largest sample size available (also corresponding to the FCR GWP group).
, sea-based and open, land-based systems (p = 0.030), between open, sea-based and closed, land-based systems (p = 0.012), as well as open, land-based and closed, land-based systems (p = 0.006).For the four FCR groups (Figure 1), each cluster validated the condition of normality and results showed clear heterogeneity of variances (p = 0.0001).The null-hypothesis of Welch was rejected (p = 0.023) indicating a significant difference between at least two of the clusters' means.Post hoc testing indicated that the difference was lying between cluster 2-3 (p = 0.009; open, land-based and open, sea-based systems) while the remaining of the clusters had statistically equal means (see supplementary data, Sections 3.3.2.1 and 3.3.2.2).

Table 1 .
List of articles and reports selected for review.French Research Institute for the Exploitation of the Sea; e UMA = University of Mohaghegh Ardabili; f SINTEF = The Company for Industrial and Technical Research; g SYKE = Finnish Environmental Institute; h NTNU = Norwegian University of Science and Technology; i GMA = Society for Marine Aquaculture; j LUKE = Natural Resource Institute Finland; a INRA = National Institute for Agricultural Research; b EarthShift Global = sustainability consultancy company; c Jacques Witford = environmental consulting engineering company; d IFREMER = k Matis = Icelandic food and biotechnology R&D institute; l SIK = Swedish Institute for Food and Biotechnology; and U = University.* Non-peer reviewed reports.

Table 2 .
List of salmon and trout life cycle processes.

Table 3 .
Main methodological choices used in salmonids life cycle assessment (LCA).
Sustainability 2019, 11, x FOR PEER REVIEW 10 of 27 be central for several aquafeed models [31,33,34,52,56,59].Like for trout, Thrane [76], Schau and Fet [74], and Tyedmers [77] are sporadically used in Ellingsen and Aanondsen [25], Pelletier et al. [31], and Ziegler et al. [48] to model capture fisheries.In recent work, Parker [57] used Cashion et al. [78] reported an FCR of 1.47 with emissions as low as 18.7 kg SO2-e.Out of 29 EP scores extracted from LCA scenarios, values ranged from 4 kg [50] to 84 [57] kg of PO 4 -e, generating a difference of factor 21. Results suggested that closed systems perform best, emitting only 17.3 kg PO 4 -e (land-based) and 26.7 (sea-based) kg PO 4 -e.Open systems generated on average 47.3 kg PO 4 -e (sea-based) and 50.6 (land-based) kg PO 4 -e into fresh and marine water bodies (Figure Sustainability 2019, 11, x FOR PEER REVIEW 12 of 27 output (Figure5).Intra-cluster heterogeneity is comparable to that of climate change impacts.Deltas of 10.1 kg SO2-e, 23.3 kg SO2-e, 33.6 kg SO2-e, and 51.5 kg SO2-e are observed between minimums and maximums of clusters 1, 2, 3, and 4 respectively.Closed land-based systems perform worst with 26.7 kg SO2-e, while closed sea-based systems outperform all other clusters by generating 15.1 kg SO2-e per FU on average.Open land-based and sea-based systems reach the 2 nd and 3 rd positions with average scores of 16.3 kg SO2-e and 18.7 kg SO2-e.Graphic representation of the FCR seems to fluctuate independently of AP impacts scores (Figure5).
are commonly found.

Table 4 .
Overview of clusters average FCR for each impact categories.

Table 4 .
Overview of clusters average FCR for each impact categories.