Microphysics and convection in the grey zone
|
|
- David Williamson
- 6 years ago
- Views:
Transcription
1 Microphysics and convection in the grey zone Luc Gerard RMIB, Av. Circulaire, B10 Brussels ABSTRACT We assess the behaviour of a limited area NWP model with different handlings of deep convection when varying the grid-box length from 8 to 1km. 1 Introduction When the grid-box length is close to the dimensions of deep-convection (DC) cells or systems, these are badly handled by the kind of DC-parametrization commonly used in coarser resolution models. On the other hand, trying to treat them explicitly (with no parametrization) is neither satisfactory. More specific approaches are required to run at these grey-zone resolutions. Here we compare the behaviour of different approaches in an academic case and in a real case. Four ways to handle deep convection are tried: no specific parametrization (nocp), the Bougeault diagnostic scheme (dia), the MT package (MT) (Gerard et al. 200) and our new complementary subgrid updraught approach (CSU). 2 Academic experiment and schemes features The academic setup is inspired from Weisman and Klemp (182) and uses their same vertical profile of temperature, moisture and wind over the entire cyclic domain. At the centre of the domain (and at a height equal to its radius) we add a potential-temperature perturbation with a vertical radius of 1400m, a horizontal radius of km, that we assume to be resolved by a 1km grid-length mesh. The perturbation is maximum (2K) at its centre, decreasing to zero at its border. The model has 85 isobar vertical levels and is run with a non-hydrostatic dynamics; the radiation scheme is switched off. The nocp 1-km experiment is taken as a reference. The nocp setup cannot produce precipitation at 8km resolution, and gives a too narrow and too intense pattern at 4km (Fig. 1, left column). The dia experiment (Fig. 2, left column) gives a (too weak) signal at 8km, but a much too wide and intense area at 4 and 2km, with an important contribution of the subgrid. The wind and the resolved circulation are further affected by the inadequate response of the DC parametrization. The MT approach was developed to produce a smoother response down to resolutions around 4km. It prevents that the subgrid (deep convection) and the resolved condensation parametrizations both feed on the same available moisture, by using a sequential organisation of the parametrizations. The deep convection scheme affects the resolved model variables through condensation and transport fluxes, and its condensation fluxes are combined with the ones from the resolved scheme, before being passed to the microphysical scheme. In addition, evolutive aspects are included by using prognostic variables for ECMWF Workshop on Parametrization of Clouds and Precipitation, 5-8 November 20 1
2 m10y0c_noc : 2010/1/1 z0:0 +2h m10y0c_de1 : 2010/1/1 z0:0 +2h max=5.0, mean= BOT t+00 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC S074WIND: 10.8<u< 0.0, 7.5<v< ; ff<10. m20y0c_noc : 2010/1/1 z0:0 +2h max=.8, mean= 4 BOT t+00 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC S074WIND: 10.4<u< 2.7, 7.7<v<.0; ff<10.5 m20y0c_de1 : 2010/1/1 z0:0 +2h max=.5, mean= 0. BOT t+00 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC S074WIND:.2<u< 2.0,.<v< 2.1; ff<10. m40y0c_noc : 2010/1/1 z0:0 +2h max=48.2, mean= 0. BOT t+00 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC S074WIND:.5<u< 0.8,.4<v< 1.; ff<10.5 m40y0c_de1 : 2010/1/1 z0:0 +2h max=8.0, mean= (e) BOT t+000 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC S074WIND: 7.2<u< 0., 2.8<v< 2.5; ff< 7.2 m80y0c_noc : 2010/1/1 z0:0 +2h max=22., mean= BOT t+000 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC S074WIND: 7.<u< 2.0, 4.<v< 0.8; ff< 7. m80y0c_de1 : 2010/1/1 z0:0 +2h max=, mean= 0.00 BOT t+00 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC S074WIND: 5.2<u< 4., 0.2<v< ; ff< max=.5, mean= 0.25 (h) BOT t+00 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC S074WIND: 7.<u< 1., 4.4<v< 1.; ff< 7. Figure 1: 2-hour accumulated parametrization, over the lower half of the symmetric domain. Left: nocp ; right: new CSU scheme. (a,b): x=1km; (c,d): 2km; (e,f): 4km; (g,h): 8km the updraught vertical velocity (non-hydrostatic vertical motion equation) and the updraught mesh fraction (prognostic closure). The microphysics considers subgrid cloud and precipitation fraction overlaps depending on the relative subgrid and resolved shares of the condensation. Finally, a memory of the convective cloud fraction (updraught + detrainment fraction) is also used from one time step to the next, to prevent a re-evaporation of not yet precipitated convective condensate that could occur if this was assumed to be distributed over the entire grid box. The MT scheme as used in the Alaro-0 model has brought a significant improvement in the operational forecast of cloud and precipitation; however, while allowing a relatively smooth evolution, the subgrid share does not diminish with decreasing grid-box 8 ECMWF Workshop on Parametrization of Clouds and Precipitation, 5-8 November 20
3 m20y0c_lsp : 2010/1/1 z0:0 +2h m20y0c_mt : 2010/1/1 z0:0 +2h max=4., mean= 0. BOT t+00 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC S074WIND: 10.<u< 0.7,.2<v< 1.7; ff<10.4 m40y0c_lsp : 2010/1/1 z0:0 +2h max=45.8, mean= 0.48 BOT t+00 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC S074WIND: 1<u< 0.8, 10.8<v< 2.5; ff<. m40y0c_mt : 2010/1/1 z0:0 +2h max=, mean= 0.0 BOT t+000 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC S074WIND:.<u< 1.4, 2.0<v< 2.2; ff<. m80y0c_lsp : 2010/1/1 z0:0 +2h max=4.5, mean= 0. BOT t+000 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC S074WIND: 8.4<u<, 2.<v< ; ff< 8.5 m80y0c_mt : 2010/1/1 z0:0 +2h max= 0.4, mean= 0.0 (e) BOT t+00 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC S074WIND: 5.2<u< 4., 0.2<v< 0.0; ff< max= 5.1, mean= 0.42 BOT t+00 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC S074WIND: 5.7<u<., 0.8<v< 0.; ff< 5.7 Figure 2: 2-hour accumulated parametrization, over the lower half of the symmetric domain. Left: dia ; right: MT scheme. (a,b): x=2km; (c,d): 4km; (e,f): 8km length. For the academic experiment, we observe (Fig. 2, right column) that the forecast is still excessive at 4km, and that the absence of an explicit triggering criterion allows an unwanted triggering (but with low precipitation) in an area around the main perturbation at the three resolutions. In order to produce a gradual extinction of the subgrid scheme while increasing the resolution, we had to develop a scheme with a set of additional features. Considering that with increasing resolution, the mean grid-box motion and resolved condensation represent an increasing fraction of the real subgrid updraught, the idea is that the subgrid scheme has to provide a decreasing contribution representing a complement to the resolved part. This complementary subgrid updraught (CSU) scheme keeps the MT sequential organisation, with a single call to the microphysical scheme, after combining subgrid and resolved condensation. The ascent properties are expressed as perturbations to the mean grid-box values, derived from the anelastic equations. The CSUperturbation ascent is then closed and confined in a single model-column, and affected by the mesh fraction and the subgrid environmental profile. Evolution equations are used for the perturbation updraught velocity and the mesh fraction, and they allow (unlike in MT) the gradual rise of the updraught top along several time steps. The CAPE closure also requires special care, since the grid-column CAPE normally decreases with increasing mesh fraction (CAPE would be zero when the updraught covers the entire grid-column area, hence we cannot write that the mesh fraction is proportional to the grid-column ECMWF Workshop on Parametrization of Clouds and Precipitation, 5-8 November 20
4 m20y0c_d1e1 : 2010/1/1 z0:0 +2h m40y0c_d1e1 : 2010/1/1 z0:0 +2h max=0., mean= 0. BOT t+00 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC max=7.8, mean= 0.42 BOT t+000 PREC EAU.CON+EAU.GEC+NEI.CON+NEI.GEC Figure : 2-hour accumulated surface precipitation, CSU scheme but using Kain 2004 s triggering criterion with fixed threshold, w 0 =0.02m/s. : x =2km; : 4km. CAPE). For this we try to evaluate the CAPE for an environment outside the updraughts. Finally a triggering criterion is introduced. Commonly used triggering criteria (e.g. Kain 2004, using a threshold of resolved vertical velocity w but also energetic criteria like Rogers and Fritsch 1, those based on CKE/TKE and CIN, or on moisture convergence) tend to produce more intense triggering with increasing resolution; and often an initial triggering has a positive feedback on further triggering in the neighbourhood, which can deteriorate the results in the vicinity of triggered ascents. In the present academic simulation, where CAPE is present over the entire domain, this phenomenon is critical. In the approach of Kain 2004, a mixed parcel is lifted up to its lifting condensation level (LCL) where it receives a buoyancy kick expressed by T v,kf = [ γ(w LCL w 0 min(1, z LCL z 0 ) ] 1/, 1 γ 0.01m s 1 K, z 0 = 2km, w 0 = 0.02m s 1. (1) If the parcel is then able to reach its LFC and continue upwards, the scheme is triggered. In (1) the buoyancy kick increases together with the mean grid-box velocity w, that should on its turn increase together with the updraught mesh fraction. Larger mesh fractions occur at higher resolution, but we can hardly find a general rule to increase the threshold w 0 in function of the resolution, because the variation of the mesh fraction depends on the granularity of the updraughts in the mesh. Fig. shows the total precipitation at 4km and 2km when using a fixed value of w 0. We developed an alternative criterion consistent with the CSU-spirit, based on on a threshold of resolved condensation, proportional to the inverse of the grid box area. This was used in the CSU experiment shown on Fig. 1 (right). In this case, the precipitation area appears consistent at the four resolutions, while the subgrid contribution is decreasing regularly with increasing resolution (not shown). Real case The model was tested with full physics at the same four resolutions on the thunderstorms of over Belgium. The runs are hydrostatic at 8 and 4km, and non-hydrostatic at 2 and 1km; there are 41 vertical hybrid levels. The nocp 1-km run is taken as a control, while the final reference is now the 1-hour accumulated radar image (Fig. 4). Using the CSU scheme (Fig. 4 a,b,d,e), the structure of the precipitation is quite similar at the four resolutions, while the grid-box averaged amounts decrease at coarser resolution. The CSU manages to suppress the spurious band of intense precipitation at the Northern edge of the domain appearing on the control. The number of kernels at 1km is closer to reality and we observe a slightly better North-South alignment for CSU on the Western side, which corresponds to the location of subgrid activity as can be seen on Fig ECMWF Workshop on Parametrization of Clouds and Precipitation, 5-8 November 20
5 t1de1 : 2005//10 z:0 +5h t2de1 : 2005//10 z:0 +5h max=8.2, mean= 1. CLSVENT: 5.5<u< 7.0, 7.<v< 7.0; ff< 8.2 max=7.4, mean= 1.02 CLSVENT:.8<u< 8.1, 5.2<v< 5.7; ff< 8.1 t4de1 : 2005//10 z:0 +5h t8de1 : 2005//10 z:0 +5h (e) t1nocxr : 2005//10 z:0 +5h max=, mean= 0. CLSVENT: 5.<u< 5.8,.2<v< 5.7; ff<.8 max=., mean= 0.1 CLSVENT: 2.5<u<.8, <v< 4.; ff< 4. max=8, mean= 1.57 CLSVENT:.5<u<.4,.<v< 7.7; ff<.8 Figure 4: Thunderstorms of over Belgium. 10-m wind and 1-hour accumulated precipitation, forecast using the CSU scheme at 1km, 2km, 4km and 8km (e) resolution. 1-hour accumulated radar image (c, same colour scale) and control simulation nocp at 1km resolution. Fig. 7 shows the results at 2, 4 and 8km resolution with no DC parametrization (nocp), with a diagnostic parametrization (dia) and with MT. At this particular range the nocp experiments does not show very intense grey-zone syndrome at 4 and 2km, but this is not the case for the whole simulation. The unwanted band of intense precipitation is still present on the nocp-2km, and there is still excessive activity in this region at 4km. The amounts at 8km are clearly exaggerated. The diagnostic parametrization produces very large areas of intense precipitation at 4km and 2km, mostly due to the DC scheme. At 8km the amounts remain much too big when considering the averaging effect. The alignment and location with the diagnostic parametrization is quite similar to the one with no parametrization. Here the MT response at 4km, while less extreme than the diagnostic scheme, is still excessive (the behaviour varies from one range to the next). The forecasters are generally happy to receive a warning even on a coarse resolution model grid, but this comes often together with an over-estimation of the precipitation that should represent an average over a large grid-box. In that case the results give no good estimation of the accumulation in hydrological basins. Fig. shows the accumulation over the entire area of Fig. 4. Unlike nocp and MT, the CSU scheme appears to produce (at the end of the event) the same total amount of accumulated water at all four resolutions. ECMWF Workshop on Parametrization of Clouds and Precipitation, 5-8 November 20 11
6 t1de1 : 2005//10 z:0 +5h t2de1 : 2005//10 z:0 +5h max=14.4, mean= 0.0 BOT t+0005 PREC EAU.CON+NEI.CON, 4 to 5 max=1., mean= 8 BOT t+0005 PREC EAU.CON+NEI.CON, 4 to 5 t4de1 : 2005//10 z:0 +5h t8de1 : 2005//10 z:0 +5h max=., mean= 0. BOT t+0005 PREC EAU.CON+NEI.CON, 4 to 5 max= 5., mean= 0.28 BOT t+0005 PREC EAU.CON+NEI.CON, 4 to 5 Figure 5: Same as Fig. 4 but subgrid share of the total precipitation, CSU scheme. 1km, 2km, 4km, 8km resolution. Figure : Evolution of accumulated surface precipitation. Red: 1km, green: 2km, blue: 4km, black: 8km. Solid line: total, dashed line: subgrid (DC-scheme) part. : nocp, : MT, : CSU simulations. 4 Conclusions It appears possible to make consistent forecasts at the grey-zone resolutions of the deep convection (between 8km and less than 1km, following the situation) by parametrizing aspects usually neglected at coarser resolution. The complementary subgrid updraught approach is specific for the triggering, ECMWF Workshop on Parametrization of Clouds and Precipitation, 5-8 November 20
7 the representation of the ascent, the evolutive aspects, the closure, the organisation with the other parametrizations. Its features make it compatible with the full range of model resolutions, from several tens of km to less than 1km. It produces a gradual transition to fully resolved convection when the resolution becomes sufficient. The complementarity with the resolved scheme can allow to maintain the total amount of precipitation feeding the hydrological basins independently of the resolution, which is an asset for the coupling with hydrological simulations. Further tests are on the way to verify the good model behaviour when changing some other aspects of the physics and in a wider set of situations. Acknowledgements This work benefited from enlightening exchanges in the frame of the COST ES005 action subsided by the European Community. References L. GERARD, J.-M. PIRIOU, R. BROŽKOVÁ, J.-F. GELEYN, and D. BANCIU. Cloud and precipitation parameterization in a meso-gamma scale operational weather prediction model. Mon.Wea.Rev, 17:0 77, November 200. J. S. KAIN. The Kain-Fritsch convective parameterization: an update. J. Appl. Meteor., 4:1 1, January Rogers, R. F. and J. M. Fritsch. A general framework for convective trigger functions. Mon. Wea. Rev.,1:8 52, 1. M. L. Weisman and J. B. Klemp. The dependence of numerically simulated convective storms on vertical wind shear and buoyancy. Mon.Wea.Rev, 110:4 520, 182. ECMWF Workshop on Parametrization of Clouds and Precipitation, 5-8 November 20 1
8 t2noc : 2005//10 z:0 +5h t4noc : 2005//10 z:0 +5h t8noc : 2005//10 z:0 +5h 8 max=7.2, mean= 1.47 CLSVENT:.1<u< 7., 7.<v< 7.4; ff< 8.0 max=47.1, mean= CLSVENT: 8.7<u< 8., 7.<v<11.4; ff<. max=2.7, mean= 1.1 CLSVENT:.7<u< 7.1, 7.2<v< 7.4; ff< 7.8 t2_lsp : 2005//10 z:0 +5h t4_lsp : 2005//10 z:0 +5h (e) t8_lsp : 2005//10 z:0 +5h max=4., mean= 2.08 CLSVENT:.<u<.8,.<v<.; ff<. max=, mean= 2.22 CLSVENT:.5<u<.5,.4<v< 7.4; ff< 7. max=8.4, mean= 1. CLSVENT: 4.<u<.,.<v< 5.1; ff< 7.2 t2_mt : 2005//10 z:0 +5h t4_mt : 2005//10 z:0 +5h (g) (h) (i) t8_mt : 2005//10 z:0 +5h max=4.7, mean= 2 CLSVENT: <u< 7.,.<v< 7.7; ff< 8. max=45.0, mean= 8 CLSVENT:.<u<.2, 5.<v< 7.0; ff<. max=., mean= 1.10 CLSVENT: 4.7<u< 5., 5.5<v< 2.7; ff<.0 Figure 7: Same as Fig. 4 but for nocp experiment at 2km, 4km, 8km, dia experiment at 2km, 4km (e), 8km and MT experiment at 2km (g), 4km (h), 8km (i) resolution. 14 ECMWF Workshop on Parametrization of Clouds and Precipitation, 5-8 November 20
Neva Pristov. LACE area leader for physics
ALARO Physics Developments Neva Pristov LACE area leader for physics Main developments TOUCANS turbulence scheme Radiation Convection Ivan Baštak Duran, Jean-Francois Geleyn, Filip Vana Jan Mašek, Radmila
More informationScale aware deep convection parameterization
Scale aware deep convection parameterization Luc Gerard Royal Meteorological Institute of Belgium 3 September 26 Convective clouds in a model grid box Quasi-Equilibrium hypothesis: Large subgrid population,
More informationALARO-0 at the RMI. Alex, Daan, Geert, Josette, Luc, Olivier, Rafiq, Piet( s first operational talk)
ALARO-0 at the RMI Alex, Daan, Geert, Josette, Luc, Olivier, Rafiq, Piet( s first operational talk) Operational use of Alaro-0 0 at RMIB The model (cy35t1) is run 4 times a day, at: 1. Resolution 4km (square
More informationImpact of different cumulus parameterizations on the numerical simulation of rain over southern China
Impact of different cumulus parameterizations on the numerical simulation of rain over southern China P.W. Chan * Hong Kong Observatory, Hong Kong, China 1. INTRODUCTION Convective rain occurs over southern
More informationSimulating orographic precipitation: Sensitivity to physics parameterizations and model numerics
Simulating orographic precipitation: Sensitivity to physics parameterizations and model numerics 2nd COPS-Meeting, 27 June 2005 Günther Zängl Overview A highly idealized test of numerical model errors
More informationMultiscale performance of the ALARO-0 model for simulating extreme summer precipitation climatology in Belgium
Multiscale performance of the ALARO-0 model for simulating extreme summer precipitation climatology in Belgium Rozemien De Troch 1,2, Rafiq Hamdi 1, Hans Van de Vyver 1, Jean-François Geleyn 2,3, Piet
More informationSami Niemelä and Carl Fortelius Finnish Meteorological Institute
6B.6 APPLICABILITY OF GRID-SIZE-DEPENDENT CONVECTION PARAMETERIZATION TO MESO-γ-SCALE HIRLAM. Sami Niemelä and Carl Fortelius Finnish Meteorological Institute. INTRODUCTION The representation of convective
More informationBio diversity and sustainable development. In European NWP
Bio diversity and sustainable development In European NWP ALADIN Consortium and development needs Some facts: We are network and so development is organized as such: We need good rules to be able to work
More informationThe prognostic deep convection parametrization for operational forecast in horizontal resolutions of 8, 4 and 2 km
The prognostic deep convection parametrization for operational forecast in horizontal resolutions of 8, 4 and 2 km Martina Tudor, Stjepan Ivatek-Šahdan and Antonio Stanešić tudor@cirus.dhz.hr Croatian
More information9D.3 THE INFLUENCE OF VERTICAL WIND SHEAR ON DEEP CONVECTION IN THE TROPICS
9D.3 THE INFLUENCE OF VERTICAL WIND SHEAR ON DEEP CONVECTION IN THE TROPICS Ulrike Wissmeier, Robert Goler University of Munich, Germany 1 Introduction One does not associate severe storms with the tropics
More informationModified PM09 parameterizations in the shallow convection grey zone
Modified PM09 parameterizations in the shallow convection grey zone LACE stay report Toulouse Centre National de Recherche Meteorologique, 02. February 2015 27. February 2015 Scientific supervisor: Rachel
More informationMEA 716 Exercise, BMJ CP Scheme With acknowledgements to B. Rozumalski, M. Baldwin, and J. Kain Optional Review Assignment, distributed Th 2/18/2016
MEA 716 Exercise, BMJ CP Scheme With acknowledgements to B. Rozumalski, M. Baldwin, and J. Kain Optional Review Assignment, distributed Th 2/18/2016 We have reviewed the reasons why NWP models need to
More informationParameterization of effects of unresolved clouds and precipitation
Parameterization of effects of unresolved clouds and precipitation eas471_cumparam.odp JDW, EAS, U. Alberta Last modified: 29 Mar. 2016 (from Physical Parameterizations in Canadian Operational Models,
More informationPrecipitation Structure and Processes of Typhoon Nari (2001): A Modeling Propsective
Precipitation Structure and Processes of Typhoon Nari (2001): A Modeling Propsective Ming-Jen Yang Institute of Hydrological Sciences, National Central University 1. Introduction Typhoon Nari (2001) struck
More informationPhysics/Dynamics Coupling at very high resolution:
Physics/Dynamics Coupling at very high resolution: permitted versus parametrized convection Sylvie Malardel ECMWF, Shinfield Park, Reading RG2 9AX, United Kingdom Sylvie.Malardel@ecmwf.int 1 Introduction
More informationTHE NON-SATURATED DOWNDRAFT IN ALARO-1
THE NON-SATURATED DOWNDRAFT IN ALARO-1 RC-LACE stay report CHMI, Prague, 5.06-3.07.2016 Simona Briceag (National Meteorological Administration, Romania) Supervisor: Radmila Brožková (CHMI) The proposed
More informationFor the operational forecaster one important precondition for the diagnosis and prediction of
Initiation of Deep Moist Convection at WV-Boundaries Vienna, Austria For the operational forecaster one important precondition for the diagnosis and prediction of convective activity is the availability
More informationRepresenting deep convective organization in a high resolution NWP LAM model using cellular automata
Representing deep convective organization in a high resolution NWP LAM model using cellular automata Lisa Bengtsson-Sedlar SMHI ECMWF, WMO/WGNE, WMO/THORPEX and WCRP WS on Representing model uncertainty
More informationINVESTIGATION FOR A POSSIBLE INFLUENCE OF IOANNINA AND METSOVO LAKES (EPIRUS, NW GREECE), ON PRECIPITATION, DURING THE WARM PERIOD OF THE YEAR
Proceedings of the 13 th International Conference of Environmental Science and Technology Athens, Greece, 5-7 September 2013 INVESTIGATION FOR A POSSIBLE INFLUENCE OF IOANNINA AND METSOVO LAKES (EPIRUS,
More informationImplementation and validation of the. ECMWF IFS convection scheme. in COSMO-CLM. Peter Brockhaus. Daniel Lüthi. Christoph Schär
Implementation and validation of the ECMWF IFS convection scheme in COSMO-CLM Peter Brockhaus Peter Bechtold Daniel Lüthi Oliver Fuhrer Christoph Schär (ETH) (ECMWF) (ETH) (MeteoSwiss) (ETH) COSMO/CLM
More informationParametrizing Cloud Cover in Large-scale Models
Parametrizing Cloud Cover in Large-scale Models Stephen A. Klein Lawrence Livermore National Laboratory Ming Zhao Princeton University Robert Pincus Earth System Research Laboratory November 14, 006 European
More informationNWP Equations (Adapted from UCAR/COMET Online Modules)
NWP Equations (Adapted from UCAR/COMET Online Modules) Certain physical laws of motion and conservation of energy (for example, Newton's Second Law of Motion and the First Law of Thermodynamics) govern
More informationECMWF Workshop on "Parametrization of clouds and precipitation across model resolutions
ECMWF Workshop on "Parametrization of clouds and precipitation across model resolutions Themes: 1. Parametrization of microphysics 2. Representing sub-grid cloud variability 3. Constraining cloud and precipitation
More information11A.2 Forecasting Short Term Convective Mode And Evolution For Severe Storms Initiated Along Synoptic Boundaries
11A.2 Forecasting Short Term Convective Mode And Evolution For Severe Storms Initiated Along Synoptic Boundaries Greg L. Dial and Jonathan P. Racy Storm Prediction Center, Norman, Oklahoma 1. Introduction
More informationSimulations with different convection parameterizations in the LM
Simulations with different convection parameterizations in the LM Linda Smoydzin Almut Gassmann Andreas Bott Marco Arpagaus (Meteo Swiss) Meteorological Institute of the University of Bonn, Germany Aims
More informationMODEL TYPE (Adapted from COMET online NWP modules) 1. Introduction
MODEL TYPE (Adapted from COMET online NWP modules) 1. Introduction Grid point and spectral models are based on the same set of primitive equations. However, each type formulates and solves the equations
More informationAn integrated package for subgrid convection, clouds and precipitation compatible with meso-gamma scales
QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY Q. J. R. Meteorol. Soc. 1: 711 7 (007) Published online in Wiley InterScience (www.interscience.wiley.com).58 An integrated package for subgrid convection,
More informationDevelopment of new products by operational forecasters
Development of new products by operational forecasters Probability of Thunder -algorithm Paavo Korpela Operational meteorologist & developer Safety Weather Centre Finnish Meteorological Institute Content
More informationPALM - Cloud Physics. Contents. PALM group. last update: Monday 21 st September, 2015
PALM - Cloud Physics PALM group Institute of Meteorology and Climatology, Leibniz Universität Hannover last update: Monday 21 st September, 2015 PALM group PALM Seminar 1 / 16 Contents Motivation Approach
More informationThe perturbation pressure, p, can be represented as the sum of a hydrostatic pressure perturbation p h and a nonhydrostatic pressure perturbation p nh
z = The perturbation pressure, p, can be represented as the sum of a hydrostatic pressure perturbation p h and a nonhydrostatic pressure perturbation p nh, that is, p = p h + p nh. (.1) The former arises
More informationAn integrated package for subgrid convection, clouds and precipitation compatible with meso-gamma scales.
QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY Published online in Wiley InterScience (www.interscience.wiley.com).06.42 An integrated package for subgrid convection, clouds and precipitation compatible
More informationThe Effect of Sea Spray on Tropical Cyclone Intensity
The Effect of Sea Spray on Tropical Cyclone Intensity Jeffrey S. Gall, Young Kwon, and William Frank The Pennsylvania State University University Park, Pennsylvania 16802 1. Introduction Under high-wind
More information1.3 HIGH-RESOLUTION MESOSCALE SIMULATIONS ON THE ROLE OF SHALLOW AND DEEP CONVECTION ON DUST EMISSION AND TRANSPORT IN A DESERT AREA.
1.3 HIGH-RESOLUTION MESOSCALE SIMULATIONS ON THE ROLE OF SHALLOW AND DEEP CONVECTION ON DUST EMISSION AND TRANSPORT IN A DESERT AREA Tetsuya Takemi Dept. of Environmental Science and Technology, Tokyo
More informationPrediction of clouds and precipitation in atmospheric models
Prediction of clouds and precipitation in atmospheric models Cumulus congestus. Sønderjylland, 2002. Foto: Flemming Byg. B. H. Sass NetFam Summer School July 2008 outline Basic definitions and characteristics
More informationNesting and LBCs, Predictability and EPS
Nesting and LBCs, Predictability and EPS Terry Davies, Dynamics Research, Met Office Nigel Richards, Neill Bowler, Peter Clark, Caroline Jones, Humphrey Lean, Ken Mylne, Changgui Wang copyright Met Office
More informationInteraction of boundary-layer turbulence and convection in the AROME model: Preliminary results
Interaction of boundary-layer turbulence and convection in the AROME model: Preliminary results Balázs Szintai Hungarian Meteorological Service COST ES0905 Workshop Cambridge, England 23 March 2010 Outline
More informationWeather Research and Forecasting Model. Melissa Goering Glen Sampson ATMO 595E November 18, 2004
Weather Research and Forecasting Model Melissa Goering Glen Sampson ATMO 595E November 18, 2004 Outline What does WRF model do? WRF Standard Initialization WRF Dynamics Conservation Equations Grid staggering
More informationA Study of Convective Initiation Failure on 22 Oct 2004
A Study of Convective Initiation Failure on 22 Oct 2004 Jennifer M. Laflin Philip N. Schumacher NWS Sioux Falls, SD August 6 th, 2011 Introduction Forecasting challenge: strong forcing for ascent and large
More informationEvolution and Maintenance of the June 2003 Nocturnal Convection
Evolution and Maintenance of the 22-23 June 2003 Nocturnal Convection Jerilyn Billings NOAA/NWS Wichita, KS August 6 th, 2011 Work Completed at North Carolina State University for MS Thesis During the
More informationANALYSIS OF THE MPAS CONVECTIVE-PERMITTING PHYSICS SUITE IN THE TROPICS WITH DIFFERENT PARAMETERIZATIONS OF CONVECTION REMARKS AND MOTIVATIONS
ANALYSIS OF THE MPAS CONVECTIVE-PERMITTING PHYSICS SUITE IN THE TROPICS WITH DIFFERENT PARAMETERIZATIONS OF CONVECTION Laura D. Fowler 1, Mary C. Barth 1, K. Alapaty 2, M. Branson 3, and D. Dazlich 3 1
More informationz g + F w (2.56) p(x, y, z, t) = p(z) + p (x, y, z, t) (2.120) ρ(x, y, z, t) = ρ(z) + ρ (x, y, z, t), (2.121)
= + dw dt = 1 ρ p z g + F w (.56) Let us describe the total pressure p and density ρ as the sum of a horizontally homogeneous base state pressure and density, and a deviation from this base state, that
More informationPV Generation in the Boundary Layer
1 PV Generation in the Boundary Layer Robert Plant 18th February 2003 (With thanks to S. Belcher) 2 Introduction How does the boundary layer modify the behaviour of weather systems? Often regarded as a
More informationA brief overview of the scheme is given below, taken from the whole description available in Lopez (2002).
Towards an operational implementation of Lopez s prognostic large scale cloud and precipitation scheme in ARPEGE/ALADIN NWP models F.Bouyssel, Y.Bouteloup, P. Marquet Météo-France, CNRM/GMAP, 42 av. G.
More informationChapter 3 Convective Dynamics
Chapter 3 Convective Dynamics Photographs Todd Lindley 3.2 Ordinary or "air-mass storm 3.2.1. Main Characteristics Consists of a single cell (updraft/downdraft pair) Forms in environment characterized
More information1. Current atmospheric DA systems 2. Coupling surface/atmospheric DA 3. Trends & ideas
1 Current issues in atmospheric data assimilation and its relationship with surfaces François Bouttier GAME/CNRM Météo-France 2nd workshop on remote sensing and modeling of surface properties, Toulouse,
More informationLarge-Eddy Simulations of Tropical Convective Systems, the Boundary Layer, and Upper Ocean Coupling
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Large-Eddy Simulations of Tropical Convective Systems, the Boundary Layer, and Upper Ocean Coupling Eric D. Skyllingstad
More informationVariational data assimilation of lightning with WRFDA system using nonlinear observation operators
Variational data assimilation of lightning with WRFDA system using nonlinear observation operators Virginia Tech, Blacksburg, Virginia Florida State University, Tallahassee, Florida rstefane@vt.edu, inavon@fsu.edu
More informationTornadoes forecasting, dynamics and genesis. Mteor 417 Iowa State University Week 12 Bill Gallus
Tornadoes forecasting, dynamics and genesis Mteor 417 Iowa State University Week 12 Bill Gallus Tools to diagnose severe weather risks Definition of tornado: A vortex (rapidly rotating column of air) associated
More informationChapter 14 Thunderstorm Fundamentals
Chapter overview: Thunderstorm appearance Thunderstorm cells and evolution Thunderstorm types and organization o Single cell thunderstorms o Multicell thunderstorms o Orographic thunderstorms o Severe
More informationNumerical simulation of marine stratocumulus clouds Andreas Chlond
Numerical simulation of marine stratocumulus clouds Andreas Chlond Marine stratus and stratocumulus cloud (MSC), which usually forms from 500 to 1000 m above the ocean surface and is a few hundred meters
More informationMesoscale meteorological models. Claire L. Vincent, Caroline Draxl and Joakim R. Nielsen
Mesoscale meteorological models Claire L. Vincent, Caroline Draxl and Joakim R. Nielsen Outline Mesoscale and synoptic scale meteorology Meteorological models Dynamics Parametrizations and interactions
More informationIncorporation of 3D Shortwave Radiative Effects within the Weather Research and Forecasting Model
Incorporation of 3D Shortwave Radiative Effects within the Weather Research and Forecasting Model W. O Hirok and P. Ricchiazzi Institute for Computational Earth System Science University of California
More informationOperational convective scale NWP in the Met Office
Operational convective scale NWP in the Met Office WSN09 Symposium. 18 st of May 2011 Jorge Bornemann (presenting the work of several years by many Met Office staff and collaborators) Contents This presentation
More informationExperiments with Idealized Supercell Simulations
Latent Heat Nudging in almo: Experiments with Idealized Supercell Simulations Daniel Leuenberger and Andrea Rossa MeteoSwiss Talk Outline Introduction Supercell Storms Characteristics Reference simulation
More informationPhysics/Dynamics coupling
Physics/Dynamics coupling Sylvie Malardel ECMWF November 8, 2010 Sylvie Malardel (ECMWF) Physics/Dynamics coupling November 8, 2010 1 / 21 Coupling between Physics and Dynamics for convection permitting
More informationParameterization of Convection in the Global NWP System GME of the German Weather Service
Parameterization of Convection in the Global NWP System GME of the German Weather Service Dmitrii Mironov and Bodo Ritter German Weather Service, Offenbach am Main, Germany (dmitrii.mironov@dwd.de, bodo.ritter@dwd.de)
More information4.4 DRIZZLE-INDUCED MESOSCALE VARIABILITY OF BOUNDARY LAYER CLOUDS IN A REGIONAL FORECAST MODEL. David B. Mechem and Yefim L.
4.4 DRIZZLE-INDUCED MESOSCALE VARIABILITY OF BOUNDARY LAYER CLOUDS IN A REGIONAL FORECAST MODEL David B. Mechem and Yefim L. Kogan Cooperative Institute for Mesoscale Meteorological Studies University
More informationPreliminary results with very high resolution COSMO model for the forecast of convective events. Antonella Morgillo. Arpa-Simc.
2 Working Group on Physical Aspects 52 Preliminary results with very high resolution COSMO model for the forecast of convective events. Antonella Morgillo Arpa-Simc amorgillo@arpa.emr.it 1 Introduction
More informationSuccesses and challenges of a seamless development of physical parameterizations
Successes and challenges of a seamless development of physical parameterizations F. Bouyssel, I. Beau, E. Bazile, Y. Bouteloup, J.F. Gueremy, J.M. Piriou, D. Pollack, R. Roehrig CNRM-GAME, Météo-France
More informationTheories on the Optimal Conditions of Long-Lived Squall Lines
Theories on the Optimal Conditions of Long-Lived Squall Lines References: Thorpe, A. J., M. J. Miller, and M. W. Moncrieff, 1982: Two -dimensional convection in nonconstant shear: A model of midlatitude
More informationCharles A. Doswell III, Harold E. Brooks, and Robert A. Maddox
Charles A. Doswell III, Harold E. Brooks, and Robert A. Maddox Flash floods account for the greatest number of fatalities among convective storm-related events but it still remains difficult to forecast
More informationWaVaCS summerschool Autumn 2009 Cargese, Corsica
Introduction Part I WaVaCS summerschool Autumn 2009 Cargese, Corsica Holger Tost Max Planck Institute for Chemistry, Mainz, Germany Introduction Overview What is a parameterisation and why using it? Fundamentals
More informationBOUNDARY LAYER ENERGY POTENTIAL (BLEP): AN OPERATIONAL TOOL TO ESTIMATE MAXIMUM SURFACE WIND SPEEDS IN CONVECTIVE STORMS?
BOUNDARY LAYER ENERGY POTENTIAL (BLEP): AN OPERATIONAL TOOL TO ESTIMATE MAXIMUM SURFACE WIND SPEEDS IN CONVECTIVE STORMS? Bill Nichols Lead Forecaster, WFO DVN Quad Cities OVERVIEW Purpose Of BLEP General
More informationPrecipitation in climate modeling for the Mediterranean region
Precipitation in climate modeling for the Mediterranean region Simon Krichak Dept. of Geophysics Atmospheric and Planetary Sciences, Tel Aviv University, Israel Concepts for Convective Parameterizations
More informationDiagnostic tools using a mesoscale NWP model for the early warning of convection
Meteorol. Appl. 5, 329 349 (1998) Diagnostic tools using a mesoscale NWP model for the early warning of convection Véronique Ducrocq, Diane Tzanos and Stéphane Sénési, GAME/CNRM, Météo-France, 31057 Toulouse
More informationRepresenting model uncertainty using multi-parametrisation methods
Representing model uncertainty using multi-parametrisation methods L.Descamps C.Labadie, A.Joly and E.Bazile CNRM/GMAP/RECYF 2/17 Outline Representation of model uncertainty in EPS The multi-parametrisation
More informationCONVECTIVE CLOUD MICROPHYSICS IN A HIGH-RESOLUTION NWP MODEL
CONVECTIVE CLOUD MICROPHYSICS IN A HIGH-RESOLUTION NWP MODEL J. Trentmann 1, A. Seifert 2, H. Wernli 1 1 Institute for Atmospheric Physics, Johannes Gutenberg University Mainz, Germany 2 German Weather
More informationAtmospheric Processes
Atmospheric Processes Atmospheric prognostic variables Wind Temperature Humidity Cloud Water/Ice Atmospheric processes Mixing Radiation Condensation/ Evaporation Precipitation Surface exchanges Friction
More informationKelly Mahoney NOAA ESRL Physical Sciences Division
The role of gray zone convective model physics in highresolution simulations of the 2013 Colorado Front Range Flood WRF model simulated precipitation over terrain in CO Front Range Kelly Mahoney NOAA ESRL
More informationSPECIAL PROJECT PROGRESS REPORT
SPECIAL PROJECT PROGRESS REPORT Progress Reports should be 2 to 10 pages in length, depending on importance of the project. All the following mandatory information needs to be provided. Reporting year
More informationEnhancing information transfer from observations to unobserved state variables for mesoscale radar data assimilation
Enhancing information transfer from observations to unobserved state variables for mesoscale radar data assimilation Weiguang Chang and Isztar Zawadzki Department of Atmospheric and Oceanic Sciences Faculty
More informationSimulation of shallow cumuli and their transition to deep convective clouds by cloud-resolving models with different third-order turbulence closures
Q. J. R. Meteorol. Soc. (2006), 132, pp. 359 382 doi: 10.1256/qj.05.29 Simulation of shallow cumuli and their transition to deep convective clouds by cloud-resolving models with different third-order turbulence
More informationSensitivity of precipitation forecasts to cumulus parameterizations in Catalonia (NE Spain)
Sensitivity of precipitation forecasts to cumulus parameterizations in Catalonia (NE Spain) Jordi Mercader (1), Bernat Codina (1), Abdelmalik Sairouni (2), Jordi Cunillera (2) (1) Dept. of Astronomy and
More informationDepartment of Geosciences San Francisco State University Spring Metr 201 Monteverdi Quiz #5 Key 150 pts.
Department of Geosciences Name San Francisco State University Spring 2013 Metr 201 Monteverdi Quiz #5 Key 150 pts. 1. Definitions. (5 points each for a total of 20 points in this section). (a) Coriolis
More informationTornadoes forecasting, dynamics and genesis. Mteor 417 Iowa State University Week 12 Bill Gallus
Tornadoes forecasting, dynamics and genesis Mteor 417 Iowa State University Week 12 Bill Gallus Tools to diagnose severe weather risks Definition of tornado: A vortex (rapidly rotating column of air) associated
More informationDiabatic processes and the structure of extratropical cyclones
Geophysical and Nonlinear Fluid Dynamics Seminar AOPP, Oxford, 23 October 2012 Diabatic processes and the structure of extratropical cyclones Oscar Martínez-Alvarado R. Plant, J. Chagnon, S. Gray, J. Methven
More informationComparison of Convection-permitting and Convection-parameterizing Ensembles
Comparison of Convection-permitting and Convection-parameterizing Ensembles Adam J. Clark NOAA/NSSL 18 August 2010 DTC Ensemble Testbed (DET) Workshop Introduction/Motivation CAMs could lead to big improvements
More informationA Large-Eddy Simulation Study of Moist Convection Initiation over Heterogeneous Surface Fluxes
A Large-Eddy Simulation Study of Moist Convection Initiation over Heterogeneous Surface Fluxes Song-Lak Kang Atmospheric Science Group, Texas Tech Univ. & George H. Bryan MMM, NCAR 20 th Symposium on Boundary
More informationTowards a new hybrid cumulus parametrization scheme for use in non-hydrostatic weather prediction models
QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY Q. J. R. Meteorol. Soc. 133: 479 49 (27) Published online in Wiley InterScience (www.interscience.wiley.com) DOI: 1.12/qj.28 Towards a new hybrid cumulus
More informationReading. What meteorological conditions precede a thunderstorm? Thunderstorms: ordinary or single cell storms, multicell storms, supercell storms
Thunderstorms: ordinary or single cell storms, multicell storms, supercell storms Reading Ahrens, Chapter 14: Thunderstorms and Tornadoes This lecture + next (Lightning, tornadoes) will cover the topic.
More informationmeso-sas, a modification of the SAS for meso-scale models Hua-Lu pan Qingfu Liu
meso-sas, a modification of the SAS for meso-scale models Hua-Lu pan Qingfu Liu Problem with the conventional mass flux schemes Most of the mass-flux schemes are based on the original Arakawa-Schubert
More informationSIMULATED EFFECTS OF AN ISOLATED SUPERCELL ON THE EVOLUTION OF A NEARBY SQUALL LINE
5.55 SIMULATED EFFECTS OF AN ISOLATED SUPERCELL ON THE EVOLUTION OF A NEARBY SQUALL LINE Jacey Wipf* and Adam French South Dakota School of Mines and Technology 1. INTRODUCTION 2. METHODOLOGY Operational
More informationThe Contribution Of Fine Scale Atmospheric Numerical Models In Improving The Quality Of Hydraulic Modelling Outputs
City University of New York (CUNY) CUNY Academic Works International Conference on Hydroinformatics 8-1-2014 The Contribution Of Fine Scale Atmospheric Numerical Models In Improving The Quality Of Hydraulic
More informationMidlatitude Frontal Clouds: GCM-Scale Modeling Implications
2729 Midlatitude Frontal Clouds: GCM-Scale Modeling Implications J. J. KATZFEY AND B. F. RYAN CSIRO Atmospheric Research, Aspendale, Victoria, Australia (Manuscript received 8 March 1999, in final form
More informationCumulus parameterization in non-convection-resolving models
Cumulus parameterization in non-convection-resolving models Given a column profile of model variables*, what convective tendencies will* occur? Hard questions: *1 is mean thermo. sounding enough information?»if
More informationTHE IMPACT OF GROUND-BASED GPS SLANT-PATH WET DELAY MEASUREMENTS ON SHORT-RANGE PREDICTION OF A PREFRONTAL SQUALL LINE
JP1.17 THE IMPACT OF GROUND-BASED GPS SLANT-PATH WET DELAY MEASUREMENTS ON SHORT-RANGE PREDICTION OF A PREFRONTAL SQUALL LINE So-Young Ha *1,, Ying-Hwa Kuo 1, Gyu-Ho Lim 1 National Center for Atmospheric
More informationExploring stochastic model uncertainty representations
Exploring stochastic model uncertainty representations with relevance to the greyzone Sarah-Jane Lock, Martin Leutbecher, Peter Bechtold, Richard Forbes Research Department, ECMWF ECMWF November 15, 2017
More informationFuture directions for parametrization of cloud and precipitation microphysics
Future directions for parametrization of cloud and precipitation microphysics Richard Forbes (ECMWF) ECMWF-JCSDA Workshop, 15-17 June 2010 Cloud and Precipitation Microphysics A Complex System! Ice Nucleation
More informationAnalysis of Cloud-Radiation Interactions Using ARM Observations and a Single-Column Model
Analysis of Cloud-Radiation Interactions Using ARM Observations and a Single-Column Model S. F. Iacobellis, R. C. J. Somerville, D. E. Lane, and J. Berque Scripps Institution of Oceanography University
More informationTornadogenesis in Supercells: The Three Main Ingredients. Ted Funk
Tornadogenesis in Supercells: The Three Main Ingredients Ted Funk NWS Louisville, KY Spring 2002 Environmental Parameters Supercells occur within environments exhibiting several wellknown characteristics
More informationPrediction of cirrus clouds in GCMs
Prediction of cirrus clouds in GCMs Bernd Kärcher, Ulrike Burkhardt, Klaus Gierens, and Johannes Hendricks DLR Institut für Physik der Atmosphäre Oberpfaffenhofen, 82234 Wessling, Germany bernd.kaercher@dlr.de
More informationDevelopment of a stochastic convection scheme
Development of a stochastic convection scheme R. J. Keane, R. S. Plant, N. E. Bowler, W. J. Tennant Development of a stochastic convection scheme p.1/44 Outline Overview of stochastic parameterisation.
More information11.B4 The influence of horizontal convective rolls on the morphology of low-level rotation in idealized simulations of supercell thunderstorms
11.B4 The influence of horizontal convective rolls on the morphology of low-level rotation in idealized simulations of supercell thunderstorms CHRISTOPHER J. NOWOTARSKI, PAUL M. MARKOWSKI, AND YVETTE P.
More informationThe Plant-Craig stochastic Convection scheme in MOGREPS
The Plant-Craig stochastic Convection scheme in MOGREPS R. J. Keane R. S. Plant W. J. Tennant Deutscher Wetterdienst University of Reading, UK UK Met Office Keane et. al. (DWD) PC in MOGREPS / 6 Overview
More informationIdealized Nonhydrostatic Supercell Simulations in the Global MPAS
Idealized Nonhydrostatic Supercell Simulations in the Global Joe Klemp, Bill Skamarock, and Sang-Hun Park National Center for Atmospheric Research Boulder, Colorado Typical characteristics: Supercell Thunderstorms
More informationWhere does the memory of convection stem from? Why can it be useful for parameterizations?
Where does the memory of convection stem from? Why can it be useful for parameterizations? D'où vient la mémoire de la convection? En quoi cela peut-il être utile pour les paramétrisations? Maxime Colin,
More informationAviation Hazards: Thunderstorms and Deep Convection
Aviation Hazards: Thunderstorms and Deep Convection TREND Empirical thunderstorm forecasting techniques Contents Necessary conditions for convection: Instability Low-level moisture Trigger mechanism Forecasting
More informationEnvironmental stability control of the intensity of squall lines under low-level shear conditions
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2007jd008793, 2007 Environmental stability control of the intensity of squall lines under low-level shear conditions Tetsuya Takemi 1 Received 12
More informationFigure 5: Comparison between SAFIR warning and radar-based hail detection for the hail event of June 8, 2003.
SAFIR WARNING : Expected risk Radar-based Probability of Hail 0915 0930 0945 1000 Figure 5: Comparison between SAFIR warning and radar-based hail detection for the hail event of June 8, 2003. Lightning
More informationJP1J.9 ASSIMILATION OF RADAR DATA IN THE MET OFFICE MESOSCALE AND CONVECTIVE SCALE FORECAST SYSTEMS
JP1J.9 ASSIMILATION OF RADAR DATA IN THE MET OFFICE MESOSCALE AND CONVECTIVE SCALE FORECAST SYSTEMS S. Ballard *, M. Dixon, S. Swarbrick,, Z. Li,, O.Stiller, H. Lean, F. Rihan 2 and C.Collier 2 Met Office
More information