Heavy-metal concentrations in small mammals from a diffusely polluted floodplain: importance of species- and location-specific characteristics.

The soil of several floodplain areas along large European rivers shows increased levels of heavy metals as a relict from past sedimentation of contaminants. These levels may pose risks of accumulation in food webs and toxicologic effects on flora and fauna. However, for floodplains, data on heavy-metal concentrations in vertebrates are scarce. Moreover, these environments are characterised by periodical flooding cycles influencing ecologic processes and patterns. To investigate whether the suggested differences in accumulation risks for insectivores and carnivores, omnivores, and herbivores are reflected in the actual heavy-metal concentrations in the species, we measured the current levels of Zn, Cu, Pb, and Cd in 199 specimens of 7 small mammal species (voles, mice, and shrews) and in their habitats in a diffusely polluted floodplain. The highest metal concentrations were found in the insectivorous and carnivorous shrew, Sorex araneus. Significant differences between the other shrew species, Crocidura russula, and the vole and mouse species was only found for Cd. The Cu concentration in Clethrionomys glareolus, however, was significantly higher than in several other vole and mouse species. To explain the metal concentrations found in the specimens, we related them to environmental variables at the trapping locations and to certain characteristics of the mammals. Variables taken into account were soil total and CaCl(2)-extractable metal concentrations at the trapping locations; whether locations were flooded or nonflooded; the trapping season; and the life stage; sex; and fresh weight of the specimens. Correlations between body and soil concentrations and location or specimen characteristics were weak. Therefore; we assumed that exposure of small mammals to heavy-metal contamination in floodplains is significantly influenced by exposure time, which is age related, as well as by dispersal and changes in foraging and feeding patterns under influence of periodic flooding.

Industrial and communal wastewater discharges and agricultural activities have caused large-scale soil contamination of a majority of the floodplains along the large European rivers (Nienhuis et al. 1998). The heavy metals Cd, Cu, Pb, and Zn are present in large amounts in these diffusely polluted floodplains (Middelkoop & Van Haselen 1999; Vink et al. 1999). Several studies have suggested that the present contaminant levels pose risks to floodplain ecosystems through accumulation of heavy metals in food webs and possible toxicologic effects in a variety of species (Hendriks et al. 1995; Van den Brink et al. 2003; Kooistra et al. 2001, 2005; Leuven et al. 2005).

Studies of contaminant levels in floodplain species have been scarce and generally focused on lower trophic levels, such as vegetation (Schröder 2005), and macro-invertebrates, such as earthworms, snails, spiders, and insects (Hobbelen et al. 2004; Notten et al. 2005; Van Vliet et al. 2005). Exceptions are studies including the common shrew Sorex araneus (Hendriks et al. 1995), the little owl Athene noctua vidalli (Van den Brink et al. 2003), and the badger Meles meles, which were expected to forage in floodplains (Van den Brink & Ma 1998). Assessments of the risk of contaminant accumulation in vertebrates (e.g., mammals and birds) in floodplains have generally been based on soil contaminant levels combined with accumulation factors, consumption rates, and life expectancies (Jongbloed et al. 1996; Pascoe et al. 1996; Kooistra et al. 2001, 2005). Except for the soil contaminant levels, these parameters for the calculation of accumulation risks are generally derived from literature data from inland area or laboratory studies. However, floodplains are highly dynamic environments, with periodical flooding affecting species distribution, life expectancy and mortality, habitat suitability patterns within the landscape, food availability, and recolonisation processes from the nonflooded areas (Robinson et al. 2002; Klok et al. 2006; Wijnhoven et al. 2006).

Small mammals (voles, mice, and shrews) play an important role in floodplain food webs by acting at different trophic levels. They include predominantly herbivorous as well as insectivorous and carnivorous species. They are prey to a whole range of predatory mammals and birds of prey (Erlinge et al. 1983; Jongbloed et al. 1996). Because these small-mammal species can be numerous in certain areas within the floodplains and are often mentioned as species at risk of toxicologic effects of the current contaminant levels in floodplains themselves, they are suitable as monitors of contaminant levels in floodplain ecosystems. It has been shown that the distribution and densities of the common small-mammal species are highly influenced by flooding events, resulting in the highest densities throughout the year occurring on and near the nonflooded areas (Wijnhoven et al. 2005, 2006). Densities were also found to be much higher at the end of summer and in autumn compared with winter and spring. This is a result of the gradual population growth after increased mortality in winter, especially during floods, which is reflected in the age distribution. Total metal concentrations in the soil are generally higher in the lower, frequently flooded areas than in the nonflooded areas, but CaCl2-extractable concentrations do not show a similar pattern (Wijnhoven et al. 2006). It is assumed that metal concentrations in small mammals are reflected not only by trophic level but also by exposure time and metal concentrations in the exposure areas (Hunter et al. 1989; Torres & Johnson 2001). Furthermore, different contaminant levels in small-mammal species can occur because of seasonal variation (Greville & Morgan 1989; Ma et al. 1991) as well as differences in size and sex (Dodds-Smith et al. 1992; Damek-Poprawa & Sawicka-Kapusta 2004). We hypothesized that of the seven common small mammal species in our research area (Wijnhoven et al. 2005), the highest average metal concentrations would be found in the insectivores and carnivores, (e.g., common shrew Sorex araneus and white-toothed shrew Crocidura russula), whereas the lowest ones were expected in predominantly herbivorous species (e.g., common vole Microtus arvalis, short-tailed field vole Microtus agrestis, wood mouse Apodemus sylvaticus, and harvest mouse Micromys minutus), with the more omnivorous bank vole Clethrionomys glareolus in between. Species-specific metal concentrations were expected to be positively related to the metal concentrations in the soil of the trapping locations. We investigated whether these patterns of species- and location-related body metal concentrations were present in a diffusely polluted floodplain. The study tried to answer the following research questions: (1) Are there differences in average metal concentrations between species, and can these differences be explained by feeding behaviour (herbivory, omnivory, and/or insectivory and carnivory)? (2) Are such differences in metal concentration between species similar for each of the investigated metals? (3) Are the possible interspecific and intraspecific differences in body metal concentrations related to soil metal concentrations or CaCl2-extractable concentrations at the trapping locations? (4) Are there intraspecific differences in metal concentrations that may be related to sex, life stage or size of the animals, the time of the year, or the positioning the traps? The article discusses the consequences of our findings for heavy-metal exposure and accumulation risks in floodplain food webs.

Materials and Methods

Data Collection

All data were collected at the Afferdensche en Deestsche Waarden (ADW), a moderately to heavily polluted embanked floodplain area along the river Waal, the main distributary of the Rhine in The Netherlands (Fig. 1). The research area consists of lands inside and outside the summer dikes. The summer dikes are the lower inner embankments protecting agricultural areas in the floodplain against summer floods. Large parts of the floodplain are periodically flooded, on average once a year, predominantly in winter. The floodplain includes areas with and without agricultural activities. Those without agriculture, feature naturally developed vegetation and offer a wide range of habitats. Detailed descriptions of the research area are given in Wijnhoven et al. (2005, 2006).

Fig. 1.

Location of the ADW floodplain area in The Netherlands and positions of the trapping locations (with 5 to 10 traps each) within the floodplain

Small mammals from the ADW floodplain were collected at 58 sites (Fig. 1) between 2001 and 2003 using Longworth live traps in lines of 5 to 10 traps at each site. The traps were baited with apple, carrot, and rinsed meat and were stuffed with hay and tissue. Specifically for this study, there were three sessions of two 3-day trapping trials, all traps were checked twice a day in August 2002 and again in June and October 2003. Furthermore, all trapping casualties from monitoring studies were included, especially specimens trapped in winter and spring because mortality in that time of the year is higher (Wijnhoven et al. 2005, 2006). Trapping locations were originally selected to monitor recolonisation of the floodplain after flood events (Wijnhoven et al. 2005), so they were chosen based on habitat characteristics (vegetation structure, soil type, and management type) without previous information on the levels of contamination. Therefore, trapping sites covered the whole range from nonflooded parts to flooded locations situated far from the nonflooded areas, and were expected to show a representative variation in contaminant levels for the study area. At each of the sites, three soil cores, prepared from three or five soil samples from a 1 m2 plot, were taken with line intervals of at least 10 m. A 5-g portion of soil from each sample was oven dried for 24 hours at 105°C. The total metal content of 0.2 mg dry weight (DW) substrate in a mixture of 3.0 ml 65% HNO3 and 1.5 ml 37% HCl was measured after microwave destruction using a MLS-1200 MEGA microwave oven (Milestone, Sorisole, Italy). The samples were topped up to 50 ml, after which the metal content was measured using inductively coupled plasma–atomic emission spectrometry (ICP-AES; Spectro Analytical Instruments, Kleve, Germany). The 0.01 M CaCl2-exchangeable fraction was determined as a measure of the potential metal solubility. A 6-g fresh-weight portion of substrate, to which 0.01 M CaCl2 had been added in a 1:10 (m[DW]/v) ratio, was mixed for 2 hours, after which the suspension was centrifuged at 12,000 rpm (5,000 x g) for 15 minutes. After the pHCaCl2 had been measured in the substrate suspension in 0.01 M CaCl2, the supernatant was filtered over a 0.45-μm pore filter. A pH of 2 was obtained by adding a few droplets of 65% HNO3, and the metal content of the sample was subsequently measured on the spectrometer.

Fresh weights (FWs) of the mammals were determined, and the liver and kidneys of each specimen were weighed (FW). Parts of these organs and flank muscles were oven dried for 24 hours at 105°C, after which DWs were measured. The metal contents of the animal tissues were measured after microwave extraction of approximately 0.01 to 0.25 mg DW with HNO3 and HCl and analysed using ICP-AES as described previously. The metal concentrations in whole animals were calculated from the concentrations in the liver, kidneys, and muscle tissue. Concentrations in muscle tissue were assumed to reflect the concentration in the animal’s remaining tissues, i.e., the entire animal minus liver and kidneys. Calculations were based on the species-specific weighed average distributions in percentage DW for liver, kidneys, and other tissues (tissue-to-total body ratios) and the tissue-specific DW-to-FW ratios derived from our own data.

Statistical Analysis

Interspecific differences in metal concentration distributions were statistically tested at P < 0.05 using the two-sample Kolmogorov–Smirnov test (KS test) in Systat for Windows 11 (Systat Software Inc., Richmond, CA, USA) because the metal content values within the species did not show a normal distribution (per one-sample KS test). To search for possible explanations for interspecific differences in metal concentrations, metal concentrations in the soil (average total and CaCl2-extractable concentrations at the locations where individuals were trapped) were also compared with the two-sample KS test. Two-sample KS tests are used because concentration distribution variation (“shape”) and average or median concentration distribution (“location”) will affect exposure of populations (Sokal & Rohlf 1995). To calculate functional relations of metal concentrations in species, regressions between metal concentrations and mammal characteristics and environmental factors were calculated making use of Microsoft Office Excel 2003 (Microsoft Corporation, Redmond, WA, USA), according to the stepwise-method, with a critical F value for the regression equations calculated at P < 0.05. A significance level of P < 0.05 for the individual regression coefficients was used. Parameters included were the total and CaCl2-extractable metal concentrations at the trapping sites; the trapping locations in either flooded or nonflooded areas; the trapping times (for which the dates were divided into the four seasons); and the sexes (male vs. female), life stages (juvenile vs. adult); and sizes and conditions of the animals (recorded as mg FW). Data of specimens were only included if all of the information for the selected parameters was available, which means that the sample size in these analyses was sometimes slightly smaller than the total number of specimens collected. Similarity in sex ratio for each of the species and the numbers of a species trapped in flooded and nonflooded areas was tested with the binomial test (with a two-tailed significance level of 5%) (Sokal & Rohlf 1995). Interspecific differences in the distribution of DW-to-FW ratios, organ-to-total body ratios, sex ratio, and life-stage composition of the populations, as well as the proportional distribution of the populations over the flooded and nonflooded areas and the numbers trapped in each season, were also tested in Systat 11 (two-sample KS test) after checking for possible normal distributions (one-sample KS test).

Results

Differences in Heavy-Metal Concentrations Between Species

In total, 199 specimens of 7 small-mammal species were collected (Table 1). The variance in organ-to-total body ratio within the species groups was not large, with the exception of S. araneus (Table 2). Comparing these ratios between species, the two shrew species showed heavier livers and kidneys than the other species. Distributions in DW-to-FW ratios of especially liver, but also muscle tissue, also differed between species. Differences for this ratio in kidneys were found only between A. sylvaticus and M. agrestis. In all species, the DW-to-FW ratios were somewhat higher for livers than for kidneys and lowest for muscle tissue.

Table 1.

Characteristics of the small mammals collected

Species	N	FWbody (g)	Sex ratio		Life stage	Flooding		Season
Apodemus sylvaticus	21	16.06 ± 2.42	1.30 ± 0.47	More male animals	1.52 ± 0.51	1.33 ± 0.48a	NS	2.95 ± 0.50a
Clethrionomys glareolus	56	19.68 ± 4.62	1.55 ± 0.50	NS	1.68 ± 0.47	1.79 ± 0.41b	More nonflooded	2.32 ± 0.64bd
Crocidura russula	11	10.97 ± 1.04	1.45 ± 0.52	NS	2	1ac	More flooded	2.82 ± 0.40ab
Microtus agrestis	9	21.58 ± 8.72	1.33 ± 0.50	NS	1.67 ± 0.50	2bde	More nonflooded	2.44 ± 0.53bc
Microtus arvalis	31	19.48 ± 6.03	1.55 ± 0.51	NS	1.71 ± 0.46	1.06 ± 0.25cf	More flooded	2.61 ± 0.50b
Micromys minutus	4	4.53 ± 0.18	1	NS	2	1adf	NS	3ab
Sorex araneus	67	8.32 ± 2.07	1.45 ± 0.50	NS	1.66 ± 0.48	1.55 ± 0.50ae	NS	2.18 ± 0.85cd

a The numbers per species are shown with average values ( ± SD) of fresh weights’ sex ratios (1 = male; 2 = female), life stages (1 = juvenile; 2 = adult), trapping locations (1 = flooded; 2 = nonflooded) and trapping seasons (1 = spring; 2 = summer; 3 = autumn; 4 = winter). Unequal intraspecies distributions of sex ratio and flooding were tested (P < 0.05) using binomial test (indicated when unequal; NS = no significant unequal distribution). Interspecies differences were tested when relevant using the two-sample Kolmogorov-Smirnov test (P < 0.05); significant differences are indicated by different superscript letters.

Table 2.

Species and tissue-specific DW to FW ratios (± SD) and species-specific tissue-to-total body ratios (± SD) based on DWs used to calculate whole-body metal concentrations of specimens

Species	N	DW-to-FW ratio			Tissue-to-total body ratio (in DW)
		Liver	Kidney	Muscle	Liver	Kidney	Remaining tissue
A. sylvaticus	21	0.316 ± 0.057ac	0.253 ± 0.027a	0.235 ± 0.057a	0.0648 ± 0.0103a	0.0152 ± 0.0027a	0.920 ± 0.011a
C. glareolus	56	0.272 ± 0.081b	0.243 ± 0.043ab	0.218 ± 0.066c	0.0701 ± 0.0123a	0.0165 ± 0.0029a	0.913 ± 0.014c
C. russula	11	0.315 ± 0.041a	0.246 ± 0.022ab	0.204 ± 0.056abc	0.0904 ± 0.0191b	0.0213 ± 0.0031b	0.888 ± 0.021b
M. agrestis	9	0.251 ± 0.038b	0.228 ± 0.034b	0.210 ± 0.044abc	0.0701 ± 0.0115a	0.0149 ± 0.0027ac	0.915 ± 0.012ac
M. arvalis	31	0.287 ± 0.034cd	0.253 ± 0.029ab	0.235 ± 0.033ac	0.0678 ± 0.0097a	0.0127 ± 0.0013c	0.919 ± 0.010ac
M. minutus	4	0.264 ± 0.012bd	0.244 ± 0.032ab	0.173 ± 0.088abc	0.0679 ± 0.0104a	0.0250 ± 0.0060b	0.907 ± 0.014abc
S. araneus	67	0.305 ± 0.056ad	0.239 ± 0.046ab	0.178 ± 0.072b	0.106 ± 0.084b	0.0241 ± 0.0239b	0.869 ± 0.107b

a Interspecies differences were tested using the two-sample Kolmogorov-Smirnov test (P < 0.05); significant differences are indicated by different superscript letters

The average whole-body Zn concentration was highest in S. araneus with 126 mg kg–1 DW, followed in descending order by C. glareolus, C. russula, A. sylvaticus, M. arvalis, M. agrestis, and M. minutus (Fig. 2). The differences between S. araneus and M. arvalis, between S. araneus and A. sylvaticus and between S. araneus and C. glareolus were significant (P < 0.05). The variation in Zn concentrations in S. araneus was large. The highest Zn concentration was, however, observed in a specimen of C. glareolus, which contained 866 mg kg–1 DW, more than twice as much as the highest concentration measured in a specimen of S. araneus (416 mg kg–1 DW).

Fig. 2.

Metal concentrations in various small-mammal species trapped in the ADW floodplain. Average metal concentrations are indicated by crosses; error bars show SDs; total range (minimum and maximum observations) is shown as a vertical line. The columns indicate the 25% and 75% percentiles, with the median in between. Different letters indicate significant differences (P < 0.05) in concentration distributions; identical letters indicate no significant differences between species

Average Cu concentrations were also highest in S. araneus, followed in descending order by M. arvalis, C. glareolus, C. russula, M. minutus, A. sylvaticus, and (lowest) M. agrestis. The average of 15.6 mg kg–1 DW in S. araneus was significantly higher than in all other species except C. russula and M. minutus. The average Cu concentration of 6.71 mg kg–1 DW in C. glareolus was also significantly higher than those found in A. sylvaticus and M. agrestis, with 4.05 and 3.38 mg kg–1 DW, respectively. The highest Cu concentration was measured in a specimen of S. araneus, with 88.9 mg kg–1 DW, but an individual of M. arvalis also had a high concentration, i.e., 62.0 mg kg–1 DW.

The order from high to low in average Pb concentrations was S. araneus, C. russula, C. glareolus, M. agrestis, M. arvalis, A. sylvaticus, and M. minutus. The average value of 32.6 mg Pb kg–1 DW in S. araneus was significantly higher than the average concentrations in A. sylvaticus and C. glareolus. For all species, only a few individuals had increased Pb levels, as was indicated by the median concentration of < 1.60 mg kg–1 DW for all species. The highest concentrations in individuals were found in C. glareolus (270 mg Pb kg–1) and S. araneus (264 mg Pb kg–1).

The average Cd concentration was significantly higher in M. agrestis than in S. araneus, C. russula, and M. arvalis; however, this is largely the result of one individual, which contained 63.9 mg Cd kg–1 DW. The average Cd concentration of 3.49 mg kg–1 in S. araneus was significantly higher than the average concentrations in all other species except C. russula. The median concentrations were highest in S. araneus, with 1.82 mg Cd kg–1 DW, followed by C. russula and M. agrestis, with median concentrations of 0.58 and 0.57 mg Cd kg–1 DW, respectively, and the other species with median concentrations < 0.20 mg Cd kg–1 DW.

Exposure of Small Mammals to Heavy Metals

Exposure of species to total metal concentrations showed similar patterns for each of the metals. The total metal concentrations at the trapping sites of C. glareolus and M. agrestis were significantly lower than those at the sites of most of the other species, depending on the investigated metal (Fig. 3). Total Zn concentrations at the trapping locations of A. sylvaticus (398 mg kg–1 DW), C. russula (392 mg kg–1 DW), and M. arvalis (377 mg kg–1 DW) were significantly higher than those at the trapping locations of C. glareolus (200 mg kg–1 DW) and M. agrestis (110 mg kg–1 DW) and higher at the trapping locations of S. araneus (289 mg kg–1 DW) than at those of M. agrestis. For Cu and Pb, significant differences were found between A. sylvaticus, M. arvalis, S. araneus, and C. russula, all of which had high exposure concentrations, and C. glareolus and M. agrestis, which had low exposure concentrations. With regard to Cd, A. sylvaticus and S. araneus were exposed to significantly higher concentrations than were C. glareolus and M. agrestis.

Fig. 3.

Total and CaCl2-extractable metal concentrations in the upper 10 cm of soil at the ADW floodplain area trapping locations. Average metal concentrations are indicated by a cross; error bars show SDs; total range (minimum and maximum observations) is shown as a vertical line. The columns indicate the 25% and 75% percentiles, with the median in between. Different letters indicate significant differences (P < 0.05) in concentration distributions; identical letters indicate no significant differences between metal concentrations in the soil

Exposure to CaCl2-extractable Zn was highest for M. agrestis and C. glareolus, which is significantly higher than for the other species (Fig. 3). Of these species, A. sylvaticus was exposed to significantly higher concentrations than M. arvalis, C. russula, and M. minutus. Differences in CaCl2-extractability of Cu and Pb were smaller, but A. sylvaticus, C. glareolus, M. arvalis, and S. araneus were exposed to higher Cu levels than M. agrestis, whereas A. sylvaticus, C. glareolus, M. arvalis, and S. araneus were exposed to higher Pb levels than C. russula and M. agrestis. Exposure to CaCl2-extractable Cd concentrations was similar for all species.

Regression of Metal Concentrations

Significant regressions between metal concentrations in the various species on one hand and species characteristics and/or environmental variables on the other hand were found for A. sylvaticus with respect to Zn, Cu, and Pb, respectively (Table 3). Of the variance in Zn concentrations in this species, 58% was explained by the parameters taken into account: season, total Zn concentration in the soil, locations being flooded or nonflooded, and CaCl2-extractable Zn concentration in the soil. Cu concentrations were related to life stages, with higher concentrations in adults. Pb concentrations were related to the sizes of the animals, with the highest concentrations in the smallest specimens. Significant regressions of metal concentrations in C. glareolus were found for Zn, with the highest concentrations found in the adults trapped in spring at locations with high CaCl2-extractable Zn concentrations, and for Cd, with the highest concentrations found in the larger animals trapped toward winter. Cu concentrations in C. russula were positively related to CaCl2-extractable concentrations and negatively related to total Cu concentrations in the soil; they were highest in the larger female animals. Cd concentrations in this species were more closely related to the season (highest toward winter) and were also positively related to CaCl2-extractable concentrations in soil; again, they were highest in female individuals. A significant regression (P < 0.05) for Zn concentrations in M. arvalis was only found with the sexes, with the highest concentrations in female animals.

Table 3.

Significant regressions (P < 0.05) between metal concentrations in small-mammal species, characteristics of the species (see Table 1), and environmental variables (critical F value calculated at P < 0.05)

Metal		F	n	R2	P
	A. sylvaticus
Zn	ln([Znbody]) = 1.41(season) + 0.00107([Zntot]) + 1.11(flood) – 0.722([ZnCaCl2]) – 1.24	4.51	18	0.581	0.017
Cu	ln([Cubody]) = 0.918(life stage) – 0.490	6.01	18	0.273	0.026
Pb	ln([Pbbody]) = −0.000610(FWbody) + 7.30	4.94	18	0.236	0.041
Cd	NS		18
	C. glareolus
Zn	ln([Znbody]) = 0.693(life stage) – 0.554(season) + 0. 295([ZnCaCl2]) +6.58	5.95	49	0.284	0.002
Cu	NS		49
Pb	NS		49
Cd	ln([Cdbody]) = 0.000236(FWbody) + 1.48(season) –10.5	6.49	49	0.220	0.003
	C. russula(s2, f1)
Zn	NS		10
Cu	ln([Cubody]) = 45.1([Cu]CaCl2) – 0.0187([Cutot]) + 0.828(sex) + 0.000247(FWbody) – 3.57	19.9	10	0.941	0.003
Pb	NS		10
Cd	ln([Cdbody]) = 4.87(season) +199([Cd]CaCl2) +1.52(sex) –19.2	17.2	10	0.896	0.002
	M. agrestis(f2)
Zn	NS		8
Cu	NS		8
Pb	NS		8
Cd	NS		8
	M. arvalis
Zn	ln([Znbody]) = 0.342(sex) + 3.73	4.65	29	0.147	0.040
Cu	NS		29
Pb	NS		29
Cd	NS		29
	S. araneus
Zn	NS		45
Cu	NS		45
Pb	NS		45
Cd	NS		45

M. minutus is not included because only four animals were trapped. [Mebody] = metal concentration in animal (mg kg−1 DW). Species characteristics and environmental variables initially included are sex = 1 for male and 2 for female; life stage = 1 for juvenile and 2 for adult; FWbody = fresh weight of animal (mg); [Metot] = total metal concentration in the soil (mg kg−1 DW); [MeCaCl2] = CaCl2-extractable metal concentration in the soil (mg kg−1 DW); flooding = 1 for flooded and 2 for nonflooded trapping locations; season = 1 for animals trapped in spring, 2 for summer, 3 for autumn, 4 for winter; s2 = all animals were adults; f1 = all animals were trapped in flooded areas; f2 = all animals were trapped in nonflooded areas. NS = nonsignificant.

Differences in Population Structure

Because differences in metal concentrations between species may result from differences in population structure or distribution, the characteristics of each species were analysed. For all species, more adults than juveniles were trapped, but for none of the species was the adults-to-juveniles ratio significantly higher than for another species. The fact that only adults were trapped from two species (C. russula and M. minutus) was reflected in the FWs of these species, which showed a relatively narrow range, with an SD <10% of the average weight. Comparable numbers of male and female individuals were trapped for nearly all species, except for A. sylvaticus, for which significantly more (i.e., twice as many) male than female animals were trapped. With regard to seasonal differences, S. araneus and C. glareolus were trapped significantly earlier, because half of those individuals had already been collected in early summer, than A. sylvaticus, of which more individuals were collected throughout fall. S. araneus was trapped earlier than M. arvalis, M. agrestis, C. russula, and M. minutus.

With regard to the numbers trapped in flooded and nonflooded areas, significantly (P < 0.05) more specimens of C. glareolus and M. agrestis were trapped in nonflooded areas, whereas significantly more specimens of C. russula and M. arvalis were trapped in flooded areas.

Discussion

Differences in average metal concentrations between species can be the result of differences in population structure between the species. The data set also showed differences in the relative numbers of each of the species trapped in each of the seasons. Nevertheless, the trapping results were in accordance with the distribution patterns of the species in the ADW floodplain found during the 2001 and 2002 monitoring programme (Wijnhoven et al. 2006). C. glareolus and S. araneus showed rather stable densities throughout the year. Other species showed density peaks during late summer and autumn. The highest densities of A. sylvaticus and M. minutus were found during late autumn and winter. The male-dominated trapping results for A. sylvaticus seem to be a common phenomenon (Randolph 1977). Kikkawa (1964) suggests that it results from a larger home range size and therefore greater trappability of male individuals.

With regard to the numbers of species trapped in flooded and nonflooded areas, the trapping results of A. sylvaticus and M. minutus in this study confirm the patterns found earlier by Wijnhoven et al. (2006), with more individuals of these species found in flooded than in nonflooded areas throughout the year. Wijnhoven et al. (2006) also found higher numbers of M. arvalis and C. russula in flooded areas. Only the fact that the largest proportion of C. glareolus was trapped in nonflooded areas is not in agreement with the 2001 and 2002 monitoring results, in which the reverse was found. We collected 56 specimens of C. glareolus for this study, which is much more than for most of the other species. This is not in line with the relatively low densities of this species found in 2001 and 2002.

Our measurements of metal concentrations in small-mammal species are representative of the values in the floodplain as a whole, and the previously mentioned differences between the data sets for the various species are generally in line with expectations for the species and the research area. Data for M. minutus should be interpreted with care because only four specimens were included in the analyses. C. glareolus might be exposed to different average metal concentrations in years of lower densities or when a larger proportion of the animals is present in the flooded areas. Variable factors—such as the frequency, timing, and duration of flooding and climatic variations and management measures, such as mowing and grazing in the area—can also influence the accumulation risks through effects on habitat suitability and connectivity and therefore on species abundance and distribution.

In accumulation and risk-assessment studies, most attention is given to shrews, in particular S. araneus, because these are the small mammals expected to show the highest accumulation of heavy metals in view of their insectivorous and carnivorous diet. However, with regard to the risks of accumulation and possible toxic effects to top predators, vole species are more relevant because they generally occur in much higher densities and are the dominant prey for several species (Jongbloed et al. 1996). In addition, the fact that a particular small-mammal species accumulates larger amounts of heavy metals does not necessarily mean that this is the species most at risk of toxic effects from pollutants. Some species could be more sensitive to heavy metals than others, and storage of heavy metals in organs, such as the liver, could be a good mechanism to cope with toxicants (Shore & Douben 1994). In our study, S. araneus did indeed show the highest heavy-metal concentrations of the seven small mammal species in the ADW floodplain area. For Cd, a maximum difference of a factor of 5.2 was found between S. araneus and the vole and mouse species. The differences are not as large as those found in inland areas or as those that have sometimes been suggested. A study by Ma et al. (1991), in polluted and relatively clean inland areas, found that estimated differences, by factors of 26 to 57, in Cd intake concentrations resulted in a difference of factors of 46 to 182 in kidney and 83 to 812 in liver concentrations. Hunter et al. (1989) recorded differences of factors of 2.1 to 12.1, 8.6 to 28.2, and 6.2 to 9.5 in kidney, liver, and muscle tissue, respectively, in control and polluted inland areas. The difference in bioaccumulation of Cd from soil to S. araneus and the voles (Microtidae) is a factor of 32 in the model described by Jongbloed et al. (1996). The estimated difference in exposure risk for shrews and voles, as calculated by Kooistra et al. (2001) for areas in the ADW floodplain, are factors of 2.1 to 9.8, which is more in line with our findings. For Pb, the difference between S. araneus and the vole and mouse species varied between 2.7 and 9.3 in our study. The difference of factors of 6.2 to 11.0 in Pb intake recorded by Ma et al. (1991), resulting in a difference of factors of 4.9 to 22.3 in kidney and 1.6 to 7.1 in liver concentrations is in line with our results. However, concentrations of Pb in C. russula and C. glareolus were not significantly different in our study. For Cu, the difference between S. araneus and the vole and mouse species varied by factors of 2.2 to 4.6 in our study. These values are similar or even slightly higher than the differences of 1.6 to 2.0, 1.6 to 4.1, and 1.2 to 1.7 in kidney, liver, and muscle tissue, respectively, as recorded by Hunter et al. (1989). The average Zn contents in S. araneus in our study did not differ by more than a factor of 2 from those in the Microtine rodent species. In our study, the metal concentrations in, for example, C. glareolus and M. agrestis, did differ less from those in S. araneus than recorded in the literature (for Cu and Cd) or were similar to them (for Zn and Pb). In our study, levels, except for Cd, in the shrew C. russula were similar to those in C. glareolus and also to those in several other vole and mouse species.

We expected that other than differences in metal concentrations between species of different trophic levels, the actual exposure concentrations would be also reflected in the metal concentrations in small-mammal species and individuals (inter and intraspecific variation). Because C. glareolus and M. agrestis were predominantly trapped in nonflooded areas, their exposure to total Zn concentrations in the soil was generally lower than that of the other species, whereas their exposure to CaCl2-extractable Zn concentrations was higher. This was, however, not reflected in the Zn concentrations in these two species. This suggests that the accumulation of Zn in these vole species is not determined by total or by CaCl2-extractable Zn fractions alone, indicating that Zn is regulated by either the food species (Heikens et al. 2001) or by the small mammals themselves (Mertens et al. 2001), which influences observed tissue values. With regard to exposure to Cu, Pb, and Cd, only exposure to total concentrations was different, i.e., lower for C. glareolus and M. agrestis, than for some of the other species depending on the investigated metal, which was not reflected in the metal concentrations in the species either.

Although exposure concentrations in the soil for C. russula were generally similar to those for other species, the metal concentrations show a trend toward concentrations being lower than in S. araneus and similar to those in voles. This could be the result of differences in feeding patterns between the two shrew species. Earthworms, which are known to be strong accumulators of heavy metals (Hobbelen et al. 2004; Van Vliet et al. 2005), are more important food items for S. araneus than for C. russula. It is also possible that macroinvertebrates are more often eaten by vole species, in particular, C. glareolus and M. agrestis, than was thought previously.

In addition to the influence of soil metal concentrations for accumulation of heavy metals, several interfering factors might be of importance. For several combinations of small mammal species and metals tested in multiple regressions, no regressions with measured parameters explaining a significant part of the variance (P < 0.05) in the observed metal concentrations in the animals were found. For M. agrestis, this can be ascribed to the small sample size and small variance, but also for S. araneus, with 45 specimens, no significant regressions were found. None of the parameters were substantially more often found to be related to the metal concentrations in the small mammals. This means that the total metal concentrations in the soil at the trapping location are not a good predictor of the metal concentrations in the small mammals found there. Using CaCl2-extractable concentrations in the soil hardly improved the relation with the accumulated metals in the small mammals. We conclude that the soil metal concentrations at the trapping locations do not necessarily reflect the exposure concentrations throughout the animals’ life history. This could be caused by (1) large variations in exposure time (or age of the animals); (2) heterogeneity of the soil concentrations, which means that exposure within an animal’s home range may not be similar to that at the trapping location; or (3) movements, dispersal, and/or shifts in feeding patterns, making such correlations irrelevant.

Life-stage and FW of the small mammals only occasionally showed significant regressions. This could reflect a poor relationship between FW and age because condition may play a role. S. araneus are known to show a decrease in body mass during winter (Ochocińska & Taylor 2003), which is probably also reflected as a larger variance in tissue-to-total body ratio (Table 2), but such a phenomenon is not assumed for the other species. Life-stage could be a poor predictor because small mammals remain juvenile for only a short period, causing the greatest variation in body metal concentrations within the group of adults. It seems most likely that the poor relations of soil metal concentrations with body metal concentrations result from the migration of species from one area to another or from shifts in diets. However, we found a relationship between the sex and body metal concentrations in three cases, with higher concentrations being found in female than in male individuals, suggesting that there are differences in feeding patterns or home-range structures between the sexes, or pregnancy effects on weight, for C. russula and M. arvalis.

It is difficult to assess the exposure risk and contaminant levels of individual small mammals in the highly dynamic floodplain environment. However, based on our results, it is nevertheless possible to calculate average risks for species within ecologic units or floodplains as a whole. It is important that the animals collected reflect the species distribution and densities throughout the year, making it crucial to combine measurements of metal accumulation with population monitoring.

Metal concentrations in wild small mammals are often measured in the target organs (e.g., liver and kidneys) (Hunter & Johnson 1982; Hunter et al. 1984; Ma et al. 1991; Damek-Poprawa & Sawicka-Kapustra 2003, 2004). Although this is relevant from a toxicologic point of view, total body concentrations may be more relevant when considering small mammals as prey animals. The observed differences in metal concentrations between small-mammal species have consequences for the risks to predators in floodplains. Species preying on S. araneus, for instance, run greater risks of heavy-metal accumulation than species preying on voles and mice. However, S. araneus is generally not the most important prey species for predators (Jongbloed et al. 1996) because the densities of other species are often higher (Wijnhoven et al. 2006). Predators will forage where their preferred prey is available in large numbers, which means that in the ADW floodplain, exposure of predators specialising in bank voles mainly occurs in nonflooded areas, whereas exposure of predators of M. arvalis occurs especially in flooded areas. This leads to differences in exposure between predators, such as the weasel (Mustela nivalis) and the Eurasian kestrel (Falco tinnunculus), as a result of their different foraging behaviour (Erlinge et al. 1983), distribution of their prey, and contaminant levels within the prey. During floods, as well as for quite some time after an area has been flooded, exposure of predators of small mammals to heavy metals will generally occur in nonflooded areas. Predators are probably sparse in the area during this time because they are forced to move away because of low prey availability, or, in the case of generalists, forced to prey on other food sources (Van den Brink et al. 2003).

Differences in metal concentration ranges between the small-mammal species seem to be small for Zn, which is probably to a certain extent regulated. Extremely high levels of metals were found in a few individual voles: Zn and Pb in C. glareolus, Cd in M. agrestis, and Cu in M. arvalis. Although these extreme values could not be directly related to soil concentrations, they can possibly be explained by point sources of pollution and by foraging in other areas than the sample sites alone. All individuals with extreme body concentrations were trapped at or near nonflooded areas. Only in a few cases were the metal concentrations in the mammals related to soil metal concentrations, either total or extractable. We know that this floodplain is diffusely polluted, which implies heterogeneous contamination patterns. However, large differences in contaminant levels over short distances are not to be expected, unless there are certain point sources in the nonflooded areas where industrial activities have taken place in the past’ or near the borders between the flooded and nonflooded parts. It is more likely that in this dynamic environment, exposure of small mammals at or near the trapping location only occurs during a part of their life history. Movements and dispersal probably interfere with the relations between soil metal concentrations and the concentrations in the bodies of the mammals. In addition, the metal concentrations in the small mammals could also be affected by shifts in feeding patterns, seasonal and flood-related aspects of food availability, habitat suitability and connectivity, and life-stage–related food preference, combined with the same variations in the metal contents in the food items themselves. Finally, exposure time, and therefore age of the animals, might be an explanatory factor, as was indicated in a few cases by the relation between metal concentrations in small mammals and life stage, FW, or season.

Conclusion

Heavy-metal concentrations in seven small-mammal species from a diffusely polluted floodplain along the river Rhine differed between species. These differences can be partly explained by the trophic level of the species because metal concentrations were highest in the carnivorous and insectivorous shrew, S. araneus. However, average metal concentrations in S. araneus differed from those in the vole and mouse species at the maximum by factors of 2.0, 4.6, 9.3, and 5.2 for Zn, Cu, Pb, and Cd, respectively, which is less than has been reported by several studies of inland areas. No significant differences were found in average Zn, Pb, and Cd concentrations between S. araneus, M. agrestis, and C. russula. The vole C. glareolus and the shrew C. russula had similar concentrations of Cu, Pb, and Zn. Although possible differences in accumulation between life-stages and sexes or size- and season-related differences were corrected for, relationships between total or CaCl2-extractable soil metal concentrations at the trapping locations and the metal concentrations in the mammals were poor or absent. We suspect that exposure time, dispersal, and changes in foraging behaviour might be important factors influencing the exposure of small mammals in highly dynamic environments, such as floodplains. The observed differences in metal concentrations between small-mammal species will have consequences in terms of risks to predators in floodplains and should therefore be considered in risk assessments.