Climate scientists have a lot of insight into the factors driving weather systems in the mid-latitudes, where the rotation of the earth is an important influence. The tropics are less well served, and this can be a problem for global climate models which don’t capture many of the phenomena observed in the tropics that well.
What we do know about the tropics however is that despite significant contrasts in sea surface temperatures (Fig. 1) there is very little horizontal temperature variation in the atmosphere (Fig. 2) – because the Coriolis force (due to the Earth’s rotation) that enables this gradient in more temperate climates is not present. We believe that the large-scale circulation acts to minimise the effect these surface contrasts have higher up. This suggests a model for vertical wind which cools the air over warmer surfaces and warms it where the surface is cool, called the Weak Temperature Gradient (WTG) Approximation, that is frequently used in studying the climate in the tropics.
Thermodynamic ideas have been around for some 200 years. Carnot, a Frenchman worried about Britain’s industrial might underpinning its military potential(!), studied the efficiency of heat engines and showed that the maximum mechanical work generated by an engine is determined by the ratio of the temperatures at which energy enters and leaves the system. It is possible to treat climate systems as heat engines – for example Kerry Emanuel has used Carnot’s idea to estimate the pressure in the eye of a hurricane. I have been building on a recent development of these ideas by Olivier Pauluis at New York University who shows how to divide up the maximum work output of a climate heat engine into the generation of wind, the lifting of moisture and a lost component, which he calls the Gibbs penalty, which is the energetic cost of keeping the atmosphere moist. Typically, 50% of the maximum work output is gobbled up by the Gibbs penalty, 30% is the moisture lifting term and only 20% is used to generate wind.
For my PhD, I have been applying Pauluis’ ideas to a modelled system consisting of two connected tropical regions (one over a cooler surface than the other), which are connected by a circulation given by the weak temperature gradient approximation. I look at how this circulation affects the components of work done by the system. Overall there is no impact – in other words the WTG does not distort the thermodynamics of the underlying system – which is reassuring for those who use it. What is perhaps more interesting however, is that even though the WTG circulation is very weak compared to the winds that we observe in the two columns, it does as much work as is done by the cooler column – in other words its thermodynamic importance is huge. This suggests that further avenues of study may help us better express what drives the climate in the tropics.
Current numerical modelling and data assimilation methods still face problems in strongly nonlinear cases, like in convective scales. A different, but interesting tool to help overcome these issues can be found in the synchronisation theory.
It all started in 1665, when Christiaan Huygens, a Dutch scientist, discovered that his two pendulum clocks were suddenly oscillating in opposite directions, but in a synchronised way. He tried to desynchronise them, by perturbing randomly one of the clocks, but surprisingly, after some time, both devices were synchronised again. He has attributed the phenomenon to the frame both clocks were sharing and after that, synchronisation field was opened to the world.
Figure 1: A drawing by Christiaan Huygens of his experiment in 1665.
Nowadays, researchers use these synchronisation concepts to reach a main goal: synchronise a model (any) with the true evolution of a system, using measurements. And even when only a reduced part of this system is observed, synchronisation between models and the true state can still be achieved. This is quite similar to what data assimilation looks for, as it aims to synchronise a model evolution with the truth by using observations, finding the best estimate of the state evolution and its uncertainty.
So why not investigate the benefits of recent synchronisation findings and combine these concepts with a data assimilation methodology?
At the start of this project, the first noticeable step that should be taken was to open up the synchronisation field to higher-dimension systems, as the experiments performed in the area were all focused on low-dimension, non-realistic systems. To this end, a first new idea was proposed: an ensemble version of a synchronisation scheme, what we are calling EnSynch (Ensemble Synchronisation). Tests with a partly observed 1000-dimension chaotic model show a very efficient correspondence between the model and the true trajectories, both for estimation and prediction periods. Figures 2 and 3 show how our estimates and the truth are on top of each other, i.e. synchronised. Note that we do not have observations for all of the variables in our system. So, it is amazing to obtain the same successful results for the observed and also for the unobserved variables in this system!
Figure 2: Trajectories of 2 variables (top:observed and bottom: unobserved). Blue lines: truth. Green lines: estimates/predictions. (Predictions start after the red lines, i.e. no data assimilation is used.)
Figure 3: Zoom in the trajectory of a variable, showing how the model matches with the truth. Blue line: truth. Red line: our model. Yellow dots: observations.
The second and main idea is to test a combination of this successful EnSynch scheme with a data assimilation method called Particle Filter. As a proper data assimilation methodology, a particle filter provides us the best estimation of the state evolution and its uncertainty. Just to illustrate the importance of data assimilation in following the truth, figure 4 compares the case of only counting on an ensemble of models running freely in a chaotic nonlinear system, with the case of a data assimilation method applied to it.
Figure 4: Trajectories of ensemble members. Blue: with data assimilation. Red: without data assimilation. Truth is in black.
Efficient results are found with the combination between the new EnSynch and the particle filters. An example is shown in figure 5, where particles (ensemble members) of an unobserved variable nicely follow the truth during the assimilation period and also during the forecast stage (after t=100).
Figure 5: Trajectory for an unobserved variable in a 1000-dimension system. Observations occur at every 10 time steps until t=100. Predictions start after t=100.
These results are motivating and the next and big step is to implement this combined system in a bigger atmospheric model. This methodology has been shown to be a promising solution for strongly nonlinear problems and potential benefits are expected for numerical weather prediction in the near future.
Rey, D., M. Eldridge, M. Kostuk, H. Abarbanel, J. Schumann-Bischoff, and U. Parlitz, 2014a: Accurate state and parameter estimation in nonlinear systems with sparse observations. Physics Letters A, 378, 869-873, doi:10.1016/j.physleta.2014.01.027.
Zhu, M., P. J. van Leeuwen, and J. Amezcua, 2016: Implicit equal-weights particle filter. Quart. J. Roy. Meteorol. Soc., 142, 1904-1919, doi:10.1002/qj.2784.
The partitioning of shortwave radiation by vegetation into absorbed, reflected, and transmitted terms is important for most biogeophysical processes including photosynthesis. The most commonly used radiative transfer scheme in climate models does not explicitly account for vegetation architectural effects on shortwave radiation partitioning, and even though detailed 3D radiative transfer schemes have been developed, they are often too computationally expensive and require a large number of parameters.
Using a simple parameterisation, we modified a 1D radiative transfer scheme to simulate the radiative balance consistently with 3D representations. Canopy structure is typically treated via a so called “clumping” factor which acts to reduce the effective leaf area index (LAI) and hence fAPAR (fraction of absorbed photosynthetically radiation, 400-700 nm). Consequently from a production efficiency standpoint it seems intuitive that any consideration of clumping can only lead to reduce GPP (Gross Primary Productivity). We show, to the contrary, that the dominant effect of clumping in more complex models should be to increase photosynthesis on global scales.
The Joint UK Land Environment Simulator (JULES) has recently been modified to include clumping information on a per-plant functional type (PFT) basis (Williams et al., 2017). Here we further modify JULES to read in clumping for each PFT in each grid cell independently. We used a global clumping map derived from MODIS data (He et al., 2012) and ran JULES 4.6 for the year 2008 both with and without clumping using the GL4.0 configuration forced with the WATCH-Forcing-Data-ERA-Interim data set (Weedon et al., 2014). We compare our results against the MTE (Model Tree Ensemble) GPP global data set (Beer et al., 2010).
Fig. 1 shows an almost ubiquitous increase in GPP globally when clumping is included in JULES. In general this improves agreement against the MTE data set (Fig. 2). Spatially the only significant areas where the performance is degraded are some tropical grasslands and savannas (not shown). This is likely due to other model problems, in particular the limited number of PFTs used to represent all vegetation globally. The explanation for the increase in GPP and its spatial pattern is shown in Fig 3. JULES uses a multi-layered canopy scheme coupled to the Farquhar photosynthesis scheme (Farquhar et al., 1980). Changing fAPAR (by including clumping in this case) has largest impacts where GPP is light limited, and this is especially true in tropical forests.
Beer, C. et al. 2010. Terrestrial gross carbon dioxide uptake: global distribution and covariation with climate. Science, 329(5993), pp.834-838.
Farquhar, G.D. et al. 1980. A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species. Planta, 149, 78–90.
He, L. et al. 2012. Global clumping index map derived from the MODIS BRDF product. Remote Sensing of Environment, 119, pp.118-130.
Weedon, G. P. et al. 2014. The WFDEI meteorological forcing data set: WATCH Forcing Data methodology applied to ERA-Interim reanalysis data, Water Resour. Res., 50, 7505–7514.
Williams, K. et al. 2017. Evaluation of JULES-crop performance against site observations of irrigated maize from Mead, Nebraska. Geoscientific Model Development, 10(3), pp.1291-1320.
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.
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.
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 .
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)
Forests play an important role in the global carbon cycle, removing large amounts of CO2 from the atmosphere and thus helping to mitigate the effect of human-induced climate change. The state of the global carbon cycle in the IPCC AR5 suggests that the land surface is the most uncertain component of the global carbon cycle. The response of ecosystem carbon uptake to land use change and disturbance (e.g. fire, felling, insect outbreak) is a large component of this uncertainty. Additionally, there is much disagreement on whether forests and terrestrial ecosystems will continue to remove the same proportion of CO2 from the atmosphere under future climate regimes. It is therefore important to improve our understanding of ecosystem carbon cycle processes in the context of a changing climate.
Here we focus on the effect on ecosystem carbon dynamics of disturbance from selective felling (thinning) at the Alice Holt research forest in Hampshire, UK. Thinning is a management practice used to improve ecosystem services or the quality of a final tree crop and is globally widespread. At Alice Holt a program of thinning was carried out in 2014 where one side of the forest was thinned and the other side left unmanaged. During thinning approximately 46% of trees were removed from the area of interest.
Using the technique of eddy-covariance at flux tower sites we can produce direct measurements of the carbon fluxes in a forest ecosystem. The flux tower at Alice Holt has been producing measurements since 1999 (Wilkinson et al., 2012), a view from the flux tower is shown in Figure 1. These measurements represent the Net Ecosystem Exchange of CO2 (NEE). The NEE is composed of both photosynthesis and respiration fluxes. The total amount of carbon removed from the atmosphere through photosynthesis is termed the Gross Primary Productivity (GPP). The Total Ecosystem Respiration (TER) is made up of autotrophic respiration (Ra) from plants and heterotrophic respiration (Rh) from soil microbes and other organisms incapable of photosynthesis. We then have, NEE = -GPP + TER, so that a negative NEE value represents removal of carbon from the atmosphere and a positive NEE value represents an input of carbon to the atmosphere. A schematic of these fluxes is shown in Figure 2.
The flux tower at Alice Holt is on the boundary between the thinned and unthinned forest. This allows us to partition the NEE observations between the two areas of forest using a flux footprint model (Wilkinson et al., 2016). We also conducted an extensive fieldwork campaign in 2015 to estimate the difference in structure between the thinned and unthinned forest. However, these observations are not enough alone to understand the effect of disturbance. We therefore also use mathematical models describing the carbon balance of our ecosystem, here we use the DALEC2 model of ecosystem carbon balance (Bloom and Williams, 2015). In order to find the best estimate for our system we use the mathematical technique of data assimilation in order to combine all our available observations with our prior model predictions. More infomation on the novel data assimilation techniques developed can be found in Pinnington et al., 2016. These techniques allow us to find two distinct parameter sets for the DALEC2 model corresponding to the thinned and unthinned forest. We can then inspect the model output for both areas of forest and attempt to further understand the effect of selective felling on ecosystem carbon dynamics.
In Figure 3 we show the cumulative fluxes for both the thinned and unthinned forest after disturbance in 2015. We would probably assume that removing 46% of the trees from the thinned section would reduce the amount of carbon uptake in comparison to the unthinned section. However, we can see that both forests removed a total of approximately 425 g C m-2 in 2015, despite the thinned forest having 46% of its trees removed in the previous year. From our best modelled predictions this unchanged carbon uptake is possible due to significant reductions in TER. So, even though the thinned forest has lower GPP, its net carbon uptake is similar to the unthinned forest. Our model suggests that GPP is a main driver for TER, therefore removing a large amount of trees has significantly reduced ecosystem respiration. This result is supported by other ecological studies (Heinemeyer et al., 2012, Högberg et al., 2001, Janssens et al., 2001). This has implications for future predictions of land surface carbon uptake and whether forests will continue to sequester atmospheric CO2 at similar rates, or if they will be limited by increased GPP leading to increased respiration.
Wilkinson, M. et al., 2012: Inter-annual variation of carbon uptake by a plantation oak woodland in south-eastern England. Biogeosciences, 9 (12), 5373–5389.
Wilkinson, M., et al., 2016: Effects of management thinning on CO2 exchange by a plantation oak woodland in south-eastern England. Biogeosciences, 13 (8), 2367–2378, doi: 10.5194/bg-13-2367-2016.
Bloom, A. A. and M. Williams, 2015: Constraining ecosystem carbon dynamics in a data-limited world: integrating ecological “common sense” in a model data fusion framework. Biogeosciences, 12 (5), 1299–1315, doi: 10.5194/bg-12-1299-2015.
Pinnington, E. M., et al., 2016: Investigating the role of prior and observation error correlations in improving a model forecast of forest carbon balance using four-dimensional variational data assimilation. Agricultural and Forest Meteorology, 228229, 299 – 314, doi: http://dx.doi.org/10.1016/j.agrformet.2016.07.006.
Heinemeyer, A., et al., 2012: Exploring the “overflow tap” theory: linking forest soil co2 fluxes and individual mycorrhizo- sphere components to photosynthesis. Biogeosciences, 9 (1), 79–95.
Högberg, P., et al., 2001: Large-scale forest girdling shows that current photosynthesis drives soil respiration. Nature, 411 (6839), 789–792.
Janssens, I. A., et al., 2001: Productivity overshadows temperature in determining soil and ecosystem respiration across european forests. Global Change Biology, 7 (3), 269–278, doi: 10.1046/j.1365-2486.2001.00412.x.
Mountains come in many shapes and sizes and as a result their dynamic impact on the atmospheric circulation spans a continuous range of physical and temporal scales. For example, large-scale orographic features, such as the Himalayas and the Rockies, deflect the atmospheric flow and, as a result of the Earth’s rotation, generate waves downstream that can remain fixed in space for long periods of time. These are known as stationary waves (see Nigam and DeWeaver (2002) for overview). They have an impact not only on the regional hydro-climate but also on the location and strength of the mid-latitude westerlies. On smaller physical scales, orography can generate gravity waves that act to transport momentum from the surface to the upper parts of the atmosphere (see Teixeira 2014), playing a role in the mixing of chemical species within the stratosphere.
Figure 1 shows an example of the resolved orography at different horizontal resolutions over the Himalayas. The representation of orography within models is complicated by the fact that, unlike other parameterized processes, such as clouds and convection, that are typically totally unresolved by the model, its effects are partly resolved by the dynamics of the model and the rest is accounted for by parameterization schemes.However, many parameters within these schemes are not well constrained by observations, if at all. The World Meteorological Organisation (WMO) Working Group on Numerical Experimentation (WGNE) performed an inter-model comparison focusing on the treatment of unresolved drag processes within models (Zadra et al. 2013). They found that while modelling groups generally had the same total amount of drag from various different processes, their partitioning was vastly different, as a result of the uncertainty in their formulation.
Climate models with typically low horizontal resolutions, resolve less of the Earth’s orography and are therefore more dependent on parameterization schemes. They also have large model biases in their climatological circulations when compared with observations, as well as exhibiting a similarly large spread about these biases. What is more, their projected circulation response to climate change is highly uncertain. It is therefore worth investigating the processes that contribute towards the spread in their climatological circulations and circulation response to climate change. The representation of orographic processes seem vital for the accurate simulation of the atmospheric circulation and yet, as discussed above, we find that there is a lot of uncertainty in their treatment within models that may be contributing to model uncertainty. These uncertainties in the orographic treatment come from two main sources:
Model Resolution: Models with different horizontal resolutions will have different resolved orography.
Parameterization Formulation: Orographic drag parameterization formulation varies between models.
The issue of model resolution was investigated in our recent study, van Niekerk et al. (2016). We showed that, in the Met Office Unified Model (MetUM) at climate model resolutions, the decrease in parameterized orographic drag that occurs with increasing horizontal resolution was not balanced by an increase in resolved orographic drag. The inability of the model to maintain an equivalent total (resolved plus parameterized) orographic drag across resolutions resulted in an increase in systematic model biases at lower resolutions identifiable over short timescales. This shows not only that the modelled circulation is non-robust to changes in resolution but also that the parameterization scheme is not performing in the same way as the resolved orography. We have highlighted the impact of parameterized and resolved orographic drag on model fidelity and demonstrated that there is still a lot of uncertainty in the way we treat unresolved orography within models. This further motivates the need to constrain the theory and parameters within orographic drag parameterization schemes.
Nigam, S., and E. DeWeaver, 2002: Stationary Waves (Orographic and Thermally Forced). Academic Press, Elsevier Science, London, 2121–2137 pp., doi:10.1016/B978-0-12-382225-3. 00381-9.
Air pollution is a major global problem, with the World Health Organisation recently linking 1 in 8 global deaths to this invisible problem. I say invisible, what air pollution may seem is an almost invisible problem. My PhD looks at some of the largest air pollutants, particulate matter PM10, which is still only 1/5th the width of a human hair in diameter!
My project looks at whether winter (December – February) UK PM10 concentration ([PM10]) exceedance events will change in frequency or composition in a future climate. To answer this question, a state of the art climate model is required. This model simulates the atmosphere only and is an iteration of the Met-Office HADGEM3 model. The climate simulation models a future 2050 under the RCP8.5 emissions scenario, the highest greenhouse-gas emission scenario considered in IPCC-AR5 (Riahi et al., 2011).
In an attempt to model PM10 in the climate model (a complex feat, currently tasked to the coupled UKCA model), we have idealised the problem, making the results much easier to understand. We have emitted chemically inert tracers in the model, which represent the key sources of PM10 throughout mainland Europe and the UK. The source regions identified were: West Poland, Po Valley, BENELUX and the UK. While the modelled tracers were shown to replicate observed PM10 well, albeit with inevitable sources of lost variability, they were primarily used to identify synoptic flow regimes influencing the UK. The motivation of this work is to determine whether the flow regimes that influence the UK during UK PM10 episodes, change in a future climate.
As we are unable to accurately replicate observed UK [PM10] within the model, we need to generate a proxy for UK [PM10] episodes. We chose to identify the synoptic meteorological conditions (synoptic scale ~ 1000 km) that result in UK air pollution episodes. We find that the phenomenon of atmospheric blocking in the winter months, in the Northeast Atlantic/ European region, provide the perfect conditions for PM10 accumulation in the UK. In the Northern Hemisphere winter, Rossby Wave Breaking (RWB) is the predominant precursor to atmospheric blocking (Woollings et al., 2008). RWB is the meridional overturning of air masses in the upper troposphere, so that warm/cold air is advected towards the pole/equator. The diagnostic chosen to detect RWB on is potential temperature (θ) on the potential vorticity = 2 Potential vorticity units surface, otherwise termed the dynamical tropopause. The advantages of using this diagnostic for detecting RWB have been outlined in this study’s first publication; Webber et al., (2016). Figure 1 illustrates this mechanism and the metric used to diagnose RWB, BI, introduced by Pelly and Hoskins (2003).
In Fig. 1 warm air is transported to the north of cold air to the south. This mechanism generates an anticyclone to the north of the centre of overturning (black circle in Fig 1) and a cyclone to the south. If the anticyclone to north becomes quasi-stationary, a blocking anticyclone is formed, which has been shown to generate conditions favourable for the accumulation of PM10.
To determine whether there exists a change in RWB frequency, due to climate change (a climate increment), the difference in RWB frequency between two simulations must be taken. The first of these is a free-running present day simulation, which provides us with the models representation of a present day atmosphere. The second is a future time-slice simulation, representative of the year 2050. Figure 2 shows the difference between the two simulations, with positive values representing an increase in RWB frequency in a future climate. The black contoured region corresponds to the region where the occurrence of RWB significantly increases UK [PM10].
RWB frequency anomalies within the black contoured region are of most importance within this study. Predominantly the RWB frequency anomaly, within the black contour, can be described as a negative frequency anomaly. However, there also exist heterogeneous RWB frequency anomalies within the contoured region. What is shown is that there is a tendency for RWB to occur further north and eastward in a future climate. These shifts in the regions of RWB occurrence influence a shift in the resulting flow regimes that influence the UK.
Climate shifts in flow regimes were analysed, however only for the most prominent subset of RWB events. RWB can be subset into cyclonic and anti-cyclonic RWB (CRWB and ACRWB respectively) and both have quite different impacts on UK [PM10] (Webber et al., 2016). ACRWB events are the most prominent RWB subset within the Northeast Atlantic/ European region (Weijenborg et al., 2012). Figure 1 represents ACRWB, with overturning occurring in a clockwise direction about the centre of overturning and these events were analysed for climate shifts in resultant flow regimes.
The analysis of climate flow regime shifts, provides the most interesting result of this study. We find that there exists a significant (p<0.05) increase in near European BENELUX tracer transport into the UK and a significant reduction of UK tracer accumulation, following ACRWB events. What we therefore see is that while in the future we see a reduction in the number of RWB and ACRWB events in a region most influential to UK [PM10], there also exists a robust shift in the resulting flow regime. Following ACRWB, there exists an increased tendency for the transport of European PM10 and decreased locally sourced [PM10] in the UK. Increased European transport may result in increased long-range transport of smaller and potentially more toxic (Gehring et al., 2013) PM2.5 particles from Europe.
Gehring, U., Gruzieva, O., Agius, R. M., Beelen, R., Custovic, A., Cyrys, J., Eeftens, M., Flexeder, C., Fuertes, E., Heinrich, J., Hoffmann, B., deJongste, J. C., Kerkhof, M., Klümper, C., Korek, M., Mölter, A., Schultz, E. S., Simpson, A.,Sugiri, D., Svartengren, M., von Berg, A., Wijga, A. H., Pershagen, G. and Brunekreef B.: Air Pollution Exposure and Lung Function in Children: The ESCAPE Project. Children’s Health Prespect, 121, 1357-1364, doi:10.1289/ehp.1306770 , 2013.
Pelly, J. L and Hoskins, B. J.: A New Perspective on Blocking. J. Atmos. Sci, 50, 743-755, doi: http://dx.doi.org/10.1175/1520- 0469(2003)060<0743:ANPOB>2.0.CO;2, 2003.
Riahi, K., Rao S., Krey, V., Cho, C., Chirkov, V., Fischer, G., Kindermann, G., Nakicenovic, N. and Rafaj, P.: RCP 8.5—A scenario of comparatively high greenhouse gas emissions. Climatic Change, 109, no. 1-2, 33-57, doi: 10.1007/s10584-011-0149-y, 2011.
Webber, C. P., Dacre, H. F., Collins, W. J., and Masato, G.: The Dynamical Impact of Rossby Wave Breaking upon UK PM10 Concentration. Atmos. Chem. and Phys. Discuss, doi; 10.5194/acp-2016-571, 2016.
Weijenborg, C., de Vries, H. and Haarsma, R. J.: On the direction of Rossby wave breaking in blocking. Climate Dynamics, 39, 2823- 2831, doi: 10.1007/s00382-012-1332-1, 2012.
Woollings, T. J., Hoskins, B. J., Blackburn, M. and Berrisford, P.: A new Rossby wave-breaking interpretation of the North Atlantic Oscillation. J. Atmos. Sci, 65, 609-626, doi: http://dx.doi.org/10.1175/2007JAS2347.1, 2008.