__________________________________________________________ Rangifer tarandus ) from northern and central Lapland, Sweden, 1983-2005 Swedish monitoring programme in terrestrial biota Time trends of metals in liver and muscle of reindeer (

(1)

Time trends of metals in liver and muscle of reindeer (Rangifer tarandus) from northern and central Lapland,

Sweden, 1983-2005

Swedish monitoring programme in terrestrial biota

__________________________________________________________

SWEDISH · MUSEUM · OF · NATURAL · HISTORY Department of Contaminant Research

P.O. Box 50007 SE-104 05 Stockholm

2007-05-11

__________________________________________________________

(2)

(3)

SAKRAPPORT

Time trends of metals in liver and muscle of reindeer (Rangifer tarandus) from northern and central Lapland, Sweden, 1983-2005

Swedish monitoring programme in terrestrial biota

Avtal nr 214 0630 Miljögifter i biota - fjäll

Department of Contaminant Research

Swedish Museum of Natural History P O Box 50007

SE-104 05 Stockholm

Department of Chemistry National Veterinary Institute

SE-750 07 Uppsala 2007-05-11

(4)

TIME TRENDS OF METALS IN LIVER AND MUSCLE OF REINDEER (Rangifer tarandus) FROM NORTHERN AND CENTRAL LAPLAND, SWEDEN, 1983-2005

Compiled by

Tjelvar Odsjö, Jannikke Räikkönen, Anders Bignert Department of Contaminant Research

Swedish Museum of Natural History

INTRODUCTION

The long-term monitoring of bio-accumulating contaminants in biota from terrestrial environments in Sweden as part of the Swedish National Environmental Monitoring

Programme is based on analysis of organs and tissues of different animal species collected in certain pristine areas of the Swedish mainland (Odsjö and Olsson 1979 a,b). In the

mountainous area of north-western Sweden, reindeer (Rangifer tarandus) is chosen as a representative indicator for the fauna living in that part of the country. Samples of reindeer have continuously been collected in three districts since the early 1980s. Later, the

Declaration on the Protection of the Arctic environment established an Arctic Monitoring and Assessment Programme (AMAP) to monitor levels and assess effects of anthropogenic

pollutants in components of the Arctic environment. The Programme recommends that collection of baseline data for heavy metals and radionuclides in caribou/reindeer should be mandatory for participating states due to the importance of that species in the diets of northern native people (AMAP 1993). The current material of reindeer from northern and central Lapland partly satisfies the Swedish participation in the AMAP programme.

The herbivorous reindeer spend the summer time in the westernmost part of the high

mountain areas. Summer diets include grasses, sedges, twigs, leaves and mushrooms. During autumn they migrate eastwards to winter grounds in the central coniferous forest areas of the country, where they primarily feed on lichens, which are noted for their ability to accumulate nutrients and contaminants from the air. Winter diets also include sedges and twigs.

The actual material of reindeer have earlier been utilised for analyses of e.g. radiocesium in a study of effects of the fallout of Cs-132 from the Chernobyl accident in 1986 (Forberg et al.

1992) and for studies of time trends of levels of HCHs and HCB (Odsjö et al. 1998).

AIM

The aims of the study are to present long-term trend series of bio-accumulated concentrations of metals for the period 1983-2005 and to summarise results from statistical analyses. The analysed elements are Al, Ca, Cd, Co, Cr, Cu, Fe, Mg, Mn, Ni, V, Zn, Hg, Pb and Mo.

The present paper reports on trend studies of concentrations of metals in muscle and liver of reindeer collected from Gabna, Lævas and Girjas Sámi Villages, the northernmost of the three collection districts presented below. The report also includes presentation of concentrations of metals in muscle and liver of reindeer collected in 1996-2005 in Ran and Gran Sámi Villages, the central of the three collection districts presented below. Reindeer from this district have so

(5)

far been analysed from nine years only. For aluminium there are data only from two years why the extended statistics are not relevant and thus reduced. Tissues of reindeer collected 1983-1995 in Ran and Gran, not yet analysed, are stored frozen in the Environmental

Specimen Bank at the Swedish Museum of Natural History and are available for retrospective analyses and trend studies. So are specimens of reindeer from the southernmost district, but have not yet been chemically analysed.

MATERIAL AND METHODS Sampling areas

Collection of specimens of reindeer aimed at chemical analyses of bio-accumulating noxious substances is annually carried out in three districts along the Swedish, easternmost mountain chain.

* The northernmost district, Gabna, Lævas and Girjas Sámi Villages, reaches from the Swedish/Norwegian border eastwards to the central forest areas in the northern part of Lapland. (Figure 1). Collection of specimens started in Rensjön and Aitejåkk in 1981.

* The central district, Ran and Gran Sámi Villages, is situated in the mountain area in central Lapland and reaches from the Swedish/Norwegian border eastwards to the central forest areas (Figure 1). Collection of specimens started in Ammarnäs in 1983.

* The southernmost district, Handölsdalens and Mittådalens Sámi Villages, covers the border areas between the provinces of Jämtland and Härjedalen from the

Swedish/Norwegian border eastwards to the forest areas in the central part of the country. (Figure 1). Collection of specimens started in Ottsjölägret in 1982 and was moved over to Ljungris in 1989 due to changed slaughter schedules and artificial feeding prior to slaughter as a consequence of deposition of cesium polluted from Chernobyl, the Ukraine in 1986.

Since the start of the monitoring programme, collection of specimens of muscle, liver, kidney and left under jaw with teeth (for age determination) from at least 50 male reindeer has been carried out annually at regular slaughtering in each of the three districts. In 1983-1986, the muscle samples were taken from the mandibles. After that, routines were changed and muscle samples were taken from the front leg tibia. The change of muscle samples convinced us to exclude analytical results from mandible muscle in the presentation below since they were not comparable.

The slaughtering was carried out between the end of August and mid-October, mostly in the mid of September. The age of the males from which samples were taken ranged mainly between 2 and 4 years. Most specimens were 3 years old. Ages were determined in the lab after slaughter. In 1998, the ordinary slaughter in the northernmost district was postponed from early September to early November, which should be noted when evaluating the time trends of levels. The fact that reindeer move from summer areas in the mountains to wintering areas in the forest means a change of diet, which might influence upon the exposition of metals via food. A later date for slaughtering means a longer foraging on winter diets that might contain higher concentrations of metals.

(6)

Sample preparation

Ten male reindeer, three years old, were selected for chemical analyses annually. Multi- element analyses of metals were carried out on muscle and liver on an individual basis.

Approximately 15-20 g tissue was prepared of both muscle and liver. To avoid irrelevant contamination, the surface layer of each tissue sample was removed by ceramic knives at preparation for analyses. For extra quality assurance, the surface layer was removed once again by use of knives made of titanium (Ti) when the tissues were finally prepared at the analytical laboratory.

Chemical analyses

For determination of Hg, continuous hydride generation by flow injection analysis and atomic absorption spectrometry was used (Galgan & Frank 1988). For pre-treatment of the samples, see: Frank 1976, 1983, 1988. Simultaneous analysis of the other 13 elements mentioned above was performed by use of a direct-current plasma-atomic emission spectrometer, DCP- AES (SpectraSpan IIIA, Applied Research Laboratories Inc., Valencia, CA, USA). For pre- treatment of the samples, see: Frank and Petersson 1983, 1985. The chemical analyses were carried out by the Department of Chemistry, National Veterinary Institute, Uppsala.

Statistical treatment and graphical presentation

Statistical treatment and graphical presentation have been carried out according to Bignert (1998).

Trend detection

One of the main purposes of the monitoring programme is to detect trends. The trend detection is carried out in three steps.

1. Log-linear regression analyses

Log-linear regression analyses is performed both for the entire investigated time period and for time series longer than ten years, also for the recent ten years. The slope of the line

describes the yearly percentage change. A slope of 5% implies that the concentration is halved in 14 years whereas 10% corresponds to a similar reduction in 7 years and 2% in 35 years.

See Table 2, below.

Table 2. The approximate number of years required to double or half the initial concentration assuming a continuous annual change of 1, 2, 3, 4, 5, 7, 10, 12, 15 or 20% a year.

1% 2% 3% 4% 5% 7% 10% 12% 15% 20%

Increase 70 35 24 18 14 10 7 6 5 4

Decrease 69 35 23 17 14 10 7 6 4 3

2. Non-parametric trend test

The regression analyses presuppose, among other thing, that the regression line gives a good description of the trend. The leverage effect of points in the end of the line is also a well- known fact. An exaggerated slope, caused 'by chance' by a single or a few points in the end of the line, increases the risk of a false significant result when no real trend exist. A non-

parametric alternative to the regression analysis is the Mann-Kendall trend test (Gilbert, 1987, Helsel & Hirsch, 1995, Swertz, 1995). This test has generally lower power than the regression analysis and does not take differences in magnitude of the concentrations into account, it only counts the number of consecutive years where the concentration increases or decreases compared with the year before. If the regression analysis yields a significant result but not the

(7)

Mann-Kendall test, the explanation could be either that the latter test has lower power or that the influence of endpoints in the time series has become unwarrantable great on the slope.

Hence, the eighth line reports Kendall's 'τ', and the corresponding p-value. The Kendall's 'τ' ranges from 0 to 1 like the traditional correlation coefficient ‘r’ but will generally be lower.

‘Strong’ linear correlation of 0.9 or above corresponds to τ-values of about 0.7 or above (Helsel and Hirsch, 1995, p. 212). EPA recommended this test for use in water quality monitoring programmes with annual samples, in an evaluation comparing several other trend tests (Loftis et al. 1989).

3. Non-linear trend components

An alternative to the regression line in order to describe the development over time would be some kind of smoothed line. The smoother applied here is a simple 3-point running mean smoother fitted to the annual geometric mean values. In cases where the regression line is badly fitted the smoothed line may be more appropriate. The significance of this line is tested by means of an Analysis of Variance where the variance explained by the smoother and by the regression line is compared with the total variance. This procedure is used at assessments at ICES and is described by Nicholson et al., 1995.

Outliers and values below the detection limit

Observations too far from the regression line considering from what could be expected from the residual variance around the line is subjected to special concern. These deviations may be caused by an atypical occurrence of something in the physical environment, a changed pollution load or errors in the sampling or analytical procedure. The procedure to detect suspected outliers in this presentation is described by Hoaglin and Welsch (1978). It makes use of the leverage coefficients and the standardised residuals. The standardised residuals are tested against a t.05 distribution with n-2 degrees of freedom. When calculating the ith

standardised residual the current observation is left out implying that the ith observation does not influence the slope nor the variance around the regression line. The suspected outliers are merely indicated in the figures and are included the statistical calculations except in a few cases, pointed out in the figures.

Values reported below the detection limit is substituted using the ‘robust’ method suggested by Helsel & Hirsch (1995) p 362, assuming a lognormal distribution within a year. N.B. a minimum of three values above LOQ (Limit of quantification) are required for this substitution, years with fewer results above LOQ are not represented in the figures.

Legend to the plots

The analytical results from each of the investigated elements are displayed in figures. A separate plot represents each site.

The plot displays the geometric mean concentration of each year (circles) together with the individual analyses (small dots) and the 95% confidence intervals of the geometric means.

The overall geometric mean value for the time series is depicted as a horizontal, thin, dashed line.

The trend is presented by one regression line (plotted if p < 0.10, two-sided regression analysis). Ten years is often too short a period to statistically detect a trend unless it is of considerable magnitude. Furthermore, the residual variance around the line compared to the residual variance for the entire period will indicate if the sensitivity have increased as a result

(8)

of e.g. an improved sampling technique or that problems in the chemical analysis have disappeared.

A smoother is applied to test for non-linear trend components. The smoothed line is plotted if p < 0.10. A broken line or a dashed line segment indicates a gap in the time series with a missing year.

The log-linear regression lines fitted through the geometric mean concentrations follow smooth exponential functions.

A cross inside a circle, indicates a suspected outlier, see above. The suspected outliers are merely indicated in the figures and are included in the statistical calculations except in a few cases, pointed out in the figures.

Each plot has a header with element, species name, tissue and sampling locality. Below the header of each plot the results from several statistical calculations are reported:

n(tot) = The first line reports the total number of analyses included together with the number of years ( n(yrs) = ).

m = The overall geometric mean value together with its 95% confidence interval is reported on the second line of the plot (N.B. d.f.= n of years - 1).

slope = reports the slope, expressed as the yearly percentage change together with its 95%

confidence interval.

SD(lr) = reports the square root of the residual variance around the regression line, as a measure of between-year variation, together with the lowest detectable change in the current time series with a power of 80%, one-sided test, α=0.05. The last figure on this line is the estimated number of years required to detect an annual change of 5% with a power of 80%, one-sided test, α=0.05.

power = reports the power to detect a log-linear trend in the time series (Nicholson & Fryer, 1991). The first figure represents the power to detect an annual change of 5% with the number of years in the current time series. The second figure is the power estimated as if the slope where 5% a year and the number of years were ten. The third figure is the lowest detectable change for a ten-year period with the current between year variation at a power of 80%.

y(05) = reports the concentration estimated from the regression line for the last year together with a 95% confidence interval, e.g. y(05)=40.6(36.6,45.0) is the estimated concentration of year 2005 where the residual variance around the regression line is used to calculate the confidence interval. Provided that the regression line is relevant to describe the trend, the residual variance might be more appropriate than the within-year variance in this respect.

r²= reports the coefficient of determination (r²) together with a p-value for a two-sided test (H0: slope = 0) i.e. a significant value is interpreted as a true change, provided that the assumptions of the regression analysis is fulfilled.

tao = reports Kendall's 'τ', and the corresponding p-value.

SD(sm) = reports the square root of the residual variance around the smoothed line. The significance of this line could be tested by means of an Analysis of Variance. The p-value is reported for this test. A significant result will indicate a non-linear trend component.

Below these nine lines are additional lines with information concerning the regression of the last ten years.

(9)

RESULT

Long-term trends of metals in liver and muscle of reindeer

The analytical and statistical results from Gabna, Lævas and Girjas Sámi Villages are

displayed in Figure 2-9, which visualise the trends of metal concentration in muscle and liver of reindeer from the period 1983-2005. The data of metal concentrations in muscle of reindeer from Gabna, Lævas and Girjas Sámi Villages, collected in 1983-1986, was excluded from the statistical calculations and is not visualised in the graphs. The reason is that muscle samples during the early part of the period were taken from the mandibles. From 1987 onwards, muscle samples were taken from the front leg tibia and we have reason to believe that concentrations in the both tissues are not fully comparable.

Significantly increasing trends were shown for the time series from Gabna, Lævas and Girjas Sámi Villages for the following elements: chromium (27% a year, muscle, 1994-2003), copper (1.7% a year, liver, 1983-2005), and mercury (9.2% a year, muscle, 1987-2005).

Significantly decreasing trends were shown for the time series from Gabna, Lævas and Girjas Sámi Villages for the following elements: calcium (-2.2% a year, muscle, 1996-2005 and - 0.57% a year, liver, 1983-2005), cadmium (-8.9% a year, muscle, 1987-2005), cobalt (-9.4% a year, muscle, 1996-2005 and –1.9% a year, liver, 1983-2005), copper (-0.59% a year, muscle, 1987-2005), magnesium (-0.36% a year, liver, 1983-2005), manganese (-1.3% a year, muscle, 1987-2005 and –2.2% a year, liver, 1983-2005), vanadium (-9.1% a year, muscle, 1987-2005 and -8.2% a year, liver, 1983-2001) and lead (-2.5% a year, liver, 1983-2005). 23 years of analyses are now available for liver and 19 years for muscle.

The analytical results from Ran and Gran Sámi Villages from 1996-2005 are displayed in Figure 10-16. The time series is continuously updated once a year when new material is collected and analysed.

Significantly increasing trends were shown for the time series from Ran and Gran Sámi Villages for the following elements: iron (3.2% a year, muscle, 1996-2005)

Significantly decreasing trends were shown for the time series from Ran and Gran Sámi Villages for the following elements: calcium (-3.0% a year, muscle, 1996-2005), iron (-3.6%

a year, liver, 1996-2005), nickel (-25% a year, liver, 1997-2000), mercury (-9.7% a year, liver, 1996-2005) and lead (-10% a year, liver, 1996-2005)

(10)

Gabna, Lævas and Girjas Sámi Villages

─────────────────────────────────────────────────────────

Aluminium (Al

) (Figure 2)

Residue levels of aluminium have been analysed in muscle of reindeer from the period 1987- 1998, and in liver from the period 1983-1998.

The aluminium concentrations in muscle and liver of reindeer show no significant log-linear trend during the period (parametric test).

The number of years required to detect an annual change of 5% was 30 years for both muscle and liver with a power to detect a 5% annual change varying between 0.10 and 0.18 for the full period.

The ANOVA test showed that the smoothed lines for concentrations of aluminium in both muscle and liver indicate a significant non-linear trend component (p<0.077, muscle;

p<0.002, liver).

The overall geometric mean value of aluminium in muscle and liver was 0.143 and 0.332 µg/g (fresh weight), respectively for the period 1987-1998 and 1983-1998.

_______________

Calcium (Ca)

(Figure 2)

The calcium concentrations in muscle of reindeer show no significant log-linear trend during the period 1987-2005 but for liver (1983-2005) (p<0.094; parametric test). The average annual decrease in liver was 0.57%. The calcium concentrations in muscle also show a significant decrease (2.2%, p<0.079) during the period 1996-2005.

The number of years required to detect an annual change of 5% was 13 and 9 years for muscle and liver respectively, with a power of 1.0 to detect a 5% annual change for the full period for both muscle and liver.

The ANOVA test showed that the smoothed line for concentrations of calcium in liver indicates a significant non-linear trend component (p<0.094).

The overall geometric mean value of calcium in muscle and liver was 38.8 and 42.3 µg/g (fresh weight) respectively, for the period 1987-2005 and 1983-2005.

_______________

Cadmium (Cd)

(Figure 3)

The cadmium concentrations in muscle of reindeer show a significant log-linear change during the period 1987-2005(p<0.008; parametric test), while the concentrations in liver show no significant change during the period (parametric test). The average annual decrease in muscle was 8.9%.

The number of years required to detect an annual change of 5% was 27 and 14 years for muscle and liver respectively, with a power to detect a 5% annual change varying between 0.27 and 1.0 for the full period.

The ANOVA test showed that the smoothed line for concentrations of cadmium in muscle indicates a significant non-linear trend component (p<0.074, muscle).

The overall geometric mean value of cadmium in muscle and liver was 0.010 and 0.440 µg/g (fresh weight) respectively, for the periods 1987-2005 and 1983-2005.

(11)

_______________

Cobalt (Co)

(Figure 3)

The cobalt concentrations in liver of reindeer show a significant log-linear change during the period 1983-2005 (-1.9% per year, p<0.009; parametric test). The concentrations in muscle decreased significantly during the period 1996-2005 (-9.4% per year, p<0.003).

The ANOVA test showed that the smoothed line for concentrations of cobalt in liver indicates a significant non-linear trend component (p<0.060).

The overall geometric mean value of cobalt in muscle and liver was 0.007 and 0.141 µg/g (fresh weight) respectively, for the periods 1987-2005 and 1983-2005.

_______________

Chromium (Cr)

(Figure 4)

The levels of chromium in muscle from 1996 and 2004 2005 were all below the detection limit and were excluded from any calculation of mean values and statistical calculations. So was the case for chromium levels in liver from 1981, 1982, 1987-89, 1992-96, 2004 and 2005.

The chromium concentrations in muscle of reindeer show a significant log-linear increase during the period 1994-2003 (27% per year, p<0.014, parametric test) most likely due to a low outlier in 1994.

The ANOVA test showed that the smoothed line for concentrations of chromium in liver indicates a significant non-linear trend component (p<0.001).

The overall geometric mean value of chromium in muscle and liver was 0.020 and 0.014 µg/g (fresh weight) respectively, for the periods 1987-2003 and 1983-2003.

_______________

Copper (Cu)

(Figure 4)

The copper concentrations in muscle and liver of reindeer show significant log-linear changes during the period 1987-2005 (muscle) and 1983-2005 (liver) (parametric test). The average annual decrease in muscle was 0.59% (p<0.076) and the average annual increase in liver was 1.7% (p<0.054).

The number of years required to detect an annual change of 5% was 8 and 15 years for muscle and liver respectively, with a power of 1.0 to detect a 5% annual change in both tissues for the full period.

The overall geometric mean value of copper in muscle and liver was 1.32 and 67.2 µg/g (fresh weight) respectively, for the periods 1987-2005 and 1983-2005.

_______________

(12)

Iron (Fe)

(Figure 5)

The iron concentrations in muscle and liver of reindeer show no significant log-linear change during the period 1983-2004 (parametric test).

The ANOVA test showed that the smoothed line for concentrations of iron in muscle indicates a significant non-linear trend component (p<0.043).

The overall geometric mean value of iron in muscle and liver was 29.4 and 140 µg/g (fresh weight) respectively, for the periods 1987-2005 and 1983-2005.

It should be noted that concentration of iron might vary with amount of blood and red blood capsules in the analysed samples.

_______________

Magnesium (Mg)

(Figure 5)

The magnesium concentrations in liver of reindeer show a significant log-linear decrease during the period 1983-2005 (0.36% per year, p<0.088, parametric test).

The number of years required to detect an annual change of 5% was 10 and 7 years for muscle and liver respectively, with a power of 1.0 for both tissues to detect a 5% annual change for the full period.

The ANOVA test showed that the smoothed lines for concentrations of magnesium in both muscle and liver indicate a significant non-linear trend component (p<0.015, muscle;

p<0.040, liver).

The overall geometric mean value of magnesium in muscle and liver was 202 and 183 µg/g (fresh weight) respectively, for the periods 1987-2005 and 1983-2005.

_______________

Manganese (Mn)

(Figure 6)

The manganese concentrations in muscle and liver of reindeer show a significant log-linear change during the period 1987-2005 and 1983-2005, respectively (-1.3% per year, p<0.053;

muscle, parametric test), (-2.2% per year, p<0.015, liver, parametric test).

The number of years required to detect an annual change of 5% was 11 and 15 years for muscle and liver respectively, with a power of 1.0 in both muscle and liver series to detect a 5% annual change for the full period.

The ANOVA test showed that the smoothed lines for concentrations of manganese in both muscle and liver indicate a significant non-linear trend component (p<0.088, muscle;

p<0.001, liver).

The overall geometric mean value of manganese in muscle and liver was 0.146 and 2.72 µg/g (fresh weight) respectively, for the periods 1987-2005 and 1983-2005.

_______________

(13)

Nickel (Ni)

(Figure 6)

The levels of nickel in muscle from 1988-1989, 1991, 2002, 2004 and 2005 were all below the detection limit and were excluded from any calculation of mean values and statistical calculations. So were levels in liver from 1987, 1989, 1992-96 and 2002-05.

The nickel concentrations in muscle and liver of reindeer show no significant log-linear change during the period (parametric test).

The number of years required to detect an annual change of 5% was 38 and 26 years for muscle and liver respectively, with a power to detect a 5% annual change of 0.08 and 0.12 for the full period.

The ANOVA test showed that the smoothed lines for concentrations of nickel in both muscle and liver indicate a significant non-linear trend component (p<0.047, muscle; p< 0.019, liver).

The overall geometric mean value of nickel in muscle and liver was 0.012 and 0.015 µg/g (fresh weight) respectively, for the periods shown in the graphs.

_______________

Vanadium (V)

(Figure 7)

The vanadium concentrations in liver of reindeer show a significant log-linear change during the period 1983-2001 (-8.2% per year, p<0.036 parametric test) and during 1992-2001 (-22%

per year, p< 0.012, parametric test). The concentrations in muscle also show a significant log- linear decrease during the period 1987-2005 (9.1% per year, p<0.049 parametric test)

The ANOVA test showed that the smoothed line for concentrations of vanadium in liver indicates a significant non-linear trend component (p< 0.003).

The number of years required to detect an annual change of 5% was 28 and 29 years for muscle and liver respectively, with a power to detect a 5% annual change varying between 0.07 and 0.14 for the full period. However, considering the large number of years missing from this data set, calculations of power are a blunt tool.

The overall geometric mean value of vanadium in muscle and liver was 0.005 and 0.007 µg/g (fresh weight) respectively during the period 1987-2005 and 1983-2001.

_______________

Zinc (Zn)

(Figure 7)

The zinc concentrations in muscle and liver of reindeer show no significant log-linear change during the period 1983-2005 (parametric test).

The number of years required to detect an annual change of 5% was 8 and 12 years for muscle and liver respectively, with a power of 1.0 to detect a 5% annual change in both tissues for the full period.

The overall geometric mean value of zinc in muscle and liver was 80.8 and 33.5 µg/g (fresh weight) respectively, for the period 1987-2005 and 1983-2005.

_______________

Mercury (Hg)

(Figure 8)

The mercury concentrations in muscle of reindeer show a significant log-linear increase during the period 1987-2005 (9.2% per year, p<0.002; parametric test).

The ANOVA test showed that the smoothed line for concentrations of mercury in liver indicates a significant non-linear trend component (p< 0.059).

(14)

The overall geometric mean value of mercury in muscle and liver was 0.001 and 0.044 µg/g (fresh weight) respectively, for the period 1987-2005 and 1983-2005.

_______________

Lead (Pb)

(Figure 8)

The lead concentrations in liver of reindeer show a significant log-linear change during the period (p<0.093, parametric test), while the concentrations in muscle show no significant log- linear trend during the period (parametric test). The average annual decrease in liver was 2.5%.

The overall geometric mean value of lead in muscle and liver was 0.013 and 0.117 µg/g (fresh weight) respectively, for the period 1987-2002 and 1983-2005.

_______________

Molybdenum (Mo)

(Figure 9)

The molybdenum concentrations in muscle and liver of reindeer show no significant log- linear change during the period (parametric test).

The ANOVA test showed that the smoothed line for concentrations of molybdenum in muscle indicates a significant non-linear trend component (p< 0.001).

The overall geometric mean value of molybdenum in muscle and liver was 0.003 and 0.628 µg/g (fresh weight) respectively, for the period 1997-2005 and 1994-2005.

(15)

Ran and Gran Sámi Villages

─────────────────────────────────────────────────────────

Aluminium (Al)

(Figure 10)

Aluminium has been analysed in muscle and liver sampled only in 1996 and 1997. Thus no statistics are reported

The geometric mean value of aluminium in muscle and liver was 0.035 and 0.042 µg/g (fresh weight) respectively, for the two only years 1996-1997.

_______________

Calcium (Ca)

(Figure 10)

The calcium concentrations in muscle of reindeer show a significant log-linear change during the period 1996-2005 (p<0.019; parametric test), while no significant log-linear change in liver was shown during the period (parametric test). The average annual decrease in muscle was 3.0%.

The number of years required to detect an annual change of 5% was 8 and 9 years

respectively for muscle and liver, with a power to detect a 5% annual change varying between 0.99 and 0.97 for the full period.

The overall geometric mean value of calcium in muscle and liver was 41.6 and 38.6 µg/g (fresh weight) respectively, for the period 1996-2005.

_______________

Cadmium (Cd)

(Figure 11)

The levels of cadmium in muscle from 1996 were all below the detection limit and were excluded from any calculation of a mean value and statistical calculations.

The cadmium concentrations in muscle and liver of reindeer show no significant log-linear change during the period (parametric test).

The overall geometric mean value of cadmium in muscle and liver was 0.007 and 0.498 µg/g (fresh weight) respectively, for the period 1997-2005 and 1996-2005.

_______________

Cobalt (Co)

(Figure 11)

The cobalt concentrations in muscle and liver of reindeer show no significant log-linear change during the period (parametric test).

The ANOVA test shows a significant non-linear trend component for the concentrations of cobalt in liver (p<0.027).

The number of years required to detect an annual change of 5% was 14 and 13 years for muscle and liver respectively, with a power to detect a 5 % annual change varying between 0.42 and 0.51 for the full period.

The overall geometric mean value of cobalt in muscle and liver was 0.004 and 0.100 µg/g (fresh weight) respectively, for the period 1996-2005.

(16)

_______________

Chromium (Cr)

(Figure 12

The levels of chromium in muscle and liver from 1996, 2004 and 2005 were all below the detection limit and were excluded from any calculation of mean values and statistical calculations.

The chromium concentrations in muscle and liver of reindeer show no significant log-linear change during the period (parametric test).

The number of years required to detect an annual change of 5% was 21 and 16 years for muscle and liver respectively, with a power varying between 0.08 and 0.11 to detect a 5%

annual change in both tissues for the full period.

The overall geometric mean value of chromium in muscle and liver was 0.028 and 0.013 µg/g (fresh weight) respectively, for the period 1997-2003.

_______________

Copper (Cu)

(Figure 12)

The copper concentrations in liver of reindeer show no significant log-linear change during the period (parametric test).

The overall geometric mean value of copper in muscle and liver was 1.35 and 57.9 µg/g (fresh weight) respectively, for the period 1996-2005.

_______________

Iron (Fe)

(Figure 13)

The iron concentrations in muscle of reindeer show a significant log-linear increase during the period (p<0.004; parametric test), while the concentrations in liver show a significant

decrease during the period (p<0.039; parametric test). The average annual increase in muscle was 3.2% and the decrease in liver 3.6%.

The ANOVA test showed that the smoothed line for concentrations of iron in muscle indicates a significant non-linear trend component (p<0.017).

The overall geometric mean value of iron in muscle and liver was 30.1 and 89.0 µg/g (fresh weight) respectively, for the period 1996-2005.

It should be noted that concentration of iron might vary with amount of blood and red blood capsules in the analysed samples.

_______________

Magnesium (Mg)

(Figure 13)

The magnesium concentrations in muscle and liver of reindeer show no significant log-linear change during the period (parametric test).

(17)

The overall geometric mean value of magnesium in muscle and liver was 198 and 179 µg/g (fresh weight) respectively, for the period 1996-2005.

_______________

Manganese (Mn)

(Figure 14)

The manganese concentrations in muscle and liver of reindeer show no significant log-linear change during the period (parametric test).

The ANOVA test showed that the smoothed line for concentrations of manganese in liver indicates a significant non-linear trend component (p<0.016).

The overall geometric mean value of manganese in muscle and liver was 0.148 and 2.55 µg/g (fresh weight), respectively for the period 1996-2005.

_______________

Nickel (Ni)

(Figure 14)

The levels of nickel in muscle from 1996, 2002 – 2005 were all below the detection limit and were excluded from any calculation of mean values and statistical calculations, as well as levels in liver from 1996 and 2001-2005.

The nickel concentrations in liver of reindeer show a significant log-linear change during the period 1997-2000, (average annual decrease 25%, p<0.051; parametric test).

The number of years required to detect an annual change of 5% was 27 and 10 years for muscle and liver respectively, with a power of 0.06 and 0.09 to detect a 5% annual change in both tissues for the full period.

The ANOVA test showed that the smoothed line for concentrations of nickel in muscle and liver indicate significant non-linear trend components (p<0.016; muscle, p<0.016; liver).

The overall geometric mean value of nickel in muscle and liver was 0.020 and 0.017 µg/g (fresh weight), respectively for the period 1997-2001, and 1997-2000.

_______________

Vanadium (V)

(Figure 15)

All levels of vanadium in muscle from the period, except for 2001 and 2005, were below the detection limit. These levels were all excluded of any calculations of mean values and statistical calculations. The levels in liver from 1996-1997, 2002 and 2004 were also below the detection limit.

The vanadium concentrations in liver of reindeer show no significant log-linear change during the period (parametric test).

The number of years required to detect an annual change of 5% was 26 years for liver, with a power to detect a 5% annual change of 0.06 for the full period.

(18)

The overall geometric mean value of vanadium in liver was 0.002 µg/g (fresh weight), for the period 1998-2005.

_______________

Zinc (Zn)

(Figure 15)

The zinc concentrations in muscle and liver of reindeer show no significant log-linear change during the period (parametric test).

The overall geometric mean value of zinc in muscle and liver was 80.7 and 29.4 µg/g (fresh weight) respectively, for the period 1996-2005.

_______________

Mercury (Hg)

(Figure16)

The mercury concentrations in liver of reindeer show a significant log-linear change during the period (p<0.058; parametric test). The annual decrease was 9.7%.

The ANOVA test shows that the smoothed lines for the concentrations of mercury in muscle and liver of reindeer indicate a non-linear trend component (p<0.021, muscle; p<0.068, liver).

The overall geometric mean value of mercury in muscle and liver was 0.001 and 0.029 µg/g (fresh weight) respectively, for the period 1996-2005.

_______________

Lead (Pb)

(Figure 16)

The levels of lead in muscle from 1996-1998, 2002 and 2005 were all below the detection limit and were excluded from any calculation of a mean value and statistical calculations.

The lead concentrations in liver of reindeer show a significant log-linear change during the period (p<0.082; parametric test). The annual decrease was 10%.

The number of years required to detect an annual change of 5% was 26 and 21 years for muscle and liver respectively, with a power between 0.06 and 0.13 to detect a 5% annual change in both tissues for the full period.

The overall geometric mean value of lead in muscle and liver was 0.005 and 0.113 µg/g (fresh weight) respectively, for the period 1999-2004 and 1996-2005.

_______________

(19)

REFERENCES

AMAP. 1993. The Monitoring Programme for Arctic Monitoring And Assessment Programme, AMAP. AMAP Report 93:3. ISBN 82-7655-140-8. Oslo, Norway.

Bignert, A. 1998. Comments Concerning the National Swedish Contaminant Monitoring Programme in Marine Biota. Report to the Swedish Environmental Protection Agency.

Frank, A. 1976. Automated wet ashing and multi-metal determination in biological material by atomic-absorption spectrometry. Z. Anal. Chem. 279, 101-102.

Frank, A. 1983. Alltid viss risk? (Always risky? - wet ashing with perchloric acid.) Kem. Tidskr. 95, 63 (In Swedish).

Frank, A. 1988. Semi-micro accessory to an automated wet digestion system for ashing small sample amounts. In: Trace Element Analytical Chemistry in Medicine and Biology, Vol. 5. Brätter, P.

and Schramel, P. (eds.). Proceedings of the Fifth International Workshop. Walter de Gruyter.

Berlin. New York. p. 78-83.

Frank, A. and Petersson, L.R. 1983. Selection of operating conditions and analytical procedure in multi-element analysis of animal tissue by d.c. plasma-atomic emission spectroscopy.

Spectrochim. Acta 38 B, 207-220.

Frank, A. and Petersson, L.R. 1985. Direct current plasma-atomic emission spectrometer as a simultaneous multi-element tool for analysis of biological materials. Kemia-Kemi 12, 426- 430.

Forberg, S., Odsjö, T. and Olsson, M. 1992. Radiocesium in muscle tissue of reindeer and pike from northern Sweden before and after the Chernobyl accident. A retrospective study on tissue samples from the Swedish Environmental Specimen Bank. The Science of the Total Environment, 115:179-189.

Galgan V. and Frank, A. 1988. Automated system for determination of selenium in biological materials. In: Trace Element Analytical Chemistry in Medicine and Biology, Vol. 5. Brätter, P. and Schramel, P. (eds.). Proceedings of the Fifth International Workshop. Walter de Gruyter. Berlin. New York. p. 84-8

Gilbert, R.O. 1987. Statistical Methods for Environmental Pollution Monitoring. Van Nostrand Reinhold, New York.

Helsel, D.R. & Hirsch, R.M. 1995. Statistical Methods in Water Resources, Studies in Environmental Sciences 49. Elsevier, Amsterdam.

Hoaglin, D.C. and Welsch, R.E. 1978. The hatmatrix in regression and ANOVA. Amer. Stat., 32:17- 22.

Loftis, J.C., Ward, R.C. and Phillips, R.D. 1989. An Evaluation of Trend Detection Techniques for Use in Water Quality Monitoring Programs. EPA/600/3-89/037. 139p.

Nicholson, M.D. & Fryer, R. 1991. The Power of the ICES Cooperative Monitoring Programme to Detect Linear Trends and Incidents. In: Anon. Report of the Working Group on Statistical Aspects of Trend Monitoring. ICES Doc CM 1991.

Nicholson, M.D., Fryer R. and Larsen, J.R. 1995. A Robust Method for Analysing Contaminant Trend Monitoring Data. Techniques in Marine Environmental Sciences. ICES.

(20)

Odsjö, T., Asplund, L., Eriksson, U., Kärsrud, A.-S. and Litzén, K. 1998. Time Trends of HCHs and HCB in Muscle of Reindeer (Rangifer tarandus) from Lapland, Northern Sweden, 1983-1995.

Proceedings from the 18th Symposium on Halogenated Environmental Organic Pollutants, Stockholm, Sweden, August 17-21, 1998. In: DIOXIN-98. Environmental Levels P35. (Eds.) N. Johansson, Å. Bergman, D. Broman, H. Håkansson, B. Jansson, E. Klasson Wehler, L.

Poellinger and B Wahlström. Organohalogen Compounds 39:351-354.

Odsjö, T. and Olsson, M. 1979a. Program och arbetsbeskrivning vid miljöprovsbanken. Rapport till Statens naturvårdsverk, 1979-03-15. (In Swedish).

Odsjö, T. and Olsson, M. 1979b. Förslag till miljögiftsprogram inom Program för övervakning av miljökvalitet, PMK. 1979-06-14. (In Swedish).

Swertz, O. 1995. Trend assessment using the Mann-Kendall test. Report of the Working Group on Statistical Aspects of Trend Monitoring. ICES CM 1995/D:2.

(21)

Figure 1. Sampling areas for reindeer in Sweden

The Swedish Monitoring Programme of Contaminants in Terrestrial Biota Sampling Areas

Reindeer

Sami villages 1) Gabna 2) Laevas 3) Girjas 4) Gran 5) Ran 6) Handlsdalen 7) Mittdalen Summer and winter habitats

1 2

3

4

5

6

7

TISS - 98.07.30 16:40, tjodx

(22)

Figure 2. Al and Ca in reindeer from Gabna, Lævas and Girjas Sámi Villages, N. Lapland. (µg/g fresh weight).

Log-linear regression on geometric means, suspected outliers indicated. Smoother: 3-point running mean, unweighted. Al in muscle

.0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8

81 86 91 96 01 n(tot)=122,n(yrs)=12 m=.143 (.087,.236) slope=-7.1%(-22,7.4) SD(lr)=.78,22%,30 yr power=.10/.08/31%

y(98)=.096 (.038,.248) r2=.11, NS

tao=-.24, NS

SD(sm)=.62, p<.077,30 % slope=-7.6%(-30,15) SD(lr)=.87,36%,32 yr power=.07/.07/36%

r2=.07, NS

Al in liver

.0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8

81 86 91 96 01 n(tot)=160,n(yrs)=16 m=.332 (.214,.514) slope=-6.6%(-16,2.6) SD(lr)=.79,14%,30 yr power=.18/.08/32%

y(98)=.202 (.091,.453) r2=.15, NS

tao=-.15, NS

SD(sm)=.42, p<.002 *,19 % slope=-6.8%(-32,19) SD(lr)=1.0,42%,35 yr power=.07/.07/42%

r2=.04, NS

Ca in muscle

0 20 40 60 80 100 120 140 160

81 86 91 96 01 n(tot)=190,n(yrs)=19 m=38.8 (34.9,43.1) slope=1.0%(-.89,2.9) SD(lr)=.22,2.7%,13 yr power=1.0/.43/7.8%

y(05)=42.5 (34.8,51.9) r2=.07, NS

tao=.11, NS

SD(sm)=.18, n.s.,8.1 % slope=-2.2%(-4.7,.34) SD(lr)=.10,3.6%,9 yr power=.98/.98/3.6%

r2=.33, p<.079

Ca in liver

0 20 40 60 80 100 120 140 160

81 86 91 96 01 n(tot)=230,n(yrs)=23 m=42.3 (40.4,44.4) slope=-.57%(-1.3,.11) SD(lr)=.10,1.0%,9 yr power=1.0/.97/3.8%

y(05)=39.7 (36.4,43.4) r2=.13, p<.094 tao=-.23, NS

SD(sm)=.09, p<.094,3.9 % slope=-.44%(-2.5,1.6) SD(lr)=.08,2.9%,8 yr power=1.0/1.0/2.9%

r2=.03, NS

pia - 07.03.02 10:03, alca2007

Department of Contaminant Research / Swedish Museum of Natural History in Stockholm and Department of Chemistry at the National Veterinary Institute in Uppsala.

(23)

Figure 3. Cd and Co in reindeer from Gabna, Lævas and Girjas Sámi Villages, N. Lapland. (µg/g fresh weight).

Log-linear regression on geometric means, suspected outliers indicated. Smoother: 3-point running mean, unweighted.

Cd in muscle

.00 .01 .02 .03 .04 .05 .06 .07 .08 .09

81 86 91 96 01 n(tot)=176,n(yrs)=17 m=.010 (.007,.016) slope=-8.9%(-15,-2.7) SD(lr)=.67,10%,27 yr power=.27/.09/26%

y(05)=.005 (.002,.009) r2=.38, p<.008 * tao=-.46, p<.011 *

SD(sm)=.53, p<.074,25 % slope=-9.5%(-29,9.7) SD(lr)=.63,30%,26 yr power=.08/.09/25%

r2=.16, NS

Cd in liver

.0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8

81 86 91 96 01 n(tot)=230,n(yrs)=23 m=.440 (.400,.484) slope=.78%(-.65,2.2) SD(lr)=.22,2.1%,14 yr power=1.0/.42/8.0%

y(05)=.479 (.399,.576) r2=.06, NS

tao=.08, NS

SD(sm)=.21, n.s.,9.3 % slope=-1.8%(-8.8,5.3) SD(lr)=.28,10%,16 yr power=.28/.28/10%

r2=.04, NS

Co in muscle

.00 .01 .02 .03 .04 .05 .06 .07 .08

81 86 91 96 01 n(tot)=193,n(yrs)=19 m=.007 (.004,.013) slope=4.2%(-8.0,16) SD(lr)=1.4,19%,43 yr power=.12/.06/62%

y(05)=.010 (.003,.036) r2=.03, NS

tao=-.19, NS

SD(sm)=1.5, n.s.,84 % slope=-9.4%(-15,-4.2) SD(lr)=.21,7.5%,13 yr power=.47/.47/7.5%

r2=.68, p<.003 *

Co in liver

.0 .1 .2 .3 .4 .5 .6 .7 .8

81 86 91 96 01 n(tot)=230,n(yrs)=23 m=.141 (.127,.158) slope=-1.9%(-3.4,-.54) SD(lr)=.22,2.0%,13 yr power=1.0/.43/7.9%

y(05)=.114 (.095,.137) r2=.28, p<.009 * tao=-.36, p<.016 *

SD(sm)=.18, p<.060,7.8 % slope=-3.0%(-7.5,1.5) SD(lr)=.18,6.4%,12 yr power=.59/.59/6.4%

r2=.23, NS

pia - 07.03.02 10:05, CdCo2007

(24)

Figure 4. Cr and Cu in reindeer from Gabna, Lævas and Girjas Sámi Villages, N. Lapland. (µg/g fresh weight).

Cr in muscle

.00 .04 .08 .12 .16 .20 .24 .28 .32 .36 .40 .44 .48

y(03)=.030 (.009,.094) r2=.05, NS

tao=.20, NS

SD(sm)=1.0, n.s.,53 % slope=27%(7.1,46) SD(lr)=.72,35%,28 yr power=.07/.08/29%

r2=.60, p<.014 *

Cr in liver

.00 .02 .04 .06 .08 .10 .12 .14 .16 .18

81 86 91 96 01 n(tot)=137,n(yrs)=13 m=.014 (.009,.022) slope=2.4%(-3.6,8.4) SD(lr)=.70,17%,28 yr power=.13/.08/28%

y(03)=.018 (.009,.036) r2=.07, NS

tao=.08, NS

SD(sm)=.26, p<.001 *,11 % slope=1.1%(-13,15)

SD(lr)=.29,21%,16 yr power=.11/.27/11%

r2=.01, NS

Cu in muscle

0 1 2 3 4 5 6 7 8

81 86 91 96 01 n(tot)=188,n(yrs)=19 m=1.32 (1.27,1.37) slope=-.59%(-1.3,.08) SD(lr)=.08,1.0%,8 yr power=1.0/1.0/2.7%

y(05)=1.25 (1.17,1.34) r2=.17, p<.076 tao=-.26, NS

SD(sm)=.07, n.s.,3.0 % slope=-.71%(-2.9,1.5) SD(lr)=.09,3.1%,8 yr power=1.0/1.0/3.1%

r2=.07, NS

Cu in liver

0 50 100 150 200 250

81 86 91 96 01 n(tot)=230,n(yrs)=23 m=67.2 (59.4,76.1) slope=1.7%(-.05,3.5) SD(lr)=.27,2.5%,15 yr power=1.0/.29/9.9%

y(05)= 81 ( 65, 102) r2=.16, p<.054

tao=.26, p<.077 SD(sm)=.26, n.s.,12 % slope=-4.7%(-11,1.7) SD(lr)=.25,9.2%,15 yr power=.33/.33/9.2%

r2=.26, NS

pia - 07.03.02 10:07, CrCu2007

(25)

Figure 5. Fe and Mg in reindeer from Gabna, Lævas and Girjas Sámi Villages, N. Lapland. (µg/g fresh weight).

Fe in muscle

0 25 50 75 100 125 150 175 200

81 86 91 96 01 n(tot)=190,n(yrs)=19 m=29.4 (28.1,30.7) slope=.02%(-.82,.85) SD(lr)=.09,1.2%,8 yr power=1.0/.99/3.4%

y(05)=29.4 (26.9,32.1) r2=.00, NS

tao=-.02, NS

SD(sm)=.07, p<.043 *,3.2 % slope=1.3%(-.76,3.4) SD(lr)=.08,3.0%,8 yr power=1.0/1.0/3.0%

r2=.21, NS

Fe in liver

0 100 200 300 400 500 600 700 800 900 1000 1100 1200

81 86 91 96 01 n(tot)=230,n(yrs)=23 m=140 (120 ,162 ) slope=-.74%(-3.1,1.6) SD(lr)=.35,3.3%,18 yr power=.99/.19/13%

y(05)= 129 ( 96, 173) r2=.02, NS

tao=-.19, NS

SD(sm)=.31, n.s.,14 % slope=-7.6%(-19,3.4) SD(lr)=.44,16%,21 yr power=.14/.14/16%

r2=.24, NS

Mg in muscle

0 100 200 300 400 500 600 700 800 900

81 86 91 96 01 n(tot)=190,n(yrs)=19 m=202 (190 ,214 ) slope=-.38%(-1.4,.69) SD(lr)=.12,1.5%,10 yr power=1.0/.91/4.3%

y(05)= 195 ( 174, 218) r2=.03, NS

tao=-.09, NS

SD(sm)=.09, p<.015 *,3.7 % slope=1.0%(-.56,2.6) SD(lr)=.06,2.2%,7 yr power=1.0/1.0/2.2%

r2=.22, NS

Mg in liver

0 100 200 300 400 500 600 700 800 900

81 86 91 96 01 n(tot)=230,n(yrs)=23 m=183 (178 ,189 ) slope=-.36%(-.78,.06) SD(lr)=.06,.60%,7 yr power=1.0/1.0/2.3%

y(05)= 176 ( 167, 186) r2=.13, p<.088

tao=-.26, p<.086

SD(sm)=.05, p<.040 *,2.2 % slope=-.82%(-2.6,1.0) SD(lr)=.07,2.6%,7 yr power=1.0/1.0/2.6%

r2=.12, NS

pia - 07.03.02 10:09, FeMg2007

(26)

Figure 6. Mn and Ni in reindeer from Gabna, Lævas and Girjas Sámi Villages, N. Lapland. (µg/g fresh weight).

Mn in muscle

.0 .4 .8 1.2 1.6 2.0 2.4

81 86 91 96 01 n(tot)=189,n(yrs)=19 m=.146 (.135,.158) slope=-1.3%(-2.6,.03) SD(lr)=.15,1.9%,11 yr power=1.0/.74/5.5%

y(05)=.130 (.113,.150) r2=.20, p<.053 tao=-.43, p<.011 *

SD(sm)=.13, p<.088,5.5 % slope=.22%(-3.7,4.2) SD(lr)=.15,5.6%,11 yr power=.71/.71/5.6%

r2=.00, NS

Mn in liver

0 1 2 3 4 5 6 7 8 9 10 11

81 86 91 96 01 n(tot)=230,n(yrs)=23 m=2.72 (2.39,3.10) slope=-2.2%(-4.0,-.46) SD(lr)=.27,2.5%,15 yr power=1.0/.30/9.8%

y(05)=2.13 (1.70,2.67) r2=.25, p<.015 * tao=-.39, p<.009 *

SD(sm)=.17, p<.001 *,7.2 % slope=-3.6%(-10,3.0) SD(lr)=.26,9.5%,15 yr power=.31/.31/9.5%

r2=.16, NS

Ni in muscle

.00 .02 .04 .06 .08 .10 .12 .14

y(05)=.028 (.005,.146) r2=.13, NS

tao=-.05, NS

SD(sm)=.74, p<.047 *,36 % slope=.35%(-22,22)

SD(lr)=.50,39%,23 yr power=.07/.12/19%

r2=.00, NS

Ni in liver

.00 .02 .04 .06 .08 .10 .12 .14

81 86 91 96 01 n(tot)=125,n(yrs)=12 m=.015 (.010,.022) slope=-1.3%(-7.6,5.1) SD(lr)=.64,18%,26 yr power=.12/.09/25%

y(05)=.013 (.005,.032) r2=.02, NS

tao=-.21, NS

SD(sm)=.36, p<.019 *,16 % slope=-3.0%(-32,26) SD(lr)=.29,47%,16 yr power=.07/.26/11%

r2=.04, NS

pia - 07.03.02 10:12, MnNi2007

(27)

Figure 7. V and Zn in reindeer from Gabna, Lævas and Girjas Sámi Villages, N. Lapland. (µg/g fresh weight).

V in muscle

.00 .01 .02 .03 .04 .05 .06 .07 .08 .09

81 86 91 96 01 n(tot)=94,n(yrs)=9 m=.005 (.002,.009) slope=-9.1%(-18,-.01) SD(lr)=.69,34%,28 yr power=.07/.08/27%

y(05)=.002 (.001,.005) r2=.44, p<.049 * tao=-.33, NS

SD(sm)=.56, n.s.,26 % SD(lr)=.98,>99%,34 yr r2=.07, NS

V in liver

.00 .01 .02 .03 .04 .05 .06 .07 .08 .09

81 86 91 96 01 n(tot)=143,n(yrs)=14 m=.007 (.004,.011) slope=-8.2%(-16,-.55) SD(lr)=.75,16%,29 yr power=.14/.08/30%

y(01)=.003 (.001,.007) r2=.31, p<.036 * tao=-.16, NS

SD(sm)=.39, p<.003 *,18 % slope=-22%(-34,-10) SD(lr)=.31,51%,17 yr power=.07/.23/11%

r2=.92, p<.012 *

Zn in muscle

0 20 40 60 80 100 120 140 160 180

81 86 91 96 01 n(tot)=190,n(yrs)=19 m=80.8 (77.8,84.0) slope= .00%(-.72,.72) SD(lr)=.08,1.1%,8 yr power=1.0/1.0/2.9%

y(05)=80.8 (74.9,87.2) r2=.00, NS

tao=.10, NS

SD(sm)=.08, n.s.,3.5 % slope=.79%(-.40,2.0) SD(lr)=.05,1.7%,6 yr power=1.0/1.0/1.7%

r2=.23, NS

Zn in liver

0 20 40 60 80 100 120 140 160 180

81 86 91 96 01 n(tot)=230,n(yrs)=23 m=33.5 (31.1,36.0) slope=.04%(-1.1,1.2) SD(lr)=.18,1.7%,12 yr power=1.0/.60/6.4%

y(05)=33.6 (29.0,38.9) r2=.00, NS

tao=-.01, NS

SD(sm)=.18, n.s.,7.9 % slope=-3.9%(-9.1,1.4) SD(lr)=.21,7.5%,13 yr power=.46/.46/7.5%

r2=.26, NS