Mercury concentrations in fish in Finland: Difference between revisions
Jump to navigation
Jump to search
(→Calculations: code for conc_mehg moved from Goherr assessment) |
(→Calculations: assumption changed: normal -> lognormal) |
||
| Line 83: | Line 83: | ||
<rcode name="conc_mehg" label="Initiate ovariable conc_mehg (for developers only)" embed=1> | <rcode name="conc_mehg" label="Initiate ovariable conc_mehg (for developers only)" embed=1> | ||
# This is code Op_en4004/conc_mehg on page [[Mercury concentrations in fish in Finland]] | |||
library(OpasnetUtils) | library(OpasnetUtils) | ||
| Line 95: | Line 96: | ||
test <- paste(Hg$Parameter, Hg$Naytelaji, Hg$Helcom.alue) | test <- paste(Hg$Parameter, Hg$Naytelaji, Hg$Helcom.alue) | ||
meanSalm <- Hg$Result[test == " | meanSalm <- Hg$Result[test == "MeanLog lohi "] | ||
meanHerrN <- Hg$Result[test == " | meanHerrN <- Hg$Result[test == "MeanLog silakka C-Selkämeri"] | ||
meanHerrS <- Hg$Result[test == " | meanHerrS <- Hg$Result[test == "MeanLog silakka H-Pohjoinen Itämeri"] | ||
sdSalm <- Hg$Result[test == " | sdSalm <- Hg$Result[test == "SDLog lohi "] | ||
sdHerrN <- Hg$Result[test == " | sdHerrN <- Hg$Result[test == "SDLog silakka C-Selkämeri"] | ||
sdHerrS <- Hg$Result[test == " | sdHerrS <- Hg$Result[test == "SDLog silakka H-Pohjoinen Itämeri"] | ||
N <- openv$N | N <- openv$N | ||
| Line 109: | Line 110: | ||
Iter = rep(1:N, 2), | Iter = rep(1:N, 2), | ||
Country = rep(c("FI", "SWE", "DK", "EST"), each = 2*N), | Country = rep(c("FI", "SWE", "DK", "EST"), each = 2*N), | ||
Result = c( | Result = 10^c( | ||
rnorm(N, meanHerrN, sdHerrN), | rnorm(N, meanHerrN, sdHerrN), | ||
rnorm(N, meanSalm, sdSalm), | rnorm(N, meanSalm, sdSalm), | ||
| Line 161: | Line 162: | ||
Hg <- Hg[!is.na(Hg$Conc),] | Hg <- Hg[!is.na(Hg$Conc),] | ||
Hg$Year <- as.numeric(as.character(substr(Hg$Alkupäivä, 7,8))) | Hg$Year <- as.numeric(as.character(substr(Hg$Alkupäivä, 7,8))) | ||
Hg$Century <- ifelse(Hg$Year < 18, 2000, 1900) | |||
Hg$Year <- Hg$Century + Hg$Year | |||
colnames(Hg)[4] <- "Naytelaji" # remove ä | colnames(Hg)[4] <- "Naytelaji" # remove ä | ||
Hg[Hg$Naytelaji == "lohi" , "Helcom-alue"] <- "" | Hg[Hg$Naytelaji == "lohi" , "Helcom-alue"] <- "" | ||
oprint(aggregate(Hg["Conc"], by = Hg[c("Helcom-alue", "Naytelaji")], FUN=length)) | oprint(aggregate(Hg["Conc"], by = Hg[c("Helcom-alue", "Naytelaji")], FUN=length)) | ||
| Line 180: | Line 182: | ||
Parameter = "SD", | Parameter = "SD", | ||
aggregate(Hg$Conc, by = Hg[c("Helcom-alue", "Naytelaji")], FUN=sd) | aggregate(Hg$Conc, by = Hg[c("Helcom-alue", "Naytelaji")], FUN=sd) | ||
), | |||
data.frame( | |||
Parameter = "MeanLog", | |||
aggregate(log10(Hg$Conc), by = Hg[c("Helcom-alue", "Naytelaji")], FUN=mean) | |||
), | |||
data.frame( | |||
Parameter = "SDLog", | |||
aggregate(log10(Hg$Conc), by = Hg[c("Helcom-alue", "Naytelaji")], FUN=sd) | |||
), | ), | ||
data.frame( | data.frame( | ||
Parameter = "N", | Parameter = "N", | ||
aggregate(Hg$Conc, by = Hg[c("Helcom-alue", "Naytelaji")], FUN= | aggregate(Hg$Conc, by = Hg[c("Helcom-alue", "Naytelaji")], FUN=length) | ||
) | ) | ||
) | ) | ||
Revision as of 12:57, 19 March 2018
| [show] |
|---|
Question
What are the current total mercury concentrations in different fish species in Finland?
Answer
{{#opasnet_base_link:Op_en4004}}
The results contain e.g. the KERTY database with thousands of mercury measurements from different parts of Finland, both inland lakes and the Baltic Sea.
Rationale
Data
The table below is data from EU-kalat 1 project (2003), and the values are about Baltic herring. (Other fish species were measured but data are not here).[1]
| Obs | Age | Size | Concentration |
|---|---|---|---|
| 1 | 1 | 1 | 0.0005 |
| 2 | 2 | 1 | 0.0005 |
| 3 | 3 | 1 | 0.0005 |
| 4 | 2 | 1 | 0.005 |
| 5 | 2 | 4 | 0.005 |
| 6 | 5 | 2 | 0.006 |
| 7 | 2 | 1 | 0.007 |
| 8 | 1 | 1 | 0.008 |
| 9 | 2 | 1 | 0.008 |
| 10 | 3 | 2 | 0.008 |
| 11 | 3 | 2 | 0.008 |
| 12 | 2 | 2 | 0.009 |
| 13 | 3 | 2 | 0.009 |
| 14 | 3 | 3 | 0.01 |
| 15 | 2 | 2 | 0.011 |
| 16 | 3 | 2 | 0.011 |
| 17 | 1 | 1 | 0.012 |
| 18 | 2 | 1 | 0.014 |
| 19 | 4 | 2 | 0.015 |
| 20 | 4 | 3 | 0.015 |
| 21 | 1 | 1 | 0.018 |
| 22 | 3 | 2 | 0.018 |
| 23 | 2 | 1 | 0.02 |
| 24 | 2 | 2 | 0.021 |
| 25 | 2 | 2 | 0.022 |
| 26 | 4 | 2 | 0.022 |
| 27 | 6 | 4 | 0.022 |
| 28 | 8 | 4 | 0.024 |
| 29 | 6 | 3 | 0.025 |
| 30 | 5 | 2 | 0.026 |
| 31 | 4 | 3 | 0.027 |
| 32 | 5 | 4 | 0.027 |
| 33 | 8 | 4 | 0.027 |
| 34 | 6 | 4 | 0.028 |
| 35 | 5 | 3 | 0.03 |
| 36 | 2 | 1 | 0.031 |
| 37 | 7 | 5 | 0.031 |
| 38 | 3 | 2 | 0.033 |
| 39 | 10 | 4 | 0.033 |
| 40 | 6 | 3 | 0.038 |
| 41 | 10 | 5 | 0.038 |
| 42 | 12 | 5 | 0.039 |
| 43 | 6 | 4 | 0.04 |
| 44 | 5 | 5 | 0.043 |
| 45 | 7 | 5 | 0.043 |
| 46 | 7 | 4 | 0.048 |
| 47 | 6 | 3 | 0.049 |
| 48 | 6 | 3 | 0.049 |
| 49 | 3 | 2 | 0.051 |
| 50 | 9 | 5 | 0.054 |
| 51 | 9 | 4 | 0.066 |
| 52 | 7 | 5 | 0.069 |
| 53 | 3 | 5 | 0.072 |
| 54 | 11 | 4 | 0.083 |
| 55 | 11 | 5 | 0.107 |
Calculations
- 6.10.2017 data.frame Hg initiated [2]
See also
Metalli_RAPORTTI Report EU-kalat about heavy metals in fish in Finland- Benefit-risk assessment of methyl mercury and omega-3 fatty acids in fish
- Miller WaterAirSoilPollut 2013
- NordicHg_B1761 Mercury in Nordic ecosystems 2007
- Instructions for Kerty data Data received from Kerty 3th Sep 2014. Further questions can be asked from Simo Salo (at ymparisto.fi)