Bless the rains down in Africa #DACEA3

An ultra-quick post today. I have been spending a lot of time lately writing a grant proposal (and occasionally tweeting about it  on the #DACEA3 hashtag).  Finally it’s in and after a celebratory beer or two at the famous Mikeller last week I have managed to get around to a very brief summary of what it’s all about… 

Around 17,000 years ago, Lake Victoria more or less completely dried out. I still find this absolutely staggering. In fact, the lake has dried out and reformed at least 3 times since it first formed about 400,000 years ago.

Lake Victoria is the largest lake in Africa and indeed the tropics, containing 2.75 cubic kilometres of water (though compared to the 2,850,000 cubic kilometres of water in the Greenland ice sheet that seems small, which merely goes to prove how much of our fresh water is locked up in the ice sheets), making it the 9th largest lake by volume in the world.

Gratuitous wildlife shot: A raft of hippos chilling out in the river. Photo: Pim Bussink
Gratuitous wildlife shot: A raft of hippos chilling out in the river.
Photo credit: Ruth Mottram

Clearly, the disappearance and later reappearance of the lake, and others in the region speaks to monumental shifts in the climate. The East African Rift Valley lakes are largely fed by the East African rains, long and short, delivered by the shifting position of the Intertropical Convergence Zone as the Earth’s seasons change bringing those life-giving rains.

This grant proposal started as idle speculation around the coffee machine (in the grand old scientific tradition) about how this was climatically possible and could it happen again? My colleague (and talented PI on the proposal) Peter Thejll had been reading a book about John Hanning Speke and Richard Burton (not that one) and their famous search for the source of the Nile and has some personal African connections, which prompted the conversation and it seemed obvious to try and find out what happens to the local circulation to allow the lake to dry out. A quick google search revealed an old friend, Dr. Sarah Davies at Aberystwyth University was researching this topic actively and it all fell into place.

Now, I can guess what you’re thinking – this is usually a glaciology or Arctic Climate blog, where on earth has all this Africa stuff come from? Well what happens in the Arctic does not necessarily stay in the Arctic.

There are a number of hypotheses as to the drivers of these changes in African rainfall, among which is the interesting observation that the periods of greater aridity correlate remarkably well with Heinrich events in the North Atlantic.

Heinrich events were first identified as layers of sediment most likely transported into the North Atlantic Ocean by icebergs, known as ice rafted debris – IRD. The southerly position of many of these layers thousands of kilometres from any ice sheets either at the present day or in the past suggests a truly extraordinary amount of icebergs and cold fresh water were discharged over a relatively short period of time, from a large ice sheet. The source of these sediments is most likely the gigantic Laurentide ice sheet of North America, but there is also some evidence of smaller contributions from the British and Fenno- Scandian ice sheets (which may or may not have been joined together across the North Sea depending on how you interpret the evidence). The physics behind this is that as the enormous amount of cold fresh water was discharged into the North Atlantic, the temperature and salinity changes were sufficient to push, or keep the ITCZ far to the south, preventing the rains one East Africa.

On the other hand, other research has linked the failure of the rains to El Niño and related phenomena such as the Indian Ocean dipole and the Walker circulation. Still other scientists have noted that these drying periods seem to correlate with orbital changes in the earth which would affect the seasonality, that is the annual cycle of seasons. It is known as orbital forcing as the Earth’s seasons are driven by changes in our orbit around the sun (have a look at the excellent Orbit documentary from the BBC for a very easy to follow and beautifully filmed introduction to the importance of our orbit around the sun if you’re not familiar with Milankovitch cycles etc).

Milankovitch cycles shown from ocean cores and an Antarctic ice core at the bottom compared with the theoretical cycles. Image: By Incredio (Own work) [CC BY 3.0 (http://creativecommons.org/licenses/by/3.0)], via Wikimedia Commons
Milankovitch cycles shown from ocean cores and an Antarctic ice core at the bottom compared with the theoretical cycles.
Image: By Incredio (Own work) [CC BY 3.0 (http://creativecommons.org/licenses/by/3.0)], via Wikimedia Commons
All of these hypotheses can be supported by correlations with palaeo evidence,  but to really disentangle the connections between different regions of the world and how they affect each other’s weather and climate, we need to use a climate model. Luckily, at DMI we have the perfect tool to hand, a global climate model including ice sheets, EC-Earth. Furthermore we also have a high resolution regional climate model, HIRHAM5, my usual tool of choice. Our friends Morten Dahl Larsen and Martin Drews at the Danish Technical University are experts in using hydrology models so the answer is obvious.

We want to use these model tools and an extensive archive of observations, helpfully curated by our project partners Sarah Davies and Henry Lamb at Aberystwyth University to test all these different ideas. As an extra spinoff from the project, the Aberystwyth group have been intensively involved with the collection and analysis of a new lake sediment core from Chew Bahir in Ethiopia, so it’s going to be pretty exciting seeing if we can get the models to replicate  these kind of records.

There is of course an extra urgency to this project. It’s not just a somewhat obscure academic question. A recent paper showed that the long rains have significantly reduced over the last decade, and about 300 million* people live in this region and rely on these rains for drinking water, hydroelectric power and agricultural production. During this period we have also seen rapid changes in the Arctic. Of course the two trends may not be connected, or may be linked via a common third factor which is why the physics of climate are so important to understand.

UPDATE 2: I had no time originally in the writing of this to add a little about our other project associate. One of the best things about doing science are the very smart and friendly people  you meet along the way. Social media has really helped here to keep in touch as it is a nomadic lifestyle. By sheer chance I noticed a familiar name in a tweet that seemed to have some direct relevance to the proposal as we were writing it.

Hycristal

John Marsham was an old friend from my student days at Edinburgh University who I had slightly lost touch with. Thanks to the efforts of facebook we were soon back in touch and he is one of the Investigators on the HyCRISTAL project, part of the hugely important Future Climates for Africa Project, funded by the Department for International Development (DFiD) and the Natural Environment Research Council (NERC) in the UK. DACIA has some really obvious parallels with this project, though where we would like to concentrate on past climates, they will be focusing on present day and future climates. We hope therefore to send our PhD student to collaborate with the HyCRISTAL and FCFA projects where our insights from palaeomodelling palaeodata can make a real difference to the way future climate change is adapted to in East Africa. It will be very nice to collaborate with John’s group at Leeds as well as the Aberystwyth group, now we just have to hope we get the money to do it..

Or, to put it another way, “bless the rains down in Africa” ** (As an aside, for years I had always heard this as “I miss the rains down in Africa”, assuming it was about someone from Africa who missed being there).

UPDATE 1: Having viewed the original pop video again, I am rather troubled by the casual racism, sexism and naked orientalism on display (yes it was the 70s but still…) so I think I prefer to post instead this particularly witty deconstruction courtesy of @spaceforpootling

*(based on a back of the envelope calculation based on population statistics from Wikipedia if you know the correct number do let me know).

**(Apologies if you now have cheesy 1970’s pop music going round your head all day… 🙂 )

Conversion Factors

The official end of the hydrological year in Greenland (1st September to 31st August) means I am rather busy writing reports to give an overview of where the ice sheet is this year and what happened. I will try to write a quick blogpost about this in the next week or so (in case you’re curious here’s a quick plot to show the entire annual SMB, see also: http://polarportal.dk/en/groenlands-indlandsis/nbsp/isens-overflade/)

Daily and accumulated surface mass budget of the Greenland ice sheet, 31st August, 2015, last day of the hydrological year
Daily and accumulated surface mass budget of the Greenland ice sheet, 31st August, 2015
Anyway, as I find I am constantly switching between Gigatonnes (or indeed Gigatons), cubic kilometres and sea level equivalent, here is a quick and handy guide to converting different units of mass, for my own use as much as anyone else.

1 gigatonne is 1 billion metric tonnes  (or 1 milliard if you like the old British style, that is one thousand million).

However, on the Polar Portal we usually reckon everything in water equivalent. This is to save having to distinguish between snow (with a density between ~100 kg/m3 when freshly fallen and ~350 kg/m3 m when settled after a few days), firn (snow that has survived a full annual cycle with a density up to ~800 kg/m3) and glacier ice (anything from ~850 kg/m3 to 900+). Water has a density (at 4C) of 1000 kg/m3

1 gigatonne of ice will still weigh 1 gigatonne when it is melted but the volume will be lower since ice expands when it freezes.

1 metric tonne of water is 1 cubic metre and 1 billion metric tonnes is 1 km3 (a cubic kilometre of water)

A cubic kilometre of ice does not however weight 1 gigatonne but about 10% less because of the density difference.

100 gigatonnes of water is roughly 0.28mm of sea level rise (on average, note there are big regional differences in how sea level smooths itself out).

Finally, 1 mm sea level rise is 360 Gt of ice (roughly the number of days in a year) 

EDIT: – thanks to ice sheet modeler Frank Pattyn and ice core specialist Tas van Ommen on Twitter for pointing out I’d missed this last handy conversion. Interestingly and probably entirely coincidentally this is very close to the amount of mass lost by the Greenland ice sheet reported by Helm et al., 2014 for the the period January 2011 – January 2014 (pdf here) of 375 +/-24 km3 per year.

Over the last 10 years or so, Greenland has lost on average around 250 Gigatonnes of ice a year (Shepherd et al., 2012), contributing a bit less than a millimetre to global sea level every year with some big interannual variability. This year looks like it will be a comparable number but we will have to wait for the GRACE satellite results in a couple of months to fill in the dynamic component of the mass budget and come up with our final number.

Of course, gigatonnes and cubic kilometres are rather hard to visualise so we have skeptical science to thank for this post that tries. And as aside, Chris Mooney wrote a nice piece in the Washington Post on the difficulties of visualising how much ice is being lost which contains the immortal  line “Antarctica is clearly losing billions of African elephants worth of ice each year”.

A question of observation?

It’s been a while since I lasted posted anything, not for want of ideas but mainly lack of time. I shall try to catch up over the next few weeks. For now I was inspired to write an ultra-quick post about a very trivial question that came up at work today. I think it really captures how observational meteorology works (or should work).

Today, a colleague, John Cappelen, (also known as Mr. Greenland observational data), happened to mention in passing that on the 15th July this year, the weather station at Summit on the Greenland ice sheet had transmitted back to us in Copenhagen, a temperature observation of 2.5°C. This was during one of the highest melt periods this summer.

Automatic weather station operating at Summit, June 2015
The automatic weather station doing it’s thing at Summit, June 2015. Photo: DMI

Bearing in mind that Summit Camp is at roughly 3,216m, this is a pretty high measured temperature. In fact it would be rather noteworthy, especially as it occurred on one of the highest melt days of the summer. Temperatures above 0°C at Summit are not unknown and the record, during the famous summer of 2012 when around 95% of the ice sheet surface experienced melt, the water sweeping away a bridge on the Watson River near Kangerlussuaq, was 3.6°C.

Now, my colleague is a very experienced and careful scientist. He had checked the observations and the temperatures before and after this measurement were well below zero, so, my colleague asked, was there any reason to believe this measurement or can we assume an instrument failure of some kind?

My office mate in the Arctic and Climate Research section and I obligingly had a quick look at our Polar Portal Greenland ice sheet surface plots (see below) and at the melt extent plots that are updated daily on the DMI website. We had to conclude there was no evidence of melt that high on the ice sheet and there was also no reason to believe that a sudden sharp warming had occurred at Summit on this day based on DMI’s own weather forecast. We then turned to check the weather plots, also on the polar portal and based on data from the European Centre for Medium Range Weather Forecasting (the ECMWF – probably the best weather forecast modellers in the world).

Again, the anomaly plots showed rather cold conditions prevailing over the ice sheet during this period, though at the same time very high melt and low surface mass balance from the ice sheet due to the clear skies.

Graphs showing area of the Greenland ice sheet experiencing melt conditions, compared with the average (dark grey line) and range of past summers (1990-2012), for more detail see the DMI website
Graphs showing area of the Greenland ice sheet experiencing melt conditions, compared with the average (dark grey line) and range of past summers (1990-2012), for more detail see the DMI website
Temperature record from Summit Camp for the last month.
Temperature record from Summit Camp for the last month.

Fortunately, due to the American Summit Camp we have access to a back-up dataset very close to this location and after a quick web search John Cappelen was able to confirm that indeed this measurement was an error as the nearby station has not seen anything like that during the period in question (see right).

This kind of thing happens all the time and is therefore not at all newsworthy or interesting enough to write a publication about. However, when a recent record high temperature in the UK can lead to 2 critical articles in the Daily Telegraph and a particularly vigorous exchange on twitter for Met Office scientist Mark McCarthy, as well as this corrective piece on the Carbon Brief blog, perhaps we should be more vocal about just how careful and critical we as scientists are about observations, including the ones we decide to discard as well as the ones we keep.

Surface mass balance of the Greenland ice sheet on the 15th July 2015. Intense melting around the margins led to very negative SMB (the red colours) during this period.
Surface mass balance of the Greenland ice sheet on the 15th July 2015. Intense melting around the margins led to very negative SMB (the red colours) during this period.

Addendum: I was alerted by this tweet from Gareth Jones, also a Met Office scientist, to some slightly strange cherry picking in the blogosphere of climate records from a couple of DMI stations in Greenland. These have apparently been used to claim no climatic warming trend in Greenland over the 20th Century (I’m not going to link to it).

Screenshot of tweet

Anyone who is really interested in the observational data could try checking these reports by Mr Greenland observations himself instead, here is a quick summary: 

Mean annual temperature in Copenhagen, Torshavn (Faeroes) and selected DMI weather stations in Greenland from 1873 - 2014. Figure from DMI
Mean annual temperature in Copenhagen, Torshavn (Faeroes) and selected DMI weather stations in Greenland from 1873 – 2014. Figure from DMI

Calling all students…

I’m off to the UK next week for a workshop at Sheffield University where we will discuss the Surface Mass Balance of the Greenland Ice Sheet. The ISMASS workshop includes all the main modelling groups and observation groups who are involved in assessing surface mass balance in Greenland. I will be representing DMI’s Greenland SMB work there (not an easy task condensing it down to a 20 minute talk!).

In the course of preparing my presentation I have been making plots and figures and really investigating some exciting results. Sadly, I very rarely get the chance to spend time on this these days and I am keen to recruit students to assist in this work. Should any potentially interested students want to discuss this at Sheffield do let me know.

At the risk of spoilers in my presentation, here for example is one showing how different ways of characterising the surface snow pack affects our estimates for surface mass balance, and how the effects of the specific changes can be very different in different years.

Surface mass balance map plots of Greenland
Surface mass balance for the hydrological year (Sep -Aug) ending in 2012 and 2013 calculated using HIRHAM5 with 2 different surface schemes. The maps on the right show the difference between the 2.

As I mentioned I rarely get enough time to analyse the output from our runs and I would be very happy to hear from any students who are interested in doing a project on our simulations. We have lots of MSc and Bachelors projects already listed on our website at DMI but we are always happy to hear new ideas from students on related topics. I have terabytes of data from simulations I would like to be properly analysed and this could be very interesting given we are talking about Greenland and the Arctic in the present day and in the future. It’s a really nice opportunity to work with some cutting edge research. I am also happy to hear from students who would like to do a summer project and for the right candidate I would be able to look into a paid “studentmedarbejderhjælper” position for a few months, especially if you are already a trained computer science candidate….

If you are an undergraduate looking into an MSc, I urge you to consider Denmark. EU citizens usually qualify for generous support grants (rare these days!) as we have a shortage of candidates wanting to study in the sciences in Copenhagen. The research and teaching are world class and done in English at MSc level. The possibilities for projects in Greenland are literally endless.

If you want any more details or to talk about any of the possibilities, do get in touch!

Changes in SW Greenland ice sheet melt

A paper my colleague Peter Langen wrote together with myself and various other collaborators and colleagues has just come out in the Journal of Climate. I notice that the Climate Lab Book regularly present summaries of their papers so here I try to give a quick overview of ours. The model output used in this run is available now for download.

The climate of Greenland has been changing over the last 20 or so years, especially in the south. In this paper we showed that the amount of melt and liquid water run off from the ice sheet in the south west has increased at the same time as the equilibrium line (roughly analogous to the snow line at the end of summer on the ice sheet) has started to move up the ice sheet. Unlike previous periods when we infer the same thing happened this can be attributed to warmer summers rather than drier winters.

Map showing area around Nuuk
The area we focus on in this study is in SW Greenland close to Nuuk, the capital. White shows glaciers, blue is sea, brown is land not covered by ice.

We focused on the area close to Nuuk, the capital of Greenland, as we had access to a rather useful but unusual (in Greenland) dataset gathered by Asiaq the Greenland survey. They have been measuring the run off from a lake near the margin of the ice sheet for some years and made this available to us in order to test the model predictions. This kind of measurement is particularly useful as it integrates melt and run-off from a wider area than the usual point measurements. As our model is run at 5.5 km resolution, one grid cell has to approximate all the properties of a 5.5 km grid cell. Imagine your house and how much land varies in type, shape and use in a 5.5 km square centred on your house and you begin to appreciate the problems of using a single point observation to assess what is essentially an area simulation! This is even more difficult in mountainous areas close to the sea, like the fjords of Norway or err, around south west Greenland (see below).

Represent this in a 5.5km grid cell, include glacier, sea and mountain...  Godthåbsfjord near Nuuk in August
The beautiful fjords near Nuuk. Represent this in a 5.5km grid cell…

The HIRHAM5 model is one of very few regional climate models that are run at sufficiently high resolution to start to clearly see the climate influences of mountains, fjords etc in Greenland, which meant we didn’t need to do additional statistical downscaling to see results that matched quite closely the measured discharge from the lake.

Graph comparing modelled versus measured discharge as a daily mean from Lake Tasersuaq near Nuuk, Greenland. The model output was summed over the Tasersuaq drainage basin and smoothed by averaging over the previous 7 days. This is because the model does not have a meltwater routing scheme so we estimated how long it takes for melt and run-off fromt he ice sheet to reach this point.
Graph comparing modelled versus measured discharge as a daily mean from Lake Tasersuaq near Nuuk, Greenland. The model output was summed over the Tasersuaq drainage basin and smoothed by averaging over the previous 7 days. This is because the model does not have a meltwater routing scheme so we estimated how long it takes for melt and run-off from the ice sheet to reach this point.

We were pretty happy to see that HIRHAM5 manages to reproduce this record well. There’s tons of other interesting stuff in the paper including a nice comparison of the first decade of the simulation with the last decade of the simulation, showing that the two look quite different with much more melt, and a lower surface mass balance (the amount of snowfall minus the amount of melt and run – off) per year in recent years.

Red shows where more snow and ice melts than falls and blue shows where more snow falls than is melted on average each year.
Red shows where more snow and ice melts than falls and blue shows where more snow falls than is melted on average each year.

Now, as we work at DMI, we have access to lots of climate records for Greenland. (Actually everyone does, the data is open access and can be downloaded). This means we can compare the measurements in the nearest location, Nuuk, for a bit more than a century. Statistically we can see the last few years have been particularly warm, maybe even warmer than the well known warm spell in the 1920s – 1940s  in Greenland.

Graphs comparing and extending the model simulation back in time with Nuuk observations
Graphs comparing and extending the model simulation back in time with Nuuk observations

There is lots more to be said about this paper, we confirm for example the role of increasing incoming solar radiation (largely a consequence of large scale atmospheric flow leading to clearer skies) and we show some nice results which show how the model is able to reproduce observations at the surface, so I urge you to read it (pdf here) but hopefully this summary has given a decent overview of our model simulations and what we can use them for.

I may get to the future projections next time…

The Present Day and Future Climate of Greenland

Regional Climate Model Data from HIRHAM5 for Greenland

In this post I am linking to a dataset I have made available for the climate of Greenland. In my day job I run a Regional Climate Model (RCM) over Greenland called HIRHAM5 . I will write a simple post soon to explain what that means in less technical terms but for now I just wanted to post a link to a dataset I have prepared based on output from an earlier simulation.

Mean annual 2m  temperature over Greenland (1989 - 2012) from HIRHAM5 forced by ERA-Interim on the boundaries
Mean annual 2m temperature over Greenland at 5km resolution (1989 – 2012) from HIRHAM5 forced by ERA-Interim on the boundaries [Yes I know it’s a rainbow scale. Sorry! it’s an old image – will update soon honest…]

This tar file gives the annual means for selected variables at 0.05degrees (5.5km) resolution over the Greenland/Iceland domain.

I am currently running a newly updated version of the model but the old run gave us pretty reasonable and could be used for lots of different purposes. I am very happy for other scientists to use it as they see fit, though do please acknowledge us, and we especially like co-authorships (we also have to justify our existence to funding agencies and governments!).

This is just a sample dataset we have lots of other variables and they are available at 3 hourly, daily, monthly, annual, decadal timescales so send me an email (rum [at] dmi [dot] dk) if you would like more/a subset/different/help with analysis of data. This one is for the period 1989 – 2012. I have now updated it to cover up to the end of 2014. The new run starts in 1979 and will continue to the present and has a significantly updated surface scheme plus different SST/sea ice forcing and a better ice mask.

I have also done some simulations of future climate change in Greenland at the same high resolution of 5km using the EC-Earth GCM at the boundaries for RCP4.5 and RCP8.5 scenarios which could be fun to play with if you are interested in climate change impacts in Greenland, Iceland and Arctic Canada.

Mean annual 2m temperature change between control period (1990 - 2010) and end of the century (2081 - 2100) under RCP45 from HIRHAM5 climate model runs forced by EC-Earth GCM at the boundaries
Mean annual 2m temperature change between control period (1990 – 2010) and end of the century (2081 – 2100) under RCP45 from HIRHAM5 climate model runs forced by EC-Earth GCM at the boundaries.  This plot shows the full domain I have data for in the simulations.

This run should be referenced with this paper:

Quantifying energy and mass fluxes controlling Godthåbsfjord freshwater input in a 5 km simulation (1991-2012), Langen, P. L., Mottram, R. H., Christensen, J. H., Boberg, F., Rodehacke, C. B., Stendel, M., van As, D., Ahlstrøm, A. P., Mortensen, J., Rysgaard, S., Petersen, D., Svendsen, K. H., Aðalgeirsdóttir, G.,Cappelen, J., Journal of Climate (2015)

http://journals.ametsoc.org/doi/abs/10.1175/JCLI-D-14-00271.1 

PDF here

Finally I should acknowledge that this work has been funded by a lot of different projects:

Picture4

The Arctic Tern

Sterna Paradisaea is the Arctic tern. This amazing little bird lives around 30 years in the wild and every year completes the longest migration in the world, flying from the Arctic, where it breeds, to the Antarctic to feed. It sees more daylight than any other creature and is not only a great endurance flier but a marvellously agile one, a true aerial acrobat like a marine version of the more familiar swallow. It is also a brave parent, chasing away much bigger animals from it’s ground nesting colonies and braving attacks by skuas and other aerial pirates to feed its young.

It has been my privilege and delight to observe these birds while doing fieldwork in Iceland, Svalbard and Greenland. Some years ago while working out of UNIS in Svalbard, I found a desk in an empty office that overlooked a colony. Over 4 short weeks, albeit of continuous daylight, I was able to watch the spectacular aerial courtship displays, the ultra brief matings, the brooding of eggs, the feeding and raising of young and the eventual fledging of new individuals. I had never known I was a nascent birdwatcher before this.

I have decided to start a blog, partly to share my opinions with the world (and thus avoid boring my long-suffering friends), but mainly to improve my writing by preparing short easy to read pieces in an every day style.  I chose the Arctic tern as a symbol of the wide range of subject areas I want to cover. I also like to think that as scientists we are in a constant chasing of the light, and in my field of climate research especially, trying to evade and chase off predators and pirates.

As a scientist, I often feel I am not very good at communicating with non-scientists and I hope this blog will provide me with some good practice and discipline, and of course I hope it will provide you, dear reader, with some interesting diversion. Please feel free to comment and provide feedback on my subjects and my writing style.