C a r b o

1
Carbo
S t r a t e g y
E u rope
R e gion a l
E x peri m e n t
ATMOPHERIC MODELLING AT THE REGIONAL SCALE DURING
CERES
Claire Sarrat, Joël Noilhan, Pierre Lacarrère,
Sylvie Donier
2
Evaluation of Meso-NH in  Operational  Mode

• CONFIGURATION
• 2 nested models
• First domain 240km - 240km, Resolution
  8km
• Second domain 60km - 60km, Resolution 2km
  
• Center of both domains La Cape Sud

• Each D day during CERES
• Forecast for D1
• Beginning of the simulation D1 at 3hTU
• End of the simulation D1 at 18hTU
• Initial fields and LBC ALADIN forecast
• Operational surface fields soil water
  content
• Meso-NH forecasted fields (for D1) were
  available daily at 8h30TU.

ALADIN forecast
27 H
00UTC
18UTC
3UTC
240 x 240 km Resolution 8 km
D1
D
D2
Meso-NH forecast
3
Evaluation of Meso-NH in  Operational  Mode
Comparisons T2M observed (squares) and simulated
(isolines)
05/27/05
06/02/05
4
Evaluation of Meso-NH in  Operational  Mode
Statistical scores calculated for 66
stations and for 19 IOP days
5
Evaluation of Meso-NH in  Operational  Mode
Comparisons from Radiosounding in La Cape Sud
IOP3 06/01/05
IOP4a 06/06/05 Lagrangian Exp.
6
CarboEurope Regional Experiment Evaluation of
Meso-NH in  Operational  Mode
Boundary layer height comparisons MESO-NH
simulation / observed RS in La Cape Sud
7
Evaluation of Meso-NH in  Operational  Mode
PIPER IN MARMANDE 1445 UTC
Theta
CO2
PIPER IN LA CAPE SUD 1400 UTC
Zi 1300m
Zi 700m
PBL MESONH 12 UTC 27.05.05 (36 h forecast)
8
Evaluation of Meso-NH in  Operational  Mode
IOP 4 a Lagangian , June, 6
MESO-NH forecast 35 H
STILT backward plume
9
Evaluation of Meso-NH in  Operational 
Mode Conclusion

• Encouraging results using Meso-NH in operational
  mode during CERES
• ? Dry and warm bias at low level in the
  morning to be analysed (ALADIN role to be
  confirmed)
• ? this bias is reduced during the day (role
  of the boundary layer development and
  surface processes?)
• ? Satisfactory simulation of the BL height
  at La Cape Sud and its regional variation
• ? Meso-NH proved very useful for the
  operations planning during CERES (Height of CO2
  flask sampling, choice of aircraft flights,
  Lagrangian trajectories)

10
The CNRM Modelling program during CERES

• Aims
• ? Atmospheric CO2 modelling in interaction with
  the surface coupling Meso-NH atmospheric CO2
  with CO2 surface flux from ISBA-A-gs
• -gt To improve surface-atmosphere exchanges
  coupling water and carbon cycles
• -gt To improve ABL modelling (vertical mixing,
  entrainment)
• -gt To estimate the regional budget of atmospheric
  CO2
• -gt To collaborate to CO2 sources and sinks
  retrieval from CO2 concentrations (collaboration
  with Ceres partners, LSCE, ( T. Louvaux PhD),
  Uni. Amsterdam, MPI, )

11
The CNRM Modelling program during CERES

• Strategy
• Implementation of CO2 in MesoNH-Isba-A-gs for
  Ceres
• ? Large scale forcing and initialisation of
  CO2 with LMDZ sensitivity to large scale
  forcing and to the spatio-temporal resolution of
  LMDZ
• ? Anthropogenic emissions from GENEMIS
  database
• ? New vegetation classification Ecoclimap
  South-West (Champeaux et al. 05)
• Modelling program
• ? further calibration of Isba A- gs for
  Ceres (0-D and 1-D)
• ? 1-D simulations analysis of the ABL-Surface
  interactions and entrainment for IOP with weak
  winds,
• ? 3-D runs at 8 km resolution for all IOPs
  sensitivity to initialisation (Arpege, Aladin,
  ECMWF), sensitivity to surface processes versus
  advection, first guess for inversion,
• ? 3-D runs, 2 ways nesting at high resolution
  2.5 km for selected days 27may, 1june, 6june,
  22june,
• ? CO2 budgets analysis

12
The CNRM Modelling program during CERES

• Surface Scheme ISBA-Standard

• Coupling Meso-NH and ISBA-A-gs

LAI
Met.Forcing MESO-NH CO2atm.
Active Biomass
ISBA-A-gs
LE, H, Rn, W, Ts
CO2 Flux
13
The CNRM Modelling program during CERES

• Resolution 8km
• Arpege initialisation and LBC
• 65 vertical levels
• Surface scheme ISBA-A-gs

14
The CNRM Modelling program during CERES
Anthropogenic CO2 emissions at 10km resolution
from Univ. of Stuttgart
Emissions on 1st of june 2000, 16H
15
The CNRM Modelling program during CERES

• ECOCLIMAP SW 62 classes of vegetation (winter
  vs. summer crops)

• Land-cover map processed from Corine 2000 and
  VEGETATION NDVI
• (Champeaux et al. 2005)

• Soon available for the CERES database
• ? Classification of 62 vegetation classes for
  CERES domain (SW)
• NDVI profiles decadal from 1999 to 2003 and
• from may
  to june 2005
• ? Mean NDVI profiles for 62 vegetation classes

16
Meso-NH / ISBA-A-gs 1-D Simulations 24may IOP2

• Simulation conditions
• Radiation
• No large scale forcing no horizontal
  advection, w 0
• 50 vertical levels including 30 below 2000m
• Surface scheme ISBA-A-gs
• Initial conditions at 500 UTC
• CO2 concentrations from Sky-Arrow (QL) at 700
  above La Cape Sud
• Meteorological parameters from RS in La Cape Sud
  at 500 UTC
• Sensitivity experiment to surface fluxes
• - with forest conditions
• - with maize conditions bare ground
• - 50 forest 50 maize

Problem in the simulation due to the lack of
subsidence estimation of
!
17
Meso-NH / ISBA-A-gs 1-D Simulations 24may
IOP2 Surface fluxes
Bowen ratio (2) observed and simulated in
agreement for both sites
18
Meso-NH / ISBA-A-gs 1-D Simulations 24may IOP2
THETA (C) Vertical profiles
19
CarboEurope Regional Experiment Meso-NH /
ISBA-A-gs Atm CO2 Modelling
3-D Modelisation pré-campaign Meso-NH et CO2
atm COCA Experiment juin 2001 Comparisons
Obs-Simu Le Bray (forest)

• RECAB calibration of parameters
• Gm 0.5 mm.s-1
• RE25 1.5 10-7 kg/kg.m.s-1

20
CarboEurope Regional Experiment Meso-NH /
ISBA-A-gs Atm CO2 Modelling

• CO2 gradient forest- agricultural areas
  simulated confirmed during the campaign

CO2
H2O
? simulated CO2 at 15H (june 2001) Advection
Assimilation vertical mixing vertical ?CO2
decrease
21
CarboEurope Regional Experiment Meso-NH in
 Operational  Mode Evaluation
Statistical scores calculated for 66 stations and
19 IOP days
22
(No Transcript)
23
(No Transcript)
24
(No Transcript)
25
Meso-NH / ISBA-A-gs 1-D Simulations 24may IOP2
Surface fluxes for both sites
Bowen ratio (2) observed and simulated in
agreement