NCEP Operational Air Quality Forecast Change Log

AQF Operational Run Forecasts

Date	Cycle	Change
		2015-2021 Changes
July 20, 2021	06 UTC	Operational Air Quality Model (AQMv6) Upgrade: NWS Service Change Notice on this upgrade. For the CONUS domain, the CMAQ Kalman Filter Analog (KFAN) bias correction system for fine Particulate Matter (PM2.5) and Ozone will be improved to use a consistent training data set, additional monitor sites, and a unified KFAN bias correction system. Over the CONUS, Alaska, and Hawaii domains, the CMAQ code will be unified with the Environmental Protection Agency (EPA) Version 5.3.1 using an updated 2016 National Emissions Inventory (NEI 2016). All domains will be driven by the NWS Global Forecast System (GFS) v16 meteorology with air quality forecasts extended from 48 to 72 hours for the 06 and 12 UTC run cycles. The Air Resources Laboratory (ARL) Fengsha dust model will be updated with improved threshold velocities by soil type. The Updated Biomass Emissions Land Database (BELD) will be updated to v5 for biogenic emissions processing. Greenness vegetation fraction (GVF) will be updated daily with a 7-day average GVF provided by the National Environmental Satellite, Data, and Information Service (NESDIS) satellite products. Leaf Area Index (LAI) will be updated from a constant value in space and time to a climatological field. The NESDIS Global Biomass Burning Emissions Product eXtended (GBBEPx) product will be used to initialize fire particulate and gas-phase emissions and its associated Fire Radiative Power, which is used to drive fire smoke plume rise. Both GVF and LAI are used for biogenic and deposition processes and should result in better land-related processes. See the NCEP OD Briefing on this implementation package for more details.
Oct 1, 2020	06 UTC	Wild Fire emissions turned off.
July 1, 2020	06 UTC	Wild Fire emissions turned on.
Oct. 1, 2019	06 UTC	All fire emissions turned off.
July 01, 2019	06 UTC	HMS fire emissions turned on (CONUS only)  but without GOES-16 fire points.
May 1,2019	06 UTC	Fire emissions turned back on automatically HMS fire points, assuming all wild-fires (24hrs), diurnal profile Anthropogenic Emissions summary for FY19
April 16, 2019	12 UTC	Point source emissions update
Dec. 18, 2018	06 UTC	AQM V5.1 (CMAQ V5.0.2) Public Release Notice Update to NEI 2014V2 emissions for area and mobile sources  Updated oil and gas sector emissions for current year  Turned of all wild fire and ag fire emissions for winter from HMS Update to the  fine Particulate Matter (PM2.5) bias correction system to use: Consistent training model predictions for the unified Kalman Filter Analog (KFAN) bias correction system. Increase the number of observation sites for model bias correction to an average over 900 monitors. Improvements to forecast extreme events by adding the difference between the current raw model forecast and historical analogs’ mean to the KFAN bias-corrected predictions. The assumption behind this correction is that the CMAQ model, when predicting an extreme event, has some skill beyond that of its historical best analogs. A measure of the exceptionality of such a rare event is given by the difference between the CMAQ forecast and its previous analogs , which is then added to the final bias-corrected forecast. Implementation of a new unified  ozone bias correction system with the same codes and configuration as for PM2.5 except the use of NOx, NOy and ozone as parameters to identify analogs Improvements to address rare ozone exceedance events Quantitative evaluation of AQM V5.1 CONUS upgrade Update the Alaska and Hawaii domain CMAQ code to the same version 5.0.2 used for CONUS : CB05 gas-phase and aero6 aerosol chemistry (155 species) Improved heterogeneous, aqueous, winter-time reactions Improved SOA and coarse mode PM
August 1, 2018	06 UTC	Use of HMS fire activity file for 06 and 12 Z cycle runs.  All fires assumed 24 hour type.  This replaced former NESDIS processing where an analyst would ensure that all fire points were associated with visible smoke plumes and  smoke emission strength would be qualitatively adjusted by adding more fire points (thereby increasing fire area) for large fires.  This implementation has decreased fire smoke emissions for large fires by a small amount. ARL bugzilla 7049 bug fix for dynamically allocated fire points array.  Also initiailized PBOT array (plume base).  Appears fire smoke plume rise is lower and nearer source.
May 1, 2018	12 UTC	CONUS V5.0.2 Point source emissions updated with DOE projections to 2018.  A detailed description of V5.0.2 point, area, moblile, smoke and dust emsissions is found here.
April 17, 2018	06 UTC	NESDIS SAB HMS fire activity processing now filters out all fire points with 2 hours burn duration or less to try and reduce smoke emissions and number of fire points.
April 11, 2018	06 UTC	NESDIS SAB HMS fire activity file (_prelim.txt) updated processing to produce automated fire points over the East of 102W and manual points over West.  More than 500 fire points were created and PREMAQ failed.             Fix implemented to increase number fire points to 10000 .
Sept 2X, 2017	12 UTC	Updated KFAN bias correction training  files to use all CMAQ V5.0.2 for 12 months before present date
Aug 24, 2017		Updated the 8 hr average ozone  file hour range labels to span a full 8 hours (eg: 0-8, rather than 1-8).  The 8hr avg files on AWIPS 227 grid affected:  aqm.t06z.ave_8hr_o3.227.grib2, aqm.t12z.ave_8hr_o3.227.grib2. Bias corrected PM25 daily file day 2 starting date inadvertantly changed to start on day2 rather than start of model run.  Will be fixed with AQM V5.1 implementation.
July 1, 2017	12 UTC	CONUS V5.0.2 Point source emissions updated with DOE projections to 2017.  A detailed description of V5.0.2 point, area, moblile, smoke and dust emsissions is found here.
		CMAQ V5.0.2  IMPLEMENTATION
June 14, 2017	12 UTC	CMAQ V5.0.2  CONUS system promoted to production Release Notification CMAQ V5.0.2 system with increased PM speciation  CB05 gas-phase and aero6 aerosol chemistry (155 species) Improved hererogenous, aqueous, winter-time reactions Improved SOA and coarse mode PM Lateral Boundary Conditions :  Static GEOS-Chem 2006 climatology Dynamic dust boundary conditions from NEMS Global Aerosol Capability (NGAC V2) A detailed description of V5.0.2 emsissions is found here. Point sources to be updated on July 1,2017: 2011 NEI projected to 2017 (with DOE estimates and CEM data) Area Sources : 2011 NEI (Canada 2012) Mobile Sources: 2005 projected to 2011 with CSAPR inventories Biogenics: BEIS 3.14 Updated smoke emissions USFS BlueSky system update to V3.5.1 24 h pre-analysis cycle to coincide with NESDIS HMS one day old fire information Dust:  ARL Fengsha emissions mode Updated Bias Correction to Kalman Filter Analog (KFAN) approach Ozone and PM2.5 files in grib2 format only
March 21, 2017	12 UTC	NAM V4 promoted to operations
Feb. 4, 2016	12 UTC	CMAQ V4.7.2 Implemented (EPA Version number 4.6) Updated emissions using EPA NEI 2011 base year values (Except Mobile emissions). Increased number of vertical levels from 22 to 35. Using Daily NEMS Global Aerosol Capability (NGAC) dust predictions at lateral boundaries. Separate bias corrected PM 2.5 output using NOAA/ESRL bias correction system.  See Djalalova, et al. (2015): Atm. Envir., 108, May,  pgs 76-87 Technical Information Notice
May 1, 2015	12 UTC	Updated emissions using EPA NEI 2011 base year values.
Jan. 29, 2015	12 UTC	V4.6.5 CB05/AERO-4 CMAQ Based on CMAQ V4.6, Byun and Schere (2006).  CB05: 51 species with  156 reactions,  (Sarwar et al., 2013 and 2008). AERO-4: Represent size distribution by log normal distributions of φ, geometric diameter of the particles: Aitken (φ <0.1 μm), accumulation (0.1<φ <2.5μm) and coarse (2.5<φ <10μm). New particle formation: gas conversion and nucleation. (Binkowski and Rosselle, 2003).  Heterogeneous hydrolysis reaction of N2O5 ( key linkage between gas and aerosol phase reactions, includes temperature and humidity impact). ISORROPIA V1.7  gas-particle partitioning:  Partitioning between inorganic gas and particulate species due to thermodynamic equilibrium (v1.7 increased stability). NTR, organic nitrate, biases high and influenced ozone production. It is now photolyzed and removed quicker (Dickerson et al., 2014) by shortening NTR lifetime by a factor of 10. Update NOx emissions with satellite data (Pan, et al., 2014). Modulate fugitive dust emission with a binary switch: whenever there is ice/snow suppress emission. Incorporate NESIDS Hazard Mapping System “observed” wild fire  near real-time smoke emissions  (24-48 hrs old fire ponts) Layer specific time step was added to speed up code. Deposition velocity calculation over trees  were adjusted as VOCs and ozone were overpredicted Increased  the minimum PBL to 50 m where high pollutants were noted esp. over marine boundary layers Implemented  a hybrid lateral boundary condition where GEOS-Chem 2006 climatological LBCs were used up to 7 km. Updated emissions using 2005 base year while extrapolating to 2013 by using: Updated CEM point source data to 2012. Updated DOE projection factors to 2014.
		2014 Changes
August 12, 2014	12 UTC	NAM  upgraded and implemented in production.  All CMAQ runs driven by upgraded NAM.
May 1, 2014	12 UTC	Point, area and mobile emissions  for CB04 chemistry mechanism continue to use recent EPA National Emissions Inventory (NEI,  2005) but projected for the current year.  EGU sources use 2012 CEM data projected for 2014. Details on Alaska Emissions updates Details on Hawaii Emissions updates
		2011-2013 Changes
April 1, 2013	12 UTC	Point, area and mobile source emissions for CONUS CB-IV run updated and extrapolated to current year.
June 1, 2012	12 UTC	Point, area and mobile source missions for CONUS CB-IV run updated and extrapolated to current year.
October 18, 2011	12 UTC	The North American Model (NAM) was upgraded to the National Environmental Modeling System (NEMS) framework Non-hydrostatic Multi-scale Model (NMMB).  The NAM provides hourly meteorological predictions at 12 km to drive CMAQ.
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. Details on Alaska Emissions updates Details on Hawaii Emissions updates
		2009-2010  Changes
July 6, 2010	12 UTC	Corrected powert plant emissions were incorporated; from our evaluations the impact is a minor reduction in predicted ground-level ozone (a few ppb)
September 30, 2010	12 UTC	EMC modified POST3 codes to produce both day 1 and day 2 ozone daily maximum predictions. Implementation of Daily Maximum Day 1 ozone GRIB2 NDGD files
May 1, 2010	12 UTC	Point, area and mobile emissions  for CB04 chemistry mechanism continue to use recent EPA National Emissions Inventory (NEI,  2005) but projected for the current year.  EGU sources use 2008 CEM data projected for 2010.
August 12, 2009	12 UTC	NAM-CMAQ production AQF model capability declared operational on NCEP IBM Power 6 Computing platform
May 01, 2009	12 Z	Point, area and mobile emissions  for CB04 chemistry mechanism upgraded based upon recent EPA National Emissions Inventory (NEI,  2005) and then projected for the current year.  EGU sources will use 2007 CEM data projected for 2009.   For mobile sources, the  EPA Office of  Transportation and Air Quality (OTAQ ) estimates will be used in addition to 2005 NEI v1 emission data sets.
		2008  Changes
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
May 13, 2008	12Z Run	Point, area and mobile emissions  upgraded based upon recent EPA National Emissions Inventory (NEI,  2005) and then projected for the current year.  EGU sources will use 2006 CEM data projected for 2008.   For mobile sources, the  EPA Office of  Transportation and Air Quality (OTAQ ) estimates will be 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.
March 31, 2008	12Z Run	Changes to the NAM WRF-NMM weather forecast model implemented.  Impact on AQ specific weather and ozone predictions are shown here. Expanded the Computational domain by 18% Gravity wave mountain drag parameterization implemented Modified horizontal advection for improved mass conservation Improved Surface Longwave radiation calculation Upgraded GSI data assimilation with NMM based backgroun error covariances Improved assimilation of AIRS radiances, GOES single field of view radiances, MODIS and surface Mesonet winds Use of NAM forecast precip to drive soil moisture calculations outside CONUS where precip observations are lacking Increased terrain smoothing
September 18, 2007	12Z run	Experimental CONUS CMAQ gas-phase only forecast system moved to NCO Operational slot: Expanded domain (5X) to cover Continental U.S. Common WRF-NMM hybrid sigma-P vertical coordinate NAM clear sky radiation prediction used to scale CMAQ photolysis Asymmetric Convective Model used to drive CMAQ moist convection mixing Asymmetric Convective Model-2 used to drive CMAQ PBL mixing Constant, static Ozone Lateral Boundary conditions for all CMAQ levels (turned off use of GFS ozone at CMAQ Lateral boundary top) Updates to CMAQ V4.6 with bug corrections to deposition, optimized advection scheme and plume rise Special California NEI 2002 Non-Road mobile emissions
Eastern U.S. Operational Run Change Log (2005-Sept. 2007)
July 3, 2007	12Z Run	CMAQ Point source emission files updated with 2007 projections.
May 16, 2007	00Z Run	CMAQ static  area emissions updated with 2007 projections.
September 5, 2006	12Z 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 22, 2006	18Z Run	CRISIS Fix:  PREMAQ array bounds corrections implemented. Very small improvement in ozone predictions.
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.  The NAM-Y impact on CMAQ forecasts can be seen from Tsidulko analysis.
August 01, 2006	12Z Run	CRISIS Fix:  Some Point  Source Emissions were still configured for 2005 .  These emissions files were updated for 2006 projections.
July 28, 2006	18Z Run	CRISIS FIX: The WRF post code interpolates WRF outputs from the WRF hybrid vertical  coordinate systems to CMAQ sigma levels and  uses simple interpolation for all  elements except temperature.  For temperature,  to be consistent  with hydrostatic balance, it first computes geopotential heights  for upper and lower boundaries (using the hydrostatic equation) and then computes temperature at desired sigma level using the hydrostatic equation again. This algorithm doesn't work in the first model layer, and in this case temperature is extrapolated from the upper layers.  In case of an inversion, since the algorithm involves several upper layers, the bottom sigma temperature happens to be higher then the first model layer temperature. To fix the problem for   loose coupling,. temperature is simply interpolated  (just like all other elements) for the lowest model layer, while using hydrostatically derived temperature for all the upper layers. Impact of these changes can be seen here.
June 20, 2006	12Z Run	Operational CMAQ Eastern U.S. ozone run driven by operational NAM-WRF. 12 km CMAQ gas w/ reduced CB4 gas  phase chemistry forecasts run twice/day(06 and 12Z) to 48 hours.  WRF outputs interpolated to CMAQ sigma verfical coordinate for loose vertical coupling. 6 hour CMAQ cycling for initial conditions Lateral Boundaries are static except at CMAQ top layer where GFS ozone is specified Above-cloud downward convective mixing turned off Point/Area Source Emissions  for 2006
April 4, 2006	18 Z cycle	Eastern U.S. (3X) run initiated with NAM-WRF parallel run.  WRF outputs interpolated to CMAQ sigma verfical coordinate for loose vertical coupling for 6 hour cycling.  Operational CMAQ NAM-Eta and NAM-WRF 3X domain runs are viewed here
August 31, 2005	00 Z cycle	NE operational forecasts replaced with Eastern U.S. (3X)  Experimental runs. Click here for more info on Eastern U.S. configuration The 3X operational forecasts were the same as the NE U.S. runs except: Larger CMAQ domain: 268x259x22 Lateral Boundaries are static except at CMAQ top layer where GFS ozone is specified Above-cloud downward convective mixing turned off
N.E. U.S. Operational Run Change Log (2004-2005)
June 10, 2005	00 Z Run	A bug was corrected in the interpolation  of ozone from the new GFS T382 to the CMAQ 12 km grid.  This has helped reduce ozone overpredictions found in runs from  May 31-June 10, 2005.
May 31, 2005	12 Z GFS Upgrade	The NCEP Global Forecast System (GFS) was upgraded to T382 (~35 km). In addition, the following implementations were changed: A new land/sea ice physics Reduced vertical diffusion Enhanced mountain blocking assimilation of AIRS satellite data GFS ozone is used to prescribe the CMAQ lateral Boundary Conditions above 6 km.
May 3, 2005	12 Z NAM-12 CMAQ Operational (1X domain)	Upgraded NAM-12.   Upgraded CMAQ  :    PREMAQ updated to use NAM-12 1 km landuse data     Minor CMAQ updates to version 4.4     2002 Point/Area Source Emissions projected for 2005     Updated Mobile 6 source emissions     GFS ozone for CMAQ LBCs above 6 km
2004
September 15, 2004	00 Z	N.E. U.S. Run declared operational
May 1, 2004	12 Z Eta-12 CMAQ Experimental	Implemented Experimental System 2002  Point/Area Source Emissions Projected for 2004 Updated Mobile 6 mobile source emissions 6 hour initial condition cycling Use of GFS ozone predictions for CMAQ upper Lateral Boundary Conditions March 2004 Eta-12 implementation 48 hour forecasts for both 06 and 12 UTC runs
May 5, 2004	12 Z Experimental	Correction to 6 hour Cycling system. Now turned on.
May 17, 2004	12 Z Experimental	Correction to GFS ozone LBC ingest.  Now working.
May 26, 2004	12z Experimental	Sfc ozone Grib files now available on /com/aqm/prod
July 20, 2004	06z	1. Limited  GFS LBC  ozone ingest to above 6 km AGL (CMAQ sigma lvl 18) 2. Specification of  a minimum mixing coefficiant (Kz) in CMAQ pbl routines    *  Allows for lower Kz in rural areas than urban at night :  reduces  nightime O3 overprediction in rural  areas    * Prevents precursor concentrations in urban areas becoming too large at night and      thereby reduces overtitration in urban areas (more mixing in urban areas at nght)

NCEP FVS verification system

Date	System Component	Change
June 1, 2005	Grid2obs VSDB records	VSDB exceedence records (FHO) now being produced for daily maximum 1 and 8 hour ozone predictions.
March 25, 2005	Grid2obs VSDB	VSDB file generation turned on for 2005 Ozone season: Backward averaged convention used for both Airnow bufr and cmaq forecasts L1L2 and FHO statistics now computed (except daily maximums)
2004 FVS VSDB File changes
May 31, 2004	Grid2obs VSDB records	VSDB records (L1,L2) begun for NE domain on 146 grid using forward averaged BUFR observations
June 9-June 20	BUFR generation	BUFR software changed at EPA to produce backward averages.  VSDB records not produced
June 21, 2004	BUFR Generation	BUFR softward changed at EPA to produce foreward averages for NCO data tanks. Bufr tanks and VSDB records for June 9-20 recreated with foreward avging
July 12, 2004	VSDB records	VSDB records are now computed for East US domain on Grid 146, 142 Subdomain records computed for all grids (eg: NEC, SEC,APL...)
July 22, 2004	3X VSDB Records	VSDB records computed for developmental 3x domain (Grid 142) over full domain, and subregions, and over 1 x domain
August 1, 2004	VSDB records	VSDB subdomains updated to new FVS regions.  MDW, GMC and  LMV now replace NMW, SMW (midwest) subregions.
August 1, 2004	FHO VSDB records	Began creating  1h & 8 h VSDB FHO records for events : 1hr ozone > 20, 40, 65, 85, 105, 125 ppb 8h avg ozone > 20, 40, 65, 85 ppb

Summer 2003 North East Domain Forecasts

Date	Model Run	Change
June 25, 2003	12 Z Eta-12 – CMAQ Experimental	Implemented Experimental System with 24 hour ozone initialization cycling. Used Spring 2003 Eta-12 implementation
July 23-24	12 Z Experimental	CMAQ Archives lost
July 31- August 4	12 Z Experimental	Cold Start Runs…No cycling
August 5	12 Z NCO run	Transfer to NCO Cold Start…no cycling
August 5	06Z Experimental	06Z implemented 6 hour cycling (6z only)
August 15	06 Z Experimental	Cold Start
August 20	06 Z NCO 12 Z NCO	Transfer of 6 hour cycling to NCO for 6 and 12 Z runs. Cold Start
September 9	06 Z, 12 Z NCO	Eta-Post:  landuse, soiltype variable packing fixed.  Previously landuse values incorrectly packed to equal 0 (Ocean) or 10.