Ero sivun ”Talvivaaran kaivoksen mittauksia” versioiden välillä

Versio 18. marraskuuta 2012 kello 16.39

Vastuuvapautus: Tämä sivu käsittelee Talvivaaran kaivosta tai siihen liittyviä asioita. Sivustoa ylläpidetään joukkoistamalla tiettyjen sääntöjen mukaisesti (mm. keskustelusäännöt), eikä kukaan yksittäinen henkilö tai organisaatio vastaa sivujen sisällöstä. Sisällöstä vastaa jokainen kirjoittaja itse omien tekstiensä osalta. Sivut eivät myöskään edusta minkään valtion laitoksen tai organisaation virallista kantaa, vaan ne löytyvät kunkin laitoksen omilta sivuilta.

Tämä sivu on tutkimus. Sivutunniste: Op_fi3280
Moderaattori:Jouni (katso kaikki) Sivun edistymistä ei ole arvioitu. Arvostuksen määrää ei ole arvioitu (ks. peer review).
Lisää dataa {{#opasnet_base_link:Op_fi3280}}

Talvivaaran kaivoksen mittauksia esittelee mittaustuloksia ja havaintoja kaivoksen lähiympäristöstä.

Kysymys

Mitä mittauksia on tehty ja mitä tuloksia saatu Talvivaaran ympäristössä?

Vastaus

Tällä sivulla oli aiemmin suuri taulukko tuloksista. Nyt se on luettavuussyistä poistettu, mutta voit nähdä tulokset eri tavoilla:

Piirrä karttoja ja kuvaajia haluamastasi altisteesta tämän sivun kohdassa Laskenta.
Mene alla olevasta linkistä tietokantakäyttöliittymään ja katso uusinta 18.11.2012 ladattua dataa (series_id 7408).

Katso tämän sivun aiempaa versiota, jossa taulukko vielä on.

Perustelut

Data

Mittauspisteiden sijainnit leveys- ja pituusasteina
Name	LA	LO	Description
Lumijärvi	63.942533	27.936602
Ylä-Lumijärvi	63.945662	27.971964
Kivijärvi	63.931893	27.907677
Lumijoki, Kivijärveä lähin piste	63.925618	27.929521
Lumijoki, keskimmäinen piste	63.936103	27.926474
Lumijoki, Lumijärveä lähin piste	63.945605	27.953424
Lampi kaivosalueen reunalla 1	63.986605	27.971706	Ei selvää, mihin lampeen viittaa? Näissä koordinaateissa sijaitsee Haukilampi.
Lampi kaivosalueen reunalla 2	63.992419	27.983251	Ei selvää, mihin lampeen viittaa? Näissä koordinaateissa sijaitsee Kärsälampi.
Puro Kalliojärven ja Salmisen välillä	64.004996	28.002756
Tuhkajoki	64.044465	28.104916
Kortelampi	63.950091	27.984538
Kärsälampi	63.992843	27.984581
Rengaskaivo, Pappila	63.968536	27.953897
Kalliojärvi
Kalliojoki
Lumijoki
Salmisenpuro

Huom: epäselvyys STUK:in kartassa näkyy kaksi lähekkäistä näytteenottopaikkaa, jotka ovat kaivospiirin reunalla ja niihin viittaa vain yhdessä nimike "Lampi". Voisi siis kuvitella, että ne ovat "Lampi kaivosalueen reunalla 1" ja "Lampi kaivosalueen reunalla 2". Pisteiden osoittamissa kohdissa ovat Haukilampi ja Kärsälampi, sivulla on kuitenkin myös mittaus, joka on nimetty Kärsälammeksi. Haukilampea ei löydy mittauksista.

Nikkeli- ja sinkkipitoisuuksia lähialueen pintavesistä (Talvivaaran omia mittauksia) [1] [2]
Uraanipitoisuuksia lähialueen vesistä ja pohjasedimenteistä (STUK:n mittauksia). [3]

Laskenta

Tämän koodin avulla voi piirtää kuvaajia mittaustuloksista.

+ Näytä koodi - Piilota koodi

library(OpasnetUtils)
library(ggplot2)
library(xtable)
library(OpasnetUtilsExt)
library(rgdal)
library(maptools)
library(RColorBrewer)
library(classInt)
library(raster)


pitoisuus <- new("ovariable",
	name = "pitoisuus",
	data = tidy(op_baseGetData("opasnet_base", "Op_fi3280", series_id = 7408), objname= "pitoisuus")
)

# Remove empty rows.
pitoisuus@data <- pitoisuus@data[pitoisuus@data$pitoisuusResult != "" , ]

# Convert values that are below detection limit to ranges 0 - detection limit.
pitoisuus@data$pitoisuusResult <- gsub("<", "0-", pitoisuus@data$pitoisuusResult)

# Convert different time formats used to POSIXct time format.
temp <- as.POSIXct(pitoisuus@data$Aika, format = "%d.%m.%Y %H:%M")
temp <- ifelse(is.na(temp), as.POSIXct(pitoisuus@data$Aika, format = "%d.%m.%Y %H.%M"), temp)
temp <- ifelse(is.na(temp), as.POSIXct(pitoisuus@data$Aika, format = "%d.%m.%Y"), temp)
temp <- as.POSIXct(temp, origin = "1970-01-01 02:00.00 UTC")
pitoisuus@data$Aika <- temp

pitoisuus <- EvalOutput(pitoisuus, N = 1)

data <- pitoisuus@output[pitoisuus@output$Altiste == altiste, ]
data.time <- data

########### Concentrations on map

locs <- data.frame(Paikka = c(
	"Lumijärvi",
	"Ylä-Lumijärvi",
	"Kivijärvi",
	"Lumijoki, Kivijärveä lähin piste",
	"Lumijoki, keskimmäinen piste",
	"Lumijoki, Lumijärveä lähin piste",
	"Lampi kaivosalueen reunalla 1",
	"Lampi kaivosalueen reunalla 2",
	"Puro Kalliojärven ja Salmisen välillä",
	"Tuhkajoki",
	"Kortelampi",
	"Kärsälampi",
	"Rengaskaivo, Pappila",
	"Kalliojärvi",
	"Kalliojoki",
	"Lumijoki",
	"Salmisenpuro"
),
LA = c(
	63.942533,
	63.945662,
	63.931893,
	63.925618,
	63.936103,
	63.945605,
	63.986605,
	63.992419,
	64.004996,
	64.044465,
	63.950091,
	63.992843,
	63.968536,
	64,
	64.01,
	64.02,
	64.03
),
LO = c(
	27.936602,
	27.971964,
	27.907677,
	27.929521,
	27.926474,
	27.953424,
	27.971706,
	27.983251,
	28.002756,
	28.104916,
	27.984538,
	27.984581,
	27.953897,
	28,
	28.01,
	28.02,
	28.03
))

data <- merge(data, locs)
result <- "pitoisuusResult"

title <- paste(altiste, "Talvivaaran läheisyydessä (", data$pitoisuusUnit[1], ")", sep = "")
coordinates(data)=c("LO","LA")

# Plot the data
proj4string(data)<-("+init=epsg:4326")
epsg4326String <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
shp<-spTransform(data,epsg4326String)

#Create blank raster
rast<-raster()

#Set raster extent to that of point data
extent(rast)<-extent(shp)

#Choose number of columns and rows
ncol(rast) <- 64
nrow(rast) <- 64

#Rasterize point data
rast2<-rasterize(shp, rast, shp[[result]], fun=mean)

start <- min(shp[[result]])
end <- max(shp[[result]])

steps <- approx(c(start,end),n=6)$y
colors <- rev(rainbow(length(steps), start=0, end=0.50))

par(mfrow=c(6,1), mar=c(3,1,0,1), cex=1.5)

colorstrip <- function(colors, labels)
{
count <- length(colors)
m <- matrix(1:count, count, 1)
image(m, col=colors, ylab="", axes=FALSE)
axis(1,approx(c(0, 1), n=length(labels))$y, labels)
}

cat("<span style='font-size: 1.2em;font-weight:bold;'>", title, "</span>\n")


colorstrip(colors, steps)

#Plot data

google.show_raster_on_maps(rast2, col=colors, style="height:500px;")

########### Concentrations in time

ggplot(data.time, aes(x = Aika, y = pitoisuusResult, colour = Paikka)) + 
	geom_point(size = 5) + 
	theme_grey(base_size = 24) +
	labs(
		title = paste(altiste, "Talvivaaran läheisyydessä"), 
		y = paste("Pitoisuus,", data$pitoisuusUnit[1]),
		x = "Aika"
	)

Katso myös

Viitteet

Aiheeseen liittyviä tiedostoja