Pot experiment

A pot experiment was conducted during 32 weeks at Wageningen, the Netherlands. The seven treatments (Table 1) in four replicates (blocks) included five doses of olivine powder (including zero). Olivine doses increased fivefold each next level. Two levels of kieserite, a highly soluble Mg sulphate fertiliser, were included. We used kieserite as a reference expected to provide its full Mg content (16.2% on mass basis) as bioavailable Mg right from the start of the experiment. Our forsterite-dominated olivine product from ‘North Cape Olivine Sand’ (Sibelco Nordic Ltd.) contained 23.4 mass % Mg, and 4.0 mass % Fe (this corresponds to a molar Fe∶Mg ratio of 1∶13.4). See also Table S1 for chemical composition. The finely-ground olivine product consisted of 7 mass-% of particles <2 µm, 66% between 2 and 50 µm, and 27% between 50 and 200 µm. The full particle size distribution curve is given in Figure S1. The olivine powder was further characterised by a BET (Brunauer, Emmett and Teller) specific surface area of 4.8 m² g⁻¹, determined by N₂ adsorption [11] at 77 K on a Micromeritics TriStar 3000 analyzer. (Prior to adsorption measurement, the olivine sample was degassed in vacuum for 16 hours at 120°C and 350°C, respectively. The two temperatures yielded the same BET areas.)

Download:

• PPT
  PowerPoint slide
• PNG
  larger image
• TIFF
  original image

Table 1. Treatments in pot experiment, Wageningen, 2009–2010.

https://doi.org/10.1371/journal.pone.0042098.t001

The olivine powder was mixed with sandy soil collected from a field under arable cultivation at Droevendaal farm near Wageningen. The soil was selected for its relatively low bioavailability (among Dutch arable soils) of Mg (44.4 mg kg⁻¹) and K (41.7 mg kg⁻¹). Further chemical characteristics of the soil are given in Table S2.

PVC pots of 10 litre, with top diameter of 250 mm, were filled with the soil-olivine mixture. Pots contained 11.0 kg dry soil, except at the highest olivine rate where 9.5 kg soil was used. We used perennial ryegrass (Lolium perenne, cv. Barata) as the test crop, for its capacity to absorb large amounts of nutrients over a long growth period, and for the elasticity of nutrient contents in its biomass. The grass was sown on 20 August 2009, at 0.3 g seed per pot. At the start, all pots were supplied with adequate initial amounts of nitrogen (N), phosphorus (P), and potassium (K), equivalent with rates of 80 kg N ha⁻¹, 19.6 kg P ha⁻¹, and 80 kg K ha⁻¹. Fertilisers were finely ground and mixed through the upper half of soil in the pots. To sustain high biomass production, extra N, P and K were supplied after each harvest as finely ground fertilisers on top of the remaining grass stubble, and flushed into the soil with irrigation water. See Table S3 for a full account on fertiliser management. To avoid the possibility of nutrient losses via downward percolation, irrigation water was supplied daily onto a tray below each pot, for a predominantly upward flow. Once every 14 days, 400 ml water was supplied on top to prevent accumulation of solutes near the surface.

Pots were kept outdoors under transparent rain shelter, until 3 November. They were then transferred to a greenhouse to sustain crop growth and soil processes during winter and early spring. The temperature regime for both periods is given in Figures S2 and S3. Always leaving the grass stubble for regrowth, we harvested aboveground plant biomass on 28 September, 2 November, 3 December, 11 January, 22 February, and 31 March. Fresh biomass yield was measured, and samples were dried at 70°C to assess dry matter yield. The experiment ended on 7 April, 2010. The soil from each pot was then thoroughly mixed and sampled for chemical analysis.

Plant samples were analysed separately (per harvest, per pot) for N, P, K, Ca, Mg, S, Ni, Fe, and Si. Soil was analysed for total and bioavailable contents of these elements, at the end of the experiment. Methods of soil destruction, element extraction (bioavailability), and elements analysis in soil and crop are given in Table S4. Soil water samples were collected only occasionally, and from selected treatments. This was done by suction using a 100 mm long synthetic microporous tube, embedded horizontally halfway between top and bottom in the centre of each pot. Water samples were analysed for pH and total alkalinity, for Mg and Ni (two occasions), and Si (once).

Estimation of the fraction of olivine weathered

We estimated olivine weathering in our experiment by two methods. Method 1 is based on the Mg balance: the amount of olivine weathered corresponds to the amount of bioavailable Mg accumulated in soil (Mg_bio,soil, g pot⁻¹) and plant biomass (Mg_plant, g pot⁻¹), in excess of that in the Control. The mass fraction of olivine weathered (F_weath,Meth1) relative to olivine applied, is then written as:(2)where the subscripts of groups in brackets refer to treatments with olivine (OLIV) and control (CON). The symbol Mg_app,_OLIV refers to the amount of Mg applied in the form of olivine (g pot⁻¹).

Alternatively (Method 2), the kieserite treatments were used as references. This enables to express an olivine dose as an equivalent kieserite dose, that is, a dose that has the same impact on a particular response variable, X. Thus, we can calculate by interpolation (between Control and KIES1; or between KIES1 and KIES2) or by extrapolation (beyond KIES2) how much Mg in kieserite form was needed to achieve the effect on X found in OLIV1. Presuming full dissolution of kieserite, then, the same amount of Mg must have been dissolved from olivine. For example, if the response variable X in treatment OLIV1 is between the X values found in treatments KIES1 and KIES2, respectively, then the mass fraction of olivine weathered (F_weath,Meth2) in OLIV1 is calculated as:(3)where subscripts again refer to treatments, and Mg_app are amounts (g pot⁻¹) of magnesium applied in the respective treatments (KIES1, KIES2, OLIV1). In contrast to Method 1, this Method 2 does not rely on a complete balance, but on similarity of fate of dissolved Mg across treatments, irrespective of its origin (olivine or kieserite). We applied Method 2 to three response variables (X in Eq.3): Mg concentration in plant biomass, total Mg offtake in plant biomass, and bioavailable soil Mg.

Statistical analysis

Analysis of variance was used to analyse the randomized block experiment. Residual plots revealed that for most response variables a log-transformation stabilized the variances. Therefore all response variables, except soil water pH, were logarithmically transformed prior to analysis. Each pair of treatments means was tested for significance at the 1% level using a Student t-test which employs the residual mean square of the analysis of variance as an estimate of the variance, see e.g. [12]. The results of these pairwise tests are presented by annotating mean values, given on the original scale, with a superscript letter such that treatments sharing the same letter are not significantly different at the 1% level, while treatment means with no letter in common are significantly different. We verified that only rare and minor discrepancies occurred between statistical analysis on the original and on the log-scale.

Harvested biomass and nutrient concentrations in plant tissue were measured per separate harvest event. Their values were first aggregated to total biomass yield and to mass-weighted average nutrient concentrations, prior to statistical analysis.

Olsen model for olivine weathering

Estimates of weathering from our pot experiment were compared with reaction kinetics based on laboratory measurements under a wider range of conditions. We used a simple model based on [10], [13], expressing olivine weathering per unit area of crystal surface in terms of the ‘carbonation rate’ r (mol m⁻² s⁻¹):(4a)(4b)

These regression equations were based by [10] on geometrical surface area for spherical particles, as opposed to BET surface area. (The two approaches give different estimates of the regression parameters, see [10] for a comparative study on a large set of laboratory data.) The dependence on temperature (T) is given by the Arrhenius equation:(5)with T_ref for reference temperature (298 K), and E_a for activation energy (63 kJ/mol; [14]). R is the Universal gas constant (8,31 J K⁻¹ mol⁻¹) and r_Tref is the carbonation rate at reference temperature. Using 4.4 10⁻⁵ m³ mol⁻¹ for the molar volume of olivine, V_m [15], the weathering rate per unit surface area is converted into a corresponding retreat of the reactive surface position (‘shrinking particle model’). Following [14], the model thus accounts for the time t_d that is required to completely dissolve a particle, in function of its size class. For particles of diameter D₀ (m), t_d (s) is approximated [10] as(6)To implement the calculations, we expressed the particle size distribution in 95 size classes. The model was applied at the lowest and highest soil water pH observed in our Control treatment, respectively. Model results are given as total mass fraction weathered, as well as contributions from selected particle size classes.

Plant biomass and plant analysis

Plant growth did not differ from the Control in any treatment, except OLIV4 which gave a 15.6% higher dry matter (DM) yield (Table S5; Treatment codes in Table 1). Expressed per unit pot surface area, DM yields corresponded to 19.2 to 23.6 10³ kg ha⁻¹. These values are roughly 50% higher than typically obtained under field conditions, light interception in pot culture being larger than in the field. Total growth duration, however, was comparable to field conditions in North-West Europe. The yields obtained are proof of vigorous growth in all treatments.

Element concentrations in plant dry matter are given in Table 2. N concentration was not affected by olivine addition. Values are considered low - relative to values between 25 and 35 g kg⁻¹ for normal production conditions. They show that fertiliser N application at the equivalent of 480 kg N ha⁻¹ (including topdressing after each harvest) was modest, relative to crop demand. P concentrations were within the range of 3.0 to 4.5 g kg⁻¹ typically found, and were somewhat reduced by olivine. Total P uptake, however, was not. K uptake, in contrast, was increased by olivine, but only at its highest dose (OLIV4). Extra K uptake by plants was possibly due to preferential adsorption of Mg⁺⁺ (over K⁺) on the soil complex, releasing K⁺ into solution. Mg, Si and Ni concentrations in plant dry matter were higher than in the Control, even at the smallest olivine dose, and effects increased with higher doses. Ca, in contrast, decreased with larger olivine doses. This is attributed to competition between Ca and Mg uptake, and shows that olivine might induce a nutritional imbalance as it does in natural systems with high inherent Mg/Ca ratios [16], [17], [18], [19]. There was a significant rise in Ni concentration in grass at all doses of olivine, with a fivefold increment from 531 (Control) to 2669 µg per kg dry biomass (OLIV4).

Download:

• PPT
  PowerPoint slide
• PNG
  larger image
• TIFF
  original image

Table 2. Element concentrations¹ in plant dry matter, mass-weighted average over all six harvests.

https://doi.org/10.1371/journal.pone.0042098.t002

Soil analysis

Bioavailability of elements in soil at the end of the experiment is listed in Table 3, only for elements that showed significant responses. Available soil Mg, Si, and Ni, and soil pH were all increased by olivine, and effects increased with larger doses. Soil pH increased from 4.89 in the Control to 5.96 at the largest olivine dose. While changes in Mg, Si and Ni may directly reflect the addition of these elements (with olivine), some other elements were affected, too. Available K increased while total soil K remained unaffected. Available P decreased slightly with higher olivine dose, but only in one extractant (0.01 M CaCl₂, Table 3), not in ammoniumlactate-acetic acid (P_AL). (P_AL in the Control was 24.8 mg per 100 g soil). These responses were possibly related to changes in soil pH. Inorganic soil N was not affected. Its value in the Control treatment was 15.00 mg per kg dry soil.

Download:

• PPT
  PowerPoint slide
• PNG
  larger image
• TIFF
  original image

Table 3. Bioavailability¹ of elements in soil at the end of the experiment.

https://doi.org/10.1371/journal.pone.0042098.t003

Increments in total soil Mg and total soil Ni reflect the total amounts of olivine added, because all olivine was dissolved by sample destruction with Aqua regia. For these and other elements with significant responses, total concentrations in soil are given in Table S6. Total N, Ca, K and S in soil were not affected by olivine. Their values in the Control treatment were 1.15 g N, 923 mg Ca, 355 mg K, and 198 mg S per kg dry soil.

Soil water analysis

Element concentrations in soil water were measured at two occasions in selected treatments (Table S7). Concentrations of Mg, Ni and Si in water were higher in OLIV2, OLIV3 and OLIV4 than in the Control. Soil water pH and alkalinity were measured in selected treatments (Table 4). All values refer to soil water at a harvest event (grass cutting), prior to dressing of new fertiliser for regrowth. Throughout the experiment, pH and alkalinity were higher in OLIV4 than in the Control treatment. Only on April 7 were soil water pH and alkalinity measured in all treatments. Soil water pH, then, was increased only in OLIV4, while alkalinity was increased in OLIV3 and OLIV4. The fluctuations in pH and alkalinity over time, and the absence of pH and alkalinity responses at lower olivine doses, are not fully understood. They may reflect buffering by the soil system. Also, the application of fertilisers, root activity, and fluctuations in soil water content affect the inorganic C balance. For these reasons we did not attempt to estimate olivine weathering or net CO₂ sequestration from changes in alkalinity.

Download:

• PPT
  PowerPoint slide
• PNG
  larger image
• TIFF
  original image

Table 4. Soil water pH¹ and alkalinity¹.

https://doi.org/10.1371/journal.pone.0042098.t004

Olivine weathering in the experiment

Olivine weathering estimated by Method 1 is shown in Table 5. At the lowest olivine dose (OLIV1), 14.8% of applied olivine Mg became bioavailable during the 32 weeks of the trial. This fraction decreased steeply with higher olivine doses. According to Method 2 (Table 6), estimates of olivine weathered in treatment OLIV1 ranged from 13% to 58%, depending on the indicator chosen. The lower value refers to Mg content in plant biomass, the upper value to bioavailable soil Mg.

Download:

• PPT
  PowerPoint slide
• PNG
  larger image
• TIFF
  original image

Table 5. Magnesium (Mg) balance terms, and ratio of bioavailable to applied Mg (Method 1).

https://doi.org/10.1371/journal.pone.0042098.t005

Download:

• PPT
  PowerPoint slide
• PNG
  larger image
• TIFF
  original image

Table 6. Fraction of olivine weathered in OLIV1, estimated via Method 2 (similarity with kieserite).

https://doi.org/10.1371/journal.pone.0042098.t006

Estimation of olivine weathering

Method 1 presumes that all Mg released from olivine is retrieved in soil and plant sampling. This, however, is challenged by the only partial retrieval of Mg from kieserite (30%, Table 5). While kieserite is considered highly soluble, Mg from kieserite then must have precipitated in a form not extracted as bioavailable (in 0.01 M CaCl₂); or it was lost from the system. The former option seems more likely, as cations are generally held by the soil's negatively charged adsorption complex and, moreover, we aimed to avoid net water percolation. The same fate, then, may have applied to Mg released from olivine, and so our estimation of weathering by Method 1 may be too low. (Amounts of Mg in soil water are ignored as a balance term. These were two orders of magnitude smaller than bioavailable soil Mg; divalent cations are strongly adsorbed by soil). For Method 2, while the two lower values in Table 6 (column f) were obtained by interpolation within the data range, the highest value was extrapolated beyond the Mg dose of KIES2. Without extrapolation, the equivalency of the olivine dose in OLIV1 must be set to 0.486 g kieserite-Mg pot⁻¹ (as in KIES2), and the corresponding estimate of the fraction weathered becomes 26%.

Weathering in experiment versus Olsen model

The Olsen model indicated that for our experiment, with monthly average temperatures as measured and with pH fixed at 5.4, a skin of thickness 0.12 µm on olivine particles would react during the trial. Model outcomes are independent of olivine dose. Given the particle size distribution of our olivine product (Fig. S1), 0.12 µm corresponds to a mass fraction of 9.0% of olivine applied. At pH = 6.7, the model predicted a layer of thickness 0.04 µm, or 4.0 mass%. Disproportionality between skin thickness and mass weathered is due to complete consumption of very fine particles. At pH 5.4, the model predicted that particles smaller than 2 µm contributed 56% of total weathering, particles smaller than 20 µm contributed 90%. At pH 6.7, these fractions were 66% and 93%, respectively, according to the model. Above model estimates for overall weathering (all particle sizes) at both pH values are below our experimental value for OLIV1 (Tables 5, 6), and are around our value for OLIV2 (Table 5). We conclude that modelled and measured weathering rates differ by less than one order of magnitude for the lower doses. But also that negative feedback occurred at high olivine doses in soil. As stated, the model ignores this and thus overestimates weathering at high olivine doses, under our conditions. Olsen's model was based on a compilation of laboratory data for (sometimes pre-treated) olivine grains in stirred buffer solutions at various pH; the contrasting conditions complicate a direct comparison with our data. Feedbacks in soil are possibly due to changes in boundary layer pH, high aqueous Si concentration, or the formation of a passivating silica layer on olivine particles [10], [20].

The model indicates that, at our particle size distribution, particles smaller than 20 µm accounted for about 90% of all weathering during the trial. If this is true, then a drop in weathering rate must be expected over longer periods as fine particles are consumed. On the other hand, we do not know the impact on olivine particles of prolonged exposure to the soil environment. Also, the large discrepancy between measured BET surface of our olivine product (4.8 m² g⁻¹), and surface area calculated from measured particle size distribution (0.39 m² g⁻¹ for spherical particles) warrants caution in extrapolating the model predictions over longer periods. The high ratio between the two estimators of specific surface area may indicate micro-porosity or micro-cracks, phenomena recently documented for olivine [21].

Carbon sequestration potential

Gross sequestration amounts to 1.2 kg CO₂ per kg of olivine weathered (Eq. 1), and we use this figure to estimate gross CO₂ sequestration associated with weathering in our experiment. By our most conservative estimate of weathering (14.8%), gross sequestration at the dose of 8 g olivine per pot (1630 kg ha⁻¹, a dust layer of about 0.1 mm thickness) was equivalent to 29 10³ kg CO₂ km⁻² during the course of the experiment. (This is calculated as 0.148 * 1630 kg ha⁻¹ * 1.2 kg CO₂ per kg olivine = 290 kg ha⁻¹ = 29 10³ kg km⁻².) The amount of olivine weathered, and hence of gross CO₂ sequestration, roughly doubled for each fivefold increment of olivine applied. (Multiply fraction weathered from Table 5 with corresponding olivine dose, Table 1.) At the highest dose (OLIV4), gross sequestration as calculated by the same rule (1.2 kg kg⁻¹) was 269 10³ kg CO₂ km⁻². To put these figures into perspective: if 29 to 269 10³ kg CO₂ km⁻² would be sequestered annually on the world's entire agricultural area, this would correspond to 1.5 to 13.9 Pg CO₂ a⁻¹ or 4.7% to 43.6% of the annual global CO₂ release from combustion of fossil fuels (emission data 2008 in [22]). Our values can also be compared to annual CO₂ sequestration by natural silicate weathering as calculated for catchments. The global average is estimated at 1.9 10³ kg C (7.0 10³ kg CO₂) km⁻²a⁻¹, based on the GEM-CO₂ model [23]. Values below 5 10³ kg CO₂ km⁻²a⁻¹ were listed for acidic formations in the humid tropics [24], whereas the basaltic Deccan Traps would capture some 55 10³ kg CO₂ km⁻²a⁻¹ [25]. As an average for the volcanic Japanese Archipelago [26], reports an intermediate value of 6.05 10³ kg C (about 22 10³ kg CO₂) km⁻²a⁻¹. So, gross CO₂ sequestration at our lowest dose was four times larger than the annual global average, and about 30% higher than the annual Japanese average. Gross sequestration in our OLIV2 treatment (54.7 10³ kg km⁻² at olivine dose of 8150 kg ha⁻¹) would be similar to the high extreme given for the Deccan Traps. This is still far below extreme rates reported for carbonation of certain mine wastes [27].

Extrapolated to the world's agricultural area, our lower and upper olivine doses would correspond to global olivine inputs of 8 and 1000 Pg. For comparison: annual global hard coal production in 2011 was about 6.2 Pg [28]. Further extrapolation must take into account that feedbacks in oceanic carbonate chemistry would reduce CO₂ sequestration efficiency by some 20% according to [6]; but also that, once applied, olivine would continue to sequester CO₂ during many years, if our relative weathering rates (Table 5) remained valid over longer periods. Net sequestration will be smaller than gross sequestration, as part of the carbon captured as bicarbonate by ‘enhanced weathering’ may first be neutralized by soil acids, to escape as CO₂ back into the atmosphere. Of course, this is no zero-effect operation because the neutralisation of soil acidity would otherwise have required fossil carbonates, releasing CO₂ into the atmosphere. Indeed, basic silicates can be used as an alternative for the common practice of liming to counter soil acidification. Contrary to liming with chalk, there is no net release of CO₂ here. Further, estimating net sequestration potential requires accounting for CO₂ equivalents spent in mining, transporting and milling olivine rock. According to [29], the total energy required for crushing and milling, including tertiary (ultra-fine, 10 µm) milling, corresponds to an emission of 174 g CO₂ per kg CO₂ sequestered, far more important than the CO₂ cost for mining and transportation. This value was calculated by [29] from primary sources [30], [31].

Notes on olivine application to agricultural systems

It is difficult to anticipate all impacts that input of fresh basic silicates may have on agricultural systems. Impacts are likely to vary widely depending on soil, crop and climate characteristics. Based on Eq. 1 as well as our observations (Tables 3, 4) it can be expected that massive application of fast-weathering silicates will induce changes in soil pH and associated chemistry, and this may affect the availability of plant nutrients. Mg-induced Ca deficiency, for example, is well known from non-adapted crop species on serpentine-derived soils [32]. Such deficiencies could possibly be corrected in well-managed high input systems, but may pose problems in some extensively managed systems. Other aspects that warrant closer inspection are (a) the possible enhancement of soil organic matter decomposition which might result in loss of soil quality and net CO₂ emission (liming of acidified soil is known to reduce soil organic carbon stocks); and (b) the possible fixation of phosphates on freshly formed Fe-complexes, if silicates rich in Fe are used.

Finally, the relatively fast release of bioavailable Ni from olivine into the food chain and the wider environment could set limits to permissible olivine doses. While our experiment revealed no negative impacts on plant growth, it seems not currently possible to set general no-effect thresholds. Bioavailability of Ni, and toxicity to plants of Ni in soils amended with soluble metal salts were studied by [33], [34], [35]. Those studies underline that while soil pH is key to bioavailability, thresholds for toxicity depend strongly on other characteristics such as cation exchange capacity (CEC) and, within the acidic pH range, organic matter content. The study by [33] showed that the growth of oats (Avena sativa L.) seedlings was reduced by 50% at thresholds for total Ni in soil that varied, depending on soil properties, from 41 to 1321 mg kg⁻¹, with lowest values for sandy soils. Ranges by [34], who used barley (Hordeum vulgare L.) root growth and tomato (Solanum lycopersicum L.) shoot growth to express Ni-induced inhibition, are of the same order of magnitude. For comparison, our maximum total soil Ni content was 284 mg kg⁻¹ in OLIV4 (Table S6), but this included all Ni in unweathered olivine, too. (Weathering in OLIV4 was estimated at 1.1% only; Table 5.) Bioavailable Ni in this treatment was 1.38 mg kg⁻¹, and had no negative impact on growth (to the contrary). Although all Ni in olivine added to soils will ultimately be released from the silicate, this does not necessarily imply a build-up of bioavailable or phytotoxic Ni. Toxicity thresholds in ‘aged soils’ (18 months after amendment with soluble Ni salts) were shown to be up to 100 times larger than in freshly amended soils [35], suggesting immobilisation of previously bioavailable Ni as time proceeds.

Ni effects on plant growth were also studied in areas that are naturally rich in Ni derived from igneous bedrock. Trace metals in pastures over basalts were studied by [36] in the French Massif Central. They reported total Ni in soil between 168 and 214 mg kg⁻¹, and values of around 2.17 mg kg⁻¹ for Ni in aboveground plant biomass. The latter is close to 2.67 mg kg⁻¹ found at the highest dose (OLIV4) in our trial. These are well below the threshold of 10 mg kg⁻¹ for Ni toxicity in plants according to [37]. Higher levels (11.1–39.3 mg kg⁻¹ in forage) were reported by [32] for pastures on serpentine-derived soils in Galicia, Spain. These were associated with high Ni accumulation in kidney tissue of grazing cattle, while Ni in liver and muscle tissues remained undetectable. For the same area in northern Spain, [38] reported Ni contents of 12 to 34 mg kg⁻¹ in foliage of various crops, with highest values for sugar beet (Beta vulgaris L.), but they found no indications of Ni toxicity. Growth inhibition was demonstrated for ryegrass (Lolium perenne L.), at shoot Ni concentrations of about 100 mg kg⁻¹ [39]. In summary, while toxicity of Ni to plants and animals remains undisputed, evidence in the cited studies occurred at much higher concentrations than those found in our trial at high olivine doses. Multi-annual studies on olivine application to different types of crops and soils are needed for a more comprehensive assessment of the risks associated with Ni from olivine. The use of other fast weathering basic silicates, low in heavy metals, might waive this issue, but such minerals are less abundant than olivine.

Conclusions

The weathering rate of finely ground olivine in our soil was substantial. Between 13% and 58% of added olivine weathered during the 32 weeks of the trial, depending on the method of estimation. This range applies to our lowest olivine dose only, and is higher than estimates made by the Olsen model, which is based on laboratory conditions. At higher doses, however, weathering (relative to dose applied) decreased steeply in our experiment, and was less than predicted by the model. This suggests negative feedbacks for weathering in the soil environment, unaccounted for by the model.

Olivine increased soil water pH and alkalinity, plant uptake of Mg, Si and Ni, as well as bioavailability of Mg and Ni in soil. Olivine clearly suppressed Ca uptake which is attributed to competition by Mg. All these effects increased with higher olivine doses. There was a slight negative effect on P content in plant biomass at high olivine doses. At the highest dose (204 10³ kg ha⁻¹), olivine increased plant biomass (+15.6%) as well as K concentration in plant biomass (+16.5%). The cause of increased growth remains unclear, and may at our plant nutrient levels be unrelated to increased K uptake.

The appreciable weathering rate and lack of evidence for negative impacts on plant growth support the feasibility of the ‘enhanced weathering’ concept. Yet, massive application of olivine in agriculture may cause imbalances in plant nutrition, notably at low Ca availability, and will bring Ni into the food chain. In our case, olivine increased grass Ni concentration already at the dose of 1630 kg olivine ha⁻¹, from 0.531 to 0.696 mg kg⁻¹. At the extreme olivine dose of 204*10³ kg ha⁻¹, grass Ni concentration was 2.67 mg kg⁻¹. Although this is below a phytotoxic threshold of 10 mg kg⁻¹, it implies that the use of olivine in agricultural systems must remain within certain limits. Long term field studies are required to assess such limits under different climatic, soil and crop conditions.