Can You Guess The Ingredients Of A Cake?

By: Amos Lawless

“Mmm this cake is lovely, what’s in it?” “Try to guess!” How often have we had that response from a friend or colleague who is proud of the cake they have just baked? And we usually try to guess the main ingredients – “I think there must be ginger or cinnamon. And can I taste lemon?”. But what if that friend persisted and asked you to try to guess all the ingredients – how many eggs they have used, how many grams of sugar are in the cake and how much butter it contains? Maybe you’d think they’d gone a bit crazy! Surely it is impossible to work out all the ingredients just by tasting it? It may sound unreasonable, but this is effectively what we try to do each day to interpret satellite measurements for our weather forecasts.

Weather satellites, besides giving us the nice pictures that we see on television, provide a wealth of other information about the atmosphere. Satellites actually measure the radiation emitted from the atmosphere at different frequencies, and these measurements depend on the properties of that part of the atmosphere that the satellite is looking at, such as its temperature, humidity and winds. It is as if these “ingredients” of the atmosphere are brought together into a “cake” that the satellite can taste. But what we are really interested in knowing is these ingredients. So how can we split the satellite measurement back into its atmospheric ingredients?

Thankfully we have a mathematical technique for doing this, which we call data assimilation. Each satellite instrument can measure at many different frequencies (as if they have many “taste buds” sensitive to different ingredients), so by combining measurements from different satellites in an intelligent way, as well as other more conventional measurements made on the ground, data assimilation helps us to build up a complete picture of the atmosphere all around the globe. This is done every day as part of modern weather forecasting, since knowing what the atmosphere is like now is essential if we are to make accurate forecasts. Most data assimilation techniques are based on finding an optimal combination of what we think is the current state of the atmosphere and our measurements, taking into account the precision of the different pieces of information we have. Writing down the theory of how to do this is fairly easy, but putting into practice is usually much harder.

Scientists of the Data Assimilation Research Centre (DARC) at the University of Reading work on a variety of problems related to data assimilation, from developing new approaches to applying it in practice. Each year, jointly with the National Centre for Earth Observation (NCEO), we organise a training course for scientists round the world to learn about the theory of data assimilation and how to apply it in practice. Lectures from DARC scientists are combined with computer practical exercises, so that participants can learn the theory of data assimilation and get a feel for how different methods perform in practice. Normally the course is held in-person, but this year there was the challenge of whether it was possible to hold it online. So it was that at the start of May our first ever training course on data assimilation using Microsoft Teams took place. Joining were 29 scientists from the UK, Belgium, Bulgaria, Denmark, Germany, Greece, Italy, Spain and the USA, working in universities, research institutes and meteorological forecasting centres.

Figure 1: Lecture by DARC scientist Dr Javier Amezcua

So how did we do it? By now we are already used to giving and listening to talks online, so the lecture part of the course was fairly straightforward. However, an important aspect of a course such as this is that it is interactive, with the possibility to ask questions. Thankfully the chat function worked well here, with participants putting questions in the chat continually and other DARC scientists responding if it wasn’t necessary to interrupt the lecture. Then computer practical exercises took place in breakout rooms, with groups of three participants working together. And during the breaks informal discussions took place using Gather.Town (a very impressive tool that I have only just discovered), including use of a virtual whiteboard to discuss further the mathematics. So what did the participants say about the online delivery? Comments included “I think the format worked really well”, “the arrangements for the remote delivery of the course were excellent”, “I think the practicals were organised well with lecturers rotating and coming to different rooms. That made me feel like I was in a classroom with having constant access to help”. Running this course certainly taught us a lot about how to teach data assimilation online, with lots of lessons learnt for the future. But everybody also realised that there are limitations to such a format. Hopefully next year we will be able to run the course in person again, with the opportunity for more informal discussions over coffee … and plenty of cake!

Figure 2: Online group photo of some of the lecturers and participants.

References

Data Assimilation Research Centre (n.d.), What is data assimilation? https://research.reading.ac.uk/met-darc/aboutus/what-is-data-assimilation/

Data Assimilation Research Centre (2019). Online lecture notes from 2019 training course.
https://research.reading.ac.uk/met-darc/training/ecmwf2019/

Lawless, A.S. (2013), Variational data assimilation for very large environmental problems. In Large Scale Inverse Problems: Computational Methods and Applications in the Earth Sciences (2013), Eds. Cullen, M.J.P., Freitag, M. A., Kindermann, S., Scheichl, R., Radon Series on Computational and Applied Mathematics 13. De Gruyter, pp. 55-90.

Nichols, N.K. (2009), Mathematical concepts of data assimilation. Preprint MPS_2009-04. Department of Mathematics, University of Reading.
http://www.reading.ac.uk/nmsruntime/saveasdialog.aspx?lID=48408&sID=90309

 

Posted in data assimilation, earth observation, Teaching & Learning | Leave a comment

Data Assimilation Improves Space Weather Forecasting Skill

By: Matthew Lang

Over the past few years, I have been working on using data assimilation methodologies that are prevalent in meteorology to improve forecasts of space weather events (Lang et al. 2017; Lang and Owens 2019). Data assimilation does this by incorporating observations from spacecraft orbiting the Sun into numerical solar wind models, allowing for estimates of the solar wind to be updated. These updated solar wind conditions are then used to drive a solar wind model that produces forecasts of the solar wind at Earth. I have shown that over the lifetime of the STEREO-B spacecraft (2007-2014), data assimilation is able to reduce errors in solar wind forecasts by about 31% compared to forecasts performed without (Lang et al. 2020). Furthermore, these data assimilated forecasts can compensate for systematic errors in forecasts produced from in-situ observations alone.

Space weather is the study of the changing environmental conditions in near-Earth space and its impacts on humans and our technologies, both in space and on Earth. One of the major drivers of space weather events is the solar wind, the constant outflow of plasma (the fourth state of matter that can be thought of as a hot, highly magnetised gas) from the Sun’s surface. The solar wind fills the solar system with particles and magnetic field and is constantly bombarding the Earth’s magnetic field.

Coronal Mass Ejections (CMEs) are huge eruptions of plasma from the Sun’s atmosphere that can travel from the Sun to Earth, through the solar wind, in as little as 18 hours and can drive the most severe space weather events. These include depletion of a part of the ionosphere that is responsible for bouncing radio signals around the planet, hence hampering long-distance communication systems.

Figure 1: Transformer damage from a CME that caused a blackout in Quebec in 1989.

Another major impact on human technologies is that the solar wind and CMEs drive changes to the Earth’s magnetic field, inducing electrical currents in the Earth’s atmosphere that have the potential to overload power systems causing transformer fires and widespread blackouts (this occurred in Quebec in 1989 (see Figure 1) and Sweden in 2003). Most of the impacts of a severe space weather event can be mitigated against if accurate forecasts are available. And that’s where data assimilation comes in.

Data assimilation is the combination of information from forecasts and observations of a system to produce an optimal estimate for the true state of that system. It is an invaluable tool in many aspects of modern life, with applications ranging from course correction during the Apollo Moon landing missions, satellite navigation in areas of poor GPS coverage and oil reservoir modelling. The most notable application for this blog, however, is its use in numerical weather prediction where it is a necessary step for producing more accurate starting points for weather forecasts. This reduces the impact of the “butterfly effect”, where a small change can lead to a vastly different outcomes in the future (the famous hypothetical example being that the titular butterfly flaps it’s wings in Japan leading to a tornado forming in the USA). By ensuring that weather forecasts are started as close to the truth as possible, the resulting forecasts will be more accurate over longer periods.

For consecutive 27-day periods (the time taken for the Sun to rotate once at its equator, relative to the Earth) between 2007 and the end of 2014, an empirical solar wind model called MAS (Magnetospherics Around a Sphere) (Linker et al. 1999) was used to generate a prediction of the solar wind conditions, which I call the prior solar wind. Data assimilation is then performed using data from the STEREO-A, STEREO-B and ACE spacecraft to generate a new set of solar wind conditions which I shall refer to as the posterior solar wind. The STEREO spacecraft orbit the Sun at approximately the same radial distance as Earth, however STEREO-A orbits at a rate of 22⁰ faster and STEREO-B 22⁰ slower. The ACE spacecraft is in near-Earth space, between the Earth and the Sun. The prior and posterior solar winds are then input into the simplified solar wind model, HUXt (Owens et al. 2020), which was developed at the University of Reading to produce forecasts for the subsequent 27-days. Finally, the prior and posterior forecasts were compared with a forecast from the STEREO-B spacecraft (the closest in-situ observation of the solar wind behind the Earth during this time), generated by assuming that the solar wind speed observed at STEREO-B will occur at Earth with a time-lag defined by the distance of the spacecraft behind Earth and the rotational speed of the Sun.

Figure 2: Plot showing the Root Mean Squared Errors (± one standard error) of the prior (blue), posterior(red) and STEREO-B corotation (orange) mean 27-day forecast of the solar wind speed over the lifetime of STEREO-B.

The results of these forecasts are summarised in Figure 2, where the mean 27-day solar wind speed forecast from the prior, posterior and STEREO-B corotation forecast are shown over the lifetime of STEREO-B. The posterior and corotation forecasts have lower Root Mean Squared Errors (RMSEs) than the prior forecasts, showing that both are good improvements over the prior forecast at all lead-times. It is understandable that the STEREO-B corotation and posterior forecast are similar, as
both use the observations from the STEREO-B spacecraft in their forecast.

Figure 3:  Solar wind speed forecasts using the HUXt mode where the Sun is in the centre and Earth is in the same location as the ACE spacecraft (black circle). The left one is initialised from the MAS empirical model without data assimilation and the right one is initialised from a data assimilation analysis, where STEREO (black triangles) and ACE observations have been assimilated. A coronal mass ejection (CME) initialised with the same characteristics and released from the Sun at the same time and propagated through the two ambient solar winds yielding very different evolutions of the CME.

A major difference between the STEREO-B corotation and the posterior forecast, however, is that the corotation produces a forecast at a single point as opposed to the posterior forecast which produces a forecast at every point in the model domain (at all radii and longitudes of interest). This is an especially useful feature as accurate specification of the solar wind can influence how CMEs evolve on their way to Earth. Figure 3 shows two CMEs initialised with the same properties (obtained from (Barnard et al. 2020)); the left one is propagated through a ‘prior’ solar wind that has not had data assimilation performed on it, compared to the ‘posterior’ on the right which does include data assimilation. The CME evolution is changed greatly by the different ambient solar winds, with the CME arriving 19 minutes earlier than it was observed at Earth with a data assimilated solar wind, compared to 41 hours late in a solar wind without data assimilation. By comparison, for the same CME, the operational solar wind model used by the Met Office arrived at Earth 10 hours before it was observed (Barnard et al. 2020). This shows that there is great potential in this field for data assimilation to improve forecasts of not only the solar wind, but also the more hazardous coronal mass ejections.

References

Barnard, L., M. J. Owens, C. J. Scott, and C. A. Koning, 2020: Ensemble CME Modeling Constrained by Heliospheric Imager Observations. AGU Adv., 1, https://doi.org/10.1029/2020av000214

Lang, M., and M. J. Owens, 2019: A Variational Approach to Data Assimilation in the Solar Wind. Sp. Weather, 17, 59–83, https://doi.org/10.1029/2018SW001857.

——, P. Browne, P. J. van Leeuwen, and M. Owens, 2017: Data Assimilation in the Solar Wind: Challenges and First Results. Sp. Weather, 15, 1490–1510, https://doi.org/10.1002/2017SW001681.

——, J. Witherington, H. Turner, M. Owens, and P. Riley, 2020: Improving solar wind forecasting using Data Assimilation. http://arxiv.org/abs/2012.06362.

Linker, J. A., and Coauthors, 1999: Magnetohydrodynamic modeling of the solar corona during Whole Sun Month. J. Geophys. Res. Sp. Phys., 104, 9809–9830, https://doi.org/10.1029/1998JA900159.

Owens, M., and Coauthors, 2020: A Computationally Efficient, Time-Dependent Model of the Solar Wind for Use as a Surrogate to Three-Dimensional Numerical Magnetohydrodynamic Simulations. Sol. Phys., 295, 43, https://doi.org/10.1007/s11207-020-01605-3.

Posted in Climate, data assimilation, space weather, Weather forecasting | Leave a comment

Cold Winter Weather: Despite or Because of Global Warming?

By: Marlene Kretschmer

This year’s winter was cold. There was heavy snowfall across the UK, Europe and parts of the United States including Texas. This severe weather came with significant societal and economic impacts.

Every time cold extremes like this occur, one can almost predict the media headlines.  On the one hand, dubious media will use a regional cold snap to sow doubt about human-made global warming by deliberately misunderstanding the difference between weather and climate. In a similarly absurd manner, other newspapers will state that climate change was responsible for the cold snap. In between, there are debates among scientists about the role of climate change in causing cold extremes. This is where it gets complicated and, hence, interesting.

Climate change manifests itself in different ways. While the increase of CO2 in the atmosphere leads to warmer temperatures globally, there might be indirect mechanisms causing opposite effects regionally. In recent years, researchers have hypothesised that the melting of Arctic sea ice – a direct result of global warming – favours winter cold extremes in the Northern Hemisphere mid-latitudes. In particular, it has been suggested that the decline in Barents and Kara sea ice weakens the stratospheric polar vortex, a band of fast blowing westerly winds circling the Arctic during winter at approximately 15-50 km altitude. Weak phases of the vortex are linked to cold winter weather in Eurasia and North America. In other words, it was proposed that climate change indirectly leads to colder weather. The polar vortex this year was extremely weak, and therefore likely to be the culprit of the cold weather. But are Arctic changes also making these weak vortex phases more likely?

Figure 1: Schematic overview of the different plausible causal mechanisms making it difficult to quantify the influence of autumn Barents and Kara sea ice concentrations (BK-SIC) on the winter stratospheric polar vortex (SPV); sea level pressure over the Ural Mountains (Ural-SLP) and over the North Pacific (NP-SLP), lower-stratospheric poleward eddy heat flux (vT), North Pacific sea ice concentrations (NP-SIC) and El Niño–Southern Oscillation/Madden–Julian Oscillation (ENSO/MJO). The arrows represent assumed causal relationships. (Taken from Kretschmer et al, 2020)

The scientific debate regarding a causal role of Arctic sea ice loss is controversial (see e.g. Cohen et al. 2020, Screen et al. 2018). Scientists face a dilemma. In observational data, a statistically significant signal has been detected. Given the large natural variations in climate data and different possible mechanisms which are difficult to disentangle, it is hard to tell if this signal reflects a causal influence (see also Fig. 1). This is further compounded by partly opposing results from climate model simulations. So far all that can be said conclusively, is that the question of whether the decline of Arctic sea ice is causing a weakening of the polar vortex cannot be answered conclusively.

But should we ignore the potential risk the decline of the Arctic holds for our future weather and climate, just because the current data do not allow a clear statement? The short answer is: No!

We explore this aspect in our latest study (see Kretschmer et al. 2020). In contrast to previous studies, which examined whether the decrease in sea ice causes a weakening of the polar vortex (and thus severe winter weather), we pose a different question. We ask: Assuming there is a causal influence of sea ice loss, what does this imply?

To address this question we use different climate model simulations of the next 100 years. All climate projections agree that sea ice will continue to melt as climate change progresses. This is a sad but unsurprising fact highlighting the need to evaluate possible consequences of a changing Arctic. Based on the model simulation data and using methods from causal inference, we further conclude that the causal effect of Arctic sea ice on the polar vortex is, if it exists, plausibly only very small. However, given that the decrease of sea ice will be huge, this small effect can have large implications. In fact, the climate models project a weakening of the polar vortex as long as the autumn sea ice in the Barents and Kara Sea melts. Whilst this is no definitive proof for a causal influence of sea ice loss, it is consistent with the initial hypothesis. Moreover, we find that once all sea ice is gone, the vortex strengthens again, suggesting there are other, poorly understood mechanisms by which global warming affects the polar vortex and thereby our weather in the mid-latitudes.

More generally, our study calls for more focus on understanding plausible climate-change related risks. Absolute statements about the regional effects of global warming are often not possible, given the complexity of the climate system and often contradictory climate predictions. This forces decision makers to act under large uncertainties. It is therefore necessary for climate scientists to evaluate different causal possibilities (such as an influence of the sea ice loss on the polar vortex) to gain a better understanding of regional climate risks. This also requires the use of different statistical tools and techniques – some of which we apply and discuss in our study.

The next time a cold snap hits Europe the same oversimplistic media headlines can be expected. Hopefully, however, the scientific debate will then have shifted towards a more conditional risk-based understanding of the plausible impacts of the changing Arctic.

References:

Cohen, J., Zhang, X., Francis, J. et al. Divergent consensuses on Arctic amplification influence on midlatitude severe winter weather. Nat. Clim. Chang. 10, 20–29 (2020). https://doi.org/10.1038/s41558-019-0662-y

Screen, J.A., Deser, C., Smith, D.M. et al. Consistency and discrepancy in the atmospheric response to Arctic sea-ice loss across climate models. Nature Geosci 11, 155–163 (2018). https://doi.org/10.1038/s41561-018-0059-y

Kretschmer, M., Zappa, G., and Shepherd, T. G.: The role of Barents–Kara sea ice loss in projected polar vortex changes, Weather Clim. Dynam., 1, 715–730.https://doi.org/10.5194/wcd-1-715-2020

 

Posted in Arctic, Climate, Climate change, Cryosphere, Polar | Leave a comment

What Did You Get For Number 9?

By: Todd Jones

A common way to check your work in school is to turn to your neighbour and ask, “What did you get for this one?”  With a little extra effort, though, students end up having productive discussions and learning to solve problems they didn’t fully understand or discovering new, clearer routes to the solution.  Even answering broadly defined questions, a group of comparisons could lead to consensus or a narrowing of the possibilities.

Particularly effective teachers encourage and schedule these comparison sessions. Scientists, ever the continual students, bring this technique to their research, aspiring to uncover solutions to challenging problems they have not previously considered by comparing their research with that of others. 

While the methods of solving 23÷1.4 with pen and paper varies little, questions about the motions of the atmosphere are not always so well constrained.  The many sensitive equations that describe these motions can often only be solved approximately, and scientists may reasonably choose from a number of approximations (with varying levels of accuracy) based on practical issues, such as how powerful their computers are.  For example, these calculations are often very large problems that require division of the atmosphere to be divided into a number of points where calculations about temperature, wind, and rain can be performed. Between the points, these parameters must be approximated with something like a “best guess.”  Each of these justifiable choices will lead to differing solutions that can generate years of classroom-style “compare and discuss” activity.

For example, we could compare solutions to simple models of the atmosphere.  One can remove complications of the real world and create close approximations that allow an easier solution.  For instance, picture a non-orbiting, non-rotating world that is entirely warm ocean, where the oscillations of night and day replaced by constant moderate sunshine.  The “world” doesn’t even have to be a sphere!  Modelling the atmosphere of this world, we would see that the atmosphere cools off gradually, radiating energy to space.  As the lower atmosphere warms from the ocean’s heat, moist convective bubbles begin to rise to then cool, forming clouds and rain.  Over time, the heat from condensation of water vapour into clouds and rain balances out the radiative cooling of the air.  We call this energetic balancing “radiative-convective equilibrium,” or RCE.  This model is a close approximation to Earth’s climate, and it can be used as a “toy Earth” to learn how the climate might change in response to parameter changes. 

Figure 1.  A scattered deep convective cloud scene from a simulation of the climate of a simplified world in radiative-convective equilibrium with an ocean constantly at 22°C in the UK Met Office model.  The back left wall shows a slice of relative humidity (hur). The back right wall shows a slice of specific humidity (water vapour concentration, hus).  The bottom surface shows the total amount of water vapour in the column of air above each point (prw).  Cloud surfaces are coloured for various levels of frozen cloud particles (cli), liquid cloud droplets (clw), and rain (plw).  Orange arrows show the velocity of the wind near the model surface.

Playing with these models over the past few decades, scientists have noticed some intriguing behaviour.  Choosing different global temperatures, we can investigate how clouds respond to global warming: will more reflective clouds spread and counter the warming?  Much of the time, the deep convective clouds that are generated in these models appear as one might guess: randomly scatted, sputtering across the little world (Figure 1).  However, when oceans are warm enough or the modeled worlds are sufficiently large, the deep convective clouds can spontaneously cluster into isolated locations, with very dry regions in between (Figure 2).  News of this phenomenon spurred tens of independent studies for comparison [1], and scientists began to uncover that phenomena like interactions between radiation and clouds can lead to this convective clustering.

Figure 2.  A clustered deep convective cloud scene from a simulation of the climate of a simplified world in radiative-convective equilibrium with an ocean constantly at 32°C in the UK Met Office model.  Features are as described in Figure 1.

For a fair test, the parameters used by each model should be the same, so a group of scientists gathered to define a specific set of parameters to test warming in RCE climates in models of many geometries and scales.  These “rules” were codified and shared [2], and volunteers reported their solutions for comparison [3].  Formally, this is an intercomparison of models simulating RCE, known as RCEMIP.  The 30+ representations varied between areas on the scale of 100 km to the full globe and between levels of detail (resolution) from 200 m to 50 km.

Though there were many small differences between the model results, there was broad agreement over the formation of aggregating clusters. Over small areas, only one model (that in Figure 2) developed convective clusters, whereas over large areas, all but a few models developed convective clusters.  The deep clouds in most models showed that, as the world is warmed, anvil tops become warmer, are located higher in the atmosphere, and cover smaller areas.  This means that the effect of high cloud tops on climate would vary little under global warming.  Instead, changes in low cloud properties and the degree of convective clustering can influence this [4].

Compared to high-resolution models, lower resolution global models show a change in clustering with global warming that indicates a smaller amount of warming for a given greenhouse gas forcing.  Because the higher resolution models tend to be more correct, it’s possible that coarser climate models have painted too rosy a picture about future warming. 

Though there is disagreement, there is much to be said for comparing solutions.  Many investigations comparing model patterns are underway, ultimately steering toward a better-understood solution of the climate system problem.

References:

[1] Wing, A. A., K. Emanuel, C. E. Holloway, and C. Muller, 2017: Convective self-aggregation in numerical simulations: A review. Surveys in Geophysics, 38 (6), 1173–1197, doi:10.1007/ s10712-017-9408-4, URL https://doi.org/10.1007/s10712-017-9408-4.

[2] Wing, A. A., K. A. Reed, M. Satoh, B. Stevens, S. Bony, and T. Ohno, 2018: Radiative– convective equilibrium model intercomparison project. Geoscientific Model Development, 11 (2), 793–813, doi:10.5194/gmd-11-793-2018, URL https://gmd.copernicus.org/articles/11/ 793/2018/.

[3] Wing, A. A., and Coauthors, 2020: Clouds and convective self-aggregation in a multimodel ensemble of radiative-convective equilibrium simulations. Journal of Advances in Modeling Earth Systems, 12 (9), e2020MS002138, doi:https://doi.org/10.1029/2020MS002138, URL https://doi.org/10.1029/2020MS002138.

[4] Becker, T., and A. A. Wing, 2020: Understanding the extreme spread in climate sensitivity within the radiative-convective equilibrium model intercomparison project. Journal of Advances in Modeling Earth Systems, 12 (10), e2020MS002 165, doi:https://doi.org/10.1029/2020MS002165, URL https://doi.org/10.1029/2020MS002165.

Posted in Climate, Convection | Leave a comment

April Flowers – A story of bluebells and frosts

By: Pete Inness

Figure 1: Bluebells in a wood near Reading on the 16th of April 2020 (left) and the same date in 2021 (right). In 2021 the flowers are yet to emerge and there are no leaves on the trees.

Bluebells regularly come out top in surveys of Britain’s favourite wildflower. From mid-April to mid-May they form carpets of lilac-coloured and strongly scented flowers in woodlands from Southern England to Scotland. Reading is particularly well placed for seeing these flowers with the beech woods of the Chilterns to the north of Reading being a favoured location for them. In fact, you don’t even need to leave the University campus as there are several good locations for them within a couple of minutes’ walk from the Meteorology Department.

Since I first came to Reading as a student in the late 80’s I’ve tried to get out into the countryside most years to see the bluebells at their best. This involves careful timing. Back in those early days I would have said that the May Day Bank Holiday weekend was the best time to catch them, but in the intervening 30 years that date has moved forward somewhat and I’d now say that going out a week or so earlier than this gives you a better chance of seeing them in their prime.

The main cause of year-to-year variations in the flowering date of bluebells is temperature variability. Whilst, like most woodland flowers, they are primed to get through much of their above-ground life cycle before the leaf canopy gets too thick and cuts down the sunlight reaching the forest floor, flowering can be accelerated or delayed by warmer or colder temperatures through the Spring months.

A few years ago, I decided to turn my interest in bluebells, and the annual cycle of nature in general, into something more productive by running undergraduate projects looking at relationships between weather patterns and the occurrence of events in the natural world.  This has been made possible by an excellent citizen science project called Nature’s Calendar which is run jointly by the Centre for Ecology and Hydrology and the Woodland Trust. This project encourages members of the public to report their sightings of a wide range of natural events such as first flowering of flowers and shrubs, first nest building of common birds, or first appearance of certain species of butterfly and other insects. Using the data recorded by this project, together with weather data such as the Met Office’s Central England Temperature record, students can explore relationships between weather and the annual cycle of the natural world and then relate them to specific weather events such as “the Beast from the East” in 2018 or longer-term changes in climate.

These studies by our students have shown that the flowering date of bluebells is sensitive to the average temperature through February and March – the months when the leaves emerge from the ground and the flower stalks and buds form. Every 1 degree Celsius rise in mean temperature across these months leads to bluebells flowering about 5 days earlier. The average temperature in April seems to have very little impact on the flowering date and this makes sense. Because bluebells produce their first flowers in mid-April (earlier in sheltered spots and in the south of the country), the temperature of the remainder of the month is immaterial to the flowering date.

2021 seems to be an exception to that rule. Whilst February 2021 was quite a bit colder than 2020, March 2021 was actually warmer than 2020. The differences in temperature between these 2 years in February and March are nowhere near large enough to explain the difference in the state of the bluebells in the pictures above, both taken on the 16th of April, in 2020 and 2021 respectively. April 2021 has been one of the coolest Aprils in recent years and in Reading has been the frostiest April since 1917. There have been 11 air frosts recorded at our Atmospheric Observatory through the month and only 5 nights in the month when there wasn’t a ground frost. To put these numbers in context, in a typical April in Reading we would expect 2 air frosts.

These frosts effectively slammed the brakes on the flowering process. Bluebells have evolved to avoid exposing the delicate reproductive parts of their flowers to frost and so during the first half of April the flower buds remained closed.  Even now, at the start of May, the bluebells in our local area are still some way behind the dense carpets of flowers that we saw in mid-April last year.

So, this year’s project students will be studying the effect of these exceptional frosts on UK wildlife and looking for their impacts on other plants, trees, birds and insects.

Posted in Climate, Phenology | Leave a comment

Some thoughts on future energy supply, such as an “Instantaneous Energy Market”

By: Peter Cook

We all know that it’s time to stop using fossil fuels, due to the greenhouse gasses emitted and the finite amount of these fuels.  Many renewable sources of energy are now being adopted but a lot of work and ingenuity will be needed for these to become the only sources of energy, and most people will need to be involved to make this happen.

A very different energy grid will be needed with multiple supplies (see figure), instead of the few large power companies at present, plus a lot of storage rather than just the National Grid balancing the load.  However, this should be seen as an opportunity not a problem.

There will be many opportunities for small companies and individuals to get involved, by generating their own electricity to sell, or by storing energy for other people, or by using energy in more efficient ways.  This could encourage a new entrepreneurial society, speeding up the adoption of new technology and the transition from fossil fuels to renewable energy.

A possible way to create the new energy grid would be to set up an “Instantaneous Energy Market”

Sources of renewable energy are often criticised for being intermittent, and their widespread adoption is dismissed as impractical because of the problems in matching energy supply to demand. These critics claim we need large-scale energy storage or backup sources of energy.  But is this way of thinking correct?  What about matching the demand to the supply instead?

Like other products, electricity can be priced according to supply and demand, and in many places, electricity is already cheaper at night than during the day.  Many of us make use of this, charging our storage heaters and running our washing machines and dishwashers at night, but this has the potential to be taken much further.  Prices could be adjusted second by second according to the instantaneous supply and demand.  Many uses such as heating, water heating and charging do not need to be on continuously and could be stopped for short periods, if demand (and price) became particularly high, without causing much inconvenience.

To do this the electricity supply would need to include a signal to show the price.  At present, the mains alternating current frequency is usually 60 Hz, but the frequency is reduced if the supply is low, so small changes to the frequency could be used to show the price.  There could also be information on how the supply, demand and price are changing in the short term, which would be used to predict the price in the very near future (minutes) to help people manage the changing price.  While on longer timescales (days) there could be electricity price forecasts, which would depend on the weather (sun and wind for supply, extra demand in cold weather), problems with supply, and large demands (during popular TV shows), which people could use to plan their electricity use and so reduce costs.

People who generate their own electricity (e.g. from solar panels) could sell their excess power, using large batteries to store electricity when it’s cheap and then sell it when the price increases.  Others could just have a large battery to buy electricity cheap and sell dear.  With this control of electricity demand and supply, adding new sources of energy would be easier, and energy suppliers would have less need for backup sources.

With many people adjusting their demand according to price, changes would be smoothed and variations in the price kept to a minimum.  When electricity is cheap, the resulting increase in energy use would lead to a price rise, whilst when electricity is expensive, the resulting drop in demand and increased electricity supply from people selling their own electricity would lead to a price fall.  People would also set their own thresholds of when to use electricity or not so that abrupt jumps in the overall demand would be avoided.  Attempts at profiteering (storing energy to raise the price) would be difficult because of the large amount of storage that would be needed.

The use of instantaneous energy pricing might work better at a local rather than at a national level, and modelling studies are required to see how it would work in practice, identify potential problems and to investigate the extent to which such a system could be scaled up.

Reference

The attached figure (but not any of the above text) is from the paper “Smart management system for improving the reliability and availability of substations in smart grid with distributed generation”, by Shady S. Refaat and Amira Mohamed, January 2019, The Journal of Engineering (17), DOI:10.1049/joe.2018.8215

 

Posted in Climate, Energy meteorology | Leave a comment

TerraMaris: Plans, Progress And Setbacks Of Atmospheric Research In Indonesia

By: Emma Howard 

To some of us weather enthusiasts, there’s nothing more exciting than a good tropical thunderstorm. For the best storms, you need a good source of humid air from a warm ocean and a hot land surface. If you can find some mountains to push air upwards and initiate convection (the intense vertical motion of air in updrafts and downdrafts which drive storms) all the better.

Figure 1: Development of convection offshore of West Papua. Photo credit: Megan Howard 

As a volcanic archipelago centred right on the equator, Indonesia has all of this and more. So it’s no surprise that Indonesia is the largest of the three major tropical convective hotspots on Earth. Local lore says that rain comes like clockwork during the wet season, occurring every day at the same time for weeks on end. This is borne out in quantitative rainfall observations, which show that after forming over mountains and land during the mid-afternoon and evening, storms tend to move offshore, with regular night-time and early-morning showers over the oceans and seas adjacent to islands. At present, most atmospheric forecast models (which parameterise atmospheric convection rather than resolving it) don’t represent these diurnally propagating systems very well. This makes it challenging to use these models to predict the timing and intensity of convection in Indonesia.

Unfortunately, some of the more intense thunderstorms can have severe impacts on local communities, particularly when associated with large-scale forcing such as Tropical Cyclone Seroja, which struck Timor-Leste and the Indonesian Nusa Tenggara provinces just two weeks ago. Beyond their immediate impact, these storms have subtle impacts further afield. By condensing water vapour into ice and liquid miles above the earth’s surface, intense storms cause latent heat to be released in the upper atmosphere. This heat source drives the Hadley and Walker cells, global scale atmospheric circulation systems which influence weather and climate across the world, including the UK. For these reasons, scientific research into the convection that occurs in thunderstorms in Indonesia is critical for our understanding of the Earth system and improving climate models.

TerraMaris is a large, collaborative research project that is furthering scientific understanding of atmospheric convection in the Indonesian region. The project involves researchers from three UK universities (East Anglia, Reading and Leeds), the UK Met Office and Indonesia’s weather and space agencies (BMKG and LAPAN). TerraMaris aims to transform our understanding of convective processes in Indonesia and their interactions with the largescale flow through an intensive observational and modelling campaign focussed on the circulation systems associated with the daily development and offshore propagation of convection.

Thankfully, the modelling component of our project hasn’t been so affected by the pandemic and is chugging away as normal. We’re generating a set of very high-resolution model simulations over the whole of Indonesia that are able to (at least partially) resolve the convective updrafts and downdrafts in the daily-repeating storms. Unlike many lower resolution models, these simulations are capable of accurately simulating offshore propagating convection. We intend to run 10 simulations, with one during the long-awaited field campaign season, covering the entire December – February rainy season. A wide range of weather conditions will be represented in this sample, and we’ll be able to study the simulated thunderstorms during all of them.

We are able to compare the role these storms play in heating the upper atmosphere to that in more conventional, lower resolution models, which aren’t able to resolve the updrafts and downdrafts and instead have to parameterise them. These models generally don’t represent Indonesian convection very well. It’s early days, but we’re finding that there’s a lot more variability in the height above the ground where heating occurs in the high-resolution models than the low resolution models. Our high-resolution models also simulate the daily formation of storms in the afternoon/evening and their overnight propagation into the oceans really well (see video).

Figure 2: Mean diurnal cycle of precipitation in early TerraMaris simulations.

Because interactions between the atmosphere and the warm tropical oceans are really important in this part of the world, we’re using a carefully designed coupled atmosphere-ocean model to run all these simulations. Full ocean models are very computationally expensive to run, so we’re using a multi-column KPP ocean model in order to simulate turbulent vertical mixing in the near-surface mixed layer. This is the oceanic process that interacts most strongly with the atmosphere, as it transports heat and freshwater fluxes from the atmosphere at the sea surface further down through the upper ocean. The role of ocean currents and other processes are represented by imposing “corrective” sources and sinks of heat and salt which ensure that in the long run, our simulated ocean matches up with observations of the real ocean.

We’re hoping that these simulations will be able to answer some really fundamental questions about how large-scale weather conditions modulate the vertical distribution of convective heating and how important the daily propagating systems are for providing the heat that drives global circulation. This will be useful for improving the representation of Indonesian convection in lower resolution models. If we can improve that, we hope that weather forecasts will improve both locally in Indonesia and globally through interactions with the Hadley and Walker cells. With any luck, by the time we finally step onto that plane, we’ll know a lot more about the storms that we’re trying to observe than we do now!

 

Posted in Atmospheric circulation, Climate, Convection, Rainfall, Thunder Storms, Tropical convection | Leave a comment

Pacific and Atlantic Conversations

By: Daniel Hodson

The Earth is a world of water – oceans spread out across much of the planet and they exert a profound influence over the climate. Ascending from the Earth, the churning waves and surf shrink away and the oceans relax into seemingly silent, passive bodies of water. But this seeming passivity belies a complex network of currents and flows hidden beneath the surface, driven by heat at the equator flowing to the colder poles, but being frustrated in doing so by the spin of the Earth.

Figure 1: The Atlantic Meridional Overturning Circulation

In the Atlantic, an immense flow of water drives northwards towards Greenland and Iceland in the top kilometre of ocean, before plunging down kilometres and returning southwards at depth, towards Antarctica (Figure 1). This is the deep Atlantic Meridional Overturning Circulation (AMOC). This circulation involves such large flows of water that oceanographers had to invent a new unit of measurement to think about the volumes involved: the Sverdrup (Sv) is a million metres cubed per second – that’s a cube of water 100m on a side, flowing past every second. This northward flowing water carries heat with it, sometimes speeding up, sometimes slowing down – bringing more or less heat as it does so, leading to warming or a cooling of the surface of the ocean. This heat can then be carried away by the atmosphere leading to warmer air temperature, or perhaps driving changes in surface wind patterns.

If, whilst orbiting over the Pacific, you tuned your eyes away from the blue of the Pacific and into the infrared, you would see what the satellites see: a vast pattern of warm and cold spread out across the expanse of the Pacific Ocean. Over the years, you would see this pattern pulse warm and then cold; in the semi-regular cycle of El Nino: the heartbeat of the climate system which dominates the tropics.Figure 2: The Pacific Decadal Oscillation pattern

El Nino is driven by complex interactions between the winds blowing over the Pacific Ocean, and the waters sloshing between Asia and the Americas. It leads to a 3-6-year cycle of warming and cooling in the equatorial Pacific Ocean. In the warm phase, large pulses of heat are released from the ocean into the atmosphere, shifting climate patterns leading to droughts and deluges across the globe. Over many decades of watching, a more widespread pattern of warming and cooling emerges across the Pacific – a pattern known as the Pacific Decadal Oscillation (PDO) (Figure 2). The connection between the PDO and El Nino remains to be fully understood.

Both the AMOC and the PDO play a key role in storing and moving heat around; their variations over time, in turn, modulate our climate system, potentially in profound ways. The way these climate features respond to external factors like changing levels of greenhouse gases or industrial pollution may affect the medium-term trajectory of anthropogenic climate change.Figure 3: The Pacific and Atlantic Oceans

Figure 4: The Tropical Walker Circulation

For a long time, it was thought that these two siblings (AMOC and PDO) continued their existence in ignorance of the other; bounded by Africa and Eurasia but divided by the Americas (Figure 3). They may hear distant echoes of each other, mediated by the turbulent Southern Ocean around Antarctica, or the icy Arctic Ocean – but signals in the ocean are ponderous, slow and noisy. New simulations with modern complex climate models suggest that they hear and feel each other’s presence over, rather than around, the wall of the Americas; mediated by the atmosphere. The Walker circulation is the large-scale pattern of ascending and descending air one encounters when travelling around the equator (Figure 4). Air heated and pushed upwards by a warm ocean in one place, must be replaced by descending air elsewhere in the tropics. This circulation seems to allow the two oceans to talk to and influence each other. Climate model simulations [1][2] seem to show that, over many decades, a warmer Atlantic can nudge a cooler Pacific Ocean, whilst a Warmer Pacific ocean can lead to a warmer Atlantic.

Whilst we are seeing a clearer picture of how these two oceans coordinate their climate modulations, challenges remain. Many decades of observations are needed to understand the slow influences of these twin oceans – but whilst the 21st-century ocean is well observed, ocean observations before 1950 are much scarcer. Remarkable efforts are underway, however, to utilize the vast datasets buried in old ships logs. We also rely on climate models to tease apart the complex interactions in the climate system. Are the models we use accurate enough? Are we doing the right experiments with these models to understand how these features of climate interact? If we can begin to understand the conversation between these two oceans better, we may be better able to predict their future influences on climate and, in turn, on us.

References:

[1] Meehl, G.A., and Coauthors, 2021: Atlantic and Pacific tropics connected by mutually interactive decadal-timescale processes. Nat. Geosci. 14, 36–42 . https://doi.org/10.1038/s41561-020-00669-x

[2] Ruprich-Robert, Y., Msadek, R., Castruccio, F., Yeager, S., Delworth, T., & Danabasoglu, G., 2017: Assessing the Climate Impacts of the Observed Atlantic Multidecadal Variability Using the GFDL CM2.1 and NCAR CESM1 Global Coupled Models, Journal of Climate, 30(8), 2785-2810. https://doi.org/10.1175/JCLI-D-16-0127.1

Posted in Climate | Leave a comment

Satellite data used to provide life-saving weather forecasts in tropical Africa

By: Peter Hill

Much of the population of tropical Africa are vulnerable to severe weather, often caused by intense storms that can generate heavy rainfall, strong winds and flooding. For instance, thousands of fishermen drown each year in Lake Victoria as a result of accidents caused by storms. As a result, improved weather forecasting systems in tropical Africa could save lives and protect livelihoods.

The Global Challenges Research Fund (GCRF) African Science for Weather Information and Forecasting Techniques (SWIFT) project aims to enable African weather forecasting services to develop such improved weather forecasting systems. A partnership between meteorologists from Senegal, Ghana, Nigeria, Kenya and the UK, including several scientists at the University of Reading, SWIFT is striving to improve forecasts from timescales of a few hours to a few weeks ahead.

Much of my work in the SWIFT project involves very short-range predictions – from 0 to 12 hours ahead – based directly on observations, something meteorologists term “nowcasting”. The simplest nowcasts take weather observations and extrapolate them forwards in time, using the assumption that the weather will continue to develop along the same trajectory as the recent past. Nowcasts can be crucial for severe weather events, providing timely information to enable authorities and the public to respond appropriately to safeguard lives and livelihoods.

One of the major obstacles to nowcasting in tropical Africa is the lack of rainfall radar observations, which are used for nowcasting in other parts of the world, including the UK. Passive satellite observations, which measure the naturally occurring energy at the top of the atmosphere, provide a less direct measure of weather systems. Yet in the absence of other observations, this satellite data can provide vital information for nowcasting purposes.

To this end, the SWIFT project has made satellite-based nowcasts for tropical Africa freely available from a new website. Figure 1 provides examples of two such products. These nowcasts are based on software provided by the European Nowcasting Satellite Applications Facility (NWCSAF). However, these products have been calibrated and validated for mid-latitude European weather systems and it is therefore necessary to evaluate how well they perform for tropical Africa.

Figure 1: Examples of two NWCSAF products over tropical Africa. (a) shows the convective rainfall rate in different regions (b) shows the rapidly developing thunderstorms convection-warning product over the Guinea Coast region.

To understand the suitability of this NWCSAF software for tropical Africa, I compared the two products shown in Figure 1 to higher quality satellite rainfall estimates that incorporate data from multiple sources including direct rainfall estimates from rain gauges at the surface. This higher quality data cannot be used for nowcasting because it is not available sufficiently quickly.

The comparison demonstrates that both NWCSAF products provide useful information, despite some limitations. For instance, the convective rain rate product has valuable skill for predictions at least 90 minutes ahead (Figure 2). The rapidly developing thunderstorms product can also identify the occurrence of heavy precipitation, correctly identifying around 60% of strong (5 mm of rain per hour) events at least one hour before they occur. These products could be used to inform flood warnings, disaster response, or provide warnings to fishermen.

Figure 2: Skill metrics for the convective rainfall rate product, compared to predictions based on the historical occurrence of rainfall events. Hit rate means the proportion of true rainfall events that are successfully identified, and false alarm ratio means the proportion of the predicted rainfall events that do not occur in reality. “Retrieval” here is the skill of the satellite products, versus higher quality data. This higher quality data is regarded as “truth” but is not available sufficiently quickly to be useful for nowcasts. The “extrapolation” is the forecast made by projecting the observed storms forward in time. The “climatology” skill is skill from assuming today’s storms can be predicted using previous years storms on the same time of day and time of year.

This analysis is crucial in providing forecasters with confidence in the products which GCRF African SWIFT has made available to them to issue warnings. It has also highlighted some aspects of both the convective rain rate and rapidly developing thunderstorm – convection warning products that could be improved upon. Future work will aim to further develop these products to provide better nowcasts for tropical Africa.

Ongoing work within the African SWIFT project is also enabling African groups to generate these products locally, as well as supporting forecasters to understand and use these products effectively to minimise adverse impacts of severe weather on lives and livelihoods in Africa.

References:

Roberts, A.J., Fletcher, J.K., Groves, J., Marsham, J.H., Parker, D.J., Blyth, A.M., Adefisan, E.A., Ajayi, V.O., Barrette, R., de Coning, E., Dione, C., Diop, A.-L., Foamouhoue, A.K., Gijben, M., Hill, P.G., Lawal, K.A., Mutemi, J., Padi, M., Popoola, T.I., Rípodas, P., Stein, T.H.M., Woodhams, B.J. 2021. Nowcasting for Africa: advances, potential and value. Weather (In press).

Hill, P. G., Stein, T. H. M., Roberts, A. J., Fletcher, J. K., Marsham, J. H. & Groves, J. 2020. How skilful are Nowcasting Satellite Applications Facility products for tropical Africa? Meteorological Applications, 27(6). DOI: https://doi.org/10.1002/met.1966

Posted in Africa, Climate, Rainfall, Remote sensing | Leave a comment

Flood forecasting for the Negro River in the Amazon Basin

By: Amulya Chevuturi

Figure 1: Photograph of the Negro River and the Amazon rainforest.

The Amazon is the largest river basin in the world, with large free-flowing rivers, draining about one-sixth of global freshwater to the ocean. The Amazonian floodplains have been long settled and used by indigenous populations, providing essential ecosystem services and natural resources for human needs (Junk et al., 2014). Increasing frequency and magnitude of floods in the last two decades has caused considerable environmental and socio-economic losses in many regions of the Amazon basin (Marengo and Espinoza, 2016). Although some studies have estimated flood risk for the Amazon basin (de Andrade et al., 2017), most towns and cities in this region still lack operational flood forecasts and integrated flood risk management plans.

The main aim of the PEACFLOW (Predicting the Evolution of the Amazon Catchment to Forecast the Level Of Water) project was to develop skilful forecasting systems for high water levels of Amazonian rivers, at sufficiently long lead time, for effective implementation of disaster risk management actions. In this project, we focused on developing forecast models for annual maximum water level for the Negro River at Manaus, Brazil (Figure 1), as a pilot case study, using a multiple linear regression approach. We used various potential predictors from preceding months: rainfall, water level, Pacific and Atlantic Ocean conditions and linear trend, all of which strongly influence the water levels in the Amazon basin. Flood levels in the Negro River occur between May and July and are strongly influenced by the rainfall during November to February, as its large floodplains delay the flood wave by months (Schöngart and Junk, 2007). This delay and the regularity between the rainfall and peak water level allows for the development of skilful statistical forecast models that can issue forecasts by March or earlier.

Figure 2: The Negro, Solimões and Madeira Rivers (blue lines) and their catchment basins (regions bounded by black lines) contributing to the river water level at Manaus (yellow circle; 3.14°S, 60.03°W).

In collaboration with Brazilian scientists, from various partner institutes, our team developed forecast models of the annual maximum water level (flood level) for the Negro River at Manaus by finding the best model fit over the training period of 1903 to 2004. For our models, rainfall over the catchment of the Negro River as well as from the catchment of nearby Solimões and Madeira Rivers (Figure 2) is the predominant predictor. We developed three models in this project, which use observations as input and can be implemented operationally to provide flood forecasts for Manaus. We compared the models developed in this project against current operational forecasts, provided by Brazilian agencies (CPRM and INPA), for the period of 2005 to 2019. The three PEACFLOW models issue forecasts of flood levels in the middle of January, February and March each year, but the skill of the models increase with decreasing lead-time (Figure 3a). Our results show that the models developed in this study can provide forecasts with the same skill as existing operational models one month in advance.

We also gained an additional month of lead-time when we relaced the observed input data with the ECMWF seasonal ensemble forecast. We developed two operational models using this data, which provide probabilistic forecasts at the beginning of January and February (Figure 3b). The probabilistic forecasts for the maximum water level, using ECMWF input, show good skill for extreme flood likelihood.

Figure 3: Comparison of models developed in PEACFLOW project and existing models (CPRM and INPA) and observed values for annual maximum water level at Manaus for models using (a) observations and (b) seasonal forecasts as input.

The methods developed in this project can also be used to develop forecast models for flood and drought levels over other regions of the Amazon basin. We provide the fully automated PEACFLOW models in a GitHub repository at https://github.com/achevuturi/PEACFLOW_Manaus-flood-forecasting. We retrospectively forecasted the annual maximum water levels at Manaus for 2020, and we are actively forecasting for 2021 (Table 1). Our forecasts this year show the maximum water levels crossing 29m, which is the extreme flood threshold for Manaus, at which the government declares emergency conditions.

Table 1: Forecasts for 2020 and 2021 using PEACFLOW models at different lead-times. Observed annual maximum water level at Manaus for 2020 was 28.52m.

 

References:

de Andrade MMN et al. (2017) Flood risk mapping in the Amazon. Flood Risk Management, 41. DOI:  https://doi.org/10.5772/intechopen.68912

Junk WJ et al. (2014) Brazilian wetlands: their definition, delineation, and classification for research, sustainable management, and protection. Aquatic Conservation: marine and freshwater ecosystems, 24, 5–22. DOI: https://doi.org/10.1002/aqc.2386

Marengo JA and Espinoza JC (2016) Extreme seasonal droughts and floods in Amazonia: causes, trends and impacts. International Journal of Climatology, 36, 1033–1050. DOI: https://doi.org/10.1002/joc.4420

Schöngart J and Junk WJ (2007) Forecasting the flood-pulse in Central Amazonia by ENSO-indices. Journal of Hydrology, 335(1),124–132. DOI: https://doi.org/10.1016/j.jhydrol.2006.11.005

Posted in Amazon, Climate, Flooding | Leave a comment