Notes for running WRF with the Aerosol-aware Thompson Scheme (mp_physics = 28)
(See updated info for V3.9 below)
When the namelist variable, 'use_aero_icbc' is set to false, the Thompson & Eidhammer (2014) scheme will assume all model horizontal grid points have the same vertical profiles of water nucleating aerosols (CCN, also known as number of water-friendly aerosols, NWFA) and ice nucleating aerosols (IN, also known as number of ice-friendly aerosols, NIFA). These profiles are controlled by parameter settings of variables at the top of 'phys/module_mp_thompson.F': naCCN0 (300 per cubic centimeter) is the near-surface value of CCN and naCCN1 (50 per cubic centimeter) is the free tropospheric value of CCN. A set of similar variables are used for IN. The vertical profile is terrain height dependent in a manner that was designed to fit the Continental U.S. in which the near-surface value is found to exist within an idealized boundary layer of approximately 200 to 1000 meters depending on starting elevation. The formulation was designed not to follow either height above ground nor height above sea level, because the so-called boundary-layer height is not the same at the top of high mountains like Pikes Peak as it is in Denver or Houston. In effect, the formulation tries to account for a very thin idealized boundary layer height of tens of meters in high terrain above 2500 meters but closer to 1000 meters for grid points near sea level. An exponential decay of aerosol number from the higher numerical value in the boundary layer to the lower free tropospheric number is used to complete the vertical profile. These settings are done once at model initial time (inside subroutine thompson_init) regardless of land versus ocean or other potential geographic information. A future version could incorporate marine versus continental differences.
During model integration, the NWFA and NIFA variables are advected and diffused exactly as other scalars (e.g., cloud ice number concentration), and a zero-gradient lateral boundary condition also follows the other scalars. A fake surface aerosol emissions/flux/tendency is computed as a 2D field (computed in subroutine thompson_init and held in variable called nwfa2d) based on horizontal grid spacing and starting aerosol number concentration for the NWFA variable. No surface emission tendency is applied for NIFA as of this writing (April 2014). The 2D tendency field is added each time step to the first model vertical level NWFA value. Future versions are expected to use more appropriate aerosol emission inventories or other available data.
When the namelist variable, use_aero_icbc is true, the Thompson & Eidhammer (2014) scheme uses an auxiliary aerosol climatology file placed into WRF through the WPS program. Aerosol number concentrations were derived from multi-year (2001-2007) global model simulations (Colarco, 2010) in which particles and their precursors are emitted by natural and anthropogenic sources and are explicitly modeled with multiple size bins for multiple species of aerosols by the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model (Ginoux et al. 2001). The aerosol input data we used included mass mixing ratios of sulfates, sea salts, organic carbon, dust, and black carbon from the 7-year simulation with 0.5-degree longitude by 1.25-degree latitude spacing. We transformed these data into our simplified aerosol treatment by accumulating dust mass larger than 0.5 microns into the ice nucleating, non-hygroscopic mineral dust mode, NIFA, and combining all other species besides black carbon as an internally-mixed cloud droplet nucleating, hygroscopic CCN mode, NWFA. Input mass mixing ratio data were converted to final number concentrations by assuming log-normal distributions with characteristic diameters and geometric standard deviations taken from Chin et al. (2002; Table 2).
Update in V3.9: Capability is added to WPS and real to take in this monthly averaged native model level data. The new data file, QNWFA_QNIFA_SIGMA_MONTHLY.dat, is available for download and can be found on the Yellowstone/Cheyenne supercomputer at /glade/p/work/wrfhelp/WPS_files.
In order to use the data, you will need to place the file inside the WPS/ directory. Before running metgrid.exe, you will need to add the following line in the &metgrid section of the namelist.wps file:
constants_name = 'QNWFA_QNIFA_SIGMA_MONTHLY.dat'
You can then run metgrid.exe. (Note that an updated METGRID.TBL is required). Twelve monthly data for two variables W_WIF and I_WIF, as well as pressure field, P_WIF will be produced in metgrid output file.
To use the data in program real.exe and model wrf.exe, add these namelists:
wif_input_opt = 1
num_wif_levels = 30
mp_physics = 28, 28, 28,
use_aero_icbc = .true.
Questions/concerns can be directed to Greg Thompson (gthompsn AT ucar.edu) or Trude Eidhammer (trude AT ucar.edu).
Thompson, G. and T. Eidhammer, 2014: A study of aerosol impacts on clouds and precipitation development in a large winter cyclone. J. Atmos. Sci., DOI: http://dx.doi.org/10.1175/JAS-D-13-0305.1.