NCEP CONUS Experimental Run Air Quality Forecast Change Log

 CB05/AERO-4  Experimental Run Forecasts 

Date	Experimental  Model Run	Change
		2011-2013 CONUS EXPERIMENTAL RUNS
April 1, 2013	12 UTC	Point, area and mobile source emissions for CONUS CB05 run updated and extrapolated to current year.
June 1, 2012	12 UTC	Point, area and mobile source emissions for CONUS CB05 run updated and extrapolated to current year.
June 1, 2011	12 UTC	Updated CEM point source data from 2008 to 2009. Updated DOE projection factors from 2010 to 2011.  Changed the number of eletricity marketing module (EMMs) from 13 to 22, consistent with the changes made in DOE models.  Removed average fire emissions from area sources in the experimental domain.
		2010 CONUS EXPERIMENTAL RUNS
July 6, 2010	12 UTC	Corrected power plant emissions were incorporated; from our evaluations the impact is a minor reduction in predicted ground-level ozone (a few ppb)
May 26, 2010	18 UTC	Post-processing upgraded for efficiency.  Now using cmaq2grib codes.  Ozone day 1 Daily Maximum predictions added to NDGD files and display.
May 1, 2010	CONUS 12 UTC	Point, area and mobile emissions were upgraded based upon recent EPA National Emissions Inventory (NEI,  2005) and then projected for the current year. EGU sources use 2008 CEM data projected for 2010.  These emissions are also consistent with the CB05 chemical mechanism.
		2009
August 12, 2009	CONUS 12 UTC	NAM-CMAQ V4.6 production AQF model capability declared operational on NCEP IBM Power 6 Computing platform. A correction to Biogenic emissions cycling from 24 to 6 hour cycling was implemented with the machine switch.
May 1, 2009	12 Z	Point, area and mobile emissions were upgraded based upon recent EPA National Emissions Inventory (NEI,  2005) and then projected for the current year. EGU sources use 2007 CEM data projected for 2009.  These emissions are also consistent with the CB05 chemical mechanism.  For mobile sources, the  EPA Office of  Transportation and Air Quality (OTAQ ) was used in addition to 2005 NEI v1 emission data sets.   A ten day test showing the minor impact of these upgrades are shown here.
		2008
December 16, 2008	12Z	Changes to the NAM  driving meteorological model forecast system: The background for the first (tm12) analysis in each NDAS run is now from the GDAS instead of the previous NDAS run (so-called "partial cycling"). Land states are still fully cycled from the previous NDAS cycle. WRF-NMM Model changes (also implemented into the DGEX): The PBL/turbulance schemes were modified to mix each hydrometeor species in the vertical. To apply vertical diffusion for separate water species, the model was changed so that (a) it can apply vertical diffusion to an arbitrary number of species, (b) the counter gradient option can be applied to some or all of the species if desired, and (c) option to set to zero some or all of the surface fluxes is also made available. In the radiation parameterization, the absorbtion coefficients for water and ice have been doubled to 1600 and 1000, respectively Changes to land-sfc physics: Let the potential evaporation decrease linearly with Bulk Richardson number under stable condition, and weighted by snow coverage. Let the slope of saturated humidity function wrt temperature decrease linearly with snow coverage. Changes to GSI analysis Use latest (1Q 2008) version of the GSI analysis code Assimilate METOP radiance data Assimilate TAMDAR/AMDAR aircraft data New version of Communitity Radiative Transfer Model Use AFWA 1/16 bedient snow depth analysis Use WPS (instead of WRF-SI) codes to process GDAS first guess input files, which are used as a first guess to the first (tm12) GSI analysis in the NDAS
August 19, 2008	12Z Cycle	The  biogenic emission files for the experimental run CB05 chemical mechanism were incorrect and did not include some secondary organic aerosols.  The point source emission files were also incorrect.  These files were corrected for the experimental run mechanism.  Retrospective testing showed that ozone forecasts were decreased by 3-4 ppb overall while developmental PM 2.5 forecast products were increased slightly.
June 12, 2008	06Z Cycle	PM2.5 developmental product predicted from the experimental run was computed incorrectly with updated CB05 emissions.  The primary organic aerosol emission was given in carbonaceous mass instead of aerosol mass.  The error resulted in an underestimate of PM2.5 since the upgraded emissions were implemented on June 10, 2008.
June 10, 2008	06Z Cycle	Specific experimental run system improvements include: The CB05 chemical mechanism with more detailed Volatile Organic Compound (VOC) chemistry was implemented.  A harmonized version of the AERO-4 aerosol chemistry mechanism was turned on. The harmonized NOx night time chemistry is configured with heterogeneous pathways turned off but gas-phase N205 hydrolysis turned on.  This was done to minimize the impact of aerosol chemistry on the ozone foreceast. Sea salt emissions were included in CMAQ. Point, area and mobile emissions were upgraded based upon recent EPA National Emissions Inventory (NEI,  2005) and then projected for the current year. EGU sources use 2006 CEM data projected for 2008.  These emissions are also consistent with the CB05 chemical mechanism.  For mobile sources, the  EPA Office of  Transportation and Air Quality (OTAQ ) was used in addition to 2005 NEI v1 emission data sets.  Use of OTAQ on-road emission estimates is a departure from the temperature dependent regression approach used in previous years. Experimental and Developmental products from these runs are now cycled every 6 hours with forecasts to 48 hours at both 06 and 12 Z cycles. Impacts of these potential changes are shown here Additional Developmental Products: These additional products are  restricted to designated scientists and  state forecasters through the EMC  restricted access AQF web page.
Summer 2007 Experimental Run Changes
September 18, 2007	12 Z run	Experimental CONUS Run moved to operational Status
July 18, 2007	12 Z run	1. An error in the dry deposition velocity calculation was identified and corrected for selected species (NO, NO2, and CO). The error is associated with the mesophyll resistance which was erroneously set to zero. The impact is a reduction in deposition velocity and a resultant increase in predicted concentrations of these species. Higher NOx concentrations then result in higher predicted O3 concentrations. Test simulations confirm this trend with widespread enhancements in predicted O3 in regions with vegetation (primarily in the East U.S.). The max 8-hr O3 increased by 2-5ppb for cases during May and June 2007. The direction of the predicted O3 change was in the correct direction, i.e., reduction in under-estimation in the current version. 2. An error was corrected in the plume-rise calculations to use the correct pressure fields estimated from the NAM Hybrid vertical coordinate. Test cases during June 2007 showed minimal impact.
June 12, 2007	12 Z run	BL mixing scheme in CMAQ has been change from RADM version 2.6 (Byun et al. 1999) using the NAM TKE PBL height to Asymmetric Convective Model-2 (ACM-2) using a Richardson calculation in PREMAQ (Jon Pleim, JAM in press). Minimum eddy diffusivity used in ACM-2 reduced from 0.5 to 0.1 m^2/s
May 29, 2007	18 Z run	Two Summer 2006 changes were not implemented on 5/25/07. These were the ACM convective and NAM clear sky radiation photolysis attenuation and  were implemented today.
May 25, 2007	12 Z run	Ozone has been underpredicted in the Conus experimental run since CMAQ was upgraded.  Therefore, Summer 2006 CMAQ configuration with RADM-2 PBL parameterization and NAM PBL height has been temporarily reinstated. Vertical diffusivity (KZL) minimum has been set to 0.1 m^2/s.
May 16, 2007	00 Z run	Mobile emissions were missing from experimental run since May 1. Therefore, mobile emissions not included in CMAQ predictions from May 1-15, 2007.  Reverted back to using the 2006 Mobile emissions for now.
May 4, 2007	12 Z run	Implemented the Asymmetric Convective Model-2  PBL parameterization into CMAQ.  This uses the flux Richardson  number calculation in PREMAQ to compute PBL height. Vertical diffusivity (KZL) minimum has been set to 0.5 m^2/s.  Static emissions for point sources updated for 2007 projections.   Evaluation of results are shown here.
May 1, 2007	12 Z run	Area source emissions updated for 2007 projection. Updated California Non-Road NEI 2002 inventories
Summer 2006 Experimental Changes
September 5, 2006	18 Z run	PREMAQ met processing calculations of vertical eddy diffusivities were mistakenly set from the 1st 22 layers from the NAM 60 layer predictions rather than for the actual  22 CMAQ layers.  This adversely impacted the vertical advection in the PREMAQ Jacobian calcluations by increasing vertical diffusion in CMAQ.  This was corrected and impact can be evaluated from Lee, Slides 7-17
September 5, 2006	12 Z run	CMAQ now driven by updates to NAM model (additional NAM-Y changes implemented) that include increase to horizontal diffusion and corrections to SST initial conditions.
August 15, 2006	12 Z run	CMAQ now driven by NAM-Y meteorology  which was implemented into NAM slot.  NAM-Y had corrections to vertical and horizontal diffusion and lake SST initialization which resulted in more realistic warmer,dryer boundary layers that increase production of ozone.  Impacts can be evaluated from Lee, Slides 7-17.
Augus 04, 2006	18 Z cycle	CRISIS Fix:  Some Point  Source Emissions were still configured for 2005 .  These emissions files were updated for 2006 projections. CRISIS FIX:   Major errors found in the tightly coupled PREMAQ emissions processor.  Isoprene biogenic emissions were never set, causing strong ozone underpredictions in the tightly coupled 5X CMAQ runs.  Also, some arrays were improperly dimensioned causing array bound errors.  For impact, see Mathur, slides 20-25
July 12, 2006	00Z cycle	Most Emissions updated with 2006 Projections
June 8, 2006	06Z cycle	WRF-CMAQ tightly coupled runs were available late for this case only at 18:40 UTC.  Removal of GFS ozone interpolation scripts(since GFS ozone no longer used for CMAQ) and optimizing job procedures are being tested for quicker turn around.
May 8, 2006	18Z cycle	WRF-CMAQ tightly coupled system transferred to NCO 4x/day  parallel testing and replaces the parallel NAM-Eta CMAQ system. 6 hour cycling turned on for CMAQ initialization. 06 and 12 UTC cycles are run to 48 hours.
April 12, 2006	12Z Cycle	WRF-CMAQ tight vertical coupling Testing CMAQ driven by NAM-WRF parallel run with tight vertical coupling (using the same Sigma-p Pressure hybrid vertical coordinate).  WRF NMM output on all 60 sigma-pressure hybrid layers now read in by PREMAQ and collapsed to 22 CMAQ model layers.  These 22 layers are a subset of the WRF 60 level system. CMAQ ACM convective parameterization Constant Static Lateral Boundary Conditions (O3=100 ppb)  Emissions for 2005 Optimized CMAQ advection scheme addition of diagnostic tracer outputs
SUMMER 2005 Developmental Run Changes
Aug. 20, 2005	12Z Cycle	Cold start.  Aug. 19 warm start rerun but not in time to initialize Aug. 20 cold start. Aug. 17 run GRIB output not created.
July 26, 2005	12Z Cycle	CMAQ CONUS Run  now using ACM  approach to parameterize deep convection and  photolysis attenuation from NAM clear sky radiation predictions and  static lateral boundary conditions (ozone at lbcs=100ppb). This was done to correct overmixing of ozone downward in high terrain.
July 16, 2005	12Z Cycle	CONUS Run  w/S3 convection modified to use static lateral boundary conditions (ozone at lbcs=100ppb). This was done to correct overmixing of ozone downward in high terrain.  Problem may be related to CMAQ mass adjustment scheme and NAM interpolated winds that cause spurious vertical motions near model top.
July 9, 2005	12Z Cycle	CONUS Run modified to replace ACM convective mixing scheme with Standard CMAQ mixing scheme except with downward convective mixing turned off (S3) This was done to correct overmixing of ozone downward in high terrain.
June 29, 2005	12Z Cycle	CONUS run was restarted.  This first run was initialized as a cold start (no atm. chemistry init fields).  However, updated emissions are used.  Successive runs(6/30/05 and on) will use 24 hour cycling to initialize CMAQ atm. chemistry with a previous runs 24 hour forecast.
June 24, 2005	12Z Cycle	NCEP Develop machine crashed.  No CONUS runs available due to this hardware problem. We expect CONUS runs to be available again soon with a cold start initialization (no cycling).
June 17, 2005	12Z  Cycle	GFS Ozone predictions are used to specify the CMAQ lateral boundaries at only the CMAQ model top level (100 mb, level 22).  This appears to help overprediction over the Rocky Mountain areas.
June 01, 2005	12 Z  Cycle	The CONUS AQF (5X) 48 hour Forecasts were begun with the following options: CONUS Emissions, with Special California Inventories Testing the Asymmetric Convective Model (ACM) in CMAQ for pollutant mixing in deep convection (S4 configuration) Use of GFS Ozone at lateral boundaries above 6 km Use of NAM clear-sky radiation to drive CMAQ photolysis rates Maximum cloud fraction is recomputed to be the global maximum cloud fraction from all possible cloud types 24 hour cycling