# Spatial and Temporal Abundance of Interacting Populations of White Clover and Grass Species as Assessed by Image Analyses

Authors: Anne Kjersti Bakken; Helge Bonesmo; Bård Pedersen
Journal: Dataset Papers in Science
 (2015)

---

## Abstract

The dataset comprises detailed mappings of two communities of interacting populations of white clover (Trifolium repens L.) and grass species under differing experimental treatments over 4-5 years. Information from digital photographs acquired two times per season has been processed into gridded data and documents the temporal and spatial dynamics of the species that followed from a wide range of spatial configurations that arose during the study period. The data contribute a unique basis for validation and further development of previously published models for the dynamics and population oscillations in grass-white clover swards. They will be well suited for estimating parameters in spatially explicit versions of these models, like neighborhood based models that incorporate both the dispersal and the local nature of plant-plant interactions.


---

## Body

## 1. Introduction

White clover (Trifolium repens L.) is a clonal legume that is commonly grown in grasslands managed for grazing and forage production. The facilitation of long term, stable coexistence of clover and grass in such a system is a challenge [1–4] and the community dynamics are still poorly understood. The community has therefore turned out to be an important model system for studying and understanding ecological interactions between plants and the consequences of these interactions for coexistence of populations of different plant species.Various mechanisms have been proposed to cause population oscillations and formation of clover patches in pastures and swards. Amongst them are facilitation and exploitation [5], balancing local extinction and invasion of clover [6], and synchronous fragmentation of ramets into small and vulnerable individuals [7]. Grass and clover populations are suggested to remain stable at a nitrogen level that balances their competitive advantages; otherwise oscillations are expected [5, 8]. An increase in nitrogen availability in the sward increases the abundance of grass, and clover is excluded by grass. Time delayed responses to changes in the availability of nitrogen may increase the probability of oscillations [5]. In order to stabilize field scale oscillations of white clover population size, it is suggested that it must have a patchy distribution, be able to invade space where it can achieve high growth rate, and avoid space where it will be outcompeted [9]. In addition, environmental and demographic stochasticity [10] affect dynamics of real ecological systems of interacting species. Such stochastic noise may also induce dynamic phenomena like oscillations, local extinctions, and spatial pattern formation, features that characterize dynamics of grass-clover systems [11, 12].Thus, both the temporal and spatial variation in clover abundance may be relevant to deduce mechanisms that cause the dynamics. The above theories suggest it is especially relevant to study the spatiotemporal dynamics of clover-grass systems under contrasting environmental conditions that differ with respect to nitrogen availability and disturbance regimes. Simulation models have been developed to explore the spatiotemporal dynamics of white clover populations [5, 6, 8, 9, 13]. There are, however, few reports from field studies contributing data on spatial dynamics [6, 14–17] and to the best of our knowledge none that cover temporal and spatial abundance of both grass and clover in contrasting environments.We have constructed an experiment with swards dominated by white clover and smooth meadow-grass (Poa pratensis ssp.pratensis L.), with different management regimes with respect to cutting frequency and nitrogen supply. The management and treatments were chosen to develop environments that gave differing competitive advantages of grass and clover.Both the spatial distribution and abundance of white clover and grass were recorded during four or five years by digital photography and image analyses. The results were further processed to gridded data appropriate as inputs to spatial statistical analyses as suggested by Pedersen et al. [18] or simulation models for population dynamics in space and time.


## 2. Methodology

Site description and experimental layout are as follows. The experiment was initiated in 2001 and undertaken at two sites in pastures/leys dominated by white clover and smooth meadow-grass. The first one was located at Mære Agricultural School (63°56′N, 11°25′E, 40 m a.s.l.) in central Norway in a pasture previously grazed by dairy cattle. The species sown in 1993 besides smooth meadow-grass and white clover were timothy (Phleum pratense L.) and meadow fescue (Festuca pratensis Hudson). In 2001 the timothy and meadow fescue had nearly disappeared and there was a low appearance of weeds. The soil at the site is a sandy loam, and the content of organic matter in the topsoil is 5%.The second experimental site was located at Kvithamar Research Centre (63°29′N 10°52′E, 40 m a.s.l.) in a ley sown in spring 2000 with white clover (variety Norstar, Graminor AS, Stange, Norway) and smooth meadow-grass (variety Entopper, Innoseeds b.b., Vlijmen, The Netherlands). Altogether 18 main plots (4 m × 8 m each) were split into small plots (0.5 m × 0.5 m) and sown either with white clover at a rate of 0.8 g m−2 or smooth meadow-grass at a rate of 1.7 g m−2 or a mixture of both species (0.5 g m−2 + 1.5 g m−2). The three types of small plots were randomly distributed within main plots and the respective numbers of them in each main plot were 20 with clover, 20 with grass, and 88 with a mixture of the two species. The leys were cut once in late summer 2000 and were not fertilized until spring 2001. The soil at the site is a silt loam (20–25% clay and 55–65% silt) and the content of organic matter in the topsoil is 8%.From spring 2001 until autumn 2004 (Mære)/spring 2005 (Kvithamar), a factorial experiment with two harvesting regimes combined with three nitrogen-fertilizer levels was conducted at both sites. At Mære the treatments were replicated four times on 4 m × 4 m plots within a total area of 480 m2 and at Kvithamar there were three replicates laid out on 4 m × 8 m plots within an area of 684 m2. Plots were bordered with 20 cm wide strips kept free from vegetation by regular spraying with glyphosate. The plots were cut at a stubble height of 5 cm either two (late June and late August) or four times (late May, late June, late July, and late August) each growing season. The first/second and second/fourth cuts were always synchronised.Nitrogen was supplied from a compound mineral fertiliser (N-P-K, 18 : 3 : 15) at rates of 0, 10, or 30 g m−2 yr−1. In the first treatment, 1.5 g P and 5.1 g K m−2 were supplied from a P-K fertilizer.In the two-cut regime, 60% of the fertiliser was distributed in early spring and 40% after the first cut, whereas in the four-cut regime 40% was supplied in spring and 20% after each of the cuts in late May, late June, and late July.Dry yields were recorded at each of the cuts at Kvithamar. For the years 2002–2005, the content of clover in plotwise subsamples of the yield was determined by Near Infrared Reflectance Spectroscopy [19]. At Mære, there were no yield recordings. On both sites and just before the cuts in late August in 2003 and in late June and August in 2004, the height of the plant canopy was determined by a rising-plate meter [20]. There were 128 and 64 recording points within each plot at Kvithamar (Figure 1) and Mære (Figure 2), respectively. Precipitation, temperature, and global radiation for the experimental period are available at Agrometeorology Norway, lmt.bioforsk.no, stations Kvithamar and Mære.Figure 1
Orientation of the maps at Kvithamar: (a) 18 plots (4 m × 8 m plots each) and (b) 128 subplots (small plots: 0.5 m × 0.5 m plots each).
(a)
(b)Figure 2
Orientation of the maps at Mære: (a) 24 plots (4 m × 4 m plots each) and (b) 64 subplots (small plots: 0.5 m × 0.5 m plots each).
(a)
(b)Plots harvested two times per season and supplied with high levels of N were, by time, invaded by tussock forming and more erect and high yielding grass species than smooth meadow-grass. They wereElytrigia repens L.,P. pratense L.,Agrostis capillaris L., andLolium perenne L. After harvest, the initial regrowth rate expressed as plant coverage of soil surface, was slower in such plots. In line with this, bare soil at the time of image acquisition (cf. next section) does not indicate low grass abundance or production.Determination of plant coverage image acquisition and processing was as follows. To determine the coverage of clover, grass, and dicotyledonous weeds in the experimental fields, they were mapped by means of digital photography about ten days after the cuts in late June and August every year from 2001 and onwards (Figure3). To ensure that the images were taken at the same positions in the fields at each acquisition, 1.5 m steel tubes with a diameter of 1.9 cm were inserted into the ground at the corners of the 4.0 m × 4.0 m and 4.0 m × 8.0 m plots. Portable rails were anchored on these tubes on each side of the plot and on the rails, across the plot, a moveable bar was placed, and a sliding 0.5 m × 0.5 m metal frame that delimited the squares to be analysed sequentially was placed upon the bar (Figure 3).Figure 3
Positioned image acquisition to map the grass-clover sward.The area to be photographed was shielded from sunlight by an opaque parasol and the internal blitz of the camera was used to eliminate shadows and ensure even illumination. The camera was mounted to a rack on the metal frame such that the distance from the camera to the ground was kept constant at 80 cm. A digital Olympus Camedia C-3040Zoom camera was used, and the images were stored in a compressed format that gave a spatial resolution of approximately 1600 × 1600 pixels m−2.The information from the digital colour photographs was processed by software (“Trifolium.exe”) specially designed for the purpose. The software classifies each pixel in the image as clover, grass, weeds, and bare ground. A thorough description of the image analysis techniques is given by Bonesmo et al. [21]. The 0.5 m × 0.5 m single images were after processing by the software (Figure 4) combined into one picture (TIF-format) (Figure 5) or data file for each plot after identifying overlap zones and removing them from one of the overlapping images. The procedure for removing overlap zones as well as the software Trifolium.exe is available from the authors. The copyright for the source code of Trifolium.exe is not owned by the authors and cannot be published here.Figure 4
Outputs at sequential steps in the analysis of (a) a digital colour image by the software Trifolium.exe, (b) grey level image, (c) edge image, and (d) classified image, where red indicates clover, blue indicates grass, and black indicates soil (published with permission from Taylor & Francis Group, and previously printed in [21]).
(a)
Image
(b)
Grey level image
(c)
Edge image
(d)
ClassificationFigure 5
An example of a 4 m × 4 m plot with patterns occurring from identification of pixels either as clover (red), grass (blue), bare ground (black), or dicotyledonous weeds (white).


## 3. Dataset Description

The dataset associated with this Dataset Paper consists of 21 items which are described as follows.1
Photos for each plot and image acquisition at Kvithamar. All photos together covered the whole plot. The photos are named according to the number of image acquisition (9 image acquisitions), plot number (18 main plots), and subplot number (128 subplots). These photos were the basis for the image analyses conducted by the software Trifolium.exe. There are altogether 20736 photos (9 image acquisitions × 18 plots/acquisition × 128 photos/plot). A few photographs are missing.
Olympus Camedia C-3040Zoom camera
Olympus Camedia master2
Photos for each plot and image acquisition at Mære. All photos together covered the whole plot. The photos are named according to the number of image acquisition (8 image acquisitions), plot number (24 main plots), and subplot number (64 subplots). These photos were the basis for the image analyses conducted by the software Trifolium.exe. There are altogether 12288 photos (8 image acquisitions × 24 plots/acquisition × 64 photos/plot). A few photographs are missing.
Olympus Camedia C-3040Zoom camera
Olympus Camedia master3
Treatments allocated to the different plots are listed together with harvest dates, data for harvested yield, and clover proportion in dried yield samples, as determined by NIRS. Yields were recorded at the site Kvithamar only. The column Site identifies whether the site is Kvithamar or Mære. The column Plot Number presents the plot number within site (1–18 at Kvithamar and 1–24 at Mære). In the column Harvesting Regime, 1 indicates two cuts per season and 2 indicates four cuts per season. In the column Nitrogen Supply Rate, 1 indicates 0 g, 2 indicates 10 g, and 3 indicates 30 g. The column Yield One presents the dry matter yield (g m−2) at harvest in spring (May/June); Yield Two, the dry matter yield (g m−2) at harvest in early summer (June/July); Yield Three, the dry matter yield (g m−2) at harvest in late summer (July/August); Yield Four, the dry matter yield (g m−2) at harvest in autumn (August/September). The column Clover One presents the proportion (%, weight/weight) of clover in yields in spring; Clover Two, the proportion (%, weight/weight) of clover in yields in early summer; Clover Three, the proportion (%, weight/weight) of clover in yields in late summer; Clover Four, the proportion (%, weight/weight) of clover in yields in autumn. The column Date One presents the date of harvest in spring; Date Two, the date of harvest in early summer; Date Three, the date of harvest in late summer; Date Four, the date of harvest in autumn. The harvest dates at Mære were 2-3 days later than respective dates at Kvithamar. A dot (·) in a cell indicates that no value was recorded for this combination of harvesting regime, N-supply rate, year, and time of season.
115
17
0
Site
Plot Number
Harvesting Regime
Nitrogen Supply Rate
Year
Yield One
Yield Two
Yield Three
Yield Four
Clover One
Clover Two
Clover Three
Clover Four
Date One
Date Two
Date Three
Date Four4
Registration of the first image acquisition at Kvithamar containing 18 matrices that contain 2400 × 1200 cells describing the 8 m × 4 m plots. Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.5
Registration of the second image acquisition at Kvithamar containing 18 matrices that contain 2400 × 1200 cells describing the 8 m × 4 m plots. Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.6
Registration of the third image acquisition at Kvithamar containing 18 matrices that contain 2400 × 1200 cells describing the 8 m × 4 m plots. Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.7
Registration of the fourth image acquisition at Kvithamar containing 18 matrices plots that contain 2400 × 1200 cells describing the 8 m × 4 m. Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.8
Registration of the fifth image acquisition at Kvithamar containing 18 matrices that contain 2400 × 1200 cells describing the 8 m × 4 m plots. Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.9
Registration of the sixth image acquisition at Kvithamar containing 18 matrices that contain 2400 × 1200 cells describing the 8 m × 4 m plots. Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.10
Registration of the seventh image acquisition at Kvithamar containing 18 matrices that contain 2400 × 1200 cells describing the 8 m × 4 m plots (plot 9 is missing). Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.11
Registration of the eighth image acquisition at Kvithamar containing 18 matrices that contain 2400 × 1200 cells describing the 8 m × 4 m plots (plot 8 is missing). Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.12
Registration of the ninth image acquisition at Kvithamar containing 18 matrices that contain 2400 × 1200 cells describing the 8 m × 4 m plots. Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.13
Registration of the first image acquisition at Mære containing 24 matrices that contain1200×1200 cells describing the 4 m × 4 m plots. Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.14
Registration of the second image acquisition at Mære containing 24 matrices that contain1200×1200 cells describing the 4 m × 4 m plots (plots 11, 19, and 20 are missing). Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.15
Registration of the third image acquisition at Mære containing 24 matrices that contain1200×1200 cells describing the 4 m × 4 m plots. Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.16
Registration of the fourth image acquisition at Mære containing 24 matrices that contain1200×1200 cells describing the 4 m × 4 m plots. Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.17
Registration of the fifth image acquisition at Mære containing 24 matrices that contain1200×1200 cells describing the 4 m × 4 m plots. Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.18
Registration of the sixth image acquisition at Mære containing 24 matrices that contain1200×1200 cells describing the 4 m × 4 m plots. Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.19
Registration of the seventh image acquisition at Mære containing 24 matrices that contain1200×1200 cells describing the 4 m × 4 m plots. Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.20
Registration of the eighth image acquisition at Mære containing 24 matrices that contain1200×1200 cells describing the 4 m × 4 m plots. Each cell in the matrices contains one digit (1, 2, 3, or 4), which identifies the character of single cells as being covered by clover (1), grass (2), bare ground (3), or dicotyledonous weeds (4). For a few plots, there are smaller areas that are not identified. Such cells are annotated with the digit 0.21
Recordings of plant canopy height at Kvithamar and Mære in late August in 2003, late June in 2004, and late August in 2004. There were 128 and 64 recording points within each plot at Kvithamar and Mære, respectively.
11521
5
0
Cut Date
Site
Plot Number
Subplot Number
Height


## 4. Concluding Remarks

The dataset presented here will contribute a unique basis for validation and further development of previously published models for the dynamics and population oscillations in grass-white clover swards. The set documents the spatial dynamics of the two species that followed from a wide range of spatial configurations that arose during the experiment. It should therefore be well suited for estimating parameters in spatially explicit versions of these models, like neighborhood based models that incorporate both the dispersal and the local nature of plant-plant interactions (e.g., [22–25]). For example, the hypothesis that a patchy distribution of clover and its ability to invade space where it is not outcompeted will damp oscillations may be tested.



---

*Source: 620164-2015-06-08.xml*