Characterising the seasonal and geographical variability in tropospheric ozone, stratospheric influence and recent changes

Email: r.s.williams@pgr.reading.ac.uk

Williams, R. S., Hegglin, M. I., Kerridge, B. J., Jöckel, P., Latter, B. G., and Plummer, D. A.: Characterising the seasonal and geographical variability in tropospheric ozone, stratospheric influence and recent changes, Atmos. Chem. Phys., 19, 3589–3620, https://doi.org/10.5194/acp-19-3589-2019, 2019.


Approximately 90 % of atmospheric ozone (O3) today resides in the stratosphere, which we know as the ozone layer (extending from ~15-35 km), where it plays a critical role in filtering out most of the harmful ultraviolet (UV) rays from the sun. The gradual formation of the ozone layer from around 600 million years ago was key in Earth’s evolutionary history, as it enabled life to flourish on land. Lesser known is the importance of the remaining ~ 10 % of atmospheric ozone, which is found in the troposphere and has implications for air quality, radiative forcing and the oxidation capacity of the troposphere. Whilst ozone is a pollutant at ground level, contributing to an estimated 6 million premature deaths globally per year, it also acts to cleanse the troposphere by breaking down a large number of pollutants, along with some greenhouse gases. Ozone is however a greenhouse gas in itself – where it has a maximum radiative forcing in the upper troposphere. It is an example of a non-well mixed gas, owing to its spatially and temporally highly varying sources and sinks, as well as its relatively short global mean tropospheric lifetime of about 3 weeks.

Figure 1 – Seasonal composites of monthly averaged 1000-450 hPa (0-5.5 km) subcolumn O3 (DU) for 2005-2010 (left to right) from (a) OMI, (b) EMAC minus OMI and (c) CMAM minus OMI. Circles denote (a) equivalent ozone-sonde derived subcolumn O3 (DU), (b) EMAC minus ozone-sonde differences and (c) CMAM minus ozone-sonde differences. All data were regridded to 2.5° resolution (~ 275 km). 1 Dobson Unit (DU) equates to a thickness of 0.01 mm if it were compressed at sea level.

A major source of tropospheric ozone is the photochemical reactions of emission precursors such as carbon monoxide (CO), nitrogen oxides (NOx) and volatile organic compounds (VOCs), which have both natural and anthropogenic sources, in addition to the natural influx of ozone-rich air from the stratosphere. The magnitude of these two competing influences has been poorly quantified until the recent advent of satellite observations and the development of comprehensive chemistry-climate models (CCMs), which simulate interactive chemistry and are stratospherically well-resolved.

Our study aimed to update and extend the knowledge of a previous key study (Lamarque et al., 1999), that investigated the role of stratosphere-troposphere exchange (STE) on tropospheric ozone, using two contemporary state-of-the-art CCMs (EMAC and CMAM) with stratospheric-tagged ozone tracers as a diagnostic. We first sought to validate the realism of the model ozone estimates with respect to satellite observations from the Ozone Monitoring Instrument (OMI), together with spatially and temporally limited vertical profile information provided from ozonesondes, which we resolved globally on a seasonal basis for the troposphere (1000-450 hPa) (Figure 1).

Whilst we found broad overall agreement with both sets of observations, an overall systematic bias in EMAC of + 2-8 DU (Dobson Units) and regionally and seasonally varying biases in CMAM (± 4 DU) can be seen in the respective difference panels (Figure 1b and 1c). A height-resolved comparison of the models with respect to regionally aggregated ozonesonde observations helped us to understand the origin of these model biases. We showed that apparent closer agreement in CMAM arises due to compensation of a low bias in photochemically produced ozone in the troposphere, resulting from the omission of a group of emission precursors in this model, by excessive smearing of ozone from the lower stratosphere due to an inherent high bias. This smearing is induced when accounting for the satellite observation geometry of OMI, necessary to ensure a direct comparison with vertically well-resolved models, which has limited vertical resolution due to its nadir field of view. The opposite was found to be the case in EMAC, with a high (low) bias in the troposphere (lower stratosphere) relative to ozonesondes. Given the similarity in the emission inventories used in both models, the high bias in this model indicates that excess in situ photochemical production from emission precursors is simulated within the interactive chemistry scheme. These findings emphasise the importance of understanding the origin of such biases, which can help prevent erroneous interpretations of subsequent model-based evaluations.

Noting these model biases, we next exploited the fine scale vertical resolution offered by the CCMs to investigate the regional and seasonal variability of the stratospheric influence. Analysis of the model stratospheric ozone (O3S) tracers revealed large differences in the burden of ozone in the extratropical upper troposphere-lower stratosphere (UTLS) region, with some 50-100 % more ozone in CMAM compared to EMAC. We postulated that CMAM must simulate a stronger lower branch of the Brewer-Dobson Circulation, the meridional stratospheric overturning circulation, since the stratospheric influence is isolated using these simulations. This has implications for the simulated magnitude and distribution of the downward flux of ozone from the stratosphere in each model. Shown in Figure 2 is the zonal-mean monthly evolution of ozone volume mixing ratio (ppbv) from ozonesondes and EMAC over the period 1980-2013 for the upper (350 hPa), middle (500 hPa) and lower (850 hPa) troposphere, together with the EMAC O3S and derived fraction of ozone of stratospheric origin (O3F) (%) evolution.

We found that the ozonesonde evolution closely resembles that of both EMAC and CMAM (not shown) throughout the troposphere. A clear correspondence in the seasonality of ozone is also evident for the EMAC O3S tracer, and in turn the O3F evolution, particularly towards the upper troposphere. Nonetheless, both models imply that over 50 % of near-surface ozone is derived from the stratosphere during wintertime in the extratropics, which is substantially greater than that estimated by Lamarque et al. (1999) (~ 10-20 %), and still considerably higher than more recent studies (~ 30-50 %) (e.g. Banarjee et al., 2016). This indicates that the stratospheric influence may indeed be larger than previously thought and is thus an important consideration when attempting to understand past, present and future trends in tropospheric ozone.

Figure 2 – Zonal-mean monthly mean evolution of ozone (O3) volume mixing ratio (ppbv) derived from (a) ozonesondes and (b) EMAC. The evolution of the (c) EMAC stratospheric ozone (O3S) tracer and (d) stratospheric fraction (O3F) (%) are additionally included over the period 1980-2010 for 350 hPa (top row), 500 hPa (middle row) and 850 hPa (bottom row).

Finally, we analysed height-resolved seasonal changes in both the model O3 and O3S between 1980-89 and 2001-10. The calculated hemispheric springtime (MAM/SON) changes in ozone are shown in Figure 3, and equivalently for O3S in Figure 4, for the upper and middle troposphere (350 and 500 hPa), as well as for the surface model level. A general increase in tropospheric ozone was found worldwide in all seasons, which is maximised overall during spring in both the Northern Hemisphere (~ 4-6 ppbv) and the Southern Hemisphere subtropics (~ 2-6 ppbv), corresponding to a relative increase of about 5-10 %. Respectively, a significant stratospheric contribution to this change of ~ 3-5 ppbv and ~ 1-4 ppbv is estimated using the model O3S tracers (~ 50-80 % of the total change), although with substantial inter-model disagreement over the magnitude and sometimes the sign of the attributable change for any given region or season from the stratosphere.

Figure 3 – Seasonal change in EMAC ozone volume mixing ratio (O3) (ppbv) between 1980-89 and 2001-10 for MAM (top) and SON (bottom) at (a) 350 hPa, (b) 500 hPa and (c) the surface model level. Stippling denotes regions of statistical significance according to a paired two-sided t-test (p < 0.05).

Although surface ozone changes are dominated by regional changes in precursor emissions between the two periods – the largest, statistically significant increases (> 6 ppbv) being over south-east Asia – the changing influence from the stratosphere were estimated to be up to 1–2 ppbv between the two periods in the Northern Hemisphere, albeit with high regional, seasonal and inter-model variability. In relative terms, the stratosphere can be seen to typically explain 25-30 % of the surface change over regions such as the Himalayas, although locally it may represent the dominant driver (> 50 %) where changes in emission precursors are negligible or even declining due to the enforcement of more stringent air quality regulations over regions such as western Europe and eastern North America in recent years.

Figure 4 – Seasonal change in EMAC stratospheric ozone volume mixing ratio (O3S) (ppbv) between 1980-89 and 2001-10 for MAM (top) and SON (bottom) at (a) 350 hPa, (b) 500 hPa and (c) the surface model level. Stippling denotes regions of statistical significance according to a paired two-sided t-test (p < 0.05). Note the scale difference between (a-b) and (c).

To summarise, our paper highlights some of the shortcomings of the EMAC and CMAM CCMs with respect to observations and we emphasise the importance of understanding model bias origins when performing subsequent model-based evaluations. Additionally, our evaluations highlight the necessity of a well-resolved stratosphere in models for quantifying the stratospheric influence on tropospheric ozone. We find evidence that the stratospheric influence may be larger than previously thought, compared with previous model-based studies, which is a highly significant finding for understanding tropospheric ozone trends.

References:
Lamarque, J. F., Hess, P. G. and Tie, X. X.: Three‐dimensional model study of the influence of stratosphere‐troposphere exchange and its distribution on tropospheric chemistry., J. Geophys. Res. Atmos., 104(D21), 26363-26372, https://doi:10.1029/1999JD900762, 1999.

Banerjee, A., Maycock, A. C., Archibald, A. T., Abraham, N. L., Telford, P., Braesicke, P., and Pyle, J. A.: Drivers of changes in stratospheric and tropospheric ozone between year 2000 and 2100., Atmos. Chem. Phys., 16, 2727-2746, https://doi.org/10.5194/acp-16-2727-2016, 2016.

The Colour of Climate

Email: Jake.J.Gristey@noaa.gov
Web: https://cires.colorado.edu/researcher/jake-j-gristey

Gristey, J.J., J.C. Chiu, R.J. Gurney, K.P. Shine, S. Havemann, J. Thelen, and P.G. Hill, 2019: Shortwave Spectral Radiative Signatures and Their Physical Controls. J. Climate, 32, 4805–4828, https://doi.org/10.1175/JCLI-D-18-0815.1

Sunlight reaching the Earth is comprised of many different colours, or wavelengths. Some of these wavelengths cannot be detected by the human eye, such as the ultraviolet (UV) wavelengths which famously cause sunburn. Fortunately for us, the most intense sunlight is found at harmless visible wavelengths and reaches the surface with relative ease, allowing us to see during the daytime. Sometimes nature aligns to dramatically separate these wavelengths, producing beautiful optical phenomena such as rainbows. More often, however, the properties of the atmosphere and surface lead to intricate differences in the wavelengths of sunlight that get reflected back to space (Fig. 1).

Fig. 1. Schematic showing how the spectral structure of reflected sunlight at the top of the atmosphere can emerge via interactions with various atmospheric/surface properties*.

Satellites have observed specific wavelengths of reflected sunlight to infer the properties and evolution of our climate system for decades. Satellites have also independently measured the total amount of reflected sunlight across all wavelengths to track energy flows into and out of the Earth system. It has been less common to make spectrally resolved measurements at many contiguous wavelengths throughout the solar spectrum. In theory, these measurements would simultaneously provide the total energy flow – by integrating over the wavelengths – and the “spectral signature” associated with all atmospheric and surface properties that determined this energy flow. Our recent study puts this theory to the test.

Almost 100,000 spectra of reflected sunlight were computed at the top-of-atmosphere under a diverse variety of conditions. Applying a clustering technique to the computed spectra (which identifies “clusters” in a dataset with similar characteristics) revealed distinct spectral signatures. When we examined the atmospheric and surface properties that were used to compute the spectra belonging to each spectral signature, a remarkable separation of physical properties was found (Fig. 2).

Fig. 2. (top row) Three of the extracted “spectral signatures” of reflected sunlight. (bottom row) Their relationship to the underlying atmospheric/surface properties. Seven others are shown in the published article.

Surprisingly, the separation of physical properties by distinct spectral signatures, as shown in Fig. 2, was found to be robust up to the largest spatial scales tested of 240 km. This is similar to the footprint size of one of the only previous satellite instruments to measure contiguous spectrally resolved reflected sunlight, the SCIAMACHY**, providing an exciting opportunity to investigate spectral signature variability in real observations. We found that the frequency of spectral signatures in real SCIAMACHY observations followed the expected behaviour during the West African monsoon very closely (Fig. 3).

Fig. 3. (left) The annual cycle of precipitation [mm/day] associated with the West African monsoon, and (right) frequency of the three “spectral signatures” shown in Fig. 2 from real satellite observations during 2010 over West Africa.

Overall, the separation of physical properties by distinct spectral signatures demonstrates great promise for monitoring evolution of the Earth system directly from spectral reflected sunlight in the future.

Funding acknowledgement: This work was supported by the Natural Environment Research Council (NERC) SCience of the Environment: Natural and Anthropogenic pRocesses, Impacts and Opportunities (SCENARIO) Doctoral Training Partnership (DTP), Grant NE/L002566/ 1, and from the European Union 7th Framework Programme under Grant Agreement 603502 [EU project Dynamics–Aerosol–Chemistry–Cloud Interactions in West Africa (DACCIWA)]

*Note several key simplifications in Fig. 1 for the purposes of visual effect: atmospheric properties are separated, but often occur simultaneously and throughout the atmosphere; the depicted path of sunlight is one option, but sunlight emerging at the top of the atmosphere will come from many different paths; sunlight reflected by the surface will need to travel back through the same gases (and likely other properties) on its way back to the top of the atmosphere, which is not shown. The spectra in Fig. 1 are generated with SBDART using a set of arbitrary but realistic atmospheric and surface properties.

** SCIAMACHY = Scanning Imaging Absorption Spectrometer for Atmospheric Chartography.

Jake completed his PhD at Reading in 2018 and now works at the NOAA Earth System Research Laboratory (ESRL) in Boulder, Colorado.

Extending the predictability of flood hazard at the global scale

Email: rebecca.emerton@reading.ac.uk

When I started my PhD, there were no global scale operational seasonal forecasts of river flow or flood hazard. Global overviews of upcoming flood events are key for organisations working at the global scale, from water resources management to humanitarian aid, and for regions where no other local or national forecasts are available. While GloFAS (the Global Flood Awareness System, run by the European Centre for Medium-Range Weather Forecasts (ECMWF) and the European Commission Joint Research Centre (JRC) as part of the Copernicus Emergency Management Services) was producing operational, openly-available flood forecasts out to 30 days ahead, there was a need for more extended-range forecast information. Often, due to a lack of hydrological forecasts, seasonal rainfall forecasts are used as a proxy for flood hazard – however, the link between precipitation and floodiness is nonlinear, and recent research has shown that seasonal rainfall forecasts are not necessarily the best indicator of potential flood hazard. The aim of my PhD research was to look into ways in which we could provide earlier warning information, several weeks to months ahead, using hydrological analysis in addition to the meteorology.

Presidente Kuczynski recorre zonas afectadas por lluvias e inund
Flooding in Trujillo, Peru, March 2017 (Photo: Presidencia Perú on Twitter)

Broadly speaking, there are two key ways in which to provide early warning information on seasonal timescales: (1) through statistical analysis based on large-scale climate variability and teleconnections, and (2) by producing dynamical seasonal forecasts using coupled ocean-atmosphere GCMs. Over the past 4.5 years, I worked on providing hydrologically-relevant seasonal forecast products using these two approaches, at the global scale. This blog post will give a quick overview of the two new forecast products we produced as part of this research!

Can we use El Niño to predict flood hazard?

ENSO (the El Niño Southern Oscillation), is known to influence river flow and flooding across much of the globe, and often, statistical historical probabilities of extreme precipitation during El Niño and La Niña (the extremes of ENSO climate variability) are used to provide information on likely flood impacts. Due to its global influence on weather and climate, we decided to assess whether it is possible to use ENSO as a predictor of flood hazard at the global scale, by assessing the links between ENSO and river flow globally, and estimating the equivalent historical probabilities for high and low river flow, to those that are already used for meteorological variables.

With a lack of sufficient river flow observations across much of the globe, we needed to use a reanalysis dataset – but global reanalysis datasets for river flow are few and far between, and none extended beyond ~40 years (which includes a sample of ≤10 El Niños and ≤13 La Niñas). We ended up producing a 20th Century global river flow reconstruction, by forcing the Camaflood hydrological model with ECMWF’s ERA-20CM atmospheric reconstruction, to produce a 10-member river flow dataset covering 1901-2010, which we called ERA-20CM-R.

elnino_flood_hazard_gif_beccalize

Using this dataset, we calculated the percentage of past El Niño and La Niña events, during which the monthly mean river flow exceeded a high flow threshold (the 75th percentile of the 110-year climatology) or fell below a low flow threshold (the 25th percentile), for each month of an El Niño / La Niña. This percentage is then taken as the probability that high or low flow will be observed in future El Niño/La Niña events. Maps of these probabilities are shown above, for El Niño, and all maps for both El Niño and La Niña can be found here. When comparing to the same historical probabilities calculated for precipitation, it is evident that additional information can be gained from considering the hydrology. For example, the River Nile in northern Africa is likely to see low river flow, even though the surrounding area is likely to see more precipitation – because it is influenced more by changes in precipitation upstream. In places that are likely to see more precipitation but in the form of snow, there would be no influence on river flow or flood hazard during the time when more precipitation is expected. However, several months later, there may be no additional precipitation expected, but there may be increased flood hazard due to the melting of more snow than normal – so we’re able to see a lagged influence of ENSO on river flow in some regions.

While there are locations where these probabilities are high and can provide a useful forecast of hydrological extremes, across much of the globe, the probabilities are lower and much more uncertain (see here for more info on uncertainty in these forecasts) than might be useful for decision-making purposes.

Providing openly-available seasonal river flow forecasts, globally

For the next ‘chapter’ of my PhD, we looked into the feasibility of providing seasonal forecasts of river flow at the global scale. Providing global-scale flood forecasts in the medium-range has only become possible in recent years, and extended-range flood forecasting was highlighted as a grand challenge and likely future development in hydro-meteorological forecasting.

To do this, I worked with Ervin Zsoter at ECMWF, to drive the GloFAS hydrological model (Lisflood) with reforecasts from ECMWF’s latest seasonal forecasting system, SEAS5, to produce seasonal forecasts of river flow. We also forced Lisflood with the new ERA5 reanalysis, to produce an ERA5-R river flow reanalysis with which to initialise Lisflood, and to provide a climatology. The system set-up is shown in the flowchart below.

glofas_seasonal_flowchart_POSTER_EGU

I also worked with colleagues at ECMWF to design forecast products for a GloFAS seasonal outlook, based on a combination of features from the GloFAS flood forecasts, and the EFAS (the European Flood Awareness System) seasonal outlook, and incorporating feedback from users of EFAS.

After ~1 year of working on getting the system set up and finalising the forecast products, including a four-month research placement at ECMWF, the first GloFAS -Seasonal forecast was released in November 2017, with the release of SEAS5. GloFAS-Seasonal is now running operationally at ECMWF, providing forecasts of high and low weekly-averaged river flow for the global river network, up to 4 months ahead, with 3 new forecast layers available through the GloFAS interface. These provide a forecast overview for 307 major river basins, a map of the forecast for the entire river network at the sub-basin scale, and ensemble hydrographs at thousands of locations across the globe (which change with each forecast depending on forecast probabilities). New forecasts are produced once per month, and released on the 10th of each month. You can find more information on each of the different forecast layers and the system set-up here, and you can access the (openly available) forecasts here. ERA5-R, ERA-20CM-R and the GloFAS-Seasonal reforecasts are also all freely available – just get in touch! GloFAS-Seasonal will continue to be developed by ECMWF and the JRC, and has already been updated to v2.0, including a calibrated version of the hydrological model.

NEW_WEB_figure1_basins
Screenshot of the GloFAS seasonal outlook at www.globalfloods.eu

So, over the course of my PhD, we developed two new seasonal forecasts for hydrological extremes, at the global scale. You may be wondering whether they’re skilful, or in fact, which one provides the most useful forecasts! For information on the skill or ‘potential usefulness’ of GloFAS-Seasonal, head to our paper, and stay tuned for a paper coming soon (hopefully! [update: this paper has just been accepted and can be accessed online here]) on the ‘most useful approach for forecasting hydrological extremes during El Niño’, in which we compare the skill of the two forecasts at predicting observed high and low flow events during El Niño.

 

With thanks to my PhD supervisors & co-authors:

Hannah Cloke1, Liz Stephens1, Florian Pappenberger2, Steve Woolnough1, Ervin Zsoter2, Peter Salamon3, Louise Arnal1,2, Christel Prudhomme2, Davide Muraro3

1University of Reading, 2ECMWF, 3European Commission Joint Research Centre

Modelling windstorm losses in a climate model

Extratropical cyclones cause vast amounts of damage across Europe throughout the winter seasons. The damage from these cyclones mainly comes from the associated severe winds. The most intense cyclones have gusts of over 200 kilometres per hour, resulting in substantial damage to property and forestry, for example, the Great Storm of 1987 uprooted approximately 15 million trees in one night. The average loss from these storms is over $2 billion per year (Schwierz et al. 2010) and is second only to Atlantic Hurricanes globally in terms of insured losses from natural hazards. However, the most severe cyclones such as Lothar (26/12/1999) and Kyrill (18/1/2007) can cause losses in excess of $10 billion (Munich Re, 2016). One property of extratropical cyclones is that they have a tendency to cluster (to arrive in groups – see example in Figure 1), and in such cases these impacts can be greatly increased. For example Windstorm Lothar was followed just one day later by Windstorm Martin and the two storms combined caused losses of over $15 billion. The large-scale atmospheric dynamics associated with clustering events have been discussed in a previous blog post and also in the scientific literature (Pinto et al., 2014; Priestley et al. 2017).

Picture1
Figure 1. Composite visible satellite image from 11 February 2014 of 4 extratropical cyclones over the North Atlantic (circled) (NASA).

A large part of my PhD has involved investigating exactly how important the clustering of cyclones is on losses across Europe during the winter. In order to do this, I have used 918 years of high resolution coupled climate model data from HiGEM (Shaffrey et al., 2017) which provides a huge amount of winter seasons and cyclone events for analysis.

In order to understand how clustering affects losses, I first of all need to know how much loss/damage is associated with each individual cyclone. This is done using a measure called the Storm Severity Index (SSI – Leckebusch et al., 2008), which is a proxy for losses that is based on the 10-metre wind field of the cyclone events. The SSI is a good proxy for windstorm loss. Firstly, it scales the wind speed in any particular location by the 98th percentile of the wind speed climatology in that location. This scaling ensures that only the most severe winds at any one point are considered, as different locations have different perspectives on what would be classed as ‘damaging’. This exceedance above the 98th percentile is then raised to the power of 3 due to damage from wind being a highly non-linear function. Finally, we apply a population density weighting to our calculations. This weighting is required because a hypothetical gust of 40 m/s across London will cause considerably more damage than the same gust across far northern Scandinavia, and the population density is a good approximation for the density of insured property. An example of the SSI that has been calculated for Windstorm Lothar is shown in Figure 2.

 

figure_2_blog_2018_new
Figure 2. (a) Wind footprint of Windstorm Lothar (25-27/12/1999) – 10 metre wind speed in coloured contours (m/s). Black line is the track of Lothar with points every 6 hours (black dots). (b) The SSI field of Windstorm Lothar. All data from ERA-Interim.

 

From Figure 2b you can see how most of the damage from Windstorm Lothar was concentrated across central/northern France and also across southern Germany. This is because the winds here were most extreme relative to what is the climatology. Even though the winds are highest across the North Atlantic Ocean, the lack of insured property, and a much high climatological winter mean wind speed, means that we do not observe losses/damage from Windstorm Lothar in these locations.

figure_3_blog_2018_new
Figure 3. The average SSI for 918 years of HiGEM data.

 

I can apply the SSI to all of the individual cyclone events in HiGEM and therefore can construct a climatology of where windstorm losses occur. Figure 3 shows the average loss across all 918 years of HiGEM. You can see that the losses are concentrated in a band from southern UK towards Poland in an easterly direction. This mainly covers the countries of Great Britain, Belgium, The Netherlands, France, Germany, and Denmark.

This blog post introduces my methodology of calculating and investigating the losses associated with the winter season extratropical cyclones. Work in Priestley et al. (2018) uses this methodology to investigate the role of clustering on winter windstorm losses.

This work has been funded by the SCENARIO NERC DTP and also co-sponsored by Aon Benfield.

 

Email: m.d.k.priestley@pgr.reading.ac.uk

 

References

Leckebusch, G. C., Renggli, D., and Ulbrich, U. 2008. Development and application of an objective storm severity measure for the Northeast Atlantic region. Meteorologische Zeitschrift. https://doi.org/10.1127/0941-2948/2008/0323.

Munich Re. 2016. Loss events in Europe 1980 – 2015. 10 costliest winter storms ordered by overall losses. https://www.munichre.com/touch/naturalhazards/en/natcatservice/significant-natural-catastrophes/index.html

Pinto, J. G., Gómara, I., Masato, G., Dacre, H. F., Woollings, T., and Caballero, R. 2014. Large-scale dynamics associated with clustering of extratropical cyclones affecting Western Europe. Journal of Geophysical Research: Atmospheres. https://doi.org/10.1002/2014JD022305.

Priestley, M. D. K., Dacre, H. F., Shaffrey, L. C., Hodges, K. I., and Pinto, J. G. 2018. The role of European windstorm clustering for extreme seasonal losses as determined from a high resolution climate model, Nat. Hazards Earth Syst. Sci. Discuss., https://doi.org/10.5194/nhess-2018-165, in review.

Priestley, M. D. K., Pinto, J. G., Dacre, H. F., and Shaffrey, L. C. 2017. Rossby wave breaking, the upper level jet, and serial clustering of extratropical cyclones in western Europe. Geophysical Research Letters. https://doi.org/10.1002/2016GL071277.

Schwierz, C., Köllner-Heck, P., Zenklusen Mutter, E. et al. 2010. Modelling European winter wind storm losses in current and future climate. Climatic Change. https://doi.org/10.1007/s10584-009-9712-1.

Shaffrey, L. C., Hodson, D., Robson, J., Stevens, D., Hawkins, E., Polo, I., Stevens, I., Sutton, R. T., Lister, G., Iwi, A., et al. 2017. Decadal predictions with the HiGEM high resolution global coupled climate model: description and basic evaluation, Climate Dynamics, https://doi.org/10.1007/s00382-016-3075-x.

The Role of the Cloud Radiative Effect in the Sensitivity of the Intertropical Convergence Zone to Convective Mixing

Email: j.f.talib@pgr.reading.ac.uk

Talib, J., S.J. Woolnough, N.P. Klingaman, and C.E. Holloway, 2018: The Role of the Cloud Radiative Effect in the Sensitivity of the Intertropical Convergence Zone to Convective Mixing. J. Climate, 31, 6821–6838, https://doi.org/10.1175/JCLI-D-17-0794.1

Rainfall in the tropics is commonly associated with the Intertropical Convergence Zone (ITCZ), a discontinuous line of convergence collocated at the ascending branch of the Hadley circulation, where strong moist convection leads to high rainfall. What controls the location and intensity of the ITCZ remains a fundamental question in climate science.

ensemble_precip_neat_thesis
Figure 1: Annual-mean, zonal-mean tropical precipitation (mm day-1) from Global Precipitation Climatology Project (GPCP, observations, solid black line) and CMIP5 (current coupled models) output. Dashed line indicates CMIP5 ensemble mean.

In current and previous generations of climate models, the ITCZ is too intense in the Southern Hemisphere, resulting in two annual-mean, zonal-mean tropical precipitation maxima, one in each hemisphere (Figure 1).  Even if we take the same atmospheric models and couple them to a world with only an ocean surface (aquaplanets) with prescribed sea surface temperatues (SSTs), different models simulate different ITCZs (Blackburn et al., 2013).

Within a climate model parameterisations are used to replace processes that are too small-scale or complex to be physically represented in the model. Parameterisation schemes are used to simulate a variety of processes including processes within the boundary layer, radiative fluxes and atmospheric chemistry. However my work, along with a plethora of others, shows that the representation of the ITCZ is sensitive to the convective parameterisation scheme (Figure 2a). The convective parameterisation scheme simulates the life cycle of clouds within a model grid-box.

Our method of showing that the simulated ITCZ is sensitive to the convective parameterisation scheme is by altering the convective mixing rate in prescribed-SST aquaplanet simulations. The convective mixing rate determines the amount of mixing a convective parcel has with the environmental air, therefore the greater the convective mixing rate, the quicker a convective parcel will become similar to the environmental air, given fixed convective parcel properties.

AEIprecipCREon
Figure 2: Zonal-mean, time-mean (a) precipitation rates (mm day-1}$) and (b) AEI (W m-2) in simulations where the convective mixing rate is varied.

In our study, the structure of the simulated ITCZ is sensitive to the convective mixing rate. Low convective mixing rates simulate a double ITCZ (two precipitation maxima, orange and red lines in Figure 2a), and high convective mixing rates simulate a single ITCZ (blue and black lines).

We then associate these ITCZ structures to the atmospheric energy input (AEI). The AEI is the amount of energy left in the atmosphere once considering the top of the atmosphere and surface energy budgets. We conclude, similar to Bischoff and Schneider, 2016, that when the AEI is positive (negative) at the equator, a single (double) ITCZ is simulated (Figure 2b). When the AEI is negative at the equator, energy is needed to be transported towards the equator for equilibrium. From a mean circulation perspective, this take place in a double ITCZ scenario (Figure 3). A positive AEI at the equator, is associated with poleward energy transport and a single ITCZ.

blog_figure_ITCZ_simulation
Figure 3: Schematic of a single (left) and double ITCZ (right). Blue arrows denote energy transport. In a single ITCZ scenario more energy is transported in the upper branches of the Hadley circulation, resulting in a net-poleward energy transport. In a double ITCZ scenario, more energy is transport equatorward than poleward at low latitudes, leading to an equatorward energy transport.

In our paper, we use this association between the AEI and ITCZ to hypothesize that without the cloud radiative effect (CRE), atmospheric heating due to cloud-radiation interactions, a double ITCZ will be simulated. We also hypothesize that prescribing the CRE will reduce the sensitivity of the ITCZ to convective mixing, as simulated AEI changes are predominately due to CRE changes.

In the rest of the paper we perform simulations with the CRE removed and prescribed to explore further the role of the CRE in the sensitivity of the ITCZ. We conclude that when removing the CRE a double ITCZ becomes more favourable and in both sets of simulations the ITCZ is less sensitive to convective mixing. The remaining sensitivity is associated with latent heat flux alterations.

My future work following this publication explores the role of coupling in the sensitivity of the ITCZ to the convective parameterisation scheme. Prescribing the SSTs implies an arbitary ocean heat transport, however in the real world the ocean heat transport is sensitive to the atmospheric circulation. Does this sensitivity between the ocean heat transport and atmospheric circulation affect the sensitivity of the ITCZ to convective mixing?

Thanks to my funders, SCENARIO NERC DTP, and supervisors for their support for this project.

References:

Blackburn, M. et al., (2013). The Aqua-planet Experiment (APE): Control SST simulation. J. Meteo. Soc. Japan. Ser. II, 91, 17–56.

Bischoff, T. and Schneider, T. (2016). The Equatorial Energy Balance, ITCZ Position, and Double-ITCZ Bifurcations. J. Climate., 29(8), 2997–3013, and Corrigendum, 29(19), 7167–7167.

 

Baroclinic and Barotropic Annular Modes of Variability

Email: l.boljka@pgr.reading.ac.uk

Modes of variability are climatological features that have global effects on regional climate and weather. They are identified through spatial structures and the timeseries associated with them (so-called EOF/PC analysis, which finds the largest variability of a given atmospheric field). Examples of modes of variability include El Niño Southern Oscillation, Madden-Julian Oscillation, North Atlantic Oscillation, Annular modes, etc. The latter are named after the “annulus” (a region bounded by two concentric circles) as they occur in the Earth’s midlatitudes (a band of atmosphere bounded by the polar and tropical regions, Fig. 1), and are the most important modes of midlatitude variability, generally representing 20-30% of the variability in a field.

Southern_Hemi_Antarctica
Figure 1: Southern Hemisphere midlatitudes (red concentric circles) as annulus, region where annular modes have the largest impacts. Source.

We know two types of annular modes: baroclinic (based on eddy kinetic energy, a proxy for eddy activity and an indicator of storm-track intensity) and barotropic (based on zonal mean zonal wind, representing the north-south shifts of the jet stream) (Fig. 2). The latter are usually referred to as Southern (SAM or Antarctic Oscillation) or Northern (NAM or Arctic Oscillation) Annular Mode (depending on the hemisphere), have generally quasi-barotropic (uniform) vertical structure, and impact the temperature variations, sea-ice distribution, and storm paths in both hemispheres with timescales of about 10 days. The former are referred to as BAM (baroclinic annular mode) and exhibit strong vertical structure associated with strong vertical wind shear (baroclinicity), and their impacts are yet to be determined (e.g. Thompson and Barnes 2014, Marshall et al. 2017). These two modes of variability are linked to the key processes of the midlatitude tropospheric dynamics that are involved in the growth (baroclinic processes) and decay (barotropic processes) of midlatitude storms. The growth stage of the midlatitude storms is conventionally associated with increase in eddy kinetic energy (EKE) and the decay stage with decrease in EKE.

ThompsonWoodworth_Fig2a_SAM_2f_BAM(1)
Figure 2: Barotropic annular mode (right), based on zonal wind (contours), associated with eddy momentum flux (shading); Baroclinic annular mode (left), based on eddy kinetic energy (contours), associated with eddy heat flux (shading). Source: Thompson and Woodworth (2014).

However, recent observational studies (e.g. Thompson and Woodworth 2014) have suggested decoupling of baroclinic and barotropic components of atmospheric variability in the Southern Hemisphere (i.e. no correlation between the BAM and SAM) and a simpler formulation of the EKE budget that only depends on eddy heat fluxes and BAM (Thompson et al. 2017). Using cross-spectrum analysis, we empirically test the validity of the suggested relationship between EKE and heat flux at different timescales (Boljka et al. 2018). Two different relationships are identified in Fig. 3: 1) a regime where EKE and eddy heat flux relationship holds well (periods longer than 10 days; intermediate timescale); and 2) a regime where this relationship breaks down (periods shorter than 10 days; synoptic timescale). For the relationship to hold (by construction), the imaginary part of the cross-spectrum must follow the angular frequency line and the real part must be constant. This is only true at the intermediate timescales. Hence, the suggested decoupling of baroclinic and barotropic components found in Thompson and Woodworth (2014) only works at intermediate timescales. This is consistent with our theoretical model (Boljka and Shepherd 2018), which predicts decoupling under synoptic temporal and spatial averaging. At synoptic timescales, processes such as barotropic momentum fluxes (closely related to the latitudinal shifts in the jet stream) contribute to the variability in EKE. This is consistent with the dynamics of storms that occur on timescales shorter than 10 days (e.g. Simmons and Hoskins 1978). This is further discussed in Boljka et al. (2018).

EKE_hflux_cross_spectrum_blog
Figure 3: Imaginary (black solid line) and Real (grey solid line) parts of cross-spectrum between EKE and eddy heat flux. Black dashed line shows the angular frequency (if the tested relationship holds, the imaginary part of cross-spectrum follows this line), the red line distinguishes between the two frequency regimes discussed in text. Source: Boljka et al. (2018).

References

Boljka, L., and T. G. Shepherd, 2018: A multiscale asymptotic theory of extratropical wave, mean-flow interaction. J. Atmos. Sci., in press.

Boljka, L., T. G. Shepherd, and M. Blackburn, 2018: On the coupling between barotropic and baroclinic modes of extratropical atmospheric variability. J. Atmos. Sci., in review.

Marshall, G. J., D. W. J. Thompson, and M. R. van den Broeke, 2017: The signature of Southern Hemisphere atmospheric circulation patterns in Antarctic precipitation. Geophys. Res. Lett., 44, 11,580–11,589.

Simmons, A. J., and B. J. Hoskins, 1978: The life cycles of some nonlinear baroclinic waves. J. Atmos. Sci., 35, 414–432.

Thompson, D. W. J., and E. A. Barnes, 2014: Periodic variability in the large-scale Southern Hemisphere atmospheric circulation. Science, 343, 641–645.

Thompson, D. W. J., B. R. Crow, and E. A. Barnes, 2017: Intraseasonal periodicity in the Southern Hemisphere circulation on regional spatial scales. J. Atmos. Sci., 74, 865–877.

Thompson, D. W. J., and J. D. Woodworth, 2014: Barotropic and baroclinic annular variability in the Southern Hemisphere. J. Atmos. Sci., 71, 1480–1493.

Climate model systematic biases in the Maritime Continent

Email: y.y.toh@pgr.reading.ac.uk

The Maritime Continent commonly refers to the groups of islands of Indonesia, Borneo, New Guinea and the surrounding seas in the literature. My study area covers the Maritime Continent domain from 20°S to 20°N and 80°E to 160°E as shown in Figure 1. This includes Indonesia, Malaysia, Brunei, Singapore, Philippines, Papua New Guinea, Solomon islands, northern Australia and parts of mainland Southeast Asia including Thailand, Laos, Cambodia, Vietnam and Myanmar.

subsetF1
Figure 1: JJA precipitation (mm/day) and 850 hPa wind (m s−1) for (a) GPCP and ERA-interim, (b) MMM biases and (c)–(j) AMIP biases for 1979–2008 over the Maritime Continent region (20°S–20ºN, 80°E–160ºE). Third panel shows the Maritime Continent domain and land-sea mask

The ability of climate model to simulate the mean climate and climate variability over the Maritime Continent remains a modelling challenge (Jourdain et al. 2013). Our study examines the fidelity of Coupled Model Intercomparison Project phase 5 (CMIP5) models at simulating mean climate over the Maritime Continent. We find that there is a considerable spread in the performance of the Atmospheric Model Intercomparison Project (AMIP) models in reproducing the seasonal mean climate and annual cycle over the Maritime Continent region. The multi-model mean (MMM) (Figure 1b) JJA precipitation and 850hPa wind biases with respect to observations (Figure 1a) are small compared to individual model biases (Figure 1c-j) over the Maritime Continent. Figure 1 shows only a subset of Fig. 2 from Toh et al. (2017), for the full figure and paper please click here.

We also investigate the model characteristics that may be potential sources of bias. We find that AMIP model performance is largely unrelated to model horizontal resolution. Instead, a model’s local Maritime Continent biases are somewhat related to its biases in the local Hadley circulation and global monsoon.

cluster2
Figure 2: Latitude-time plot of precipitation zonally averaged between 80°E and 160°E for (a) GPCP, (b) Cluster I and (c) Cluster II. White dashed line shows the position of the maximum precipitation each month. Precipitation biases with respect to GPCP for (d) Cluster I and (e) Cluster II.

To characterize model systematic biases in the AMIP runs and determine if these biases are related to common factors elsewhere in the tropics, we performed cluster analysis on Maritime Continent annual cycle precipitation. Our analysis resulted in two distinct clusters. Cluster I (Figure 2b,d) is able to reproduce the observed seasonal migration of Maritime Continent precipitation, but it overestimates the precipitation, especially during the JJA and SON seasons. Cluster II (Figure 2c,e) simulate weaker seasonal migration of Intertropical Convergence Zone (ITCZ) than observed, and the maximum rainfall position stays closer to the equator throughout the year. Tropics-wide properties of clusters also demonstrate a connection between errors at regional scale of the Maritime Continent and errors at large scale circulation and global monsoon.

On the other hand, comparison with coupled models showed that air-sea coupling yielded complex impacts on Maritime Continent precipitation biases. One of the outstanding problems in the coupled CMIP5 models is the sea surface temperature (SST) biases in tropical ocean basins. Our study highlighted central Pacific and western Indian Oceans as the key regions which exhibit the most surface temperature correlation with Maritime Continent mean state precipitation in the coupled CMIP5 models. Future work will investigate the impact of SST perturbations in these two regions on Maritime Continent precipitation using Atmospheric General Circulation Model (AGCM) sensitivity experiments.

 

 

References:

Jourdain N.C., Gupta A.S., Taschetto A.S., Ummenhofer C.C., Moise A.F., Ashok K. (2013) The Indo-Australian monsoon and its relationship to ENSO and IOD in reanalysis data and the CMIP3/CMIP5 simulations. Climate Dynamics. 41(11–12):3073–3102

Toh, Y.Y., Turner, A.G., Johnson, S.J., & Holloway, C.E. (2017). Maritime Continent seasonal climate biases in AMIP experiments of the CMIP5 multimodel ensemble. Climate Dynamics. doi: 10.1007/s00382-017-3641-x

Should we be ‘Leaf’-ing out vegetation when parameterising the aerodynamic properties of urban areas?

Email: C.W.Kent@pgr.reading.ac.uk

When modelling urban areas, vegetation is often ignored in attempt to simplify an already complex problem. However, vegetation is present in all urban environments and it is not going anywhere… For reasons ranging from sustainability to improvements in human well-being, green spaces are increasingly becoming part of urban planning agendas. Incorporating vegetation is therefore a key part of modelling urban climates. Vegetation provides numerous (dis)services in the urban environment, each of which requires individual attention (Salmond et al. 2016). However, one of my research interests is how vegetation influences the aerodynamic properties of urban areas.

Two aerodynamic parameters can be used to represent the aerodynamic properties of a surface: the zero-plane displacement (zd) and aerodynamic roughness length (z0). The zero-plane displacement is the vertical displacement of the wind-speed profile due to the presence of surface roughness elements. The aerodynamic roughness length is a length scale which describes the magnitude of surface roughness. Together they help define the shape and form of the wind-speed profile which is expected above a surface (Fig. 1).

blogpostpic

Figure 1: Representation of the wind-speed profile above a group of roughness elements. The black dots represent an idealised logarithmic wind-speed profile which is determined using the zero-plane displacement (zd) and aerodynamic roughness length (z0) (lines) of the surface.

For an urban site, zd and z0 may be determined using three categories of methods: reference-based, morphometric and anemometric. Reference-based methods require a comparison of the site to previously published pictures or look up tables (e.g. Grimmond and Oke 1999); morphometric methods describe zd and z0 as a function of roughness-element geometry; and, anemometric methods use in-situ observations. The aerodynamic parameters of a site may vary considerably depending upon which of these methods are used, but efforts are being made to understand which parameters are most appropriate to use for accurate wind-speed estimations (Kent et al. 2017a).

Within the morphometric category (i.e. using roughness-element geometry) sophisticated methods have been developed for buildings or vegetation only. However, until recently no method existed to describe the effects of both buildings and vegetation in combination. A recent development overcomes this, whereby the heights of all roughness elements are considered alongside a porosity correction for vegetation (Kent et al. 2017b). Specifically, the porosity correction is applied to the space occupied and drag exerted by vegetation.

The development is assessed across several areas typical of a European city, ranging from a densely-built city centre to an urban park. The results demonstrate that where buildings are the dominant roughness elements (i.e. taller and occupying more space), vegetation does not obviously influence the calculated geometry of the surface, nor the aerodynamic parameters and the estimated wind speed. However, as vegetation begins to occupy a greater amount of space and becomes as tall as (or larger) than buildings, the influence of vegetation is obvious. Expectedly, the implications are greatest in an urban park, where overlooking vegetation means that wind speeds may be slowed by up to a factor of three.

Up to now, experiments such as those in the wind tunnel focus upon buildings or trees in isolation. Certainly, future experiments which consider both buildings and vegetation will be valuable to continue to understand the interaction within and between these roughness elements, in addition to assessing the parameterisation.

References

Grimmond CSB, Oke TR (1999) Aerodynamic properties of urban areas derived from analysis of surface form. J Appl Meteorol and Clim 38:1262-1292.

Kent CW, Grimmond CSB, Barlow J, Gatey D, Kotthaus S, Lindberg F, Halios CH (2017a) Evaluation of Urban Local-Scale Aerodynamic Parameters: Implications for the Vertical Profile of Wind Speed and for Source Areas. Boundary-Layer Meteorology 164: 183-213.

Kent CW, Grimmond CSB, Gatey D (2017b) Aerodynamic roughness parameters in cities: Inclusion of vegetation. Journal of Wind Engineering and Industrial Aerodynamics 169: 168-176.

Salmond JA, Tadaki M, Vardoulakis S, Arbuthnott K, Coutts A, Demuzere M, Dirks KN, Heaviside C, Lim S, Macintyre H (2016) Health and climate related ecosystem services provided by street trees in the urban environment. Environ Health 15:95.

The onset and end of wet seasons over Africa

Email: c.m.dunning@pgr.reading.ac.uk

For many Africans, the timing of the wet season is of crucial importance, especially for those reliant upon subsistence agriculture, who depend on the seasonal rains for crop irrigation. In addition, the wet season recharges lakes, rivers and water storage tanks which constitute the domestic water supply in some areas. The timing of the wet season also affects the availability of energy from hydroelectric schemes, and has impacts upon the prevalence of certain disease carrying vectors, such as mosquitoes.

Climate change is already threatening many vulnerable populations, and changes in the timing or intensity of the wet season, or increasing uncertainty in the timing of the onset, may lead to significant socio-economic impacts. But before we consider future projections or past changes in the seasonality, we need to go back a few steps.

The first step is to find a method for determining when the wet season starts and ends (its ‘onset’ and ‘cessation’). In order to look at large-scale shifts in the timing of the wet season and relate this to wider-scale drivers, this method needs to be applicable across the entirety of continental Africa. Most previous methods for determining the onset focus on the national to regional scale, and are dependent on the exceedance of a certain threshold e.g. the first week with at least 20mm of rainfall, with one rainfall event of more than 10mm, and no dry spell of more than 10 days after the rain event for the next month. While such definitions work well at a national scale they are not applicable at a continental scale where rainfall amounts vary substantially. A threshold suitable for the dry countries at the fringes of the Sahara would not be suitable in the wetter East African highlands.

In addition to a vast range of rainfall amounts, the African continent also spans multiple climatic regimes. The seasonal cycle of precipitation over continental Africa is largely driven by the seasonal progression of the ITCZ and associated rain belts, which follows the maximum incoming solar radiation. In the boreal summer, when the thermal equator sits between the equator and the Tropic of Cancer, the ITCZ sits north of the equator and West Africa and the Sahel experience a wet season. During the boreal autumn the ITCZ moves south, and southern Africa experiences a wet season during the austral summer, followed by the northward return of the ITCZ during the boreal spring. As a consequence of this, central African regions and the Horn of Africa experience two wet seasons per year – one as the ITCZ travels north, and a second as the ITCZ travels south. A method for determining the onset and cessation at the continental scale thus needs to account for regions with multiple wet seasons per year.

In our paper (available here) we propose such a method, based on the method of Liebmann et al (2012). The method has three steps:

  • Firstly, determine the number of seasons experienced per year at the location (or grid point) of interest. This is achieved using harmonic analysis – the amplitude of the first and second harmonic were computed, using the entire timeseries and their ratio compared. If the ratio was greater than 1.0, i.e. the amplitude of the second harmonic was greater than the amplitude of the first harmonic then the grid point was defined as having two wet seasons per year (biannual), if the ratio was less than one then it was defined as having an annual regime. Figure 1 shows the ratio for one African rainfall dataset (TARCATv2). Three regions are identified as biannual regions; the Horn of Africa, an equatorial strip extending from Gabon to Uganda and a small region on the southern West African coastline.

    blog_fig1
    Figure 1: Location of regions with one and two seasons per year, determined using harmonic analysis. Yellow indicates two seasons per year, while pink/purple indicates one season per year. Computed from TARCATv2 data.
  • Secondly the period of the year when the wet season occurs was determined. This was achieved by looking for minima and maxima in the climatological cumulative daily rainfall anomaly to identify one or two seasons.
  • The third and final stage is to calculate the onset and cessation dates for each year. This is done by looking for the minima and maxima in the cumulative daily rainfall anomaly, calculated for each season.

Figure 2 shows the seasonal progression of the onset and cessation, with the patterns observed in agreement with those expected from the driving physical mechanisms, and continuous progression across the annual/biannual boundaries. Over West Africa and the Sahel, Figure 2a-b shows zonally-contiguous progression patterns with onset following the onset of the long rains and moving north, and cessation moving southward, preceding the end of the short rains. Over southern Africa Figure 2c-d shows the onset over southern Africa starting in the north-west and south-east, following the onset of the short rains, reaching the East African coast last, and cessation starting at the Zimbabwe, Mozambique, South Africa border and spreading out radially into the cessation of the long rains.

As well as testing the method for compatibility with known physical drivers of African rainfall, agreement across multiple satellite-based rainfall estimates was also examined. In general, good agreement was found across the datasets, particularly for regions with an annual regime and over the biannual region of East Africa.

blog_fig2
Figure 2: Southward and northward progression of the onset and cessation across the annual/biannual boundaries, computed using GPCP daily rainfall data 1998-2013.

The advantage of having a method that works at the continental scale is the ability to look at the impact of large-scale oscillations on wider-scale variability. One application of this method was to investigate the impact of El Niño upon both the annual rains and short rains (Figure 3). In Figure 3 we see the well-documented dipole in rainfall anomaly, with higher rainfall totals over 0–15°S and the Horn of Africa in El Niño years and the opposite between 15°S and 30°S.  This anomaly is stronger when we use this method compared with using standard meteorological seasons. We can also see that while the lower rainfall to the south is colocated with later onset dates and a consequentially shorter season, the higher rainfall over the Horn of Africa is associated with later cessation of the short rains, with only small differences in onset date.

blog_fig3
Figure 3: a-c) Composite of onset, cessation and wet season rainfall in El Niño years for annual rains and short rains, minus the mean over 1982-2013, computed using CHIRPS data d) Oct-Feb rainfall anomaly in  years (CHIRPS).

In addition to using this method for research purposes, its application within an operational setting is also being explored. Hopefully, the method will be included within the Rainwatch platform, which will be able to provide users with a probabilistic estimate of whether or not the season has started, based on the rainfall experienced so far that year, and historical rainfall data.

For more details, please see the paper detailing this work:

Dunning, C.M., E Black, and R.P. Allan (2016) The onset and cessation of seasonal rainfall over Africa, Journal of Geophysical Research: Atmospheres, 121 11,405-11,424, doi: 10.1002/2016JD025428

References:

Liebmann, B., I. Bladé, G. N. Kiladis, L. M. Carvalho, G. B. Senay, D. Allured, S. Leroux, and C. Funk (2012), Seasonality of African precipitation from 1996 to 2009, J. Clim.25(12), 4304–4322.

Can we really use El Niño to predict flooding?

R. Emerton, H. Cloke, E. Stephens, E. Zsoter, S. Woolnough, F. Pappenberger (2017). Complex picture for likelihood of ENSO-driven flood hazard. Nature Communications. doi: 10.1038/NCOMMS14796

Email: r.e.emerton@pgr.reading.ac.uk

When an El Niño is declared, or even forecast, we think back to memorable past El Niños (such as 1997/98), and begin to ask whether we will see the same impacts. Will California receive a lot of rainfall? Will we see droughts in tropical Asia and Australia? Will Peru experience the same devastating floods as in 1997/98, and 1982/83?

news_headlines_elnino

El Niño and La Niña, which see changes in the ocean temperatures in the tropical Pacific, are well known to affect weather, and indeed river flow and flooding, around the globe. But how well can we estimate the potential impacts of El Niño and La Niña, and how likely flooding is to occur?

This question is what some of us in the Water@Reading research group at the University of Reading have been looking to answer in our recent publication in Nature Communications. As part of our multi- and inter-disciplinary research, we work closely with the Red Cross / Red Crescent Climate Centre (RCCC), who are working on an initiative called Forecast-based Financing (FbF, Coughlan de Perez et al.). FbF aims to distribute aid (for example providing water purification tablets to prevent spread of disease, or digging trenches to divert flood water) ahead of a flood, based on forecasts. This approach helps to reduce the impact of the flood in the first place, rather than working to undo the damage once the flood has already occurred.

Photo credit: Red Cross / Red Crescent Climate Centre

In Peru, previous strong El Niños in 1982/83 and 1997/98 had resulted in devastating floods in several regions. As such, when forecasts in early 2015 began to indicate a very strong El Niño was developing, the RCCC and forecasters at the Peruvian national hydrological and meteorology agency (SENAMHI) began to look into the likelihood of flooding, and what FbF actions might need to be taken.

Typically, statistical products indicating the historical probability (likelihood [%] based on what happened during past El Niños) of extreme precipitation are used as a proxy for whether a region will experience flooding during an El Niño (or La Niña), such as these maps produced by the IRI (International Research Institute for Climate and Society). You may also have seen maps which circle regions of the globe that will be drier / warmer / wetter / cooler – we’ll come back to these shortly.

These rainfall maps show that Peru, alongside several other regions of the world, is likely to see more rainfall than usual during an El Niño. But does this necessarily mean there will be floods? And what products are out there indicating the effect of El Niño on rivers across the globe?

For organisations working at the global scale, such as the RCCC and other humanitarian aid agencies, global overviews of potential impacts are key in taking decisions on where to focus resources during an El Niño or La Niña. While these maps are useful for looking at the likely changes in precipitation, it has been shown that the link between precipitation and flood magnitude is nonlinear (Stephens et al.),  – more rain does not necessarily equal floods – so how does this transfer to the potential for flooding?

The motivation behind this work was to provide similar information, but taking into account the hydrology as well as the meteorology. We wanted to answer the question “what is the probability of flooding during El Niño?” not only for Peru, but for the global river network.

To do this, we have taken the new ECMWF ERA-20CM ensemble model reconstruction of the atmosphere, and run this through a hydrological model to produce the first 20th century global hydrological reconstruction of river flow. Using this new dataset, we have for the first time estimated the historical probability of increased or decreased flood hazard (defined as abnormally high or low river flow) during an El Niño (or La Niña), for the global river network.

elnino_flood_hazard_gif_beccalize
Figure 1: The probability of increased (blue) or decreased (red) flood hazard during each month of an El Nino. Based on the ensemble mean of the ERA-20CM-R 20th century river flow reconstruction.

The question – “what is the probability of flooding during El Niño?”, however, remains difficult to answer. We now have maps of the probability of abnormally high or low river flow (see Figure 1), and we see clear differences between the hydrological analysis and precipitation. It is also evident that the probabilities themselves are often lower, and much more uncertain, than might be useful – how do you make a decision on whether to provide aid to an area worried about flooding, when the probability of that flooding is 50%?

blob_floodhazard_comparison_map
Figure 2: Historical probability of increased / decreased flood hazard map for February, with overlay showing the typical impact map for winter during an El Nino. This highlights the complexity of the link between El Nino and flooding compared to the information usually available.

The likely impacts are much more complex than is often perceived and reported – going back to the afore-mentioned maps that circle regions of the globe and what their impact will be (warmer, drier, wetter?) – these maps portray these impacts as a certainty, not a probability, with the same impacts occurring across huge areas. For example, in Figure 2, we take one of the maps from our results, which indicates the probability of increased or decreased flood hazard in one month during an El Niño, and draw over this these oft-seen circles of potential impacts. In doing this, we remove all information on how likely (or unlikely) the impacts are, smaller scale changes within these circles (in some cases our flood hazard map even indicates a different impact), and a lot of the potential impacts outside of these circles – not to mention the likely impacts can change dramatically from one month to the next. For those organisations that take actions based on such information, it is important to be aware of the uncertainties surrounding the likely impacts of El Niño and La Niña.

“We conclude that while it may seem possible to use historical probabilities to evaluate regions across the globe that are more likely to be at risk of flooding during an El Niño / La Niña, and indeed circle large areas of the globe under one banner of wetter or drier, the reality is much more complex.”

PS. During the winter of 2015/16, our results estimated an ~80% likelihood of increased flood hazard in northern coastal Peru, with only ~10% uncertainty surrounding this. The RCCC took FbF actions to protect thousands of families from potentially devastating floods driven by one of the strongest El Niños on records. While flooding did occur, this was not as severe as expected based on the strength of the El Niño. More recently, during the past few months (January – March 2017), anomalously high sea surface temperatures (SSTs) in the far eastern Pacific (known as a “coastal El Niño” in Peru but not widely acknowledged as an El Niño because central Pacific SSTs are not anomalously warm) have led to devastating flooding in several regions and significant loss of life. And Peru wasn’t the only place that didn’t see the impacts it expected in 2015/16; other regions of the world, such as the US, also saw more rainfall than normal in places that were expected to be drier, and California didn’t receive the deluge they were perhaps hoping for. It’s important to remember that no two El Niños are the same, and El Niño will not be the only influence on the weather around the globe. While El Niño and La Niña can provide some added predictability to the atmosphere, the impacts are far from certain.

Presidente Kuczynski recorre zonas afectadas por lluvias e inund
Flooded areas of Trujillo, Peru, March 2017. Photo credit: Presidencia Peru, via Floodlist

Full reference:

R. Emerton, H. Cloke, E. Stephens, E. Zsoter, S. Woolnough, F. Pappenberger (2017). Complex picture for likelihood of ENSO-driven flood hazard. Nature Communications. doi: 10.1038/NCOMMS14796

Press Release: