Mercury concentrations in fish in Finland: Difference between revisions
Jump to navigation
Jump to search
| (3 intermediate revisions by the same user not shown) | |||
| Line 93: | Line 93: | ||
conc_mehg <- Ovariable( | conc_mehg <- Ovariable( | ||
"conc_mehg", | "conc_mehg", | ||
dependencies = data.frame(Name="Hg", Ident="Op_en4004/ | dependencies = data.frame(Name="Hg", Ident="Op_en4004/Hg2"), | ||
formula = function(...) { | formula = function(...) { | ||
Hg <- Hg[!(is.na(Hg$MeanLog) | is.na(Hg$SDLog)),] | |||
out <- data.frame() | |||
out | for(i in 1:nrow(Hg)) { | ||
for(i in 1:nrow( | |||
if(openv$N>1) { | if(openv$N>1) { | ||
tmp <- data.frame( | |||
Iter=1:openv$N, | Iter=1:openv$N, | ||
Result=10^rnorm(openv$N, | Result=10^rnorm(openv$N, Hg$MeanLog[i], Hg$SDLog[i]) | ||
) | ) | ||
} else { | } else { | ||
tmp <- data.frame(Result = 10^Hg$MeanLog[i]) | |||
} | } | ||
out <- rbind(out, | |||
data.frame(Area=as.character( | data.frame(Area=as.character(Hg$Area[i]), | ||
Kala=as.character( | Kala=as.character(Hg$Kala[i]), | ||
tmp | |||
) | ) | ||
) | ) | ||
} | } | ||
out <- Ovariable(output=out, marginal=colnames(out) != "Result") * | |||
out <- Ovariable(output= | |||
Ovariable( | Ovariable( | ||
output = data.frame( | output = data.frame( | ||
| Line 268: | Line 264: | ||
==== Updated Hg code ==== | ==== Updated Hg code ==== | ||
There was a unit conversion error in the Hg code, so an update was made on 2021-02-24 [http://en.opasnet.org/en-opwiki/index.php?title=Special:RTools&id= | There was a unit conversion error in the Hg code, so an update was made on 2021-02-24 [http://en.opasnet.org/en-opwiki/index.php?title=Special:RTools&id=6b4wREivb5piSXFs]. | ||
<rcode name="Hg2" label="Initiate data.frame Hg"> | <rcode name="Hg2" label="Initiate data.frame Hg"> | ||
| Line 278: | Line 274: | ||
"Op_en4004", subset = "Kerty database")[c(8,9,12,15,16,17,20,21,22)] | "Op_en4004", subset = "Kerty database")[c(8,9,12,15,16,17,20,21,22)] | ||
Hg <- Hg[Hg$Näyteosa=="lihas",] | Hg <- Hg[Hg$Näyteosa=="lihas",] | ||
Hg$Tulos <- as.numeric(Hg$Tulos) # 140 results has type L0.05, T0.1 or similar and are turned into NA. | Hg$Tulos <- as.numeric(as.character(Hg$Tulos)) # 140 results has type L0.05, T0.1 or similar and are turned into NA. | ||
Hg$Tulos <- ifelse(Hg$Yksikkö=="µg/kg", Hg$Tulos * 0.001, Hg$Tulos) | Hg$Tulos <- ifelse(Hg$Yksikkö=="µg/kg", Hg$Tulos * 0.001, Hg$Tulos) | ||
Hg$Yksikkö[Hg$Yksikkö=="µg/kg"] <- "mg/kg" | Hg$Yksikkö[Hg$Yksikkö=="µg/kg"] <- "mg/kg" | ||
summary(Hg) | oprint(summary(Hg)) | ||
aggregate(Hg$Tulos, by = Hg[c("Helcom-alue")], FUN=length) | oprint(aggregate(Hg$Tulos, by = Hg[c("Helcom-alue")], FUN=length)) | ||
aggregate(Hg$Tulos, by = Hg[c("Näytelaji")], FUN=length | oprint(aggregate(Hg$Tulos, by = Hg[c("Näytelaji")], FUN=length)) | ||
# summary(lm(Tulos ~ `Helcom-alue` + Näytelaji, data = Hg)) | # summary(lm(Tulos ~ `Helcom-alue` + Näytelaji, data = Hg)) | ||
| Line 316: | Line 311: | ||
) | ) | ||
) | ) | ||
Hg$Compound <- "Mercury (µg/g)" | Hg$Compound <- "Mercury (µg/g)" | ||
Hg <- reshape(Hg, idvar = c("Area","Näytelaji","Compound"), timevar = "Parameter", v.names = " | Hg <- reshape(Hg, idvar = c("Area","Näytelaji","Compound"), timevar = "Parameter", v.names = "x", direction = "wide") | ||
colnames(Hg)[colnames(Hg)=="Näytelaji"] <- " | colnames(Hg)[colnames(Hg)=="Näytelaji"] <- "Kala" | ||
colnames(Hg) <- gsub(" | colnames(Hg) <- gsub("x\\.", "", colnames(Hg)) | ||
oprint(Hg) | oprint(Hg) | ||
Latest revision as of 15:13, 14 March 2021
| [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
conc_mehg for PFAS disease burden
This code was developed for op_fi:PFAS-yhdisteiden tautitaakka. It contains all fish species from Hg table but only Gulf of Finland and Bothnian Bay as areas.
conc_mehg for Goherr assessment
This code contains only Baltic herring and salmon as fish species.
- 6.10.2017 data.frame Hg initiated [2]
Updated Hg code
There was a unit conversion error in the Hg code, so an update was made on 2021-02-24 [3].
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)