<<

Clonal architecture and patch formation of perfoliatus L. in response to environmental conditions

Susanne R. Wolfer

Promotor: Prof.dr. M. Scheffer Hoogleraar Aquatische ecologie en waterkwaliteitsbeheer Departement Omgevingswetenschappen Wageningen Universiteit

Co-promotoren: Dr.ir. E. H. van Nes Universitair docent, leerstoelgroep Aquatische ecologie en waterkwaliteitsbeheer Departement Omgevingswetenschappen Wageningen Universiteit

PD Dr. D. Straile Limnological Institute University of Constance

Promotiecommissie: Prof.dr. F. Berendse, Wageningen Universiteit Prof.dr. H.J. Lindeboom, Wageningen Universiteit Prof.dr. J.T.A. Verhoeven, Universiteit Utrecht Dr.ir. G.J. van Geest, Deltares, Utrecht Dr. M.S. van den Berg, Rijkswaterstaat, Waterdienst (RIZA), Lelystad

Clonal architecture and patch formation of L. in response to environmental conditions

Susanne R. Wolfer

Proefschrift ter verkrijging van de graad van doctor op gezag van de rector magnificus van Wageningen Universiteit, Prof.dr. M.J. Kropff, in het openbaar te verdedigen op maandag 20 oktober 2008 des namiddags te vier uur in de Aula.

Clonal architecture and patch formation of Potamogeton perfoliatus L. in response to environmental conditions / S.R. Wolfer Thesis Wageningen University - With ref. - With summary in Dutch ISBN 978-90-8585-211-7

Contents

List of figures...... 6

List of tables ...... 7

Chapter 1 Introduction...... 9

Chapter 2 Lake and macrophyte characteristics...... 15

Chapter 3 Spatio-temporal dynamics and plasticity of clonal architecture in Potamogeton perfoliatus... 23

Chapter 4 Density control in Potamogeton perfoliatus L and Potamogeton pectinatus L...... 35

Chapter 5 Effect of sediment structure and nutrients on the growth of Potamogeton perfoliatus ...... 45

Chapter 6 To share or not to share - clonal integration in ramets of Potamogeton perfoliatus...... 55

Chapter 7 Modelling the clonal growth of the rhizomatous macrophyte Potamogeton perfoliatus ...... 65

Chapter 8 Synthesis and discussion...... 85

References ...... 95

Summary ...... 111

Samenvatting ...... 113

Acknowledgement ...... 115

List of publications ...... 117

Curriculum vitae ...... 119

List of figures

Fig. 1.1 Research at different scales ...... 10 Fig. 2.1 Location of Lake Constance in Europe and geography of Lake Constance ...... 15 Fig. 2.2 Lake Constance, satellite image...... 16 Fig. 2.3 Daily water levels at Konstanz (Upper Lake) for the years 1998 and 1999 ...... 17 Fig. 2.4 P. perfoliatus (a) with 17 ramets and 8 branchings (b) turions...... 18 Fig. 2.5 Annual growth cycle of P. perfoliatus in Lake Constance ...... 18 Fig. 2.6 P. perfoliatus patches, Lower Lake Constance...... 18 Fig. 2.7 (a) 3-D representation of distribution and shoot length (b) density (c) changes in patch growth ...... 20 Fig. 2.8 Patch patterns and patch dynamics over a period of 3 years, north-west Reichenau...... 21 Fig. 3.1 Apical section of P. perfoliatus plant...... 25 Fig. 3.2 Reconstructed spatio-temporal growth dynamics of P. perfoliatus at poor growth and rich growth site. ...27 Fig. 3.3 (a) Temporal development of ramet numbers and (b) average shoot length of P. perfoliatus...... 28 Fig. 3.4 (a) Frequency distribution of spacer length (b) spacer length between primary and higher order shoots..28 Fig. 3.5 (a) Shoot growth rate and (b) growth rate per plant of P. perfoliatus...... 29 Fig. 4.1 (a) Average number of ramets produced per site (b) average number of ramets recruited per initial shoot (c) average total dry weight production per initial shoot of P. perfoliatus...... 37 Fig. 4.2 Average shoot, rhizome and root biomass allocation (%) of P. perfoliatus ...... 38 Fig. 4.3 (a) Average rhizome spacer length (b) relationship between total rhizome length and final number of ramets per site (c) average total rhizome length per plant of P. perfoliatus...... 38 Fig. 4.4 Average number of rhizome branchings (a) per ramet (b) per plant of P. perfoliatus ...... 39 Fig. 4.5 (a) Average number of tubers (b) fresh weight frequency distribution of tubers of P. pectinatus...... 40 Fig. 5.1 Mesocosm, Limological Institute, University of Constance ...... 48 Fig. 5.2 Ramet number, shoot length, and total biomass of P. perfoliatus in unfertilized and fertilized quadrats ...49 Fig. 5.3 Biomass allocation and spacer length of P. perfoliatus without and with fertilization...... 50 Fig. 5.4 Average rhizome spacer length of P. perfoliatus in relation to position within rhizome...... 50 Fig. 5.5 Ramet number, ramet length and patch size of P. perfoliatus on unfertilized and fertilized sediment ...... 51 Fig. 6.1 (a) Aquarium (b) sketch of an experimental unit ...... 57 Fig. 6.2 Total shoot biomass of both offspring ramets of P. perfoliatus in the 2-ramet versus 3-ramet system ...... 58 Fig. 6.3 Shoot biomass of the offspring ramets O1 and O2 in the different treatments ...... 58 Fig. 6.4 Biomass ratio of second offspring and first offspring of P. perfoliatus grown in aquaria ...... 59 Fig. 6.5 (a) Shoot (b) rhizome and (c) root allocation of 3-ramet system and 2-ramet system ...... 60 Fig. 6.6 Spacer length between offspring ramets in relation to first offspring shoot biomass...... 60 Fig. 7.1 Schematic overview of the factors that determine the growth of each shoot ...... 69 Fig. 7.2 Scheme of the clonal growth of P. perfoliatus...... 71 Fig. 7.3 Schematic representation of a clone and the transport of energy from the ramets to the turions ...... 71 Fig. 7.4 Example of a generated pattern (clone growth during one year) ...... 73 Fig. 7.5 20 years run, (a) year 2009 (b) 1990 - 2010 ...... 74 Fig. 7.6 Relation between the mean number of ramets and biomass at day 249 ...... 75 Fig. 7.7 Modelling results on the influence of branching probability after a simulation period of 20 years...... 76 Fig. 7.8 Patch density and expansion in dependence to rhizome growth direction relative to the previous year....77 Fig. 7.9 Simulation of the response of P. perf. biomass to increasing and subsequently decreasing turbidity ...... 78 Fig. 7.10 Cluster analysis of the sensitivity parameters of the model ...... 79 Fig. 7.11 Time course of the relative sensitivity coefficients of various parameters on model outcomes ...... 80 Fig. 8.1 Factors that influence patch characteristics...... 85 Fig. 8.2 Simplified responses of patch growth to light and nutrients as found in this study ...... 88

6 List of tables

Tab. 1.1 Environmental parameters which influence plant growth and / or clonal architecture of submerged macrophytes...... 11 Tab. 2.1 Lake Constance in brief...... 16 Tab. 3.1 Biomass allocation (% dry weight) of P. perfoliatus in poor growth and fertilized poor growth quadrats ..30 Tab. 3.2 N and P contents (% dry weight), C:N, C:P and N:P weight ratios of P. perfoliatus leaf tissue in poor growth, fertilized poor growth, and rich growth quadrats...... 30 Tab. 4.1 Calculated patch sizes and spatial densities of P. perfoliatus ...... 39 Tab. 5.1 Sediment characteristics and nutrient concentrations in porewater of sites with or without vegetation cover ...... 49 Tab. 5.2 Tissue nutrient contents at non-fertilized and fertilized sites ...... 49 Tab. 6.1 ANOVA with shoot biomass as dependent variable and severing, shaded ramet order (first or second offspring) and aquarium as independent variables...... 59 Tab. 7.1 An overview of the main differences between non-clonal macrophyte models and the clonal growth model CLOMO ...... 67 Tab. 7.2 Default parameters for the clonal model of P. perfoliatus...... 68 Tab. 7.3 Model results of growth parameters in comparison to in-situ and mesocosm assessments...... 73 Tab. 8.1 Simplified responses of clonal architecture on the growth factors nutrients, light and clonal integration..86 Tab. 8.2 Simplified relationships between clonal architecture/growth factors and patch size...... 88

7 8 Introduction

Chapter 1. INTRODUCTION

Submersed macrophytes usually occur in patches which differ in size, height, shape and shoot density. Also Potamogeton perfoliatus L. (perfoliate pondweed), a species common to many Central European Lakes, forms such patches ranging from semi-continuous to widely dispersed stands over a wide range of scales (Lehmann 1997). In many lakes, Potamogeton patches change location, size, shape and density in the long run (Schmieder 1997), whereas in other lakes they remain spatially stable in the medium range (Gafny & Gasith 1991, Walser 1996) or disappear and reappear from year to year (Scheffer et al. 1992). Macrophyte patches may be round (Valta-Hulkkonen et al. 2003), elliptical (Gusak 1983), ring-shaped (Wikberg & Svensson 2003) or irregular. They may consist of a few shoots or extend to several kilometres and comprise millions of shoots.

Since P. perfoliatus is a modular plant consisting of reiterations of ramets along that are anchored in the sediment, the visible spatial patterns provide only a partial picture of the complicated belowground patterns. The proper spatial pattern is hidden to the researcher: neighbouring shoots may stem from two different , while one plant might spread over several square meters.

The knowledge about the formation of spatial patterns is important for a variety of reasons and applications. Presence or absence of macrophytes as well as size, shape and density of vegetation patches are important for both biotic and abiotic processes in the ecosystem. Spatial patterns play a vital role in structuring of the littoral zone of lakes (Jeppesen et al. 1997). Macrophytes provide a refuge for herbivores (Timms & Moss 1984), shelter for prey (Dvorak & Best 1982), refuge for young fish (Walser & Löffler 1993), substrate for epiphytes (Dvorak & Best 1982, Wetzel 1983), and food for waterfowl, mammals and invertebrates (Moen & Cohen 1989). Moreover, macrophyte composition and patch characteristics affect the number and diversity of invertebrates (De Vos 1954, Webster et al. 1998) influence the distribution and age-characteristics of littoral fish (Weaver et al. 1997), and modify littoral food webs (Beklioglu & Moss 1996, Lauridsen et al. 1996), as well as host-epiphyte relationships (Levin & Mathieson 1991).

In dependence of their spatial characteristics, patches influence nutrient cycles (Barko & James 1997), reduce current velocity (Madsen & Warncke 1983, Fonseca 1989, Losee & Wetzel 1993) and consequently enhance sediment stability, increase sedimentation of particles (Sand-Jensen 1998), stabilize clear water by preventing resuspension of sediments (Meijer et al. 1999) and reducing nutrient levels (Van Donk et al. 1989, Sand-Jensen 1998), lower sediment pH, raise sediment redox potential, release oxygen through their roots, increase P retention in the sediment (Jaynes & Carpenter 1986), and release organic carbon which might influence the bacterial activity or the microbial biomass in the sediment as well as the denitrification potential within deeper root zones (Karjalainen et al. 2001). The P released by decomposing aquatic plants is an important source of P for pelagic phytoplankton (Barko & Smart 1980, Smith & Adams 1986). Plant patches thus substantially modify their own environment and that of the littoral biocoenosis.

Despite the recognition of the importance of spatial patterns of submersed macrophytes for the littoral ecosystem, the rules of pattern formation are not well understood (Lehmann 1997). Why does a fast growing species like P. perfoliatus not cover larger uniform lake areas after a certain number of years, even if competitors are absent? Since seeds are considered as relatively unimportant for propagation (Chambers & Kalff 1985), the knowledge of its architectural programme at the individual plant level is the key to understanding patch patterns and patch dynamics. As in other clonal plants, the growth patterns of Potamogeton are assumed to be non-random (Dale 1999). Clonal plant architecture comprises growth rules such as number of ramets, branching patterns (frequency and

9 Chapter 1

angle) of the rhizome, biomass allocation to rhizomes and shoots, and the rhizome spacer length in between consecutive shoots. These rules are thought to be plastic to a certain degree allowing for adaptive responses to a spatially heterogeneous environment (Marbà & Duarte 1998).

Whereas clonal architecture and spatial patterns of terrestrial plant species have received considerable attention (Huber & Stuefer 1997, Klimes et al. 1997), there are less studies on clonal growth of aquatic macrophytes. The aqueous medium requires a completely different physiology of aquatic plants compared to terrestrial forms. Moreover, life cycle, reproduction, morphology, clonal architecture and growth pattern are species-specific and can not simply be transferred to other species. If we want to model clonal aquatic macrophytes, we therefore need more consistent studies. Most existing spatial studies on aquatic macrophytes concern above-sediment patch patterns (Vidondo et al. 1997, Bradley & Stolt 2006). In-situ research on the architecture of rhizomes and roots is rare for difficulties with assessment under water and below the sediment (Ervin & Wetzel 1997). Previous work was mainly published on descriptive architecture of seagrasses (Olesen & Sand- Jensen 1994b, Olesen & Sand-Jensen 1994a, Duarte et al. 1996, Duarte et al. 1998, Marbà & Duarte 1998, Klimes 2000, Molenaar et al. 2000, Nielsen & Pedersen 2000, Pagès et al. 2000, Paling & McComb 2000), and Reusch et al. 1998 excavated some seagrass systems in the context of a genetic study. Freshwater macrophyte clonal morphology included studies on emergent Scirpus (Clevering & Hundscheid 1998), Phragmites australis (Vretare et al. 2001), Potamogeton pectinatus (Idestam- Almquist & Kautsky 1995), and (Brux et al. 1987). In general, architectural studies are neither spatially concise nor related to patch characteristics, and spatial studies exclude clonal architecture. In both fields, comprehensive studies on P. perfoliatus are lacking.

The aim of this thesis is to describe the spatial patterns of P. perfoliatus, to search for the mechanisms that generated the patterns, and to understand the influence of plants on the patterns. Our hypothesis was that in case of the submersed macrophyte species P. perfoliatus, (1) architectural growth rules pre-determine spatial patterns, and (2) environmental factors influence architecture, i.e. the species shows morphological plasticity. In order to understand the spatial growth patterns, a combination of in-situ observations, aquarium experiments and mesocosm experiments was used at small, medium and landscape-size scales (Fig. 1.1).

m cm a water 25 100 20 51 b90 30 50 c 38 18 20 80 88 land 46 66 56 92 80 70 47 17 57 44 70 78 62 15 60 63 87 29 75 50 36 37 42 16 53 72 15 74 10 40 28 71 3 26 52 30 8 33 14 24 32 7 84 85 31 40 5 20 12 1 6 83 13 22 5 41 100 m 10 11 21 39 2 67 0 0 10 50 m 0 5 10 15 20 25 m 0 102030405060708090100 species 1 2 3 4 cm

Fig. 1.1 Research at different scales: (a) landscape scale: aerial image (DLR 2001), Reichenau-Niederzell, Lower Lake (b) medium scale: plan-view of lake vegetation at 2.50 m depth, 1 m² resolution (c) small scale: plan-view, 1 m², map of underground structures; lines represent rhizomes, numbers represent shoots (method see Chapter 3)

10 Introduction

The core task of the research project was to assess growth patterns. As a result. the thesis provides a detailed architectural library of clonal architecture including biomass allocation, shoot length, and variation of spacer length, branching frequency and branching angle of P. perfoliatus in Lake Constance along with detailed maps of below-sediment rhizome networks (Chapter 3).

A number of environmental factors is thought to be responsible for the variability and plastic growth of aquatic macrophytes (Tab. 1.1). In our study, we investigated the effects of nutrient and light conditions on both architecture and pattern. In order to test for nutrient effects, we fertilized small quadrats in-situ, analysed clonal architecture and plant tissue samples at sites with different vegetation growth (Chapter 3), and carried out a mesocosm experiment with different nutrient levels (Chapter 5). In order to test for light effects, we carried out a self-shading experiment with three different stand densities and compared clonal architecture and patch characteristics (Chapter 4). An aquarium experiment focusing on the question of resource transfer between ramets also included a light component (Chapter 6).

Additionally, we were interested in clonal integration, i.e. the support of younger ramets by older ramets (or vice versa) by the transfer of nutrients, signals or photosynthetates. The question was whether clonal integration enhances survival, overall performance, and patch expansion in P. perfoliatus, and in which manner it affects clonal architecture (Chapter 6).

Tab. 1.1 Environmental parameters which influence plant growth and / or clonal architecture of submerged macrophytes factors related to light Canfield et al. 1985, Chambers & Kalff 1985, Dale 1986, Duarte & Kalff 1990b, attenuation sensu latu Abernethy et al. 1996, Agami et al. 2001 water depth Duarte & Kalff 1990b, Rea et al. 1998 slope Duarte & Kalff 1986, Collins et al. 1987, Duarte & Kalff 1990b, Rea et al. 1998, Mallik & Rasid 1993 hydrodynamics, waves, Keddy 1982, Keddy 1983, Chambers 1987, Kautsky 1987, Duarte & Kalff currents, fetch 1990b, Chambers et al. 1991, Idestam-Almquist & Kautsky 1995, Strand & Weisner 1996, Lehmann et al. 1997, Fonseca & Bell 1998, Rea et al. 1998, Barrat-Segretain 2001 water temperature Barko & Smart 1981a, Dale 1986, Chambers 1987, Pip 1989, Spencer & Ksander 1992 water chemistry Spence 1967, Casper & Krausch 1980, Lachavanne 1985, Canfield et al. 1985, Jackson & Charles 1988, Melzer et al. 1988, Husák et al. 1989, Lachavanne 1991, Royle & King 1991, Walser 1996, Demars & Harper 1998 sediment Aiken & Picard 1980, Mantai & Newton 1982, Barko & Smart 1983, 1986b, Chambers & Kalff 1985, Anderson & Kalff 1986, 1988, Chambers 1987, Barko et al. 1988, 1991, Duarte 1995, Walser 1996, Lehmann et al. 1997, Koch 1999, Xiao et al. 2007b competition Agami & Waisel 1985, Chambers & Prepas 1990, Abernethy et al. 1996 herbivory Armellina 1999, Jonzén et al. 2002

11 Chapter 1

All findings were used to create a clonal model based on the previous non-clonal growth model Charisma (Van Nes 2002, Van Nes et al. 2003), a ramet-based spatially explicit model (Chapter 7). Our model is the first detailed model including architectural, spatial, environmental, physiological and demographic components and simulates both individual clonal plant growth and patch growth. By simulating different clonal growth architecture, we built the bridge between clonal architecture and patch patterns.

In the synthesis and discussion (Chapter 8), the most important results are presented, conclusions are drawn, and different theories of architectural responses to habitat quality are discussed. Ultimately, we hope that the information on clonal growth rules can be used for the understanding and management of recolonization and restoration of littoral biotopes.

12 Introduction

13 Chapter 1

14 Lake and macrophyte characteristics

Chapter 2. LAKE AND MACROPHYTE CHARACTERISTICS

Lake characteristics

Lake Constance is a large water body in Central Europe shared by Germany, Switzerland and Austria (Figs. 2.1, 2.2). The perialpine lake has been formed by the action of glaciers about 15000 years ago. Some of its characteristics are listed below (Tab. 2.1). The lake consists of two separate basins which are connected by a 4 km stretch of the river Rhine. The larger basin, called ‘Obersee’ (Upper Lake), is deep and monomictic. It has steep banks, a relatively well-mixed epilimnion, a true pelagic zone, and a relatively small littoral zone of 15 % (IGKB 2004). The littoral banks of Lake Überlingen are, for example, about 20-30 m wide. The smaller basin, called ‘Untersee’ (Lower Lake), consists of three individual dimictic basins and has an extensive littoral area (28 %).

Fig. 2.1 Location of Lake Constance in Europe and geography of Lake Constance (adapted from Straile & Geller 1998, Pinnel 2007); 1: sampling area near the island of Reichenau, Lower Lake Constance (Chapter 3), 2: experimental area at Upper Lake Constance and location of mesocosm (Chapters 4, 5)

Lake Constance is situated in the warm-temperate climate zone with atlantic and continental influences and an average annual air temperature of 8-9 °C. The lake surface has not been frozen since the year 1963.

15 Chapter 2

Fig. 2.2 Lake Constance, satellite image (Landesvermessungsamt Baden-Württemberg)

Tab. 2.1 Lake Constance in brief (IGKB 2004)

Lake Constance Upper Lake Lower Lake location 47 n. lat., 9° e. long. elevation 395 m a.s.l. average air temperature 8 - 9 °C average air temperature, July 18°C average air temperature, January -1°C average water temperature 5.7 °C average precipitation 805 - 1380 mm/a main wind direction W, S/W catchment area 11500 km² volume 48.4 km³ 47.6 km³ 0.8 km³ surface area 535 km² 472 km² 62 km² shore length 273 km 186 km 87 km maximum length 63 km maximum width 14 km maximum depth 253 m 40 m mean depth 101 m 13 m

The hydrograph of Lake Constance is characterized by low water in January, high water in June after snow-melt in the Alps and a mean annual water level fluctuation of about 2 m (Fig. 2.3). These fluctuations are very important for biological and chemical processes in the littoral zone. During the last years, the annual fluctuations tended to increase and the absolute water levels tended to fall (IGKB 2004).

16 Lake and macrophyte characteristics

Fig. 2.3 Daily water levels at Konstanz (Upper Lake Constance) for the years 1998 and 1999, mean daily values and standard deviation (grey band) for the period 1817–2002 (Joehnk et al. 2004)

After an eutrophication period in the 1970ies and a slow recovery, the Upper Lake is now again meso-oligotrophic, the Lower Lake mesotrophic. In comparison to the Upper Lake, the Lower Lake has warmer water, higher water and sediment nutrient levels, and higher sedimentation rates (IGKB 2004). Most in-situ experiments in the framework of this study were carried out at the Lower Lake. The sediment for the mesocosm experiments was also taken from the Lower Lake, while lake water for the aquaria and mesocosm experiments originated from the Upper Lake.

Macrophyte characteristics

The lake has extensive stands of submersed macrophytes which underwent a significant change in abundance, species composition and zonation during the last 50 years. The eutrophication period of the 1970s coincided with an increase of eutraphentic species and a decrease of most other species. During the re-oligotrophication period of the 1980s, most disappeared species returned back to the lake. The present state is similar to the one before eutrophication. For a detailed analysis of species, see Schmieder (1998).

Potamogeton perfoliatus - the species investigated in this study - is the second most common submersed species in Lake Constance after Potamogeton pectinatus (Schmieder 1997). Plants consist of primary ramets (i.e. leafy shoot with root), rhizomes and a varying number of secondary respectively higher order ramets (Fig. 2.4). P. perfoliatus has a vegetation period from May until September, vanishes in autumn due to senescence and autumn storms and survives the winter in the form of turions in the sediment (Fig. 2.5). P. perfoliatus occupies depths between 2 and 5 meters, develops distinct patches, has shoots up to 4.5 m in length (Fig. 2.6) and develops flowers that may or may not reach the water surface depending on the water level fluctuations. In some years, canopies are formed on the surface. During the phase of eutrophication, P. perfoliatus decreased considerably in the Lower Lake, but increased again after re-oligotrophication.

17 Chapter 2

a

b

Fig. 2.4 P. perfoliatus (a) plant with 17 ramets and 8 branchings (b) turions

Fig. 2.5 Annual growth cycle of P. perfoliatus in Lake Constance; left axis/squares: no of shoots/m² (rich growth patch, Chapter 3); right axis/curve: av. daily water temperature at 50 cm depth (°C) at Pegel Bregenz (www.vorarlberg.at)

a b

Fig. 2.6 P. perfoliatus patches, Lower Lake Constance (photos: J. van Schie, Rijkswaterstaat Waterdienst)

18 Lake and macrophyte characteristics

Typical medium to landscape scale macrophyte composition at the Lower Lake

The main project area was located at the Lower Lake, in a bay W-SW of Niederzell on the island of Reichenau (Fig. 2.1). This bay has an extensive shallow littoral area with relatively strong wind exposure. The main species are Potamogeton perfoliatus, Potamogeton pectinatus, Chara spp. and Najas marina ssp. intermedia. Medium to landscape scale data was collected with permanent grids and airborne multispectral scanners.

In order to gain a first overview, a typical landscape-scale grid (25 x 25 m, squares of 1 m²) was established by fixed lines (average depth: 2.50 m) and plants were mapped by snorkelling along the lines. The area was mapped on 21.07.2000 and on 17.05.2001 respectively 25.07.2001 to analyse patch patterns and to look at seasonal changes (Fig. 2.7).

The remotely sensed data of the sampling area was acquired for the years 2000-2002 by a multispectral sensor onboard an aircraft. The Modular Inversion & Processing System (MIP) was used as a processing tool to account for aerosols, sun glitter correction, atmosphere, water surface and water column corrections with a physically based process chain (Heege et al. 2003). The final results are shown as percentage coverage of bare sediment, tall macrophytes (mainly P. perf. and P. pect.) and low macrophytes (Characeae) up to a depth of 4.5 m by fitting with ground truth (Fig. 2.8). Whereas this earlier automatic classification did not discriminate between P. perfoliatus and P. pectinatus, the improved method now includes species discrimination. This is achieved by implementing data from a hyperspectral scanner and utilizing the optimal wavelengths from a comprehensive spectral library of eight macrophyte species, acquired by a submersible spectroradiometer (Pinnel 2007).

Both methods show that the plants form species-specific patches in this area (Figs. 2.7, 2.8). Chara spec. covers extensive bottom areas, especially in form of near-shore bands. In the examined area, P. perfoliatus has its main distribution at the shoreward side of the shallow zone and forms distinct patches. P. perfoliatus shoots reach 2.5 m (the water surface), density is frequently highest in the centre of the patch (Fig. 2.7b).

The comparison of processed scanner images from three different years shows that medium and high growth forms occur roughly at the same sites each year, but patch size and shoot lengths differ (Fig. 2.8). Short species (Chara ssp.) show more variability in bottom coverage. All differences are partly due to different development states at the assessment day due to different weather conditions during earlier growth phases. The P. perfoliatus patch mapped by snorkelling showed not much spatial expansion during one season (Fig. 2.7c). Although at some places, the plants extended by 3 m during the vegetation period, the net area of P. perfoliatus decreased by 1.9 %. Patch receding occurred mainly at the lakeward side with deeper water and may be partially due to light attenuation after increasing water levels.

19 Chapter 2

a

Potamogeton perfoliatus Potamogeton pectinatus Chara Najas marina

]

m

[ 2.5

h

t

g 2.0

n

e

l

1.5

t

o

o 1.0

h

S 0.5 25 0.0 25 20 20 15 15 10 m m 10 E 5 5 R O 0 SH

b m 25

20

15

10

5

0 0 5 10 15 20 25 m

density low medium high

c m 2 5

2 0

1 5

1 0

5

0 0 5 10 15 20 25 m both dates first visit second visit

Fig. 2.7 (a) 3-D representation of distribution and shoot length and plan view of (b) density (21.07.2000) (c) changes in patch growth of P. perfoliatus between 17.05.2001 and 25.07.2001

20 Lake and macrophyte characteristics

essed by a physically

r

e t

me ersität Hohenheim). s

te u

s

hy

nt na p e s oot

o t s

r h t o m he c o o s di o a ic e m h sh u m s t s r d e di g Re

2 r o

f e n ge sh m lo ba r o e 200 nd y l bm a l u Is Ju S

r

e t me s e t u y

s t t

na o o oph

ts

s

r h t he o men c o o di o ic h e h e

um s t s s r d ma e s

R di g o h ge s me lon bar of r 001 2 nd y l bme a l u Is Ju S

r

e t me es u yt

s a t

h n p oot

o ts s

h r o men he i c o d oot a ic

h e h um s m i t s r e s d

Re 0 o f e ng s

h s med lo bar o Patch patterns and patch dynamics over a period of 3 years, north-west Reichenau. Data acquired by an airborne scanner proc erg 200 nd m y l b a l Is Ju Su based algorithm. The ‘golf club’ is a shallow area connected to the mainland by man made earth dam (graphic: A. Woithon, Univ Fig. 2.8

21 Chapter 2

22 Spatio-temporal dynamics

Chapter 3. SPATIO-TEMPORAL DYNAMICS AND PLASTICITY OF CLONAL ARCHITECTURE IN POTAMOGETON PERFOLIATUS

Susanne R. Wolfer & Dietmar Straile

Aquatic Botany 78, 307-318 (2004)

Abstract

Aboveground and belowground clonal growth architecture was compared for two neighbouring patches of Potamogeton perfoliatus L. in Lake Constance, representing typical sparse and dense growth types. A detailed map of individual ramets and their corresponding rhizome network was produced and the seasonal development of ramet sprouting and rhizome growth was reconstructed. Rhizomes extended at rates ranging between 40 cm y-1 and 63 cm y-1, and added a new shoot for every 1 cm to 20 cm of rhizome produced. The main rhizome axis grew in a semi-linear fashion, with deviation means of 15° and 24°, and developed 0 to 0.5 branches per plant, with an insertion angle from 15° to 90° degree. Total rhizome length of all plants amounted to 5 m m-² and 11 m m-² at the two sites. The total biomass produced at the two sites differed ten-fold. The neighbouring patches also showed strikingly different allocation patterns which are interpreted as foraging.

Contrary to the prediction of the foraging hypothesis, there was only a small difference in mean spacer length between sites. On the other hand, the spacer length was surprisingly plastic relatively to total biomass. We propose that spacer length is not an independent part of the foraging strategy, but rather a result of overall productivity and biomass allocation.

Phosphorus content of plant tissue at the more productive site was two-fold higher than that at the less productive site which suggests that the different growth types are due to different nutrient availability in the sediment. In-situ fertilization of a less productive site increased the P content to the levels of the productive site.

In-situ fertilization also resulted in higher ramification, lower root allocation and decreased spacer length and confirmed the foraging capability of Potamogeton perfoliatus.

Keywords: macrophyte; clonal architecture; rhizome; biomass allocation; foraging, spatial growth

23 Chapter 3

Introduction

Macrophytes usually occur in patches which play a vital role in structuring of the littoral zone of lakes (Jeppesen et al. 1997). Potamogeton perfoliatus L. forms such spatial patterns ranging from semi- continuous to widely dispersed patches over a wide range of scales (Walser 1996, Lehmann 1997). In many lakes, Potamogeton patches change location, size, shape and density in the long run (Schmieder 1997) but remain spatially stable in the medium range (Walser 1996, Gafny & Gasith 1991), whereas in other lakes, they disappear and reappear from year to year (Scheffer et al. 1992). Presence or absence of macrophytes, as well as size, shape and density of vegetation patches influence element cycles (Barko & James 1997) and littoral food webs (Lauridsen et al. 1996). Despite the recognition of the importance of spatial patterns for the littoral ecosystem, it is not yet understood how these patterns are generated (Lehmann 1997) and which role can be attributed to architectural growth rules, morphological plasticity in response to spatially heterogeneous environmental factors, and competition between species or genets.

The concept of plant architecture suggests that each plant species has its own growth form (Bell 1984). Compared to non-clonal plants, clonal plant architecture comprises additional growth rules such as number of ramets, branching patterns (frequency and angle) of the rhizome, and the rhizome spacer length in between consecutive shoots (Callaghan et al. 1990, Huber et al. 1999). These rules are plastic to a certain degree allowing for adaptive responses to a spatially heterogeneous environment (Marbà & Duarte 1998).

Since seeds are considered as relatively unimportant for propagation of Potamogeton (Walser 1996), the knowledge of its architectural programme is essential for understanding the different rates and patterns at which this plant occupies space. Information on clonal growth rules is also important for the understanding and management of recolonization and restoration of littoral biotopes. Whereas descriptive mappings of aquatic vegetation patches have a long tradition, the importance of modular structure in studies of macrophyte stands has only recently been acknowledged. Most previous work on aquatic plant architecture has been published for seagrasses (e.g. Marbà & Duarte 1998, Molenaar et al. 2000) and emergent macrophytes (e.g. Clevering & Hundscheid 1998, Klimes 2000), whereas the clonal morphology of submersed freshwater macrophytes has received less attention (but see Idestam-Almquist and Kautsky 1995). In general, architectural studies lack spatial information, and vice versa, furthermore investigations on temporal aspects of clonal architecture and spatial pattern formation are rare. In all above-mentioned fields, comprehensive studies on P. perfoliatus are missing.

Below, we present detailed information on the characteristics of clonal growth architecture (length of rhizome spacers, branching intensity, branching angle) of P. perfoliatus, temporal dynamics, and natural variability of two typical growth types in Lake Constance. One site was characterized by lush patch growth with shoots reaching the water surface (referred to as ‘rich growth’ below) whereas at the other site, the shoots remained low and sparse (referred to as ‘poor growth’ below). The conclusions of the study are supported by an in-situ fertilization experiment and will be discussed in the light of foraging theory.

Materials and methods

Study site and general definitions The growth architecture of P. perfoliatus was investigated at the ‘Lower Basin’ (‘Untersee’) of Lake Constance (9°18'E, 47°39'N), a large, subalpine lake on the borders of Germany, Switzerland and

24 Spatio-temporal dynamics

Austria. The sub-basin has a surface area of 62 km², a mean depth of 13 m and a maximum depth of 40 m (Wessels 1998) and holds an extensive littoral area. P. perfoliatus is the second most common submersed macrophyte species in the area (Schmieder 1997), has a local vegetation period from approximately May to September and survives the winter in the form of subterraneous turions. The patches with sampling sites A and B (1 m² each) were located in a bay south of the north-western edge of the island of Reichenau, at a depth of approximately 2.2 m ± 0.2 m throughout the growth period. The distance between the two monospecific stands was 30 m.

For description of patch characteristics, the following terminology will be used: 'Plant' is defined as a complete unit of ramets connected by rhizomes, originating from a single primary shoot. 'Ramet' is a single module of a clonal plant, consisting of shoot, rhizome and roots. Spacer length is the rhizome length between two consecutive shoots of the same plant (Fig. 3.1).

ramet

main rhizome turion axis

branch

roots spacer length

Fig. 3.1 Apical section of P. perfoliatus plant. Rhizome diameter increases in growth direction. Turions are formed at main axis and side branches. The younger ramets are more stout than the older ramets

In-situ growth monitoring

On May 01 2000, when the first primary shoots were visible, two square frames with 1 m side length and a 10 x 10 cm grid subdivision were placed on the lake sediment at site A and B. The position of individual shoot bases within the grid and their corresponding lengths were noted on PVC-boards by snorkelers on four occasions in May, June and August 2000. The quadrats represent open systems: plants that had their origin outside could grow into the quadrats and plants from inside were allowed to grow out of the quadrats.

Architectural analysis Final harvest was made by SCUBA divers in August 2000. Each shoot was labelled at the top as well as close to the sediment. A transparent board (1 m²) placed above the sampling quadrat was used to draw a map of shoot positions. The shoots were cut off directly above the sediment and harvested with the upper set of tags before the sediment was removed by hand up to a depth of 15 cm. Rhizome networks were removed as completely as possible with the lower set of tags attached.

25 Chapter 3

In the lab, shoot length, ramifications and number of ramets per plant were recorded. All plants were cleaned and separated into shoots, rhizomes and roots. Turions were discarded because owing to their deep position in the sediment, they could not be harvested quantitatively. Plant fractions were dried at 105 °C until constant weight. After allowing for cooling down to room temperature, dry weights were determined on an analytic scale. Ramet number, branching frequency, branch angle, ramet length, a subset of spacer lengths, total rhizome length, specific rhizome weight, and biomass allocation values were determined for all plant material included in the m². For calculations of plant size, plant rhizome length, and radial angle, only those plants growing completely within the quadrat were considered. Spacer length was calculated as Euclidian distance between the nearest sibling ramets. Calculated spacer lengths compared adequately with a subset of measured spacer lengths (r = 0.98, p < 0.001).

The drawings of shoot positions and rhizome networks were scaled down by transformed coordinate systems. The mapping of shoots over the investigation period and the additional information about connections between these shoots made it possible to assign daughter shoots to parent shoots and to reconstruct the spatio-temporal growth of individual rhizome systems. A rough estimate of the temporal dynamics of biomass allocation was based on the number of shoots, mean shoot length and the assumption of constant weight-specific allocation into shoots and ramets. Angles were determined to the nearest degree on the reconstructed map.

Fertilization experiment and tissue analysis On July 02 2001, five typical poor growth areas of 0.25 m² were fertilized (‘fertilized type A’) with 250 ® g slow-release N-P-K-fertilizer Plantacote Depot 4 M (14 % N, 9 % P2O5, 15 % K2O), additional three low growth sites were assessed as control (‘type A’). All sites had an initial density of eight shoots per m² and were located in the same lake area as A and B, at a maximum distance of 50 m between each other. On August 15 2001, the complete plants including rhizomes were removed from the sediment and processed in the same way as those from locations A and B. Leaf tissue of the fertilized and poor growth sites as well as additional rich growth sites (‘type B’) was analysed for carbon and nutrient content. C was determined by Autoanalyzer (ThermoQuest NCS-2500) after combustion. P and N values were determined by AutoanalyzerII (Bran & Luebbe) following digestion with potassium persulphate. For the statistical analysis of percentages, data were arcsin transformed. Statistical analyses were performed with SAS software (SAS Institute Inc., Cary, NC, USA 1999-2001).

Results

General growth characteristics At the first field visit, ramet numbers were identical, but during the following two months, ramet numbers were always higher at site B than at site A. The sites also showed a strikingly different development in several other aspects (Fig. 3.2). At site A, biomass at harvest was more than ten-fold lower (17 g DW m-2) than at site B (200 g DW m-2). The maximum density was 90 ramets at site A compared to 210 at site B (Fig. 3.3a). The calculated number of plants was 18 with 5 ± 1 ramets each at site A, compared to 25 plants with 8 ± 3 ramets each at site B. On average, shoots measured 12 cm ± 7 cm at site A compared to 66 cm ± 75 cm at site B (Fig. 3.3b). At site B, many shoots had reached the water surface and formed a canopy (285 cm max. shoot length), whereas at site A the longest shoot measured only 48 cm.

26 Spatio-temporal dynamics

Fig. 3.2 Reconstructed spatio-temporal growth dynamics of P. perfoliatus at poor growth (A) and rich growth (B) site. Points represent ramets, lines show rhizome connections, arrows indicate growth direction, double arrows show sites of turion formation

27 Chapter 3

Fig. 3.3 (a) Temporal development of ramet numbers and (b) average shoot length ± S.E. of P. perfoliatus at poor growth (A) and rich growth (B) site

Growth architecture Primary P. perfoliatus ramets developed plagiotropic stems (lateral rhizomes) initiated at the ramet base which produced further ramets on every second rhizome node and adventitious roots at every node (Fig. 3.1). Overwintering organs, also called turions, were formed at the end of the main rhizome as well as on most branch tips. Rhizome branching was acropetal and occurred only at the base of ramets. The maximum level of branching encountered at the research site was 1st order. Remarkably, there were some primary ramets bearing two rhizome axes at two different vertical nodes. Rhizome branch growth was subordinate to the growth of the main rhizome axis. No true branches were found at site A, whereas at site B every second plant (every 17th ramet) had developed a branch. Initial rhizomatous growth direction appeared to be random. The mean radial angle was 15° ± 12° at site A and 24° ± 25° at site B, and the mean insertion angle (branch angle) was 68° ± 15° at site B (15° to 90°). Especially at site B, rhizomes of different plants crossed each other frequently.

Fig. 3.4 (a) Frequency distribution of spacer length, values on x-axis indicate midpoints of intervals (b) spacer length ± S.E. between primary and higher order shoots (main axis only) of P. perfoliatus at poor growth (A) and rich growth (B) site

28 Spatio-temporal dynamics

The spacer length was highly variable at both sites ranging from 1 cm to 20 cm of rhizome produced between two shoots (Fig. 3.4a). However, average spacer lengths were similar with 8.0 ± 3.2 cm at site A and 7.6 ± 3.4 cm at site B. Part of the spacer length variation can be attributed to the order of shoots (Fig. 3.4b). At site B, spacing increased in the direction of new ramet recruitment. Average distance between ramet 6 and 7 (11 cm) was almost three times the distance between primary and secondary ramets. In contrary, spacer length at site A approached a constant already between the secondary and tertiary ramets. Spacer order (ANOVA: F1,162= 48.5, p < 0.0001) and the square of spacer order (F1,162 = 24.1, p < 0.0003) explained 37 % of the variability of spacer length. The significant contribution of the square root of spacer length indicates a nonlinear relationship between the length and order of spacers. After accounting for the effects of spacer order, analysis of covariance shows that the two sites differ significantly in spacer length (F1,160= 12.3, p < 0.001) but also in the increase of spacer length with spacer order (F1,160= 10.4, p < 0.001). The site effect, however, increases the adjusted r2 of the model by only 4 %, showing its secondary importance compared to the effects of spacer order.

Final total rhizome length was 5.4 m at site A compared to 11 m at site B. The lower weight-length relationship of rhizomes at site A (6.9 mg cm-1) compared to site B (7.7 mg cm-1) was probably due to a smaller rhizome diameter at site A. Mean rhizome growth of individual plants ranged between 40 cm yr-1 (site A) and 63 cm yr-1 (site B) and determines the maximum extension of the outer perimeter of a patch.

Fig. 3.5 (a) Shoot growth rate ± S.E. and (b) rhizome growth rate ± S.E. per plant of P. perfoliatus at poor growth (A) and rich growth (B) site

Biomass allocation At the time of harvesting, sites A and B differed considerably in biomass allocation. Shoot-, rhizome- and root biomass amounted to 73 %, 22 % and 5 % of total dry weight at site A compared to 95 %, 4 % and 1 % at site B. The estimation of temporal dynamics of biomass allocation (Fig. 3.5a,b) shows that the initial shoot growth rates were similar. At the second field visit, the shoot growth rate at site B was several times higher than at site A (Fig. 3.5a) which implies a higher total photosynthetic capacity. Despite the large difference in photosynthetic tissue, the rhizome growth rates showed no significant difference between the two sites (Fig. 3.5b).

29 Chapter 3

Fertilization experiment and tissue analysis Fertilization of plants resulted in longer shoots, almost three-fold biomass (t-test, p < 0.08) and a doubling in shoot number (t-test, p < 0.08). Mean rhizome length was shorter in the fertilized treatment (5.1 cm) than in the control (6.4 cm) (t-test, p < 0.06). The unfertilized plants had no branches, whereas under fertilization, an average of 3.6 branchings per quadrat was found. Shoot allocation was higher (t-test, p < 0.03) and rhizome and root allocation were lower (t-test, p > 0.05 and p < 0.01) in the fertilized quadrats than in the controls (Tab. 3.1).

Tab. 3.1 Biomass allocation (% dry weight) of P. perfoliatus in poor growth (type A, n=3) and fertilized poor growth (fertilized type A, n= 5) quadrats

Shoot ± S.D. (%) Rhizome ± S.D. (%) Root ± S.D. (%)

Type A 64 ± 4 29 ± 1 7 ± 3

Fertilized type A 80 ± 10 17 ± 9 3 ± 1

N content (% dry weight) was lower in type A plants than in fertilized type A plants (Tukey’s HSD test, p < 0.005), P content (% dry weight) was lower in type A plants than in fertilized type A plants and rich growth type B plants (Tukey’s HSD test, p < 0.001), and C:N, C:P and N:P weight ratios of leaf tissue were higher in type A plants than in fertilized type A plants and rich growth type B plants (Tukey’s HSD test, all p < 0.006) (Tab. 3.2).

Tab. 3.2 N and P contents (% dry weight), C:N, C:P and N:P weight ratios of P. perfoliatus leaf tissue in poor growth (type A, n=3), fertilized poor growth (fertilized type A, n=5), and rich growth (type B, n=5) quadrats

% N % P C:N C:P N:P

Type A 1.7 ± 0.4 0.14 ± 0.03 22.4 ± 4.4 268.4 ± 55.3 12.2 ± 2.5

Fertilized type A 3.1 ± 0.8 0.36 ± 0.11 13.8 ± 3.9 120.5 ± 39.8 8.7 ± 0.9

Type B 2.5 ± 0.3 0.35 ± 0.11 16.1 ± 1.8 120.7 ± 29.1 7.5 ± 1.6

Discussion

The results of our study show a common architectural base-programme at two neighbouring stands of P. perfoliatus, but also striking variability in some growth characteristics including ten-fold total biomass, two-fold density, five-fold shoot length, two-fold total rhizome length, higher branching frequency, and higher biomass allocation to shoots and less to roots at site B compared to site A.

The unequal nutrient concentrations of plant tissue at the two sites is a strong indication that growth differences were due to differences in sediment nutrient availability (Langeland et al. 1983, Hilbert 1990) which is known to show considerable spatial variability (Downing & Rath 1988). Causes such as water column nutrients, light conditions and currents can be ruled out because the two patches are located close by. Genotypic differences are unlikely because Potamogeton propagates vegetatively in Lake Constance (Walser 1996). Also the increase in tissue P in the fertilized quadrats suggests that

30 Spatio-temporal dynamics

type A plants were P limited. In fact, the P content of type A plants (0.14 %) is close to that considered as limiting for aquatic macrophytes (0.13 %) by Gerloff and Krombholz (1966).

Foraging theory considers changes in branching frequency, spacer length, and biomass allocation as mechanisms for the preferential location of modules in favourable microsites, and as means of an improved exploitation of resources such as light or nutrients (Hutchings 1997). Intraspecific clonal variability due to nutrient availability has been described for terrestrial plants (Cheplik 1995) and seagrasses (Pérez et al. 1994). As predicted by foraging theory, we found increased branching intensity at the rich site B as well as in the fertilized quadrats. Branching intensity is positively related to nutrient supply in many terrestrial clonal plants (Hutchings & De Kroon 1994), aquatic freshwater plants (Evans 1988), and seagrasses (Terrados et al. 1997a) and is thought to maximize the number of ramets at a suitable site.

As predicted by foraging theory, plants at the nutrient-deficient site allocated considerably more biomass to roots than plants at the rich site, and plants in the controls of the fertilization experiment allocated more biomass to roots than the fertilized plants. This allocation rule has been confirmed by a variety of experiments with clonal terrestrial plants (review of Hutchings 1997) and submersed macrophytes (Kautsky 1991, Reusch et al. 1994). Recently, it has been proposed that allocation differences under different nutrients budgets are rather a consequence of allometry than of morphological plasticity, as smaller plants frequently show a higher allocation to belowground structures (Gedroc et al. 1996, Müller et al. 2000). Site A ramets, however, have ceased height growth in favour of producing longer rhizomes and additional modules and therefore represent an architectural model different from that of site B. Hence, the observed differences in biomass allocation are apparently not the result of allometric differences, but rather the result of different growth strategies.

Mean spacer length did not differ strongly between the two sites. According to the foraging hypothesis, however, plants in good habitats (B) should have shorter spacers in order to place a maximum amount of shoots in the favourable habitat while plants in less favourable habitats (A) should grow long spacers in order to escape from the bad site. Although some investigations confirm the hypothesis for rhizomatous species, other researchers found indifferent or even longer rhizomes in more favourable habitats (De Kroon & Hutchings 1995). To explain the latter results, Hutchings and De Kroon (1994) proposed the growth hypothesis (see also Stoll et al. 1998), according to which spacer length should increase at favourable sites as they allow a higher productivity. The increase in spacer length with the order of ramets at both stands is, for example, in accordance with the prediction of the 'growth hypothesis' and can be interpreted as a result of an increasing overall productivity of a rhizome system with increasing amount of photosynthetically active and clonally integrated ramets. The lack of distinct differences between the mean spacer length at sites A and B, however, can neither be explained by growth hypothesis nor foraging theory alone. We therefore suggest that a synthesis of both theories is necessary for a complete understanding of spacer length in Potamogeton. Despite having only one tenth of photosynthetically active biomass, plants at the poor site produced the same average spacer length and more than sixty percent of the total rhizome length compared to site A plants. Principally, longer rhizomes at inhospitable sites can be achieved by means of two mechanisms: (a) by an increased biomass allocation into rhizome tissue, as seen in our study, and (b) by reducing the tissue-mass density and thereby increasing spacer length. Only in case (b), spacer length can be considered as an independent part of the foraging behaviour of the plant, whereas in our study spacer length is mainly an effect of biomass allocation into rhizomes. Such foraging sensu latu can easily be overseen when biomass is ignored. We therefore emphasize that

31 Chapter 3

studies which investigate the foraging capacity of clonal plants need to include measurements of spacer length, biomass and biomass allocation. In case of the fertilization experiment, foraging was more obvious, because due to comparable biomass, mean spacer length tended to be shorter and biomass allocation into rhizomes was lower in the fertilized plants than in the controls.

Overall, plants must function as a balanced system in terms of resource uptake and use, and resources may not be allocated independently from each other (e.g. Hirose 1986; Ågren & Ingestad 1987). On the assumption of different nutrient regimes at the two sites, the mechanisms for biomass allocation may be interpreted in the light of stochiometric theory. Under nutrient limitation, increased photosynthesis as a consequence of higher biomass can not increase biomass production because of missing nutrients. Consequently, it would be better for the plant to invest into roots to capture more nutrients, and rhizomes in order to grow faster out of a limiting site. Only if a high nutrient uptake rate makes C the limiting factor, it is beneficial for the plant to produce aboveground photosynthetic tissue instead of belowground biomass. At our study site, the allocation response of P. perfoliatus to sediment nutrient limitation allows for a patch expansion of roughly two thirds of that under thriving conditions, although the biomass is only one tenth. However, these patches have much shorter shoots and are less dense due to the lower number of plants per area, lower number of ramets per plant, and specific clonal architecture of the species (marginally longer spacers and less branchings). Hence our study shows how plasticity of clonal architecture in P. perfoliatus may result in two patch types that constitute two distinctively different habitats within the littoral biocoenosis.

32 Spatio-temporal dynamics

33 Chapter 3

34 Density control

Chapter 4. DENSITY CONTROL IN POTAMOGETON PERFOLIATUS L. AND POTAMOGETON PECTINATUS L.

Susanne R. Wolfer & Dietmar Straile

Limnologica 34, 98-104 (2004)

Abstract

The effect of density on the growth, recruitment of new ramets, biomass allocation, rhizome spacer length and rhizome branching of the submersed macrophyte Potamogeton perfoliatus L. was experimentally evaluated in a mesocosm with three different initial shoot densities. The findings suggest that the number of primary shoots available in the beginning of the season can strongly influence patterns of growth and clonal reproduction. In contrast to many studies that found decreasing production parameters with plant density, number of ramets and total biomass of P. perfoliatus plants were highest at medium plant density, and lower at low density and at high density. This is probably due to negative feed-back through crowding at high density and unclear negative effects at low density. Shoot allocation tended to increase with density, rhizome allocation tended to decrease with density and root biomass remained unchanged. At low density P. perfoliatus produced longer rhizome spacers and more branchings than at higher densities. Shorter rhizome spacers at high plant density probably restrict patch expansion and cause discrete patch shapes. A likely mechanism for plastic changes in clonal architecture is increasing competition for light at increasing density, but other density-dependent factors may also contribute.

Investigations on propagule numbers of tubers and turions show that self-regulation of ramet number was associated with self-regulation of propagule number in both P. perfoliatus and Potamogeton pectinatus, with higher densities producing less propagules per plant.

Key words: Potamogeton, submersed macrophyte, density, self-regulation, propagules

Introduction

Potamogeton perfoliatus L. and Potamogeton pectinatus L. are common clonal macrophytes in many parts of the world. Natural stands of P. perfoliatus have densities from one shoot / m² to several hundred shoots / m², Potamogeton pectinatus densities may reach more than 1000 shoots / m² (own observation). Since Potamogeton propagates mainly clonally, shoot density is determined by (a) the number of sprouting propagules (= plant density), (b) the number of ramets per plant, and (c) the rhizome length between shoots. The number of sprouting propagules is a product of shoot growth and biomass allocation into reproductive structures of the previous growth season, and propagule survival during winter. In aquatic macrophytes, morphological characteristics like the number of ramets per plant, spacer length and branching may react plastic on environmental factors such as water depth (Strand & Weisner 2001), light (Barko et al. 1982), nutrient availability (Barko et al. 1991), or wave action (Idestam-Almquist & Kautsky 1995).

The shoot density of macrophytes is important for processes in the littoral ecosystem. Macrophyte density or patchiness determine the distribution and age-characteristics of littoral fish (Weaver et al.

35 Chapter 4

1997) and can affect macro-invertebrates (Webster et al. 1998). Density determines the competitive ability of aquatic clonal plants in mixed stands (Davis & Fourqurean 2001), and fluctuations in density can facilitate species coexistence by delaying the exclusion phenomenon (Bonis et al. 1995). Density may have intraspecific feed-backs such as light and nutrient modification, or regulation of flowering (Thompson et al. 1990).

Most studies of macrophyte density are confined to shoot numbers but ignore details on larger hierarchical levels, e.g. clone fragments or clones, as well as belowground dynamics. We experimentally evaluated the effect of initial ‘propagule’ density on the growth, recruitment of new ramets, biomass allocation, spacer length and branching of plants of the submersed macrophyte P. perfoliatus in a mesocosm with three monospecific primary shoot densities (corresponding to plant number). Turions of Potaomgeton perfoliatus and tubers of Potamogeton pectinatus were also investigated, because the number of propagules determines the initial density in the following year and is therefore crucial for the modelling and prediction of patch dynamics. Furthermore the weight of propagules influences growth and survival of the produced shoots (Idestam-Almquist & Kautsky 1995, Santamaría & Rodriguez-Girones 2002).

The main hypotheses to be tested were: (1) Ramet production decreases with density in P. perfoliatus (self-regulation hypothesis) (2) Propagule production decreases with density in P. perfoliatus and P. pectinatus.

Materials and methods

The mesocosm experiments were conducted in a basin of 10 m x 5 m x 1.2 m on the premises of the Limnological Institute (University of Constance), nearby Lake Constance, Germany. The distribution of about 8 m³ of lake sediment originating from the lower lake basin created a sediment layer of approximately 15 cm depth.

On May 31, 2000, P. perfoliatus and P. pectinatus shoots from nearby natural stands (P. pectinatus were collected near Egg, P. perfoliatus near Unteruhldingen) were planted in three monocultural densities (5, 15 and 40) per 0.25 m² in three replicates. The terms ‘low’, ‘medium’ and ‘high’ density are used in a relative sense on the potential scale for the number of initial shoots (= plant density). All shoots had a shoot length of 15 – 20 cm, and a maximum rhizome length of 10 cm. The quadrats were located within 3 m distance of each other to avoid interaction. Treatments were randomly located within the mesocosm and plants were randomly rooted within the quadrats.

After planting, the mesocosm was carefully filled with lake water to 1 m depth. A pump of a capacity of 3 l/s maintained a constant flow-through with a renewal time of roughly one day. Since the lake water originated from 15 m depth, excessive heating of the water was avoided. Initial mortality of transplants was < 1 % after 8 days. Shoots grew well and at the end of the experiment, part of the shoots from all treatments topped at the water surface.

During the growth period, average length was estimated by measuring a subsample of ramets. On August 4, the water was drained and the rhizome systems were carefully dug out by hand. Turions were present and showed that the growth phase was already completed. Primary, secondary and higher order ramets were counted, spacer length was measured. Plants were separated into shoots, rhizomes and roots, and dried at 105 °C to constant weight. The fractions were weighed on an analytic scale.

36 Density control

Turions could not be harvested quantitatively, because they were located deeply in the sediment and broke off frequently during digging. Since turions are built at all main axes and nearly all branch tips (Wolfer & Straile 2004b), the number of plants and the number of branchings was used to estimate the number of turions instead.

Since the patch perimeter was extremely irregular with single runner shoots far from the centre, patch size was estimated as (a + (2 * r * b))² where a is the side length of the quadrat in meters, r is the average total rhizome length of the plants, and b is a factor correcting for branching (b low density =

2.8, b medium density = 2.5, b high density = 1.8). ‘Spatial ramet density’ (the number of ramets per square meter) was calculated as final ramet number divided by patch size.

Clones of Potamogeton pectinatus were not harvested, but tubers were collected by sieving the sediment through a mesh of width = 1 mm. One replicate of the density 15 treatment had to be omitted because of harvesting problems.

Statistical analyses were performed with SAS software (SAS Institute Inc., Cary, NC, USA 1999- 2001). Proportional data were transformed to the arcsine of the square root before analysis.

Results

Length growth dynamics

The first P. perfoliatus plants reached the water surface after 49 days which is an average growth rate of 16 mm per day. In the beginning, shoots at higher density grew faster in length than at low density (results not shown). At harvest, average length per treatment was not density dependent any more (Anova, F2,8 = 1.34; p < 0.33).

Production parameters

P. perfoliatus produced between 30 and 343 ramets per site, variability within treatments was considerable (Fig. 4.1a). Ramet recruitment (Fig. 4.1b) was highest at medium plant density (Tukey’s

HSD-test, F1,8 = 8.7, p < 0.02). Total biomass production per initial shoot was also highest at medium plant density (Fig. 4.1c). The average specific shoot dry weight of 0.85 g m-1 was not different between treatments (Anova, F2,8 = 0.5, p < 0.6).

Fig. 4.1 (a) Average number of ramets produced per site ± S.E. (b) average number of ramets recruited per initial shoot ± S.E. (c) average total dry weight production per initial shoot ± S.E. of P. perfoliatus grown in a mesocosm at different densities

37 Chapter 4

Biomass allocation and spacer length

Shoot allocation of P. perfoliatus tended to increase with initial shoot density and rhizome allocation tended to decrease with density. The final shoot density can, however, better explain shoot allocation (r² = 0.63, n = 9, p < 0.01) and rhizome allocation (r² = 0.74, n = 9, p < 0.003) than initial density. Root allocation was not affected by final density (r² = 0.001, n = 9, p < 0.93) (Fig. 4.2).

Fig. 4.2 Average shoot, rhizome and root biomass allocation (%) ± S.E. of P. perfoliatus grown in a mesocosm at different densities

The decrease in rhizome allocation with increasing initial density goes along with a non-significant decrease in spacer length (Anova, F2,8 = 1.03, p < 0.4) and a non-linear relationship between shoot number and total rhizome length (r² = 0.98, n = 9, p < 0.0001, quadratic term: p < 0.04) (Fig. 4.3a,b). Since medium plant density produced more ramets, average rhizome length produced per initial shoot was highest at medium plant density, medium at low plant density and lowest at high plant density, all were significantly different (Anova, F2,8 = 16.5, p < 0.004) (Fig 4.3c). The average specific rhizome dry -1 weight of 0.5 gm was not different between treatments (Anova, F2,8 = 0.8, p < 0.5).

Fig. 4.3 (a) Average rhizome spacer length ± S.E. (b) relationship between total rhizome length and final number of ramets per site (c) average total rhizome length per plant ± S.E. of P. perfoliatus grown in a mesocosm at different densities

38 Density control

Branchings

Rhizome branching probability per ramet and number of rhizome branchings per plant of P. perfoliatus were affected by density (Anova, F2,8 = 4.24, p < 0.07, respectively F2,8 = 7.4, p < 0.02) (Fig. 4.4). The number of rhizome branchings per plant was significantly lower at high density than at low and medium density (Tukey’s HSD-test, F1,8 = 8.9, p < 0.02) (Fig. 4.4b).

Fig. 4.4 Average number of rhizome branchings (a) per ramet ± S.E. (b) per plant ± S.E. of P. perfoliatus grown in a mesocosm at different densities

Patch size and spatial density

Patch size of P. perfoliatus was higher at medium plant density than at medium and at low plant density (Tukey’s HSD-test, F2,8 = 16.7, p < 0.004) (Tab. 4.1). Spatial density was higher at high density than at medium and low density (Tukey’s HSD-test, F2,8 = 15.7, p < 0.004) (Tab. 4.1).

Tab. 4.1 Calculated patch sizes and spatial densities of P. perfoliatus plant density patch size spatial ramet turion numbers tuber numbers per of P. (m²) of density (m-²) of per P. perfoliatus P. pectinatus plant perfoliatus P. perfoliatus P. perfoliatus plant

5 2.2 20 2.1 ± 0.75 8.0 ± 3.1

15 4.8 42 2.0 ± 0.10 4.6 ± 0.2

40 1.9 138 1.3 ± 0.25 2.0 ± 0.9

Propagules

Turion numbers per plant of P. perfoliatus were affected by density (Anova, F1,8 = 5.4, p < 0.05), and were lowest at high plant density and higher at low and at medium plant density. Tuber production of P. pectinatus decreased with density (Fig. 4.5a). Numbers produced at low density were significantly higher than those at high density (Tukey’s HSD-test, F1,7 = 6.51, p < 0.04). Propagule production per plant was roughly 4-times higher at low density, and 2-times higher at medium and high density in P.

39 Chapter 4

pectinatus compared to P. perfoliatus (Tab. 4.1). Average fresh weight per tuber was not significantly different between treatments (Anova, F2,7 = 0.99, p < 0.44). All fresh weight frequency curves were skewed to the right and leptokurtic (Fig. 4.5b).

Fig. 4.5 (a) Average number of tubers ± S.E. (b) fresh weight frequency distribution of tubers of P. pectinatus plants grown in a mesocosm at different densities (n5, 40 = 3; n15 = 2).

Discussion

The findings suggest that initial density can strongly influence patterns of growth and clonal reproduction of P. perfoliatus, by changing (a) relative ramet recruitment (b) biomass allocation (c) branching and (d) spacer length.

The hypothesis of negatively density-dependent growth was only partly confirmed by the results. In accordance with the hypothesis, ramet production and total biomass production per initial shoot were lower at high density than at medium plant density. Reduced rates of ramet recruitment and growth are the most common effects of increased density in terrestrial plant species (De Kroon & Kwant 1991, De Kroon 1993) and higher aquatic plant species (Moen & Cohen 1989, Groth et al. 1996). Regulation can take place via mortality control (Lovett-Doust 1981) or natality control (Briske & Butler 1989, Thompson et al. 1990). The investigation of rhizome nodes of our experimental plants revealed no loss of ramets, suggesting a predominantly natality-regulated ramet recruitment in P. perfoliatus. Hutchings (1979) proposed intraclonal regulation as underlying reason for self-regulation of ramet density, whereas other authors favour a negative feed-back of deteriorating conditions (e.g. self- shading, nutrient depletion) due to overcrowding (De Kroon & Kwant 1991).

Contrary to the self-regulation hypothesis we found less ramet production and growth at low density than at medium density. Cain et al. (1995) also found highest performance of a terrestrial clonal plant at medium density and Worm and Reusch (2000) report that reduced plant density had negative effects on eelgrass shoot growth during the early stages of recolonization. Reduced fitness at low population size (‘allee effect’) is often caused by genetic deficiencies (Fischer et al. 2000), but low population sizes may also involve growth impairment by edge effects (Gardner 1998), increased levels of epiphytes (Reed 1990), excessive irradiance (Scrosati & Dewreede 1998), or stochastic mortality (Dennis 2002). We propose that relative ramet recruitment and biomass production in P. perfoliatus increases at low ramet densities, reaches an optimum where benefits and disadvantages of density- related factors counterbalance each other, and decreases at high densities.

40 Density control

Density did not only influence growth but also the clonal architecture of P. perfoliatus: Shoot allocation tended to increase with density and rhizome allocation tended to decrease with density. The same results have been reported for terrestrial plants (Bazzaz & Grace 1997). In our study, shoot allocation lay well in the range reported for P. perfoliatus in the literature (58-95%, Ozimek et al. 1976) and from in-situ investigations in Lake Constance (Wolfer & Straile 2004b). Although researchers often found lower root allocation of Potamogeton with increasing density (Moen & Cohen 1989, Sand- Jensen & Vindbaek Madsen 1991), our results show no effect of density on root allocation. Root allocation probably indicates that nutrient supply is adequate, otherwise the plants should react plastic on nutrient shortage (Aerts et al. 1991). The reason for the observed differences in shoot and rhizome allocation can be strategies towards different shading conditions (e.g. Ryser & Eek 2000). In case of increased self-shading due to higher density it would be better for the plant to invest (a) in rhizomes in order to grow faster out of the shaded site (foraging sensu strictu) at the cost of shoots or (b) shoots to catch more light (foraging in a wider sense; ‘light foraging’) at the cost of rhizomes. Strategy (a) might be problematic because shoots in high density stands already have to cope with a reduced net photosynthesis as a result of reduced irradiance (Scrosati & Dewreede 1998), and might therefore lack the energy to built longer rhizomes. In fact, the plants in our study showed strategy (b). This is in contradiction to the foraging hypothesis which proposes that plants in ‘good’ habitats (= low density, no shade, nutrients) should have shorter spacers in order to establish a maximum amount of shoots in the favourable habitat while plants in less favourable habitats (= higher density, self shading, possibly nutrient limitation) should grow long spacers in order to escape from the bad site (e.g. Sutherland & Stillman 1988). Field investigations, however, showed that only some plants forage by means of spacer length (Hartnett & Bazzaz 1985, Hutchings & De Kroon 1994), and that foraging might be plastic and depends on the relative levels of resources available (Wolfer & Straile, unpublished data).

Density affected branching, where plants of the low density treatment showed higher branching frequencies per ramet than plants of the medium and high density treatment. Experiments have shown that self-shading reduces branching and light induces branching in terrestrial plant species (Ellison & Niklas 1987, Hutchings & De Kroon 1994) and seagrasses (Vermaat & Verhagen 1996, Marbà et al. 1996). Branching frequency is also interpreted as foraging behaviour, since increased branching maximizes the number of ramets at a suitable site (Sutherland & Stillman 1988). In our experiment, however, reduced biomass production at low density compared to medium density shows that patch quality was not optimal at low density, and nevertheless, branching was higher. As a consequence of different trends in ramet numbers and branching frequencies, low and medium density stands reach similar number of branchings per plant, whereas high density plants had less branchings. This leads to self-regulation of density and has consequences for the production of propagules (see below).

Patch expansion is partly ruled by the density-regulated architecture of clonal growth. Our estimations show lower patch expansion at low plant density (due to increased branching and decreased growth) and at higher density (due to shorter rhizome length and decreased growth), and highest patch expansion at medium density. Also the spatial ramet density was highest at high initial density. High density probably restricts patch expansion and causes compact patch shapes.

Besides self-regulation with regard to ramet production and branching, P. perfoliatus also showed self-regulation with regard to propagule formation. Since turions are formed on branches and are therefore a direct consequence of plant architecture, low and medium density plants produced more propagules per plant whereas high density plants produced less propagules. In P. pectinatus, the tuber production was negatively density dependent throughout. Since the tuber weight did not differ with density, the self-regulation of tuber numbers will lead to better starting conditions in the following

41 Chapter 4

year. The higher relative propagule production in P. pectinatus compared to P. perfoliatus is probably required because P. pectinatus tubers are subject to a high mortality due to bird foraging, fungus etc. (e.g. Hangelbroek et al. 2002) whereas turions are buried deeper in the sediment at water depths probably not reached by water fowl.

Summarizing, P. perfoliatus and P. pectinatus plants can lessen low density effects and overcrowding by self-regulation of ramet recruitment, growth and propagule production associated with changes in clonal architecture.

42 Density control

43 Chapter 4

44 Nutrient experiments

Chapter 5. EFFECT OF SEDIMENT STRUCTURE AND NUTRIENTS ON THE GROWTH OF POTAMOGETON PERFOLIATUS

Abstract

A series of assessments and experiments was carried out to test for effects of sediment nutrients on the growth of Potamogeton perfoliatus:

(1) Bare sediments and sediments with differing vegetation cover at Lower Lake Constance were tested for differences in grain size and nutrients. There was no relationship between sediment composition, nutrient content and macrophyte type. We suggest that this is due to sampling methods, lack of knowledge about the appropriate nutrient fractions to be measured in porewater, and the distinct seasonal characteristics of nutrient cycles.

(2) The growth substrate of P. perfoliatus was nutrient enriched with slow-release fertilizer in a mesocosm. The results were compared with in-situ enrichment results described in Chapter 3. Whereas in-situ fertilization increased ramet recruitment and biomass growth, fertilization of the mesocosm sediments showed no significant effect. As predicted by the foraging theory, ramets allocated less biomass to roots when fertilized and rhizome spacer length was significantly shorter in the fertilized plants. Shoot / rhizome allocation was not significantly different.

(3) P. perfoliatus ramets were planted in areas without previous macrophyte growth at Upper Lake Constance, with or without nutrient addition. Only fertilized stands could establish, showing that the occurrence of P. perfoliatus in parts of the Upper Lake is at least partly limited by nutrient availability.

The different responses of growth to fertilization are attributed to the different sediments used: when sediments were naturally nutrient poor, fertilization was beneficial, whereas on nutrient rich sediment, fertilization was indifferent.

Introduction

P. perfoliatus is a species which is known to react to changes of water quality. In Lake Constance, an originally oligotrophic lake, nutrient loads of macronutrients (N and P) increased until the end of the 1970ies due to agricultural runoff and increasing wastewater discharge. Since the construction of area-wide treatment plants, nutrient levels have been decreasing steadily and are now back to pre- eutrophication levels.

While most parts of Upper Lake Constance show little macrophyte growth, large stands are found at the influents of nutrient-rich rivers such as the Schussen. At Lower Lake Constance, P. perfoliatus has been a typical plant since the beginning of the 20th century. During the phase of eutrophication, the species decreased considerably but increased again after re-oligotrophication. In Upper Lake Constance, there was not much change over the years (Schmieder 1997).

Since macrophytes may take up nutrients from the water column over their surface (Sculthorpe 1967, review of Agami & Waisel 1986), water nutrient concentrations and macrophyte biomass are sometimes correlated (Carbiener et al. 1995, Bini et al. 1999). However, root uptake is dominant (Denny 1980, Barko & Smart 1986a, Barko & Smart 1986b) and suggests predominant relationships between sediment nutrients and biomass. Whereas water chemistry may change quickly, the sediment

45 Chapter 5

nutrient status reflects the recent chemical history of the water body. This means that nutrient concentrations in the interstitial water are commonly a better indicator of macrophyte growth.

Nutrients do not only influence biomass but may also affect the morphological characteristics of aquatic macrophytes. Under varying nutrient conditions, there may be variations in plant size, leaf number, leaf morphology, leaf and root turnover rate and storage capability (Aiken & Picard 1980, Hutchings & De Kroon 1994). Moreover, clonal architecture such as spacer length, rhizome allocation and branching may be affected (Wolfer & Straile 2004b; Chapter 3). These changes may be beneficial when they allow for a better usage of the resources available in the growth environment. The ability of stoloniferous and rhizomatous plant to concentrate ramets in nutrient rich patches or to avoid nutrient deficient patches by preferential placement of modules has been interpreted as ‘foraging’ behaviour analagous to the search path of a foraging animal. Rhizome spacer length is supposed to decrease, and branching is supposed to increase with patch quality (Bell 1984, Schmid 1986, Slade & Hutchings 1987a, Slade & Hutchings 1987b, Slade & Hutchings 1987c, Evans 1988, De Kroon & Schieving 1990, De Kroon & Knops 1990, Kelly 1990, Evans 1991, Oborny 1991, Hutchings & De Kroon 1994, De Kroon & Hutchings 1995).

There is ample literature on experiments concerning nutrient uptake and utilization by aquatic macrophytes (Barko & Smart 1979, Best & Mantai 1987, Barko & Smart 1981b, Carignan 1985, Christiansen et al. 1985, Anderson & Kalff 1986, Barko et al. 1991, Beklioglu & Moss 1996, Bini et al. 1999, Clarke & Wharton 2001). Fertilization experiments in the aquatic environment have been conducted for both seagrasses (Agawin et al. 1996, Alcoverro et al. 1997, Worm 2000, Worm et al. 2000, Peralta et al. 2003) and freshwater macrophytes (Daldorph & Thomas 1991, Demars & Harper 1998) but the effect on the clonal architecture of submersed macrophytes has hardly been studied.

The aim of the study was to assess the effect of sediment nutrients on presence or absence, growth and architecture of P. perfoliatus. The assessments and experiments comprised:

(1) grain size and nutrient analysis of bare lake sediments as well as lake sediments with differing vegetation cover;

(2) nutrient enrichment of growth substrate with a slow-release fertilizer in a mesocosm; analysis of clonal growth parameters and tissue nutrient concentrations;

(3) planting of P. perfoliatus ramets in lake areas without significant macrophyte growth, either without or with nutrient addition.

Materials and methods

Sediment analysis

All sites were situated in a bay W-SW of Niederzell, island of Reichenau, within an area of 50 m x 25 m, at a depth between 2.00 m and 2.50 m (see Fig. 2.1).

Three to five replicates of sediment at growth sites of each of the following 5 vegetation types were taken: (a) P. perfoliatus tall and dense (> 2.50 m; > 50 shoots/m²), identical to rich growth sites of Chapter 3, (b) P. perfoliatus tall and sparse (>2.50 m; < 50 shoots /m²), (c) P. perfoliatus low (< 50 cm max. length), identical to poor growth sites of Chapter 3, (d) Chara sp., and (e) bare sediment.

46 Nutrient experiments

The top layer of sediment (10 cm) was collected by means of hand-driven corers (diameter 12.5 cm). This depth was chosen to include all roots. Overlying water was immediately removed by the principle of communicating tubes. Porewater was extracted from the sediment by centrifugation at 4500 rpm for 10 minutes. Subsamples of the supernatants were filtered (0.45 µm, prewashed Millipore), filtrated and unfiltrated samples were immediately deep frozen and stored until further analysis.

For particle size analysis of the sediments, all visible plant fragments were removed. The samples were wet-sieved through 63 µm and both fractions were dried until constant weight and weighed. The coarse fraction (> 63 µm) was dry sieved through meshes of different sizes according to DIN 18123 and weighed.

Sediment water content was determined by drying 5 ml samples of known weight at 100 °C for 24- 48 h (until constant weight), and weighing of the dried sediment.

Total organic content of the sediment was determined as LOI (loss of ignition) by drying the sample at 100 ° C and determining the weight loss after combustion at 500 °C for 12 h in a muffle furnace. The total carbon and total nitrogen content were determined with a CHN analyser (Hereaeus).

Concentrations of nitrate and SRP in interstitial water were determined using ionchromatography (Dionex). Concentration of total P was determined photometrically by flow-injection analysis (Eppen- dorf) following digestion with K2S2O8.

Proportional data were transformed to the arcsine of the square root before analysis. Statistical analyses were performed with SAS software (SAS Institute Inc., Cary, NC, USA, 1999–2001).

Mesocosm fertilization and tissue analysis

A basin of 5 m width, 10 m length and 1.2 m depth (Fig. 5.1) was filled with a 15 cm layer of lake sediment from Lower Lake Constance. Flow-through of lake water was regulated to about 10-20 l/min, i.e. an approximate retention time of 56 hours. P. perfoliatus shoots of 15-20 cm length were collected near Unteruhldingen on 13.06.2001 and immediately planted at densities of 8 shoots per 0.25 m². The experimental set-up consisted of three treatments and 4 randomly distributed replicates. The test quadrats were left untreated, fertilized with 1 kg m-² slow-release N-P-K-fertilizer Plantacote Depot -1 - 4M® (14% N, 9% P2O5, 15% K2O) (= 1 g P l sediment), or fertilized with 4 kg m ² Plantacote Depot 4M® (= 4 g P l-1 sediment). This is in the range of in-situ fertilization experiments listed in the review by Worms et al. 2000. After 14 weeks, the complete clonal plants were removed from the sediment. Number and biomass of shoots, length, spacer length, biomass and number of branchings of rhizomes, as well as the root biomass were determined. Poor plant growth in the treatment with 4 kg P m-2 was attributed to overfertilization and results are not presented here.

Tissue samples were taken to test for differences in nutrient contents. Three to five healthy looking apical leaves without any sign of senescence were ground in a mortar. C was determined by Autoanalyzer (ThermoQuest NCS-2500) after combustion. P and N values were determined by AutoanalyzerII (Bran&Luebbe) following digestion with potassium persulfate. Proportional data were transformed to the arcsine of the square root before analysis. Statistical analyses were performed with SAS software (SAS Institute Inc., Cary, NC, USA, 1999–2001).

47 Chapter 5

Fig. 5.1 Mesocosm, Limological Institute, University of Constance

Planting experiment

Shoots of P. perfoliatus were collected at the main experimental site near the island of Reichenau (see Chapter 3) on 05.07.2001. On 06.07.2001, the shoots were planted at previously vegetation-free sites of Upper Lake Constance south of the island of Mainau (depth = 4 m), east of Staad (depth = 2.50 m and 4 m) and in the ‘Litoralgarten’ (depth = 4 m) by SCUBA divers. At each location, 6 patches (8 ramets/0.25 m²) were planted, half of them without fertilization, the other half with fertilization of 1 g -1 P l sediment slow release fertilizer Plantacote Depot 4M® (14% N, 9% P2O5, 15% K2O).

After one growth season, on 25.09.2001, ramets were counted and measured, and patch size was determined by SCUBA divers. Since ramet numbers, shoot lengths and patch sizes did not differ between sites under the same treatment, sites were pooled for analysis. Means of the sediment parameters were compared using one-way ANOVA at a level of significance of p < 0.05. One year later, on 06.08.2002, half of the sites were re-visited and re-assessed by Scuba divers according to the same protocol.

Results

Sediment analysis

Grain size, H2O-content and organic C of sediment were not related to vegetation cover type (Tab. 5.1). Furthermore, no significant relationships were found between vegetation cover and porewater nutrients such as SRP, TN and NO3-N (Tab. 5.1). P values ranged from oligotrophic conditions (< 10 µg l−1 SRP (P. perf. sparse, P. perf. poor) to mesotrophic conditions (10 and 30 µg l−1) (P. perf. dense, Chara).

48 Nutrient experiments

Tab. 5.1 Sediment characteristics and nutrient concentrations in porewater of sites with or without vegetation cover

P. perf. P. perf. P. perf. Unit Chara sp. Sediment Pr > F dense sparse poor n 5 3 3 4 5

H2O % 35.64 ± 3.43 34.17 ± 2.03 33.91 ± 4.72 39.10 ± 1.91 31.64 ± 1.61 < 0.43 organic C % 3.96 ± 0.46 3.84 ± 0.34 4.64 ± 0.14 4.20 ± 0.49 3.76 ± 0.07 < 0.52 < 63µm % 44.52 ± 6.02 45.67 ± 12.86 54.70 ± 10.47 35.43 ± 14.01 58.81 ± 4.69 < 0.42 > 63µm % 43.17 ± 5.53 45.83 ± 11.92 16.63 ± 5.62 52.81 ± 16.32 30.94 ± 5.23 < 0.15 > 63mm % 12.31 ± 4.22 8.50 ± 2.21 28.66 ± 16.09 11.79 ± 5.10 10.25 ± 2.15 < 0.30 -1 TP µg l 82.98 ± 42.36 54.03 ± 14.50 122.73 ± 15.16 116.38 ± 46.74 82.84 ± 27.71 < 0.75 -1 SRP µg l 18.46 ± 6.77 6.57 ± 1.83 9.23 ± 1.60 20.08 ± 8.23 11.20 ± 3.37 < 0.45 -1 TN mg l 2.58 ± 0.90 2.23 ± 0.69 2.16 ± 0.10 2.87 ± 0.64 2.65 ± 0.41 < 0.95 -1 NO3-N µg l 3.60 ± 1.86 28.00 ± 23.59 95.67 ± 76.21 722.50 ± 722.5 18.80 ± 4.91 < 0.48 N:P 43.91 39.68 18.10 33.82 42.24

Mesocosm fertilization and tissue analysis

In comparison to controls, nutrient addition did not produce significantly different growth results concerning the number of ramets (ANOVA, F1,7 = 2.39, p < 0.1734), average shoot length (ANOVA,

F1,7 = 0.05, p < 0.8300) and total biomass (ANOVA, F1,7 = 2.29, p < 0.1809) (Fig. 5.2).

n shoot length (cm) total biomass (g DW) 180 45 40 160 30 140 30 120 100 25 20 80 60 20 40 10 20 0 15 0 no nutrients nutrients no nutrients nutrients no nutrients nutrients

Fig. 5.2 Ramet number, shoot length, and total biomass of P. perfoliatus in unfertilized and fertilized quadrats in a mesocosm ± S.E.

Likewise, plant tissues on fertilized and unfertilized sites showed no significant differences with regard to % N (ANOVA, p < 0.81) and % C (p < 0.62) (Tab. 5.2).

Tab. 5.2 Tissue nutrient contents at non-fertilized and fertilized sites

n % N % C C:N no fertilizer 19 4.79 ± 0.13 40.01 ± 0.20 8.44 ± 0.19 fertilizer 15 4.74 ± 0.19 40.22 ± 0.42 8.65 ± 0.29

49 Chapter 5

However, root allocation was significantly lower in the fertilized treatment (ANOVA, F1,7 = 15.64, p <

0.008) while biomass allocation into aboveground (F1,7 = 1.08, p < 0.34) and belowground biomass

(F1,7 = 1.08, p < 0.34) was not significantly different between treatments (Fig. 5.3). Rhizome spacer length was significantly shorter in the nutrient treatment (ANOVA, F1,7 = 16.38, p < 0.056). The spacer length differed not only significantly between but also within treatments (Fig. 5.3). Part of this variation can be attributed to the order of shoots (Fig. 5.4). Rhizome length was shortest between primary and secondary ramets (4 cm) and increased in the direction of new ramet recruitment. The average branch spacer was longer than the average main axis spacer (i.e. the axis originating from the primary ramet) and increased the later the branch was initiated.

shoots (%) rhizomes (%) root biomass (%) spacer length (cm) 88 17 4 10 87 16 9 86 15 3 85 14 8 84 13 83 7 2 12 82 6 81 11 80 10 1 5 no nutrients nutrients no nutrients nutrients no nutrients nutrients

Fig. 5.3 Biomass allocation and spacer length of P. perfoliatus ± S.E. without and with fertilization

spacer length (cm) spacer length (cm) 15 20

15 10

10

5 5

0 0 1 2 3 4 5 6 7 m b1 b2 b3 b4 spacer number main axis (m) / branch number (b) Fig. 5.4 Average rhizome spacer length of P. perfoliatus in relation to position within rhizome ± S.E. (a) rhizome spacer length of main axis, 1 is the oldest, 2 is the youngest spacer (b) spacer length of main axis and branches, b1 is the oldest branch, b4 is the youngest branch

50 Nutrient experiments

Planting experiment

After termination of the experiment, ramet number and ramet length were significantly larger in the fertilized treatments than in the control (t-test, p < 0.001, Fig. 5.5). At unfertilized sites, there was little growth at all. Moreover, on 33 % of unfertilized sites all plants had died.

In 2002, when sites were re-visited, patch extension on fertilized sites had increased at an average of 6.5-fold, whereas on unfertilized sites, patch area had either decreased, or patches had completely disappeared.

ramet number length (cm) 70 70

60 60

50 50

40 40

30 30

20 20

10 10

0 0 control fertilizer control fertilizer

Fig. 5.5 Ramet number, ramet length and patch size of P. perfoliatus planted on unfertilized and fertilized sediment ± S.E.

Discussion

Sediment analysis

Although previous experiments have clearly shown that growth performance and morphology of submersed macrophytes may significantly be related to sediment characteristics such as grain size, density, organic matter, and nutrients (Anderson & Kalff 1986, Barko & Smart 1986a, Walser 1996), we found no differences of sediment composition and nutrient content at sites with different or no vegetation growth in Lake Constance. However, also other authors have failed to relate vegetation distribution and growth to sediment parameters in-situ (Barko & Smart 1979, Anderson & Kalff 1986, Lehmann 1997).

We suggest that the lack of correlation could be due to sampling problems: Reproducibility depends on both the precision of the technique and the spatial heterogeneity of the sediments. However, although previous papers have presented replicate porewater profiles (Carignan 1984, Carignan et al. 1994), the reproducibility has seldom been quantitatively examined (Urban et al. 1997). One reason for our insignificant results was the high variability of values leading to a high standard error. In fact, the average SRP value at poor growth sites was only 50% of the value at rich growth sites and more replicates could possibly have led to a significant result. Furthermore, nutrients were measured based on volume in porewater, Barko & Smart (1986a) report, however, that growth was rather correlated with sediment nutrient concentrations based on sediment volume.

In addition, there is a lack of knowledge which nutrient fraction is most representative for macrophyte growth. For example, N is more easily taken up as NH4, but more often taken up as NO3

(Agami & Waisel 1986). Measurements reported in literature include total P, SRP, RP, PO4-P, o-P,

51 Chapter 5

total N, NO3-N, NH4-N, NH3-N etc. Moreover, these parameters may be assessed either in the sediment or in the porewater. This is problematic because P in the porewater constitutes only a very small fraction of the total P bound in the sediment (Søndergaard 2007), but there is hardly any information about bioavailability of sediment P (Golterman 2004). While most authors found relationships between porewater nutrient concentration and nutrient status of plant tissue (Fourqurean 1992, Agawin et al. 1996, Alcoverro et al. 1997), the results are more inconsistent for relationships between porewater nutrients and biomass, plant growth or vegetation cover: There are reports of positive correlations between biomass and PO3 (Chambers & Fourqurean 1991), NH4 (Chambers & Fourqurean 1991, Cízková 2001), SRP (Fourqurean 1992, Cízková 2001), P, TP and N, (Cízková 2001). In other cases, no relations were found between biomass and P (Chambers & Fourqurean

1991), NH4 (Fourqurean 1992), TN (Chambers & Fourqurean 1991, Cízková 2001), or total dissolved solids (Cízková 2001). Similarily, only poor relations were found between shoot nutrient content and growth (Barko & Smart 1986a).

Furthermore, it has to be taken into account that N and P in the interstitial water are characterized by a high spatial and temporal variability (Carignan 1985) and a pronounced seasonal cycle with a minimum in summer (Madsen 1988, McComb et al. 1998). The same applies to tissue nutrient concentrations (Nichols & Keeney 1976). Therefore, it is likely that more replicates have to be taken throughout the growth season to provide meaningful results.

Mesocosm fertilization and planting experiment

Production parameters and tissue analysis

In contrary to the expectation that fertilization will increase growth, we found no significant effects of fertilization on biomass in the mesocosm. Nevertheless, tissue nutrient levels were rather high - %N was 92 % higher than in-situ and C:N was 42 % lower than in-situ (Wolfer & Straile 2004b, Chapter 3) - indicating that nitrogen tissue content is not always correlated to growth in macrophytes as proposed by Ågren G.I. (1985). Similarly, there are reports about insignificant effects of fertilization of seagrass meadows, even though raised tissue contents indicate increased uptake of nutrients (Erftemejier et al. 1994, Worms et al. 2000). In contrary, results presented in Chapter 3 show that in-situ fertilization of P. perfoliatus plants resulted in increased shoot numbers, longer shoots, and higher biomass on poor growth sites, and planted ramets grew better when nutrients were added. Also Alcoverro et al. (1997) found different responses at different sites to fertilization. These results may be due to the naturally different sediments. The fact that fertilization showed no increase of growth parameters with regard to ramet and biomass production and the high tissue nutrient concentrations in the mesocosm indicated that the sediment utilized was already rich in nutrients. Since productive sites are dominant in the Lower Lake Constance, we assumed that the poor growth sites were a product of either sediment disturbance or nutrient depletion by previous plant uptake. On the other hand, the establishment of healthy stands of P. perfoliatus at fertilized sites at Upper Lake Constance shows that the lack of macrophytes in extensive littoral areas of the Upper Lake Constance might be due to N or P nutrient limitation, whereas sediment grain size, currents, herbivores, other nutrients etc. play a minor role.

Root allocation, rhizome / shoot allocation and spacer length

Lower root allocation in fertilized treatment or nutrient rich sites as found in our study is commonly reported (Barko & Smart 1981b, Ervin & Wetzel 1997, Wertz & Weisner 1997, Peralta 2003) and has been experimentally proven to maximize plant growth (Hirose 1987, Mooney et al. 1988). Resource

52 Nutrient experiments

allocation requires balancing of the competing requirements of photosynthetically productive tissue and the nutrient uptake apparatus (Madsen 1991).

On the other hand, biomass allocation into aboveground and belowground biomass was not significantly different between treatments. This is in contrast to other studies (review of Hutchings 1997) where rhizome allocation decreased with increasing nutrients. The reason is probably due to the fact that in our study, biomass of fertilized and unfertilized sites was equal, whereas in other studies, biomass increased due to fertilization.

As predicted by the foraging theory, rhizome spacer length was significantly shorter in the nutrient treatment. Also in-situ (Chapter 3), mean rhizome spacer length was shorter in the fertilized treatment than in the control. According to the foraging theory the benefit of shorter rhizomes consists of a concentration of shoots in favourable places, whereas longer spacers ‘move’ the ramets away from unfavourable sites (for a more detailed discussion of foraging see Chapter 8). This seems to contradict the recent rejection of the foraging theory by many scientists including the original authors (De Kroon & Hutchings 1995), mainly because in rhizomatous species effects of higher nutrients on rhizome spacer length are less consistent (Granéli et al. 1992, Grace 1993, Cain 1994, Dong & De Kroon 1994, Hutchings & De Kroon 1994) than in stoloniferous species. Our findings are also in contrast to the growth hypothesis which suggests that rhizomatous species increase rather than decrease the length of their rhizomes with increasing resource-availability (De Kroon & Knops 1990, Hutchings & Mogie 1990, Hutchings & De Kroon 1994, Slade & Hutchings 1987, Stoll et al. 1998). The fact that rhizome spacer length decreased significantly under fertilization, whereas other growth parameters including rhizome allocation remained the same shows that in this case rhizomes became thinner and longer and that the change in spacer length was an independent part of the foraging behaviour. As another example of the influence of nutrients on the clonal architecture of P. perfoliatus, we could also confirm increased branching at higher nutrient availability in-situ (Chapter 3).

Summarizing,

(a) we found no correlation between sediment characteristics and vegetation cover. On the other hand, in the planting experiment only fertilized ramets could establish, indicating that the plants are nutrient limited.

(b) in contrary to in-situ measurements, we found no correlation between tissue nutrient concentration and fertilization in the mesocosm; this is most likely due to the fact that the original nutrient level was high and plants were not nutrient limited in the mesocosm.

(c) we found shorter rhizomes and increased branching at fertilized sites which is in line with the foraging theory

The effect of nutrients on P. perfoliatus has implications on its distribution in Lake Constance: We could confirm that P. perfoliatus occurs in meso- and eutrophic water bodies and has a high tolerance against raised phosphate and nitrate concentrations (Pietsch 1982) as fertilization increased growth in- situ, and levels of fertilization in the mesocosm were rather high. Therefore, we assume that the decline of P. perfoliatus during eutrophication was partly due to secondary effects of high nutrient loads, such as shading by 5 to 10-fold higher phytoplankton production (Müller 1977) and periphyton growth.

53 Chapter 5

54 Clonal integration

Chapter 6. TO SHARE OR NOT TO SHARE – CLONAL INTEGRATION IN RAMETS OF POTAMOGETON PERFOLIATUS

Abstract

The ability of clonal plant species to share resources has been studied in many experiments. Potamogeton produces interconnected ramets within short (time) intervals and hence may or may not share resources with ramets growing in less favourable microhabitats.

From a genet point of view sharing with ramets growing under less favourable conditions might not be an optimal strategy when photosynthates could be used to establish other ramets growing under more favourable conditions.

In order to analyse the plasticity in clonal integration of a submerged macrophyte (Potamogeton perfoliatus), we set up a factorial aquaria experiment with unshaded or shaded recipient ramets (offspring) which were connected to or separated from donor ramets (parents).

Increased biomass production of offspring in parent-offspring-systems compared to severed offspring in both light and shade showed that ramets share resources via clonal integration.

The relative translocation to the first and second offspring generation was influenced by habitat quality: if first offspring ramets (i.e. the ramets that develop after the mother ramet) were shaded by nets, second offspring ramets clearly profited. This may be partly due to the fact that resources are shifted from first offspring to second offspring ramets.

This complex sharing behaviour might be relevant when plants produce ramets within a dense patch of macrophytes, where support of a shaded ramet might not pay off.

Introduction

Clonal plants are characterized by the reiteration of potentially independent modules, called ramets which consist of shoots, rhizomes or stolons, and roots. Clonal integration involves resource sharing through rhizomes and plays an important role in the regulation of shoot growth. The transport of water, nutrients and photosynthetates has been shown to increase the capacity of plants to tolerate resource heterogeneity, to colonize different microhabitats, and to recover from herbivory (Ong & Marshall 1979, Schmid et al. 1988, Alpert 1999). The degree of resource sharing of a clonal plant species is under both genetical and environmental control (Alpert 1999, Van Kleunen et al. 2000). The primary motor behind clonal integration might either be the resource export from parents acting as a source (push model), or the demand of offspring acting as a sink (pull model) (Pitelka & Ashmun 1985, Marshall & Price 1997). Hormones such as auxins, or cytokinins have been proposed as regulating mechanisms (Alpert et al. 2002).

Among aquatic plants, clonal integration is known from emergent macrophytes (Hester et al. 1994, Amsberry et al. 2000), floating stoloniferous species (Alpert et al. 1991), submerged macrophytes (Xiao et al. 2007a) and marine seagrasses (Tomasko & Dawes 1989, Marbà et al. 2002). We therefore expected that the submerged rhizomatous freshwater species Potamogeton perfoliatus L. is also capable of clonal integration. In-situ surveys in Lake Constance revealed that each P. perfoliatus plant sprouts from a turion in spring, and produces horizontal rhizomes of up to 1.5 m length, bearing up to

55 Chapter 6

15 and more ramets during one season (Wolfer & Straile 2004a,b). Whereas existing experiments studied overall transfer within a multi-ramet clone or sharing within a 2-ramet system (parent and offspring), we investigated the trade-off between the parental support of two ramet generations growing under different habitat conditions. Optimality modelling suggests that in this case, biomass allocation depends strongly on microhabitat heterogeneity (Gardner & Mangel 1999).

In a 2 x 2 factorial aquarium experiment, we tested for the effects of shading of one offspring on the growth of two offspring generations. The first offspring was either left connected to or severed from the parental ramet. Hence, we analysed the relative performance of two ramets in a 3-ramet system (connected) and in a 2-ramet system (severed). At the start of the experiment, the 2nd of the offspring ramets was not yet established. Therefore, in the 2-ramet system where no support from the parent is possible, the success of the clonal fragment will depend on the photosynthesis of the shaded ramet. In contrast, in the three-ramet system, the shaded ramet is not crucial for the success of the clonal fragment as the parent ramet can share resources with the unshaded offspring. Hence, we tested the following hypotheses:

(1) growth of offspring ramets in the 3-ramet (parent-offspring) system is enhanced as compared to growth of offspring ramets in the 2-ramet system (severed offspring system), i.e. there is clonal integration between parent and offspring.

(2) growth differences between the two offspring ramets will be different in the 2- versus the 3-ramet system.

(3) shading influences the growth differences of the two offspring ramets, i.e., the relative translocation of resources to different offspring generations depends on the habitat quality of the individual ramets.

Materials and methods

Origin and pre-cultivation of plant material

P. perfoliatus shoots were collected at the Lower Basin of Lake Constance, a large, meso- oligotrophic lake in central Europe (9°18'E, 47°39'N). All shoots originated from the same patch with a diameter of approximately 15 m and had an intermediate developmental age. Since they were cut off above the sediment, they had no rhizomes and no roots when planted. The shoots were planted across three aquaria (length: 80 cm, width: 40 cm, height: 50 cm), filled with 10 cm of natural sediment from Lake Constance and with 125 l of filtered lake water. The water was exchanged twice a week; light was provided by pairs of white and plant-grow tubes at 14 hours a day.

Experimental design

The experiment was started when all planted shoots (parent ramets) were well established and had produced one offspring ramet (O1), i.e. the first offspring ramet. 50 % of the offspring ramets were left connected to their parent, 50 % were severed from their parent. In addition, habitat heterogeneity was introduced in 50 % of the connected and severed observation units by shading O1 with cylinder- shaped nets (Agroflor, height: 40 cm, diameter: 6 cm, shading effect: 63 %) (Fig. 6.1). Eight replicates of each treatment were distributed randomly across the three aquaria. At the beginning of the experiment, neither shoot lengths of P ramets (31 cm ± 8 cm), nor shoot lengths of O1 ramets (12 cm ± 5 cm) differed between treatments.

56 Clonal integration

a

b PO1O2

SO1

Fig. 6.1 (a) Aquarium (b) sketch of an experimental unit: P = parent ramet, O1 = first offspring ramet, O2 = second offspring ramet, SO1 = spacer length between offspring ramets. Rhizomes between P and O1 were either left connected or severed midway between the shoots. O1 ramets were either kept unshaded or shaded

The experiment was ended after three weeks, when plants had produced the second offspring ramet (O2). The plants were carefully removed from the sediment with their rhizomes and roots, and thoroughly washed. After the measurement of shoot lengths and lengths of the rhizomes between O1 and O2 - called spacer length (SO1) below - shoot, rhizome and root fractions were dried at 105°C, cooled down and weighed on an analytic scale.

Resource sharing was inferred from increased shoot length and biomass in connected, compared to severed offspring ramets (Tomasko & Dawes 1989, Amsberry et al. 2000), but results are only presented for biomass since analyses based on shoot measurements yielded consistent results. Effects of treatments and their interactions on plant growth were tested by ANOVA. As we were not able to use one aquarium for every replicate, we included the factor ’aquarium’ in the ANOVA to correct for overall growth differences between the three aquaria. Proportional data were transformed to the arcsine of the square root before analysis. Statistical analyses were performed with SAS software (SAS Institute Inc., Cary, NC, USA 1999-2001).

Results

Clonal integration

After three weeks of growth, severing had reduced total shoot biomass of offspring ramets (O1 +

O2) strongly (63 % reduction, ANOVA; F1,28 = 25.9, p < 0.0001, Fig. 6.2), whereas shading of the first offspring (O1) had a smaller effect (30 % reduction, F1,28 = 5.5, p < 0.05). Shading did not influence the response of offspring ramet shoot biomass to being severed (severing x shading interaction: F1,28 = 0.02, p < 0.8).

57 Chapter 6

2.40 ) ) 2.15 ght i g dw we

m 1.90 ( 10 g o

Shoot 1.65 (l

1.40 3-ramet 2-ramet system system Fig. 6.2 Total shoot biomass of both offspring ramets (mean ± S.E) of P. perfoliatus in the 2-ramet versus 3- ramet system. Solid lines indicate ‘light’, hatched lines indicate ‘shade’ treatments

Offspring ramets were differently affected by severing and by shading (Fig. 6.3). Biomass of O1 was highest in the 3-ramet - light treatment as compared to all other treatments. When shaded, there was no significant biomass difference between the 2- and 3-ramet systems, meaning that O1 performance was not significantly increased by the presence of a parent plant. In contrast, O2 did profit from the 3- ramet system under both light and shade conditions (of O1). Furthermore, its biomass increase due to clonal integration was higher when O1 was shaded than when O1 was unshaded. Shading of O1 did not influence the biomass of O2 in the 3-ramet system. As a consequence, biomass of O2 was similar to the biomass of O1 in the 3-ramet - shade treatment, but was less than O1 in all other treatments (Fig. 6.3).

) AAAB 100 dw g (m

10 weight

Shoot 1 abcc 2 2 3 3 S L S L Fig. 6.3 Shoot biomass of the offspring ramets O1 (filled circles) and O2 (open circles) in the different treatments. 2 and 3 refer to the 2-, respectively 3-ramet systems, and S (shade) and L (light) refer to the two light conditions. There are significant treatment differences for O1 biomass (ANOVA: F3,28 = 9.41, p < 0.0002) and for O2 biomass (ANOVA: F3,28 = 15.2, p < 0.0001). Capital letters at the top of the graph indicate treatment differences in O2 biomass, lower case letters at the bottom of the graph show treatment differences in O1 biomass. Different letters indicate significant differences between treatments (Duncan’s post-hoc test, p < 0.05)

This complex response pattern is supported when analysing the full data set with an ANOVA and distinguishing the two offspring ramets with the factor ’ramet order’. Both O1 and O2 were negatively influenced by severing and shading (significant main effects: severing, shading, and ramet order, Tab. 6.1). Severing had a stronger effect on O2 biomass as compared to O1 biomass (significant severing x ramet order interaction, Tab. 6.1), whereas there was no indication that shading of O1 affected O1 and

58 Clonal integration

O2 differently (non-significant shading x ramet order interaction, Tab. 6.1). Finally, offspring ramets differed in their response to shading in the 2- versus 3-ramet system (significant three way interaction: severing * shading * ramet order, Tab. 6.1). There was also a small aquaria effect, but the inclusion or omission of the factor ’Aquarium’ in the ANOVA neither affected the significance levels nor the effect sizes of the studied main effects or interactions.

Tab. 6.1 ANOVA with shoot biomass (log transformed) as dependent variable and severing, shaded ramet order (first or second offspring) and aquarium as independent variables. The latter was included to account for effects due to different growth conditions in the three aquaria. Model R2 = 0.74, n = 58, F = 15.0, p <0.0001.

DF Type III SS F Value Pr > F Severing 1 5.55 56.53 <.0001 Shading 1 1.28 13.09 0.0007 Severing * shading 1 0.08 0.83 0.3672 Ramet order 1 2.91 29.60 <.0001 Severing * ramet order 1 1.47 14.94 0.0003 Shading * ramet order 1 0.06 0.66 0.4201 Severing * shading * ramet order 1 0.63 6.39 0.0148 Aquarium 2 0.77 3.92 0.0264

As a consequence of different responses of O1 and O2 to severing and shading, the growth of O2 relatively to O1 (O2 / O1) was significantly influenced by severing (F1,28 = 22.82, p < 0.0001, Fig. 6.4) and the interaction between severing and shading (F1,28 = 8.73, p < 0.01): Shading of O1 increased the O2 / O1 biomass ratio when O1 was connected to a parent ramet, but decreased it when the connection was severed. Biomass of parent ramets was neither affected by severing nor by shading

(F1,28 = 0.33, p < 0.6, and F1,28 = 0.1, p < 0.75).

2.000 1 /O

2 1.000 O

o ti

a 0.500 s r s a 0.250 m o i B 0.125 3-ramet 2-ramet system system Fig. 6.4 Biomass ratio of second offspring and first offspring (O2 / O1) (mean ± S.E.) of P. perfoliatus grown in aquaria. Solid lines indicate ‘light’, hatched lines indicate ‘shade’ treatments

59 Chapter 6

Biomass allocation and spacer length

Root allocation of O1 was affected by severing (F1,28 = 3.53, p < 0.08), shading (F1,28 = 3.75, p <

0.07) as well as the interaction between the two (F1,28 = 7.93, p < 0.009): Shading increased root allocation of connected O1 but reduced it in severed O1 (Fig. 6.5a). No significant treatment effects were observed for shoot allocation (Fig. 6.5b), and rhizome allocation (Fig. 5c). Shoot biomass was significantly related to shoot (r = 0.36, p < 0.05) and rhizome (r = -0.47, p < 0.006) allocation but not to root allocation (r = 0.03, p > 0.5). Consequently, treatment effects on root allocation were also significant when including shoot biomass as a covariate into the statistical model (shading: F1,27 =

4.67, p < 0.05, severing: F1,27 = 4.45, p < 0.05, shading x severing: F1,27 = 6.00, p < 0.03, shoot biomass: F1,27 = 0.93, p < 0.4).

50 100 ) 50 ) )

a b % c ( % %

( (

40 90 n 40 o i on t on ti ti a a c a 30 80 30 c c o l o o l l l l

20 a 70 20 e al a t t o oo om o

10 60 z 10 R i Sh

0 50 Rh 0 3-ramet 2-ramet 3-ramet 2-ramet 3-ramet 2-ramet system system system system system system

Fig. 6.5 (a) Shoot (b) rhizome and (c) root allocation ± S.E. of 3-ramet system and 2-ramet system in P. perfoliatus first offspring (O1). Solid lines indicate ‘light’, hatched lines indicate ‘shade’ treatments

Spacer length between offspring ramets (SO1) was reduced in the 2-ramet system (F1,28 = 6.09, p <

0.03), whereas shading had no significant effect on SO1 (F1,28 = 0.17, p < 0.7). There was also no influence of shading on the response of SO1 to severing (severing x light interaction F1,28 = 2.13, p < 0.16) (Fig. 6.6).

)

m 8 (c h t

ng 6 le er

c 4 a p s

e 2 om iz

h 0 R 0.0 0.1 0.2

Shoot biomass (mg dw)

Fig. 6.6 Spacer length between offspring ramets (SO1) in relation to first offspring (O1) shoot biomass. Filled symbols indicate 3-ramet system treatments, open symbols indicate 2-ramet system treatments

60 Clonal integration

SO1 was positively related to O1 shoot biomass (r = 0.57, p < 0.007). However, if O1 in the 2-ramet systems was analysed separately, SO1 was negatively related to shoot biomass (r = -0.51, p < 0.05), whereas in 3-ramet systems, SO1 was positively related to shoot biomass (r = 0.71, p < 0.005) (Fig. 6.6). This suggest that the 3-ramet system O1 responded to growth conditions - which were at least partially related to shading - by decreasing rhizome length, whereas the 2-ramet system O1 responded to growth conditions by increasing rhizome length.

Discussion

Clonal integration

The substantially higher biomass of offspring in the 3-ramet systems as compared to the 2-ramet systems shows that parent ramets of P. perfoliatus were capable of sharing resources acropetally via clonal integration. However, clonal integration was not uniform, but rather sensitive to the growth conditions of the first offspring ramet: O1 ramets benefited substantially from integration only when in the light but were not supported when shaded. In contrast, many previous studies have found physiological integration especially where recipient ramets experienced stress through resource limitation in heterogeneous environments (Alpert & Stuefer 1997, Hutchings & Wijesinghe 1997, Alpert 1999). For example, connected shaded shoots of the seagrass Thalassia testudinum achieved growth rates equal to non-shaded controls (Tomasko & Dawes 1989). Likewise, severely shaded ramets of Lolium perenne continued to grow and produced new leaves as a result of physiological integration (Ong & Marshall 1979), and shaded shoots of Eichhornia crassipes suffered less when connected to parent shoots growing in the light (Methy et al. 1990). However, in our study the proximal offspring (O1) benefited from the 3-ramet system, i.e., from the connection to its parental ramet only under full light, but was not supported when shaded. In the latter situation, the shaded ramet O1 was passed by and resources were directed to the distal offspring O2. Indeed, O2 ramets of P. perfoliatus benefited from resource sharing relatively more strongly when O1 ramets were shaded. This shows that clonal integration in P. perfoliatus can also occur between non-adjacent shoots (see also (Terrados et al. 1997) and that resources are shared with distal offspring depending on the microhabitat, i.e. growth conditions, of the proximate offspring. Likewise, resource sharing has often been observed with ramets suffering from herbivory (Marshall & Sagar 1965, Schmid et al. 1988, Olson & Wallander 1999). However, lack of support for damaged ramets of a perennial herb was recently demonstrated and attributed to competition between sibling ramets (Hellström et al. 2006). Competition within branches of a plant has also been demonstrated in pine trees (Honkanen & Haukioja 1994) and pea plants (Novoplansky et al. 1989) and has been termed the branch- competition hypothesis (Sachs & Novoplansky 1997). This hypothesis predicts that a plant module which is inferior due to e.g., herbivory damage or microhabitat unsuitability should be left out of support when more viable sinks are available. Our results are in line with this prediction. However, shaded O1 in the 3-ramet system showed even strong signs of chlorosis which was not the result of shading per se, as shaded O1 in the 2-ramet system appeared vigorous and retained green leaves. This observation suggests that the branch- ompetition hypothesis is unlikely to give a complete explanation of our results: At the end of the experiment, shoot biomass of O1 and O2 were similar (Fig. 6.3). Assuming a faster growth rate of unshaded O2 relative to shaded O1 suggests that during most of the experimental time, biomass of O1 was larger than biomass of O2. It is hence difficult to believe that despite these biomass differences, competitive superiority of O2 was large enough to result into chlorosis of O1. Rather, the large differences in performance of O1 suggest controlled senescence of shaded O1 in the 3-ramet

61 Chapter 6

treatment possibly associated with remobilization of resources (Ong & Marshall 1979, Stapel & Hemminga 1997) towards O2. Hence, integration seems to have qualitatively altered the response of shaded O1 by inducing a novel response, i.e. chlorosis (De Kroon et al. 2005).

Biomass allocation and spacer length

Shading increased root allocation of connected offspring but reduced root allocation in severed offspring. Our statistical analyses suggest that the differences in biomass allocation are not due to allometric growth rules which predict that smaller plants show a higher biomass allocation to belowground structures (Müller et al. 2000). Despite their small biomass, O1 in 2-ramet - shade systems had the lowest root allocation. This might result from a shortage of carbohydrates and the need to invest in shoot biomass (Alcoverro et al. 1997). The lower root allocation in the severed shaded offspring compared to higher root allocation in connected shaded offspring is also in line with the foraging hypothesis proposing that single plants specialize in the most limiting resource (here: light), and integrated plants specialize in the most abundant resource (here probably nutrients) (Stuefer et al. 1996, Stuefer 1998). In a clonal ramet system with ‘division of labour’ (Hutchings & Wijesinghe 1997) ramets may continue to take-up nutrients by roots even when they are non- photosynthetic (Jónsdottir & Callaghan 1990) and in our experiment, shaded O1 may still contribute to plant growth by supplying nutrients to unshaded O2.

As in some terrestrial plants (Wijesinghe & Handel 1994, Van Kleunen et al. 2000), severing significantly reduced spacer lengths of P. perfoliatus offspring. This might be attributed to the stronger effect of severing on offspring shoot biomass, and to the overall positive relationship between shoot biomass and rhizome length. In-situ, spacer length of P. perfoliatus strongly increases with distance from the primary ramet, possibly as a consequence of an increasing biomass and production of an interconnected clonal fragment with the number of produced ramets (Wolfer & Straile 2004b). As a result of shading, O1 in the 3-ramet system responded to the unfavourable growth conditions by decreasing rhizome length, whereas O1 the 2-ramet system responded to shade conditions by increasing rhizome length, even though shoot biomass was slightly decreased. This is in line with the predictions of the foraging hypothesis: shaded shoots are expected to produce longer rhizomes in order to 'escape' from the unfavourable habitat (Hartnett & Bazzaz 1983, Sutherland & Stillman 1988). Possibly, foraging is only expressed in the growth patterns when growth is strongly impaired by shading, and integration is not possible (but see De Kroon & Hutchings (1995) for a critical discussion of the plant foraging hypothesis).

To conclude, our experiment has shown that (1) P. perfoliatus parent ramets support their clonal offspring ramets via translocation of resources, i.e. there is clonal integration within a genet, and (2) the relative translocation of resources to different offspring generations depends on the habitat quality of the individual ramets: ramets in unfavourable microhabitats are not integrated when support of other ramets provides higher benefits for the genet. This behaviour might be highly relevant when plants produce new ramets within dense patches of macrophytes (Wolfer & Straile 2004a, Wolfer & Straile 2004b). In such a case it might not benefit the plant to support a severely shaded ramet, but rather to invest in rhizome growth and new ramets at the outer perimeter of the patch, where microhabitats are more suitable.

62 Clonal integration

63 Chapter 6

64 Modelling

Chapter 7. MODELLING THE CLONAL GROWTH OF THE RHIZOMATOUS MACROPHYTE POTAMOGETON PERFOLIATUS

Susanne R. Wolfer, Egbert H. van Nes & Dietmar Straile

Ecological Modelling 192, 67-82 (2006)

Abstract

Macrophytes play a crucial role in the functioning of lake ecosystems. Until now most macrophyte models neglected the fact that the majority of macrophyte species expand clonally during the growing season. Inclusion of a detailed description of clonal growth in models can facilitate our understanding of space occupation and patch expansion and predict future macrophyte development. ‘CLOMO’ is an individual-based model which includes a detailed, spatially explicit description of rhizome formation and clone expansion as well as a realistic description of photosynthesis including light limitation and temperature. The model also accounts for transfers of energy or resources between different parts of the clone (‘clonal integration’).

Although the clonal growth of macrophytes is complex and poorly known, the first model results for the macrophyte species Potamogeton perfoliatus were promising and compared well with the field data. The model can produce growth networks very similar to those found in the field. A Monte-Carlo sensitivity analysis showed systematically which parameters have the largest effect on the architecture and expansion of the clones.

The application of the model provided new insights into growth dynamics and patch development:

(1) The model showed that a lack of branching will lead to the extinction of the clonal fragment after a certain number of years. This is due to the fact that the reproduction organs (turions) are formed at the end of a branch and even a small turion mortality will cause a reduction in surviving plant numbers.

(2) The growth of rhizome axes relative to those in the previous year determines the patch density and patch expansion rate. Reversing rhizomes lead to compact patch growth whereas continuing rhizomes lead to loose aggregates.

65 Chapter 7

Introduction

The majority of submersed macrophytes form clones which are connected by belowground rhizomes. Observed from above, a macrophyte patch appears like a simple collection of shoots, but a closer look into the sediment reveals a complex network of ‘ramets’ (= potentially independent units with leafy shoots and roots) interconnected by ‘rhizomes’ (= horizontal shoots lacking chlorophyll). One plant of P. perfoliatus can, for example, consist of more than 15 ramets connected by more than 1 m of rhizome (Wolfer & Straile 2004a). Experiments have shown that ramets exchange energy or resources through these rhizomes (Hartnett & Bazzaz 1983), a process called ‘clonal integration’. Furthermore, rhizomes are mobility units, as a clone can move slowly by creating new rhizomes. The architectural growth of main rhizome axes, rhizome branchings, rhizome angles and rhizome ‘spacers’ (= rhizome connection between two neighbouring shoots) follows species-specific clonal rules (Callaghan et al. 1990, Evans & Cain 1995). Nonetheless, spacer lengths and branch angles also show considerable variation (Cain et al. 1995, Wolfer & Straile 2004a,b), sometimes in response to their growth conditions (De Kroon & Hutchings 1995), aiming at the effective exploitation of local resources such as light and nutrients (Callaghan et al. 1990).

Clonal growth architecture is important for macrophyte fitness because it determines propagation, growth and survival, space occupation and patch expansion. Also the sharing of resources such as carbohydrates or nutrients increases fitness (Wijesinghe & Handel 1994, Stuefer et al. 1996, Hutchings & Wijesinghe 1997).

Our objective was to develop a model which can improve our understanding of mechanisms that could cause the observed field patterns of clonal architecture and patch expansion. It does not aim at replacing existing models of aquatic macrophytes, but rather adds a supplementary application. On the long run, it could be used to predict future macrophyte development and the impact of environmental variation.

Although there are many detailed dynamical models of macrophytes (Collins & Wlosinski 1989, Scheffer et al. 1993, Hootsmans 1994, Chen & Coughenour 1996, Muhammetoglu & Soyupak 2000, Van Nes et al. 2003), the clonal growth of submerged macrophytes has, to our knowledge, never been modelled in detail, including architectural, spatial, environmental, temporal and demographic components. For terrestrial plants, only a few of such detailed models exists (Herben & Suzuki 2002, Oborny & Kun 2002). The facts that photosynthesis in water has fundamentally other restrictions than in air and that light submission in water is far more complex require separate models for submersed species.

Here we present the individual-based model ‘CLOMO’ (‘clonal module’), a spatially explicit model that combines a stochastic description of clonal expansion with a ramet-based calculation of primary production and respiration using physiological rules of the non-clonal model Charisma (Van Nes et al. 2002, 2003). Clonal integration was described as simple as possible, since transport between ramets is a poorly understood process. We applied the model to describe the spatial architecture of P. perfoliatus L. in Lake Constance.

Special analyses were performed with regard to branching frequency and rhizome growth direction. Furthermore, a Monte Carlo sensitivity analysis (Klepper et al. 1994) was applied to show systematically which parameters have the largest effect on the architecture and expansion of the clones.

66 Modelling

Methods

Model description

Overview

CLOMO is an extension of the model Charisma (Van Nes et al. 2002, 2003). Although CLOMO uses the biomass growth unit of Charisma, its new components make it fundamentally different from the earlier version (see Tab. 7.1). While the biomass model Charisma is non-clonal and uses super- individuals to cope with large numbers of individuals (Scheffer et al. 1995), the current model is truly individual-based and operates on a smaller spatial scale. Clones are modelled as a set of ramets (nodes with shoots) interconnected by rhizomes. The model is explicitly spatial, describing patches of plants growing on a grid.

Tab. 7.1 An overview of the main differences between non-clonal macrophyte models (e.g. Collins & Wlosinski 1989, Scheffer et al. 1993, Hootsmans 1994, Chen & Coughenour 1996, Muhammetoglu & Soyupak 2000, Van Nes et al. 2003) and the clonal growth model CLOMO

Non-clonal macrophyte models Clonal macrophyte model

Shoots are the basic unit for photosynthesis, (self-) Ramets are the basic unit for photosynthesis shading is calculated for each shoot

From each seed or tuber emerges a maximum of From each tuber emerges one ramet in spring; one shoot in spring; no further shoot emergence additional ramets are formed during the growth during the growth season season if there is enough biomass accumulated

Individual shoots, no clones Several ramets form a clone

(Almost) all shoots are assumed to be of equal age Each ramet has a different age and length and length

Spacer lengths and number of ramets per clone are Spacer lengths and number of ramets per clone are ignored important components

No patch expansion within a season Position of new ramets is determined stochastically and leads to patch expansion within a season

Transfer of energy between shoots and Additional transfer of energy between ramets roots/reproductive organs only

The basic units of the calculations are the ramets. Each ramet has an exact x and y coordinate and is associated with one grid cell. A clone can expand over more than one grid cell. The biomass growth of each ramet is calculated with time steps of 1 day and depends on photosynthesis and local environmental conditions. During their growth, the shoots reserve an increasing part of their net production for creating rhizomes and new ramets. The positioning of new ramets is determined by the length of the rhizome and a stochastic component. Above a certain shoot length, part of the production will be transferred to neighbouring ramets (clonal integration). During one season, apical growth of rhizomes and ramets and branching lead to the origin of a large clonal system (‘plant’). At the end of the growth season, part of the biomass is redirected to overwintering organs (‘turions’) in the sediment, and the rest of the aboveground biomass dies. The turions will develop new sprouts at a preset day in spring. Many of the default parameters have been derived from Lake Constance field data (Tab. 7.2).

67 Chapter 7

Each grid cell has environmental variables associated with it (vertical light attenuation, water level and nutrients). The grid dimensions can be defined by the user.

In the following, we will describe the model focusing on the new features and summarizing the parts that are taken from the earlier version of the model (Van Nes et al. 2002, 2003). The model is implemented in Delphi 5.5, an object oriented version of Pascal. It is freely available on: http://www.dow.wau.nl/aew/charisma/.

Tab. 7.2 Default parameters for the clonal model of P. perfoliatus. Source: a - calibrated, b - estimated from field observations (Wolfer, unpublished results), c - assumed to be the same as P. pectinatus (Scheffer et al. 1993, Van Nes et al. 2003), d - Wolfer & Straile 2004b, e - Wolfer & Straile 2004a

Parameter Description Unit Value HAge Half saturation for senescence d 100 b PAge Exponent in Hill function for senescence 2 b pMax Maximal gross photosynthesis hr-1 0.006 a Q10 Q10 for maintenance respiration 2 c Resp20 Respiration at 20 °C day-1 0.00206 c MaxLength Maximum shoot length m 5 b MaxWeightLenRatio Mean weight of 1 m young shoot g m-1 0.8 e RootShootRatio Proportion of shoot allocated to the roots fraction 0.06 b FracPeriphyton Fraction of light reduced by periphyton fraction 0.2 c hPhotoLight Half-saturation light intensity (PAR) for photosynthesis µE m-2 s-1 30 c hPhotoTemp Half-saturation temperature for photosynthesis °C 14 c pPhotoTemp Exponent in temperature effect (Hill function) for photosynthesis 3 c sPhotoTemp Scaling of temperature effect for photosynthesis 1.35 c HTurbReduction Half saturation biomass of light attenuation reduction g m-2 100 c pTurbReduction Power in Hill function of light attenuation reduction 1 c PlantK Light attenuation of plant tissue m2 g-1 0.02 c cTuber Fraction of turion weight lost daily when sprouts start growing fraction 0.1 c MaxTurionWeight Maximum weight of turions g 0.2 b MeanTurionWeight Mean (initial) weight of turions g 0.2 b MinTurionWeight Minimum weight of turions g 0.1 b TurionFraction Fraction of biomass allocated to turions fraction 0.05 a TurionGerminationDay Day of turion sprouting 114 b TurionMortality Annual mortality of turions yr-1 0.1 b TurionReproDay Formation day of turions 250 b Alpha Parameter that determines the biomass allocation to subsequent 0.6 a rhizomes and ramets MeanRhiAngle Mean rhizome angle rad 0 d NewRametLength Length of the ramet at which it creates a new rhizome and ramet m 0.15 a PBranchingLong Branching probability of rhizomes fraction 0.1 d RhizomeWeightPerM Average weight of 1 m rhizome g m-1 0.5 e SDRhiAngle Standard deviation of rhizome angle rad π/9 d TurionAngle Angle of the turion with the parent rhizome rad π d TurionAngleRange Range in the angle of the turion with the parent rhizome rad 0 d

68 Modelling

Production and respiration

Most rules for production and respiration are taken from the Charisma model (Van Nes et al. 2003), but bicarbonate is not assumed to be limiting, and senescence is described in a more detailed way.

The daily growth of each shoot (∆W) depends on the photosynthesis (P), the maintenance respiration (Rm), the import from previous ramets (T1) and export to subsequent ramets (T2) (Fig. 7.1):

∆Ws = Ws P −W Rm +T1 −T2 (1)

temperature light limitation

primary production respiration integrated over day senescence and shoot length

previous ramets net production rhizomes, integrated over day or turion and shoot length next ramets

shoot biomass shoot length

root biomass

Fig. 7.1 Schematic overview of the factors that determine the growth of each shoot

Maintenance respiration (Rm) is arbitrarily taken in the middle of the range of values published by -1 Madsen and Adams (1989) and Ikusima (1970) for miscellaneous submerged plants (r20 = 0.024 g g -1 d at 20 ºC). Temperature dependence of the respiration is formulated using a Q10 formulation (default value 2):

T −20 10 Rrm = 20 Q10 (2)

Only the shoots take part in the primary production (P). The maximum photosynthesis (Pmax) is limited by the in-situ light intensity at plant leaves (I), temperature (T), and the age of the ramet (A):

I 1.35T 3 H 2 A (3) P = Pmax 3 3 2 2 I + H I T +14 H A + A

o The parameter Pmax represents the specific daily production of the plant top at 20 C not having light limitation. Light limitation is described by a Monod function, HI is the half-saturation constant. Irradiance follows a daily as well as a yearly cycle and light is attenuated in the water column both by the water and the vegetation (Van Nes et al. 2003).

Temperature dependence of photosynthesis was fitted to values of Potamogeton pectinatus (Scheffer et al. 1993). As field observations indicate that older ramets are often covered with

69 Chapter 7

periphyton that reduces irradiance and consequently production, a factor for aging was added as a Hill function (HA is the half-saturation constant of aging).

In the present analyses we assumed that the plants are not limited by phosphorus, nitrogen or carbon in the water, but the model can optionally account for such limitations.

As light varies during the day and with water depth, we integrated photosynthesis over the day and over the length of each shoot by three-point Gaussian integration (Goudriaan 1986). The in-situ light conditions were averaged for each grid cell taking self-shading into account (see below). As accounting for self-shading can be very computer intensive and as the biomass changes are relatively slow, we calculated the light attenuation in each grid cell once in 7 days (at 15 depths distributed evenly over the water column) and interpolated linearly between the calculated points.

For vertical light attenuation, temperature and water level we used 10 years of interpolated data from Lake Constance. For more theoretical questions we averaged these years to avoid stochastic differences between years.

We also included water clarification by macrophytes (Scheffer 1998, Van Nes et al. 2003). In Lake Constance this effect is probably not important as the water is usually rather clear, but this effect is essential for the existence of alternative stable states (Scheffer 1998; Van Nes et al. 2002).

In most simulations we used a grid of 10 x 10 cells of 0.2 m x 0.2 m. To avoid edge effects, ramets that grew beyond the border of the grid were coupled to the opposite grid cells.

Biomass allocation, clonal growth and clonal integration (transfer)

The produced biomass is allocated to shoots, rhizomes and roots. The root is a fixed proportion of the shoot that does not take part in photosynthesis. The height of the young shoots increases proportionally with their biomass until the water surface or a predefined maximum height has been reached. If the net production of a shoot is positive, a part is allocated to the rhizomes. The proportion of the production that is allocated depends on the height of the plant: very small shoots do not attribute, fully-grown shoots add all their production to the rhizomes. Between these extremes, the allocated proportion increases linearly with shoot height. The biomass is assumed to be transported acropetally from older to younger ramets (Evans & Whitney 1992, Marbà et al. 2002). Each younger shoot adds a part (α) of the transferred biomass to its net production. The remainder goes to the next shoot (Fig. 7.2). Thus, the transport to the shoot of ramet i (Ti) sums up to:

i k−1 i > 1 (4) Ti = ∑ Si−k α (α −1) k =1

in which Si is the surplus production of shoot i (i = 1 is the oldest shoot). For simplicity we assume that there is no cost for transport. Biomass that is transferred to full-grown shoots, will not be used but immediately be transferred further to the next ramet. If there is a branch, the transported biomass is divided equally over both branches. The surplus production of the last ramet goes into the growth of the new rhizome which has a fixed weight-length ratio. If the last shoot reaches a certain fixed height, the growth of the rhizome is ceased and a new daughter ramet is created.

The main rhizome may have an (usually small) angle with the previous rhizome which is drawn from a normal distribution. Each ramet has a certain probability of branching which is drawn from a Bernoulli distribution. The angle of the branch with the main rhizome is also drawn from a predefined normal distribution.

70 Modelling

ramet

shoot

S 1 α S 2 α S3 α S4 T 2 T3 T4

growth root rhizome spacer apex

Fig. 7.2 Scheme of the clonal growth of P. perfoliatus. Si is the surplus production of shoot i, Ti is the transfer to shoot i, α is a parameter that determines which part of the transfer (Ti) goes to the next shoot

Mortality

The model accounts for (a) a fixed background mortality and (b) mortality at the end of the season. Background mortality results either in a weight loss or in a stochastic loss of the complete shoot depending on whether it has already reached the surface or not (Van Nes et al. 2003). At a pre-set day in autumn, all shoots die, and their biomass is reallocated to the rhizomes. At the last (youngest) ramet of each branch of the clone, a turion is created (Wolfer & Straile 2004b). Since it has been shown that turion production often depends on plant biomass (Spencer et al. 1993), the turion receives a proportional part of the accumulated energy in each ramet (Fig. 7.3). The amount of energy (biomass) that is transported to each turion is a fixed part of the rhizome biomass up to a fixed maximum. The turion is only created if the rhizome biomass exceeds a certain minimum weight. The growth direction of the clone can be continuing, reversed or stochastic in the next season.

Turion

(optional) reversion of growth direction next year

Turion

(optional) reversion of growth direction next year

Fig. 7.3 Schematic representation of a clone and the transport of energy from the ramets (circles) to the turions (triangles). At the end of the season the turions get a proportional part (same color) of the accumulated biomass in the rhizomes and shoots. If enough biomass is accumulated, the turion produces a new clone in the next season in which the growth direction is optionally reversed

71 Chapter 7

Scanning of asymptotic regimes

In the previous version of the model we showed that there can be two alternative stable states in the model due to a feedback of vegetation on their light climate (Van Nes et al. 2002, Van Nes et al. 2003). In the current more detailed model, we show whether the model still has alternative stable states, by analysing of the effect of increasing and decreasing light attenuation on the equilibrium biomass of P. perfoliatus. Vertical light attenuation was slowly increased in small steps while the model was not reset. After a period of stabilising, the biomass at the end of the growth season was plotted for 5 years. When the water had reached maximum vertical light attenuation and the vegetation had disappeared, the same procedure was repeated backwards (i.e. vertical light attenuation was reduced). A small import of turions prevented total extinction of the plants. If the model has alternative stable states it will show a hysteresis in the response, i.e. the vegetation will recolonize at a lower turbidity than the turbidity at which they disappeared.

Monte Carlo sensitivity analysis

We applied a Monte Carlo sensitivity analysis to select the parameters which have the strongest impact on clonal architecture. We generated 20 000 sets of parameters, drawing all parameters randomly and independently from uniform probability distributions within ranges of ±10 % around the default values. Three years were simulated with each parameter setting, and the model results (mean shoot length, mean spacer length, fraction belowground biomass, mean number of ramets per clone, approximate expansion area, mean expansion per clone) were stored at three dates per year (at day 189, 219 and 249) of each year. At the end of each simulation, the model was reset to the starting number of turions per grid cell (5 m-2). Sensitivity coefficients were defined by linear regression between the parameter values and each model output value, scaled by the ranges used for each parameter (Klepper 1989). Cluster analysis (average linkage) was used to form groups of parameters that had the same or opposite effect on the qualitative model results. As similarity measure the absolute sine of the angle between the vectors of sensitivity coefficients was used. As measure of the total sensitivity the length of this vector was used (Klepper 1989).

In a Monte-Carlo analysis there is a probability that a parameter has a positive sensitivity coefficient because it is related to the model outcomes by chance. To determine the significance level of the sensitivity coefficients, we added 100 dummy parameters that had no effect on the model (Van Nes et al. 2003). These dummy parameters were arbitrarily set to 1. As the real model parameters, the dummy parameters were drawn from uniform distributions, and their ‘sensitivity coefficients’ were calculated. The 0.02 significance level of the sensitivity coefficients was estimated by the 98% percentile of the sensitivity coefficients of the dummy parameters. In the cluster analyses we only included parameters with a significant sensitivity coefficient.

Results

Architecture and growth parameters

The simulation produces complex architectural patterns of shoot and rhizome connections (Fig. 7.4) which compare well with field data from Lake Constance (Wolfer & Straile 2004b). For example, the observed increase in length of younger rhizomes during one season (Wolfer & Straile 2004a, 2004b) is reproduced by the model.

72 Modelling

Shoot length (m) 3

2.00 1.6 1.2 Y 0.8 2.0 0.4 1.6 1.2 0.0 0.8 0.4 0.0 X

Fig. 7.4 Example of a generated pattern (clone growth during one year)

A 20 years simulation starting with 5 propagules m-2 yielded the following results: Annual curves of most growth parameters (Fig. 7.5a) start in May, increase more or less exponentially, peak in August, and fall back to zero after September. With regard to area, there is a net patch growth each year (difference between end-of-season and before-season turion occurrence). In 20 years (Fig. 7.5b), biomass, ramet number and total patch size initially increase and then stabilize after 10 - 15 years; ramet number per plant and rhizome length decrease slightly, spacer length, shoot length, belowground biomass and limitation stay about the same.

The values of simulated growth parameters such as shoot length, rhizome spacer length, rhizome length etc. are in the range of the field data (Wolfer & Straile 2004b; Tab. 7.3). In general, growth parameters of macrophytes are extremely variable and also the cited reference only covers a small part of the repertoire. The simulated biomass and ramet densities occur frequently in-situ (Wolfer, unpublished observation); for the cited field work, lower densities were deliberately chosen. The simulated number of ramets per plant is rather low, but this is partly due to the high equilibrium density of ramets, implying strong light limitation. At lower ramet densities the number of ramets per plant is quite realistic (Fig. 7.6).

Tab. 7.3 Model results of growth parameters in comparison to in-situ and mesocosm assessments (average site values)

Parameter Model result In-situ (Wolfer and Mesocosm (Wolfer and Straile, 2004b) Straile, 2004a) biomass 280 g dw m–2 17 - 200 g dw m-² na ramets 350 m–2 90 - 210 na ramets per plant 5.5 5 - 8 7 - 13 shoot length 100 cm 12 - 66 cm na rhizome spacer length 8.5 cm 7.6 - 8 cm 7 - 10 cm rhizome length 0.4 m 0.4 - 0.63 m 0.6 - 1.25 m

73 Chapter 7

. 300 300 ) (g

s 200 200 as om i 100 100 b

0 0 r

e 400 400 b

m 300 300 nu t

e 200 200 m

ra 100 100 0 0

t 6.0 12 n a l 4.5

r p 8

pe 3.0 s t

e 4

m 1.5

ra 0.0 0

) 10.0 15 m c ( 7.5 h

t 10

ng 5.0 e l r 5 e 2.5 ac p

s 0.0 0

) 0.4 0.8 m (

h 0.3 0.6 t

ng 0.2 0.4 le e 0.1 0.2 om z i h

r 0.0 0.0

) 1.2 15 m ( h t 0.8 10 g n e t l

o 0.4 5 o h s 0.0 0 36 40

35 30 ²)

(m 34 20 ea

ar 33 10 32 0 JFMAMJJASOND 90 92 94 96 98 00 02 04 06 08 10 month year

Fig. 7.5 20 years run, (a) year 2009 (b) 1990 - 2010.

74 Modelling

16

12 one cl r

pe 8 s t e m

a 4 R

0 0 100 200 300 400 500 Biomass (g m-2)

Fig. 7.6 Relation between the mean number of ramets and biomass at day 249. Data based on 10 replicates of 20 years of simulation. The gap between biomass 200 and 300 g m-2 is due to the exponential phase of the growth in macrophytes

Stabilization of biomass and ramet number after an initial colonization phase is common in aquatic macrophytes (Duarte & Sand-Jensen 1990). Ramet number per plant and total rhizome length decrease because of higher density and the resulting light limitation (Wolfer & Straile 2004a).

Average spacer length, shoot length and belowground fractions do not change much during the 20 simulated years. The behaviour of spacer length does not comply with observations which show that average spacer length decreases with increasing light limitation (Wolfer & Straile 2004a). This decrease must be an active light foraging strategy of the plants (shoot growth for light capturing at the expense of rhizome growth) which could not be reproduced by the simplified plant strategy in the model.

As in Wolfer & Straile 2004a, we found lower patch expansion at low plant density and at higher density, and highest patch expansion at medium density. Long-term stabilization of patch sizes under model conditions (no disturbance, homogeneous environment) is probably related to field observations, where patches appear at the same site and in a very similar size from year to year (Walser 1996; Gafny & Gasith 1991). In-situ, patches may even keep their distinct shape over long time periods, although, according to expectations, they should expand until they fill up homogeneous lake compartments. The growth curves suggest that this could be a result of the stabilization of biomass and ramet number which in turn are due to light limitation in a dense patch.

Branchings

A variation of the branching probability per plant from 0 to 1 reveals strong consequences for equilibrium summer biomasses and total number of ramets (Fig. 7.7). Plants with zero branching die. Due to the clonal growth architecture of P. perfoliatus, each plant with zero branching produces exactly one turion at the end of the rhizome axis (Wolfer & Straile 2004b). Under the assumption of zero mortality, the number of plants will remain the same forever. However, turions are usually subjected to waterbird foraging, fungus infection, or mortality of stochastic reasons. Additionally, it is assumed that turions also fail to propagate if their weight is too low (for example due to bad shoot

75 Chapter 7

biomass growth in the previous season). The model shows that, if there is no branching, even the lowest turion mortality will eventually lead to the extinction of the plant.

)

-2 350 300 (a) 250 200 150 100 50 0 Summer biomass (g m 0 0,2 0,4 0,6 0,8 1 Probability of branching

400 350 (b) )

-2 300 m

( 250 s

er 200 b 150 100 Num 50 0 0 0,2 0,4 0,6 0,8 1 Probability of branching

Fig. 7.7 Modelling results on the influence of branching probability on (a) summer biomass and (b) number of ramets in P. perfoliatus after a simulation period of 20 years

Maximum summer biomass and maximum density are achieved at the (relatively low) branching intensity of 0.1 - 0.2. At higher branching rates, biomass and number of ramets decrease again. The optimum branching probabilities found by the model are equal to those found in the field in patches with favourable plant growth (Wolfer & Straile 2004b), and lower than those found under experimental conditions (Wolfer & Straile 2004a), probably because of lower water level and more light in the experimental water basin. The decreasing biomass at higher branching is due to our assumption that energy reserves at the end of the growing season are equally divided up into the number of turions. Branching increases the number of turions, therefore individual turions will have not enough reserves to sprout successfully in the next season.

Spatial expansion with and without reversing turions

The main rhizomes axes of P. perfoliatus usually grow relatively straight and terminate with a turion at the end of the season. During winter, the rhizomes decompose and only the turion remains intact. In the following year, the growth of the new rhizome can theoretically (a) be random, (b) continue in the

76 Modelling

direction of the previous year, or (c) reverse its direction. In practice however, the growth direction of the new rhizome is probably not random but rather determined by the position of the meristem along the clonal growth axis (Watson & Cook 1982). There is some indication of growth reversal, since turions grow basipetally in the sediment in a hook-like fashion (Wolfer & Straile 2004b).

The model allowed us to compare the influence of random, continuing, and reverse rhizome growth. There are dramatic effects on patch density and expansion after a 20 year simulation (Fig. 7.8). Reversal leads to a very compact patch with high density (Fig. 7.8b). Continuation of growth into last years direction creates a larger, sparser patch (Fig. 7.8d) and a random direction lies in between the two (Fig. 7.8c).

20 (a) (b)

15

10

5

m 0 e t 20 e (c) (d) r 15

10

5

0

0 5 10 15 20 0 5 10 15 meter

Fig. 7.8 Patch density and expansion in dependence to rhizome growth direction relative to the previous year. (a) common starting condition (b) growth reversal (c) random growth (d) continuing growth (simulation period 20 years)

These model results show that not only extrinsic growth conditions but also intrinsic architectural growth rules are of major importance for the patch characteristics. There are intra-specific feed-backs of patch density and structure such as modification of light, nutrients, and sediment detritus (Cebrián & Duarte 2001), regulation of flowering (Thompson et al. 1990), or susceptibility against wave attack (Coops et al. 1991). The observed differences in patch size, pattern and density will have impacts on processes in the littoral ecosystem, such as element cycles (Barko & James 1997), as well as on

77 Chapter 7

littoral food webs (Lauridsen et al. 1996) and the distribution and age-characteristics of littoral fish (Weaver et al. 1997) and macro-invertebrates (Webster et al. 1998).

Scanning of asymptotic regimes

Although the new model ‘CLOMO’ is much more complex than the original model Charisma, the results were quite similar (Van Nes & Scheffer 2005). Starting from the turbid state, P. perfoliatus growth tolerates less turbidity, than starting from the clear state. We can therefore show that Potamogeton influences its own environment through feedback mechanisms and that the system can have two alternative stable states at higher turbidity: one with vegetation and one with little or no vegetation. The main difference is that the zone with alternative equilibria is much smaller in the clonal model (Fig. 7.9) which is due to the fact that P. perfoliatus is assumed to have a much smaller effect on water clarity than P. pectinatus.

Biomass (g m-2) 500

400

300

200

100

0 0 1 2 3 4 Vertical light attenuation (m-1)

Fig. 7.9 Simulation of the response of P. perfoliatus biomass to increasing and subsequently decreasing turbidity without resetting the model

Monte Carlo sensitivity analysis of clonal architecture

Cluster analysis of the sensitivity coefficients shows that there are only two clear clusters of parameters that determine the clonal architecture. The strongest effect is due to a group of parameters with a strong effect on photosynthesis. This group includes hPhotoTemp, pMax, sPhotoTemp and hAge (Fig. 7.10). These parameters have a strong effect on the expansion area and the number of ramets per clone (Fig. 7.11). In an earlier analysis, they also had the strongest effect on biomass and numbers in the model Charisma (Van Nes et al. 2003).

78 Modelling

Dendrogram using Average linkage

0.624 RootShootRatio 0.431 PBranchingLong 0.968 MaxWeightLenRatio 0.675 MaxTurionWeight 0.661 pPhotoTemp 1.586 NewRametLength 5.245 Alpha 1.907 RhizomeWeightPerM 0.574 PlantK 0.413 Resp20 2.037 HAge 6.518 hPhotoTemp 1.329 hPhotoLight 1.370 PAge 5.979 pMax 5.796 sPhotoTemp 7.245 TurionGerminationDay

0.0 0.2 0.4 0.6 0.8 1.0 Abs. sine distance

Fig. 7.10 Cluster analysis of the sensitivity parameters of the model. Parameters with the same or opposite effect on the model results of 4 subsequent years at 3 points in time (day 189, 219 and 249). Sensitivity was based on the following model output: mean spacer length, mean shoot length, fraction belowground biomass, mean number of ramets per clone, approximate expansion area, mean expansion per clone. All parameters shown have a significant effect on model results (p > 0.01)

A second smaller cluster (Fig. 7.10) includes the parameter Alpha which determines the fraction of transported biomass that is used for the growth of subsequent ramets (Fig. 7.2). The parameters determining the length/weight relationship of the rhizomes (RhizomeWeightPerM) and the length of the shoots at which a new ramet is created (NewRametLength) are also included in this cluster (Fig. 7.10). These parameters also have a strong effect on the clonal expansion because they determine total rhizome length, but also on rhizome spacer length (Fig 7.11). Furthermore there are some parameters that belong to no cluster, but also have a significant effect. It seems plausible that the day of turion germination (TurionGerminationDay) influences the growth of the clones because any additional day of growth increases the biomass. The importance of earlier germination is also confirmed by experiments and in-situ observations (Spencer & Rejmanek 1989, Spencer et al. 2000). In the model, it strongly affects all variables in the first year but has much less effect the next years.

79 Chapter 7

0.5 a. hPhotoTemp Approx. expansion area (m2)

0 t Fraction root and rhizome n e i c i f

f Mean expansion per clone (m)

e -0.5 o c

y Mean length (m) it v i

it -1

s Mean no. of ramets per clone n e S -1.5 Mean spacer length (m)

-2 91/01 92/01 93/01 94/01 95/01

2 0.5 b. pMax c. Alpha

1.5 t t n n 0 e e i c ci i i f f f e ef 1 o c y y co -0.5 t it i v iv i t

0.5 it s n e ensi S S -1 0

-0.5 -1.5 91/01 92/01 93/01 94/01 95/01 91/01 92/01 93/01 94/01 95/01

3.5 0.5 d. TurionGerminationDay e. MaxWeightLenRatio 3 t t n 2.5 n e e i i c c i i f f f 2 f e e o o c c y 1.5 y 0 it it iv iv it it

s 1 s n n e Se 0.5 S

0

-0.5 -0.5 91/01 92/01 93/01 94/01 95/01 91/01 92/01 93/01 94/01 95/01

Fig. 7.11 Time course of the relative sensitivity coefficients of various representative parameters (expressed as sensitivity coefficients) on model outcomes

80 Modelling

Discussion

Comparison with other models

The models available for clonal plant growth can be classified into empirical models and mechanistic models (Carr et al. 1997).

(1) The majority of aquatic plant models are mechanistic biomass models that calculate plant growth from physiological processes such as photosynthesis but disregard clonal architecture. These models can be (a) non-spatial (Collins & Wlosinski 1989, Scheffer et al. 1993, van Dijk & Janse 1993, Davis & McDonnell 1997, Hootsmans 1999, Asaeda & Karunaratne 2000, Calado & Duarte 2000, Best et al. 2001) or (b) spatial (Wortmann et al. 1997, Van Nes et al. 2003). Compared to such models our model needs many extra parameters that have to be assessed or calibrated. Complex rules have to be set up which make the model complex and prone to uncertainties. The advantage of our model is that it gives more insight into plant growth strategies than non-architectural mechanistic models. The model allows us to study the effects and implications of specific growth rules for growth characteristics. For example: is branching necessary for successful propagation? It also allows for insights into patch formation and growth.

(2) In contrast to mechanistic models, there are some empirical ‘design models’ of clonal architecture (Bell 1976, Bell et al. 1979), including submersed species (Molenaar et al. 2000). Empirical models use the ranges and standard deviations of rhizome lengths, branching frequency, branching angles etc. and generate simulations of growth patterns. They provide insights into possible variability of patterns but can not explain the underlying mechanisms nor simulate the patterns under different environmental conditions. In contrast, our rhizome architecture is growth related and consequently ruled by environmental factors.

(3) Deterministic and stochastic spatio-temporal models (cellular automata) are also used for clonal plants in competition models (Colasanti & Grime 1993, Chiarello & Barrat-Segretain 1997, Balzter et al. 1998) but they describe clonal expansion on a grid basis using simple empirical rules, or random processes. They can be used to predict patterns on a larger scale, but being empirical, they do not provide mechanistic insight.

(4) A range of models, mainly for terrestrial plants, takes into account both mechanistic and empirical aspects (Cain & Cook 1988, Callaghan et al. 1990, Cain 1994, Oborny 1994, Cain et al. 1995, Evans & Cain 1995, Cain et al. 1996, Klimes 2000), but neither of these models deals with submersed macrophytes nor clonal integration.

(5) Models of clonal integration exist for theoretical plants and are usually non-spatial (Stuefer et al. 1998, Chesson & Peterson 2002, Suzuki 2002). Only a few spatial models include clonal integration (Oborny et al. 2000, Herben & Suzuki 2002, Oborny & Kun 2002), but in contrast to our model these models are not truly mechanistic, ignore clonal architecture or have a rigid architecture, i.e. do not assume any plastic adjustment of clonal architecture to resource availability.

To our knowledge there is no other model that includes all following components: submersed species, mechanistic biomass growth (dependent on environment), individually based, spatially concise, non-rigid clonal architecture, and clonal integration.

81 Chapter 7

Controlling complexity of the model

Being complex and detailed, our model has also several drawbacks (Van Nes & Scheffer 2005). The model includes many unknown parameters and processes, resulting in a high uncertainty of the results. Although our results were qualitatively correct, we cannot be sure that we get good results for the correct reason. Therefore, it is not sure whether this model can also be applied to other lakes and sites.

In particular, the results of the simulation of branching variability are sensitive to some parameters for which we do not have good field data, such as fraction of biomass allocated to turions or minimum turion weight for sprouting. The description of poorly known processes was kept as simple as possible in the model, ignoring more flexible plant strategies like omitting small side branches or forming a few big turions to overcome the problem of too small turions. Moreover, the model does not include biomass allocation strategies and foraging by means of rhizome length control, as we lack enough knowledge of these processes.

We tried to manage the complexity of the model by implementing detailed visualisation of the results which we regard as essential for keeping track of the complexity (Grimm 2002). Our modelling environment could visualise the clones and the various features of the ramets, including shoot length, roots and rhizomes, transport between shoot and rhizomes and between rhizomes and patch expansion in 2D and 3D graphs. Furthermore, it was essential that the generated patterns could be compared with real data of P. perfoliatus in Lake Constance at the same level of detail (Wolfer & Straile 2004a,b). Keeping track of several properties of the clones (spacer length, biomass, number of ramets) simultaneously, we could reduce the probability of getting realistic results for wrong reasons. Finally, we compared our results with a simpler version of the model (Van Nes & Scheffer 2005), namely the original Charisma model. This way we could show that some behaviour of the model, like the existence of alternative stable states, is not sensitive to the complexity of clonal growth.

The benefits of our model approach are detailed ecological insights through qualitatively correct results. The model has also provided first insights into the understanding of growth architecture and patch expansion of submersed macrophytes. The strong graphical user interface makes it possible to analyse the existing patterns and find the underlying mechanisms. In future, the model will be used for more detailed analyses and to generate hypotheses about macrophyte growth rules which can be tested experimentally. It can easily be extended by additional parameters and processes, for example wave mortality, or heterogeneous grid cells in order to provide an improved picture of the growth of submersed macrophytes.

82 Modelling

83 Chapter 7

84 Synthesis and discussion

CHAPTER 8. SYNTHESIS AND DISCUSSION

Our study on Potamogeton perfoliatus provides the first insights into clonal architecture of this species and its plasticity in response to the environmental factors light and nutrients. The results show that architectural growth rules govern spatial patterns, but that local environmental factors may also affect architecture. Based on our experimental results and the literature, we are able to construct models that simulate patch growth quite well. Nevertheless we are still far from understanding all details of patch growth found in-situ. A number of factors that potentially influence patch patterns have not been investigated in this study and provide possibilities for future research (Fig. 8.1). Below, a summary of our findings with respect to architectural features and the effects of nutrients, light and clonal integration will be presented. Subsequently, some other factors that seem likely to be influential in the field will be highlighted.

species factors physiology clonal integration

clonal architecture spacer length habitat parameters branching frequency light branching angle nutrients growth direction sediment (texture) size, density exposure and shape of water chemistry macrophyte water temperature patch

growth biomass

biotic factors herbivory epiphytes competition

Fig. 8.1 Factors that influence patch characteristics; bold font: investigated in this study; standard font: studies available; italic font: further research required

Architecture

The investigation of small-scale quadrats (1 m²) of P. perfoliatus by divers in Lower Lake Constance (Chapter 3) has provided the first detailed data on architectural features such as number of ramets per clonal fragment, rhizome spacer length, branching angles of rhizomes, frequency of branching and turion formation. We found recurring architectural rules, among them a particularly interesting feature that has not been described before for this species: spacer length tended to increase considerably with the order of rhizome spacers i.e. the younger the spacer, the longer. We attributed this to the larger amount of photosynthetic activity of both individual shoots and the whole plant and to the transfer of resources to younger sections of the clonal system. Spacing between ramets of the same

85 Chapter 8

rhizome system was larger than spacing between nearest (non-clonal) neighbours which has interesting implications for competition.

The architectural information collected was used as input in the new simulation model ‘CLOMO’ (Chapter 7). Using parameterization of data from our experiments, the model generated distinct patches similar to those in the field, rather than all-covering ‘plant carpets’. The simulations also revealed that the direction of rhizome growth relative to that in the previous year determines the patch density and patch expansion rate. Moreover, the model showed that a lack of branching will lead to the extinction of the clonal fragment if the turion (i.e. the propagation unit) is non-productive. The reason for this is the fact that only one single turion is formed at the tip of each unbranched rhizome axis, whereas each branch may produce one additional turion.

Besides architectural mechanisms, the explanations for different patch sizes and shapes of P. perfoliatus lie also in various environmental factors that cause plastic growth of plants. As summarized in the following sections, our results demonstrate that clonal growth and patch architecture are affected by nutrients, light conditions and clonal integration (Tab 8.1).

Nutrients

Our experiments have shown that nutrients control parameters of clonal growth morphology such as branching, rhizome spacer length and biomass allocation of P. perfoliatus as well as total productivity, density and size of patches. The lack of P. perfoliatus patches in parts of Upper Lake Constance as well as the poor growth patches we found next to good growth patches in Lower Lake Constance can be explained by insufficient nutrients in the sediment porewater (Chapters 3, 5). Only P. perfoliatus shoots planted with additional nutrient supply (slow-release N-P-K fertilizer) survived in the Upper Lake, whereas unfertilized controls vanished eventually. Moreover, poor growth of plant patches in Lower Lake Constance showed increased growth after fertilization. An associated tissue analysis confirmed the positive relationship between nutrient availability, growth and tissue nutrient concentration. By contrast, in the mesocosm, fertilization did not enhance productivity on naturally nutrient-rich sediment of Lower Lake Constance. Since productive sites are dominant in the Lower Lake, we assumed that the unproductive sites were a product of either sediment disturbance or nutrient depletion. As an example of the influence of nutrients on the clonal architecture of P. perfoliatus, we showed increased branching at higher nutrient availability (e.g. Chapter 3). Furthermore, reduced allocation to rhizomes and (relatively) shorter rhizomes were found at higher nutrient concentrations (Chapters 3, 5). Under nutrient poor conditions, in contrast, P. perfoliatus allocated more biomass to roots and rhizomes - primarily by means of a reduction of shoot length - and developed (relatively) longer rhizomes.

Tab. 8.1 Simplified responses of clonal architecture on the growth factors nutrients, light and clonal integration as found in this study. + : positive relationship, - : negative relationship, = : indifferent, n.t. : not tested; integration may also change the local effects of light and nutrients on architecture (not shown in the table)

branching spacer length rhizome allocation root allocation nutrients + - / = - - light + + + = clonal integration n.t. - / + = / + - / +

86 Synthesis and discussion

Light / shade

Planting density had a clear effect on total productivity and parameters of clonal growth morphology such as branching, rhizome spacer length and biomass allocation of P. perfoliatus in a mesocosm. Nutrient or light competition seem obvious potential explanations for these effects. However, since fertilization did hardly affect growth on the same sediment (Chapter 5), we suppose that light effects are dominant. Remarkably, production of P. perfoliatus was highest at medium initial plant density. It remains uncertain if the lower production at low initial density is a result of physiological damage by excessive irradiation or other allee effects. Lower production at higher initial density might be due to light competition. The findings imply that the number of primary shoots i.e. propagules available in the beginning of the season can strongly influence patterns of growth and clonal reproduction. Higher light levels significantly increased the number of branchings, while medium and lower light levels resulted in reduced branching, implying that branching also works indirectly as self-regulation of ramet numbers (Chapter 4). Shading stimulated vertical growth of shoots at the expense of rhizome allocation and rhizome spacer length.

Clonal integration

An experiment in which parts of the growing clones were disconnected showed that translocation of resources between different ramets of the same plant may modify the effects of shading (Chapter 6). We showed that P. perfoliatus ramets were capable of sharing resources acropetally (i.e. from the primary shoot in direction of the younger shoots) with recipient ramets. As a consequence, integrated ramets improve nutrient supply of younger ramets and influence the productivity of the patch in a positive way. However, overall dynamics of redistribution of resources may be quite variable. For instance, if ramets occupy patchy environments, support by other ramets may be unavailable to ramets at unfavourable micro-sites but may be shifted to ramets at better micro-sites. Besides, rhizome spacer length of integrated ramets decreased in the shade (Chapters 6, 4), whereas in single ramets, spacer length increased in the shade (Chapter 6). Moreover, root allocation was affected by integration and shading as well as the interaction between the two. This shows that integration is a complex process that likely works towards optimum performance of the total plant. Implementation of integration in our simulation model may therefore be important to improve realism.

Patch size and density

The reactions of biomass and rhizome spacer lengths of individual plants to environmental factors influence the characteristics of the patch as a whole. Patch sizes increase with decreasing branching, increasing rhizome spacer length, increasing rhizome allocation, or straight rhizome growth direction and decrease with increasing branching, decreasing rhizome spacer length or reversing rhizome growth direction (Tab. 8.2). The trends found in the course of our nutrient and light experiments had the following implications (Fig 8.2): Nutrient supply decreased rhizome spacer lengths and resulted in dense and compact patches with smaller patch expansion and discrete patch shapes whereas a lack of nutrients stimulated horizontal expansion and larger (but low density and low height) patches. Besides, increased nutrient supply also increased branching. Since higher branching rates decrease patch size, the reactions of spacer lengths and branching to increased nutrients both worked in the same direction, i.e. a reduction of patch size. Responses to the factor light were more complicated. On one hand, increasing rhizome spacer lengths with increasing light led to larger patches and reduced rhizome spacer lengths in the shade led to reduced horizontal growth and resulted in dense and

87 Chapter 8

compact patches with smaller patch expansion and discrete patch shapes. On the other hand, increased branching with light resulted in reduced patch sizes. Therefore, reactions of rhizome spacer lengths and branching to different light levels may partially cancel each other out.

Tab. 8.2 Simplified relationships between clonal architecture / growth factors and patch size as found in this study. + : positive relationship, - : negative relationship factor patch size branching - spacer length + rhizome allocation + reversing growth direction -

+ - + N, P - N, P

Fig. 8.2 Simplified responses of patch growth to light and nutrients as found in this study

Several aspects were not addressed in our study, but might nevertheless have an important effect on clonal growth. In the following, the potential role of exposure, periphyton, interspecific competition and feedback between plants and their environment will be briefly highlighted.

Exposure

Among the environmental factors that influence patch growth, exposure has received considerable attention. Wave action and currents affect submersed macrophyte stands directly or indirectly. The most obvious direct effect in shallow water zones is mechanical disturbance through uprooting (Keddy 1982, Lehmann 1997). Although P. perfoliatus is often present at exposed sites (Luther 1951, Schmieder 1995), waves can prevent patch establishment (Coops et al. 1991), impact patch growth, or destroy parts of patches (Walser 1996). It has been shown that exposure leads to more compact outlines of macrophyte patches (Fonseca & Bell 1998). Indirect effects of exposure on macrophytes may include light attenuation by particle resuspension or partial removal of periphyton (Coops et al. 1991, Fonseca 1996). This may affect species composition (Keddy 1983), biomass (Kautsky 1987, Duarte & Kalff 1990a) and percent cover of submerged macrophyte species (Fonseca & Bell 1998) as well as clonal architecture, patch characteristics and aspects such as shoot lengths (Kautsky 1987),

88 Synthesis and discussion

root number and length (Koch 1999), shoot branching and shoot internode length (Kautsky 1987), shoot density (Lundholm & Simser 1999) and stolon length (Puijalon & Bornette 2006). Importantly, plants also influence current velocity (Fonseca et al. 1983, Losee & Wetzel 1993, Koch 1996), implying a two-way interaction between plants and their environment.

Epiphyton and periphyton

Since epiphytes take up the same nutrients as macrophytes and additionally reduce light and CO2 availability (Jones et al. 2000), they may consequently decrease macrophyte chlorophyll contents, reduce macrophyte growth (Ozimek et al. 1991), or even kill the host plant (Jones et al. 2000). Blindow (1987) could show that epiphyte colonization differed between plant parts, sites and macrophyte species. As epiphyte species also affect the amount of marl encrustation (i.e. carbonate deposits), they may considerably influence patch characteristics. Since epiphyte layers increase during the growing season, fast growth is favourable. Moreover, it may well be that shading by epiphytes has similar effects to self-shading or experimental shading: reduced rhizome spacer lengths, less roots, more shoots, and reduced branching.

Interspecific competition

In the context of this study, we have tried to work on mono-specific stands where ever possible. However, P. perfoliatus occurs also in intermixed stands. In Lake Constance, for example, important associated species are Chara sp., Potamogeton pectinatus and Najas marina. In such intermixed stands, there is likely to be some competition for space occupancy, light and nutrients which might influence patch characteristics. Competitive reactions of P. perfoliatus include elongation and concentration of photoreceptive biomass at the water surface (light foraging, see below). Since simulations show that spatial aggregation has a dramatic effect in reducing the rate at which stronger competitors are able to exclude weaker ones (Silvertown 1992), we would expect shorter rhizome spacer lengths, and increased branching at higher levels of competition. Furthermore, in case of nutrient competition, foraging by root systems should increase (Campbell et al. 1991). Further competitive actions involve the change of sediment characteristics or the release of allelopathic substances and might affect the positioning of macrophyte patches.

Effects on own environment

Once a plant has established itself on a formerly uncolonized area, it will influence its own environment. After the sprouting of ramets and building up of biomass, patches act like sediment traps. During the vegetation period, suspended allochtonous and autochtonous organic matter (Sand- Jensen 1998, Lee et al. 2001) or sand (Sánchez et al. 2001) accumulates within the patch thereby changing both density of the sediment (Barko & Smart 1986b) and nutrient status (Szcepanska & Szcepanski 1973). In Lake Constance, detritus at macrophyte sites can amount to 10 cm thickness, but after the growing season, the site might be completely eroded by autumn storms (Walser 1996). Plants on old growth sites may profit from increased nutrients or might be negatively influenced by high concentrations of nutrients and / or toxic decomposition products. Indeed, the lowest density is sometimes found at the centre and patches may take on a ring-shaped appearance (Sánchez et al. 2001). These could, however, also be explained by nutrient depletion in the centre due to proliferate growth and nutrient uptake.

89 Chapter 8

Theories of architectural responses to habitat quality

Clonal responses and plastic morphological changes in architecture of modular plants as found in this study have received much attention and are often interpreted as strategies to optimize space occupation and utilization. Here we discuss whether our results are in line with the predictions from two main theories in this field: foraging theory and allometric growth rules.

The ‘foraging theory’ considers changes in branching frequency, biomass allocation and spacer length as mechanism for the increase of ramet density in favourable microsites, and as means of improving exploitation of resources such as light or nutrients (Hutchings 1988, Sutherland & Stillman 1988, De Kroon & Knops 1990, Vermaat & Verhagen 1996, Marbà et al. 1996).

Some definitions of foraging

Grime (1979): ‘the ability of plants to project leaves and roots into patches of high resource supply within the environment’ (the word foraging as such was not used; pertaining non-clonal species ).

Bell (1984): clonal plants consist of ‘feeding sites’ and spacers that connect the feeding sites

Slade and Hutchings (1987b): ‘plasticity [in stolons] enables clones to consolidate occupation of favourable sites, through intensive foraging, and to grow through less favourable sites’

De Kroon & Hutchings (1995): ‘the processes whereby an organism searches, or ramifies within its habitat, which enhance its acquisition of essential resources’.

Oborny & Cain (1997): ‘Plant foraging occurs when a plant exhibits any morphological plasticity that is selectively advantageous for resource acquisition at a particular spatiotemporal distribution of the resource.’

As predicted by the foraging theory, we found reduced branching at low nutrient availability (Chapter 5) or low light levels (Chapter 4). Such reduced branching allows for larger spatial extension and for searching for better sites whereas increased branching at higher nutrient and high light availability allows for denser stands and placement of more shoots at the favourable site (Oborny 1994, Sutherland & Stillman 1988). As in other plant species, branching patterns are thought to be regulated by hormones such as auxins, cytokinins and gibberellins which strengthen apical dominance of rhizomes in case of low nutrient availability (see Jackson et al. 1985).

In recent years, the foraging theory has been rejected by many authors including the original authors, especially with regard to spacer length of rhizomatous species (De Kroon & Hutchings 1995). The main reason is that many studies did not find shorter spacer lengths in favourable environments (e.g. Hartnett & Bazzaz 1985). Furthermore, experiments showed that plasticity in spacer length did not succeed in concentrating ramets in favourable habitat (De Kroon & Hutchings 1995). On the other hand, some rhizomatous species definitely showed foraging (Slade & Hutchings 1987a, De Kroon & Knops 1990, Dong et al. 1997, Klimes et al. 1997) and simulation models suggested that foraging can succeed in concentrating ramets in favourable patches (Sutherland & Stillman 1988, Oborny 1994, Cain et al. 1996). According to a literature review (Miller 1999), 15 % of examined rhizomatous species (all terrestrial) showed foraging. Unfortunately, the presented growth parameters generally suffer from a lack of completeness and comparability. In our study, as predicted by the foraging theory, we found longer rhizome spacers at non-fertilized sites (Chapter 5) which resulted in sparsely vegetated patches with a relatively large horizontal expansion rate. By contrast, the foraging theory

90 Synthesis and discussion

seems not to be appropriate to explain some of our other results such as decreasing rhizome spacer lengths with shade (Chapter 4) or equal spacer lengths of good and bad growth sites in-situ (Chapter 3).

The main problem of the foraging theory is that relevant studies mostly compare absolute rhizome spacer lengths at ‘good’ and ‘bad’ sites, disregarding that plants may differ considerably in biomass, shoot length, number of ramets and weight-density relationships. However, when resources increase, plant growth increases including both shoots and rhizomes. This is the main prediction of the ‘growth hypothesis’ according to which spacer length should increase at favourable sites as they allow a higher productivity (Hutchings & De Kroon 1994, De Kroon & Hutchings 1995, Stoll et al. 1998). However, if rhizomes increase at a slower rate than shoots, this foraging sensu latu might be easily overseen. In this case only a closer look at biomass allocation will reveal foraging, following the newer definitions of De Kroon & Hutchings (1995) or Oborny & Cain (1997). Principally, longer rhizomes can be achieved by means of increased biomass allocation into rhizome tissue, a reduced rhizome diameter or reduced tissue-mass density. The latter aspects have hardly received any attention by researchers and are not covered by this study. In our study, reactions of spacer length to changing resources were usually related to changes in biomass allocation. In-situ, rhizome spacers had the same length whereas rhizome allocation was five-fold higher on unfavourable sites compared to favourable sites (Chapter 3). Increased allocation to rhizomes under bad conditions, however, is only possible at the expense of shoot growth because the total biomass can at most be equal to that at a favourable site. In fact, we would expect total biomass to be lower as reduced resources decrease photosynthesis. Indeed, our study showed that the ramets on unfavourable sites had less biomass, were much smaller and had higher rhizome allocation than those on favourable sites. This is also in line with the ‘allometry theory’ which argues that allocation differences are a consequence of allometric growth rules rather than of morphological plasticity, as smaller plants frequently show a higher biomass allocation to rhizomes and roots (Duarte 1991, Gedroc et al. 1996, Müller et al. 2000). Allometric rules should not be confounded with morphological plasticity (Coleman & McConnaughay 1995, Huber & Stuefer 1997, Huber et al. 1999, Preston & Ackerly 2004) and have been complemented with the concept of ‘developmental reaction norms’ (DRN) to include the ontogenetic stage of plants (Schlichting & Pigliucci 1998).

We cannot determine whether spacer length in P. perfoliatus follows the foraging or the allometry theory. This is because both theories lead to the same predictions: favourable sites in the context of the foraging theory should produce large plants in the context of the allometry theory, and both result in shorter rhizomes respectively less allocation into belowground parts. Foraging and allometry theory differ only in terms of causation which is an active ‘search’ in the former and an ontogenetic growth rule which may also have been evolved through natural selection in the latter (Preston & Ackerly 2004). If we associate foraging with morphological plasticity, we have to compare responses of spacer length / belowground allocation to a resource at the same ontogenetic level which is rather difficult. Allometric relationships differ considerably between species (Duarte 1991, Müller et al. 2000) and adequate data for P. perfoliatus are missing. However, it seems very unlikely that, at the course of their ontogenetic growth, ramets might grow by more than 50 cm without any changes in rhizomes (Chapter 3). Therefore it seems likely that there is indeed some kind of foraging in P. perfoliatus and that its effects add to allometric rules of ontogenetic growth. The systematic increase of younger rhizome spacer lengths within the season suggests, however, that P. perfoliatus does not ‘search’ actively by sensing a favourable site and placing an individual ramet, as the earlier definitions of foraging would suggest (Slade & Hutchings 1987a,b,c) because in that case we would expect to find a

91 Chapter 8

more diverse picture. As the growing rhizome is initially rootless until a new ramet is established, the sensor for microsite quality remains unclear. We suggest that ramet placement is rather influenced by the conditions of the older ramets than the new ramet. This would be beneficial in environments with a scale of heterogeneity in the range of plant size where a systematically smaller or larger increase of all rhizome spacers of a plant may succeed in increasing overall fitness of the genet fragment. In addition, localized morphological responses of biomass allocation are highly important (De Kroon & Hutchings 1995).

With respect to foraging, the mechanism for different biomass allocation may also be interpreted in the light of stochiometric theory. Under nutrient limitation, for example, nutrients and CO2 invested in photosynthetic tissue can not further increase biomass production because of missing nutrients. Consequently, it would be better for the plant to invest nutrients either in rhizomes in order to grow faster out of the bad site, or in roots to capture more nutrients. Only when CO2 becomes the limiting factor, it is beneficial for the plant to produce aboveground photosynthetic active biomass. Likewise, allometry is based on ‘economic’ allocation for optimal resource acquisition (Bloom et al. 1985). Overall, plants must function as a balanced system in terms of resource uptake and use, and resources may not be allocated independently from each other (Hirose 1986, Ågren & Ingestad 1987). Higher allocation to roots in smaller plants might also be due to the need for anchorage, the lower possibility to take up additional nutrients over their surface and the faster uptake due to a higher growth rate.

The decreasing rhizome allocation with shade (Chapter 4) seems to contradict the allometry theory since mean shoot length and mean shoot weight remained the same. In case of shading, however, investment into aboveground biomass at the expense of rhizomes enabled shoots to escape the shading of neighbours by growing up to the light. The decrease of spacer length with shade can definitely be interpreted as light foraging if the more general definitions of Grime (1979) and Oborny & Cain (1997) are applied instead of the original definition of shorter rhizomes at good sites. When the foraging rule was originally formulated, research was done on terrestrial plants with rosette growth forms or on grasses which have the strategy to grow sideward out of the shade (Slade & Hutchings 1987a). As light diminishes sharply with depth, the better strategy is to develop tall shoots and to grow upwards out of the shade, partly by etiolation.

Besides appropriate placement of ramets, a clonal plant must optimally exploit the resources, e.g. take-up nutrients by roots. As predicted by the foraging theory, plants at the nutrient-deficient site allocated considerably more biomass to roots than plants at the rich site, and plants in the controls of the in-situ fertilization experiment allocated more biomass to roots than the fertilized plants (Chapter 3) whereas in the mesocosm, root allocation was not affected by fertilization due to the fact that biomass was not affected either (Chapter 5). This allocation rule has been confirmed by a variety of experiments with clonal terrestrial plants (review of Hutchings 1997) and submersed macrophytes (Barko & Smart 1981b, Kautsky 1991, Reusch et al. 1994, Wertz & Weisner 1997). Such shifts in allocation patterns under changing environments have been proven to maximize plant growth (Hirose 1987) and can be interpreted in the light of both foraging and allometry theories analogous to rhizome allocation (see above).

In summary, the dynamics of clonal architecture we found in P. perfoliatus are in line with predictions by allometric rules as well as foraging theory. While this is a poorly explored area of research in aquatic ecology, the responses we demonstrated may be the key to explaining the ecological success of clonally growing submerged plants.

92 Synthesis and discussion

93 Chapter 8

94 References

References

Abernethy V.J., Sabbatini M.R. & Murphy K.J. (1996) Response of Michx. and Myriophyllum spicatum L. to shade, cutting and competition in experimental culture. Hydrobiologia 340, 219-224. Aerts R., Boot R.G.A. & van der Aart P.J.M. (1991) The relation between aboveground and belowground biomass allocation patterns and competitive ability. Oecologia (Berlin) 87, 551-559. Agami M. & Waisel Y. (1985) Interrelationships between Najas marina and 3 other species of aquatic macrophytes. Hydrobiologia 126, 169-174. Agami M. & Waisel Y. (1986) The ecophysiology of roots of submerged vascular plants. Physiol. Veg. 24, 607- 624. Agami M., Beer S. & Waisel Y. (2001) Seasonal variations in the growth capacity of Najas marina as a function of various water depths at the Yarkon Springs Israel. Aquat. Bot. 19, 45-52. Agawin N.S.R., Duarte C.M. & Fortes M.D. (1996) Nutrient limitation of Philippine seagrasses (Cape Bolinao, NW Philippines): In situ experimental evidence. Mar. Ecol. Prog. Ser. 138, 233-243. Ågren G.I. (1985) Theory for growth of plants derived from the nitrogen productivity concept. Physiol. Plant. 64, 17-28. Ågren G.I. & Ingestad T. (1987) Root:shoot ratio as a balance between nitrogen productivity and photosynthesis. Plant Cell Environ. 10, 579-586. Aiken S.G. & Picard R.R. (1980) The influence of substrate on the growth and morphology of Myriophyllum exalbescens and Myriophyllum spicatum. Can. J. Bot. 58, 1111-1118. Alcoverro T., Romero J., Duarte C.M. & Lopez N.I. (1997) Spatial and temporal variations in nutrient limitation of seagrass Posidonia oceanica growth in the NW Mediterranean. Mar. Ecol. Prog. Ser. 146, 155-161. Alpert P. & Stuefer J.F. (1997) Division of labour in clonal plants. In: De Kroon H. & Van Groenendael J. (eds.) The Ecology and Evolution of Clonal Plants. Backhuys Publishers, Leiden, 137-154. Alpert P. (1999) Clonal integration in Fragaria chiloensis differs between populations: ramets from grassland are selfish. Oecologia 120, 69-76. Alpert P., Warembourg F.R. & Roy J. (1991) Transport of carbon among connected ramets of Eichhornia crassipes (Pontederiaceae) at normal and high levels of carbon dioxide. Am. J. Bot. 78, 1459-1466. Alpert P., Holzapfel C. & Benson J.M. (2002) Hormonal modification of resource sharing in the clonal plant Fragaria chiloensis. Funct. Ecol. 16, 191-197. Amsberry L., Baker M.A., Ewanchuk P.J. & Bertness M.D. (2000) Clonal integration and the expansion of Phragmites australis. Ecol. Appl. 10, 1110-1118. Anderson M.R. & Kalff J. (1986) Regulation of submerged aquatic plant distribution in a uniform area of a weedbed. J. Ecol. 74, 953-962. Anderson M.R. & Kalff J. (1988) Submerged aquatic macrophyte biomass in relation to sediment characteristics in ten temperate lakes. Freshw. Biol. 19, 115-121. Armellina A.D., Bezic C.R. & Gajardo O.A. (1999) Submerged macrophyte control with herbivorous fish in irrigation channels of semiarid Argentina. Hydrobiologia 415, 265-269. Asaeda T. & Karunaratne S. (2000) Dynamic modeling of the growth of Phragmites australis: model description. Aquat. Bot. 67, 301-318. Balzter H., Braun P.W. & Köhler W. (1998) Cellular automata models for vegetation dynamics. Ecol. Mod. 107, 113-125. Barko J.W. & James W.F. (1997) Effects of submerged aquatic macrophytes on nutrient dynamics, sedimentation, and resuspension. In: Jeppesen E., Søndergaard M. & Christoffersen K. (eds.) The Structuring Role of Submerged Macrophytes in Lakes. Springer, New York, 197-214. Barko J.W. & Smart R.M. (1979) The nutritional ecology of Cyperus esculentus, an emergent aquatic plant, grown on different sediments. Aquat. Bot. 6, 13-28.

95 References

Barko J.W. & Smart R.M. (1980) Mobilization of sediment phosphorus by submersed freshwater macrophytes. Freshw. Biol. 10, 229-238. Barko J.W. & Smart R.M. (1981a) Comparative influences of light and temperature on the growth and metabolism of selected submersed freshwater macrophytes. Ecol. Monogr. 51, 219-236. Barko J.W. & Smart R.M. (1981b) Sediment-based nutrition of submersed macrophytes. Aquat. Bot. 10, 339-352. Barko J.W. & Smart R.M. (1983) Effects of organic matter additions to sediment on the growth of aquatic plants. J. Ecol. 71, 161-176. Barko J.W. & Smart R.M. (1986a) Sediment-related mechanisms of growth limitation in submersed macrophytes. Ecology 67, 1328-1340. Barko J.W. & Smart R.M. (1986b) Effects of sediment composition on growth of submersed aquatic vegetation. Tech. Rep. U.S.Army Eng.Waterways Exp. Stn. Barko J.W., Hardin D.G. & Matthews M.S. (1982) Growth and morphology of submersed freshwater macrophytes in relation to light and temperature. Can. J. Bot. 60, 877-887. Barko J.W., Smart R.M. & McFarland D.G. (1991) Interactive effects of environmental conditions on the growth of submersed aquatic macrophytes. J. Freshw. Ecol. 6, 199-208. Barko J.W., Smart R.M., Mcfarland D.G. & Chen R.L. (1988) Interrelationships between the growth of verticillata (L.f.) Royle and sediment nutrient availability. Aquat. Bot. 32, 205-216. Bazzaz F.A. & Grace J. (eds.) (1997) Plant Resource Allocation. Academic Press, San Diego, London, Boston, N.Y., Sydney, Tokyo, Toronto. Beklioglu M. & Moss B. (1996) Mesocosm experiments on the interaction of sediment influence, fish predation and aquatic plants with the structure of phytoplankton and zooplankton communities. Freshw. Biol. 36, 315- 325. Bell A.D. (1976) Computerized vegetative mobility in rhizomatous plants. In: Lindenmayer A. & Rozenberg G. (eds.) Automata, Languages, Development. North-Holland, Amsterdam, 3-14. Bell A.D. (1984) Dynamic morphology: a contribution to plant population ecology. In: Dirzo R. & Sarukhan J. (eds.) Perspectives on Plant Population Ecology. Sinauer Associates, Sunderland, 48-65. Bell A.D., Roberts D. & Smith A. (1979) Branching patterns: Simulation of plant architecture. J. Theor. Biol. 81, 351-375. Best M.D. & Mantai K.E. (1978) Growth of Myriophyllum: Sediment or lake water as the source of nitrogen and phosphorous. Ecology 59, 1075-1080. Best E.-P.H., Buzzelli C.P., Doyle R.D., Wetzel R.L., Boyd W.A. & Campbell K.R. (2001) Modeling submersed macrophyte growth in relation to underwater light climate: modeling approaches and application potential. Hydrobiologia 4441, 43-70. Bini L.M., Thomaz S.M., Murphy K.J. & Camargo A.F.M. (1999) Aquatic macrophyte distribution in relation to water and sediment conditions in the Itaipu Reservoir, Brazil. Hydrobiologia 415, 147-154. Black M.A., Maberly S.C. & Spence D.H.N. (1981) Resistances to carbon dioxide fixation in four submerged freshwater macrophytes. New Phytol. 89, 557-568. Blindow I. (1987) The composition and density of epiphyton on several species of submerged macrophytes - the neutral substrate hypothesis tested. Aquat. Bot. 29, 157-168. Bloom A.J., Chapin F.S. & Mooney H.A. (1985) Resource limitation in plants: an economic analogy. Ann. Rev. Ecol. Syst. 16, 363–392. Bonis A., Lepart J. & Grillas P. (1995) Seed bank dynamics and coexistence of annual macrophytes in a temporary and variable habitat. Oikos 74, 81-92. Briske D.D. & Butler J.L. (1989) Density-dependent regulation of ramet populations within the bunchgrass Schizachyrium scoparium interclonal versus intraclonal interference. J. Ecol. 77, 963-974. Brux H., Todeskino D. & Wiegleb G. (1987) Growth and reproduction of Potamogeton alpinus Balbis growing in disturbed habitats. Arch. Hydrobiol. Beih. Ergebn. Limnol. 27, 115-127.

96 References

Caffrey J.M. (1991) Seasonal and spatial patterns of oxygen production, respiration and root-rhizome release in Potamogeton perfoliatus L. and Zostera marina L. Aquat. Bot. 109, 128. Cain M.L. (1994) Consequences of foraging in clonal plant species. Ecology 75, 933-944. Cain M.L. & Cook R.E. (1988) Growth in Medeola virginiana clones. II Stochastic simulation of vegetative spread. Am. J. Bot. 75, 730-736. Cain M.L., Dudle D.A. & Evans J.P. (1996) Spatial models of foraging in clonal plant species. Am. J. Bot. 83, 76- 85. Cain M.L., Pacala S.W., Silander J.A.Jr. & Fortin M.-J. (1995) Neighborhood models of clonal growth in the white clover Trifolium repens. Am. Nat. 145, 888-917. Calado G. & Duarte P. (2000) Modelling growth of cirrhosa. Aquat. Bot. 68, 29-44. Callaghan T.V., Svensson B.M., Bowman H., Lindley D.K. & Carlsson B.Å. (1990) Models of clonal plant growth based on population of dynamics and architecture. Oikos 57, 257-269. Campbell B.D., Grime J.P. & Mackey J.M.L. (1991) A trade-off between scale and precision in resource foraging. Oecologia 87, 532-538. Canfield D.E.jr., Langeland K.A., Linda S.B. & Haller W.T. (1985) Relations between water transparency and maximum depth of macrophyte colonization in lakes. J. Aquat. Plant Manage. 23, 25-28. Carbiener R., Tremolieres M. & Muller S. (1995) Vegetation of running waters and water quality: Thesis, debates and prospects. Acta Bot. Gall. 142, 489-531. Carignan R. (1985) Nutrient dynamics in a littoral sediment colonized by the submersed macrophyte Myriophyllum spicatum. Can. J. Fish. Aquat. Sci. 42, 1303-1311. Carignan R., Pierre S.St. & Gachter R. (1994) Use of diffusion samplers in oligotrophic lake sediments: Effects of free oxygen in sampler material. Limnol. Oceanogr. 39, 468-474. Carr G.M., Duthie H.C. & Taylor W.D. (1997) Models of aquatic plant productivity: A review of the factors that influence growth. Aquat. Bot. 59, 195–215. Casper S.J. & Krausch H.-D. (1980) Süßwasserflora von Mitteleuropa. 1. Teil: Pteridohyta und Anthophyta. Gustav Fischer Verlag, Stuttgart. Cebrián J. & Duarte C.M. (2001) Detrital stocks and dynamics of the seagrass Posidonia oceanica (L.) Delile in the Spanish Mediterranean. Aquat. Bot. 70, 295-309. Chambers P.A. (1987) Light and nutrients in the control of aquatic plant community structure II. In-situ observations. J. Ecol. 75, 621-628. Chambers P.A. & Fourqurean J.W. (1991) Alternative criteria for assessing nutrient limitation of a wetland macrophyte Peltandra virginica L. KUNTH. Aquat. Bot. 40, 305-320. Chambers P.A. & Kalff J. (1985) The influence of sediment composition and irradiance on the growth and morphology of Myriophyllum spicatum. Aquat. Bot. 22, 253-264. Chambers P.A. & Prepas E.E. (1990) Competition and coexistence in submerged aquatic plant communities: The effects of species interactions versus abiotic factors. Freshw. Biol. 23, 541-550. Chambers P.A., Prepas E.E., Hamilton.H.R. & Bothwell M.L. (1991) Current velocity and its effect on aquatic macrophytes in flowing waters. Ecol. Appl. 1, 249-257. Chen D.X. & Coughenour M.B. (1996) A mechanistic model for submerged aquatic macrophyte photosynthesis:

Hydrilla in ambient and elevated CO2. Ecol. Mod. 89, 133-146. Cheplick G.P. (1995) Genotypic variation and plasticity of clonal growth in relation to nutrient availability in Amphibromus scabrivalvis. J. Ecol. 83, 459-468. Chesson P. & Peterson A.G. (2002) The quantitative assessment of the benefits of physiological integration in clonal plants. Evol. Ecol. Res. 4, 1153-1176. Chiarello E. & Barrat-Segretain M.-H. (1997) Recolonization of cleared patches by macrophytes: Modelling with point processes and random mosaics. Ecol. Mod. 96, 61-73.

97 References

Christiansen R., Friis N.J.S. & Søndergaard M. (1985) Leaf production and nitrogen and phosphorus tissue content of Littorella uniflora in relation to nitrogen and phosphorus enrichment of the sediment in oligotrophic Lake Hampen Denmark. Aquat. Bot. 23, 1-12. Cízková H., Pechar L., Husák S., Kvet J., Bauer V., Radová J. & Edwards K. (2001) Chemical characteristics of soils and pore waters of three wetland sites dominated by Phragmites australis: relation to vegetation composition and reed performance. Aquat. Bot. 69, 235-249. Clarke S.J. & Wharton G. (2001) Sediment nutrient characteristics and aquatic macrophytes in lowland English rivers. Science Total Env. 266, 103-112. Clevering O.A. & Hundscheid M.P.J. (1998) Plastic and non-plastic variation in growth of newly established clones of Scirpus (Bolboschoenus) maritimus L. grown at different water depths. Aquat. Bot. 62, 1-17. Colasanti R.L. & Grime J.P. (1993) Resource dynamics and vegetation processes: a deterministic model using two-dimensional cellular automata. Funct. Ecol. 7, 169-176. Coleman J.S & McConnaughay K.D.M. (1995) A non-functional interpretation of a classical optimal-partitioning example. Funct. Ecol. 9, 951-954. Collins C.D. & Wlosinski J.H. (1989) A macrophyte submodel for aquatic ecosystems. Aquat. Bot. 33, 191-206. Collins C.D., Sheldon R.B. & Boylen C.W. (1987) Littoral zone macrophyte community structure: distribution and association of species along physical gradients in Lake George, New York, U.S.A. Aquat. Bot. 29, 177-194. Coops H., Boeters R. & Smit H. (1991) Direct and indirect effects of wave attack on helophytes. Aquat. Bot. 41, 333-352. Daldorph P.-W.G. & Thomas J.D. (1991) The effect of nutrient enrichment on a freshwater community dominated by macrophytes and molluscs and its relevance to snail control. J. Appl. Ecol. 28, 685-702. Dale H.M. (1986) Temperature and light: The determining factors in maximum depth distribution of aquatic macrophytes in Ontario, Canada. Hydrobiologia 133, 73-77. Dale M.R.T. (1999) Spatial Pattern Analysis in Plant Ecology. Cambridge University Press, Cambridge. Davis B.C. & Fourqurean J.W. (2001) Competition between the tropical alga, Halimeda incrassata, and the seagrass, Thalassia testudinum. Aquat. Bot. 71, 217-232. Davis J.F. & McDonnell A.J. (1997) Development of a partitioned-biomass model for rooted macrophyte growth. Aquat. Bot. 56, 265-276. De Kroon H. (1993) Competition between shoots in stands of clonal plants. Plant Spec. Biol. 8, 85-94. De Kroon H. & Hutchings M.J. (1995) Morphological plasticity in clonal plants: The foraging concept reconsidered. J. Ecol. 83, 143-152. De Kroon H. & Knops J. (1990) Habitat exploration through morphological plasticity in two chalk grassland perennials. Oikos 59, 39-49. De Kroon H. & Kwant R. (1991) Density-dependent growth responses in two clonal herbs: regulation of shoot density. Oecologia (Berlin) 86, 298-304. De Kroon H. & Schieving F. (1990) Resource partitioning in relation to clonal growth strategy. In: Van Groenendael J. & De Kroon H. (eds.) Clonal Growth in Plants: Regulation and Function. SPB Academic, The Hague, 113-130. De Kroon H., Huber H., Stuefer J.F. & Van Groenendael J.M. (2005) A modular concept of phenotypic plasticity in plants. New Phytol. 166, 73-82. De Vos A.P.C. (1954) The littoral fauna of Lake IJssel [Dutch]. In: De Beaufort L.F. (ed.) Veranderingen in de flora en fauna van de Zuiderzee (thans IJsselmeer) na de afsluiting in 1932. De Boer, Den Helder, 268-276. Demars B.O.L. & Harper J.K. (1998) The aquatic macrophytes of an English lowland river system: Assessing response to nutrient enrichment. Hydrobiologia 384, 75-88. Dennis B. (2002) Allee effects in stochastic populations. Oikos 96, 389-401. Denny P. (1980) Solute movement in submerged angiosperms. Biol. Rev. 55, 65-92.

98 References

DLR (2001) Aerial image, Deutsches Zentrum für Luft- und Raumfahrt, Oberpfaffenhofen, Germany. Dong M. & De Kroon H. (1994) Plasticity in morphology and biomass allocation in Cynodon dactylon, a grass species forming stolons and rhizomes. Oikos 70, 99-106. Dong M., During H.J. & Werger M.J.A. (1997) Clonal plasticity in response to nutrient availability in the pseudoannual herb, Trientalis europaea L. Plant Ecol. 131, 233-239. Downing J.A. & Anderson M.R. (1985) Estimating the standing biomass of aquatic macrophytes. Can. J. Fish. Aquat. Sci. 42, 1860-1869. Downing J.A. & Rath L.C. (1988) Spatial patchiness in the lacustrine sedimentary environment. Limnol. Oceanogr. 33, 447-458. Duarte C.M. (1991) Allometric scaling of seagrass form and productivity. Mar. Ecol. Prog. Ser. 77, 289-300. Duarte C.M. (1995) Submerged aquatic vegetation in relation to different nutrient regimes. Ophelia 41, 87-112. Duarte C.M. & Kalff J. (1986) Littoral slope as a predictor of the maximum biomass of submerged macrophyte communities. Limnol. Oceanogr. 31, 1072-1080. Duarte C.M. & Kalff J. (1990a) Biomass density and the relationship between submerged macrophyte biomass and plant growth form. Hydrobiologia 196, 17-24. Duarte C.M. & Kalff J. (1990b) Patterns in the submerged macrophyte biomass of lakes and the importance of the scale of analysis in the interpretation. Can. J. Fish. Aquat. Sci. 47, 357-363. Duarte C.M. & Sand-Jensen K. (1990) Seagrass colonization: Biomass development and shoot demography in Cymodocea nodosa patches. Mar. Ecol. Prog. Ser. 67, 97-103. Duarte C.M., Hemminga M.A. & Marbà N. (1996) Growth and population dynamics of Thalassodendron ciliatum in a Kenyan back-reef lagoon. Aquat. Bot. 55, 1-11. Duarte C.M., Merino M., Agawin N.S.R., Uri J., Fortes M.D., Gallegos M.E., Marbà N. & Hemminga M.A. (1998) Root production and belowground seagrass biomass. Mar. Ecol. Prog. Ser. 171, 97-108. Dvorak J. & Best E.P.H. (1982) Macro invertebrate communities associated with the macrophytes of Lake Vechten The Netherlands structural and functional relationships. Hydrobiologia 95, 115-126. Ellison A.M. & Niklas K.J. (1987) Branching patterns of Salicornia europaea (Chenopodiaceae) at different successional stages: a comparison of theoretical and real plants. Am. J. Bot. 75, 501-512. Erftemeijer P.L.A., Stapel J., Smekens M.J.E. & Drossaert W.M.E. (1994). The limited effect of in situ phosphorus and nitrogen additions to seagrass beds on carbonate and terrigenous sediments in South Sulawesi, Indonesia. J. Exp. Mar. Biol. Ecol. 182, 123-140. Ervin G.N. & Wetzel R. (1997) Shoot:root dynamics during growth stages of the rush Juncus effusus L. Aquat. Bot. 59, 63-73. Evans J.P. (1988) Nitrogen translocation in a clonal dune perennial Hydrocotyle bonariensis. Oecologia 77, 64- 68. Evans J.P. (1991) The effect of resource integration on fitness related traits in a clonal dune perennial Hydrocotyle bonariensis. Oecologia 86, 268-275. Evans J.P. & Cain M.L. (1995) A spatially explicit test of foraging behavior in a clonal plant. Ecology 76, 1147- 1155. Evans J.P. & Whitney S. (1992) Clonal integration across a salt gradient by a nonhalophyte, Hydrocotyle bonariensis (Apiaceae). Am. J. Bot. 79, 1344-1347. Fischer M., Van Kleunen M. & Schmid B. (2000) Genetic Allee effects on performance, plasticity and developmental stability in a clonal plant. Ecol. Letters 3, 530-539. Fonseca M.S. & Bell S.S. (1998) Influence of physical setting on seagrass landscapes near Beaufort, North Carolina, USA. Mar. Ecol. Prog. Ser. 171, 109-121. Fonseca M.S. (1989) Sediment stabilization by Halophila decipiens in comparison to other seagrasses. Est. Coast. Shelf Sci. 29, 501-507.

99 References

Fonseca M.S. (1996) Scale dependence in the study of seagrass systems. In: Kuo J., Phillips R.C., Walker D.I. & Kirkman H. (eds.) Seagrass Biology: Proceedings of an International Workshop, Rottnest Island, Western Australia, 25 - 29 Jan. 1996, Sciences UWA, Nedlands,. 95-104. Fonseca M.S., Zieman J.C., Thayer G.W. & Fisher J.S. (1983) The role of current velocity in structuring eelgrass (Zostera marina L) meadows. Est. Coast. Shelf Sci. 17, 367-380. Fourqurean J.W., Zieman J.C. & Powell G.V.N. (1992) Relationships between porewater nutrients and seagrasses in a subtropical carbonate environment. Mar. Biol. (Berlin) 114, 57-65. Gafny S. & Gasith A. (1991) Spatially and temporally sporadic appearance of macrophytes in the littoral zone of Lake Kinneret, Israel taking advantage of a window of opportunity. Aquat. Bot. 62, 249-267. Gardner J.L. (1998) Experimental evidence for edge-related predation in a fragmented agricultural landscape. Aust. J. Ecol. 23, 311-321. Gardner S.N. & Mangel M. (1999) Modeling investments in seeds, clonal offspring, and translocation in a clonal plant. Ecology 80, 1202-1220. Gedroc J.J., McConnaughay K.D.M. & Coleman J.S. (1996) Plasticity in root/shoot partitioning: Optimal, ontogenetic, or both? Funct. Ecol. 10, 44-50. Gerloff G.C. & Krombholz P.H. (1966) Tissue analysis as a measure of nutrient availability for the growth of angiosperm aquatic plants. Limnol. Oceanogr. 11, 529-537. Golterman H.L. (2004) The Chemistry of Phosphate and Nitrogen Compounds in Sediments. Kluwer, Dordrecht. Goudriaan J. (1986) A simple and fast numerical method for the computation of daily totals of crop photosynthesis. Agric. Forest Meteorol. 38, 249-254. Grace J.B. (1993) The adaptive significance of clonal reproduction in angiosperms: an aquatic perspective. Aquat. Bot. 44, 159-180. Granéli W., Weisner S.E.B. & Sytsma M.D. (1992) Rhizome dynamics and resource storage in Phragmites australis. Wetland Ecol. Manage. 1, 239-247. Grime J.P. (1979) Plant Strategies and Vegetation Processes. John Wiley, Chichester. Grimm V. (2002) Visual debugging: a way of analyzing, understanding and communicating bottom-up simulation models in ecology. Nat. Res. Mod. 15, 23-38. Groom M. (1998) Allee effects limit population viability on an annual plant. Am. Nat. 151, 487-496. Groth A.T., Lovett-Doust L. & Lovett-Doust J. (1996) Population density and module demography in Trapa natans (Trapaceae), an annual, clonal aquatic macrophyte. Am. J. Bot. 83, 1406-1415. Gusak S. (1983) Contribution to the biology and coenology of Potamogeton perfoliatus L. Hydrobiol. J. 19, 54-58. Hangelbroek H.H., Ouborg N.J., Santamaría L., Schwenk K. (2002) Clonal diversity and structure within a population of the pondweed Potamogeton pectinatus foraged by Bewick's swans. Molec. Ecol. 11, 2137-2150. Hartnett D.C. & Bazzaz F.A. (1983) Physiological integration among intraclonal ramets in Solidago canadensis. Ecology 64, 779-788. Hartnett D.C. & Bazzaz F.A. (1985) The regulation of leaf, ramet and genet densities in experimental populations of the rhizomatous perennial Solidago canadensis. J. Ecol. 73, 429-443. Heege T., Bogner A. & Pinnel N. (2003) Mapping of submerged aquatic vegetation with a physically based processing chain. In: Kramer E. (ed.) SPIE-The International Society for Optical Engineering 5233, Proc. Remote Sensing 2003, September 8-12 2003, Barcelona. Hellström K., Kytöviita H.-M., Tuomi J. & Rautio P. (2006) Plasticity of clonal integration in the perennial herb Linaria vulgaris after damage. Funct. Ecol. 20, 413-420. Herben T. & Suzuki J. (2002) A simulation study of the effects of architectural constraints and resource translocation on population structure and competition in clonal plants. Evol. Ecol. 15, 403-423. Hester M.W., McKee K.L., Burdick D.M., Koch M.S., Flynn K.M., Patterson S. & Mendelssohn I.A. (1994) Clonal integration of Spartina patens across a nitrogen and salinity gradient. Can. J. Bot. 72, 767-770.

100 References

Hilbert D.W. (1990) Optimization of plant root:shoot ratios and internal nitrogen concentration. Ann. Bot. 66, 91- 99. Hirose T. (1986) Nitrogen uptake and plant growth. II. An empirical model of vegetative growth and partitioning. Ann. Bot. 58, 487-496. Hirose T. (1987) A vegetative plant growth model: Adaptive significance of phenotypic plasticity in matter partitioning. Funct. Ecol. 1, 195-202. Hofstra D.E., Adam K.D. & Clayton J.S. (1995) Isozyme variation in New Zealand populations of Myriophyllum and Potamogeton species. Aquat. Bot. 52, 121-131. Honkanen T. & Haukioja E. (1994) Why does a branch suffer more after branch-wide than after tree-wide defoliation? Oikos 71, 441-450. Hootsmans M.J.M. (1999) Modelling Potamogeton pectinatus: For better or for worse. Hydrobiologia 415, 7-11. Hootsmans M.J.M. (1994) A growth analysis model for Potamogeton pectinatus L. In: Van Vierssen W., Hootsmans M. & Vermaat J. (eds.) Lake Veluwe. A Macrophyte-Dominated System under Eutrophication Stress, Kluwer Academic Publishers, Dordrecht, 250-286. Huber H. & Stuefer J.F. (1997) Shade-induced changes in the branching pattern of a stoloniferous herb: functional response or allometric effect? Oecologia 110, 478-486. Huber H., Lukács S. & Watson M.A. (1999) Spatial structure of stoloniferous herbs: an interplay between structural blue-print, ontogeny and phenotypic plasticity. Plant Ecol. 141, 107-115. Husák S., Sladecek V. & Sladeckova A. (1989) Freshwater macrophytes as indicators of organic pollution. Acta hydrochim. hydrobiol. 17, 731-735. Hutchings M. (1979) Weight-density relationships in ramet populations of clonal perennial herbs, with special reference to the -3/2 power law. J. Ecol. 67, 21-33. Hutchings M.J. (1988) Differential foraging for resources and structural plasticity in plants. Trends Ecol. Evol. 3, 200-204. Hutchings M.J. (1997) Resource allocation patterns in clonal herbs and their consequences for growth. In: Bazzaz F.A. & Grace J. (eds.) Plant Resource Allocation. Academic Press, San Diego, 161-189. Hutchings M. & Mogie M. (1990) The spatial structure of clonal plants: control and consequences. In: Van Groenendael J. & De Kroon H. (eds.) Clonal Growth in Plants: Regeneration and Function. SPB Academic Publishing, The Hague, 57-76. Hutchings M.J. & De Kroon H. (1994) Foraging in plants: the role of morphological plasticity in resource acquisition. Adv. Ecol. Res. 25, 159-238. Hutchings M.J. & Wijesinghe D.K. (1997) Patchy habitats, division of labour and growth dividends in clonal plants. Trends Ecol. Evol. 12, 390-394. Idestam-Almquist J. & Kautsky L. (1995) Plastic responses in morphology of Potamogeton pectinatus L. to sediment and above-sediment conditions at two sites in the northern Baltic proper. Aquat. Bot. 52, 205-216. IGKB (2004) Der Bodensee. Zustand - Fakten - Perspektiven. IGKB (Internationale Gewässerschutzkommission für den Bodensee), Bregenz. Ikusima I. (1970) Ecological studies on the productivity of aquatic plant communities IV. Light condition and community photosynthetic production. Bot. Magaz. Tokyo 83, 330-341. Jackson J.B.C., Buss L.W. & Cook R.E. (1985) Population Biology and Evolution of Clonal Organisms. Yale University Press. Jackson S.T. & Charles D.F. (1988) Aquatic macrophytes in Adirondack (New York) lakes: Patterns of species composition in relation to environment. Can. J. Bot. 66, 1449-1460. Jaynes M.L. & Carpenter S.R. (1986) Effects of vascular and nonvascular macrophytes on sediment redox and solute dynamics. Ecology 67, 875-882. Jeppesen E., Søndergaard M. & Christoffersen K. (1997) The Structuring Role of Submerged Macrophytes in Lakes. Ecol. Studies 131. Springer, New York.

101 References

Joehnk K.D., Straile D. & Ostendorp W. (2004) Water level variability and trends in Lake Constance in the light of the 1999 centennial flood. Limnologica 34, 15-21. Jones J.I., Eaton J.W. & Hardwick K. (2000) The influence of periphyton on boundary layer conditions: a pH microelectrode investigation. Aquat. Bot. 67, 191-206. Jonzén N., Nolet B.A., Santamaría L. & Svensson M.G.E. (2002) Seasonal herbivory and mortality compensation in a swan-pondweed system. Ecol. Mod. 147, 209-219. Jónsdottir I.S. & Callaghan T.V. (1990) Intraclonal translocation of ammonium and nitrate nitrogen in Carex bigelowii Torr. ex Schwein. using nitrogen-15 and nitrate reductase assays. New Phytol. 114, 419-428. Karjalainen H., Stefansdottir G., Tuominen L. & Kairesalo T. (2001) Do submersed plants enhance microbial activity in sediment? Aquat. Bot. 69, 1-13. Kautsky L. (1987) Life-cycles of three populations of Potamogeton pectinatus L. at different degrees of wave exposure in the asko area northern Baltic proper. Aquat. Bot. 27, 177-186. Kautsky L. (1991) In situ experiments on interrelationships between six brackish macrophyte species. Aquat. Bot. 39, 159-172. Keddy P.A. (1982) Quantifying within lake gradients of wave energy: interrelationships of wave energy, substrate particle size and shoreline plants in Axe Lake, Ontario. Aquat. Bot. 14, 41-58. Keddy P.A. (1983) Shoreline vegetation in Axe Lake, Ontario: effects of exposure on zonation patterns. Ecology 64, 331-344. Kelly C.K. (1990) Plant foraging a marginal value model and coiling response in Cuscuta subinclusa. Ecology 71, 1916-1925. Klepper O. (1989) A model of carbon flows in relation to macrobenthic food supply in the Oosterschelde estuary (S.W. Netherlands). PhD Thesis, Wageningen University, The Netherlands. Klepper O., Van der Tol M.W., Scholten H. & Herman P.M. (1994) Smoes - A simulation model for the Oosterschelde ecosystem.1.Description and uncertainty analysis. Hydrobiologia 283, 437-451. Klimes L. (2000) Phragmites australis at an extreme altitude: rhizome architecture and its modelling. Fol. Geobot. Phytotaxon. 35, 403-417. Klimes L., Klimesova J., Hendriks R. & Van Groenendael J. (1997) Clonal plant architecture: a comparative analysis of form and function. In: De Kroon H. & Van Groenendael J. (eds.) The Ecology and Evolution of Clonal Plants. Backhuys Publishers, Leiden, 1-29. Koch E.W. (1999) Preliminary evidence on the interdependent effects of currents and porewater geochemistry on Thalassia testudinum Banks ex König seedlings. Aquat. Bot. 63, 95-102. Koch E.W. (1996) Hydrodynamics of a shallow Thalsassia testudinum bed in Florida, USA. In: Kuo J., Phillips R.C., Walker D.I. & Kirkman H. (eds.) Seagrass Biology: Proceedings of an International Workshop, Rottnest Island, Western Australia, 25 - 29 Jan. 1996, Sciences UWA, Nedlands, 105-110. Lachavanne J.-B. (1985) The influence of accelerated eutrophication on the macrophytes of Swiss lakes: Abundance and distribution. Verh. Int. Verein. Limnol. 22, 2950-2955. Lachavanne J.-B. (1991) The consequences of water oligotrophication on macrophytic vegetation of Swiss lakes. Verh. Int. Verein. Limnol. 24, 943-948. Langeland K.A., Sutton D.L., Canfield D.E.jr. (1983) Growth response of Hydrilla to extractable nutrients in prepared substrates. J. Freshw. Ecol. 2, 263-272. Lauridsen T.L., Pedersen L.J., Jeppesen E. & Søndergaard M. (1996) The importance of macrophyte bed size for cladoceran composition and horizontal migration in a shallow lake. J. Plankton Res. 18, 2283-2294. Lee S.Y., Fong C.W. & Wu R.S.S. (2001) The effects of seagrass (Zostera japonica) canopy structure on associated fauna: a study using artificial seagrass units and sampling of natural beds. J. Exp. Mar. Biol. Ecol. 259, 23-50. Lehmann A. (1997) Apports des systèmes d'information géographique à l'étude des plantes aquatiques submergées en milieu lacustre. PhD Thesis, University of Geneva, Geneva.

102 References

Lehmann A., Castella E. & Lachavanne J.-B. (1997) Morphological traits and spatial heterogeneity of aquatic plants along sediment and depth gradients, Lake Geneva, Switzerland. Aquat. Bot. 55, 281-299. Levin P.S. & Mathieson A.C. (1991) Variation in a host-epiphyte relationship along a wave-exposure gradient. Mar. Ecol. Prog. Ser. 77, 271-278. Losee R.F. & Wetzel R.G. (1993) Littoral flow rates within and around submersed macrophyte communities. Freshw. Biol. 29, 7-17. Lovett-Doust L. (1981) Population dynamics and local specialization in a clonal perennial (Ranunculus repens) I. The dynamics of ramets in contrasting habitats. J. Ecol. 69, 743-755. Lundholm J.T. & Simser W.L. (1999) Regeneration of submerged macrophyte populations in a disturbed lake Ontario Coastal marsh. J. Great Lakes Res. 25, 395-400. Luther H. (1951) Verbreitung und Okologie der hoheren Wasserpflanzen im Brackwasser der Ekenäs-Gegend in Südfinnland. II. Spezieller Teil. Acta Bot. Fennica 50, 1-350. Lynch M., Conery J. & Bürger R. (1995) Mutation accumulation and the extinction of small populations. Am. Nat. 146, 489-518. Madsen J.D. (1991) Resource allocation at the individual plant level. Aquat. Bot. 41, 67-86. Madsen J.D. & Adams M.S. (1989) The light and temperature dependence of photosynthesis and respiration in Potamogeton pectinatus L. Aquat. Bot. 36, 23-32. Madsen T.V. & Warncke E. (1983) Velocities of currents around and within submersed aquatic vegetation. Arch. Hydrobiol. 97, 389-394. Mallik A.U. & Rasid H. (1993) Root-shoot characteristics of riparian plants in a flood control channel: implications for bank stabilization. Ecol. Eng. 2, 149-158. Mantai K.E. & Newton M.E. (1982) Root growth in Myriophyllum: A specific plant response to nutrient availability? Aquat. Bot. 13, 45-55. Marbà N. & Duarte C.M. (1998) Rhizome elongation and seagrass clonal growth. Mar. Ecol. Prog. Ser. 174, 269- 280. Marbà N., Hemminga M.A., Mateo M.A., Duarte M.D., Mass Y.E.M., Terrados J. & Gacia E. (2002) Carbon and nitrogen translocation between seagrass ramets. Mar. Ecol. Prog. Ser. 226, 287-300. Marbà N., Cebrián J., Enríquez S., Duarte C.M. (1996) Growth patterns of Western Mediterranean seagrasses: species-specific responses to seasonal forcing. Mar. Ecol. Prog. Ser. 133, 203-215. Marshall C. & Price E.A.C. (1997) Sectoriality and its implications for physiological integration. In: De Kroon H. & Van Groenendael J. (eds.) The Ecology and Evolution of Clonal Plants. Backhuys Publishers, Leiden, 79-107. Marshall C. & Sagar G.R. (1965) The influence of defoliation on the distribution of assimilates in Lolium multiflorum Lam. Ann. Bot. 29, 365-370. McComb A.J., Qiu S., Lukatelich R.J. & McAuliffe T.F. (1998) Spatial and temporal heterogeneity of sediment phosphorus in the Peel Harvey Estuarine System. Est. Coast.Shelf Sci. 47, 561-577. Meijer M.-L., de Bois I., Scheffer M., Portielje R. & Hosper H. (1999) Biomanipulation in shallow lakes in The Netherlands: An evaluation of 18 case studies. Hydrobiologia 408-409, 13-30. Melzer A. & Greiner K. (1985) Die Beeinflussung der Netto-Photosyntheserate makrophytischer Wasserpflanzen durch Außenfaktoren. In: Melzer A. & Hamm A. (eds.) Die Bedeutung von Makrophyten für die Gewässerökologie 1, Bayerische Landesanstalt für Wasserforschung, München, 148-164. Melzer A., Held K., Harlacher R. & Vogt E. (1988) Die Makrophytenvegetation des Ammer-, Wörth- und Pilsensees sowie des Weißlinger Sees. Informationsberichte 1, Bayrisches Landesamt für Wasserwirtschaft, München. Methy M., Alpert P. & Roy J. (1990) Effects of light quality and quantity on growth of the clonal plant Eichhornia- crassipes. Oecologia 84, 265-271. Miller B. (1999) Clonal Plant Foraging: a Review and Simulation Study. B.A. Thesis, Reed College, Portland, USA, pp. 105 (digital).

103 References

Moen R.A. & Cohen Y. (1989) Growth and competition between Potamogeton pectinatus L. and Myriophyllum exalbescens Fern. in experimental ecosystems. Aquat. Bot. 33, 257-270. Molenaar H., Barthélémy D., De Reffye P., Meinesz A. & Mialet I. (2000) Modelling architecture and growth patterns of Posidonia oceanica. Aquat. Bot. 66, 85-99. Mooney H.A., Küppers M., Koch G.W., Gorham J., Chu C.C. & Winner W.E. (1988) Compensating effects to

growth of carbon partitioning changes in response to SO2-induced photosynthetic reduction in radish. Oecologia 72, 502-506. Muhammetoglu A.B. & Soyupak S. (2000) A three-dimensional water quality-macrophyte interaction model for shallow lakes. Ecol. Mod. 133, 161-180. Müller G. (1977) Schadstoff-Untersuchungen an datierten Sedimentkernen aus dem Bodensee. Z. Naturforsch. 32, 920-925. Müller I., Schmid B. & Weiner J. (2000) The effect of nutrient availability on biomass allocation patterns in 27 species of herbaceous plants. Perspect. Plant Ecol. Evol. Syst. 3, 115-127. Nichols D.S. & Keeney D.R. (1976) Nitrogen nutrition of Myriophyllum spicatum: variation of plant tissue nitrogen concentration with season and site in Lake Wingra. Freshw. Biol. 6, 137-144. Nielsen S.L. & Pedersen M.F. (2000) Growth, photosynthesis and nutrient content of seedlings and mature plants of Cymodocea nodosa - the importance of clonal integration. Aquat. Bot. 68, 265-271. Novoplansky A., Cohen D. & Sachs T. (1989) Ecological implications of correlative inhibition between plant shoots. Physiol. Plant. 77, 136-140. Oborny B. (1991) Criticisms on optimal foraging in plants: a review. Abstr. Bot. 15, 67-76. Oborny B. (1994) Growth rules in clonal plants and environmental predictability - a simulation study. J. Ecol. 82, 341-351. Oborny B. & Cain M.L. (1997) Models of spatial spread and foraging in clonal plants. In: De Kroon H. & Van Groenendael J. (eds.) The Ecology and Evolution of Clonal Plants. Backhuys Publishers, Leiden, 155-183. Oborny B. & Kun Á. (2002) Fragmentation of clones: how does it influence dispersal and competitive ability? Evol. Ecol. 15, 319-346. Oborny B., Kun Á., Tamás C. & Szilárd B. (2000) The effect of clonal integration on plant competition for mosaic habitat space. Ecology 81, 3291-3304. Olesen B. & Sand-Jensen K. (1994a) Demography of shallow eelgrass (Zostera marina) populations: Shoot dynamics and biomass development. J. Ecol. 82, 379-390. Olesen B. & Sand-Jensen K. (1994b) Patch dynamics of eelgrass Zostera marina. Mar. Ecol. Prog. Ser. 106, 147- 156. Olson B.E. & Wallander R.T. (1999) Carbon allocation in Euphorbia esula and neighbours after defoliation. Can. J. Bot. 77, 1641-1647. Ong C.K. & Marshall C. (1979) The growth and survival of severly-shaded tillers in Lolium perenne L. Ann. Bot. 43, 147-155. Ozimek T., Prejs A. & Prejs K. (1976) Biomasss and distribution of underground parts of Potamogeton perfoliatus L. and P. lucens L. in Mikolajskie lake, Poland. Aquat. Bot. 2, 309-316. Ozimek T., Pieczynska E. & Hankiewicz A. (1991) Effects of filamentous algae on submerged macrophyte growth: a laboratory experiment. Aquat. Bot. 41, 309-315. Pagès L., Asseng S., Pellerin S. & Diggle A. (2000) Modelling root system growth and architecture. In: Smit A.L., Bengough A.G., Engels C., Van Noordwijk M., Pellerin S. & van de Geijn S.C. (eds.) Root Methods. A Handbook. Springer, Berlin, Heidelberg, 113-146. Paling E.I. & McComb A.J. (2000) Autumn biomass, below-ground productivity, rhizome growth at bed edge and nitrogen content in seagrasses from Western Australia. Aquat. Bot. 67, 207-219. Peralta G., Bouma T. J., van Soelen J., Pérez-Lloréns J. L. & Hernández I. (2003) On the use of sediment fertilization for seagrass restoration: a mesocosm study on Zostera marina L. Aquat. Bot. 75, 95-110.

104 References

Pérez M., Duarte C.M., Romero J., Sand-Jensen K., Alcoverro T. (1994) Growth plasticity in Cymodocea nodosa stands: the importance of nutrient supply. Aquat. Bot. 47, 249-264. Pietsch W. (1982) Makrophytische Indikatoren für die ökochemische Beschaffenheit der Gewässer. In: VEB, Gustav Fischer Verlag, Jena, 67-88. Pinnel N. (2007) A Method for Mapping Submerged Macrophytes in Lakes Using Hyperspectral Remote Sensing. PhD Thesis, Technische Universität München. Pip E. (1989) Water temperature and freshwater macrophyte distribution. Aquat. Bot. 34, 367-374. Pitelka L. & Ashmun J. (1985) Physiology and integration of ramets in clonal plants. In: Jackson J.B.C., Buss L.W. & Cook R.E. (eds.) Population Biology and Evolution of Clonal Organisms. Yale Univ. Press, New Haven, Connecticut, 399-436. Preston K.A. & Ackerly D.D. (2004) Allometry and evolution in modular organisms. In: Pigliucci M. & Preston K.A. (eds.) Modularity and Phenotypic Complexity. Oxford University Press, Oxford, 80-106. Puijalon S. & Bornette G. (2006) Phenotypic plasticity and mechanical stress: biomass partitioning and clonal growth of an aquatic plant species. Am. J. Bot. 93, 1090-1099. Rea T.E., Karapatakis D.J., Guy K.K., Pinder J.E.I. & Mackey H.E.jr. (1998) The relative effects of water depth, fetch and other physical factors on the development of macrophytes in a small southeastern US pond. Aquat. Bot. 61, 289-299. Reed D.C. (1990) An Experimental Evaluation of Density Dependence in a Subtidal Algal Population. Ecology 71, 2286-2296. Reusch T.B.H., Chapman A.R.O., Groeger J.P. (1994) Blue mussels, Mytilus edulis do not interfere with eelgrass, Zostera marina but fertilize shoot growth through biodeposition. Mar. Ecol. Prog. Ser. 108, 265-282. Reusch T.B.H., Stam W.T. & Olsen J.L. (1998) Size and estimated age of genets in eelgrass, Zostera marina, assessed with microsatellite markers. Mar. Biol. 133, 519-525. Royle R.N. & King R.J. (1991) Aquatic macrophytes in Lake Liddell New South Wales Australia biomass nitrogen and phosphorus status and changing distribution from 1981 to 1987. Aquat. Bot. 41, 281-298. Ryser P. & Eek L. (2000) Consequences of phenotypic plasticity vs. interspecific differences in leaf and root traits for aquisition of aboveground and belowground resources. Am. J. Bot. 87, 402-411. Sachs T. & Novoplansky A. (1997) What does aclonal organization suggest concerning clonal plants? In: De Kroon H. & Van Groenendael J (eds.) The Ecology and Evolution of Clonal Plants. Backhuys Publishers, Leiden, 55-78. Sánchez J.M., SanLeon D.G. & Izco J. (2001) Primary colonisation of mudflat estuaries by Spartina maritima (Curtis) Fernald in Northwest Spain: vegetation structure and sediment accretion. Aquat. Bot. 69, 15-25. Sand-Jensen K. (1998) Influence of submerged macrophytes on sediment composition and near-bed flow in lowland streams. Freshw. Biol. 39, 663-679. Sand-Jensen K. & Vindbaek Madsen T. (1991) Minimum light requirements of submerged freshwater macrophytes in laboratory growth experiments. J. Ecol. 79, 749-764. Santamaría L. & Rodriguez-Girones M.A. (2002) Hiding from swans: Optimal burial depth of sago pondweed tubers foraged by Bewick's swans. J. Ecol. 90, 303-315. Scheffer M. (1998) Ecology of Shallow Lakes. Chapman & Hall, London. Scheffer M., De Redelijkheid M.R. & Noppert F. (1992) Distribution and dynamics of submerged vegetation in a chain of shallow eutrophic lakes. Aquat. Bot. 42, 199-216. Scheffer M., Bakema A.H. & Wortelboer F.G. (1993) MEGAPLANT: a simulation model of the dynamics of submerged plants. Aquat. Bot. 45, 341-356. Scheffer M., Baveco J.M., DeAngelis D.L. & Van Rose K.A. (1995) Super-individuals a simple solution for modelling large populations on an individual basis. Ecol. Mod. 80, 161-170. Schlichting C.D. & Pigliucci M. (1998) Phenotypic Evolution: a Reaction Norm Perspective. Sinauer, Sunderland, USA.

105 References

Schmid B. (1986) Spatial dynamics and integration within clones of grassland perennials with different growth form. Proc. R. Soc. Lond. B 228, 173-186. Schmid B., Puttick G.M., Burgess K.H. & Bazzaz F.A. (1988) Clonal integration and effects of simulated herbivory in old-field perennials. Oecologia 75, 465-471. Schmieder K. (1995) Submersed macrophytes of two habitats in Lake Constance (Untersee) and their relations to chemical composition of surface and sediment interstitial water. Acta Bot.Gall. 142, 651-657. Schmieder K. (1997) Littoral zone-GIS of Lake Constance: A useful tool in lake monitoring and autoecological studies with submersed macrophytes. Aquat. Bot. 58, 333-346. Schmieder K. (1998) Submerse Makrophyten der Litoralzone des Bodensees 1993 im Vergleich mit 1978 und 1967. Ber. Int. Gewässerschutzkomm. Bodensee 46, 1-169. Scrosati R. & Dewreede R.E. (1998) The impact of frond crowding on frond bleaching in the clonal intertidal alga Mazzaella cornucopiae (Rhodophyta, Gigartinaceae) from British Columbia, Canada. J. Phycol. 34, 228-232. Sculthorpe C.D. (1967) The Biology of Aquatic Vascular Plants. Edward Arnold, London. Silvertown J. (1992) Cellular automaton models of interspecific competition for space - the effect of pattern on process. J. Ecol. 80, 527-534. Slade A.J. & Hutchings M.J. (1987a) Clonal integration and plasticity in foraging behaviour in Glechoma hederacea. J. Ecol. 75, 1023-1036. Slade A.J. & Hutchings M.J. (1987b) The effects of light intensity on foraging in the clonal herb Glechoma hederacea. J. Ecol. 75, 639-650. Slade A.J. & Hutchings M.J. (1987c) The effects of nutrient availability on foraging in the clonal herb Glechoma hederacea. J. Ecol. 75, 95-112. Smith C.S. & Adams M.S. (1986) Phosphorus transfer from sediments by Myriophyllum spicatum. Limnol. Oceanogr. 31, 1312-1321. Søndergaard M. (2007) Nutrient Dynamics in Lakes - with Emphasis on Phosphorous, Sediment and Lake Restoration. PhD Thesis, National Environmental Research Institute, University of Aarhus, Denmark. Spence D.H.N. (1967) Factors controlling the distribution of freshwater macrophytes with particular reference to the lochs of Scotland. J. Ecol. 55, 147-170. Spencer D.F. & Ksander G.G. (1992) Influence of temperature and moisture on vegetative propagule germination of Potamogeton species: Implications for aquatic plant management. Aquat. Bot. 43, 351-364. Spencer D.F. & Rejmanek M. (1989) Propagule type influences competition between two submersed aquatic macrophytes. Oecologia 81, 132-137. Spencer D.F., Van Vierssen W., Ryan F.J. & Ksander G.G. (1993) Influence of photoperiod and plant weight on tuber production by Potamogeton pectinatus L. J. Freshw. Ecol. 8, 1-11. Spencer D.F., Ksander G.G., Madsen J.D. & Owens C.S. (2000) Emergence of vegetative propagules of , Potamogeton pectinatus, Vallisneria americana, and Hydrilla verticillata based on accumulated degree-days. Aquat. Bot. 67, 237-249. Suzuki J. (2002) Shoot growth dynamics and the mode of competition of two rhizomatous Polygonum species in the alpine meadow of Mt. Fuji. Folia Geobot. Phytotaxon. 29, 203-216. Stapel J. & Hemminga M.A. (1997) Nutrient resorption from seagrass leaves. Mar. Biol. 128, 197-206. Stoll P., Egli P. & Schmid B. (1998) Plant foraging and rhizome growth patterns of Solidago altissima in response to mowing and fertilizer application. J. Ecol. 86, 341-354. Straile D. & Geller W. (1998) The response of Daphnia to changes in trophy and weather patterns: a case study from a large lake (L. Constance). ICES J. Mar. Sci. 55, 775-782. Strand J.A. & Weisner S.E.B. (1996) Wave exposure related growth of epiphyton: implications for the distribution of submerged macrophytes in eutrophic lakes. Hydrobiologia 325, 113-119.

106 References

Strand J.A. & Weisner S.E.B. (2001) Morphological plastic responses to water depth and wave exposure in an aquatic plant (Myriophyllum spicatum). J. Ecol. 89, 166-175. Stuefer J.F. (1998) Two types of division of labour in clonal plants: Benefits, costs and constraints. Perspect. Plant Ecol. Evol. Syst. 1, 47-60. Stuefer J.F., De Kroon H. & During H.J. (1996) Exploitation of environmental heterogeneity by spatial division of labour in a clonal plant. Funct. Ecol. 10, 328-334. Stuefer J.F., During H.J. & Schieving F. (1998) A model on optimal root-shoot allocation and water transport in clonal plants. Ecol. Mod. 111, 171-186. Sutherland W.J. & Stillman R.A. (1988) The foraging tactics of plants. Oikos 52, 239-244. Szcepanska W. & Szcepanski A. (1973) Emergent macrophytes and their role in wetland ecosystems. Polskie Archiwum Hydrobiologii 20, 41-50. Terrados J., Duarte C.M., Kenworthy W.J. (1997a) Experimental evidence for apical dominance in the seagrass Cymodocea nodosa. Mar. Ecol. Prog. Ser. 148, 263-268. Terrados J., Duarte C.M. & Kenworthy W.J. (1997b) Is the apical growth of Cymodocea nodosa dependent on clonal integration. Mar. Ecol. Prog. Ser. 158, 103-110. Thompson J.D., Gray A.J. & Mcneilly T. (1990) The effects of density on the population dynamics of Spartina anglica. Acta Oecol. 11, 669-682. Timms R.M. & Moss B. (1984) Prevention of growth of potentially dense phytoplankton populations by zooplankton grazing in the presence of zooplanktivorous fish in a shallow wetland ecosystem. Limnol. Oceanogr. 29, 472-486. Tomasko D.A. & Dawes C.J. (1989) Evidence for physiological integration between shaded and unshaded short shoots of Thalassia testudinum. Mar. Ecol. Prog. Ser. 54, 299-305. Urban N.R., Dinkel C. & Wehrli B. (1997) Solute transfer across the sediment surface of a eutrophic lake: I. Porewater profiles from dialysis samplers. Aquat. Sci. 59, 1-25. Valta-Hulkkonen K., Pellikka P., Tanskanen H., Ustinov A. & Sandman O. (2003) Digital false colour aerial photographs for discrimination of aquatic macrophyte species. Aquat. Bot. 75, 71-88. Van Dijk G.M. & Janse J.H. (1993) Modeling resource allocation in Potamogeton pectinatus L. J. Aquat. Plant Manage. 31 (Jan spec. ed), 128-134. Van Kleunen M., Fischer M. & Schmid B. (2000) Clonal integration in Ranunculus reptans: By-product or adaptation? J. Ev. Biol. 13, 237-248. Van Nes E.H. (2002) Controlling Complexity in Individual-Based Models of Aquatic Vegetation and Fish Communities. PhD Thesis, Wageningen University, Wageningen. Van Nes E.H. & Scheffer M. (2005) A strategy to improve the contribution of complex simulation models to ecological theory. Ecol. Mod. 185, 153-164. Van Nes E.H., Scheffer M., Van Den Berg M.S. & Coops H. (2002) Dominance of charophytes in eutrophic shallow lakes - when should we expect it to be an alternative stable state? Aquat. Bot. 72, 275-296. Van Nes E., Scheffer M., Van Den Berg M. & Coops H. (2003) Charisma: a spatial explicit simulation model of submerged macrophytes. Ecol. Mod. 159, 103-116. Vermaat J.E. & Verhagen F.C.A. (1996) Seasonal variation in the intertidal seagrass Zostera noltii Hornem.: coupling demographic and physiological patterns. Aquat. Bot. 52, 259-281. Vretare V., Weisner S.E.B., Strand J.A. & Granéli W. (2001) Phenotypic plasticity in Phragmites australis as a functional response to water depth. Aquat. Bot. 69, 127-145. Walser R. (1996) Zur Rolle der Makrophytenbestände im Bereich von Flußmündungen am Bodensee-Obersee. PhD Thesis, Universität Hohenheim, Stuttgart. Walser R. & Löffler H. (1993) Ein air-lift Netz zum Fang von Jungfischen in dichten Unterwasser- Pflanzenbeständen. Österr. Fisch. 46, 17-18.

107 References

Watson M.A. & Cook C.S. (1982) The development of spatial pattern in clones of an aquatic plant, Eichhornia crassipes Solms. Am. J. Bot. 69, 248-253. Weaver M.J., Magnuson J.J. & Clayton M.K. (1997) Distribution of littoral fishes in structurally complex macrophytes. Can. J. Fish. Aquat. Sci. 54, 2277-2289. Webster P.J., Rowden A.A. & Attrill M.J. (1998) Effect of shoot density on the infaunal macro-invertebrate community within a Zostera marina seagrass bed. Estuar. Coast. Shelf Sci. 47, 351-357. Wertz I. & Weisner S.E.B. (1997) Potamogeton pectinatus and Myriophyllum spicatum response to sediments from a calcareous, shallow, eutrophic lake. J. Freshw. Ecol. 12, 1-10. Wessels M. (1998) Geological history of the Lake Constance area. Arch. Hydrobiol. Spec. Iss. Adv. Limnol. 53, 1- 12. Wetzel R.G. (1983) Limnology, 2nd edn. Saunders, Philadelphia. Wijesinghe D.K. & Handel S.N. (1994) Advantages of clonal growth in heterogeneous habitats: An experiment with Potentilla simplex. J. Ecol. 82, 495-502. Wikberg S. & Svensson B.M. (2003) Ramet demography in a ring-forming clonal sedge. J. Ecol. 91, 847-854. Wolfer S.R. & Straile D. (2004a) Density control in Potamogeton perfoliatus L. and Potamogeton pectinatus L. Limnologica 34, 98-104. Wolfer S.R. & Straile D. (2004b) Spatio-temporal dynamics and plasticity of clonal architecture in Potamogeton perfoliatus L. Aquat. Bot. 78, 307-318. Worm B. & Reusch T.B.H. (2000) Do nutrient availability and plant density limit seagrass colonization in the Baltic Sea? Mar. Ecol. Prog. Ser. 200, 159-166. Worm B., Reusch T.B.H. & Lotze H.K. (2000) In situ nutrient enrichment: Methods for marine benthic ecology. Int. Rev. Hydrobiol. 85, 359-375. Wortmann J., Hearne J.W. & Adams J.B. (1997) A mathematical model of an estuarine seagrass. Ecol. Mod. 98, 137-149. Xiao K., Yu D., Xu X. & Xiong W. (2007a) Benefits of clonal integration between 16 interconnected ramets of Vallisneria spiralis in heterogeneous light environments. Aquat. Bot. 86, 76-82. Xiao K., Yu D. & Wu Z. (2007b) Differential effects of water depth and sediment type on clonal growth of the submersed macrophyte Vallisneria natans. Hydrobiologia 589, 265-272.

108 References

109 References

110 Summary

Summary

Clonal growth governs the aboveground growth pattern of submerged clonal macrophytes. This study provides detailed data on the architectural growth rules of Potamogeton perfoliatus L. (perfoliate pondweed) such as rhizome spacer lengths, biomass allocation, branching frequencies, branching angles, as well as maps of rhizome networks and their seasonal development at Lake Constance.

Besides architectural reasons, the causes for different sizes and shapes of P. perfoliatus patches are dependent on nutrients, light conditions, and clonal integration. Branching was positively related to both irradiation and nutrient supply. Rhizome spacer lengths were negatively related to nutrient supply and positively related to irradiation. Where absolute rhizome lengths did not respond to the habitat conditions, relative rhizome lengths did. Rhizome allocation tended to increase with irradiation and shoot allocation tended to decrease with irradiation. Root allocation was higher at low nutrient supply. Clonal growth is further complicated because the species shows complex sharing of of resources between older and younger ramets, aiming at optimal resource partitioning.

We present an individual-based clonal model which shows that architectural growth affects patch characteristics: both rhizome spacer lengths and growth of rhizome axes relative to those in the previous year considerably determine patch density and patch expansion rates. Furthermore, the model is able to simulate seasonal growth of Potamogeton patches.

All results on clonal architecture are discussed in the light of foraging theory and allometric rules.

111 Summary

112 Samenvatting

Samenvatting

De bovengrondse groeipatronen van ondergedoken klonaal groeiende waterplanten worden in grote mate bepaald door de processen die klonale groei bepalen. Deze studie geeft gedetailleerde informatie over de factoren die de klonale groei van Potamogeton perfoliatus L. (Doorgroeid fonteinkuid) bepalen. Onder andere wordt de lengte tussen de scheuten van de wortelstokken (rhizomen), de biomassa verdeling, vertakkingsfrequenties en -hoeken precies beschreven. Bovendien zijn gedurende een seizoen gedetailleerde kaarten van de ondergrondse en bovengrondse delen van P. perfoliatus in de Bodensee bepaald.

De groeivormen en -partonen van P. perfoliatus zijn ook afhankelijk van omgevingsfactoren zoals nutriënten (meststoffen), licht en voeding vanuit andere delen van de plant via de wortelstokken. De vertakkingsfrequentie was positief gerelateerd aan zowel de lichtbeschikbaarheid en de nutriënten. De lengte van de rhizomen tussen de scheuten was negatief gerelateerd aan de beschikbaarheid van nutriënten en positief met de hoeveelheid licht. Terwijl de totale lengte van de rhizomen niet afhankelijk was van het habitat, was de relatieve lengte tussen de scheuten dat wel. De allocatie van biomassa aan rhizomen nam meestal toe met de hoeveelheid lichtbeschikbaarheid, en de allocatie van biomassa aan de scheuten nam meestal af met licht. Klonale groei is complex omdat voedsel op een complexe manier uitgewisseld wordt tussen oude en jonge scheuten. Experimenteel is het belang van de rhizomen ook vastgesteld.

We presenteren een individu-gebaseerd ruimtelijk model dat de klonale groei van P. perfoliatus over meerdere seizoenen beschrijft. Het model toont aan hoe fotosynthese en groeiregels van scheuten ruimtelijke patronen in een macrofytenveld kan genereren. Zowel de lengte tussen de scheuten en de groeirichting van een rhizoom ten opzichte van het vorige jaar bepaalt in grote mate hoe dicht de scheutdichtheidzal zijn en hoe snel een macrofytenveld zich kan uitbreiden.

De resultaten van deze studie worden in het licht van verschillende theorieën besproken. Met name de ‘foraging’ theorie en allometrische regels.

113 Samenvatting

114 Acknowledgement

Acknowledgement

Everything began when my supervisor Dietmar Straile proposed me as candidate for the project. He was always willing to support me with ideas and to help me with solutions in the areas of experimental set-up, analysis, and article publishing. Special thanks to Marten Scheffer, my main supervisor, who was very encouraging throughout the project and gave me the right directions during manuscript writing, and to Egbert van Nes, my co-supervisor, who taught me about modelling, gave me a warm welcome in Wageningen and supported me in every way. Thanks also to Hugo Coops and Marcel van den Berg for their discussions and comments on manuscripts. The project could be realised due to financial support by SFB454, University of Constance and RIZA, Lelystad.

Data collection would have been impossible without the research divers of the Limnological Institute at the University of Constance: Philipp Fischer, Thomas Jankowski, Martin Mörtl, Kirsten Pohlmann, Claudia Tesch, Oliver Walenciak, and Martin Wolf. In particular the sampling campaigns lead by Martin Mainberger and Addi Müller provided valuable data. John van Schie, RWS, took underwater photographs and kindly contributed photos for the thesis text and cover.

I also want to thank the student assistants who worked hard and reliably for me: Andreas Becker, Christoph Berron, Oliver Dehus, Patrick Fink, Armin Glaschke, Doris Greilinger, Robert Greither, Reiner Hoegg, Michael Merz, Sven Seibert, and Christian Weiss. The mesocosm experiments could be realised thanks to my colleagues Andreas Becker, Josef Halder, Thomas Jankowski, Klaus Jöhnk and Dietmar Straile who at one point spontaneously agreed to shovel 8 m³ of lake sediment (provided by the community of Moos) from a truck into the experimental basin.

Many other people contributed to the success of the project: Martin Wolf’s technical skills proved to be indispensable for the construction of sampling and measuring apparatus. Sandra Röck and Josef Halder navigated the boat with great skill and patience, Alfred Sulger taught me to operate a boat by myself. I gratefully acknowledge the lab assistance of Claudia Feldbaum, Christine Gehbauer, Petra Merkel and Gisela Richter. Ursula Haake looked after my computer and managed to get it running again after each breakdown.

Although only partly related to the topic of my thesis, I very much enjoyed the cooperation with the DLR Oberpfaffenhofen under the leadership of Thomas Heege. Due to the fact that aerial images can only be taken when the wheather is sunny, I spent beautiful days under blue skies on the lake with his group. Thanks also to Nicole Pinnel with who I share memories of challenging and adventurous sampling campaigns. Much of my knowledge about identification of aquatic macrophytes I owe to Michael Dienst and Klaus Schmieder.

I would like to thank my friend Ramesh K. Maskey who gave me a crash course in diving (at the pool) after I was accepted for the research post. In fact, that was the first time in my life I ever put my head below the water surface. Thanks to all family members and friends who urged me to finish after a lengthy break, especially to my parents who accepted my interests, visited frequently and helped with everything. The development of my propositions profited due to frequent discussions with my ‘paranymphs’ Heinz-Dieter Reinhardt and my husband. I am thankful about the support of my husband, who always reminded me that problems were never as big as they looked at first glance.

Last but not least, I would like to include Lake Constance in my acknowledgement because it is not just a basin filled with water but an extremely beautiful and inspiring place on earth.

115 Acknowledgement

116 List of publications

List of publications

Wolfer S.R., van Nes E.H. & Straile D. (2006) Modelling the clonal growth of rhizomatous macrophytes. Ecological Modelling 192, 67-82.

Wolfer S.R. & Straile D. (2004) Density control in Potamogeton perfoliatus L. and Potamogeton pectinatus L. Limnologica 34, 98-104.

Wolfer S.R. & Straile D. (2004) Spatio-temporal dynamics and plasticity of clonal architecture in Potamogeton perfoliatus. Aquatic Botany 78, 307-318.

Prinz D. & Wolfer S. (2002) Global Deforestation: Scale & Impacts. Proceedings, International Seminar on ‘Irrigation Water Management in an Intercultural Context, with Special Reference to Indonesia and Germany’. IWK, Universität Karlsruhe (TH), October/November 2002, 56-65.

Prinz D. & Wolfer S. (1999) Traditional Techniques of Water Management to Cover Future Irrigation Water Demand. Zeitschrift für Bewässerungswirtschaft 34, 41-60.

Prinz D. & Wolfer S. (1998) Water Harvesting in Trockengebieten. Möglichkeiten zur Deckung des Wasserbedarfs. Geographische Rundschau 50, 422-427.

Prinz D., Wolfer S., Siegert K. & Dessalegne Agaze T. (1997) Water Harvesting for Crop Production. FAO Training Course. FAO, Rome. 300 pp.

Geuppert S.R. & Chutter F.M. (1993) The Rapid Biological Monitoring of Water Quality. In: Proceedings, WISA Conference, 24.05.-27.05.1993. Durban.

117 List of publications

118 Curriculum vitae

Curriculum vitae

Name: Susanne R. Wolfer (nee Geuppert)

University education: Geography and biology, Universität des Saarlandes, Saarbrücken, Germany (1986 - 1991) and Simon Fraser University, Burnaby, Canada Main subjects: biogeography, vegetation geography, limnology, nature conservation, spatial planning, environmental science Thesis: ‘Die Entfernung von Stickstoffverbindungen aus dem spezifischen Emulsionsspaltwasser eines metallverarbeitenden Industriebetriebs’ (i.e. ‘The removal of nitrogen compounds from specific industrial waste water’) Graduation in Geography (Diplom)

04/1992 - 04/1993 Project leader, CSIR (Council for Industrial and Scientific Research), Division of Water Quality Information Systems, Pretoria, South Africa: - ecological monitoring of running waters

04/1994 - 08/1999 Coordinator of the MSc course ‘Resources Engineering’ (sustainable use of water, soil and biological resources), Universität Karlsruhe (TH), Germany

05/1999 - 05/2003 Researcher, Limnological Institute, Universität Konstanz, Germany: - clonal growth and patch dynamics of macrophytes

06/2003 up to date Freelance ecologist and GIS expert: - macrophyte mapping - geographic data updating, analysis and programming

119 Curriculum vitae

120