Mercury concentrations in fish in Finland: Difference between revisions

From Opasnet
Jump to navigation Jump to search
 
(5 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/Hg"),
   dependencies = data.frame(Name="Hg", Ident="Op_en4004/Hg2"),
   formula = function(...) {
   formula = function(...) {
     out <- Hg[(Hg$Helcom.alue %in% c("C-Selkämeri","F-Suomenlahti") |
     Hg <- Hg[!(is.na(Hg$MeanLog) | is.na(Hg$SDLog)),]
              Hg$Naytelaji %in% c("lohi","kirjolohi")) & !is.na(Hg$Result), -5]
     out <- data.frame()
     out <- merge(out[out$Parameter=="MeanLog",], out[out$Parameter=="SDLog",], by=c("Helcom.alue","Naytelaji"))
     for(i in 1:nrow(Hg)) {
    tmp <- data.frame()
     for(i in 1:nrow(out)) {
       if(openv$N>1) {
       if(openv$N>1) {
         tmp2 <- data.frame(
         tmp <- data.frame(
           Iter=1:openv$N,
           Iter=1:openv$N,
           Result=10^rnorm(openv$N, out$Result.x[i], out$Result.y[i])
           Result=10^rnorm(openv$N, Hg$MeanLog[i], Hg$SDLog[i])
         )
         )
       } else {
       } else {
         tmp2 <- data.frame(Result = 10^out$Result.x[i])
         tmp <- data.frame(Result = 10^Hg$MeanLog[i])
       }
       }
       tmp <- rbind(tmp,
       out <- rbind(out,
                   data.frame(Area=as.character(out$Helcom.alue[i]),
                   data.frame(Area=as.character(Hg$Area[i]),
                               Kala=as.character(out$Naytelaji[i]),
                               Kala=as.character(Hg$Kala[i]),
                               tmp2
                               tmp
                   )
                   )
       )
       )
     }
     }
    tmp$Area[tmp$Area=="C-Selkämeri"] <- "Suomi"
     out <- Ovariable(output=out, marginal=colnames(out) != "Result") *  
    tmp$Area[tmp$Area %in% c("","F-Suomenlahti")] <- "Porvoo"
     out <- Ovariable(output=tmp, marginal=colnames(tmp) != "Result") *  
       Ovariable(
       Ovariable(
         output = data.frame(
         output = data.frame(
Line 259: Line 255:
colnames(Hg)[colnames(Hg)=="x"] <- "Result"
colnames(Hg)[colnames(Hg)=="x"] <- "Result"
Hg$Compound <- "Mercury (µg/g)"
Hg$Compound <- "Mercury (µg/g)"
oprint(Hg)
objects.store(Hg)
cat("Data.frame Hg stored.\n")
</rcode>
==== 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=6b4wREivb5piSXFs].
<rcode name="Hg2" label="Initiate data.frame Hg">
# This is code Op_en4004/Hg2 on page [[Mercury concentrations in fish in Finland]]
library(OpasnetUtils)
Hg <- opbase.data(
  "Op_en4004", subset = "Kerty database")[c(8,9,12,15,16,17,20,21,22)]
Hg <- Hg[Hg$Näyteosa=="lihas",]
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$Yksikkö[Hg$Yksikkö=="µg/kg"] <- "mg/kg"
oprint(summary(Hg))
oprint(aggregate(Hg$Tulos, by = Hg[c("Helcom-alue")], FUN=length))
oprint(aggregate(Hg$Tulos, by = Hg[c("Näytelaji")], FUN=length))
# summary(lm(Tulos ~ `Helcom-alue` + Näytelaji, data = Hg))
# The different species (Näytelaji) vary generally more than +-0.35 mg/kg, while
# the different Helcom areas vary less than +-0.19
# Because we are mainly interested in Porvoo vs other Finland, we will categorise is accordingly.
Hg$Area <- ifelse(Hg$`Helcom-alue`=="F-Suomenlahti", "Gulf of Finland","Finland")
Hg <- Hg[!is.na(Hg$Tulos),]
Hg$Tulos[Hg$Tulos==0] <- 0.002 # Approximate level of detection
Hg <- rbind(
  data.frame(
    Parameter = "Mean",
    aggregate(Hg$Tulos, by = Hg[c("Area", "Näytelaji")], FUN=mean)
  ),
  data.frame(
    Parameter = "SD",
    aggregate(Hg$Tulos, by = Hg[c("Area", "Näytelaji")], FUN=sd)
  ),
  data.frame(
    Parameter = "MeanLog",
    aggregate(log10(Hg$Tulos), by = Hg[c("Area", "Näytelaji")], FUN=mean)
  ),
  data.frame(
    Parameter = "SDLog",
    aggregate(log10(Hg$Tulos), by = Hg[c("Area", "Näytelaji")], FUN=sd)
  ),
  data.frame(
    Parameter = "N",
    aggregate(Hg$Tulos, by = Hg[c("Area", "Näytelaji")], FUN=length)
  )
)
Hg$Compound <- "Mercury (µg/g)"
Hg <- reshape(Hg, idvar = c("Area","Näytelaji","Compound"), timevar = "Parameter", v.names = "x", direction = "wide")
colnames(Hg)[colnames(Hg)=="Näytelaji"] <- "Kala"
colnames(Hg) <- gsub("x\\.", "", colnames(Hg))


oprint(Hg)
oprint(Hg)

Latest revision as of 15:13, 14 March 2021



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]

Mercury in Baltic herring(mg /kg f.w.)
ObsAgeSizeConcentration
1110.0005
2210.0005
3310.0005
4210.005
5240.005
6520.006
7210.007
8110.008
9210.008
10320.008
11320.008
12220.009
13320.009
14330.01
15220.011
16320.011
17110.012
18210.014
19420.015
20430.015
21110.018
22320.018
23210.02
24220.021
25220.022
26420.022
27640.022
28840.024
29630.025
30520.026
31430.027
32540.027
33840.027
34640.028
35530.03
36210.031
37750.031
38320.033
391040.033
40630.038
411050.038
421250.039
43640.04
44550.043
45750.043
46740.048
47630.049
48630.049
49320.051
50950.054
51940.066
52750.069
53350.072
541140.083
551150.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.

+ Show code

conc_mehg for Goherr assessment

This code contains only Baltic herring and salmon as fish species.

+ Show code


+ Show code

  • 6.10.2017 data.frame Hg initiated [2]

+ Show code

Updated Hg code

There was a unit conversion error in the Hg code, so an update was made on 2021-02-24 [3].

+ Show code

See also

References

Retrieved from "https://en.opasnet.org/index.php?title=Mercury_concentrations_in_fish_in_Finland&oldid=43778"