# The Solid Phase Distribution and Bioaccessibility of Arsenic, Chromium, and Nickel in Natural Ironstone Soils in the UK

Authors: Joanna Wragg; Mark Cave; Sean Gregory
Journal: Applied and Environmental Soil Science
 (2014)

---

## Abstract

Thirty soil samples (12 residential gardens and 18 allotments) were collected from the Cherwell District of north Oxfordshire in south-central England. The underlying parent geology of the area is dominated by Jurassic ironstone. The samples were analysed for their total contents of As, Cr, and Ni by X-ray fluorescence spectroscopy and for the bioaccessible fractions of these elements using a physiologically based extraction test. Four soils (two residential soils and two allotment soils) were chosen for further determination of their element solid phase distribution. The study showed that whilst total concentrations of As, Cr, and Ni are elevated due to the soil parent material, the bioaccessibility test showed that only a small proportion of the total concentration is available for absorption into the human body (<15%). The sequential extraction test showed that the nonmobile forms of the elements are strongly sorbed on to iron oxides. Parent material geology has a significant effect on the total element concentrations and the bioaccessibility of potentially harmful element (PHE). Land use does not show such a large effect but the allotment bioaccessibility data show a bigger spread and possibly higher values for As and Cr which may be due to agronomic (cultivation) practices such as addition of fertilisers and organic matter.


---

## Body

## 1. Introduction

The mobility of potentially harmful elements (PHEs) present in soils is influenced by chemical reactions and both physical and biological processes occurring within soils. Commonly encountered soil PHEs include arsenic (As), chromium (Cr), and nickel (Ni), all of which can be harmful to human health via oral ingestion, inhalation, or dermal uptake [1, 2]. Naturally occurring soil PHE in the environment results from weathering of parent rock material, with As being associated with igneous rocks, but more commonly the highest concentrations tend to be observed in sedimentary rocks (e.g., sandstones, shales, and mudstones) and in heavily sulphidic mineralised areas [3]. Natural As enrichment of a parent material, such as a sedimentary ironstone, is thought to be due to the abundance of iron-bearing minerals, especially iron oxyhydroxides, as these have a high affinity for As [4]. Soil Cr is abundant in basic, ultrabasic, and igneous rocks and shales and in topsoils mirrors the clay content [5]. Weathering of basic igneous rocks results in a large proportion of soil Ni, although lower concentrations are found in association with sedimentary rocks including shales, clays, limestones, and sandstones [3, 6]. All three elements of interest in this study are found, often colocated, in a broad band of clay that includes the Jurassic ironstones and stretches from Bristol to the east of Nottingham and up to Scunthorpe [5]. Industrial activity in the UK, often related to the mining and smelting of mineralised PHE, has resulted in elevated soil concentrations of these elements in urban areas ranging from Swansea in south Wales to Derby and Nottingham in the Midlands and further the industrial belt of Yorkshire and Lancashire [7].Mean soil element concentrations for England and Wales have recently been published in the “Advanced soil geochemical atlas of England and Wales” [5]. This document presents distribution maps for a total of 53 elements, which include 17 from the original study by the National Soil Inventory (NSI) and the remainder were generated by reanalysis of the original collected soil samples. The ranges of soil concentrations reported were 0–820 mg kg−1, 5.1–1141 mg kg−1, and 0.26–459 mg kg−1 for As, Cr, and Ni, respectively. The mean concentrations of As, Cr, and Ni were 20 mg kg−1, 68 mg kg−1, and 23 mg kg−1.A number of researchers have shown the association of PHEs such as As, Cr, and Ni with Jurassic ironstone derived soils, in particular the Frodingham and Claxby ironstones [8–11]. Soil contaminant mobility research [8–11] has shown that the presence of PHEs in soils may not be mobile and hence does not present a risk to human health. Such studies have linked contaminant mobility, via the measurement of in vitro bioaccessibility, to the source apportionment of soil contaminants and/or properties of soils [12], where bioaccessibility is defined as the fraction of a contaminant that is soluble in the gastrointestinal environment and is available for absorption [13]. To date, a large proportion of this research has focussed on As, as this PHE is found to be a widespread naturally occurring element in UK soils and is associated with the Jurassic ironstone geology from areas such as Lincolnshire and Northamptonshire [10, 11, 14]. But there is little work available detailing studies on other PHE such as Cr and Ni from ironstones from other areas of the UK. The area around the town of Banbury in the UK shows particularly elevated concentrations of As, Cr, and Ni compared to the surrounding areas [5]. The purpose of this work was to be the first study to assess the potential availability of As, Cr, and Ni present in ironstone derived soil samples from the Cherwell District, (Oxfordshire, UK) via the oral ingestion exposure route using a combination of geochemical tests: sequential extraction to identify the fractionation of the PHE within the geochemical components that make up the soil and in vitro bioaccessibility measurements to identify which components are the source(s) of the bioaccessible fraction.


## 2. Materials and Methods

All reagents were of analytical grade or better. Reagents for the analysis of extracts and digests were supplied by Merck and all reagents used for the PBET extraction were supplied by Sigma-Aldrich Ltd.


### 2.1. Sample Locations, Collection, and Preparation

A total of 30 soil samples were collected from seven locations across the Cherwell District of north Oxfordshire in south central England in July 2008 (Table1). The underlying parent geology of the area is dominated by the Northampton Sand Formation [15]. The samples were stratified by underlying geology (Lias n=18 or Inferior Oolite Group n=12), by location, and by the type of site from which they were collected (residential garden (n=12) or allotment (n=18)), shown in Figure 1. Samples were collected for the assessment of the potential human exposure to soil contamination via the ingestion pathway according to the standard British Geological Survey (BGS) approach for geochemical surveying work, collecting composite of material from hand auger flights taken from five holes of 0–15 cm depth [16]. For this work the support area was one square meter. The composite soil materials were stored in Kraft paper bags and uniquely labelled with the date and the GPS coordinates of the sample location. After collection, the samples were dried at 35 ± 2°C on metal trays for at least 12 hours until they were visually dry, at which time each sample was gently disaggregated, to ensure the breakage of aggregates but retention of clasts, and sieved to extract the <250 μm fraction, using a nylon sieve mesh. This fraction is the upper bound of size fraction which is considered to adhere to the hands of the “at risk” receptor (children) considered in this study [17].Table 1
Summary of geological and sampling information and chemical measurement data.
Sample
Geology
Sampling location
Type
Organic matter (%)
Soil pH
Total As (mg kg−1)
Bioaccessible As (mg kg−1)
As BAF (%)
Total Cr (mg kg−1)
Bioaccessible Cr (mg kg−1)
Cr BAF (%)
Total Ni (mg kg−1)
Bioaccessible Ni (mg kg−1)
Ni BAF (%)
GBA  number 1
L
GB
A
12.7
6.86
153
12.8
8.3
275
2.0
0.74
108
3.3
3.1
GBA  number 2
L
GB
A
14.6
7.18
112
5.4
4.8
302
1.1
0.4
114
2.6
2.3
GBA  number 3
L
GB
A
13.4
6.70
82
4.9
6.0
203
1.4
0.7
88
2.9
3.3
GBR  number 1
L
GB
R
10.4
6.29
84
4.2
5.1
264
1.4
0.5
97
4.0
4.1
GBR  number 2
L
GB
R
10.8
6.41
92
3.8
4.1
274
1.4
0.5
98
4.4
4.5
SA  number 11
L
S
R
13.0
6.85
337
8.1
2.4
379
1.1
0.3
185
3.9
2.1
SA  number 2
L
S
A
12.4
6.25
348
7.7
2.2
402
0.9
0.2
190
4.4
2.3
SA  number 3
L
S
A
11.0
7.01
305
5.6
1.8
382
0.9
0.2
176
3.2
1.8
SR  number 1
L
S
R
13.8
5.83
392
6.0
1.5
406
0.9
0.2
207
4.4
2.1
SR  number 21
L
W
R
12.9
5.97
417
5.0
1.2
447
0.8
0.2
218
4.2
1.9
WA  number 1
L
W
A
12.4
7.12
289
7.4
2.6
321
1.1
0.3
159
3.6
2.3
WA  number 2
L
W
A
14.3
7.11
230
7.7
3.3
273
1.1
0.4
136
2.9
2.1
WA  number 3
L
W
A
12.0
7.09
249
6.2
2.5
293
1.1
0.4
147
3.0
2.1
WR  number 1
L
W
R
15.8
6.98
250
5.9
2.4
342
1.0
0.3
153
3.6
2.4
WR  number 2
L
W
R
13.9
7.11
280
6.3
2.3
359
1.1
0.3
158
3.6
2.3
LHA  number 1
IO
LH
A
7.71
7.25
17
1.9
11.0.
65
0.5
0.8
25
1.4
5.7
LHA  number 2
IO
LH
A
8.52
7.06
14
1.6
11.4
69
0.6
0.9
28
1.6
5.8
LHA  number 3
IO
LH
A
9.21
7.27
18
2.1
11.6
68
0.6
0.9
31
1.6
5.3
LHR  number 1
IO
LH
R
11.3
7.43
38
4.0
10.5
96
0.7
0.7
60
1.3
2.1
LHR  number 2
IO
LH
R
9.05
7.54
29
2.9
10.1
83
0.6
0.7
46
1.2
2.6
SAA  number 1
IO
SA
R
5.15
7.11
31
3.6
11.7
51
0.7
1.3
17
1.2
6.9
SAA  number 2
IO
SA
A
6.29
6.72
50
3.6
7.1
53
1.1
2.1
21
1.5
7.0
SAA  number 3
IO
SA
A
8.23
6.79
23
2.9
12.4
57
1.3
2.2
21
1.4
6.7
SAR  number 1
IO
SA
R
6.63
6.92
15
2.2
14.7
52
0.8
1.6
18
1.6
8.7
SAR  number 2
IO
SA
R
4.90
7.16
18
2.5
14.1
54
0.7
1.3
19
1.6
8.2
DA  number 1
L
D
A
16.1
7.24
217
7.4
3.4
340
1.4
0.4
154
4.3
2.8
DA  number 21
L
D
A
13.9
7.28
201
6.4
3.2
364
1.4
0.4
148
4.4
3.0
DA  number 3
L
D
A
16.0
7.18
178
7.0
3.9
321
1.5
0.5
142
4.5
3.2
UHR  number 11
IO
UH
R
12.0
6.75
20
3.1
15.3
98
1.1
1.1
60
3.2
5.3
UHR  number 2
IO
UH
R
9.42
6.88
16
2.4
14.8
80
0.9
1.1
44
2.4
5.4
L: Lias; IO: Inferior Oolite Groups; GB: Great Bourton; S: Shenington; W: Wroxton; LH: Lower Heyford; SA: Steeple Aston; D: Drayton; UH: Upper Heyford; A: allotment; R: residential sampling site.  Total and Bioaccessible PHE concentrations as mg kg−1; BAF denotes bioaccessible fraction (%). pH uncertainty is ±0.05 units and organic matter is ±1%. 1indicates samples used for sequential extraction.Figure 1
Local geology for the study site.


### 2.2. Total and Bioaccessible As, Cr, and Ni Extractions

Total As, Cr, and Ni were determined on a pressed powder pellet by Phillips wavelength dispersive WD-XRF, fitted with automatic sample changer, 60 kV generator, and controlled by dedicated PC running PANalytical Super XRF software. The pressed powder pellets were prepared by grinding 12.0 g of the <250μm soil material and 3.00 g of binder (EMU120FD, a styrene copolymer, and Ceridust 3620, a micronised polyethylene wax in a ratio of 9 : 1) in an agate planetary ball mill for 30 minutes followed by pressing into 40 mm diameter pellets at a 25 ton load. The instrument was calibrated using 1000 mg kg−1 single element standards in a silica matrix; background factors, where applicable, were calculated by angular difference or from regression values from “high purity” oxide blanks and line overlap factors were calculated from high concentration, multielement standards of the interfering analyte [18]. Uncertainties on the XRF data are ±10%.A physiologically based extraction test (PBET), a modified version of that first described by Ruby et al. [19], was applied to the test soils to determine the proportion of As, Cr, and Ni that may be available for uptake via accidental ingestion. The PBET utilised was the BGS PBET, which has been previously fully described by Cave et al. [9], Palumbo-Roe et al. [10], and Wragg et al. [11], but which, in essence, is a two stage sequential extraction that simulates the physicochemical gastrointestinal human environment. Each soil sample (1 g) was incubated at body temperature (37°C) in an end-over-end rotator at 30 rpm with 100 mL of simulated gastric solution (1.25 g pepsin, 0.50 g sodium malate, 0.50 g sodium citrate, 420 μL lactic acid, and 500 μL acetic acid made up to 1 L with freshly prepared deionised water, adjusted to pH 2.5 with concentrated HCl) for one hour to simulate stomach conditions in the fasted state. After a 1 hour incubation period, a 5 mL sample was collected, via disposable syringe, and the sample filtered through a 0.45 μm cellulose acetate syringe filter. To ensure that the solid : solution ratio was maintained, the filter was reversed and backflushed with a fresh 5 mL aliquot of gastric solution, thereby returning both the original volume of fluid and mass of solid. To simulate the small intestine, the conditions in each extraction vessel were then altered by titration against saturated NaHCO3 until a pH of 7.0 was obtained and the addition of 175 mg bile salts and 50 mg pancreatin. These conditions were maintained at 37°C for 4 hours. After an initial 2-hour period under small intestinal conditions, a 5 mL sample was collected via disposable syringe and cellulose acetate syringe filter and then again after 2 hours (4 hours in total) to complete the extraction. All samples were stored for analysis by Inductively Coupled Plasma Atomic Emission Spectrometry (ICP-AES) at 0–8°C. Uncertainties of the PBET method are ±15%.


### 2.3. Solid Phase Distribution of As, Cr, and Ni

Four soils (two residential soils and two allotment soils), with total As, Cr, and Ni concentrations covering the range found in in all soil samples, were chosen for further determination of their PHE solid phase distribution. The solid phase distributions of these elements were determined by applying a nonspecific sequential extraction to the soils, known as the Chemometric Identification of Substrates and Element Distributions (CISED). The methodology, developed by Cave et al. [20], applied increasing concentrations of aqua regia to the solid material under investigation, which was supported on a 0.45 μm polytetrafluoroethylene (PTFE) membrane filter tube insert housed inside a “Whatman Vectaspin 20” centrifuge tube. 10 mL of the nonspecific reagent was added, via pipette, to a 2 g aliquot of the test material and the extraction vessel was centrifuged at 1034 G for 10 minutes, and the extract was removed and stored in a clean, dry, uniquely labelled 16.5 mL screw top polypropylene conical tube and stored at 0–8°C prior to analysis by ICP-AES. The process was repeated sequentially until 14 extracts for each solid material were produced. Each reagent was applied to the soil in duplicate. The extractants applied were deionised water followed by 0.01, 0.05, 0.10, 0.50, 1.00, and 5.00 M aqua regia.


### 2.4. Analysis of Bioaccessible and CISED

All extracted elements of interest, As, Cr, and Ni for the PBET solutions and the major and trace cations (Al, As, B, Ba, Ca, Cd, Co, Cr, Cu, Fe, K, Li, Mg, Mn, Na, Ni, P, Pb, S, Si, Sr, V, and Zn) of interest in the CISED extracts were determined directly by a Varian Vista AX CCD simultaneous instrument with dedicated SPS-5 autosampler and PC running ICP Expert Version 3.1b as supplied by the manufacturer. All instrument operating conditions have been previously reported [21, 22]. Instrument calibration was carried out using a series of a minimum of four mixed element matrix matched standards, covering a linear range up to 10 mg L−1 for Al, As, B, Ba, Cd, Co, Cr, Cu, Li, Ni, P, Pb, Si, Sr, V, and Zn; 50 mg L−1 for K; 100 mg L−1 for Fe, Mn, and S; 500 mg L−1 for Mg and Ca; and 0–1000 mg L−1 for Na. For each sample type, a solution of 1% CsCl2 was introduced into the sample introduction stream to increase the linear range of the calibration, therefore reducing the number of analytical dilutions required. All reported measurements, in mg L−1, are based on the average of three ten-second replicate analyses. All element concentrations quoted have been converted into mg kg−1 extracted from the solid matrix and the PBET value used for further interpretation is the highest value from either the stomach or one of the two intestinal compartments.


### 2.5. Statistical Analysis and Modelling of the CISED Extraction Data

Data processing of the ICP-AES analysis of the CISED extraction data was carried out using a modified version of the original self-modelling mixture resolution (SMMR) method described by Cave et al. [20], using the Matlab© programming language. Full details of the modifications and their validation have been described by Cave [23]. Data processing was carried out to determine the number of soil components or substrates in each test material, the composition of the identified components, and the distribution of PHE between of the identified components. Previous applications of the methodology to soils with elevated PHE concentrations derived from Jurassic ironstone formations and mineralised geologies have been reported by Cave et al. [9], Palumbo-Roe et al. [10], Wragg [22], and Wragg et al. [11]. Four samples selected on the basis of being representative of the range of PHE found in the main set of test soils were subjected to the CISED extraction. In order to categorise the soil host substrates and their PHE distributions across the test materials in a consistent manner, a clustering methodology and associated data visualisation technique as previously used by Wragg [22] and Cox et al. [24] were applied to the SMMR data. Briefly, this involves setting up a data matrix consisting of the percentage element composition (Al, As, Ca, Cr, Fe, K, Mg, Mn, Na, Ni, P, S, Si) combined with the amount of material extracted for each substrate (mg kg−1) in each soil for all soils under investigation. The matrix was subjected to hierarchical clustering, where the data were mean centred and scaled with Euclidean distance and linkage using Ward’s method. The resulting colourmap was then used to identify chemical composition of the host substrates and the locations at which they were associated.The R programming language [25] was used to carry out the t-tests used to assess whether the total element and the bioaccessible concentrations of As, Cr, and Ni in the soils were significantly different when partitioned by geology and by land use.


### 2.6. Quality Control

The quality of the XRF and ICP-AES analysis and bioaccessibility and CISED extractions was checked according to standard operating procedures [18]. This included the regular analysis of quality control samples, reference materials, blanks, and extraction and analysis of sample duplicates. For the bioaccessibility testing, the materials that were used as quality control soils are one National Institute of Standards and Technology (NIST) Standard Reference Material (SRM), 2711, and one bioaccessibility guidance soil, BGS 102. The mean bioaccessibility repeatability for all duplicate extractions was <10% for both the stomach and the intestinal phases of the PBET. The bioaccessible values obtained (the highest value obtained for the two phases extracted) for As and Cr in the BGS 102 reference material [26] (n=4) were compared against current consensus values for this material generated using the Unified BARGE method (UBM) [26–28], a similar but more complex in vitro bioaccessibility method, developed since the undertaking of this study. Both As and Cr are in agreement, within error limits, with reported values of 3.77 ± 0.49 mg kg−1 and 9.15 ± 0.40 mg kg−1 compared to consensus values of 5.38 ± 2.39 mg kg−1 and 9.60 ± 2.00 mg kg−1 for As and Cr, respectively. However, there are currently no bioaccessibility values available for Ni. All blank extractions (n=4) returned As, Cr, and Ni values below the limit of quantification for this extraction type.For the CISED extraction, one sample was extracted in duplicate to assess the extraction repeatability and one blank extraction was carried out to assess any contribution made from the reagent blank. The mean repeatability across the 23 elements analysed was 4.65%. The blank extraction data indicates that the reagent blank did not contribute any elements of interest to the dataset as all results were below the limit of quantification for this extraction type.


## 2.1. Sample Locations, Collection, and Preparation

A total of 30 soil samples were collected from seven locations across the Cherwell District of north Oxfordshire in south central England in July 2008 (Table1). The underlying parent geology of the area is dominated by the Northampton Sand Formation [15]. The samples were stratified by underlying geology (Lias n=18 or Inferior Oolite Group n=12), by location, and by the type of site from which they were collected (residential garden (n=12) or allotment (n=18)), shown in Figure 1. Samples were collected for the assessment of the potential human exposure to soil contamination via the ingestion pathway according to the standard British Geological Survey (BGS) approach for geochemical surveying work, collecting composite of material from hand auger flights taken from five holes of 0–15 cm depth [16]. For this work the support area was one square meter. The composite soil materials were stored in Kraft paper bags and uniquely labelled with the date and the GPS coordinates of the sample location. After collection, the samples were dried at 35 ± 2°C on metal trays for at least 12 hours until they were visually dry, at which time each sample was gently disaggregated, to ensure the breakage of aggregates but retention of clasts, and sieved to extract the <250 μm fraction, using a nylon sieve mesh. This fraction is the upper bound of size fraction which is considered to adhere to the hands of the “at risk” receptor (children) considered in this study [17].Table 1
Summary of geological and sampling information and chemical measurement data.
Sample
Geology
Sampling location
Type
Organic matter (%)
Soil pH
Total As (mg kg−1)
Bioaccessible As (mg kg−1)
As BAF (%)
Total Cr (mg kg−1)
Bioaccessible Cr (mg kg−1)
Cr BAF (%)
Total Ni (mg kg−1)
Bioaccessible Ni (mg kg−1)
Ni BAF (%)
GBA  number 1
L
GB
A
12.7
6.86
153
12.8
8.3
275
2.0
0.74
108
3.3
3.1
GBA  number 2
L
GB
A
14.6
7.18
112
5.4
4.8
302
1.1
0.4
114
2.6
2.3
GBA  number 3
L
GB
A
13.4
6.70
82
4.9
6.0
203
1.4
0.7
88
2.9
3.3
GBR  number 1
L
GB
R
10.4
6.29
84
4.2
5.1
264
1.4
0.5
97
4.0
4.1
GBR  number 2
L
GB
R
10.8
6.41
92
3.8
4.1
274
1.4
0.5
98
4.4
4.5
SA  number 11
L
S
R
13.0
6.85
337
8.1
2.4
379
1.1
0.3
185
3.9
2.1
SA  number 2
L
S
A
12.4
6.25
348
7.7
2.2
402
0.9
0.2
190
4.4
2.3
SA  number 3
L
S
A
11.0
7.01
305
5.6
1.8
382
0.9
0.2
176
3.2
1.8
SR  number 1
L
S
R
13.8
5.83
392
6.0
1.5
406
0.9
0.2
207
4.4
2.1
SR  number 21
L
W
R
12.9
5.97
417
5.0
1.2
447
0.8
0.2
218
4.2
1.9
WA  number 1
L
W
A
12.4
7.12
289
7.4
2.6
321
1.1
0.3
159
3.6
2.3
WA  number 2
L
W
A
14.3
7.11
230
7.7
3.3
273
1.1
0.4
136
2.9
2.1
WA  number 3
L
W
A
12.0
7.09
249
6.2
2.5
293
1.1
0.4
147
3.0
2.1
WR  number 1
L
W
R
15.8
6.98
250
5.9
2.4
342
1.0
0.3
153
3.6
2.4
WR  number 2
L
W
R
13.9
7.11
280
6.3
2.3
359
1.1
0.3
158
3.6
2.3
LHA  number 1
IO
LH
A
7.71
7.25
17
1.9
11.0.
65
0.5
0.8
25
1.4
5.7
LHA  number 2
IO
LH
A
8.52
7.06
14
1.6
11.4
69
0.6
0.9
28
1.6
5.8
LHA  number 3
IO
LH
A
9.21
7.27
18
2.1
11.6
68
0.6
0.9
31
1.6
5.3
LHR  number 1
IO
LH
R
11.3
7.43
38
4.0
10.5
96
0.7
0.7
60
1.3
2.1
LHR  number 2
IO
LH
R
9.05
7.54
29
2.9
10.1
83
0.6
0.7
46
1.2
2.6
SAA  number 1
IO
SA
R
5.15
7.11
31
3.6
11.7
51
0.7
1.3
17
1.2
6.9
SAA  number 2
IO
SA
A
6.29
6.72
50
3.6
7.1
53
1.1
2.1
21
1.5
7.0
SAA  number 3
IO
SA
A
8.23
6.79
23
2.9
12.4
57
1.3
2.2
21
1.4
6.7
SAR  number 1
IO
SA
R
6.63
6.92
15
2.2
14.7
52
0.8
1.6
18
1.6
8.7
SAR  number 2
IO
SA
R
4.90
7.16
18
2.5
14.1
54
0.7
1.3
19
1.6
8.2
DA  number 1
L
D
A
16.1
7.24
217
7.4
3.4
340
1.4
0.4
154
4.3
2.8
DA  number 21
L
D
A
13.9
7.28
201
6.4
3.2
364
1.4
0.4
148
4.4
3.0
DA  number 3
L
D
A
16.0
7.18
178
7.0
3.9
321
1.5
0.5
142
4.5
3.2
UHR  number 11
IO
UH
R
12.0
6.75
20
3.1
15.3
98
1.1
1.1
60
3.2
5.3
UHR  number 2
IO
UH
R
9.42
6.88
16
2.4
14.8
80
0.9
1.1
44
2.4
5.4
L: Lias; IO: Inferior Oolite Groups; GB: Great Bourton; S: Shenington; W: Wroxton; LH: Lower Heyford; SA: Steeple Aston; D: Drayton; UH: Upper Heyford; A: allotment; R: residential sampling site.  Total and Bioaccessible PHE concentrations as mg kg−1; BAF denotes bioaccessible fraction (%). pH uncertainty is ±0.05 units and organic matter is ±1%. 1indicates samples used for sequential extraction.Figure 1
Local geology for the study site.


## 2.2. Total and Bioaccessible As, Cr, and Ni Extractions

Total As, Cr, and Ni were determined on a pressed powder pellet by Phillips wavelength dispersive WD-XRF, fitted with automatic sample changer, 60 kV generator, and controlled by dedicated PC running PANalytical Super XRF software. The pressed powder pellets were prepared by grinding 12.0 g of the <250μm soil material and 3.00 g of binder (EMU120FD, a styrene copolymer, and Ceridust 3620, a micronised polyethylene wax in a ratio of 9 : 1) in an agate planetary ball mill for 30 minutes followed by pressing into 40 mm diameter pellets at a 25 ton load. The instrument was calibrated using 1000 mg kg−1 single element standards in a silica matrix; background factors, where applicable, were calculated by angular difference or from regression values from “high purity” oxide blanks and line overlap factors were calculated from high concentration, multielement standards of the interfering analyte [18]. Uncertainties on the XRF data are ±10%.A physiologically based extraction test (PBET), a modified version of that first described by Ruby et al. [19], was applied to the test soils to determine the proportion of As, Cr, and Ni that may be available for uptake via accidental ingestion. The PBET utilised was the BGS PBET, which has been previously fully described by Cave et al. [9], Palumbo-Roe et al. [10], and Wragg et al. [11], but which, in essence, is a two stage sequential extraction that simulates the physicochemical gastrointestinal human environment. Each soil sample (1 g) was incubated at body temperature (37°C) in an end-over-end rotator at 30 rpm with 100 mL of simulated gastric solution (1.25 g pepsin, 0.50 g sodium malate, 0.50 g sodium citrate, 420 μL lactic acid, and 500 μL acetic acid made up to 1 L with freshly prepared deionised water, adjusted to pH 2.5 with concentrated HCl) for one hour to simulate stomach conditions in the fasted state. After a 1 hour incubation period, a 5 mL sample was collected, via disposable syringe, and the sample filtered through a 0.45 μm cellulose acetate syringe filter. To ensure that the solid : solution ratio was maintained, the filter was reversed and backflushed with a fresh 5 mL aliquot of gastric solution, thereby returning both the original volume of fluid and mass of solid. To simulate the small intestine, the conditions in each extraction vessel were then altered by titration against saturated NaHCO3 until a pH of 7.0 was obtained and the addition of 175 mg bile salts and 50 mg pancreatin. These conditions were maintained at 37°C for 4 hours. After an initial 2-hour period under small intestinal conditions, a 5 mL sample was collected via disposable syringe and cellulose acetate syringe filter and then again after 2 hours (4 hours in total) to complete the extraction. All samples were stored for analysis by Inductively Coupled Plasma Atomic Emission Spectrometry (ICP-AES) at 0–8°C. Uncertainties of the PBET method are ±15%.


## 2.3. Solid Phase Distribution of As, Cr, and Ni

Four soils (two residential soils and two allotment soils), with total As, Cr, and Ni concentrations covering the range found in in all soil samples, were chosen for further determination of their PHE solid phase distribution. The solid phase distributions of these elements were determined by applying a nonspecific sequential extraction to the soils, known as the Chemometric Identification of Substrates and Element Distributions (CISED). The methodology, developed by Cave et al. [20], applied increasing concentrations of aqua regia to the solid material under investigation, which was supported on a 0.45 μm polytetrafluoroethylene (PTFE) membrane filter tube insert housed inside a “Whatman Vectaspin 20” centrifuge tube. 10 mL of the nonspecific reagent was added, via pipette, to a 2 g aliquot of the test material and the extraction vessel was centrifuged at 1034 G for 10 minutes, and the extract was removed and stored in a clean, dry, uniquely labelled 16.5 mL screw top polypropylene conical tube and stored at 0–8°C prior to analysis by ICP-AES. The process was repeated sequentially until 14 extracts for each solid material were produced. Each reagent was applied to the soil in duplicate. The extractants applied were deionised water followed by 0.01, 0.05, 0.10, 0.50, 1.00, and 5.00 M aqua regia.


## 2.4. Analysis of Bioaccessible and CISED

All extracted elements of interest, As, Cr, and Ni for the PBET solutions and the major and trace cations (Al, As, B, Ba, Ca, Cd, Co, Cr, Cu, Fe, K, Li, Mg, Mn, Na, Ni, P, Pb, S, Si, Sr, V, and Zn) of interest in the CISED extracts were determined directly by a Varian Vista AX CCD simultaneous instrument with dedicated SPS-5 autosampler and PC running ICP Expert Version 3.1b as supplied by the manufacturer. All instrument operating conditions have been previously reported [21, 22]. Instrument calibration was carried out using a series of a minimum of four mixed element matrix matched standards, covering a linear range up to 10 mg L−1 for Al, As, B, Ba, Cd, Co, Cr, Cu, Li, Ni, P, Pb, Si, Sr, V, and Zn; 50 mg L−1 for K; 100 mg L−1 for Fe, Mn, and S; 500 mg L−1 for Mg and Ca; and 0–1000 mg L−1 for Na. For each sample type, a solution of 1% CsCl2 was introduced into the sample introduction stream to increase the linear range of the calibration, therefore reducing the number of analytical dilutions required. All reported measurements, in mg L−1, are based on the average of three ten-second replicate analyses. All element concentrations quoted have been converted into mg kg−1 extracted from the solid matrix and the PBET value used for further interpretation is the highest value from either the stomach or one of the two intestinal compartments.


## 2.5. Statistical Analysis and Modelling of the CISED Extraction Data

Data processing of the ICP-AES analysis of the CISED extraction data was carried out using a modified version of the original self-modelling mixture resolution (SMMR) method described by Cave et al. [20], using the Matlab© programming language. Full details of the modifications and their validation have been described by Cave [23]. Data processing was carried out to determine the number of soil components or substrates in each test material, the composition of the identified components, and the distribution of PHE between of the identified components. Previous applications of the methodology to soils with elevated PHE concentrations derived from Jurassic ironstone formations and mineralised geologies have been reported by Cave et al. [9], Palumbo-Roe et al. [10], Wragg [22], and Wragg et al. [11]. Four samples selected on the basis of being representative of the range of PHE found in the main set of test soils were subjected to the CISED extraction. In order to categorise the soil host substrates and their PHE distributions across the test materials in a consistent manner, a clustering methodology and associated data visualisation technique as previously used by Wragg [22] and Cox et al. [24] were applied to the SMMR data. Briefly, this involves setting up a data matrix consisting of the percentage element composition (Al, As, Ca, Cr, Fe, K, Mg, Mn, Na, Ni, P, S, Si) combined with the amount of material extracted for each substrate (mg kg−1) in each soil for all soils under investigation. The matrix was subjected to hierarchical clustering, where the data were mean centred and scaled with Euclidean distance and linkage using Ward’s method. The resulting colourmap was then used to identify chemical composition of the host substrates and the locations at which they were associated.The R programming language [25] was used to carry out the t-tests used to assess whether the total element and the bioaccessible concentrations of As, Cr, and Ni in the soils were significantly different when partitioned by geology and by land use.


## 2.6. Quality Control

The quality of the XRF and ICP-AES analysis and bioaccessibility and CISED extractions was checked according to standard operating procedures [18]. This included the regular analysis of quality control samples, reference materials, blanks, and extraction and analysis of sample duplicates. For the bioaccessibility testing, the materials that were used as quality control soils are one National Institute of Standards and Technology (NIST) Standard Reference Material (SRM), 2711, and one bioaccessibility guidance soil, BGS 102. The mean bioaccessibility repeatability for all duplicate extractions was <10% for both the stomach and the intestinal phases of the PBET. The bioaccessible values obtained (the highest value obtained for the two phases extracted) for As and Cr in the BGS 102 reference material [26] (n=4) were compared against current consensus values for this material generated using the Unified BARGE method (UBM) [26–28], a similar but more complex in vitro bioaccessibility method, developed since the undertaking of this study. Both As and Cr are in agreement, within error limits, with reported values of 3.77 ± 0.49 mg kg−1 and 9.15 ± 0.40 mg kg−1 compared to consensus values of 5.38 ± 2.39 mg kg−1 and 9.60 ± 2.00 mg kg−1 for As and Cr, respectively. However, there are currently no bioaccessibility values available for Ni. All blank extractions (n=4) returned As, Cr, and Ni values below the limit of quantification for this extraction type.For the CISED extraction, one sample was extracted in duplicate to assess the extraction repeatability and one blank extraction was carried out to assess any contribution made from the reagent blank. The mean repeatability across the 23 elements analysed was 4.65%. The blank extraction data indicates that the reagent blank did not contribute any elements of interest to the dataset as all results were below the limit of quantification for this extraction type.


## 3. Results and Discussion

### 3.1. Total PHE Concentrations

Table1 summarises the geology, lithology, sampling locations, organic matter content, soil pH, and the total and bioaccessible PHE concentrations of the test soils sampled from the Cherwell District. The organic matter content of the test soils ranged between 4.90 and 16.1% with a range of soil pHs from 5.83 to 7.54. The total As, Cr, and Ni ranged from 14.0 to 417 mg kg−1, from 51.0 to 447 mg kg−1, and from 17.0 to 218 mg kg−1, respectively. Mean As, Cr, and Ni concentrations were 150, 226, and 102 mg kg−1, an order of magnitude above the mean soil concentrations for these PHE reported by the UK soil herbage survey. Figures 2 and 3 are Box and Whisker plots summarizing the total PHE concentrations measured for this study, stratified by underlying geology and land use type, where “IO” and “LG” denote the Inferior Oolite and Lias Groups, respectively, and “A” and “R” denote allotment and residential, respectively. Figure 2 shows that the total concentrations of each PHE associated with the Lias and inferior Oolite Groups of the underlying geology is significantly different (one way ANOVA P values <0.001 for each element); in each case the PHE associated with the Lias group has the highest concentrations. Table 1 summarises the individual samples associated with each geological grouping and identifies Upper and Lower Heyford and Steeple Aston as overlaying the Inferior Oolite with associated comparatively low total PHE concentrations ranging between 14.0 and 50.0 mg kg−1, 51.0 and 98.0 mg kg−1, and 17.0 and 60.0 mg kg−1 for As, Cr, and Ni, respectively. In contrast, the total PHE concentrations measured for Drayton, Great Bourton, Shenington, and Wroxton sites (overlaying the Lias Group) are higher and have As, Cr, and Ni concentrations ranging from 82.0 to 417 mg kg−1, from 203 to 447 mg kg−1, and from 88.0 to 218 mg kg−1 (Figure 2). Total As, Cr, and Ni concentrations are not significantly different in the soil types when stratified by allotment or residential land use (Figure 3).Box and Whisker plots for total As, Cr, and Ni concentrations stratified by geology; “IO” represents the Inferior Oolite Group and “LG” represents the Lias Group. The bar shows the median value.
(a)
(b)
(c)Box and Whisker plots for total As, Cr, and Ni concentrations stratified by location type; “A” represents allotment land use and “R” represents residential land use. The bar shows the median value.
(a)
(b)
(c)The total As concentrations for the Cherwell District are in the same general concentration range of those previously measured by Palumbo-Roe et al. [10] who detected total As concentrations of up to 400 mg kg−1 and by Breward [29] who reported As concentrations of up to 340 mg kg−1 for the Claxby ironstones. However, the As concentrations in ironstone derived soils reported by Nathanail et al., Breward (Frodingham ironstones), and Wragg were markedly lower compared to the ranges reported here [22, 29–31]. Maximum total concentrations in the order of circa 100 mg kg−1 As were reported. In a study by Nathanail et al. [30] on As in Wellingborough allotments, the authors reported a range of total As concentrations between 39.0 and 113 mg kg−1 and a mean of 87.0 mg kg−1. For Cr, Breward [29] reported mean total concentrations of 100 mg kg−1 in topsoils in both the Claxby and Frodingham formations, with ranges of 16.0–254 mg kg−1 and 39.0–250 mg kg−1, respectively. Wragg [22] reported both total Cr and Ni concentrations, in ranges of 66.0–254 mg kg−1 and 20.0–88.0 mg kg−1, respectively. Palumbo-Roe et al. [10] summarised total Cr and Ni concentrations in soils from Scunthorpe, overlaying Jurassic ironstone, as containing 27.0–177 mg kg−1 and 31.0–140 mg kg−1, respectively. Comparison of the total Cr (51.0–447 mg kg−1) and Ni (17.0–218 mg kg−1) concentrations from this study shows that the highest total Cr concentrations measured are greater (circa double) than previously reported [10, 22, 29]. Similarly, the total Ni values determined for the Cherwell District are greater than those previously determined.


### 3.2. Bioaccessible PHE Concentrations

The PBET method used in this study has, more recently, been superseded by the Unified BARGE Method (UBM) which has been validated for As, Cd, and Pb against an animal model [32]. The justifications for using the PBET data for this specific investigation (the UBM had not been fully developed at the time of this study) were twofold as follows.(i)
The PBET test gave As and Cr bioaccessibility values for the BGS102 reference soil that were equivalent, within measurement error, with the consensus values for the reference soil that were determined using the UBM method. In addition, BGS 102 is an ironstone soil from the UK which is very similar in nature to the soil types being measured in Banbury area.(ii)
Recently published work on As bioaccessibility in the UK [33] has also compared the PBET with the UBM for bioaccessibility measurements on soil and found that they give comparable results for As in the UK.Whilst no validated methods for Ni and Cr are available, at the time of writing this paper, the results for bioaccessible Ni and Cr, as measured by the PBET method, are being used as guidance values for how the bioaccessibility of these two elements is likely to behave.The choice of which compartment (i.e., gastric or intestine) to report as containing the “bioaccessible fraction” has been discussed widely in the literature [14, 32, 34, 35]. The general consensus of opinion is that in order to maintain a conservative estimate, the highest value from the available compartments should be used. This approach has been adopted in this study and the following discussion is therefore aimed at a geochemical understanding of this higher amount.Figures4 and 5 are Box and Whisker plots summarizing the bioaccessible PHE concentrations, stratified by underlying geology and land use type, respectively. The mean As bioaccessibility obtained from all samples in this study (n=30) was 5.02 ± 2.48 mg kg−1. Broadly similar mean As bioaccessibilities have been reported for ironstone soils by workers, for example, Wragg et al. [11, 22] (mean 3.01 ± 1.37 mg kg−1; n=20), Nathanail et al. [30] (mean 2.95 ± 1.90 mg kg−1, n=70), and Palumbo-Roe et al. [10] (mean 3.85 ± 1.54 mg kg−1, n=70). Figure 4 indicates that, in general, there was a significant difference (P values for all elements <0.001) in the bioaccessible PHE concentrations observed between samples collected from the Inferior Oolite (n=12) and the Lias (n=18), with higher bioaccessible PHE concentrations associated with the Lias Group. Figure 5 shows that differences in land use did not indicate a significant difference in PHE bioaccessibility for any of the PHEs under investigation (P values for all three elements >0.05). Cr bioaccessibility for the residential sites only ranged from 0.60 to 1.38 mg kg−1, lower than that measured previously [22], where bioaccessible Cr was in the range of 1.74–6.00 mg kg−1 for other ironstone derived soils. The Ni bioaccessibility determined for the residential soils collected from the Cherwell District was similar to that determined for the Wellingborough garden soils [22] with ranges of 1.21–4.41 mg kg−1 and 1.59–3.36 mg kg−1, respectively.Box and Whisker plots for bioaccessible As, Cr, and Ni stratified by geology; “IO” represents the Inferior Oolite Group and “LG” represents the Lias Group. The bar shows the median value.
(a)
(b)
(c)Box and Whisker plots for bioaccessible As, Cr, and Ni stratified by location type; “A” represents allotment land use and “R” represents residential land use. The bar shows the median value.
(a)
(b)
(c)


### 3.3. Solid Phase Distribution of PHE

The CISED method has significant advantages over conventional extraction schemes. The use of a nonspecific extraction medium (aqua-regia) and the use of chemometric methods to separate out the different geochemical components significantly reduce the problems associated with the traditional methods which include the following: the nonspecificity of reagents; methods being operationally defined; analytical difficulties associated with the analysis of the extraction reagents; and the long reaction times and associated potential readsorption problems [20, 36]. There are a number of assumptions associated with the CISED test, which include the following.(i)
The relative solubility of the different constituents in the soils by the mineral acid extractant depends not only on the chemical make-up of the component (e.g., CaCO3, Fe2O3, MnO) but also on their physical form, for example, their particle size, degree of crystallinity, purity, and so forth. So, for example, crystalline forms of Fe2O3 may require higher acid strength for dissolution than fine grained Fe2O3; therefore, the components are defined by both their chemical and physical properties and hence are referred to as physicochemical components. The material under study consists of a mixture of these discrete physicochemical components with distinct major element compositions and the trace metals of interest are distributed amongst these components.(ii)
The physicochemical components will dissolve to different degrees and as the reagent strength increases, each solution will contain differing proportions of each of the components of the test material.(iii)
Within any given physicochemical component, all of the elements are dissolved congruently.(iv)
The CISED method will only identify geochemical components that are soluble in aqua regia.Whilst the complex mixture of natural soils will not exactly conform to these assumptions, the CISED offers very important improvements over conventional methods. Point (iii) is an assumption that is made in all sequential extraction schemes and has its limitations as it cannot always be guaranteed. For the CISED method, this would not be met for the bulk aluminosilicate minerals of the soil (these do not dissolve well in HCl/HNO3 mixtures); however, the method is designed to investigate the more reactive surface coatings on the soil particulates, for example, Ca carbonate, Fe oxides, and Mn oxides, where the most mobile forms of the metals are held. These more reactive physicochemical components are well solubilised in the HCl/HNO3 acid medium and, therefore, the assumption in point (iii) can be taken to be true for the CISED method.Four of the test samples were subjected to the CISED extraction methodology to determine the physicochemical soil components present across the range of test locations, for example, carbonates, clays, exchangeable phases, and so forth, and the solid phase distribution of As, Cr, and Ni between each identified component. Extraction recoveries for the CISED extraction method ranged from 5.71 to 8.43%, from 4.24 to 9.07%, and from 9.32 to 23.1% for As, Cr, and Ni, respectively. The recovery data alone provides an indication of the mobility of each PHE and shows that Ni is more available to acid extraction (up to 5.0 molar aqua regia) than As or Cr, which are extracted. Application of the SMMR algorithm [37] to the soil extraction data resulted in the determination of a number of physicochemical components for each soil under test. Hierarchical cluster analysis of the mean, centred, and scaled soil extraction data resulted in the production of a colourmap, a visual aid to identify the physicochemical components in the test soils (Figure 6). Each row of Figure 6 represents a physicochemical component found within a given soil; the soil sample is indicated by the first letters which refer to the name codes given in Table 1 and the following element names indicate which elements make up 10% or more of the physicochemical component composition. The element columns of Figure 6 (right-hand side of the figure) represent the relative amounts of each element present in the component and the E1 to E14 columns represent the relative amounts of the component extracted in each extraction step (increasing acid strength from E1 to E14). Figure 6 indicates the presence of 7 distinct physicochemical clusters ordered in relation to their ease of extractability from the test soils: Organic/Residual pore-water; Fertiliser; Carbonate; Mn oxide; High Carbonate; Al oxide; and Fe oxide.Figure 6
Colourmap and associated clustergram for CISED extract data for a selection of the test soils. The dendrogram along the left-hand side shows how individual components link together. The horizontal line divides the colourmap into clusters and the vertical line separates the element compositions data on the left-hand side from the extraction number data on the right. A high proportion of each component and an indication of its composition are shown by a white or pale grey colouration with a low proportion as dark grey or black.Organic/Residual pore-water: the cluster identified as Organic/Residual pore-water was extracted, primarily, during the initial extraction steps of the CISED, steps (1) and (2), using deionised water. This cluster, which was dominated by the presence of Na, K, and Si, is highly mobile and can be derived from soluble pore-water salts in the soils and possibly soluble humic material which is also water soluble and can contain these elements [38]. This component is associated with three of the four samples under investigation, the allotment and residential samples from Shenington and the allotment sample from Drayton.Fertiliser: this cluster was extracted by deionised water predominantly under low acid concentrations (0.01–0.05 M aqua regia) in a narrow window of extraction and associated with the two allotment samples (Shenington and Drayton) and one of the two residential samples (Upper Heyford). This cluster was dominated by the presence of a number of the primary and secondary plant nutrients used in fertilisers K, P, and Mg, respectively.Carbonate: a Ca dominated cluster, which also had a contribution from the presence of Mg, was extracted on the first addition of aqua regia and over a relatively narrow window (0.01–0.05 M). This cluster was identified with all samples subjected to the CISED sequential extraction methodology.Mn oxide: a cluster dominated by the presence of Mn was extracted in a distinct window of extraction, associated with three of the four samples under investigation (Shenington (residential and allotment) and Drayton) and observed at the point where H2O2 was first added into the test system. It is known that Mn oxides are readily dissolved in H2O2 [39], suggesting that this cluster was probably derived from the dissolution of a naturally occurring Mn oxide.High carbonate: a cluster that was distinctly different from that identified as carbonate was named high carbonate. The differences between the two carbonate clusters were the increased amount of high carbonate identified (7000 mg kg−1), the difference in its composition (Ca only) compared to the carbonate phase, its different window of extraction (0.10 M–5.0 M), and its association with only one of the test samples (Upper Heyford). This soil component could be associated with the presence of calcite in this individual soil, which has been similarly reported by previous studies of Jurassic ironstones [10, 11, 22].Al oxide: similarly to the Mn oxide cluster, a cluster dominated by the presence of Al is first seen on or around the first addition of H2O2, a reagent known to aid the dissolution of Al oxides [39].Fe oxide: these clusters are most probably derived from relatively pure Fe oxides, because of their high Fe contents and because high acid concentrations are required to extract them from the soils [40]. All soils under investigation are associated with this cluster, which is extracted at medium to high acid concentrations (0.50–5.0 M aqua regia) and over a broad window of extraction.


### 3.4. Relationship between Distribution of PHEs and Bioaccessibility

The CISED procedure has now identified the main easily extractable physicochemical components in the soils and the algorithm also gives the amount of each PHE associated with these components. In a similar manner to previous studies [11, 24, 36], the bioaccessibility of the PHE of interest has been related to the distribution of the PHE in the individual test soils subjected to the CISED in order to gain an understanding about the sources of contaminant bioaccessibility. Figure 7 summarises the results showing the cumulative amount of each PHE extracted going from the most easily extracted to those only extracted at higher acid concentration with the bioaccessible fraction marked as a horizontal dashed line (the total concentration of each metal is also noted on each plot but not shown graphically as these values are very much higher than the bioaccessible and CISED extracted values). In Figure 7, the cumulative amount of PHE extracted up to the point at which the bioaccessibility line intersects the cumulative mass extracted indicates which physic-chemical components are likely to be associated with the bioaccessible fraction of the PHE. Examination of Figure 7 shows that each element follows the same relationship between cumulative mass extracted and the bioaccessible fraction in each sample:(i)
for As, the bioaccessible fraction comes from the organic, fertiliser, carbonate, Mn oxide, and Al oxide fractions with some partial dissolution of the Fe oxide;(ii)
for Cr, the bioaccessible fraction comes from the organic, fertiliser, carbonate, Mn oxide, and Al oxide fractions but not from the Fe oxide fraction;(iii)
for Ni, the bioaccessible fraction comes almost exclusively from the carbonate fraction.Comparisons of the cumulative mass of As, Cr, and Ni extracted in the CISED identified physicochemical components and its relationship to the bioaccessible fractions. Solid line represents the cumulative mass extracted; the dashed line represents the bioaccessible fraction. Total concentrations of elements are noted on the plot.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)


## 3.1. Total PHE Concentrations

Table1 summarises the geology, lithology, sampling locations, organic matter content, soil pH, and the total and bioaccessible PHE concentrations of the test soils sampled from the Cherwell District. The organic matter content of the test soils ranged between 4.90 and 16.1% with a range of soil pHs from 5.83 to 7.54. The total As, Cr, and Ni ranged from 14.0 to 417 mg kg−1, from 51.0 to 447 mg kg−1, and from 17.0 to 218 mg kg−1, respectively. Mean As, Cr, and Ni concentrations were 150, 226, and 102 mg kg−1, an order of magnitude above the mean soil concentrations for these PHE reported by the UK soil herbage survey. Figures 2 and 3 are Box and Whisker plots summarizing the total PHE concentrations measured for this study, stratified by underlying geology and land use type, where “IO” and “LG” denote the Inferior Oolite and Lias Groups, respectively, and “A” and “R” denote allotment and residential, respectively. Figure 2 shows that the total concentrations of each PHE associated with the Lias and inferior Oolite Groups of the underlying geology is significantly different (one way ANOVA P values <0.001 for each element); in each case the PHE associated with the Lias group has the highest concentrations. Table 1 summarises the individual samples associated with each geological grouping and identifies Upper and Lower Heyford and Steeple Aston as overlaying the Inferior Oolite with associated comparatively low total PHE concentrations ranging between 14.0 and 50.0 mg kg−1, 51.0 and 98.0 mg kg−1, and 17.0 and 60.0 mg kg−1 for As, Cr, and Ni, respectively. In contrast, the total PHE concentrations measured for Drayton, Great Bourton, Shenington, and Wroxton sites (overlaying the Lias Group) are higher and have As, Cr, and Ni concentrations ranging from 82.0 to 417 mg kg−1, from 203 to 447 mg kg−1, and from 88.0 to 218 mg kg−1 (Figure 2). Total As, Cr, and Ni concentrations are not significantly different in the soil types when stratified by allotment or residential land use (Figure 3).Box and Whisker plots for total As, Cr, and Ni concentrations stratified by geology; “IO” represents the Inferior Oolite Group and “LG” represents the Lias Group. The bar shows the median value.
(a)
(b)
(c)Box and Whisker plots for total As, Cr, and Ni concentrations stratified by location type; “A” represents allotment land use and “R” represents residential land use. The bar shows the median value.
(a)
(b)
(c)The total As concentrations for the Cherwell District are in the same general concentration range of those previously measured by Palumbo-Roe et al. [10] who detected total As concentrations of up to 400 mg kg−1 and by Breward [29] who reported As concentrations of up to 340 mg kg−1 for the Claxby ironstones. However, the As concentrations in ironstone derived soils reported by Nathanail et al., Breward (Frodingham ironstones), and Wragg were markedly lower compared to the ranges reported here [22, 29–31]. Maximum total concentrations in the order of circa 100 mg kg−1 As were reported. In a study by Nathanail et al. [30] on As in Wellingborough allotments, the authors reported a range of total As concentrations between 39.0 and 113 mg kg−1 and a mean of 87.0 mg kg−1. For Cr, Breward [29] reported mean total concentrations of 100 mg kg−1 in topsoils in both the Claxby and Frodingham formations, with ranges of 16.0–254 mg kg−1 and 39.0–250 mg kg−1, respectively. Wragg [22] reported both total Cr and Ni concentrations, in ranges of 66.0–254 mg kg−1 and 20.0–88.0 mg kg−1, respectively. Palumbo-Roe et al. [10] summarised total Cr and Ni concentrations in soils from Scunthorpe, overlaying Jurassic ironstone, as containing 27.0–177 mg kg−1 and 31.0–140 mg kg−1, respectively. Comparison of the total Cr (51.0–447 mg kg−1) and Ni (17.0–218 mg kg−1) concentrations from this study shows that the highest total Cr concentrations measured are greater (circa double) than previously reported [10, 22, 29]. Similarly, the total Ni values determined for the Cherwell District are greater than those previously determined.


## 3.2. Bioaccessible PHE Concentrations

The PBET method used in this study has, more recently, been superseded by the Unified BARGE Method (UBM) which has been validated for As, Cd, and Pb against an animal model [32]. The justifications for using the PBET data for this specific investigation (the UBM had not been fully developed at the time of this study) were twofold as follows.(i)
The PBET test gave As and Cr bioaccessibility values for the BGS102 reference soil that were equivalent, within measurement error, with the consensus values for the reference soil that were determined using the UBM method. In addition, BGS 102 is an ironstone soil from the UK which is very similar in nature to the soil types being measured in Banbury area.(ii)
Recently published work on As bioaccessibility in the UK [33] has also compared the PBET with the UBM for bioaccessibility measurements on soil and found that they give comparable results for As in the UK.Whilst no validated methods for Ni and Cr are available, at the time of writing this paper, the results for bioaccessible Ni and Cr, as measured by the PBET method, are being used as guidance values for how the bioaccessibility of these two elements is likely to behave.The choice of which compartment (i.e., gastric or intestine) to report as containing the “bioaccessible fraction” has been discussed widely in the literature [14, 32, 34, 35]. The general consensus of opinion is that in order to maintain a conservative estimate, the highest value from the available compartments should be used. This approach has been adopted in this study and the following discussion is therefore aimed at a geochemical understanding of this higher amount.Figures4 and 5 are Box and Whisker plots summarizing the bioaccessible PHE concentrations, stratified by underlying geology and land use type, respectively. The mean As bioaccessibility obtained from all samples in this study (n=30) was 5.02 ± 2.48 mg kg−1. Broadly similar mean As bioaccessibilities have been reported for ironstone soils by workers, for example, Wragg et al. [11, 22] (mean 3.01 ± 1.37 mg kg−1; n=20), Nathanail et al. [30] (mean 2.95 ± 1.90 mg kg−1, n=70), and Palumbo-Roe et al. [10] (mean 3.85 ± 1.54 mg kg−1, n=70). Figure 4 indicates that, in general, there was a significant difference (P values for all elements <0.001) in the bioaccessible PHE concentrations observed between samples collected from the Inferior Oolite (n=12) and the Lias (n=18), with higher bioaccessible PHE concentrations associated with the Lias Group. Figure 5 shows that differences in land use did not indicate a significant difference in PHE bioaccessibility for any of the PHEs under investigation (P values for all three elements >0.05). Cr bioaccessibility for the residential sites only ranged from 0.60 to 1.38 mg kg−1, lower than that measured previously [22], where bioaccessible Cr was in the range of 1.74–6.00 mg kg−1 for other ironstone derived soils. The Ni bioaccessibility determined for the residential soils collected from the Cherwell District was similar to that determined for the Wellingborough garden soils [22] with ranges of 1.21–4.41 mg kg−1 and 1.59–3.36 mg kg−1, respectively.Box and Whisker plots for bioaccessible As, Cr, and Ni stratified by geology; “IO” represents the Inferior Oolite Group and “LG” represents the Lias Group. The bar shows the median value.
(a)
(b)
(c)Box and Whisker plots for bioaccessible As, Cr, and Ni stratified by location type; “A” represents allotment land use and “R” represents residential land use. The bar shows the median value.
(a)
(b)
(c)


## 3.3. Solid Phase Distribution of PHE

The CISED method has significant advantages over conventional extraction schemes. The use of a nonspecific extraction medium (aqua-regia) and the use of chemometric methods to separate out the different geochemical components significantly reduce the problems associated with the traditional methods which include the following: the nonspecificity of reagents; methods being operationally defined; analytical difficulties associated with the analysis of the extraction reagents; and the long reaction times and associated potential readsorption problems [20, 36]. There are a number of assumptions associated with the CISED test, which include the following.(i)
The relative solubility of the different constituents in the soils by the mineral acid extractant depends not only on the chemical make-up of the component (e.g., CaCO3, Fe2O3, MnO) but also on their physical form, for example, their particle size, degree of crystallinity, purity, and so forth. So, for example, crystalline forms of Fe2O3 may require higher acid strength for dissolution than fine grained Fe2O3; therefore, the components are defined by both their chemical and physical properties and hence are referred to as physicochemical components. The material under study consists of a mixture of these discrete physicochemical components with distinct major element compositions and the trace metals of interest are distributed amongst these components.(ii)
The physicochemical components will dissolve to different degrees and as the reagent strength increases, each solution will contain differing proportions of each of the components of the test material.(iii)
Within any given physicochemical component, all of the elements are dissolved congruently.(iv)
The CISED method will only identify geochemical components that are soluble in aqua regia.Whilst the complex mixture of natural soils will not exactly conform to these assumptions, the CISED offers very important improvements over conventional methods. Point (iii) is an assumption that is made in all sequential extraction schemes and has its limitations as it cannot always be guaranteed. For the CISED method, this would not be met for the bulk aluminosilicate minerals of the soil (these do not dissolve well in HCl/HNO3 mixtures); however, the method is designed to investigate the more reactive surface coatings on the soil particulates, for example, Ca carbonate, Fe oxides, and Mn oxides, where the most mobile forms of the metals are held. These more reactive physicochemical components are well solubilised in the HCl/HNO3 acid medium and, therefore, the assumption in point (iii) can be taken to be true for the CISED method.Four of the test samples were subjected to the CISED extraction methodology to determine the physicochemical soil components present across the range of test locations, for example, carbonates, clays, exchangeable phases, and so forth, and the solid phase distribution of As, Cr, and Ni between each identified component. Extraction recoveries for the CISED extraction method ranged from 5.71 to 8.43%, from 4.24 to 9.07%, and from 9.32 to 23.1% for As, Cr, and Ni, respectively. The recovery data alone provides an indication of the mobility of each PHE and shows that Ni is more available to acid extraction (up to 5.0 molar aqua regia) than As or Cr, which are extracted. Application of the SMMR algorithm [37] to the soil extraction data resulted in the determination of a number of physicochemical components for each soil under test. Hierarchical cluster analysis of the mean, centred, and scaled soil extraction data resulted in the production of a colourmap, a visual aid to identify the physicochemical components in the test soils (Figure 6). Each row of Figure 6 represents a physicochemical component found within a given soil; the soil sample is indicated by the first letters which refer to the name codes given in Table 1 and the following element names indicate which elements make up 10% or more of the physicochemical component composition. The element columns of Figure 6 (right-hand side of the figure) represent the relative amounts of each element present in the component and the E1 to E14 columns represent the relative amounts of the component extracted in each extraction step (increasing acid strength from E1 to E14). Figure 6 indicates the presence of 7 distinct physicochemical clusters ordered in relation to their ease of extractability from the test soils: Organic/Residual pore-water; Fertiliser; Carbonate; Mn oxide; High Carbonate; Al oxide; and Fe oxide.Figure 6
Colourmap and associated clustergram for CISED extract data for a selection of the test soils. The dendrogram along the left-hand side shows how individual components link together. The horizontal line divides the colourmap into clusters and the vertical line separates the element compositions data on the left-hand side from the extraction number data on the right. A high proportion of each component and an indication of its composition are shown by a white or pale grey colouration with a low proportion as dark grey or black.Organic/Residual pore-water: the cluster identified as Organic/Residual pore-water was extracted, primarily, during the initial extraction steps of the CISED, steps (1) and (2), using deionised water. This cluster, which was dominated by the presence of Na, K, and Si, is highly mobile and can be derived from soluble pore-water salts in the soils and possibly soluble humic material which is also water soluble and can contain these elements [38]. This component is associated with three of the four samples under investigation, the allotment and residential samples from Shenington and the allotment sample from Drayton.Fertiliser: this cluster was extracted by deionised water predominantly under low acid concentrations (0.01–0.05 M aqua regia) in a narrow window of extraction and associated with the two allotment samples (Shenington and Drayton) and one of the two residential samples (Upper Heyford). This cluster was dominated by the presence of a number of the primary and secondary plant nutrients used in fertilisers K, P, and Mg, respectively.Carbonate: a Ca dominated cluster, which also had a contribution from the presence of Mg, was extracted on the first addition of aqua regia and over a relatively narrow window (0.01–0.05 M). This cluster was identified with all samples subjected to the CISED sequential extraction methodology.Mn oxide: a cluster dominated by the presence of Mn was extracted in a distinct window of extraction, associated with three of the four samples under investigation (Shenington (residential and allotment) and Drayton) and observed at the point where H2O2 was first added into the test system. It is known that Mn oxides are readily dissolved in H2O2 [39], suggesting that this cluster was probably derived from the dissolution of a naturally occurring Mn oxide.High carbonate: a cluster that was distinctly different from that identified as carbonate was named high carbonate. The differences between the two carbonate clusters were the increased amount of high carbonate identified (7000 mg kg−1), the difference in its composition (Ca only) compared to the carbonate phase, its different window of extraction (0.10 M–5.0 M), and its association with only one of the test samples (Upper Heyford). This soil component could be associated with the presence of calcite in this individual soil, which has been similarly reported by previous studies of Jurassic ironstones [10, 11, 22].Al oxide: similarly to the Mn oxide cluster, a cluster dominated by the presence of Al is first seen on or around the first addition of H2O2, a reagent known to aid the dissolution of Al oxides [39].Fe oxide: these clusters are most probably derived from relatively pure Fe oxides, because of their high Fe contents and because high acid concentrations are required to extract them from the soils [40]. All soils under investigation are associated with this cluster, which is extracted at medium to high acid concentrations (0.50–5.0 M aqua regia) and over a broad window of extraction.


## 3.4. Relationship between Distribution of PHEs and Bioaccessibility

The CISED procedure has now identified the main easily extractable physicochemical components in the soils and the algorithm also gives the amount of each PHE associated with these components. In a similar manner to previous studies [11, 24, 36], the bioaccessibility of the PHE of interest has been related to the distribution of the PHE in the individual test soils subjected to the CISED in order to gain an understanding about the sources of contaminant bioaccessibility. Figure 7 summarises the results showing the cumulative amount of each PHE extracted going from the most easily extracted to those only extracted at higher acid concentration with the bioaccessible fraction marked as a horizontal dashed line (the total concentration of each metal is also noted on each plot but not shown graphically as these values are very much higher than the bioaccessible and CISED extracted values). In Figure 7, the cumulative amount of PHE extracted up to the point at which the bioaccessibility line intersects the cumulative mass extracted indicates which physic-chemical components are likely to be associated with the bioaccessible fraction of the PHE. Examination of Figure 7 shows that each element follows the same relationship between cumulative mass extracted and the bioaccessible fraction in each sample:(i)
for As, the bioaccessible fraction comes from the organic, fertiliser, carbonate, Mn oxide, and Al oxide fractions with some partial dissolution of the Fe oxide;(ii)
for Cr, the bioaccessible fraction comes from the organic, fertiliser, carbonate, Mn oxide, and Al oxide fractions but not from the Fe oxide fraction;(iii)
for Ni, the bioaccessible fraction comes almost exclusively from the carbonate fraction.Comparisons of the cumulative mass of As, Cr, and Ni extracted in the CISED identified physicochemical components and its relationship to the bioaccessible fractions. Solid line represents the cumulative mass extracted; the dashed line represents the bioaccessible fraction. Total concentrations of elements are noted on the plot.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)


## 4. Conclusions

A combination of total element analysis, bioaccessibility measurements, and sequential extraction has proved to be a powerful combination of techniques to understand the mobility of PHE in soils derived from Jurassic ironstone in the UK. Whilst the total concentration of As, Cr, and Ni is elevated due to the soil parent material, bioaccessibility testing, which mimics the conditions in the human gastrointestinal system, shows that only a small proportion of the total concentration is available for absorption into the human body. The CISED test shows that the nonmobile forms of the PHEs are strongly sorbed to iron oxides and are unreactive. The study also shows how parent material geology has a significant effect on the total element concentrations and the bioaccessibility of PHE. Land use (allotment versus residential) does not show such a large effect but the allotment bioaccessibility data show a bigger spread and possibly higher values for As and Cr which may be due to cultivation practices such as addition of fertilisers and organic matter.



---

*Source: 924891-2014-03-05.xml*