“Fossil Insect Study Suggests That Los Angeles Climate Has Been Relatively Stable for at Least 50,000 Years”

So sayeth the press release. But what about the paper, and the 182 beetles sampled from La Brea tar pits?

Fossil preservation in the tar pits is exceptional, but the constant stream of gas through the tar deposits mixes the fossils – there is no stratigraphy in the tar. Holden et al have to radiocarbon date each and every beetle they analyse. Naturally, this limits the number of beetles they can analyse, both because of the financial cost and the destruction of the beetles. It also limits the number of species than can be analysed (to seven), as only species with abundant fossils can be vaporised  (we are not told what other species are found – I would have liked this information).

Contrast the situation with the tar pits with a more typical site for beetle analysis, say a section through a peat bog, where large volumes of sediment in stratigraphic order can be collected, with dozens or more beetle fossils from many species in each sample, and only a few dates needed to constrain the chronology. La Brea is a challenging site.

These are dates of the seven beetle species at La Brea.


Holden et al figure 2. Median calibrated age of each beetle with 2-sigma ranges. Cases where only the error bar is shown are “greater than” ages where the results were very close to background and only a lower limit could be specified. The low quality of the image is because Elsevier are hopeless.

The first thing to note is that no beetles were dated to the last glacial maximum. Holden et al ascribe this to either the lack of insect collections from the pits containing LGM mammal fossils, or a cooler LGM climate making the tar less sticky so fewer insect were caught. If the former explanation is correct, the press release claim that “Los Angeles Climate has been relatively stable for at least 50,000 years” cannot be substantiated as there is no evidence from a critical interval. If the latter explanation is correct, the press release is refuted. A third explanation not considered by Holden et al is that the climate changed such that the seven beetles they use were not present at La Brea. Whatever the reason for the lack of LGM beetles, the headline of the press release is wrong.

There is also a beetle gap in the early Holocene thermal maximum. Again, the lack of evidence precludes a conclusion that the climate was stable.

It is always going to be difficult to reconstruct climate from a just seven species of beetles, selected in part because they were common. So it is not greatly surprising the reconstructed climate, for the intervals where there are data, is similar to modern. Holden et al report “mean summer temperatures within ±5 °C of today’s conditions”. A range wide enough that only the LGM (for which there are no data) could reasonably be expected to exceed.

The method for reconstructing climate is unclear. The paper appears to assume that the assemblage composition is constant through time and hence that the climatic conditions must have remained similar to modern. I would have liked to see a figure showing the assemblage in each time window. Something like this.


Number of beetles by time interval. ?Modern beetles have an age of < 200 cal BP, and so might represent accidental modern contamination of the tar pit.

There are distinct shifts in the species composition through time that Holden et al do not explore. I don’t know if these shifts have possible climatic interpretations. Although Holden et al collate modern records for their species, they don’t present the results in an easy-to-interpret way (they present violin plots, each species in a separate file, and scatter plots).

While the species distribution modelling could certainly have been done better, and the method for climate reconstruction made explicit, it is always going to be difficult to work on the stratigraphically-mixed deposits at La Brea, and without a huge amount of money for dating, and a willingness to atomise any beetles for dates, it will be impossible to get the data and reconstruction that could be expected from a more typical site. But that is no excuse for a press release that is unsupported by the paper.





Posted in climate, Peer reviewed literature | Tagged , | Leave a comment

Are Tibetan chironomids mesmerised by solar variability?

It’s been a while since I examined a paper that purports to present palaeoecological evidence for a climate response to solar variability. But last night, flicking through the recently published papers in Quaternary Science Reviews, I came across Zhang et al who suggest in their abstract that “solar activity could be an important mechanism driving the centennial-scale variability“.

Zhang et al present a high-resolution (50-yr) Holocene chironomid stratigraphy from the small Lake Tiancai on the south-east margin of the Qinghai-Tibetan Plateau from which they use a transfer function to infer Holocene temperature evolution. Unlike many authors, they present a range of transfer function diagnostics. These suggest there may be some problems with poor analogues for the fossil samples, especially in the late glacial/early Holocene, in the modern calibration set.

So far, so good. But what of the solar variability?

Alas, it is a correlation-by-eye.


Zhang et al Fig. 5. Multi-proxy records compared from Tiancai Lake (a) the percentage abundance of a pollen genus Tsuga; (b) December insolation at 30 °N; (c) chironomid-based mean July temperature reconstruction from this study (solid black line) with 5-point-running average (red solid line) highlighting the overall trend and (d) the percentage abundance of a diatom species Aulacoseira alpigena. Highlighted areas in blue correspond to the recognized Hallstatt solar cycle minima centered at 8200, 5500, 2500 and 500 cal yr BP.

Although the total solar irradiance reconstruction from Steinhilber et al. is shown in figure 6, figure 5 shows the totality of evidence for a link between chironomid-inferred temperatures and solar activity at Tiancai Lake. No statistical tests are used to test if this relationship is better than expected by chance.

The 5-point running-average has obviously been miscalculated – it should never reach the extremes of the data – it appears to be joining every fifth sample instead (it is done correctly in figure 6). This makes comparison of the smoothed record with the solar minima more difficult. Only one of the solar minima convincingly aligns with a temperature minimum and other large temperature minima do not align.

I would argue that this paper provides absolutely no evidence for a solar-temperature link, and that the purported link is a distraction to an otherwise good paper.

So, fulfilling Betteridge’s law of headlines, no, the chironomids were not mesmerised by solar variability. But the authors, reviewers and editors might well have been.



Posted in Peer reviewed literature, solar variability, transfer function | Tagged , | Leave a comment

Downloading Polish weather data

I needed some temperature data from Poland for some work I am doing on the pheonology of understorey plants in Białowieża Forest.

The easy source of data is the GHCN which can be accessed through the rnoaa package. This code loads the necessary packages and finds stations in GHCN within 200 km of Białowieża Forest.

#Download Bialowieza and regional temperature data

##regional climate
#find nearby sites
dat <- read.table(header = TRUE, sep = ",", text = "
longitude, latitude, id
23.894614, 52.744313,Białowieża"

stations <- meteo_nearby_stations(dat, lat_colname = "latitude", lon_colname = "longitude", station_data = ghcnd_stations(), var = "all",  year_min = 1960, year_max = 2000, radius = 200)

mp <- map_data("world", xlim = c(16, 30), ylim = c(48, 58))

ggplot(stations$Białowieża, aes(x = longitude, y = latitude, label = name)) +
geom_map(data = mp, mapping = aes(map_id = region), map = mp, fill = "grey80", colour = "black", inherit.aes = FALSE) +
geom_point() +
geom_label_repel() +
geom_point(data = dat, aes(x = longitude, y = latitude), colour = "red", size = 3, inherit.aes = FALSE)

Climate stations near Bialowieza Forest (red dot)

Now I can download the data for the closest few station

#download data
regionalData <- stations$Białowieża %>%
filter(distance < 70) %>% # 100 km radius
group_by(name, id) %>%
do(ghcnd_search(.$id, var = "TAVG")$tavg) %>% #download
mutate(tavg = tavg/10,
variable = "tavg") %>%
rename(value = tavg)

#plot regional data
g <- regionalData %>% filter(year(date) == 2000) %>%
ggplot(aes(x = date, y = value, colour = name, group = name)) +
geom_line() +
labs(x = "Date", y = "Daily mean temperature °C", colour = "Station")


Mean daily temperature for the year 2000 for the four nearest GHCN stations

Unfortunately, only a small proportion of the Polish weather data are available through GHCN. Until last year, these data were difficult to access. Now they are available to download from https://dane.imgw.pl/.

It is only possible to download seven days’ data at once. This would become tedious if you wanted data for several years. I wanted 50 years’ data, so I wrote an R function to hit the server ~2500 times slowly (server limit of 1000 queries per client per ten minutes). To use the script, you need to know the site ID for the weather station (map), the variable name (hope your Polish is good), and a registration. What I haven’t found are meta-data showing which stations have which variables and for how long.

## Bialowieza data from

#authentic <- "richard.telford@uib.no:Pa55w0rd"#not real password
#save(authentic, file = "data/authentic.Rdata")
startDate <- as.Date("2000-1-1")
endDate <- as.Date("2000-12-31")
siteCode <- "252230120"

BialowiezaDaily <- get_Polish_weather_data(
siteCode = siteCode,
variableCode = "B100B007CD",
startDate = startDate,
endDate = endDate,
authentic = authentic

BialowiezaDaily <- BialowiezaDaily %>%
mutate(name = "Białowieża", variable = "tavg")

And a plot to compare the Białowieża data with the GHCN data

## compare Białowieża with regional data
all_temperatures <- bind_rows(BialowiezaDaily, regionalData)

g %+% filter(all_temperatures, year(date) == 2000)

The four nearest GHCN stations and Bialowieza

Not surprisingly, all the data are in good agreement, but note that the GHCN-daily data are not homogenised to account for station moves etc. I suspect the IMGW data are not homogenised either.

Posted in Uncategorized | Tagged , | Leave a comment

Chironomid vs pollen: Holocene climate change in southern Europe

Pollen-inferred summer temperature reconstructions from southern Europe show cool early-Holocene summers and warmer late-Holocene summers (Davis et al 2003, Mauri et al 2015). In contrast,  warm early-Holocene summers are reconstructed elsewhere in Europe and most of the mid-high latitude Northern Hemisphere, as expected from the increased summer insolation due to orbital forcing.

Since one important use of palaeoclimate reconstructions is to validate climate models by testing whether model simulations of past climate match the climate reconstructed from proxies, it is important to understand whether these cool early-Holocene summer reconstructions are valid. Uncritical use of palaeodata to validate models is not going to be useful.

Time to call in the cavalry chironomids.

Samartin et al 2017 (which we covered in our reading group last week – not all ideas here are my own) reports chironomid-inferred summer temperature reconstructions from two lakes in the Italian mountains. Both (yes there is replication) lakes show a warm early Holocene, consistent with the insolation changes, marine proxy records and climate models, but not the pollen-based reconstructions.


Samartin et al figure 3. Chironomid-inferred mean July air temperature for Lakes Gemini and Verdarolo compared with other palaeotemperature records.

Samartin et al use a combined Swiss-Norwegian chironomid calibration set and report that all the fossil samples have fair/good analogues in this calibration set. I would have liked to have seen a plot of the analogue quality in the supplementary material. Otherwise, I like this paper.

Samartin et al suggest that the pollen-based reconstructions are so different because temperature is not the only control on vegetation in southern Europe. Moisture availability is important, and several millenia of disturbance by humans have replaced the natural vegetation with “humanized vegetation types” such as Garrigue. Thus pollen assemblages from the early Holocene may lack appropriate modern analogues and so pick inappropriate analogues from higher elevations where there are remnant forests. (In June, I am going to a PAGES workshop on improving pollen-based reconstruction. I must remember that in some cases the best reconstruction is no reconstruction.)

There is one problem with Samartin et al. Nature has strict guidelines about the availability of data etc.

A condition of publication in a Nature journal is that authors are required to make materials, data, code, and associated protocols promptly available to readers without undue qualifications.

A data availability statement is to be included in papers which

will report the availability of the ‘minimal data set’ necessary to interpret, replicate and build on the findings reported in the paper.

This is the full data availability statement from Samartin et al:

The chironomid-inferred mean July air temperature records from Gemini and Verdarolo are available at the NOAA National Centers for Environmental Information webpage (https://www.ncdc.noaa.gov/paleo/study/21030). Data from reference 1 (shown in Fig. 3) have been digitized from the original publication. Data from ref. 3 (online resource http://ncdc.noaa.gov/paleo/study/18317), ref. 21 (IGBP PAGES/World Data Center for Paleoclimatology Data Contribution Series no. 2006-106) and ref. 46 (IGBP PAGES/World Data Center for Paleoclimatology Data Contribution Series no. 92-007) have been downloaded from the NOAA National Centers for Environmental Information webpage.

Climate model data from refs 36,37 (shown in Fig. 3) are available at the NCAR climate data webpage (https://www.earthsystemgrid.org/dataset/ucar.cgd.ccsm3.trace.html). Climate model data from ref. 6 can be obtained from H.R. on request.

Notice anything missing?

There is no mention of where the chironomid calibration set and fossil data are archived or how they can be accessed. I do hope this omission was just an oversight.

There is no better explanation of why the data should be made available than this sentence from the statement.

Data from reference 1 (shown in Fig. 3) have been digitized from the original publication.

Posted in Peer reviewed literature, transfer function | Tagged , , , , , | Leave a comment

Resilience workshop in Finse: deadline soon

At the end of March, there is a workshop on Measuring Components of Resilience in Long-term Ecological Datasets

The workshop is at Finse, high in the mountains between Bergen and Oslo. Along with debate about how ecosystem resilience can estimated from palaeoenvironmental data, and work to implement our ideas in R, I can promise there will be lots of snow, sauna and skiing. If you have size 39 feet, I can lend you some skis (other skis will be available). Northern lights are possible if the weather is good (we saw them last year). The food is excellent, and the train line to Finse from Bergen or Oslo is one of the most scenic anywhere.


Not just ptarmigan and arctic hares & foxes at Finse

There are still a few days left to apply. All travel within Norway (i.e. trains to Finse and back), food and accommodation at the Finse Research Station, is provided. Some funds are available to support travel of early-career researchers.

Posted in Uncategorized | Tagged | Leave a comment

2016 in numbers

Highest temperature: 120 °C

At least that is what the (possibly broken) thermometer in the wood-fired sauna near Tallinn read when I was taken there and beaten with a viht by PhD students Annika and Agata. After a quick plunge into a hole in the ice-covered Viitna Pikkjärv and a discussion of the lake’s palaeoecological record, I repaid the vihtlema, concentrating on the soles of their feet.


The viht and the steam

Longest hiatus between submission and publication: 8 years

The diatom stratigraphy of Gårdsfjärden is now published.

New R packages learnt: three and a half

I started the year having used ggplot2only once, under duress. I now use it for almost all my figures. I have never used dplyr or tidyr and thought the pipe syntax (%;gt;%) looked rather odd. I’m now a convert, finding that dplyr makes it easy to write understandable code. I also learnt to use plyr a bit. I find it very useful for some jobs, but wish I wouldn’t get such unfathomable bugs when I accidentally load dplyr and plyr at the same time.

Coldest temperature: -13 °C

It seemed like a good idea to go skiing before breakfast at the Palaeo-Resilience workshop at Finse. Beautiful but cold before the sun came up. I should have worn a third pair of gloves.

You can get your chance to think about ecosystem resilience while watching the northern lights as you cool down from a sauna this March at the workshop “Measuring Components of Resilience in Long-term Ecological Datasets“. Alas, no viht.


Dawn at Finse

Most giraffes:  ~15

In 2015, 15 Rothschild’s giraffes were translocated to Lake Mburo from other Ugandan National Parks in the hope that this red-listed species would re-establish and help to control the acacia  that is converting the savannah into woodland.

I went to Uganda for Perpetra Akite’s, who was my PhD student, PhD defence. Managed to get away from Kampala for a couple of days to Lake Mburo National Park where I saw most, perhaps all of the herd. And lots of zebra.


Steffi & co

Number of blog posts on solar variability: 0

I found a few papers reporting links between palaeoecological data and solar variability that I wanted to discuss, but never found the time. I’ll try to cover more this year as I expect to write fewer posts about Polish lakes.

Number of email from researchgate.com: 431

At least that is how many are in my email trash folder. It seems like many more than that.

Most bison: 0

There are a lot of trees in Białowieża Forest, and it only takes one to hide a bison. Fortunately, I’m going back at the end of January to work on a 52-year phenological data-set (understorey plants, not bison), and the bison should be more visible when they come out of the forest to graze.

Posted in Uncategorized | Tagged , , , | Leave a comment

Performance of chironomid-temperature transfer functions

Chrironomid (non-biting midges) are sensitive to temperature (Eggermont and Heiri 2011), and because their head-capsules preserve well in lake sediments, their fossil assemblages can be use to reconstruct past temperature with a transfer function trained on the chironomid-temperature relationship in a modern calibration set of chironomid assemblages and associated temperatures. They are a often used to provide a independent temperature reconstruction against which pollen or macrofossil-inferred vegetation changes can be compared (e.g., Samartin et al. 2012).

The performance of transfer functions for reconstructing past temperatures from chironomid assemblages in lake sediments varies enormously. The root-mean-squared-error of prediction (RMSEP) for transfer functions in a compilation by Eggermont and Heiri (2011), augmented by a few recent studies1, ranges between 0.55 and 2.37 °C. The lowest RMSEP suggests that chironomids give very precise reconstructions; the highest suggests that chironomids have limited utility for anything other than reconstruction glacial-interglacial changes. What causes this 4.3-fold difference in RMSEP? Understanding this might suggest where efforts to improve chironomid-temperature, and other, transfer functions should be focused.

Continue reading

Posted in transfer function | Tagged , | 2 Comments