Le simulazioni del Progetto MEGAPOLI:
confronti con le simulazioni di MINNI
A. D’Allura, S. Finardi, P. Radice, C. Silibello
1-2 Nov 2010 - 2nd Annual MEGAPOLI Meeting (Hamburg, Germany)
Messaggi
Le nostre simulazioni producono una quantità di dati che spesso non
utilizziamo e che forniscono, invece, informazioni molto utili. Alcuni esempi
forniti dalle simulazioni effettuate da ARIANET all’interno del progetto
MEGAPOLI;
Utilizzo di dati di AOD (dati MODIS/rete AERONET) per una più completa
valutazione delle nostre simulazioni. Un esempio dal progetto MEGAPOLI;
“Data assimilation” ?
Il confronto delle simulazioni MINNI/Nord-Italia e MEGAPOLI/Po Valley
suggerisce
l’ipotesi
che
simulazioni
innestate
EUROPA/ITALIA/MACROREGIONI possano fornire stime più corrette
relativamente al PM
Le nostre simulazioni
producono una
quantità di dati che
spesso non utilizziamo
e che forniscono,
invece, informazioni
molto utili. Alcuni
esempi forniti dalle
simulazioni effettuate
da all’interno del
progetto MEGAPOLI
Po Valley simulations
Nested domains
G1
G2
Po Valley simulations
(June, December 2005)
Runs over Po Valley considering two nested domains.
FARM:
•SAPRC99 gas phase chemical mechanism / aero3 aerosol module
•Rosenbrock solver– ROS3
•Space resolution: 16 and 4 km (horizontal); 16 levs, up to 10 km (vertical)
Emissions:
•Anthropogenic: TNO data set ( 7 km resolution); Italian data set (ISPRA - Province
level); Regional data sets (INEMAR - Municipal level);
•BIOGENIC: SURFPro Isoprene and Terpenes from vegetation, PM (fine and coarse) from
aeolian resuspension and sea salts, wind influence;
Meteorology:
•RAMS;
IC/BC:
•MPI MATCH – Global scale: gaseous species
•GOCART – Global scale: Climatological aerosols
Po Valley simulations
Emission inventories
Italy - ISPRA
(Province level)
Regions - INEMAR
(Municipal level)
TNO data set
( 7 km resolution)
PROCESS CONTRIBUTIONS
Computed from hourly balance terms computed
runtime over inner model grid
Po Valley simulations
Process contribution to PM2.5 [g]
In/Out flows
4.E+10
Gain/Inflow
3.E+10
2.E+10
1.E+10
0.E+00
-3.E+10
Loss/Outflow
-4.E+10
Hor_inflow
June 2005
Hor_outflow
Top_flow
Net_flow
Net_flow = Hor_inflow + Hor_outflow + Top flow
179
-2.E+10
177
175
173
171
169
167
165
163
161
159
157
155
153
151
-1.E+10
Po Valley simulations
Process contribution to PM2.5 [g]
Other processes
2.0E+10
1.5E+10
BIOGENIC SOA
1.0E+10
5.0E+09
0.0E+00
Emitted
June 2005
Dep
Aerosol
179
-2.0E+10
177
-1.5E+10
175
-1.0E+10
173
171
169
167
165
163
161
159
157
155
153
151
-5.0E+09
Po Valley simulations
Process contribution to PM2.5 [g]
Biogenic emissions (terpenes) -> SOA
1.8E+08
1.6E+08
TRP1 emissions
1.4E+08
1.2E+08
1.0E+08
8.0E+07
6.0E+07
4.0E+07
2.0E+07
0.0E+00
180
178
176
174
172
170
168
166
164
162
160
158
156
154
152
3.00E+09
2.50E+09
BIOGENIC SOA
2.00E+09
1.50E+09
1.00E+09
5.00E+08
180
178
176
174
172
170
168
166
164
162
160
158
156
154
152
June 2005
0.00E+00
1.50E+09
300
1.25E+09
295
1.00E+09
290
7.50E+08
285
5.00E+08
280
2.50E+08
275
0.00E+00
Isoprene emissions: during
daylight between April and
October;
Terpenes emissions: full year
180
178
305
3.00E+09
300
2.50E+09
295
2.00E+09
290
1.50E+09
285
1.00E+09
280
5.00E+08
275
0.00E+00
180
178
176
174
Emission
172
170
168
166
164
162
160
158
156
154
Tavg
Terpenes Emissions [g]
Emission
152
June 2005
176
174
172
170
Isoprene/Terpene emissions
(daily integral over all
domain)
168
Tavg: ground level average
temperature at noon;
166
164
162
160
158
156
154
152
Tavg
Temperature dependence;
Isoprene Emissions [g]
305
T [K]
T [K]
Po Valley simulations
Biogenic emissions vs Temperature
Po Valley simulations
Process contribution to O3
Chemistry vs Temperature
2.0E+10
T [K]
300
1.5E+10
295
1.0E+10
290
5.0E+09
0.0E+00
285
180
178
176
174
172
170
168
June 2005
166
164
162
160
158
156
154
152
Tavg
Chemistry
Tavg: ground level average temperature at noon, O3 Chemistry: daily integral
O3 Chemistry [g]
2.5E+10
305
Po Valley simulations
Process contribution to CO [g]
CO production from Isoprene/terpenes ozonolysis
{190} ISOPRENE + O3 = 0.266 OH + 0.066 RO2_R + 0.008
RO2_N + 0.126 R2O2+ 0.192 MA_RCO3 + 0.275 CO +
0.592 HCHO + 0.1 PROD2 + 0.39 METHACRO + 0.16
MVK + 0.204 HCOOH + 0.15 RCO_OH
{194} TERP + O3 = 0.567 OH + 0.033 HO2 + 0.031 RO2_R
+ 0.18 RO2_N + 0.729 R2O2+ 0.123 CCO_O2 + 0.201
RCO_O2 + 0.157 CO + 0.235 HCHO + 0.205RCHO +
0.13ACET + 0.276PROD2 + 0.001GLY + 0.031 BACL +
0.103 HCOOH + 0.189 RCO_OH + TRP1AER (->AORB)
Po Valley simulations
Process contribution to CO [g]
Other processes
Chemistry ?
1.5E+09
1.0E+09
5.0E+08
0.0E+00
-1.5E+09
Emission
June 2005
Chemistry
Net_flow
Dry deposition is negligible and wet deposition is not considered for CO
180
-1.0E+09
179
178
177
176
175
174
173
172
171
170
169
168
167
166
165
164
163
162
161
160
159
158
157
156
155
154
153
152
-5.0E+08
Po Valley simulations
Process contribution to CO [mol]
Chemistry/interpretation
Chemistry [mol]
8.0E+06
6.0E+06
4.0E+06
2.0E+06
0.0E+00
June 2005
ISOP
TRP1
180
CO
179
-8.0E+06
178
-6.0E+06
177
-4.0E+06
176
175
174
173
172
171
170
169
168
167
166
165
164
163
-2.0E+06
Particle dry deposition velocity
Po Valley simulations
PM: Process contribution/dimensional analysis
Net flows [g]
Accumulation mode particles
may be transported outside
Po Valley region (because of
lower dry deposition).
Inflow
179
177
175
173
171
169
167
165
163
161
159
157
155
153
151
5.E+09
4.E+09
3.E+09
2.E+09
1.E+09
0.E+00
-1.E+09
-2.E+09
-3.E+09
-4.E+09
-5.E+09
Outflow
Aitken
Accumulation
Coarse
0
179
177
175
173
171
169
167
165
163
161
159
-50
-100
-150
-200
-250
Accumulation
Coarse
Aitken
Aitken mode
157
155
153
June 2005
0.00
-1.00
-2.00
-3.00
-4.00
-5.00
-6.00
-7.00
-8.00
-9.00
-10.00
151
Negligible outflow from Po
Valley region, for Aiken and
Coarse mode particles, due to
dry deposition process.
Accumulation/Coarse
modes
Dry deposition / Emission
Po Valley simulations
Process contribution to HNO3 [g]
Other processes
HNO3 nigh time production
due to chemistry
less dissolution of HNO3 in
the particles?
2.5E+08
2.0E+08
1.5E+08
1.0E+08
5.0E+07
0.0E+00
180
Dep
179
June 2005
Chemistry
178
Aerosol
177
-2.0E+08
176
-1.5E+08
175
-1.0E+08
174
173
172
171
170
169
168
167
166
165
164
163
162
161
160
159
158
157
156
155
154
153
152
-5.0E+07
Po Valley simulations
Process contribution to HNO3 [mol]
Aerosol processes
2.5E+06
2.0E+06
1.5E+06
1.0E+06
5.0E+05
0.0E+00
180
179
178
177
176
175
174
173
172
171
170
169
168
167
166
165
164
163
162
161
160
159
158
157
156
155
154
153
152
-5.0E+05
-1.0E+06
-1.5E+06
-2.0E+06
-2.5E+06
170
169
168
5.0E+06
0.0E+00
-5.0E+06
167
5.0E+05
0.0E+00
-5.0E+05
166
1.5E+07
1.0E+07
165
1.5E+06
1.0E+06
164
June 2005
2.5E+07
2.0E+07
163
The decrease in the
PM
liquid
water
content leads to less
dissolution of HNO3
in the particles
2.5E+06
2.0E+06
-1.0E+06
-1.5E+06
-1.0E+07
-1.5E+07
-2.0E+06
-2.5E+06
-2.0E+07
-2.5E+07
HNO3
Nitrate (aerosol)
Water (aerosol)
Po Valley simulations
Process contribution to HNO3 [g]
Other processes
1.5E+08
1.0E+08
5.0E+07
0.0E+00
0
364
363
362
Dep
361
Dec. 2005
Chemistry
360
Aerosol
359
-1.5E+08
358
-1.0E+08
357
356
355
354
353
352
351
350
349
348
347
346
345
344
343
342
341
340
339
338
337
336
335
-5.0E+07
Po Valley simulations
Process contribution to HNO3 [g]
Aerosol processes
1.5E+08
1.0E+08
5.0E+07
0.0E+00
0
364
363
362
361
360
359
358
357
356
355
354
353
352
351
350
349
348
347
346
345
344
343
342
341
340
339
338
337
336
335
-5.0E+07
-1.0E+08
3.5E+05
6.0E+07
3.0E+05
5.0E+07
2.5E+05
4.0E+07
2.0E+05
3.0E+07
1.5E+05
2.0E+07
1.0E+05
1.0E+07
5.0E+04
0.0E+00
0.0E+00
362
360
358
SO4
356
354
NO3
352
350
NH4
348
346
344
342
340
338
336
334
Dec. 2005
7.0E+07
SO4
No further HNO3 reacts
with ammonia to form
volatile salts
NH4, NO3
-1.5E+08
Aerosol types
RSO4
NH 


4

4
SO
from [Nenes et al., 1998 ] …based on the value of this ratio [in the
simulations we do not include sodium and chloride], following types of
aerosols are defined:
• Sulfate rich (free acid): This is when RSO4 < 1. The sulfates are in
abundance and part of it is in the form of free sulfuric acid. In this case,
there is always a liquid phase, because sulfuric acid is extremely
hygroscopic (i.e., DRH is 0%);
• Sulfate rich (non free acid): This is when 1 RSO4 < 2. There is enough
ammonia to partially (but not fully) neutralize the sulfates. The sulfates
are a mixture of bisulfates and sulfates, the ratio of which is
determined by thermodynamic equilibrium;
• Sulfate poor: RSO4  2: There is enough ammonia to fully neutralize the
sulfates. In this case, excess ammonia can react with the other species
(HNO3, HCl) to form volatile salts.
June 2005
RSO4

NH 


4

4
SO
Sulfate rich
(1 RSO4 < 2)
Scen_50
Scen_0
Dec. 2005
RSO4

NH 


4

4
SO
Sulfate poor
(RSO4  2)
Scen_50
Scen_0
NO-3
Scen_0
Scen_50
ANO3
1.00
Absolute differences
with the “Base case"
-50%
0.75
0.50
50% reduction of emis.
leads to a reduction of
conc. greater than 50%
0.25
0.00
2.00
Annihiliation (-100%)
1.50
1.00
0.50
0.00
June 2005
NO-3
Scen_0
Scen_50
ANO3
2.5
50% reduction of emis. leads
to a reduction of conc. lower
than 50%
2.0
Absolute differences
with the “Base case"
-50%
1.5
1.0
0.5
0.0
5
Annihiliation (-100%
4
3
2
1
0
Dec. 2005
Po Valley simulations
Po-Valley tracers emissions
“CO” and “PM fine” tracers are emitted only in the Po Valley region (green); are transported and
deposited (dry and wet); chemical and aerosol precesses are not considered (chemically not
reactive).
Po Valley simulations
“CO” tracer
CO
CO tracer
Base case – Scen_0
Base case
Absolute differences
June 2005
Base case concentrations
Po Valley simulations
“CO” tracer
CO
CO tracer
Base case – Scen_0
Base case
Absolute differences
Dec. 2005
Base case concentrations
Utilizzo di dati di AOD (dati MODIS/rete
AERONET) per una più completa valutazione
delle nostre simulazioni.
Un esempio dal progetto MEGAPOLI
AOD, Jan
FARM
Terra (MODIS)
I campi di AOD satellitari vengono confrontati con I campi medi mensili prodotti da FARM relativi
alle ore 12 (media mobile esaoraria). I campi medi osservati potrebbero non essere congruenti
con quelli stimati a seguito delle possibile mancanza di dati satellitari.
AOD, Jun
FARM
Terra (MODIS)
AERONET Stations
AOD at 550 nm
0.6
Palencia
El Arenosillo
Dec-05
Nov-05
Oct-05
Dec-05
Nov-05
Oct-05
Sep-05
Aug-05
Sep-05
Aug-05
Jul-05
Apr-05
May-05
Jun-05
Jul-05
Avignon
Barcelona
Jun-05
May-05
Apr-05
Mar-05
0.6
Feb-05
Jan-05
Dec-05
Nov-05
0.7
0.6
Mar-05
AOD at 550 nm
0.7
0.6
Feb-05
Jan-05
Nov-05
Oct-05
Sep-05
Aug-05
Jul-05
Jun-05
May-05
Apr-05
Mar-05
Feb-05
Dec-05
Nov-05
Oct-05
lug-05
Sep-05
Aug-05
Jul-05
Jun-05
giu-05
May-05
Apr-05
Mar-05
Feb-05
Jan-05
mag-05
apr-05
gen-05
Jan-05
Feb-05
Mar-05
feb-05
Apr-05
May-05
Jun-05
Jul-05
mar-05
Aug-05
Sep-05
Oct-05
Nov-05
Dec-05
0
0
0.5
AOD at 550 nm
AOD at 550 nm
0.7
MODIS
AERONET
FARM
Oct-05
Sep-05
Aug-05
Jul-05
Jun-05
May-05
Apr-05
Mar-05
AOD at 550 nm
0
Dec-05
0.1
0
Feb-05
Jan-05
0.1
0.1
0
0.1
Jan-05
0.5
0.2
0.2
0.1
0.1
0.1
0.3
0.3
0.3
0.3
0.4
0.4
0.2
0.7
0.3
0
0
0.6
Cabauw
Palaiseu (Paris)
0.4
0.4
0.4
0.2
0.2
0.5
0.5
0.3
0.3
0.6
0.6
0.4
0.4
0.7
Cabo da Roca
0.7
0.5
0.5
0.2
0.2
0.7
0.5
AOD at 550 nm
AOD: FARM vs MODIS and Aeronet, data
AOD at 550 nm
Mainz
Karlsuhe
Dec-05
Nov-05
Oct-05
Sep-05
Aug-05
Jul-05
Belsk
SMHI
Jun-05
Dec-05
Nov-05
Oct-05
Sep-05
Aug-05
Jul-05
Jun-05
May-05
Dec-05
Nov-05
lug-05
Oct-05
Sep-05
Aug-05
Jul-05
Jun-05
giu-05
May-05
Apr-05
Mar-05
Feb-05
gen-05
Jan-05
Feb-05
Mar-05
feb-05
Apr-05
May-05
Jun-05
Jul-05
mar-05
Aug-05
Sep-05
Oct-05
Nov-05
Dec-05
Jan-05
mag-05
0
apr-05
0.1
0
Apr-05
0.6
Mar-05
0.7
0.6
May-05
0
Apr-05
0.1
0
Mar-05
AOD at 550 nm
0.6
Feb-05
Jan-05
Dec-05
Nov-05
0.5
AOD at 550 nm
AOD at 550 nm
0.7
0.6
Feb-05
Jan-05
Dec-05
Nov-05
Oct-05
Sep-05
Aug-05
Jul-05
Jun-05
May-05
Apr-05
Mar-05
0.2
0.2
0.2
0.2
AOD at 550 nm
0.7
MODIS
AERONET
FARM
Oct-05
Sep-05
Aug-05
Jul-05
Jun-05
May-05
Apr-05
Mar-05
Feb-05
Jan-05
0.1
0
Erdemli
Minsk
0.4
0.4
0.4
0.3
0.3
0.1
Feb-05
Jan-05
0.5
0.5
0.5
0.6
0.6
0.2
0.2
0.7
Cabo da Roca
0.7
0.3
0.3
0
0
0.4
0.4
0.4
0.4
0.2
0.7
0.3
0.1
0.6
0.1
0.1
0.5
0.5
0.3
0.3
0.7
0.5
AOD at 550 nm
AOD: FARM vs MODIS Aeronet, data
AOD at 550 nm
Venice
Ispra
Dec-05
Nov-05
Oct-05
Nov-05
Dec-05
Nov-05
Dec-05
Sep-05
Oct-05
lug-05
Sep-05
Aug-05
Oct-05
Aug-05
Jul-05
Jun-05
giu-05
Jul-05
Jun-05
May-05
Apr-05
Mar-05
Feb-05
Jan-05
mag-05
0
gen-05
Jan-05
Feb-05
feb-05
Mar-05
Apr-05
May-05
Jun-05
Jul-05
mar-05
Aug-05
Sep-05
Oct-05
Nov-05
Dec-05
apr-05
0.1
0
Sep-05
0
May-05
0.6
Aug-05
0
Apr-05
0.7
0.6
Jul-05
Rome
Modena
Jun-05
0.1
Mar-05
0
May-05
0.1
Feb-05
Apr-05
May-05
Jun-05
Jul-05
Aug-05
Sep-05
Oct-05
Nov-05
Dec-05
Jan-05
0.1
0
Apr-05
AOD at 550 nm
0.6
Mar-05
0.5
AOD at 550 nm
AOD at 550 nm
0.7
0.6
Feb-05
AOD at 550 nm
0.7
MODIS
AERONET
0.2
0.2
Jan-05
Dec-05
Nov-05
Oct-05
Sep-05
Aug-05
Jul-05
Jun-05
May-05
Apr-05
Mar-05
Feb-05
Jan-05
0.1
0
Mar-05
Lampedusa
Lecce
0.4
FARM
0.3
0.3
0.1
Feb-05
Jan-05
0.5
0.2
0.2
0.4
0.4
0.2
0.7
0.3
0.1
0.6
0.4
0.4
0.2
0.2
0.5
0.5
0.3
0.3
0.6
0.6
0.4
0.4
0.7
Cabo da Roca
0.7
0.5
0.5
0.3
0.3
0.7
0.5
AOD at 550 nm
AOD: FARM vs MODIS Aeronet, data (ITALY)
EMEP PM10 Stations
PM10 - daiy average (June) [ g m-3]
120
120
100
100
80
80
Predicted
Predicted
PM10 - daiy average (January) [ g m-3]
60
60
40
40
20
20
0
0
120
100
Observed
80
60
40
20
0
120
100
80
60
40
20
0
Observed
O3 EMEP Network
Ozone - max 1 h [ g m-3]
240
200
Predicted
160
120
80
40
0
240
200
160
120
80
40
0
Observed
FARM
SILAM
1.7.2009, NO2 near-surface concentrations from FARM-ARIANET and SILAM-FMI
O3
Comparison with some Norvegian observations
http://tarantula.nilu.no/projects/ccc/emepdata.html
NO0001 Birkenes
NO0039 Kårvatn
NO0041R Osen
NO0043 Prestebakke
NO0048 Voss
NO0052 Sandve
NO0056 Hurdal
-3
g m
250
250
200
200
Prestebakke
Birkenes
18/07/2009 00:00
19/07/2009 00:00
20/07/2009 00:00
21/07/2009 00:00
22/07/2009 00:00
23/07/2009 00:00
24/07/2009 00:00
25/07/2009 00:00
26/07/2009 00:00
27/07/2009 00:00
28/07/2009 00:00
29/07/2009 00:00
30/07/2009 00:00
31/07/2009 00:00
17/07/2009 00:00
18/07/2009 00:00
19/07/2009 00:00
20/07/2009 00:00
21/07/2009 00:00
22/07/2009 00:00
23/07/2009 00:00
24/07/2009 00:00
25/07/2009 00:00
26/07/2009 00:00
27/07/2009 00:00
28/07/2009 00:00
29/07/2009 00:00
30/07/2009 00:00
31/07/2009 00:00
12/07/2009 00:00
13/07/2009 00:00
14/07/2009 00:00
15/07/2009 00:00
16/07/2009 00:00
17/07/2009 00:00
16/07/2009 00:00
15/07/2009 00:00
14/07/2009 00:00
Hurdal
Sandve
13/07/2009 00:00
11/07/2009 00:00
0
12/07/2009 00:00
10/07/2009 00:00
0
11/07/2009 00:00
09/07/2009 00:00
50
10/07/2009 00:00
08/07/2009 00:00
07/07/2009 00:00
06/07/2009 00:00
05/07/2009 00:00
04/07/2009 00:00
03/07/2009 00:00
02/07/2009 00:00
50
01/07/2009 00:00
100
01/07/2009 00:00
02/07/2009 00:00
03/07/2009 00:00
04/07/2009 00:00
05/07/2009 00:00
06/07/2009 00:00
07/07/2009 00:00
08/07/2009 00:00
09/07/2009 00:00
10/07/2009 00:00
11/07/2009 00:00
12/07/2009 00:00
13/07/2009 00:00
14/07/2009 00:00
15/07/2009 00:00
16/07/2009 00:00
17/07/2009 00:00
18/07/2009 00:00
19/07/2009 00:00
20/07/2009 00:00
21/07/2009 00:00
22/07/2009 00:00
23/07/2009 00:00
24/07/2009 00:00
25/07/2009 00:00
26/07/2009 00:00
27/07/2009 00:00
28/07/2009 00:00
29/07/2009 00:00
30/07/2009 00:00
31/07/2009 00:00
100
09/07/2009 00:00
200
08/07/2009 00:00
200
07/07/2009 00:00
250
06/07/2009 00:00
30/07/2009 00:00
29/07/2009 00:00
28/07/2009 00:00
27/07/2009 00:00
26/07/2009 00:00
25/07/2009 00:00
24/07/2009 00:00
23/07/2009 00:00
22/07/2009 00:00
21/07/2009 00:00
20/07/2009 00:00
19/07/2009 00:00
18/07/2009 00:00
17/07/2009 00:00
16/07/2009 00:00
15/07/2009 00:00
14/07/2009 00:00
13/07/2009 00:00
12/07/2009 00:00
11/07/2009 00:00
10/07/2009 00:00
09/07/2009 00:00
08/07/2009 00:00
07/07/2009 00:00
06/07/2009 00:00
250
05/07/2009 00:00
0
05/07/2009 00:00
300
04/07/2009 00:00
0
04/07/2009 00:00
300
03/07/2009 00:00
50
03/07/2009 00:00
-3
300
02/07/2009 00:00
50
02/07/2009 00:00
g m
-3
g m
300
01/07/2009 00:00
100
31/07/2009 00:00
100
01/07/2009 00:00
-3
150
150
150
150
g m
O3
Comparison with some Norwegian observations
Data assimilation
Sistema MINNI
Run “NI1” 2005 a 4 km
Test assimilazione medie annuali BRACE
Sistema MINNI
Concentrazioni medie annuali di NO2
Sistema MINNI
Medie annuali NO2 - 2005
Run “IT1” a 4 km (inventario ISPRA)
52
48
44
40
36
32
28
24
20
16
12
8
4
0
[ug/m3]
Sistema MINNI
Medie annuali NO2 - 2005
Run “IT1” a 4 km (inventario ISPRA) + assimilazione
52
48
44
40
36
32
28
24
20
16
12
8
4
0
[ug/m3]
Sistema MINNI
Confronto speditivo con database BRACE
In via preliminare sono state prese in considerazione tutte le stazioni BRACE
disponibili, senza selezione e controlli di coerenza; anche la tipologia
fondo/industriale/traffico e rurale/suburbana/urbana è stata mutuata direttamente
da quella presente nella base dati.
Sistema MINNI
Medie annuali NO2 - 2005
Run “IT1” a 4 km (inventario ISPRA)
Confronto su tutte le stazioni di fondo (122)
100
MILANO - JUVARA
80
CORMANO
MILANO - P.CO
LAMBRO
calc. [g/m3]
ARESE
SARONNO - SANTUARIO
60
C.SO FIRENZE
TO_1272_TO_LINGOTTO
BUSTO ARSIZIO - ACCAM
LIMITO
MONTE SAN PANTALEONE
MAGENTA
BRESCIAVF
- VIA ZIZIOLA
QUARTIERE ITALIA
OSIO
SOTTO
NO_3106_VERDI
ACQUASOLAARCONATE
MANDRIA
- VIA MEUCCI
PAVIA -BERGAMO
VIA FOLPERTI
VIMERCATE
PARCO BISSUOLA
TO_1309_VINOVO
PRCTDLA
CANTU' - VIAMODENA
MEUCCI - NONANTOLANA
MOTTA
VISCONTI
TREZZO
D'ADDA
BERGAMO
-TO_1171_ORBASSANO
VIA
GOISIS - XX SETTEMBRE
MODENA
LONATO
CASON
CREMA
VIA
INDIPENDENZA
BZ1 VIA- VIA
AMBA
ALAGI
S.OSVALDOCAIROLI
VIALE RISORGIMENTO
CREMA
XI FEBBRAIO
TRENTO GAR
TRENTO
VARESE
-SAN
VIA VIDOLETTI
ROCCO TRENTO
ALPSC
PORTO
CN_4003_ALBA
VEN
VENEZIA
SACCA
FISOLA
AT_5005_DACQUISTO
CN_4078_CUNEO
ERBAMODENA - CARPI 2
ROVERETO LGP
QUARTONO_3156_VERBANIA
SALORNO
GIOVI
LUCINICO
CORTE
DEI
CORTESI
MAGGIOLINA
FOLLO
SAREZZO - VIA MINELLI
CHIAPPA
AT_5012_BUTTIGLIERA
CENGIO1
BELLUNO-CITTA'
GAMBARA
BI_2012_BIELLA1
CN_4201_SALICETO
LUGAGNANO
RIVA GAR
BI_2046_COSSATO
VIPITENO
MOCOMERO
TO_1099_MANDRIA
SONDRIO - VIA MERIZZI
VC_2016_BORGOSESIA
VARENNA
RENON
BORGO VAL
NO_3118_PIEVEVERGONTE
COLICO
Febbio
LACES
LA THUILE
MONTE
ZONCOLAN
- SUTRIO
CHIAVENNA
MONTE
GAZA
BORMIO
40
20
0
0
20
40
60
meas. [g/m3]
80
100
Sistema MINNI
Medie annuali NO2 - 2005
Run “IT1” a 4 km (inventario ISPRA) + assimilazione
Confronto su tutte le stazioni di fondo (122)
100
80
MILANO - JUVARA
calc. [g/m3]
CORMANO
ARESE
TO_1272_TO_LINGOTTO
60
BRESCIA - VIA ZIZIOLA
VIALE RISORGIMENTO
BERGAMO
VIA MEUCCI
MILANO
- -P.CO
LAMBRO
NO_3106_VERDI
TO_1309_VINOVO
SARONNO
- SANTUARIO
BUSTO
ARSIZIO
- ACCAM
TO_1171_ORBASSANO
MONTE SAN PANTALEONE
MAGENTA VF
C.SO FIRENZE
PAVIA - VIA FOLPERTI
OSIO SOTTOMANDRIA
ARCONATE
BERGAMO -QUARTIERE
VIA GOISIS ITALIA
VIMERCATE
PRCTDLA
TRENTO GAR
CANTU' - VIA MEUCCI
ACQUASOLA
MOTTA
VISCONTI
MODENA
NONANTOLANA
AT_5005_DACQUISTO
TRENTO
PSC
PARCO CREMA
BISSUOLA
- SAN
VIA INDIPENDENZA
ROCCO TRENTO
AL
PORTO
MODENA
- XX
SETTEMBRE
VEN
CN_4003_ALBA
CASON
CREMA
XI FEBBRAIO
BZ1 VIA- VIA
AMBA
ERBA ALAGI
LONATO
VARESE
- VIA
VIDOLETTI
TREZZO
CN_4078_CUNEO
D'ADDA
ROVERETO LGP
NO_3156_VERBANIA
LIMITO
LUCINICO
CAIROLI
MODENA
- CARPI 2
LUGAGNANO
BI_2012_BIELLA1
RIVA
CORTEGIOVI
DEI
CORTESI
VENEZIA
SACCA GAR
FISOLA
S.OSVALDO
GAMBARA
SALORNO
BI_2046_COSSATO
SAREZZO
- VIA MINELLI
CENGIO1
BELLUNO-CITTA'
MOCOMERO QUARTO
CN_4201_SALICETO
AT_5012_BUTTIGLIERA
VC_2016_BORGOSESIA
FOLLO
MAGGIOLINA
VARENNA
TO_1099_MANDRIA
CHIAPPA
VIPITENO
SONDRIO
- VIA MERIZZI
COLICO
RENON
BORGO VAL
NO_3118_PIEVEVERGONTE
Febbio
LACES
MONTE
GAZA
CHIAVENNA
THUILE
MONTELA
ZONCOLAN
- SUTRIO
BORMIO
40
20
0
0
20
40
60
meas. [g/m3]
80
100
Sistema MINNI
Concentrazioni medie annuali di PM10
Sistema MINNI
Medie annuali PM10 - 2005
Run “IT1” a 4 km (inventario ISPRA)
52
48
44
40
36
32
28
24
20
16
12
8
4
0
[ug/m3]
Sistema MINNI
Medie annuali PM10 - 2005
Run “IT1” a 4 km (inventario ISPRA) + assimilazione
52
48
44
40
36
32
28
24
20
16
12
8
4
0
[ug/m3]
Sistema MINNI
Confronto speditivo con database BRACE
In via preliminare sono state prese in considerazione tutte le stazioni BRACE
disponibili, senza selezione e controlli di coerenza; anche la tipologia
fondo/industriale/traffico e rurale/suburbana/urbana è stata mutuata direttamente
da quella presente nella base dati.
Sistema MINNI
Medie annuali PM10 - 2005
Run “IT1” a 4 km (inventario ISPRA)
Confronto su tutte le stazioni di fondo (122)
100
MILANO - JUVARA
calc. [g/m3]
80
60
ARESE
SARONNO - SANTUARIO
MANDRIA
LIMITO
MODENA - XX SETTEMBRE
MAGENTA
VFDI NOVARA
TREVISO - VIA
LANCIERI
VIALE RISORGIMENTO
VILLA ADA
MANTOVA
VIA
VENEZIA SACCA FISOLA ARIOSTO
MODENA - CARPI
2 - VIA
CREMA
INDIPENDENZA
OSIO
SOTTO
VIMERCATE
CANTU'AL
- VIA
MEUCCI
CASON
SAN
ROCCO
PORTO
BERGAMO
- VIA
MEUCCI
40
GHERARDI
20
ERBA
TORCHIAROLO
SAREZZO
- VIA MINELLI
CENVS1S. Pancrazio
Salentino
CENSG1
SONDRIO - VIA MERIZZI
VIPITENO
FONTECHIARI
LACES
BORMIO
0
0
20
40
60
meas. [g/m3]
80
100
Sistema MINNI
Medie annuali PM10 - 2005
Run “IT1” a 4 km (inventario ISPRA) + assimilazione
Confronto su tutte le stazioni di fondo (122)
100
80
calc. [g/m3]
MILANO - JUVARA
60
ARESE
MAGENTA
MANTOVA
- VIA VF
ARIOSTO
MANDRIA
SARONNO
- SANTUARIO
OSIO SOTTO
TREVISO - VIA LANCIERI DI NOVARA
CASON - VIA MEUCCI
BERGAMO
MODENA - XX SETTEMBRE
CANTU' - VIA MEUCCI
VIALE RISORGIMENTO
SAN ROCCO AL PORTO
40
VIMERCATE
MODENA - CARPI
2 - VIA INDIPENDENZA
CREMA
VENEZIA
SACCA
FISOLA
LIMITO
ERBA
GHERARDI
SAREZZO - VIA MINELLI
20
SONDRIO - VIA MERIZZI
VIPITENO
BORMIO
LACES
0
0
20
40
60
meas. [g/m3]
80
100
Il confronto delle simulazioni MINNI/Nord-Italia e
MEGAPOLI/Po Valley suggerisce l’ipotesi che
simulazioni innestate
EUROPA/ITALIA/MACROREGIONI possano fornire
stime più corrette relativamente al PM
MEGAPOLI-Po Valley simulations
Nested domains
•
•
•
•
•
MEGAPOLI
2-way nesting (Central EUROPE-Po
Valleey);
Meteo: RAMS (venti piu’ intensi)
SAPRC99 KPP
Emissions: TNO, INEMAR, ISPRA 2005
BC: MPIMATH (gas), GOCART (aerosol)
MODELS:
FARM (MEGAPOLI and MINNI)
EMEP
MPI-MATCH
GOCART (aerosol, climatological)
Monthly averages (June)
PM2.5
MINNI
MEGAPOLI
EMEP
GOCART
PM10
MINNI
MEGAPOLI
EMEP
GOCART
PM2.5-10
MINNI
MEGAPOLI
EMEP
GOCART
O3
MINNI
MEGAPOLI
EMEP
MPIMATCH
MEGAPOLI-Po Valley simulations
Influence of BCs:
Nested domains vs 1-way nesting
1-5/6/2005: 5 days average
RUNS
• Nested domains: 2-way nesting
(Central EUROPE-Po Valley) using
MPIMATH
(gas),
GOCART
(aerosol) fields for BCs
• 1-way nesting: EMEP fields for
BCs
PM2.5
(average 1-5/6/2005)
Nested domains
1-way nesting (EMEP_BCs)
PM10
(average 1-5/6/2005)
Nested domains
1-way nesting (EMEP_BCs)
PM2.5-10
(average 1-5/6/2005)
Nested domains
1-way nesting (EMEP_BCs)
O3
(average 1-5/6/2005)
Nested domains
1-way nesting (EMEP_BCs)
Messaggi
La versione “OpenSource” di FARM potrebbe essere applicata alla
campagna della TRISAIA per verificarne le potenzialità (es. meccanismo
chimico più aggiornato e utilizzo del modulo di trasferimento radiativo).
Nelle prossime simulazioni MINNI utilizzare il file di output “balance.dat”.
Contiene informazioni utili alla interpretazione dei risultati!!
L’utilizzo dei dati di AOD e più in generale dei campi prodotti dai satelliti
fornisce informazioni che possono supportare la fase di valutazione del
sistema modellistico (confronto tra campi).
il modulo di assimilazione implementato in FARM consente di integrare
ancor più compiutamente il “sistema MINNI”, gli inventari delle emissioni ed
il dataset BRACE, al fine di produrre mappe di qualità dell’aria sul territorio
nazionale e sulle macroregioni sempre più realistici (non appropriato per
studi di scenario).
L’inserimento di un dominio a scala europea potrebbe determinare una
migliore stima dei livelli di inquinamento sul territorio nazionale.
Scarica

Po Valley simulations