{'fraction_non_alphanumeric': 0.026963153937050167, 'fraction_numerical': 0.010769148136212268, 'mean_word_length': 4.821009848366422, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 1, 'https://': 2, 'lorem ipsum': 0, 'www.': 14, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Parasite induced mortality is context dependent in Atlantic salmon: insights from an individual-based model

An individual-based model was parameterized to explore the impact of a crustacean ectoparasite (sea louse, Lepeophtheirus salmonis & Caligus spp.) on migrating Atlantic salmon smolt. The model explores how environmental and intrinsic factors can modulate the effect of sea lice on survival, growth and maturation of Atlantic salmon at sea. Relative to other effects, the parasite infestation pressure from fish farms and the encounter process emerge as the most important parameters. Although small variations in parasite-induced mortality may be masked by variable environmental effects, episodes of high infestation pressure from fish farms should be observable in wild populations of Atlantic salmon if laboratory studies accurately reflect the physiological effects of sea lice. Increases in temperature in the model negatively influenced fish survival by affecting the development time of the parasite at a rate that was not compensated for by the growth of the host. Discharge from rivers was parameterized to increase migration speed and influenced parasite induced mortality by decreasing time spent in areas with increased infestation pressure. Initial size and growth of the host was inversely related to the impact of the parasite because of size-dependent parasite-induced mortality in the early phase of migration. Overall, the model illustrates how environmental factors modulate effects on the host population by impacting either the parasite load or the relative effect of the parasite. The results suggest that linking population-level effects to parasite infestation pressure across climatic and environmental gradients may be challenging without correctly accounting for these effects.

The role of parasites in regulating their host population and the role of parasite-induced mortality in population dynamics has been widely debated in parasitology and ecology 1,2 . Historically, it was thought that parasites had little effect on host survival except in rare occasions, and that individuals with high levels of parasites would ultimately succumb to factors such as predation or starvation rather than be killed by the parasites themselves. This was often attributed to the fact that fitness of the parasite is dependent on the host survival, and that the stable state equilibrium between host and parasites depends on the cost of the parasite to the host and the cost of developing different strategies to combat the parasite (e.g. immunocompetence or behavioral avoidance). In the wake of the pivotal papers by Anderson and May in 1978 1,2 , the role of parasites as a regulating factor is now recognized. However, whether mortality from parasites is additive (i.e. comes in addition to other sources of mortality) or compensatory (i.e. targets or culls weak individuals) is important to understand natural dynamics in parasite-host systems [3][4][5] .
In a natural system, the balance between hosts and parasites can vary between years, occasionally creating sporadic episodic outbreaks 6 . This balance can also be affected by humans through effects on climate 7 , introduction of invasive species that function as co-hosts 8 or removal of predators that culls diseased individuals 9 . The role of disease interaction from fish farming with surrounding marine wildlife is an anthropogenic modification that has received scrutiny recently 10 . A large proportion of studies related to fish farming has focused on the salmon and salmon lice (or sea lice in general encompassing both Lepeophtheirus salmonis & Caligus spp., Krøyer, 1837) host-parasite complex [11][12][13][14] . The role of fish farms on the epidemiology of wild fish can be summed up into three major concerns: (1) amplification of parasite load decoupled from host density 15-17 ; (2) temporal alteration of purpose. The purpose of the model is to explore how louse encounters during the nearshore migration of Atlantic salmon smolts affects their marine survival, size, and maturation. The model was designed to simulate release group experiments in which groups of smolts are released simultaneously and half of the individuals are randomly assigned to a prophylactic treatment that alleviates the fish of parasites. Specifically, the model investigates how the relationship between infestation pressure of sea louse and marine survival are affected by environmental and intrinsic factors impacting the host and parasites. Finally, the model results are compared to the risk assessment suggested by Taranger, et al. 10 and the role of marine growth and mortality in determining the fraction of overall mortality attributable to parasite effects is discussed (Fig. 1). entities, state variables and scales. The entities of the models are individual salmon. The state variables include release date, age (in days from release), length (mm), weight (gram), migration speed (km × d −1 ), variation in sea growth (unitless), infection levels (including one state variable per sea louse life stage: copepodite, process overview and scheduling. The simulated salmon were released in marine water. Each simulated smolt was assigned a random length from a normal distribution (µ = 135 mm, σ = 10) representative of wild smolt from the Vosso river in Norway 32 , and randomly assigned (50/50) to a treated or an untreated group. The migration speed of the fish was defined as a function of the body length and river discharge according to Vollset, et al. 33 . The number of days (i.e. time steps) within the nearshore environment (~100 km) was therefore dependent on the migration speed. Following each time step, the fish grew according to a temperature dependent growth function 28 . Simultaneously, the sea lice on the fish (which are tracked individually in their own IBM, Fig. 2) developed from copepodites to chalimi to preadults and finally to adults over the course of the days available for growth as a function of the temperature also experienced by its host. This means that sea lice are individual entities in their own model, but the sum of lice on each salmon can be viewed as a state variable in the salmon model. Temperature was set as a fixed number in each model run. Impacts of changing the temperature throughout the season was explored, but had relatively small impact on the outcome of the model. The likelihood of host survival was a function of host length and number of preadult and adult lice per fish gram. At the end of the simulation (415 days), the model evaluated whether the fish was mature or not based on a simple Euler-Lotka equation described in Hutchings and Jones 31 .
Design concept. Emergence. Survival (and consequently estimated mortality attributed to sea lice -i.e. attributable fraction) and maturation emerged as a result of the combined stochastic processes of louse encounter, marine growth and mortality. Sensing. The individual fish can sense temperature and encounter louse according to the infestation pressure from fish farms defined in a spatial raster. The temperature affects the fish length and the development of louse on the fish. In addition, discharge from the river impacts their migration speed through the fjord system and therefore the exposure to lice.
Interaction. There are no interactions between individuals.
Stochasticity. Initial length distribution (which is also converted to mass) is drawn from a normal distribution rendering different growth trajectories based on the temperature dependent growth model. In addition, the louse encounter and survival is stochastic in the sense that it is drawn from random negative-binomial and binomial probability functions, respectively.
Collectives. Individuals are assigned randomly to either a treated or untreated fish within each release group.
Observations. The main observations from the model are mass, survival and maturation by the end of the year, and the mortality attributed to sea lice based on the relative survival between the treated and untreated individuals. This was calculated by first running the model for the first 150 days with daily updates, and then the remaining time until river return, which was defined as July 5, in one time step (i.e. 265 days). In some instances, the mass and lice encounters of all individuals (surviving and not surviving) during the first two months was logged to illustrate how development of the lice, initial mass, growth and infestation pressure affects the host probability of mortality during the initial phase of migration.
initialization. Initially, individuals were assigned a release date, release length (total length, L T ) and a treatment group (treated / not-treated). Mass was derived from a linear mass-length relationship fitted to data from the Vosso system To increase and decrease the growth rate the second parameter was multiplied by a factor of 2 and 0.2. This changed the growth rate by approximately 10% up or down throughout the thermal range.
The effect of sea lice on the specific growth rate is not known, but Susdorf, et al. 35 suggested that a 10% reduction in growth may be valid. Also, Skilbrei, et al. 23 showed that the size of returning after one winter at sea (1SW) salmon was ~100 gram larger when a fish received treatment against sea lice as smolt. Therefore, a simple linear relationship between sea lice infestations and reduction in growth was applied as follows: in which LPG is mobile lice (preadult and adult) per gram fish. This suggests that there will be 10% reduction in sgr when an individual has 0.3 lice per gram. This is the same level of reduction in growth used by Susdorf, et al. 35 .
There is however, little data to quantify the exact impact of sea lice on growth. The resulting sgr for the impacted fish will consequently be reduction × sgr.
Migration. To model progression rate of salmon, the model from Vollset, et al. 33 was applied, which estimated daily progression rate based on release and recapture of coded-wired-tagged smolt in and outside the Osterfjord in Western Norway. A simplified version of their linear model was applied in which the migration speed (in kilometers per day) is solely a function of length of the fish and discharge from the Vosso river.
T in which v is the daily migration rate of smolts in km × d −1 , L T is the total length of the fish in mm, and D is the discharge from the river Vosso measured as − m s 3 1 . This was used to calculate the number of days spent in the nearshore environment by rounding v divided by migration distance (in km) to the nearest whole number. To (2019) 9:17377 | https://doi.org/10.1038/s41598-019-53871-2 www.nature.com/scientificreports www.nature.com/scientificreports/ increase or decrease the migration speed, the intercept was multiplied by a factor of 3 and 0.5. This changed the swimming speed by approximately 10% up or down throughout its range.
Sea lice encounter. Simulated salmon smolts encountered parasites according to the parameterized model by Kristoffersen,et al. 36 , which linked infestation pressure from fish farms to lice attached to smolts in sentinel cages during a period of 20-25 days. The infestation pressure (IP) is based on the estimated total external infestation pressure from fish farms on the Norwegian coast. The estimated IP from the model is a unitless value (on a log scale) and normally varies between 5-17. The expected number of lice encountered by a fish during one day is calculated based on the following equation were TP is defined as the transmission pressure, and was used draw a random number from a negative binomial distribution, parameterized such that mean(x) = µ and var(X) = µ + µ 2 /theta, with µ = TP and dispersion parameter theta = 2.733 (as quantified in Kristoffersen, et al. 36 ). By varying theta from 1 to 5, it was possible evaluate to what degree patchiness impacted the final outcome of the model. To make comparable model runs, the infestation pressure from week 22 in 2016 from the migration route of the Vosso salmon was extracted (Fig. 2). The data extracted were provided by the veterinary institute, and can be viewed online here: http://apps.vetinst.no/lusekart. This step was taken because the migration speed model is validated from this region 33 , and detailed knowledge about salmon migration time in this region is known from trap netting 33 . High and low infestation pressures were modelled by multiplying the infestation pressure extracted as explained above (i.e. baseline infestation pressure) by a factor of 1.2 or 0.8, yielding values comparable to high and low infestation pressure estimated from Kristoffersen, et al. 36 . Histograms resulting from modelled encounters with lice for the three infestation pressures are provided in the appendix. These values are comparable to the lice numbers observed on trawled post-smolts during sea lice surveys along the Norwegian coast 25 .
Sea lice development. Sea lice are modelled individually and each individual louse has a temperature dependent development. The developmental time from attached copepodids through chalimus to mobile (preadult and adult) developmental stages, was taken from Stien, et al. 37 and modified as explained in Groner et al. (2013). The equation used to model development is the modified Belahrádek equation, and has been used in multiple publications on salmon lice: where τ ij is the minimum required development time for salmon lice i in stage j at temperature T. β 1 is a shape parameter, and β − 2 2 is the average τ at 10 °C. As in Groner et al. (2013), a stage-specific constant, ν j , was added to this temperature-dependent estimate, to represent additional time beyond the minimum development time. Lice development from eggs to chalimus attached copepodit are integrated in the infestation pressure estimate as explained in the previous methods. Thus, parameters are only needed for the developmental time for chalimus and preadult (Table 1). In addition, each louse had a daily probability of surviving (δ) based on the parameters from Stien et al. (2005). The mortality of adult lice on pelagic swimming post-smolts is assumed to be low (~0%), based on the assumption that the mortality seen in laboratory studies is mainly due to contact with physical barriers not present in the wild. Either way, the mortality parameter of adult lice did not have a large impact on the overall outcome of the model. in which Mp is the daily natural mortality and DW is the individual dry weight of the fish. DW is assumed to be 25% of the wet weight as logged in the model (M). In an earlier version of the model, the mortality function used in Piou and Prévost 39 was applied and yielded similar results. The alternative version from McGurk 40 was also tested, but was clearly not parametrized for the size range of salmon smolts and was therefore abandoned. Lice dependent mortality was based on the threshold limits presented in Taranger, et al. 10 , fitted to a logistic regression (see appendix for details), rendering the following equation: − www.nature.com/scientificreports www.nature.com/scientificreports/ where p is the probability of survival and LPG is lice per gram fish. This was used to calculate the likelihood of surviving parasite infections (p) if LPG was >0, and set to 1 if LPG = 0. In the model p is interpreted as a daily survival rate. In reality, the values presented in Taranger, et al. 10 are not daily rates, but rather overall risk of death during the smolt migration such that p = 1-m, where m is the probability of surviving given an LPG. Finstad, et al. 25 demonstrated that mortality increased after the lice developed to preadult stages. Lice per gram is therefore number of mobile lice divided by the mass of the fish. Lice per gram and consequently mortality rate (m) therefore increases after lice develop into preadults (Fig. 3). No one has previously attempted to calculate the daily mortality rate according to sea lice (Ml). Attempts were made to extract data from publications to calculate the daily mortality rates, but the details in the publication were scarce for such estimates. For simplicity, I divided the estimate of m by five assuming that the mortality takes place during a 5 day period such that Ml = m/5. During a standard run with baseline values this would estimate an average attributable fraction of ~50% for individuals with 4-6 lice, which is comparable to the risk model suggested by Taranger, et al. 10 , and the values used by Kristoffersen, et al. 36 . The total daily probability of survival (Sp) following Eqs. 7 and 8 is as follows where, for ages x between maturity α and death ω, l x represents the survival probability from birth to the beginning of the spawning season at age x, and m x represents the number of zygotes produced by an individual spawning at age x, while r represents the population growth rate. By turning the equation it is therefore possible to calculate r for different maturation strategies, i.e. maturing after one year or multiple years at sea (1SW versus MSW). I use the parameters extracted for the Vosso salmon for this study. This entails calculating l x as the product www.nature.com/scientificreports www.nature.com/scientificreports/ of the survival from egg to smolt set to 0.15% and the survival from smolt to 1SW (which is calculated as the product of daily mortality in Eqs 8, 9 and 10), and calculating m x (i.e. fecundity) as a function of the length of the fish based on the equation where F is fecundity and L T is total length. The estimate of r maturing after 1 year at sea (r1) is then compared to the r for an individual maturing after to year at sea (r2). To be able to calculate r2, the growth and survival the second year at sea was set to fixed as 40% increase and 55% survival, which is comparable to the values used in  36 . To explore how the IBM compares with the current risk assessment, I used the estimated mortality calculated from the number of lice encountered in the infestation module of the model combined with the risk assessment model by Taranger, et al. 10 , and compared it with the estimated mortality from the IBM.

Sensitivity analysis.
A simple one-at-a-time (OAT) sensitivity analysis was conducted with the outcome of interest being mortality attributed to sea lice (Attributable fraction, AF), untreated group survival, marine growth, and maturation. All model parameters, description and their range is given in Table 1. The variables that were tested and the maximum and minimum values are given in Table 2.

Results
Examples of development of lice on fish. Figure 3 illustrates a random sample of a salmon post-smolt encountering 9 lice during the early marine migration and the temporal development of lice, the mass of the fish and the resulting lice per gram. The highest value of lice per gram is around day 45, which is almost 1.5 months after the post-smolt has left the river (Fig. 3). Temperature impacted the timing of peak virulence of sea lice on migrating post-smolt for a period of 1.5 months, which is much longer than the short period of time that is spent migrating in the coastal environment. Thus, measurements of temperature in the coastal environment will not be sufficient to understand how temperature impacts lice development on wild fish.

Sensitivity analysis.
Resulting attributable fraction, overall marine survival, final weight and proportion maturing as 1SW are given in Tables 1 and 2 for the selected parameters. In Table 3, the model is parameterized with standard baseline values. In Table 4, the model is parameterized with an assumed effect of lice on weight as explained in section "Growth" (Eq. 3). Estimates are given for baseline values of infestation levels (upper tables), in addition to high infestation levels (lower tables) to be able to evaluate the effect of the different parameters during episodes of high lice loads on the post-smolt (see section Sea lice encounter for a description of high infestation pressure). The resulting AF is also plotted in Fig. 4. The results from the sensitivity analysis suggest that the www.nature.com/scientificreports www.nature.com/scientificreports/ most sensitive parameters were the functions linking parasite load to mortality and the encounter rate function, indicating that the modelled encounter with lice and resulting mortality were important parameters. Indeed, this result has been demonstrated before by Kristoffersen, et al. 36 , which showed that changing the mortality function   Table 3. Sensitivity analysis in which the model is parameterized with no effect of sea lice on growth rate on their salmon host. The values for the parameters at minimum and maximum values are given in Table 2. AF signifies attributable fraction, SC survival in untreated group, FW final weight and PM proportion of population matured after one winter at sea. www.nature.com/scientificreports www.nature.com/scientificreports/ and the link between infestation pressure and modelled lice had a large effect on the estimated mortality. Other important parameters were temperature, initial size of the fish, the growth and swimming functions and the discharge. Patchiness had a small effect on the outcome of the model. The natural mortality function (size-selective mortality) had a strong effect on the overall survival and maturation of the individuals, but had seemingly no measurable effects on the effect of parasites on survival of salmon (see appendix for the three parameterizations of size selective mortality). Finally, comparing the two tables it is evident that a parameterized effect of lice on growth generally increases the impact of lice on survival and maturation.

AF (%) SC (%) FW PM AF (%) SC FW PM
temperature. The effect of parasites increase with temperature (Fig. 5a). This is a result of the temperature dependent growth and size dependent mortality. The untreated group impacted by parasites was also impacted by temperature through the temperature dependent developmental time of the parasite. This yielded an increased impact of lice at high temperature and high lice infestation levels. The faster growth rate of salmon smolts cannot compensate for the decreased developmental time for the parasite. This is also illustrated in the overall survival of the untreated group, which shows that the survival increases with temperature in groups experiencing low or baseline infestation levels, but decreases with temperature when the model is parameterized with a high infestation pressure (Fig. 5b). Parameterizing an effect of lice on growth increased the attributable fraction due to parasite mortality, but decreased the effect of temperature (i.e. decreased the slope between temperature and AF). The ultimate size of the returning salmon increased with temperature, but there were almost no observable differences between groups from low or high infestation pressure in size of surviving individuals (Fig. 6a). Proportion maturing as 1SW increased with temperature in low and baseline infestation levels. For high infestation pressure, proportion of 1SW did not increase, but was low throughout the temperature range (Fig. 6b). In addition, parameterizing an effect of lice on growth decreased the proportion maturing as 1SW slightly.
Size of released fish. The effect of parasites decreased with increasing size of fish (Fig. 5c). This was particularly evident in groups with high infestation pressure, but was only evident in baseline infestation pressure settings when the model included an effect of lice on growth. This effect was relatively small in a natural range of wild salmonids (120-150 mm) but was evidently large when compared with sizes that are sometimes observed  Table 3. The values are also given in Tables 3 and 4 www.nature.com/scientificreports www.nature.com/scientificreports/ among hatchery reared fish (200 mm). Similar to temperature, initial size was important for growth (Fig. 5, mid left) and age at maturity (Fig. 6d).
Discharge and swimming speed. Increased discharge yielded faster progression speed through the fjords, reducing the days in areas with increased infestation pressure from fish farms. The effect was relatively small compared to the effects of temperature and initial size (Fig. 5e), and was very similar to increasing the swimming speed of the fish. However, the effect should be observable if studies are conducted in extreme conditions of spring floods versus periods of very low discharge. Discharge had very little observable effect on growth (Fig. 6,e) and maturation (Fig. 6f).
Growth. Increasing growth rate of the salmon smolt (for example in a scenario were prey were abundant) had an opposite effect to increasing growth of salmon smolts by increasing temperature. This occurs because increased growth allows the fish to outgrow the negative impact of sea lice. This was the same effects observed when increasing initial size. This was more important when infestation pressure was low compared to when the infestation pressure was high, because the relative effect of size was more important when lice loads was low.
comparison model to risk model. The resulting risk estimated mortality estimates from the low, baseline and high infestation pressure using the model from Kristoffersen, et al. 36 was 0.012%, 1.5% and 46.5%. Results from the OAT model suggests that the risk model generally underestimates the mortality compared to the IBM

Discussion
Manipulative studies have demonstrated that parasites such as sea lice affect marine survival 13,22,23,43 . However, linking infestation pressure from fish farms or parasite abundance observed on fish in situ to measurable effects of the parasites on fish populations has been challenging. For example, in an attempt to model the impact of sea lice from fish farms on wild fish, Vollset, et al. 13 used an infestation pressure model described in Jansen, et al. 44 , but could not account for heterogeneity in the effect of antiparasitic treatment on release groups of salmon. Instead they found that the heterogeneity in treatment was found to be strongly linked to the baseline survival of the release group. In other words, they found a positive effect of antiparasitic treatment on survival of salmon when marine survival was relatively low, but a smaller effect when marine survival was high. According to the authors the effect of the parasite is therefore likely modulated by other factors operating on the survival of salmon. Accounting for all possible factors in the wild that may modulate parasite induced mortality is extremely difficult; simulation studies assist in bridging these knowledge gaps and bring context to observations of fish in the wild. The results from the model presented in this study demonstrate that several environmental factors can impact the effect of parasites populations of salmon and can be summarized into three main novel findings: (1) temperature mediated growth indicates that effects of sea lice should be stronger with higher marine water temperature. This can have implications both for the effect of parasite-induced mortality as the season progress, as well as the Figure 6. Size of salmon after 1 year at sea (415 days) and proportion maturating after 1 year at sea from the individual-based model parameterized with low infestation pressure levels (circles and solid line),with baseline infestation pressure levels (triangles and dashed line), and high infestation pressure (crosses and dotted line). Each value are plotted against temperature (a,b), initial size (c,d), and discharge (e,f). Grey and black lines indicate models parameterized with or without an effect of sea lice on the growth rate of the salmon host.
impacts of different populations along latitudinal gradients and in the face of climate change; (2) initial size or growth can modulate the effect of sea lice-induced mortality; and (3) conditions that decrease the time spent in the nearshore areas (such as high river discharge or swimming speed) will reduce impacts of sea lice on salmon.
The model suggested that temperature can have an impact on the severity of parasite-induced mortality. This impact can be mainly attributed to the fact that the host (salmon) and the parasite (sea lice) respond very differently to temperature 28,37 . Compared to the shorter-lived sea lice 37 salmon develop at a much slower rate 28 . This means that at an equal sea lice infestation, lice reach the mobile stage, which is known to potentially do damage to the host 25 , at a much larger parasite:host size ratio than at colder temperatures. Warming sea temperatures are therefore predicted to shift the balance in the parasite-host complex of sea lice and salmonids. If this prediction holds, it should also be expected that sea lice is more harmful for salmon populations at the southern boundaries of their distribution. Another important aspect is that parasite distribution among their hosts is usually strongly skewed such that a few individuals host many parasites, whereas a large proportion of potential hosts has small levels of infestation. It is therefore important to understand the physiological effects of warming waters on parasite-host interactions when abundance of parasites are low in order to make ecologically accurate predictions which are valid on a population level.
Size at release was also an important parameter in the model, indicating that the effect of lice increase curvilinearly with decreasing size. In tagging and cultivation programs, size is often used to compensate for the lower returns of cultivated fish compared to wild fish 45 . Consequently, usually large fish are used in studies releasing treated and untreated salmon smolt 13 . The model suggests that future trials should attempt to specifically test the interactive effect of release size and parasite-induced mortality, something that is revealed to be important, but lacking in empirical studies. This effect can also be relevant for the understanding of the impacts of wild fish as size-at-sea-entry may vary between populations, latitude, and among years.
Similar to initial size, decreased growth rate was related to increased parasite-induced mortality. In this model, growth was purely modulated by temperature. In reality, however, growth in wild salmonids is strongly impacted by their encounter with suitable prey. However, as of now no model exists that can correctly model the post-smolt growth of Atlantic salmon, although analysis of scales has given us insights into annual and seasonal trends in growth and how it may vary with climate 46 . The general trend in the North Atlantic has been declining growth and survival in southern populations, which has been attributed to large scale ecosystem regime shifts 47 . In relation to parasite-induced mortality, variable marine growth, coupled to poor survival, may be the underlying mechanism behind the strong correlation observed between survival and parasite-induced mortality 13 . From a salmonid ecology point of view, we are still scraping the surface of the mechanisms that govern the marine survival and growth of salmon.
Discharge and migration length both had an intermediate effect on the impact of sea lice. This emerges in the model because of the longer exposure time due to lowered progression rates through the fjord. If fish migrate during peak flow, as shown in some acoustic telemetry studies 48 , it will strongly reduce the effect of sea lice. Peak flows are projected to be more contracted with climate change because of flashier river dynamics associated with less snowfall and more rapid melting. Mismatched timing between spring freshets and productive ocean conditions could imperil smolt runs as they either enter the ocean early when conditions are poor or late when lice pose a greater threat. Discharge may, in addition to affecting the progression rate of smolts through estuaries and fjords, affect the salinity that the smolt encounter during migration. Sea lice are known to avoid low salinity water 49 , and this effect may increase the importance of river discharge during smolt migration. Although this model does not account for this affect, the relationships among river discharge, estuary salinity, sea lice abundance, and smolt survival should be explored in future studies.
The encounter rate function had a relatively strong effect on the simulation outcome. This was also demonstrated in Kristoffersen, et al. 36 . Similarly, the effect of the function that relates lice per gram fish to mortality was important. Together with the infestation pressure, these variables seemed to be the most sensitive parameters in the model, reinforcing that efforts must foremost minimize uncertainty around these parameters. Various institutions in Norway and internationally have attempted to model the distribution of sea lice larvae using hydrodynamic models to be able to predict encounter rates with sea lice in space and time as post smolts embark on their marine migration [50][51][52] . Modelling encounter rates are key ecological concepts, and various model frameworks 53,54 exist that can be applied to the study of how salmon and sea lice interact. Less work has been put into repeating studies that quantify the level of sea lice that is lethal for a post-smolt. This is clearly a pivotal area of study that needs more attention. One method of collecting data on of infestation levels, proposed by Vollset et al. (2017), is experimental infestation of individual fish, tagging them (e.g. with passive integrated transponders) and releasing them into the wild. When these are then recaptured as adults the effect of varying degree of infestation levels can be assessed.
Even though sea lice was parameterized to affect the growth of the salmon, this had no effect on the ultimate size of the returning after one winter at sea. This occurred because reduction in size was also linked to decreased likelihood of survival. This could explain why the effect of sea lice on the size of returning salmon is generally small 23 . If the model predictions are correct, it is perhaps not surprising that we do not see a larger difference between treated and untreated fish when they return as adults. In contrast, the effect on age at maturation was strongly impacted by sea lice. This was mainly because the maturation equation takes into consideration the likelihood of survival and growth. Consequently, the effect on age at maturation is even stronger when sea lice were parameterized to affect the growth of the salmon. In this aspect, the model seems to deviate from the release groups studies. For example, Vollset, et al. 24 demonstrated a significant, albeit relatively small, effect of treatment on proportion maturing after one year at sea. This could either reflect that the model exaggerates the effect of sea lice on age at maturation, or that the release group studies indeed reflect that these groups in general experience a generally low infestation pressure. (2019) 9:17377 | https://doi.org/10.1038/s41598-019-53871-2 www.nature.com/scientificreports www.nature.com/scientificreports/ The impact of sea lice on the marine survival of salmon is one of the most debated topics in salmonid ecology, because of the complexity in acquiring unambiguous data to document the effect (if any). Marine survival of Atlantic salmon is variable and complex to study and most likely consists of several interacting factors acting together in time and space. This complexity is not uncommon when studying marine survival (or recruitment). However, when implementing management systems that link a management action (in this case keeping lice levels on fish farms low) to a specific survival of a target species, evidence-based decision making is crucial. In comparison with the risk evaluation suggested by Taranger, et al. 10 , and applied to modelled post-smolts in Kristoffersen,et al. 36 , the model predicts a larger effect when infestation pressure is low, while similar predictions (on average) are obtained when infestation pressure is high. Perhaps the biggest difference between this model and the risk evaluation of salmon used in various publications is that (1) the individuals are subject to size selective mortality during and after the process of encountering lice, and (2) that the host and parasite development and growth are affected by temperature. This simulation has yielded important results including previously undetected temperature-mediated risks associated with lice parasitism that could alter the balance between parasite and host in the wild.

conclusions
The model predicts a variety of patterns that can be tested: (1) Temperature mediated growth indicates that effects of sea lice should be stronger with warming water. This can have implications for the effect of parasite-induced mortality during the season, with latitudinal differences. (2) Initial size or growth can modulate the effect of sea lice-induced mortality. (3) Conditions that decrease the time spent in the nearshore areas (such as high river discharge or swimming speed) will reduce the encounter with sea lice.