A failed attempt to reproduce two ordinations

To examine the millennium-long chironomid-inferred air temperature reconstruction from Seebergsee (Larocque-Tobler et al 2012) is, after having shown that the calibration-in-time reconstruction for the upper section of the core (Larocque-Tobler et al 2011) has no skill, to flog the proverbial dead parrot.

But there is one aspect I wish to examine: Figure 5, which shows ordinations of the millennium-long chironomid stratigraphies from Seebergsee and Lake Silvaplana.

 

1-s2-0-s0277379112001230-gr5

Larocque-Tobler et al (2012) Fig. 5. Two dimensional non-metric multi-dimensional scaling in the fossil chironomid assemblages of Seebergsee and Lake Silvaplana. In both cores the samples after ca 1950 CE (black circles) have the highest chord distances with samples pre-1950 CE (white circles).

Both ordinations show a pronounced shift in community composition at 1950 CE, with all fossil samples from after this date distinct from those before. The ordinations are non-metric multidimensional scaling fitted and plotted using Primer (so the axes are scaled correctly).

One might expect such the strong pattern shown in the ordinations to be clearly visible in the underlying community data: they are not.

1-s2-0-s0277379112001230-gr3

Larocque-Tobler et al (2012) Fig. 3. Changes in the chironomid taxa (in percentages). Only the taxa from more than two samples are presented.

If there is a switch in community composition in Seebergsee, it would appear to be at about 1970 CE, where there is a zone boundary, rather than 1950 CE. It is unclear from the paper whether the zones were defined using the “Zone program” (Methods) or “based on the PCA scores of axis 1 and 2” (Results).

silvaplana_qsr.png

Larocque-Tobler et al (2010) Fig. 3. a) Number of head capsules per taxon found in each sample of the Lake Silvaplana cores. b) Percentages of taxa in merged samples of the Lake Silvaplana cores. Samples were merged to obtain a number of at least 30 head capsules (see Larocque et al., 2009).

The millennium-long chironomid stratigraphy from Lake Silvaplana (Larocque-Tobler et al 2010) does not have a zone boundary at 1950 CE; the assemblage change at ~1770 CE appears more important.

Both stratigraphies appear to have many more samples than shown in the ordination.

Figure 5 from Larocque-Tobler et al (2012) is the only published ordination of the Seebergsee chironomid stratigraphy, but there are several ordinations of the Lake Silvaplana stratigraphy. None show any indications of a large community change at 1950. For example, here is a correspondence analysis of the last 150 years from Larocque et al (2009).

silvaplana_2009_fig4

Figure 4 from Larocque et al (2009) Correspondence analysis. The numbers in the graph are the years AD of the samples. The variance explained by each axis is in brackets. The axes are incorrectly scaled.

Since I have the data for at least the last century from both lakes, I can try to reproduce the Figure 5.

sbs_nmds_fig-1

Non-metric multidimensional scaling of the chironomid stratigraphy from Seebergsee.

The NMDS of the last 150 years at Seebergsee (Larocque-Tobler et al 2011) shows no distinct split at 1950 or any other time.

sil_nmds_fig-1

Non-metric multidimensional scaling of the chironomid stratigraphy from Lake Silvaplana.

The NMDS of the > 400 year long chironomid stratigraphy from Lake Silvaplana (Larocque-Tobler et al 2009) shows a fairly distinct split into two groups. However, the split is at 1760 not 1950 (I had to use the Bray-Curtis distance as the chord distance used Larocque-Tobler et al (2012) found the last sample in the data to be an extreme outlier).

Neither ordination in Larocque-Tobler et al (2012) can be reproduced.

I had hoped this would be my last post exploring Seebergsee, but, alas, I found some further oddities while preparing this post. After discussing those, I have a another post about the chironomid-inferred reconstructions from Lake Zabinskie and, in the unlikely event that I feel inspired, the lakes at Abisko. Soon, I hope to resume my usual diet of papers reporting solar correlations with palaeoecological data.

Advertisements

About richard telford

Ecologist with interests in quantitative methods and palaeoenvironments
This entry was posted in Peer reviewed literature and tagged , . Bookmark the permalink.

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s