ecosystems and global change

New multi-lake experiment launched

It has been a busy summer for our NERC-funded project RELATED (that stands for Restoring Ecosystems by Linking Aquatic and Terrestrial Ecological Dynamics).  The project aims to test experimentally whether the productivity of aquatic food webs increases with the quantity and quality of terrestrial organic matter deposited in nearshore delta habitats.  It builds on our previous work that showed forests fuel the growth of juvenile fish by subsidizing the base of the aquatic food web.  RELATED will also feature much more work on unlocking the microbial ‘black box’ at the base of the food web as well as better understanding how greenhouse gas emissions might change with surrounding vegetation.  Inland waters are major sources of atmospheric carbon and predicting their responses to future change is of major interest (some great recent work here, here, and here).

Experiment behind the beautiful Living with Lakes Centre

New experiment behind the beautiful Living with Lakes Centre

We’ve now just finished launching the experimental platform behind RELATED.  Over the last 8 weeks, we’ve had a team of 8+ working tirelessly to submerge artificial lake sediments in 3 lakes.  This has involved collecting, mulching, and sifting organic and inorganic materials, mixing these at an industrial scale, and outfitting nearly 300 mesocosms with the appropriate sampling gear.  By replicating our experiment in 3 different lakes, we’ll be able to study terrestrial-aquatic linkages along gradients in eutrophication and climate change – the main drivers of change in the world’s inland waters.  We’ve also had excellent student help from Laurentian University’s School of Architecture work to design and build a network of sampling platforms that will allow us to work without disturbing our sediments.

Artificial sediments with increasing organics from left to right.

Artificial sediments with increasing organics from left to right.

Here’s the main team celebrating the deployment of the experiment in Ramsey Lake:


Explaining the origins of species diversity

The July 2015 issue of New Phytologist has come out and its a doozy!  The entire issue features work on evolutionary plant radiations, drawing together a range of papers that summarize the current state of knowledge about “where, when, why, and how” plant radiations happened.  An editorial by Colin Hughes, Reto Nyffeler Peter Linder outlines the content of what will surely be a landmark issue for years to come!

Many, if not all, of these papers were presented at a symposium organized in Zurich that we attended with our New Zealand collaborators Professors Bill and Daphne Lee in June 2014:

We're somewhere in the crowd

We’re somewhere in the crowd …

For me, this was possibly one of the most intellectually stimulating meetings that I’ve ever attended.  Talks drew major figures in plant ecology, evolution, and systematics and really pushed the boundaries out on ‘diversification’ research.  The field itself is still arguably quite new, with many of the key questions synthesized in a 2008 paper by Peter Linder.  In fact, we have a PhD studentship available to follow-up some of these questions and build on what we talked about at the meeting.

You might have noticed that we even have a contribution in the New Phytologist Special Issue.  Our paper tests the mechanisms by which plant evolutionary radiations emerge and influence ecological dynamics.  We apply more of our expertise in structural equation modelling to focus on 16 species-rich genera in the alpine zone of New Zealand.

Diversity in New Zealand’s alpine.  Celmisia, Chionochloa, Dracophyllum, and Veronica all appear in this photo.

One of the most exciting aspects of our paper is that we’ve tried to reconstruct the niche space that each genus has occupied over the last 20 million years.  This is fairly ambitious and has involved tasks like reconstructing sea surface temperatures through the Cenozoic from isotopic measurements of foraminifera deposited in marine sediment cores, and then using these to estimate past land temperatures.  We’ve also had to consider that the alpine zone has grown immensely over time with uplift of the Southern Alps.  To do so, we collated radiometric dates of rocks and tried to infer their rate of uplift since the Miocene.

Overall, our results suggest that genera that colonized New Zealand earlier encountered more ‘vacant’ environmental space, which promoted species diversification and further occupancy of the environment.  Genera that occupied more environmental space were subsequently more dominant in present-day vegetation plots.  The key message is that time not only explains why diversity arises, but how this diversity influences ecological dynamics.  The Special Issue has many other fabulous papers, so do check it out!

Mutualisms – a little help from friends?

It’s been a while since we’ve last blogged but we’re hoping to pick it up more regularly again over the next few months now that the group is quickly growing.

For now, I thought I’d comment on an interesting study published this month in Applied Vegetation Science by Laura Burkle and Travis Belote.  It links up nicely to work we’re doing on priority effects, which is the idea that species that arrive earlier into a habitat influence the interspecific interactions of later arriving colonists.  Priority effects are hardly surprising. Succession theory long predicted that communities develop along different trajectories depending on the species that are initially present.  But generalizing the direction and strength of priority effects remains challenging.  This is where Sgt. Pepper’s Lonely Hearts Club Band can offer some advice:

Continue reading

Talking allochthony in Cambridge and Lille

At the start of December, Andrew organized a three-day workshop in Cambridge, bringing over Brian Kielstra, John Gunn, Nikki Craig and Chris Solomon from Canada, Michael Pace, Grace Wilkinson and Stuart Jones from the USA, Jan Karlsson and Martin Breggen from Sweden, and Jon Grey from the UK. The aim was to see how terrestrially-derived organic matter (tOM) contributes to secondary production in aquatic ecosystems, by synthesizing global data collected from 594 observations on C, N and H, in over 10 zooplankton groups from 158 lakes in the northern hemisphere. Ultimately, the group would like to build a model for each consumer by lake. The ten limnologists worked all day (and all night!) but still had time to experience Cambridge, with time spent at the legendary Eagle pub (where great minds meet) and an exquisite dinner at Peterhouse, the University’s oldest college.


A few days later, while Jon, Martin, Jan, Nikki and Stuart traveled home, the rest of the group and I traveled to London to catch the Eurostar to Lille, where we attended the British Ecological Society (BES) and Société Française d’Ecologie (SFE) joint meeting. The meeting was the first of its kind and brought about 1,200 ecologists to the Lille Grand Palais convention center – mainly French and British, but also a lot of attendees from other European countries and the rest of the world.

Continue reading

Lakes turning to jelly

Our latest paper in the Proceedings of the Royal Society on the “jellification” of temperate lakes has gotten an impressive amount of on-line attention.  At the time of writing this blog, Altmetric scores it as the 22nd highest ranking paper ever published in the journal.  You can read summaries from the Washington Post, New York Times, Daily Mail, CBC, CBC Radio, and Yahoo, among others.  I’ve also done four separate interviews this week with BBC radio stations (BBC Radio 5, BBC Wales, BBC Cambridgeshire).  You can catch the latest, with the BBC World Service from the 26th of Nov, below:


The main finding of the paper is that a small planktonic animal named Holopedium glacialis has been dramatically increasing in two very different lake regions of Canada as the keystone grazer in these lakes, the water flea (Daphnia spp.), has been disappearing.  Our results show that this is mainly driven by declines in lake water [Ca].  Daphnia need large amounts of Ca to build their body shell, while Holopedium surround themselves in a gelatinous polysaccharide “bubble”:

An individual Holopedium with the jelly capsule clearly visible. (Photo: Ian Gardiner / E-Fauna BC)

This jelly also protects Holopedium from predators.  By contrast, Daphnia are increasingly susceptible to predators at low [Ca] because their ability to induce evolved defences is also impaired. Our analyses show how vanishing Daphnia have now left more algae uneaten for their competitors to exploit, allowing them to multiply in number.  Many media reports have picked up on this as Holopedium liking ‘pollution’, with low [Ca] somehow being the result of this.  But it is more in fact a legacy of pollution.  While we have curbed industrial emissions and reduced acid rain, the historical depletion of base cations from the thin soils of the boreal shield, have left behind much lower [Ca] than present prior to industrial activity.  Ca concentrations have consequently been falling across much of North America and Europe.

Continue reading

Fretting over R2 and outliers

Two things I’ve been thinking about lately.

First, outliers. A classical example of these arises in what is known as Anscombe’s quartet, four datasets with almost nearly identical properties (mean, variance, x-y correlation) yet very different when plotted:

Variation in Anscombe’s quartet despite nearly identical summary statistics.

John Kruschke at “Doing Bayesian Data Analysis” has a terrific example of how to fit robust regressions through these lines using BUGS, along with links to additional code from Rasmus Bååth.  It is definitely worth taking a look at these.  One of the key points is that for small sample sizes, where the population SD is really unknown, it is worth modelling data from a Student’s t distribution as opposed to a normal distribution, as it allows for more flexibility in outliers (such as in y3-y4 above).  These are easy to implement in BUGS and Stan.


Second, how worried should we be about R2 values?  I’ve been thinking briefly about this over the past week as I’m not getting my usual impressive 0.9 values.  Does this matter though?  What does a R2 really tell us?  Well, consider the example below:

x1 <- rnorm(10000)
y1 = 5 + 1*x1 + 1*rnorm(10000)
y2 = 5 + 1*x1 + 3*rnorm(10000)

summary(lm(y1 ~ x1))
summary(lm(y2 ~ x1))

If you run the code in R you will see that both models recover the estimated effects of x1 on y1 but that y2 has 3X larger error estimates around this effect.  The resulting R2 is ca. 1/5 that of the y1 model.  Does this mean the model is a poor fit?  Not really…  It speaks much more to the fact that the predictive power of the y2 model is diminished and there are wider prediction intervals, rather than there being anything wrong with the mean estimate of x1 per se.  I’d be curious to hear what others make.


Finally, I thought we should drop a mention of the student blog PLANeT, run by our Part II undergraduates.  It’s worth checking them out! They regularly post well written and engaging articles on the relationship between plants and society.

Summer of field work

Joanna Wolstenholme, a third year NatSci undergraduate, has just wrapped up seven weeks helping our field campaign in Canada.  She authored this entry, describing her experience.

Sudbury, on first inspection, is a rather spread-out mining town, inhabited by many trucks (most of them blue).  However the more you explore, the more remarkable the town becomes.  It is one of the few areas of the world where remediation has really worked, and the next generation will inherit a greener and cleaner city than the one that their parents inherited.  This remarkable change, from a barren ‘moonscape’ caused by years of acid rain (Sudbury was once the world’s largest point source of sulphur dioxide emissions, thanks to large-scale nickel and copper mining), to an area with burgeoning forest cover and recovering lakes, is a great success story that the area can be immensely proud of.

With this backstory, Sudbury, with its 330+ lakes, makes an ideal experimental location for a group dealing in ecosystems and global change. Our study lake, Daisy Lake, is perfectly set up for studying the effects of terrestrial influences on aquatic ecosystems.  Along its length, the shores and wetlands have recovered to various degrees. One catchment has even been limed – covered with calcium carbonate to neutralise the acidic soils, and so plant growth is relatively lush. Other areas, closer to the smelter at the north end of the lake, are far more barren; bare, stained rock predominates, with a few stunted trees.

One of the streams we study

One of the streams we study

In Daisy, we were studying eight stream deltas, each with very different personalities.  At each site Erik and I measured algae, sediment, and water. This all sounds very easy in theory, but in practice (as with any fieldwork, as I came to learn) things were far harder and more complicated… and often involved some rather novel solutions. If nothing else, this placement has certainly given me plenty of opportunities to stretch my problem solving skills!

My first job was to build algae-collectors, which were plastic tubes with cut up swim floats attached from which 6 microscope-slides dangled from fluorescent string. These floated on the surface, but we also sank clay pot holders as another surface for algae to grow on. We left these in the lake (on a beautiful sunny day) at each of the deltas and then returned to collect them 3 weeks later.  On a more high-tech note, we also made use of two chlorophyll fluorometers to characterise the algal species found in the water column and benthic layer. After several dry runs measuring the amount of algae on Erik’s office floor, we took them out to the lake, and used them at each of the deltas. The unseasonal amount of rain that Sudbury was experiencing, however, complicated things, and meant that in some sites Erik had to swim with the fluorometers, as we couldn’t reach the sediment from the boat.

This slideshow requires JavaScript.

As well as working on Daisy with Erik, I also helped Andrew collect additional data for his survey of terrestrial resource use by aquatic organisms. This meant going out to six other lakes around Sudbury, and six down in the Muskokas, to collect water samples, use fluorometers, and deploy and collect the microscope slide contraptions. Key to the project was collecting clean leaf and algal samples, to go off for stable isotope analysis, to allow Andrew to calculate the influence of the terrestrial systems on the lake ecosystems.

In order to grow clean algal samples without the influence of terrestrial DOM, we collected water from each of the lakes, then filtered it into jars and re-inoculated each jar with a small amount of unfiltered lake water, from which we hoped the algae would regrow. This seemed simple in theory, but involved hours of standing by a vacuum pump watching water drip through a filter. One night, we actually filtered water outside a hotel, so as not to set the fire alarms off! Safe to say we got many odd looks. However, the field trip down to the Muskokas was one of the best perks of the summer. We went down in September, almost at the peak of the colours changing, and had two lovely dry but crisp days. Driving down dirt tracks through beautiful forest, to find beautiful lakes to paddle out into was great fun, and a real adventure! It definitely offset the tedium of filtering.

This slideshow requires JavaScript.


The 'Hollywood sign' of Sudbury

The ‘Hollywood sign’ of Sudbury

At the end of my seven weeks here I am very sad to be leaving. It was a great experience, with plenty of messing about on boats, exploring new places, and making new friends. I have learnt a lot about the complications of fieldwork, how to solve problems on the fly with limited supplies, and just what really goes on behind those simple sounding ‘Materials and Methods’.


This article is also published on PLANeT.  Joanna’s trip was supported by funding from the Department of Plant Sciences and NERC.

Drooling over plants (literally!)

It is really nice to be able to write about our latest paper, which just came out in Biology Letters to a fair bit of fanfare.  It got picked up by Science, New Scientist – including in their 23 July print issue – and even the Royal Society had a piece about it.  BBC Radio recently had me out to talk to them:

Why so much attention?  Well the paper is one of those quirky but interesting scientific discoveries.  And it is based on drool.  Yup, you read that right, the saliva of large mammals.

Continue reading

Many measurements and few observations

Our University colleague Professor Sir David Spiegelhalter has written a brief opinion piece in the latest issue of Science on the future of probabilistic models, particularly for big datasets (think images or genomes).

Two points jumped out at me:

(1) Statistical problems have shifted from many observations (large n) and few parameters (small p) to small n and large p, creating pitfalls when testing large numbers of hypotheses.  This is because the standard “p-value”, which we’ve griped about in the past here, will declare 1 in 20 non-existent relationships “significant” simply by chance.  So procedures are needed to reduce false discoveries.  The bit that I didn’t really follow was why even bother minimizing false discoveries?  Wouldn’t an interpretation of effect sizes be more meaningful?

(2) Inferring causation from observational data will continue to be a challenge, especially when n gets cheap and p remains large.  Statistical theory to deal with causality will be needed more than ever, and thankfully, it is improving.  This is something we’re quite fond of having thought a fair bit about causality in the context of path analysis, structural equation modelling, and directed acyclic graphs (see our J Appl Ecol paper that just came out).  The problem, however, is that these approaches don’t come easy and I struggle to see how they can be used by non-statisticians (the models in our paper took years of faffing!).  Finding ways to make causal inference more accessible is going to be critical in the future.