october 1-3, 2007 1 6th annual cmas meeting comparison of cmaq and camx for an annual simulation...

October 1-3, 2007October 1-3, 2007 116th Annual CMAS Meeting6th Annual CMAS Meeting

Comparison of CMAQ and CAMx Comparison of CMAQ and CAMx for an Annual Simulation over the for an Annual Simulation over the

South Coast Air BasinSouth Coast Air Basin Jin LuJin Lu11, Kathleen Fahey, Kathleen Fahey11, Bruce Jackson, Bruce Jackson11, Paul , Paul

LivingstoneLivingstone11, Daniel Chau, Daniel Chau11, and Ajith Kaduwela, and Ajith Kaduwela1,21,2

1Planning and Technical Support DivisionAir Resources Board

California Environmental Protection Agency

2Department of Land, Air and Water ResourcesUniversity of California at Davis

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 22

DISCLAIMERDISCLAIMER

The content of this presentation does not necessarily reflect the views and policies of the California Air Resources Board, nor does mention of trade names or commercial products constitute endorsement or recommendation for use.

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 33

South Coast Air BasinSouth Coast Air Basin

One of the most polluted regions and the most studied region in the world.

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 44

South Coast Air BasinSouth Coast Air Basin

116×80 5 km2 grid cells

San LuisObispo Co.

Santa

Barbara Co.

Kern Co.

VenturaCo.

Orange Co.

ImperialCo.

LosAngeles Co.

Riverside Co.

San Diego

Co.

San

Bernardino

Co.

Northern Mexico

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 55

Simulation Setup for 2005Simulation Setup for 2005• MM5 (South Coast AQMD)• Emissions – CARB/SCAQMD• Measurements:

Rubidoux – Riverside (RUB) and Los Angeles (LA)– PM species concentrations: STN database

PM2.5, NO3, SO4, NH4, OC, EC– Hourly gaseous species concentrations: AIRS/AQS

databaseO3, CO, NO, NO2, SO2

• Model predictions– PM - 9 grid cells 24-hr average concentration– O3, - 1 grid cell hourly Meteorology

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 66

Model ConfigurationModel Configuration

Options CMAQ CAMxHorizontal advection PPM BOTT

Vertical advection PPM BOTT

Vertical diffusion Eddy diffusion Eddy diffusion

Mini. vert. diffusivity 0.1 m2s-1 0.1 m2s-1

Gas phase chemistry CB4 CB4

Gas phase chem. solver EBI CMC

Aqueous chemistry RADM RADM

Aerosol chemistry AE3 EFC

Aerosol representation 3 lognormal modes 2 sections

Vertical layers 9 8

Spin-up days 7 8

CMAQ Version 4.6 vs. CAMx Version 4.2

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 77

Species Mapping: CMAQ vs. CAMxSpecies Mapping: CMAQ vs. CAMx

Compound CMAQ CAMxNO3 ANO3I + ANO3J NO3F

SO4 ASO4I + ASO4J SO4F

NH4 ANH4I + ANH4J NH4F

EC AECI + AECJ ECF

OC AORGAI + AORGAJ + AORGPAI + AORGPAJ + AORGBI + AORGBJ

SOA1 + SOA2 + SOA3 + SOA4 + SOA5 + POMF

PM2.5 ANO3I + ANO3J + ASO4I + ASO4J + ANH4I + ANH4J + AORGAI + AORGAJ + AORGPAI + AORGPAJ + AORGBI + AORGBJ + AECI + AECJ + A25J

NO3F + SO4F + NH4F + ECF + SOA1 + SOA2 + SOA3 + SOA4 + SOA5 + POMF + OTRF

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 88

Outline of Comparison and AnalysisOutline of Comparison and Analysis• Scatter Plots

– CMAQ vs. Observation– CAMx vs. Observation– CMAQ vs. CAMx– PM species, annual plots at Rubidoux– O3, quarterly plots at Rubidoux

• Box Plots– PM species, weekly – Gas species, 24-hour

• StatisticsFBIAS - Fractional Bias FERR Fractional ErrorNMB - Normalized Mean Bias NME Normalized Mean Error

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 99

PMPM2.52.5 Scatter Plots Scatter Plots

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 1010

PMPM2.52.5 Box Plots Box Plots

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 1111

PM Species Scatter PlotsPM Species Scatter Plots

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 1212

OO33 Scatter Plots Scatter Plots

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 1313

OO33 Box Plots Box Plots

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 1414

Statistics – Gaseous SpeciesStatistics – Gaseous Species

O3 CO NO NO2 SO2

% CMAQ CAMx CMAQ CAMx CMAQ CAMx CMAQ CAMx CMAQ CAMx

RUB

FBIAS -27.7 -42.0 2.3 2.1 -116.6 -110.4 13.9 -6.6 -43.2 -44.9

FERR 33.5 43.7 53.9 55.8 129.7 122.5 50.8 49.0 126.8 128.1

NMB -23.0 -31.4 -17.2 -18.8 -77.3 -78.0 26.4 0.7 -68.1 -69.4

NME 28.1 33.0 51.7 53.8 83.5 82.7 55.6 46.6 81.9 82.4

LA

FBIAS -38.7 -61.8 27.3 42.7 -36.3 -5.6 45.1 38.0 58.6 57.9

FERR 45.6 63.2 53.7 63.3 84.0 81.0 51.0 47.6 135.4 135.7

NMB -26.2 -41.0 10.5 32.4 -36.3 -15.5 57.1 43.3 0.6 -2.3

NME 33.1 42.6 53.1 67.0 70.2 75.0 65.0 55.5 114.1 114.5

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 1515

Statistics – PM speciesStatistics – PM species

PM2.5 SO4 NO3 NH4 OC EC

% CMAQ CAMx CMAQ CAMx CMAQ CAMx CMAQ CAMx CMAQ CAMx CMAQ CAMx

RUB

FBIAS 12.4 32.4 -19.7 23.9 31.4 14.6 18.5 50.5 -52.9 -25.3 -10.2 -4.6

FERR 45.9 54.0 63.4 68.7 56.2 60.5 58.5 70.8 75.0 59.0 50.1 51.4

NMB -1.5 20.9 -17.9 22.2 7.5 -12.9 -7.2 36.0 -48.8 -34.5 -25.5 -23.9

NME 39.5 49.7 48.1 61.5 46.2 49.3 44.7 63.4 57.1 49.3 45.4 46.9

LA

FBIAS 26.0 45.6 3.9 44.0 15.1 -2.2 21.6 47.2 -29.6 1.3 -30.3 18.0

FERR 40.4 52.9 60.0 74.4 63.3 69.5 63.4 71.3 53.7 38.0 54.2 35.8

NMB 19.7 47.4 -16.3 25.8 -12.6 -29.2 -11.6 22.2 -33.0 -10.8 -33.4 7.0

NME 44.4 62.8 58.5 70.4 54.9 60.3 51.5 61.3 42.5 33.2 42.8 31.4

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 1616

Statistics at RubidouxStatistics at Rubidoux RUB soccer plot

0

20

40

60

80

100

-100 -50 0 50 100

NMB

NM

E

cmaq_o3cmaq_cocmaq_nocmaq_no2cmaq_so2cmaq_pm25cmaq_no3cmaq_so4cmaq_nh4cmaq_occmaq_eccamx_o3camx_cocamx_nocamx_no2camx_so2camx_pm25camx_no3camx_so4camx_nh4camx_occamx_ec

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 1717

ConclusionsConclusions

• Annual simulations over South Coast Air Basin performed with both CMAQ and CAMx

• CMAQ and CAMx model performance evaluated at Roubidoux and Los Angeles

• Model predictions correlate well with each other• Model predictions correlate reasonably well with

observations

October 1-3, 2007October 1-3, 2007 6th Annual CMAS Meeting6th Annual CMAS Meeting 1818

International Aerosol Modeling International Aerosol Modeling Algorithms (I AM A) ConferenceAlgorithms (I AM A) Conference

• Bi-annual in-depth look at the inner-workings of aerosol models (complements the annual CMAS Meeting)

• Wednesday, December 5th – Friday, December 7th, 2007 (just before the AGU meeting in San Francisco)

• University of California at Davis• Impressive list of U.S. and international invited

speakers• Information at http://airquality.ucdavis.edu/