Mesoscale

Air-sea heat fluxes at the oceanic mesoscale: the impact of the ratio of ocean-to-atmosphere grid resolution

sophiaashbyLeave a comment

Sophia Moreton – s.moreton@pgr.reading.ac.uk

Sea surface temperature (SST) anomalies are vital for regulating the earth’s weather and climate.  The generation and reduction of these SST anomalies are largely determined by air-sea heat fluxes, particularly turbulent heat fluxes (latent and sensible).

The turbulent heat flux feedback (THFF) is a critical parameter, which measures the change in the net air-sea turbulent heat flux in response to a 1 K change in SST. So far in current research, this feedback is well known at large scales, i.e. over the whole ocean basin. However, a quantification of this feedback at much smaller spatial scales (10-100km) over individual mesoscale ocean eddies remains absent.

Why do we care about air-sea feedbacks at the oceanic mesoscale?

Both heat and momentum air-sea exchanges at the mesoscale impact the local and large-scale atmosphere (e.g. shifting storm tracks) and alter the strength of western boundary currents and the large-scale ocean gyre circulation. However, research into this field to date is hindered by the lack of high spatial resolution in observational data at the air-sea interface.

Therefore our study uses three high-resolution configurations from the UK Met Office coupled climate model (HadGEM3-GC3). We provide the first global estimate of turbulent heat flux feedback (α) over individually tracked and composite-averaged coherent mesoscale eddies, which ranges between 35 to 45 Wm-2K-1 depending on eddy amplitude.

Estimates of the turbulent heat flux feedback (THFF) are split, depending if the feedback is calculated using SST on the ocean grid (α0) or after regridding SST to the atmosphere (αA). An example of αA using regridded SST anomalies (SSTA) is given in Fig.1 for large-amplitude eddies in the highest ocean-atmosphere resolution available (a 25km atmosphere coupled to a 1/12° ocean, labelled ‘N512-12’).

Figure 1: A scatter plot of the relationship (THFF, αA) between regridded SST (SSTA) and THF anomalies. αA is the gradient of the linear regression line (black) +/- the 95% confidence interval (shown by the text). The data is from eddy snapshots averaged over 1 year, denoted by ‘< >’. Only large-amplitude eddies in the N512-12 configuration (25km atmosphere – 1/12° ocean) are plotted.

Why is the feedback so sensitive to the ratio of grid resolution?

In high-resolution coupled climate models, the atmospheric resolution is typically coarser than in its ocean component although, to date, a quantification of what the ocean-atmosphere ratio of grid resolution should be remains absent.

We prove increasing the ratio of atmosphere-to-ocean grid resolution in coupled climate models can lead to a large underestimation of turbulent heat flux feedback over mesoscale eddies, by as much as 75% for a 6:1 resolution ratio, as circled in Fig. 2 from a 60km atmosphere coupled to a 1/12° ocean. An underestimation of the feedback is consistent across all eddy amplitudes (A) and all three model configurations shown (Fig. 2); it suggests SST anomalies within these eddies are likely to be not reduced enough by air-sea fluxes of heat, and consequently will remain too large.

The underestimation stems from the calculation of the air-sea heat fluxes in the HadGEM3-GC3.1 model on the coarser atmospheric grid, instead of the finer ocean grid. Many other climate models do the same. At present, for the long spin-ups needed for climate simulations, it is unrealistic to expect the atmospheric resolution to match the very fine (10km) ocean resolution in coupled climate models, i.e. to create a one-to-one grid ratio. Therefore, to minimise this underestimation in the feedback at mesoscales, we advise air-sea heat fluxes should be computed on the finer oceanic grid.

Figure 2: Estimates of the turbulent heat flux feedback (THFF) across different eddy amplitudes (A) for α0 (lighter colours) and αA (darker colours, using regridded SST) for three model configurations: N512-12, N216-12 and N216-025. The ocean and atmosphere resolutions are added in red for each. Increasing the ratio of grid resolution, underestimates the THFF (as α0 differs from αA). The horizontal bars indicate the width of the eddy amplitude bins, and the vertical error bars indicate 95% confidence intervals.

Correctly simulating the air-sea heat flux feedback over mesoscale eddies is fundamental to realistically represent their interaction with the local and large-scale atmosphere and feedback on the ocean, to improve our predictions of the earth’s climate.

For a full analysis of the results, including a decomposition of the turbulent heat flux feedback, the reader is referred to Moreton et al., 2021, Air-Sea Turbulent Heat Flux Feedback over Mesoscale Eddies, GRL (in review).

Manuscript available: https://doi.org/10.1002/essoar.10505981.1

This work lays the foundation for my current work, evaluating how mesoscale air-sea heat fluxes feedback and alter the strength of large-scale ocean gyre circulation, using the MIT general circulation model (MITgcm).

This work is funded by a NERC CASE studentship with the Met Office, UK.

Quantifying the skill of convection-permitting ensemble forecasts for the sea-breeze occurrence

Carlo CafaroLeave a comment

Email: carlo.cafaro@pgr.reading.ac.uk

On the afternoon of 16th August 2004, the village of Boscastle on the north coast of Cornwall was severely damaged by flooding (Golding et al., 2005). This is one example of high impact hazardous weather associated with small meso- and convective-scale weather phenomena, the prediction of which can be uncertain even a few hours ahead (Lorenz, 1969; Hohenegger and Schar, 2007). Taking advantage of the increased computer power (e.g. https://www.metoffice.gov.uk/research/technology/supercomputer) this has motivated many operational and research forecasting centres to introduce convection-permitting ensemble prediction systems (CP-EPSs), in order to give timely weather warnings of severe weather.

However, despite being an exciting new forecasting technology, CP-EPSs place a heavy burden on the computational resources of forecasting centres. They are usually run on limited areas with initial and boundary conditions provided by global lower resolution ensembles (LR-EPS). They also produce large amounts of data which needs to be rapidly digested and utilized by operational forecasters. Assessing whether the convective-scale ensemble is likely to provide useful additional information is key to successful real-time utilisation of this data. Similarly, knowing where equivalent information can be gained (even if partially) from LR-EPS using statistical/dynamical post-processing both extends lead time (due to faster production time) and also potentially provides information in regions where no convective-scale ensemble is available.

There have been many studies on the verification of CP-EPSs (Klasa et al., 2018, Hagelin et al., 2017, Barret et al., 2016, Beck et al., 2016 amongst the others), but none of them has dealt with the quantification of the skill gained by CP-EPSs in comparison with LR-EPSs, when fully exploited, for specific weather phenomena and for a long enough evaluation period.

In my PhD, I have focused on the sea-breeze phenomenon for different reasons:

  1. Sea breezes have an impact on air quality by advecting pollutants, on heat stress by providing a relief on hot days and also on convection by providing a trigger, especially when interacting with other mesoscale flows (see for examples figure 1 or figures 6, 7 in Golding et al., 2005).
  2. Sea breezes occur on small spatio-temporal scales which are properly resolved at convection-permitting resolutions, but their occurrence is still influenced by synoptic-scale conditions, which are resolved by the global LR-EPS.

Blog_Figure1
Figure 1: MODIS visible of the southeast of Italy on 6th June 2018, 1020 UTC. This shows thunderstorms occurring in the middle of the peninsula, probably triggered by sea-breezes.
Source: worldview.earthdata.nasa.gov

Therefore this study aims to investigate whether the sea breeze is predictable by only knowing a few predictors or whether the better representation of fine-scale structures (e.g. orography, topography) by the CP-EPS implies a better sea-breeze prediction.

In order to estimate probabilistic forecasts from both the models, two different methods have been applied. A novel tracking algorithm for the identification of sea-breeze front, in the domain represented in figure 2, was applied to CP-EPSs data. A Bayesian model was used instead to estimate the probability of sea-breeze conditioned on two LR-EPSs predictors and trained on CP-EPSs data. More details can be found in Cafaro et al. (2018).

Cafaro_Fig2
Figure 2: A map showing the orography over the south UK domain. Orography data are from MOGREPS-UK. The solid box encloses the sub-domain used in this study with red dots indicating the location of synoptic weather stations. Source: Cafaro et al. (2018)

The results of the probabilistic verification are shown in figure 3. Reliability (REL) and resolution (RES) terms have been computed decomposing the Brier score (BS) and Information gain (IGN) score. Finally, scores differences (BSD and IG) have been computed to quantify any gain in the skill by the CP-EPS. Figure 3 shows that CP-EPS forecast is significantly more skilful than the Bayesian forecast. Nevertheless, the Bayesian forecast has more resolution than a climatological forecast (figure 3e,f), which has no resolution by construction.

Cafaro_Fig11
Figure 3: (a)-(d) Reliability and resolution terms calculated for both the forecasts (green for the CP-EPS forecast and blue for LR-EPSs). (e) and (f) represent the Brier score difference (BSD) and Information gain (IG) respectively. Error bars represent the 95th confidence interval. Positive values of BSD and IG indicate that CP-EPS forecast is more skilful. Source: Cafaro et al. (2018)

This study shows the additional skill provided by the Met Office convection-permitting ensemble forecast for the sea-breeze prediction. The ability of CP-EPSs to resolve meso-scale dynamical features is thus proven to be important and only two large-scale predictors, relevant for the sea-breeze, are not sufficient for skilful prediction.

It is believed that both the methodologies can, in principle, be applied to other locations of the world and it is thus hoped they could be used operationally.

References:

Barrett, A. I., Gray, S. L., Kirshbaum, D. J., Roberts, N. M., Schultz, D. M., and Fairman J. G. (2016). The utility of convection-permitting ensembles for the prediction of stationary convective bands. Monthly Weather Review, 144(3):1093–1114, doi: 10.1175/MWR-D-15-0148.1

Beck,  J., Bouttier, F., Wiegand, L., Gebhardt, C., Eagle, C., and Roberts, N. (2016). Development and verification of two convection-allowing multi-model ensembles over Western europe. Quarterly Journal of the Royal Meteorological Society, 142(700):2808–2826, doi: 10.1002/qj.2870

Cafaro C., Frame T. H. A., Methven J., Roberts N. and Broecker J. (2018), The added value of convection-permitting ensemble forecasts of sea breeze compared to a Bayesian forecast driven by the global ensemble, Quarterly Journal of the Royal Meteorological Society., under review.

Golding, B. , Clark, P. and May, B. (2005), The Boscastle flood: Meteorological analysis of the conditions leading to flooding on 16 August 2004. Weather, 60: 230-235, doi: 10.1256/wea.71.05

Hagelin, S., Son, J., Swinbank, R., McCabe, A., Roberts, N., and Tennant, W. (2017). The Met Office convective-scale ensemble, MOGREPS-UK. Quarterly Journal of the Royal Meteorological Society, 143(708):2846–2861, doi: 10.1002/qj.3135

Hohenegger, C. and Schar, C. (2007). Atmospheric predictability at synoptic versus cloud-resolving scales. Bulletin of the American Meteorological Society, 88(11):1783–1794, doi: 10.1175/BAMS-88-11-1783

Klasa, C., Arpagaus, M., Walser, A., and Wernli, H. (2018). An evaluation of the convection-permitting ensemble cosmo-e for three contrasting precipitation events in Switzerland. Quarterly Journal of the Royal Meteorological Society, 144(712):744–764, doi: 10.1002/qj.3245

Lorenz, E. N. (1969). Predictability of a flow which possesses many scales of motion.Tellus, 21:289 – 307, doi: 10.1111/j.2153-3490.1969.tb00444.x

Modelling windstorm losses in a climate model

Matthew PriestleyLeave a comment

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.

Future of Cumulus Parametrization conference, Delft, July 10-14, 2017

Mark MuetzelfeldtLeave a comment

Email: m.muetzelfeldt@pgr.reading.ac.uk

For a small city, Delft punches above its weight. It is famous for many things, including its celebrated Delftware (Figure 1). It was also the birthplace of one of the Dutch masters, Johannes Vermeer, who coincidentally painted some fine cityscapes with cumulus clouds in them (Figure 2). There is a university of technology with some impressive architecture (Figure 3). It holds the dubious honour of being the location of the first assassination using a pistol (or so we were told by our tour guide), when William of Orange was shot in 1584. To this list, it can now add hosting a one-week conference on the future of cumulus parametrization, and hopefully bringing about more of these conferences in the future.

Delftware_display

Figure 1: Delftware.

Vermeer-view-of-delft

Figure 2: Delft with canopy of cumulus clouds. By Johannes Vermeer, 1661.

Delft_AULA

Figure 3: AULA conference centre at Delft University of Technology – where we were based for the duration of the conference.

So what is a cumulus parametrization scheme? The key idea is as follows. Numerical weather and climate models work by splitting the atmosphere into a grid, with a corresponding grid length representing the length of each of the grid cells. By solving equations that govern how the wind, pressure and heating interact, models can then be used to predict what the weather will be like days in advance in the case of weather modelling. Or a model can predict how the climate will react to any forcings over longer timescales. However, any phenomena that are substantially smaller than this grid scale will not be “seen” by the models. For example, a large cumulonimbus cloud may have a horizontal extent of around 2km, whereas individual grid cells could be 50km in the case of a climate model. A cumulonimbus cloud will therefore not be explicitly modelled, but it will still have an effect on the grid cell in which it is located – in terms of how much heating and moistening it produces at different levels. To capture this effect, the clouds are parametrized, that is, the vertical profile of the heating and moistening due to the clouds are calculated based on the conditions in the grid cell, and this then affects the grid-scale values of these variables. A similar idea applies for shallow cumulus clouds, such as the cumulus humilis in Vermeer’s painting (Figure 2), or present-day Delft (Figure 3).

These cumulus parametrization schemes are a large source of uncertainty in current weather and climate models. The conference was aimed at bringing together the community of modellers working on these schemes, and working out which might be the best directions to go in to improve these schemes, and consequently weather and climate models.

Each day was a mixture of listening to presentations, looking at posters and breakout discussion groups in the afternoon, as well as plenty of time for coffee and meeting new people. The presentations covered a lot of ground: from presenting work on state-of-the-art parametrization schemes, to looking at how the schemes perform in operational models, to focusing on one small aspect of a scheme and modelling how that behaves in a high resolution model (50m resolution) that can explicitly model individual clouds. The posters were a great chance to see the in-depth work that had been done, and to talk to and exchange ideas with other scientists.

Certain ideas for improving the parametrization schemes resurfaced repeatedly. The need for scale-awareness, where the response of the parametrization scheme takes into account the model resolution, was discussed. One idea for doing this was the use of stochastic schemes to represent the uncertainty of the number of clouds in a given grid cell. The concept of memory also cropped up – where the scheme remembers if it had been active at a given grid cell in a previous point in time. This also ties into the idea of transitions between cloud regimes, e.g. when a stratocumulus layer splits up into individual cumulus clouds. Many other, sometimes esoteric, concepts were discussed, such as the role of cold pools, how much tuning of climate models is desirable and acceptable, how we should test our schemes, and what the process of developing the schemes should look like.

In the breakout groups, everyone was encouraged to contribute, which made for an inclusive atmosphere in which all points of view were taken on board. Some of the key points of agreement from these were that it was a good idea to have these conferences, and we should do it more often! Hopefully, in two years’ time, another PhD student will write a post on how the next meeting has gone. We also agreed that it would be beneficial to be able to share data from our different high resolution runs, as well as to be able to compare code for the different schemes.

The conference provided a picture of what the current thinking on cumulus parametrization is, as well as which directions people think are promising for the future. It also provided a means for the community to come together and discuss ideas for how to improve these schemes, and how to collaborate more closely with future projects such as ParaCon and HD(CP)2.

Sting Jet: the poisonous (and windy) tail of some of the most intense UK storms

Ambrogio Volonté2 Comments

Email: a.volonte@pgr.reading.ac.uk

IDL TIFF file
Figure 1: Windstorm Tini (12 Feb 2014) passes over the British Isles bringing extreme winds. A Sting Jet has been identified in the storm. Image courtesy of NASA Earth Observatory

It was the morning of 16th October when South East England got battered by the Great Storm of 1987. Extreme winds occurred, with gusts of 70 knots or more recorded continually for three or four consecutive hours and maximum gusts up to 100 knots. The damage was huge across the country with 15 million trees blown down and 18 fatalities.

case_study_great_storm_fig011
Figure 2: Surface wind gusts in the Great Storm of 1987. Image courtesy of UK Met Office.

The forecast issued on the evening of 15th October failed to identify the incoming hazard but forecasters were not to blame as the strongest winds were actually due to a phenomenon that had yet to be discovered at the time: the Sting Jet. A new topic of weather-related research had started: what was the cause of the exceptionally strong winds in the Great Storm?

It was in Reading at the beginning of 21st century that scientists came up with the first formal description of those winds, using observations and model simulations. Following the intuitions of Norwegian forecasters they used the term Sting Jet, the ‘sting at the end of the tail’. Using some imagination we can see the resemblance of the bent-back cloud head with a scorpion’s tail: strong winds coming out from its tip and descending towards the surface can then be seen as the poisonous sting at the end of the tail.

Conceptual+model+of+storm+development
Figure 3: Conceptual model of a sting-jet extratropical cyclone, from Clark et al, 2005. As the cloud head bends back and the cold front moves ahead we can see the Sting Jet exiting from the cloud tip and descending into the opening frontal fracture.  WJ: Warm conveyor belt. CJ: Cold conveyor belt. SJ: Sting jet.

In the last decade sting-jet research progressed steadily with observational, modelling and climatological studies confirming that the strong winds can occur relatively often, that they form in intense extratropical cyclones with a particular shape and are caused by an additional airstream that is neither related to the Cold nor to the Warm Conveyor Belt. The key questions are currently focused on the dynamics of Sting Jets: how do they form and accelerate?

Works recently published (and others about to come out, stay tuned!) claim that although the Sting Jet occurs in an area in which fairly strong winds would already be expected given the morphology of the storm, a further mechanism of acceleration is needed to take into account its full strength. In fact, it is the onset of mesoscale instabilities and the occurrence of evaporative cooling on the airstream that enhances its descent and acceleration, generating a focused intense jet (see references for more details). It is thus necessary a synergy between the general dynamics of the storm and the local processes in the cloud head in order to produce what we call the Sting Jet .

plot_3D_sj ccb_short
Figure 4: Sting Jet (green) and Cold Conveyor Belt (blue) in the simulations of Windstorm Tini. The animation shows how the onset of the strongest winds is related to the descent of the Sting Jet. For further details on this animation and on the analysis of Windstorm Tini see here.

References:

http://www.metoffice.gov.uk/learning/learn-about-the-weather/weather-phenomena/case-studies/great-storm

Browning, K. A. (2004), The sting at the end of the tail: Damaging winds associated with extratropical cyclones. Q.J.R. Meteorol. Soc., 130: 375–399. doi:10.1256/qj.02.143

Clark, P. A., K. A. Browning, and C. Wang (2005), The sting at the end of the tail: Model diagnostics of fine-scale three-dimensional structure of the cloud head. Q.J.R. Meteorol. Soc., 131: 2263–2292. doi:10.1256/qj.04.36

Martínez-Alvarado, O., L.H. Baker, S.L. Gray, J. Methven, and R.S. Plant (2014), Distinguishing the Cold Conveyor Belt and Sting Jet Airstreams in an Intense Extratropical Cyclone. Mon. Wea. Rev., 142, 2571–2595, doi: 10.1175/MWR-D-13-00348.1.

Hart, N.G., S.L. Gray, and P.A. Clark, 0: Sting-jet windstorms over the North Atlantic: Climatology and contribution to extreme wind risk. J. Climate, 0, doi: 10.1175/JCLI-D-16-0791.1.

Volonté, A., P.A. Clark, S.L. Gray. The role of Mesoscale Instabilities in the Sting-Jet dynamics in Windstorm Tini. Poster presented at European Geosciences Union – General Assembly 2017, Dynamical Meteorology (General session)

Understanding the dynamics of cyclone clustering

Matthew Priestley2 Comments

Priestley, M. D. K., J. G. Pinto, H. F. Dacre, and L. C. Shaffrey (2016), Rossby wave breaking, the upper level jet, and serial clustering of extratropical cyclones in western Europe, Geophys. Res. Lett., 43, doi:10.1002/2016GL071277.

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

Extratropical cyclones are the number one natural hazard that affects western Europe (Della-Marta, 2010). These cyclones can cause widespread socio-economic damage through extreme wind gusts that can damage property, and also through intense precipitation, which may result in prolonged flood events. For example the intensely stormy winter of 2013/2014 saw 456mm of rain fall in under 90 days across the UK; this broke records nationwide as 175% of the seasonal average fell (Kendon & McCarthy, 2015). One particular storm in this season was cyclone Tini (figure 1), this was a very deep cyclone (minimum pressure – 952 hPa) which brought peak gusts of over 100 mph to the UK. These gusts caused widespread structural damage that resulted in 20,000 homes losing power. These extremes can be considerably worse when multiple extratropical cyclones affect one specific geographical region in a very short space of time. This is known as cyclone clustering. Some of the most damaging clustering events can result in huge insured losses, for example the storms in the winter of 1999/2000 resulted in €16 billion of losses (Swiss Re, 2016); this being more than 10 times the annual average.

figure-1
Figure 1. A Meteosat visible satellite image at 12 UTC on February 12th 2014 showing cyclone Tini over the UK. Image credit to NEODAAS/University of Dundee.

Up until recently cyclone clustering had been given little attention in terms of scientific research, despite it being a widely accepted phenomenon in the scientific community. With these events being such high risk events it is important to understand the atmospheric dynamics that are associated with these events; and this is exactly what we have been doing recently. In our new study we attempt to characterise cyclone clustering in several different locations and associate each different set of clusters with a different dynamical setup in the upper troposphere. The different locations we focus on are defined by three areas, one encompassing the UK and centred at 55°N. Our other two areas are 10° to the north and south of this (centred at 65°N and 45°N.) The previous study of Pinto et al. (2014) examined several winter seasons and found links between the upper-level jet, Rossby wave breaking (RWB) and the occurrence of clustering. RWB is the meridional overturning of air in the upper troposphere. It is identified using the potential temperature (θ) field on the dynamical tropopause, with a reversal of the normal equator-pole θ gradient representing RWB. This identification method is explained in full in Masato et al. (2013) and also illustrated in figure 2. We have greatly expanded on this analysis to look at all winter clustering events from 1979/1980 to 2014/2015 and their connection with these dynamical features.

schematic_box1
Figure 2. Evolution of Rossby waves on the tropopause. RWB occurs when these waves overturn by a significant amount. H: High potential temperature; L: Low potential temperature (Priestley et al., 2017).

We find that when we get clustering it is accompanied with a much stronger jet at 250 hPa than in the climatology, with average speeds peaking at over 50 ms-1 (figures 3a-c). In all cases there is also a much greater presence of RWB in regions not seen from the climatology (Figure 3d). In figure 3a there is more RWB to the south of the jet, in figure 3b there is an increased presence on both the northern and southern flanks, and finally in figure 3c there is much more RWB to the north. The presence of this anomalous RWB transfers momentum into the jet, which acts to strengthen and extend it toward western Europe.

figure-2
Figure 3. The dynamical setup for clustering occurring at (a) 65°N; (b) 55°N; and (c) 45°N. The climatology is shown in (d). Coloured shading is the average potential temperature on the tropopause, black contours are the average 250 hPa wind speeds and black crosses are where RWB is occurring.

The location of the RWB controls the jet tilt; more RWB to the south of the jet acts to angle it more northwards (figure 3a), there is a southward deflection when there is more RWB to the north of the jet (figure 3c). The presence of RWB on both sides extends it along a more central axis (figure 3b). Therefore the occurrence of RWB in a particular location and the resultant angle of the jet acts to direct cyclones to various parts of western Europe in quick succession.

In our recently published study we go into much more detail regarding the variability associated with these dynamics and also how the jet and RWB interact in time. This can be found at http://dx.doi.org/10.1002/2016GL071277.

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

References

Della-Marta, P. M., Liniger, M. A., Appenzeller, C., Bresch, D. N., Köllner-Heck, P., & Muccione, V. (2010). Improved estimates of the European winter windstorm climate and the risk of reinsurance loss using climate model data. Journal of Applied Meteorolo

Kendon, M., & McCarthy, M. (2015). The UK’s wet and stormy winter of 2013/2014. Weather, 70(2), 40-47.

Masato, G., Hoskins, B. J., & Woollings, T. (2013). Wave-breaking characteristics of Northern Hemisphere winter blocking: A two-dimensional approach. Journal of Climate, 26(13), 4535-4549.

Pinto, J. G., Gómara, I., Masato, G., Dacre, H. F., Woollings, T., & Caballero, R. (2014). Large‐scale dynamics associated with clustering of extratropical cyclones affecting Western Europe. Journal of Geophysical Research: Atmospheres, 119(24).

Priestley, M. D. K., J. G. Pinto, H. F. Dacre, and L. C. Shaffrey (2017). The role of cyclone clustering during the stormy winter of 2013/2014. Manuscript in preparation.

Swiss Re. (2016). Winter storm clusters in Europe, Swiss Re publishing, Zurich, 16 pp., http://www.swissre.com/library/winter_storm_clusters_in_europe.html. Accessed 24/11/16.

From foehn to intense rainfall: the importance of Alps in influencing the regional weather

Ambrogio VolontéLeave a comment

Email: a.volonte@pgr.reading.ac.uk

dsc_0196
Figure 1: View from Monte Lema (Italy-Switzerland) looking West. The Lake Maggiore region and the southern Alpine foothills are visible in the foreground whereas Monte Rosa and the Pennine Alps behind them are partially hidden by a characteristic foehn wall.  (A. Volonté, 4 January 2017)

The interaction between atmospheric flow and topography is at the origin of various important weather phenomena, as we have already seen in Carly Wright’s blog post. When a mountain range is particularly high and extended it can even block or deflect weather systems, as it happens with the Alps. For example, in Figure 1 we can see the main Alpine range with its over-4000m-high peaks blocking a cold front coming from the north. The main ridge acts as a wall, enhancing condensation and precipitation processes on the upstream side (stau condition) and leaving clear skies on the downstream lee side, where dry and mild katabatic foehn winds flow. The contrast is striking between sunny weather on Lake Maggiore and snowy conditions over Monte Rosa, just a few miles apart. The same phenomenon is shown in Figure 2 with a satellite image that highlights how a cold front coming from northwest gets blocked by the Alpine barrier. A person enjoying the sunny day in the southern side of the Alps, if unaware of this mechanism, would be very surprised  to know that the current weather is so different on the other side of the range.

poplex-2014295-terra-1km
Figure 2: Satellite image (MODIS-NASA) over the Alps and Po Valley on 22 October 2014

poplex-2016348-terra-1km
Figure 3: same as Figure 1 but on 13 December 2016

A comparison with Figure 3 helps to notice that in Figure 2 the shape of the cloud band closely mirrors the mountain range. As an additional remark,  this comparison shows that foehn bring clear skies even in the Po Valley, having blown away the typical mist/fog occurring in the region in Autumn and Winter months in high pressure regimes. The  stau/foehn dynamics is actually very fascinating, and you can read more about it in Elvidge and Renfrew (2015 ) and in Miltenberger et al. (2016), among others. Unfortunately, the interaction of weather systems with the Alps can often trigger very damaging phenomena, like heavy and long-lasting precipitation on one side of the slope, and this is what the rest of this post will be focused on. In fact, the most recent event of this kind just happened at the end of November, with intense and long-lasting rain affecting the southern slope of the Alps  and causing floods particularly in the Piedmont region, in northwestern Italy ( Figure 4).

tanaro2
Figure 4: River Tanaro flooding in the town of Garessio, 24 November 2016 (Piedmont, Italy). Source: http://www.corrierenazionale.it

arpa_piemonte
Figure 5: rainfall accumulated between 21 and 26 November 2016 in the Piedmont region. Source: Regional Agency for the protection of the Environment – Piedmont

Figure 5 shows that the accumulated rainfall in the event goes over 300 mm in a large band that follows the shape of the southern Alpine slope in the region (see map of Piedmont, from Google Maps), reaching even 600 mm in a few places. This situation is the result of moist southerly flow being blocked by the Alps and thus causing ascent and consequent precipitation to persist on the same areas for up to five days. It is quite common to see quasi-stationary troughs enter the Mediterranean region during Autumn months causing strong and long-lasting moist flows to move towards the Alps. Hence, it is crucial to understand  where the heaviest precipitation will occur. In other words, will it rain the most on top of the ridge or on the upstream plain? What processes are controlling the location of heavy precipitation with respect to the slope?

The study published by Davolio et al. (2016), available here and originated from my master degree’s thesis, tackles this issue focusing on northeastern Italy. In fact, the analysis includes three case studies in which heavy and long-lasting rain affected the eastern Alps and other three case studies in which intense rainfall was mainly located on the upstream plain. Although all the events showed common large-scale patterns and similar mesoscale settings, characterised by moist southerly low-level flow interacting with the Alps, the rainfall distribution turned out to be very dissimilar. The study highlights that the two precipitation regimes strongly differ in terms of interaction of the flow with the mountain barrier. When the flow is able to go over the Alps the heaviest rain occurs on top of the ridge. When the flow is instead blocked and deflected by the ridge (flow around), creating a so-called barrier wind, intense convection is triggered on the upstream plain (Figure 6) .

qj2731
Figure 6: Schematic diagram of the key mechanisms governing the two different wind and precipitation patterns over NE Italy. (a) Blocked low-level flow, barrier wind, convergence and deep convection over the plain, upstream the orography. (b) Flow over conditions with orographic lifting and precipitation mainly over the Alps. From Davolio et al. (2016)

convection
Figure 7: cross section going from the Adriatic Sea to the Alps in one of the events simulated. Equivalent potential temperature is shaded, thick black lines indicate clouds while arrows show tangent wind component. See Davolio et al. (2016)

The key mechanism that explains this different evolution is connected to the thermodynamic state of the impinging flow. In fact, when the southerly moist and warm air gets close to the Alpine barrier it is lifted above the colder air already present at the base of the orography. It can be said that the colder air behaves as a first effective mountain for the incoming flow. If this lifting process triggers convection, then the persistence of a blocked-flow condition is highly favoured (see Figure 7). On the contrary, if this initial lifting process does not trigger convection the intense moist flow will eventually be able to go over the ridge, where a more substantial ascent will take place, causing heavy rain on the ridge top. This study also looks at numerical parameters used in more idealised analyses (like in Miglietta and Rotunno (2009)), finding a good agreement with the theory.

To summarise, we can say that the Alpine range is able to significantly modify weather systems when interacting with them. Thus, an in-depth understanding of the processes taking place during the interaction, along with a coherent model is necessary to capture correctly the effects on the local weather, being either a rainfall enhancement, the occurrence of foehn winds or various other phenomena.

References

Davolio, S., Volonté A., Manzato A., Pucillo A., Cicogna A. and Ferrario M.E. (2016), Mechanisms producing different precipitation patterns over north-eastern Italy: insights from HyMeX-SOP1 and previous events. Q.J.R. Meteorol. Soc., 142 (Suppl 1): 188-205. doi:10.1002/qj.2731

Elvidge A. D., Renfrew, I. A. (2015). The causes of foehn warming in the lee of mountains. Bull. Am. Meteorol. Soc. 97: 455–466, doi:10.1175/BAMS-D-14-00194.1.

Miglietta M. and Rotunno R., (2009) Numerical Simulations of Conditionally Unstable Flows over a Mountain Ridge. J. Atmos. Sci., 66, 1865–1885, doi: 10.1175/2009JAS2902.1. 

Miltenberger, A. K., Reynolds, S. and Sprenger, M. (2016), Revisiting the latent heating contribution to foehn warming: Lagrangian analysis of two foehn events over the Swiss Alps. Q.J.R. Meteorol. Soc., 142: 2194–2204. doi:10.1002/qj.2816

Stationary Orographic Rainbands

Carly Wright1 Comment

Email: c.j.wright@pgr.reading.ac.uk

Small-scale rainbands often form downwind of mountainous terrain. Although relatively small in scale (a few tens of km across by up to ~100 km in length), these often poorly forecast bands can cause localised flooding as they can be associated with intense precipitation over several hours due to the anchoring effect of orography (Barrett et al., 2013).   Figure 1 shows a flash flood caused by a rainband situated over Cockermouth in 2009.  In some regions of southern France orographic banded convection can contribute 40% of the total rainfall (Cosma et al., 2002).  Rainbands occur in various locations and under different synoptic regimes and environmental conditions making them difficult to examine their properties and determine their occurrence in a systematic way (Kirshbaum et al. 2007a,b, Fairman et al. 2016).  My PhD considers the ability of current operational forecast models to represent these bands and the environmental controls on their formation.

blogfig1
Figure 1: Flash flood event caused by a rainband situated over Cockermouth, Cumbria, UK in 2009

 

What is a rainband?

  • A cloud and precipitation structure associated with an area of rainfall which is significantly elongated
  • Stationary (situated over the same location) with continuous triggering
  • Can form in response to moist, unstable air following over complex terrain
  • Narrow in width ~2-10 km with varying length scales from 10 – 100’s km

 

blogfig2
Figure 2: Schematic showing the difference between cellular and banded convection

To examine the ability of current operational forecast models to represent these bands a case study was chosen which was first introduced by Barrett, et al. (2016).  The radar observations during the event showed a clear band along The Great Glen Fault, Scotland (Figure 3).  However, Barrett, et al. (2016) concluded that neither the operational forecast or the operational ensemble forecast captured the nature of the rainband.  For more information on ensemble models see one of our previous blog posts by David Flack Showers: How well can we predict them?.

blogfig3
Figure 3: Radar observations of precipitation accumulation over a six hour period (between 3-9 am) showing a rainband located over The Great Glen Fault, Scotland on 29 December 2012.

Localised convergence and increased convective available potential energy along the fault supported the formation of the rainband.  To determine the effect of model resolution on the model’s representation of the rainband, a forecast was performed with the horizontal gird spacing decreased to 500 m from 1.5 km.  In this forecast a rainband formed in the correct location which generated precipitation accumulations close to those observed, but with a time displacement.  The robustness of this forecast skill improvement is being assessed by performing an ensemble of these convection-permitting simulations.  Results suggest that accurate representation of these mesoscale rainbands requires resolutions higher than those used operationally by national weather centres.

Idealised numerical simulations have been used to investigate the environmental conditions leading to the formation of these rainbands.  The theoretical dependence of the partitioning of dry flow over and around mountains on the non-dimensional mountain height is well understood.  For this project I examine the effect of this dependence on rainband formation in a moist environment.  Preliminary analysis of the results show that the characteristics of rainbands are controlled by more than just the non-dimensional mountain height, even though this parameter is known to be sufficient to determine flow behaviour relative to mountains.

This work has been funded by the Natural Environmental Research Council (NERC) under the project PREcipitation STructures over Orography (PRESTO), for more project information click here.

References

Barrett, A. I., S. L. Gray, D. J. Kirshbaum, N. M. Roberts, D. M. Schultz, and J. G. Fairman, 2015: Synoptic Versus Orographic Control on Stationary Convective Banding. Quart. J. Roy. Meteorol. Soc., 141, 1101–1113, doi:10.1002/qj.2409.

— 2016: The Utility of Convection-Permitting Ensembles for the Prediction of Stationary Convective Bands. Mon. Wea. Rev., 144, 10931114, doi:10.1175/MWR-D-15-0148.1.

Cosma, S., E. Richard, and F. Minsicloux, 2002: The Role of Small-Scale Orographic Features in the Spatial Distribution of Precipitation. Quart. J. Roy. Meteorol. Soc., 128, 75–92, doi:10.1256/00359000260498798.

Fairman, J. G., D. M. Schultz, D. J. Kirshbaum, S. L. Gray, and A. I. Barrett, 2016: Climatology of Banded Precipitation over the Contiguous United States. Mon. Wea. Rev., 144,4553–4568, doi: 10.1175/MWR-D-16-0015.1.

Kirshbaum, D. J., G. H. Bryan, R. Rotunno, and D. R. Durran, 2007a: The Triggering of Orographic Rainbands by Small-Scale Topography. J. Atmos. Sci., 64, 1530–1549, doi:10.1175/JAS3924.1.

Kirshbaum, D. J., R. Rotunno, and G. H. Bryan, 2007b: The Spacing of Orographic Rainbands Triggered by Small-Scale Topography. J. Atmos. Sci., 64, 4222–4245, doi:10.1175/2007JAS2335.1.