Media: improvement of the source term treatment F. Bompay, L. Borrel Mefeo France, 6'CEM/06V16', Tb^/o^^e, France
|
|
- Gervais Jenkins
- 5 years ago
- Views:
Transcription
1 Media: improvement of the source term treatment F. Bompay, L. Borrel Mefeo France, 6'CEM/06V16', Tb^/o^^e, France ABSTRACT MEDIA is an Eulerian three dimensional model of atmospheric dispersion, used at the present time by METEO-FRANCE to forecast the behaviour at middle and large scale, of an atmospheric release of pollutant. In this model, the source is described on the horizontal plan by a Gaussian distribution between the point of release and the 12 closest grid points.this scheme does not take into account the wind. In this paper we shall study six possible modifications of the source-mesh, in order to reduce the present drawbacks. the first four will consist in reducing the number of points that modelize the source so that computing by MEDIA (which takes account of the meteorological fields) will step in as soon as possible. In the last two ones we shall take into account the wind in order to reduce the concentration of pollutant upstream from the source, in a way which corresponds as much as possible.to the reality. 1-INTRODUCTION: MEDIA is an Eulerian three dimensional model of atmospheric dispersion, which was developped in METEO-FRANCE untill 1990 (JP Piedelievre and al,1990) to improve the forecasting of the behaviour of plums, particularly in case of atmospheric release of pollutant. Actually the Trajectory model used in the French Weather Service until 1990 was quite efficient to locate the source of the accident, but drawing up the map of the concentrations needed several trajectories and so a long time of computation even when a rapid decision was necessary. MEDIA uses fields (of wind, temperature...) computed by an operational model : it takes advantage of improvements in numerical forecasting. The source is described by a Gaussian distribution on the first 12 points around the release, and does not take the wind or any other meteorological field into account : points located upstream from the wind are contaminated the same way as if they were downstream. In these conditions, MEDIA was tested during the ATMES experiment, (W Klug et al, 1992) and gave good results. Nevertheless it was interesting to see if a change in the source mesh numerical treatment could lead to actual better results. We decided to value the modifications with the hypothesis of ATMES experiment for the source release. In this paper we shall first recall the basic equations and hypothesis of the MEDIA model. In a second part we shall present the different experiments : - in the four first ones, we shall reduce the number of points of the source so that the wind is taken into account by MEDIA as soon as possible.
2 114 Computer Simulation - in the two last experiments we shall use the wind field since the first points to reduce the concentration of pollutant upstream from the source. We shall then present the results from these experiments in the third part. 2 EQUATIONS REMINDER: The concentration of a passive chemical element in the atmosphere is described by the laws of a continuous medium, especially mass conservation : C is the pollutant concentration at a given node in time t ; V is the wind vector predicted by a numerical weather model ; K the REYNOLDS tensor ; S, the source terms ; Si the sink term. The pollutant concentration is advected by the wind field as predicted by the operational coupling model. In order to simplify the procedure and to be coherent with the coupling model, the diffusion is modeled using exchange coefficients. This choice means that the pollutant is assumed to be diffused in the same way as the water vapor. Source term description : In the initial code, the concentration of each of the 12 points located around the source as described on the following figure is given by a Gaussian distribution : r -</ Where: d is the distance from the source; Q the pollutant flux (UCs-0; H the vertical thickness of the pollutant cloud (m), and <jn the surface of the mesh which includes the source This method does not include the effect of the wind close to the source ; on the other hand, the loss of mass is quite small : the sum of the 12 values of concentration reaches 96 % of the emission. Coordinates system : We chose the a coordinate (a- P I P^) on a latitude longitude grid in accordance with the evolution of the data bases of the numerical weather forecasting models we are using (Model from the European Center for Medium Weather Forecasts -ECMWF -in that case).
3 The equation for the concentration evolution can then be written : Computer Simulation 115 i a Where: (/) is the latitude; A the longitde; R the earth's radius; TV the virtual temperature in Kelvin; a = da/dt and Ra = the dry air constant. 3 MODIFICATIONS OF THE SOURCE TREATMENT STUDIED : At the present time the model does not include the effect of the wind for the first 12 points around the source and nodes located upstream from the source are contaminated the same way as if they were downstream. The first 4 experiments will consist in reducing this number of points so that it will be taken account of the wind, by MEDIA as soon as possible. A - First experiment: The Gaussian distribution is only applied to the first four grid points around the source (c.f the following diagram); this method induces a 32 % loss of mass which is corrected by the addition of the following term : 0.32 x Q x CTn where : Q is the pollutant flux (UCs-i); H is the vertical thickness of the pollutant cloud (m), and On is the surface of the mesh which includes the source..ai is a value of surface as described by the following diagram. a-4 c/3 So the concentration on point i gets :
4 116 Computer Simulation where :d is the distance from the source; B - Second experiment : We consider that the whole mass of pollutant is concentrated on the grid point which is the closest from the source. C - Third experiment : We consider that every of the first four nodes around the source receives the concentration : c/c D - Fourth experiment : The first mesh around the source is devided in 4 parts as shown by the figure. The concentration on every point is then computed as following : dt exp... _ +0,32 x()x 2cTn} ~ otr where : a'l a'3 In the last two experiments the wind is taken into account in order to reduce the concentration of pollutant downstream of the source, in a way which corresponds, as much as possible to the reality.
5 E - Fifth experiment : Computer Simulation 117 We consider the 12 points around the source, the concentration of which is first calculated by a Gaussian distribution as in the initial version of the model. We then apply a symmetry about the plan which is orthogonal to the wind and includes the source : the upstream from the source has its concentration transfered on the downstream part of the field. This quantity of pollutant is there distributed among the different points according to their distance from the source. F - Sixth experiment : We consider the grid point which is the closest from the source.its coordinates are (Is,Js). We then discretize each cell of the field (Is-2,Is+2) ; (Js-2,Js+2) in 8 elements on each axis, which leads to a (33,33,15) networking (15 is the number of levels on the vertical axis). The concentration of all the points of the domain is computed with MEDIA and the assumption that the whole mass is located on the grid point which is the nearest to the source, in the new network (this point can be different from the one we called (Is,Js) at the beginning ). The fields used at this level are the outputs of the meteorological model, interpolated on the (33,33,15) domain. The concentration of the 9 points belonging to the area (is-l,is+l) /(Js-l,Js+l) is then computed as the average value of the points located in the four adjacent meshes. MEDIA is at that time started from these 9 points of grid, with the usual mesh size.
6 118 Computer Simulation 4 RESULTS We compared daily air concentrations predicted by MEDIA with measurements collected from the REM data bank (F. Raes et al, 1990). Participants in the final ATMES workshop were provided with these data. For each of the six different simulations, we computed : Mean value and variance of the sample Bias defined as B = ]T (C(i)-O(i)) where C(i) is the predicted value and 0(i) the observed value * Normalized Mean Square Error, defined as NMSE = V #Y co where C is the predicted mean value and 0 the observed mean value. Factors FA2 ans FAS, which represent the percentage of predicted values which are within factor 2 or 5 of the observed value. The measurements sample can be described by its mean value 0 = 0,664 Bq/nf and its variance 0,733. Table 1 shows the values of the parameters for each simulation and for the ATMES sample (MEDIA results during the ATMES experiment). ATMES Expl Exp2 Exp3 Exp4 ExpS Exp6 Mean value 0,423 0,429 0,326 0,313 0,317 0,430 0,579 Variance 0,870 0,747 0,788 0,520 0,625 1,086 1,026 Bias -0,24-0,23-0,33-0,35-0,35-0,23-0,085 NMSE 2,99 2,55 4,33 3,92 4,02 3,57 2,43 FA2 26% 27% 19% 21% 22% 25% 36% FAS 57% 60% 48% 48% 48% 57% 65% Table 1 : statistical coefficients for each simulation (Total sample : 156 values) It appears that we can gather the results of the different simulations in 3 classes. Class A (experiment 2,3 and 4) represents the simplest simulations of the source term. Class B stands for the gaussian (function of the distance) representation of the source term (experiments 1, 5 and ATMES). Class C includes a complex simulation of the source term, taking into account the meteorological conditions within the source mesh, by coupling a subgrid eulerian model to MEDIA The mean values are about 0,3 Bq/nf for class A, 0,43 Bq/m^ for class B and nearly 0,6 Bq/nf for class C, while the bias are near -0,3, -0,2 and -0,1 for class A, B, C. Coefficient FA2 (and FAS) shows that 20% (48%) of the predicted sample is within factor 2 (5) of the observed sample for class A, 26% (57% to 60%) for class B and 36% (65%) for experiment 6.
7 Computer Simulation 119 As a conclusion, the overall behaviour of experiments 2,3 and 4 is not as good as the atmes simulation. Experient 5 gives results similar to those obtained with the ATMES simulation. Experiment 1 is slightly better than atmes simulation and experiment 6 seems to provide a significant improvement of the results. We then selected 5 locations : Praha, Attikis, Berlin, Innsbruck and Helsinki and compared the predited values to the observed ones. We define the total coverage TC as TC = rap(i) = O if not. with rap(i) = if C(i) is lower than 0(i) or Values of TC (in percentage) for the 5 selected locations and for the whole sample are presented in table 2. ATMES Expl Exp2 Exp3 Exp4 Exp5 Exp6 Attikis Berlin Helsinki Innsbruck Praha Total sample Attikis (see figure 1) : MEDIA overestimates air concentrations for all the simulations. Experiment 3 gives the best results with a total coverage of 63%. The differences between all the simulations are not very important, almost all the predicted values are within factor 5 of the observed data Experiments 3, 1, 4 and 6 improve the results of the ATMES simulation.
8 120 Computer Simulation ATTIKIS 06/5 07/5 08/5 09/5 Figure 1 : Time evolution of daily average air concentration. Comparison between measurements and the 6 simulations from to Berlin (see figure 2) : MEDIA underestimates air concentrations for all the simulations. Experiment 5 gives the best results with a coverage of 46%. Except for experiment 2, there is no great difference between the simulations. Experiment 1 improves slightly the result of the ATMES simulation. BERLIN fli D ATMES expl exp2 D exp3 D exp4 B exps (II expg 29/4 30/4 01/5 02/5 03/5 04/5 05/5 06/5 07/5 08/5 Figure 2 : Time evolution of daily average air concentration. Comparison between measurements and the 6 simulations from to
9 Helsinki (see figure 3) : Computer Simulation 121 MEDIA slightly overestimates air concentrations for experiments L, 3,4 and 5. Except for experiment 6, the simulations are very close to each other. No experiment improves the results of the ATMES simulation. HELSINKI E ff 29/4 30/4 01/5 02/5 03/5 04/5 05/5 06/5 07/5 08/5 09/5 Figure 3 : Time evolution of daily average air concentration. Comparison between measurements and the 6 simulations from to Innsbruck (see figure 4) : MEDIA generally overestimates air concentrations by a slight amount for all the simulations. Experiment 6 gives the best results. Differences between experiments are significant over this location. Experiment 1 gives quite the same results as the ATMES simulation.
10 122 Computer Simulation INNSBRUCK i -' DATMES I _ expl i exp2! Oexp3, I ; U exp4 S exp5 ; dexp6 I 02/5 03/5 04/5 0*5/5 06/5 07/5 LJL Figure 4 : Time evolution of daily average air concentration. Comparison between measurements and the 6 simulations from to Praha (see figure 5) : Except for experiment 6, MEDIA underestimates air concentrations. The difference between experiments are not vey important. Experiment 1, 5 and 6 gives similar results to those derived from the ATMES simulation. PRAHA LJ ATNtS exp! exp2 D exp3 Dexp4 Sexp5! Ml exp6! 30/4 01/5 02/5 03/5 04/5 05/5 06/5 07/5 08/5 09/5 Figure 5 : Time evolution of daily average air concentration. Comparison between measurements and the 6 simulations from to
11 Computer Simulation 123 Figure 6 shows a chart of air concentration contour level above 0.1 Bequerel per cubic meter for the different simulations. The contours are very similar. The shapes of the polluted cloud look almost identical at long range. The differences appear near the source and on the eastern part of the clouds. Figure 6 : Air concentration on at 00 UTC (Significant value : 0.1 Bq/m3) 5 CONCLUSION AND FUTURE The comparison presented in the paper stresses that the simplest representations of the source term (putting the whole mass at the nearest grid point, contamining the 4 grid points surrounding the location of the source with quarter of the mass or proportionnaly with the distance) do not improve the results of the ATMES simulation. The Gaussian distribution, as a function of the distance from the source, used during the atmes experiment gave almost the same results when reduccing the number of polluted points (experiment 1) or when trying to only contamine points located downstream from the source (experiment 5). Nevertheless experiment 1 slightly improves the overall estimation of air concentrations. The spread of the predicted sample is similar to the spread of the observed one, the mean value is a little closer to the observed one, the bias is smaller and FA2, FA5 are higher.
12 124 Computer Simulation The complex simulation of the source term (MEDIA within MEDIA) gives a significant improvement of the results for all the criteria studied except for the spread, which is higher than the observed one. We are now trying to get more measured data in order to increase the sample with locations near the source (especially RIGA, VILNIUS, BARYSHEVKA, RAKHOV, MINSK and ST PETERSBURG). It should be emphasized that all the measurements that we get, are located far from Chernobyl; the nearest one is about 1200 km. At this distance, the differences between all simulations are not very important. But discrepancies could be worse at shorter range. On another hand, a comparison over one case does not lead to a definite conclusion, it only helps us in our way to improve the operational code MEDIA by giving some informations over different simulations of the source term. On major point is that we succeed in coupling MEDIA (large scale) with a subgrid model. In the future we will develop a lagrangian model for the source mesh and compare MEDIA(Langrangian source mesh) with the observed sample. The computing time was quite the same for experiments 1, 2, 3, 4 and for ATMES simulation. Experiment 5 required 5% more time and experiment 6, 35% more time. REFERENCES : Klug, W,G.Graziani, G.Grippa, D.Pierce & C.Tassone. Evaluation Of Long Range Atmospheric, transport Models Using Environmental Radioactivity Data From The Chernobyl Accident - The ATMES Report, Elsevier Applied Science, London and New York. PiedelievreJ.P, L Musson-Genon & F.Bompay. MEDIA An Euleran Model of Atmospheric dispersin:first Validation on the Chernobyl Relcasejournal of Applied Meteorology,Vol 29, N 12, December 1990,American Meteorological Society. Raes F,G Graziani, L Grossi, L.Marciano,D. Peirs, B.Pedersen,D. Stanners & N.Zarimpas. Radioactivity Measurements in Europe after the Chernobyl Accident,Commission of the European Communities, Nuclear Science and Technology, Directorate-General for Science, research and Development-joint research centre-ispra Site
RSMC MONTRÉAL USERS' INTERPRETATION GUIDELINES ATMOSPHERIC TRANSPORT MODEL OUTPUTS. Version 9
RSMC MONTRÉAL USERS' INTERPRETATION GUIDELINES ATMOSPHERIC TRANSPORT MODEL OUTPUTS Version 9 Environmental Emergency Response Section RSMC Montréal Canadian Meteorological Centre Meteorological Service
More information2. The french National Centre for Meteorological Forecasts (CNP)
RSMC TOULOUSE USERS INTERPRETATION GUIDELINES ATMOSPHERIC TRANSPORT MODEL OUTPUTS March 2016 1. Introduction Météo-France is a French executive agency with legal status and financial autonomy. Its responsibilities
More informationRSMC WASHINGTON USER'S INTERPRETATION GUIDELINES ATMOSPHERIC TRANSPORT MODEL OUTPUTS
RSMC WASHINGTON USER'S INTERPRETATION GUIDELINES ATMOSPHERIC TRANSPORT MODEL OUTPUTS -Version 2.0- (January 2007) 1. Introduction In the context of current agreements between the National Oceanic and Atmospheric
More informationQUALITY CONTROL OF WINDS FROM METEOSAT 8 AT METEO FRANCE : SOME RESULTS
QUALITY CONTROL OF WINDS FROM METEOSAT 8 AT METEO FRANCE : SOME RESULTS Christophe Payan Météo France, Centre National de Recherches Météorologiques, Toulouse, France Astract The quality of a 30-days sample
More information1.21 SENSITIVITY OF LONG-TERM CTM SIMULATIONS TO METEOROLOGICAL INPUT
1.21 SENSITIVITY OF LONG-TERM CTM SIMULATIONS TO METEOROLOGICAL INPUT Enrico Minguzzi 1 Marco Bedogni 2, Claudio Carnevale 3, and Guido Pirovano 4 1 Hydrometeorological Service of Emilia Romagna (SIM),
More informationApplication and verification of ECMWF products 2015
Application and verification of ECMWF products 2015 Hungarian Meteorological Service 1. Summary of major highlights The objective verification of ECMWF forecasts have been continued on all the time ranges
More informationCFD calculations of the test 2-4 experiments. Author: G. de With
CFD calculations of the test 2-4 experiments Author: G. de With 34. Model setup and boundary conditions Dimensions CFD model: x=1000m / y=100m / z=2000m. CFD Model: Transient simulation, with steady-state
More information1 INTRODUCTION 2 DESCRIPTION OF THE MODELS. In 1989, two models were able to make smog forecasts; the MPA-model and
The national smog warning system in The Netherlands; a combination of measuring and modelling H. Noordijk Laboratory of Air Research, National Institute of Public Health and Environmental Protection (WFM;,
More informationApplication and verification of ECMWF products 2009
Application and verification of ECMWF products 2009 Hungarian Meteorological Service 1. Summary of major highlights The objective verification of ECMWF forecasts have been continued on all the time ranges
More informationLessons Learned in Responding to the Accident at the Fukushima Daiichi Nuclear Power Stations: The Way Forward. World Meteorological Organization
Lessons Learned in Responding to the Accident at the Fukushima Daiichi Nuclear Power Stations: The Way Forward Geoff Love Director, Weather and Disaster Risk Reduction Services World Meteorological Organization
More informationApplication and verification of the ECMWF products Report 2007
Application and verification of the ECMWF products Report 2007 National Meteorological Administration Romania 1. Summary of major highlights The medium range forecast activity within the National Meteorological
More informationRegional Production Quarterly report on the daily analyses and forecasts activities, and verification of the CHIMERE performances
Regional Production Quarterly report on the daily analyses and forecasts activities, and verification of the CHIMERE performances December 2015 January 2016 February 2016 This document has been produced
More informationEvaluation of the operational IDC ATM products by comparing HYSPLIT BA, FA and OMEGA FOR during Level 4 events in Scandinavia
Evaluation of the operational IDC ATM products by comparing HYSPLIT BA, FA and OMEGA FOR during Level 4 events in Scandinavia Gerhard Wotawa, Martin B. Kalinowski Comprehensive Nuclear-Test-Ban Treaty
More informationModeling Radiological Consequences of Sever Accidents in BWRs: Review of Models Development, Verification and Validation
Modeling Radiological Consequences of Sever Accidents in BWRs: Review of Models Development, Verification and Validation Mohamed GAHEEN Department of Research on Nuclear Safety and Radiological Emergencies,
More informationApplication of a Three-Dimensional Prognostic Model During the ETEX Real-Time Modeling Exercise: Evaluatin of Results (u)
WSRC-MS-96-0766 COdF- 9 7 0 9 1 9 m - 9 Application of a Three-Dimensional Prognostic Model During the ETEX Real-Time Modeling Exercise: Evaluatin of Results (u) by D. P. Griggs Westinghouse Savannah River
More informationLessons learned in responding to the accident at the Fukushima Daiichi nuclear power stations: THE WAY FORWARD
WMO World Meteorological Organization Working together in weather, climate and water Lessons learned in responding to the accident at the Fukushima Daiichi nuclear power stations: THE WAY FORWARD WMO:
More informationValidation of 2-meters temperature forecast at cold observed conditions by different NWP models
Validation of 2-meters temperature forecast at cold observed conditions by different NWP models Evgeny Atlaskin Finnish Meteorological Institute / Russian State Hydrometeorological University OUTLINE Background
More informationAssimilating AMSU-A over Sea Ice in HIRLAM 3D-Var
Abstract Assimilating AMSU-A over Sea Ice in HIRLAM 3D-Var Vibeke W. Thyness 1, Leif Toudal Pedersen 2, Harald Schyberg 1, Frank T. Tveter 1 1 Norwegian Meteorological Institute (met.no) Box 43 Blindern,
More informationA VALIDATION EXERCISE ON THE SAFE-AIR VIEW SOFTWARE. Joint Research Centre NDFM Ispra, Italy 2
A VALIDATION EXERCISE ON THE SAFE-AIR VIEW SOFTWARE F. D Alberti 1, F. d Amati 1, E. Canepa 2, G. Triacchini 3 1 Joint Research Centre NDFM Ispra, Italy 2 CNR INFM CNISM Department of Physics, University
More informationOn the impact of the assimilation of ASAR wave spectra in the wave model MFWAM
On the impact of the assimilation of ASAR wave spectra in the wave model MFWAM Lotfi Aouf, Jean-Michel Lefèvre Météo-France, Toulouse SEASAR 2012, 4 th International workshop on Advances in SAR oceanography,
More informationErik L. Swanberg 1 and Steven G. Hoffert 2. Veridian Systems 1, Autometric 2. Sponsored by Defense Threat Reduction Agency
USING ATMOSPHERIC 137 CS MEASUREMENTS AND HYSPLIT TO CONFIRM CHERNOBYL AS A SOURCE OF 137 CS IN EUROPE Erik L. Swanberg 1 and Steven G. Hoffert 2 Veridian Systems 1, Autometric 2 Sponsored by Defense Threat
More informationAssimilation Experiments of One-dimensional Variational Analyses with GPS/MET Refractivity
Assimilation Experiments of One-dimensional Variational Analyses with GPS/MET Refractivity Paul Poli 1,3 and Joanna Joiner 2 1 Joint Center for Earth Systems Technology (JCET), University of Maryland Baltimore
More informationRegional Production, Quarterly report on the daily analyses and forecasts activities, and verification of the CHIMERE performances
ECMWF COPERNICUS REPORT Copernicus Atmosphere Monitoring Service Regional Production, Quarterly report on the daily analyses and forecasts activities, and verification of the CHIMERE performances September
More informationMasahiro Kazumori, Takashi Kadowaki Numerical Prediction Division Japan Meteorological Agency
Development of an all-sky assimilation of microwave imager and sounder radiances for the Japan Meteorological Agency global numerical weather prediction system Masahiro Kazumori, Takashi Kadowaki Numerical
More informationApplication and verification of ECMWF products: 2010
Application and verification of ECMWF products: 2010 Hellenic National Meteorological Service (HNMS) F. Gofa, D. Tzeferi and T. Charantonis 1. Summary of major highlights In order to determine the quality
More information17th International Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes 9-12 May 2016, Budapest, Hungary
17th International Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes 9-12 May 2016, Budapest, Hungary USING METEOROLOGICAL ENSEMBLES FOR ATMOSPHERIC DISPERSION
More informationCurrent Limited Area Applications
Current Limited Area Applications Nils Gustafsson SMHI Norrköping, Sweden nils.gustafsson@smhi.se Outline of talk (contributions from many HIRLAM staff members) Specific problems of Limited Area Model
More informationApplication and verification of ECMWF products 2009
Application and verification of ECMWF products 2009 Danish Meteorological Institute Author: Søren E. Olufsen, Deputy Director of Forecasting Services Department and Erik Hansen, forecaster M.Sc. 1. Summary
More informationApplication and verification of ECMWF products 2013
Application and verification of EMWF products 2013 Hellenic National Meteorological Service (HNMS) Flora Gofa and Theodora Tzeferi 1. Summary of major highlights In order to determine the quality of the
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 informationApplication and verification of ECMWF products 2017
Application and verification of ECMWF products 2017 Finnish Meteorological Institute compiled by Weather and Safety Centre with help of several experts 1. Summary of major highlights FMI s forecasts are
More informationUsing the Mercator ocean forecasting system to compute coastal maritime pollution risk indicators on the Atlantic European coasts
Environmental Problems in Coastal Regions VI 437 Using the Mercator ocean forecasting system to compute coastal maritime pollution risk indicators on the Atlantic European coasts S. Besnard 1, E. Dombrowsky
More informationA Global Atmospheric Model. Joe Tribbia NCAR Turbulence Summer School July 2008
A Global Atmospheric Model Joe Tribbia NCAR Turbulence Summer School July 2008 Outline Broad overview of what is in a global climate/weather model of the atmosphere Spectral dynamical core Some results-climate
More informationApplication and verification of ECMWF products 2010
Application and verification of ECMWF products 2010 METEO-FRANCE Bruno Lacroix 1. Summary of major highlights Implementation of a limited area model ALADIN configuration coupled with IFS over at 00 and
More informationApplication and verification of ECMWF products 2015
Application and verification of ECMWF products 2015 METEO- J. Stein, L. Aouf, N. Girardot, S. Guidotti, O. Mestre, M. Plu, F. Pouponneau and I. Sanchez 1. Summary of major highlights The major event is
More informationRegional Production, Quarterly report on the daily analyses and forecasts activities, and verification of the CHIMERE performances
ECMWF COPERNICUS REPORT Copernicus Atmosphere Monitoring Service Regional Production, Quarterly report on the daily analyses and forecasts activities, and verification of the CHIMERE performances June
More informationImpact of sea surface temperature on COSMO forecasts of a Medicane over the western Mediterranean Sea
Impact of sea surface temperature on COSMO forecasts of a Medicane over the western Mediterranean Sea V. Romaniello (1), P. Oddo (1), M. Tonani (1), L. Torrisi (2) and N. Pinardi (3) (1) National Institute
More informationRSMC WASHINGTON USER'S INTERPRETATION GUIDELINES ATMOSPHERIC TRANSPORT MODEL OUTPUTS
RSMC WASHINGTON USER'S INTERPRETATION GUIDELINES ATMOSPHERIC TRANSPORT MODEL OUTPUTS (January 2015) 1. Introduction In the context of current agreements between the National Oceanic and Atmospheric Administration
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 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 informationOSI SAF Sea Ice products
OSI SAF Sea Ice products Lars-Anders Brevik, Gorm Dybkjær, Steinar Eastwood, Øystein Godøy, Mari Anne Killie, Thomas Lavergne, Rasmus Tonboe, Signe Aaboe Norwegian Meteorological Institute Danish Meteorological
More information2 1-D Experiments, near surface output
Tuning CBR A.B.C. Tijm 1 Introduction The pre 6.2 reference Hirlam versions have a positive bias in the wind direction under stable conditions (night and winter conditions), accompanied by too strong near
More informationHYSPLIT RESULTS: HYSPLIT:-
HYSPLIT RESULTS: HYSPLIT:- Lagrangian particle models compute trajectories of a large number of so-called particles to describe the transport and diffusion of tracers in the atmosphere. The main advantage
More informationSwedish Meteorological and Hydrological Institute
Swedish Meteorological and Hydrological Institute Norrköping, Sweden 1. Summary of highlights HIRLAM at SMHI is run on a CRAY T3E with 272 PEs at the National Supercomputer Centre (NSC) organised together
More informationECMWF global reanalyses: Resources for the wind energy community
ECMWF global reanalyses: Resources for the wind energy community (and a few myth-busters) Paul Poli European Centre for Medium-range Weather Forecasts (ECMWF) Shinfield Park, RG2 9AX, Reading, UK paul.poli
More informationA study on the spread/error relationship of the COSMO-LEPS ensemble
4 Predictability and Ensemble Methods 110 A study on the spread/error relationship of the COSMO-LEPS ensemble M. Salmi, C. Marsigli, A. Montani, T. Paccagnella ARPA-SIMC, HydroMeteoClimate Service of Emilia-Romagna,
More informationIntroduction to atmospheric dispersion modelling. M.Sofiev
Introduction to atmospheric dispersion modelling M.Sofiev Content Main parts of the atmosphere Basic terms atmospheric tracer temporal and spatial scales life time in the atmosphere life cycle of atmospheric
More informationEvaluation of the Global Atmospheric Multi-Layer Transport (GAMUT) Model using the European Tracer Experiment (ETEX)
Poster 451 Evaluation of the Global Atmospheric Multi-Layer Transport (GAMUT) Model using the European Tracer Experiment (ETEX) Stephen E. Masters* and Mark A. Kienzle, ENSCO, Inc. 1. Introduction The
More informationNicolas Duchene 1, James Smith 1 and Ian Fuller 2
A METHODOLOGY FOR THE CREATION OF METEOROLOGICAL DATASETS FOR LOCAL AIR QUALITY MODELLING AT AIRPORTS Nicolas Duchene 1, James Smith 1 and Ian Fuller 2 1 ENVISA, Paris, France 2 EUROCONTROL Experimental
More informationWind and turbulence experience strong gradients in vegetation. How do we deal with this? We have to predict wind and turbulence profiles through the
1 2 Wind and turbulence experience strong gradients in vegetation. How do we deal with this? We have to predict wind and turbulence profiles through the canopy. 3 Next we discuss turbulence in the canopy.
More informationKnow and Respond AQ Alert Service. Paul Willis SCOTTISH AIR QUALITY DATABASE AND WEBSITE ANNUAL SEMINAR Stirling 30 th March 2011
Know and Respond AQ Alert Service Paul Willis SCOTTISH AIR QUALITY DATABASE AND WEBSITE ANNUAL SEMINAR Stirling 30 th March 2011 Outline: Know and Respond AQ Alert Service Introduction. The Website. How
More informationCredibility of climate predictions revisited
European Geosciences Union General Assembly 29 Vienna, Austria, 19 24 April 29 Session CL54/NP4.5 Climate time series analysis: Novel tools and their application Credibility of climate predictions revisited
More information350 Int. J. Environment and Pollution Vol. 5, Nos. 3 6, 1995
350 Int. J. Environment and Pollution Vol. 5, Nos. 3 6, 1995 A puff-particle dispersion model P. de Haan and M. W. Rotach Swiss Federal Institute of Technology, GGIETH, Winterthurerstrasse 190, 8057 Zürich,
More information4.13 EVALUATION OF NEAR FIELD DISPERSION LAGRANGIEN MODELLING WITH KRYPTON 85 MEASUREMENTS AROUND COGEMA LA HAGUE NUCLEAR REPROCESSING PLANT (FRANCE)
4.13 EVALUATION OF NEAR FIELD DISPERSION LAGRANGIEN MODELLING WITH KRYPTON 85 MEASUREMENTS AROUND COGEMA LA HAGUE NUCLEAR REPROCESSING PLANT (FRANCE) Christine Lac 1, Denis Maro 2, Didier Hebert 2, François
More information2D.4 THE STRUCTURE AND SENSITIVITY OF SINGULAR VECTORS ASSOCIATED WITH EXTRATROPICAL TRANSITION OF TROPICAL CYCLONES
2D.4 THE STRUCTURE AND SENSITIVITY OF SINGULAR VECTORS ASSOCIATED WITH EXTRATROPICAL TRANSITION OF TROPICAL CYCLONES Simon T. Lang Karlsruhe Institute of Technology. INTRODUCTION During the extratropical
More informationCommissariat à l Energie Atomique, Département Analyse, Surveillance, Environnement Bruyères-le-Châtel, France 2
ANALYSIS OF ATMOSPHERIC RADIOXENON ACTIVITIES MEASURED BY A RADIONUCLIDE GAS STATION LOCATED IN FRANCE: SIMULATION OF THE ATMOSPHERIC TRANSPORT WITH A MESOSCALE MODELLING SYSTEM Patrick Armand 1, Pascal
More informationBias correction of satellite data at Météo-France
Bias correction of satellite data at Météo-France É. Gérard, F. Rabier, D. Lacroix, P. Moll, T. Montmerle, P. Poli CNRM/GMAP 42 Avenue Coriolis, 31057 Toulouse, France 1. Introduction Bias correction at
More informationRegional Production, Quarterly report on the daily analyses and forecasts activities, and verification of the MATCH performances
ECMWF COPERNICUS REPORT Copernicus Atmosphere Monitoring Service Regional Production, Quarterly report on the daily analyses and forecasts activities, and verification of the MATCH performances September
More informationSupplement to the. Final Report on the Project TRACHT-MODEL. Transport, Chemistry and Distribution of Trace Gases in the Tropopause Region: Model
Anhang 2 Supplement to the Final Report on the Project TRACHT-MODEL Transport, Chemistry and Distribution of Trace Gases in the Tropopause Region: Model H. Feldmann, A. Ebel, Rheinisches Institut für Umweltforschung
More informationThe WMO Aviation Research & Demonstration Project (AvRDP) at Paris-CDG airport. Pauline Jaunet Météo-France
The WMO Aviation Research & Demonstration Project (AvRDP) at Paris-CDG airport Pauline Jaunet Météo-France AvRDP focus points Aviation Research Demonstration Project Joint effort between CAS/WWRP and CaeM
More informationNUCLEAR EMERGENCY RESPONSE ACTIVITIES COORDINATION GROUP Original: ENGLISH EC-JRC ENSEMBLE APPROACH TO LONG RANGE ATMOSPHERIC DISPERSION
WORLD METEOROLOGICAL ORGANIZATION COMMISSION FOR BASIC SYSTEMS CBS/CG-NERA/Doc. (number) (26.04.2006) OPAG on DPFS NUCLEAR EMERGENCY RESPONSE ACTIVITIES COORDINATION GROUP Original: ENGLISH VIENNA, AUSTRIA,
More informationDiagnostics of the prediction and maintenance of Euro-Atlantic blocking
Diagnostics of the prediction and maintenance of Euro-Atlantic blocking Mark Rodwell, Laura Ferranti, Linus Magnusson Workshop on Atmospheric Blocking 6-8 April 2016, University of Reading European Centre
More informationNaka-Gun, Ibaraki, , Japan
Examination of Atmospheric Dispersion Model s Performance - Comparison with the Monitoring Data under the Normal Operation of the Tokai Reprocessing Plant - M. Takeyasu 1, M. Nakano 1, N. Miyagawa 1, M.
More informationPortugal. Instituto de Meteorologia
WWW TECHNICAL PROGRESS REPORT ON THE GLOBAL DATA- PROCESSING AND FORECASTING SYSTEM (GDPFS), AND THE ANNUAL NUMERICAL WEATHER PREDICTION (NWP) PROGRESS REPORT FOR THE YEAR 2005 Portugal Instituto de Meteorologia
More informationApplication and verification of ECMWF products 2010
Application and verification of ECMWF products Hydrological and meteorological service of Croatia (DHMZ) Lovro Kalin. Summary of major highlights At DHMZ, ECMWF products are regarded as the major source
More informationP1.6 Simulation of the impact of new aircraft and satellite-based ocean surface wind measurements on H*Wind analyses
P1.6 Simulation of the impact of new aircraft and satellite-based ocean surface wind measurements on H*Wind analyses Timothy L. Miller 1, R. Atlas 2, P. G. Black 3, J. L. Case 4, S. S. Chen 5, R. E. Hood
More informationTHREE-DIMENSIONAL FINITE DIFFERENCE MODEL FOR TRANSPORT OF CONSERVATIVE POLLUTANTS
Pergamon Ocean Engng, Vol. 25, No. 6, pp. 425 442, 1998 1998 Elsevier Science Ltd. All rights reserved Printed in Great Britain 0029 8018/98 $19.00 + 0.00 PII: S0029 8018(97)00008 5 THREE-DIMENSIONAL FINITE
More informationFollow this and additional works at:
International Congress on Environmental Modelling and Software Brigham Young University BYU ScholarsArchive 6th International Congress on Environmental Modelling and Software - Leipzig, Germany - July
More informationA mechanistic model for seed dispersal
A mechanistic model for seed dispersal Tom Robbins Department of Mathematics University of Utah and Centre for Mathematical Biology, Department of Mathematics University of Alberta A mechanistic model
More informationApplication and verification of ECMWF products 2016
Application and verification of ECMWF products 2016 Hellenic National Meteorological Service (HNMS) Flora Gofa and Panagiotis Skrimizeas 1. Summary of major highlights In order to determine the quality
More informationApplication and verification of ECMWF products at the Finnish Meteorological Institute
Application and verification of ECMWF products 2010 2011 at the Finnish Meteorological Institute by Juhana Hyrkkènen, Ari-Juhani Punkka, Henri Nyman and Janne Kauhanen 1. Summary of major highlights ECMWF
More informationMETEOSAT cloud-cleared radiances for use in three/fourdimensional variational data assimilation
METEOSAT cloud-cleared radiances for use in three/fourdimensional variational data assimilation G. A. Kelly, M. Tomassini and M. Matricardi European Centre for Medium-Range Weather Forecasts, Reading,
More informationDeutscher Wetterdienst
Deutscher Wetterdienst Modelling the Volcanic Ash Episode: Experiences with COSMO-ART Detlev Majewski (FE1) Bernhard Vogel, Heike Vogel (KIT) Thomas Hanisch, Jochen Förstner (FE13), Ulrich Pflüger (FE15)
More informationClustering Forecast System for Southern Africa SWFDP. Stephanie Landman Susanna Hopsch RES-PST-SASAS2014-LAN
Clustering Forecast System for Southern Africa SWFDP Stephanie Landman Susanna Hopsch Introduction The southern Africa SWFDP is reliant on objective forecast data for days 1 to 5 for issuing guidance maps.
More information11/8/2018. Spatial Interpolation & Geostatistics. Kriging Step 1
(Z i Z j ) 2 / 2 (Z i Zj) 2 / 2 Semivariance y 11/8/2018 Spatial Interpolation & Geostatistics Kriging Step 1 Describe spatial variation with Semivariogram Lag Distance between pairs of points Lag Mean
More informationReynolds Averaging. Let u and v be two flow variables (which might or might not be velocity components), and suppose that. u t + x uv ( ) = S u,
! Revised January 23, 208 7:7 PM! Reynolds Averaging David Randall Introduction It is neither feasible nor desirable to consider in detail all of the small-scale fluctuations that occur in the atmosphere.
More informationApplication and verification of ECMWF products 2016
Application and verification of ECMWF products 2016 RHMS of Serbia 1 Summary of major highlights ECMWF forecast products became the backbone in operational work during last several years. Starting from
More informationNumerical Modelling for Optimization of Wind Farm Turbine Performance
Numerical Modelling for Optimization of Wind Farm Turbine Performance M. O. Mughal, M.Lynch, F.Yu, B. McGann, F. Jeanneret & J.Sutton Curtin University, Perth, Western Australia 19/05/2015 COOPERATIVE
More information6 th INTERNATIONAL WORKSHOP ON SAND/DUSTSTORMS AND ASSOCIATED DUSTFALL 7-9 September 2011, Athens, Greece
6 th INTERNATIONAL WORKSHOP ON SAND/DUSTSTORMS AND ASSOCIATED DUSTFALL Motivations Importance of Numerical Prediction Models to mineral dust cycle evaluation of dust effects over Italian region Identify
More informationREPORT ON THE ACTIVITIES OF THE EUMETSAT-ESSL CONVECTION WORKING GROUP
REPORT ON THE ACTIVITIES OF THE EUMETSAT-ESSL CONVECTION WORKING GROUP Marianne König EUMETSAT, Eumetsat Allee 1, 64295 Darmstadt, Germany Abstract The focus of the Convection Working Group is to have
More informationNerushev A.F., Barkhatov A.E. Research and Production Association "Typhoon" 4 Pobedy Street, , Obninsk, Kaluga Region, Russia.
DETERMINATION OF ATMOSPHERIC CHARACTERISTICS IN THE ZONE OF ACTION OF EXTRA-TROPICAL CYCLONE XYNTHIA (FEBRUARY 2010) INFERRED FROM SATELLITE MEASUREMENT DATA Nerushev A.F., Barkhatov A.E. Research and
More informationWINDS FROM METEOSAT AND GMS WV-IMAGERY USING VARIATIONAL TECHNIQUES. M. Rohn. Institut für Meteorologie, Freie Universität Berlin Berlin, F.R.G.
WINDS FROM METEOSAT AND GMS WV-IMAGERY USING VARIATIONAL TECHNIQUES M. Rohn Institut für Meteorologie, Freie Universität Berlin 14195 Berlin, F.R.G. ABSTRACT A new method is described where motion analysis
More informationStability in SeaWinds Quality Control
Ocean and Sea Ice SAF Technical Note Stability in SeaWinds Quality Control Anton Verhoef, Marcos Portabella and Ad Stoffelen Version 1.0 April 2008 DOCUMENTATION CHANGE RECORD Reference: Issue / Revision:
More informationRosemary Munro*, Graeme Kelly, Michael Rohn* and Roger Saunders
ASSIMILATION OF METEOSAT RADIANCE DATA WITHIN THE 4DVAR SYSTEM AT ECMWF Rosemary Munro*, Graeme Kelly, Michael Rohn* and Roger Saunders European Centre for Medium Range Weather Forecasts Shinfield Park,
More informationThe impact of the assimilation of altimeters and ASAR L2 wave data in the wave model MFWAM
The impact of the assimilation of altimeters and ASAR L2 wave data in the wave model MFWAM Lotfi Aouf 1, Jean-Michel Lefèvre 1 1) Météo-France, Toulouse 12 th Wave Hindcasting and Forecasting, Big Island
More informationDelia Arnold 1, Arturo Vargas 1,Milagros Montero 2, Alla Dvorzhak 2 and Petra Seibert 3
COMPARISON OF THE DISPERSION MODEL IN RODOS-LX AND MM5-V3.7-FLEXPART(V6.2). A CASE STUDY FOR THE NUCLEAR POWER PLANT OF ALMARAZ Delia Arnold 1, Arturo Vargas 1,Milagros Montero 2, Alla Dvorzhak 2 and Petra
More informationFeature-tracked 3D Winds from Satellite Sounders: Derivation and Impact in Global Models
Feature-tracked 3D Winds from Satellite Sounders: Derivation and Impact in Global Models David Santek 1, A.-S. Daloz 1, S. Tushaus 1, M. Rogal 1, W. McCarty 2 1 Space Science and Engineering Center/University
More informationEfficient Forecasting of Volcanic Ash Clouds. Roger P Denlinger Hans F Schwaiger US Geological Survey
Efficient Forecasting of Volcanic Ash Clouds Roger P Denlinger Hans F Schwaiger US Geological Survey Two basic questions addressed in this talk: 1. How does uncertainty affect forecasts of volcanic ash
More informationA simple operative formula for ground level concentration from a point source
Air Pollution XX 3 A simple operative formula for ground level concentration from a point source T. Tirabassi, M. T. Vilhena & D. Buske 3 Institute of Atmospheric cience and Climate (IAC-CNR), Bologna,
More informationDATA FUSION NOWCASTING AND NWP
DATA FUSION NOWCASTING AND NWP Brovelli Pascal 1, Ludovic Auger 2, Olivier Dupont 1, Jean-Marc Moisselin 1, Isabelle Bernard-Bouissières 1, Philippe Cau 1, Adrien Anquez 1 1 Météo-France Forecasting Department
More informationSpatial Interpolation & Geostatistics
(Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Lag Mean Distance between pairs of points 1 y Kriging Step 1 Describe spatial variation with Semivariogram (Z i Z j ) 2 / 2 Point cloud Map 3
More informationClimate Modeling Dr. Jehangir Ashraf Awan Pakistan Meteorological Department
Climate Modeling Dr. Jehangir Ashraf Awan Pakistan Meteorological Department Source: Slides partially taken from A. Pier Siebesma, KNMI & TU Delft Key Questions What is a climate model? What types of climate
More informationCurrent best practice of uncertainty forecast for wind energy
Current best practice of uncertainty forecast for wind energy Dr. Matthias Lange Stochastic Methods for Management and Valuation of Energy Storage in the Future German Energy System 17 March 2016 Overview
More informationSIMULATION OF ATMOSPHERIC STATES FOR THE CASE OF YEONG-GWANG STORM SURGE ON 31 MARCH 2007 : MODEL COMPARISON BETWEEN MM5, WRF AND COAMPS
SIMULATION OF ATMOSPHERIC STATES FOR THE CASE OF YEONG-GWANG STORM SURGE ON 31 MARCH 2007 : MODEL COMPARISON BETWEEN MM5, WRF AND COAMPS JEONG-WOOK LEE 1 ; KYUNG-JA HA 1* ; KI-YOUNG HEO 1 ; KWANG-SOON
More informationAtmospheric Transport Modelling BY PETER CHEN AND RICHARD HOGUE
Atmospheric Transport Modelling BY PETER CHEN AND RICHARD HOGUE AT M O S P H E R I C Atmospheric Transport Modelling BY PETER CHEN AND RICHARD HOGUE FACT BOX Atmospheric Transport Modelling (ATM) is concerned
More informationImplementation of Modeling the Land-Surface/Atmosphere Interactions to Mesoscale Model COAMPS
DISTRIBUTION STATEMENT A: Approved for public release; distribution is unlimited. Implementation of Modeling the Land-Surface/Atmosphere Interactions to Mesoscale Model COAMPS Dr. Bogumil Jakubiak Interdisciplinary
More informationM.Sc. in Meteorology. Numerical Weather Prediction
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Cehtre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. Text for the Course
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 informationMeteorological monitoring system for NPP Kozloduy. Hristomir Branzov
Meteorological monitoring system for NPP Kozloduy Hristomir Branzov National Institute of Meteorology and Hydrology 66 Tzarigradsko schaussee, Sofia 784, Bulgaria Tel. (+359 2) 975 35 9, E-mail: Hristomir.Branzov@meteo.bg
More informationMAP Re-Analysis with IFS/ARPEGE/ALADIN
MAP Re-Analysis with IFS/ARPEGE/ALADIN Stjepan Ivatek-Šahdan RC LACE Data Manager Croatian Meteorological and Hydrological Service ivateks@cirus.dhz.hr Contents - Why to care about MAP-SOP Period? - Modification
More information