Interactive comment on “ Assessing the CAM 5 physics suite in the WRF-Chem model : implementation , evaluation , and resolution sensitivity ”

This is an excellent paper summarizing a major accomplishment, which is to port climate model physics into a mesoscale model and make this available to a broader community. The paper documents this process well and will serve as a valuable reference for those wishing to test CAM5 physics over regions and at different resolutions, especially the chemistry aspects. The paper clearly discusses the technical challenges of the work. Overall I recommend this work for publication and include some minor suggestions that may help its clarity.


Introduction
Global climate models (GCMs) have been used to simulate and understand the mean state and the associated variability of climate, long-term trends of past climate (e.g., Chou et al., 2013;Deser et al., 2012), large-scale climatic response to various forcings (e.g., Ganguly et al., 2012a, b;Gent and Danabasoglu, 2011;Gettelman et al., 2012;5 Teng et al., 2012), and future climate under different forcing scenarios (e.g., Kravitz et al., 2011;Meehl et al., 2005Meehl et al., , 2012. However, it has been challenging for GCMs to accurately capture climate variability such as extreme weather events at regional and local scales (Kang et al., 2002;Zhang et al., 2013). This model deficiency can be attributed not only to the coarse resolution of global models, but also to their treatments 10 of physical processes. For example, increasing spatial resolution has been shown to improve simulated climate (Roeckner et al., 2006), precipitation (Giorgi andMarinucci, 1996;Li et al., 2011), and tracer transport (Rind et al., 2007). Yet, some model biases such as the diurnal cycle and spatial pattern of precipitation appear insensitive to spatial resolution, and are attributed to deficiencies in the formulation of physical processes 15 (Iorio et al., 2004;Marshall et al., 1997).
The typical grid size for current GCMs is still rather coarse and ranges from 0.5 to 4 • for the atmosphere (Taylor et al., 2012). For example, the Community Atmosphere Model (CAM) version 5 (CAM5) (Neale et al., 2010), the atmospheric component of the Community Earth System Model (CESM) version 1 (Hurrell et al., 2013) typically uses 20 a grid spacing of 1 to 2 • . GCMs are expected to take advantage of growing computational resources and run at much higher resolutions in the future. However, the applicability of climate model treatments of physical, chemical, and dynamical processes in a high-resolution setting has not been fully tested, and the resolution dependence of model physics and model biases is not well understood. Rapid development and eval- The continuous increase of computing power has enabled the scientific community to run GCMs at higher resolutions and to run regional models both at higher resolutions and over larger domains. Establishing a framework to compare the physics suite implemented in a GCM with other representations implemented in a RCM using systematic and consistent methodology is highly desirable for exploring the strengths and 10 weaknesses of different parameterizations across scales. To take advantage of some of these attributes of RCMs, we have transferred a nearly complete physics suite (including the treatment of deep and shallow convection, cloud microphysics, turbulent boundary layer, aerosols, and fractional clouds) from CAM5 to the WRF model with chemistry (WRF-Chem) (Grell et al., 2005). This modeling framework allows exploration of 15 the parameterization suite at high resolutions with a lower cost than a global model, allows direct comparison of the parameterizations commonly used in cloud/mesoscale models with those used in GCMs, and provide an internally consistent methodology to evaluate various treatments of physics, chemistry, and feedback processes for both types of models. The WRF-Chem model with the CAM5 physics suite can be used in -Perform self-consistent dynamical downscaling simulations by using the same physics in the GCM and RCM so that inconsistencies across the RCM's lateral boundary are greatly reduced -Advance process-level understanding through a systematic comparison between the CAM5 physics and other process representations by utilizing WRF's multiple 5 physics capability -Provide insights into the reformulation of parameterizations towards resolutionawareness needed for higher or variable-resolution next generation GCMs by extracting information of subgrid variability from high-resolution simulations (e.g., to make the parameterizations resolution-aware) In addition to describing how the CAM5 physics suite has been implemented in WRF-Chem, we use this new modeling framework to explore and demonstrate the resolution dependence of simulated cloud properties and aerosol concentrations associated with synoptic conditions that transport anthropogenic and natural aerosols from Asia towards the Arctic. The performance of the CAM5 physics suite is also compared to an- 15 other set of parameterizations available in the WRF-Chem model to determine whether there are significant differences in simulated clouds and aerosols associated with the parameterizations developed by the global and regional modeling communities.

Model implementation
Most of the physics parameterizations from CAM5 (Neale et al., 2010) were trans-20 ferred to WRF. The suite of parameterizations illustrated in Fig. 1 is introduced through "interface routines" that connect the parent model infrastructure to each parameterization. These interface routines serve the purpose of converting the model state and other variables into the form expected by each parameterization. They determine the shape of arrays to send into the parameterization such as individual columns or 3- 25 dimensional volumes, flip the vertical ordering of arrays, calculate derived variables 6163 the CAM5 physics modules were implemented in WRF. These modules were released in April 2013 as part of WRF (and WRF-Chem) version 3.5, and can be downloaded from the WRF model users website at http://www.mmm.ucar.edu/wrf/users/. The configuration of WRF-Chem running with the CAM5 physics suite has passed a regression test that the model produces exact same results when using different number of pro-10 cessors.
The specific CAM5 parameterizations ported to WRF include (1) the diagnosed turbulent kinetic energy based first order K-diffusion moist boundary layer scheme (Bretherton and Park, 2009), (2) the convective inhibition closure based shallow convection scheme (Park and Bretherton, 2009), (3) the consumption of convective available potential energy (CAPE) based Zhang-McFarlane deep convection scheme (Zhang and Mcfarlane, 1995) with modifications to use a dilute plume to calculate convection depth and CAPE (Neale et al., 2008) and convective momentum transport (Richter and Rasch, 2008), (4) the two-moment cloud microphysics scheme Morrison and Gettelman, 2008), and (5) (Emmons et al., 2010) simulation. CAM5 includes a cloud macrophysics scheme  that treats fractional cloudiness and condensation/evaporation rates within stratiform clouds. The cloud parameterization in WRF typically do not account for fractional clouds, even when 5 used at coarse resolutions that would benefit from such a treatment. The condensation/evaporation rates in WRF are determined by assuming an instantaneous adjustment to saturation. The CAM5 parameterization uses a rather complex treatment of fields produced by other model components in its treatment of these processes. Because the time discretization is quite different between CAM and WRF, and the order 10 of parameterization and atmospheric dynamic updates are also different, it made porting the CAM5 macrophysics parameterization to WRF quite difficult. Hence, we implemented a simplified module for the treatment of cloud macrophysics. For stratiform clouds, we implemented the subgrid scale cloud fraction, condensation, and evaporation using the same triangular probability density function (PDF) formulation as the 15 CAM5 cloud macrophysics scheme. Tunable parameters were set to standard CAM5 values for a 2-degree grid spacing: threshold relative humidity of 88.75 % for low clouds (reduced to 78.75 % for low clouds over land without snow), 80 % for high clouds, and interpolated thresholds for the mid-clouds. The parameterization produces condensation and evaporation for subgrid clouds with changing cloud fraction 20 Rasch and Kristjansson, 1998), even if the cell is sub-saturated. The deep and shallow convective cloud fractions are diagnosed from deep and shallow convective mass fluxes, respectively. The treatment of convective detrainment of liquid and ice cloud condensates (both mass and number) follows CAM5, but the detrainment rates in WRF are applied after the convection scheme is called instead of in a separate macrophysics 25 module. When the CAM5 microphysics scheme is selected, its cloud fraction is also used in the radiative transfer calculation by overwriting the standard WRF cloud fraction values calculated in the radiative transfer interface routine (Hong et al., 1998). This simpler macrophysics parameterization produces clouds that are very similar to CAM5, Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | but the implementation was much easier, and it requires minimal code modifications to switch between allowing a continuous cloud fraction (as in CAM5) and binary cloud fraction (as in the traditional WRF) to explore the consequences of the treatment of subgrid clouds (Gustafson et al., 2013). Changes in cloud fraction between two time steps also results in activation of interstitial aerosol or resuspension of cloud-borne 5 aerosol, respectively (Abdul-Razzak and Liu et al., 2012;Ovtchinnikov and Ghan, 2005), when running WRF-Chem with the MAM3 option. CAM5 physics integration uses the time splitting approach (Williamson, 2002) that parameterizations are called sequentially and the state variables are updated after each parameterization is called. However, WRF updates the model state by calling 10 parameterization simultaneously to return a tendency, and updates the state from the tendency sum (parallel process updates, with the exception of cloud microphysics) within the first stage of its Runge-Kutta time integration scheme. We have evaluated the effect of the difference between the call sequence in CAM5 and WRF for a domain encompassing the Pacific storm track by implementing a new code infrastructure that 15 allows WRF to call parameterizations and apply tendencies sequentially, as in CAM5. The results show that the difference of the simulated meteorology, clouds, and aerosols between the sequential and parallel calculation is very small (< 3 %), because WRF uses very small time steps (e.g., 1 min or less) compared to CAM5 (30 min).
We have embedded the MAM3 as one option of the WRF-Chem, and aerosol pro-20 cesses and their interactions with radiation and clouds follow CAM5 formulations Liu et al., 2012) with a few minor modifications. A slightly different version of the Rapid Radiative Transfer Model for general circulation models (RRTMG) (Iacono et al., 2008;Mlawer et al., 1997) is used, and the existing aerosol optical property module in WRF-Chem (Barnard et al., 2010;Fast et al., 2006)  can be used for configurations of WRF and for WRF-Chem simulations when aerosolcloud interactions are disabled. In this configuration, the cloud parameterization uses prescribed aerosol numbers and masses for cloud droplet nucleation, with the aerosol mass derived from the prescribed aerosol number and size distribution of each mode. Multiple aerosol surface deposition velocity calculations are available in WRF-Chem 5 when running with MAM3, and the CAM5 formulation that includes sedimentation, turbulent settling, and molecular adhesivity at the lowest model layer (Zhang et al., 2001) is one available option. However, CAM5 also calculates sedimentation velocities above the lowest model layer but this is not treated in WRF-Chem. The omission of this process is expected to have very little effect on the lifecycle of fine particles (i.e., aerosols 10 in the Aitken and accumulation modes), but may have some effect on the lifecycle of the largest particles (i.e., coarse-mode aerosols). The full effect of including sedimentation of aerosols in a RCM requires further investigation. The MAM aerosol package distinguishes between interstitial and cloud-borne aerosols. In CAM5, advection of cloud-borne aerosols is neglected (Ghan and Easter, 15 2006). This may produce some error at higher resolution so we advect both interstitial and cloud-borne aerosols in WRF-Chem.
The convective relaxation time scale (τ) in the deep convection scheme is usually set to be 3600 s in CAM5. With increasing spatial resolution, the model dynamics can produce CAPE at a rate that cannot be removed by the convection with a large τ, 20 producing fallacious "grid-scale storms" (Williamson, 2013). To address this issue, we have introduced a simple formula in WRF version 3.5 to determine τ as a function of grid spacing with a lower bound of 600 s: including τ, at the same values as used in CAM5 (e.g., τ = 3600 s) for the purpose of exploring the resolution dependence of the parameterizations.

From CAM5 to WRF-Chem: the dynamical downscaling procedure
We have developed a dynamical downscaling procedure for this study that (1) minimizes inconsistencies between the parameterizations in the global model CAM5 and 5 the regional model WRF-Chem; (2) facilitates a comparison of simulations between CAM5 and WRF-Chem with CAM5 physics; and (3) produces simulations which agree closely with observed meteorological events. We run CAM5 at a low climate model resolution of 1.9 by 2.5 • (nominally, "2-degree") grid spacing with 56 levels in the vertical as an "offline model" Ma et al., 2013b;Rasch et al., 1997), 10 where the model's winds, temperature, and pressure fields are constrained to agree with time interpolated fields from the ERA-Interim reanalysis (Dee et al., 2011) and there is a "wind-mass adjustment" made to the wind fields to make them consistent with the time evolution of the surface pressure field. The offline methodology has been routinely used for studying the atmospheric tracer transport problems (e.g., Ginoux 15 et al., 2001;Jacob et al., 1997;Lawrence et al., 1999;Liu et al., 2009;Ma et al., 2013a). In the offline CAM5 configuration, water substances and aerosols are allowed to evolve freely according to the CAM5 parameterization suite. The model simulation is archived at 6 h intervals. Meteorological fields (winds, temperature, pressure, and humidity), surface fields (winds, temperature, pressure, latent and sensible heat flux, 20 sea surface temperature, and snow height), and tracers (aerosols, trace gases, water vapor, and liquid and ice cloud condensates) are extracted from the archive and used as the initial and boundary conditions for the WRF-Chem simulations. A MOZART simulation is used to provide the initial and boundary conditions for the additional trace gas species in CBMZ. In these simulations, surface latent and sensible heat fluxes from the Introduction scheme coupled with a land model such as the Community Land Model of the CESM, which is embedded in WRF version 3.5. For the same reason, the effect of topography due to different resolutions is also removed in this study by replacing the standard WRF-Chem topography with the 2-degree CAM5 topography. The vertical sigma coordinate used in WRF is configured to match the sigma-pressure hybrid coordinate used 5 in CAM5 at the initial time step, which has 45 levels from the surface to 20 hPa. All model configurations use the aerosol emissions from the Polar Study using Aircraft, Remote Sensing, Surface Measurements and Models, of Climate, Chemistry, Aerosols, and Transport (POLARCAT) Model Intercomparison Project (POLMIP). This inventory was compiled by Louisa Emmons, and is available for download at ftp://acd.ucar.edu/user/emmons/EMISSIONS/arctas_streets_finn/. It contains anthropogenic emissions from David Streets' inventory for Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) (Jacob et al., 2010), biogenic emissions (Granier et al., 2011), and aerosol emissions from fires (Wiedinmyer et al., 2011). However, it does not include some recently identified aerosol sources such as gas flaring and domestic emissions over the Arctic that might result in under-prediction of aerosol concentration in the Arctic (Stohl et al., 2013). The POLMIP emission inventory includes mass emissions of primary aerosols, aerosol precursor gases, and other trace gases that are used by CBMZ. The emission data processing follows Liu et al. (2012). A small fraction (2.5 %) of the sulfur sources are emitted as primary 20 sulfate. Anthropogenic sulfate emissions are partitioned 85 %/15 % between the accumulation and Aitken mode, and volcanic sulfate emissions are partitioned 50 %/50 % in those modes. Secondary organic aerosol emissions are produced using fixed yields from isoprene, terpene, toluene, and higher molecular weight alkanes and alkenes precursors. Black carbon, primary organic matter, sulfate, and sulfur emissions from 25 fires are vertically distributed in accordance with the spatially and temporally varied vertical profiles described in the AEROsol model interCOMparison project (AeroCom) (Dentener et al., 2006;Textor et al., 2006). Aerosol number emissions are calculated from mass emissions using the assumed sizes listed in Liu et al. (2012) of aerosols and trace gases were first apportioned to the CAM5 grid and then regridded to all WRF-Chem grids so all emissions have the same spatial detail. The MAM3 sea salt and dust emissions are computed online in CAM5 as a function of surface wind speed, and similar parameterizations exist in WRF-Chem as well. However, for the purpose of this paper, we prescribe sea salt and dust emissions in the WRF-Chem 5 with the 6 hourly instantaneous emissions from CAM5 simulation to eliminate the effect of resolution dependency of these aerosol sources.

The aerosol model testbed case
There have been numerous field campaigns in different regions of the world that could be used to evaluate the performance of the CAM5 physics modules within the regional 10 WRF model framework. We first focus on the high-latitudes because it is a vulnerable region for climate change (e.g., Screen and Simmonds, 2010;Serreze et al., 2009), and there are still large uncertainties regarding the role of aerosols in the Arctic (north of 66.5 • N) with documented deficiencies in aerosol transport into the Arctic simulated by GCMs (Lee et al., 2013;Textor et al., 2006;Wang et al., 2013). To better understand 15 the factors controlling changes in atmospheric composition and climate over the Arctic, several field campaigns were conducted in the vicinity of Alaska during the Arctic haze season (Law and Stohl, 2007;Quinn et al., 2007)  The data collected from these Arctic field campaigns have been merged into a single data set for the AMT, as described in Fast et al. (2011). The AMT consists of the 15 WRF-Chem model and a suite of tools that can be used to evaluate the performance of atmospheric process modules via comparison with a wide variety of field measurements. The "Analysis Toolkit" software (available to the community) extracts simulated variables from model history files in a manner compatible with the available measurements using "instrument simulators". These instrument simulators are used to evaluate 20 the simulated meteorology, clouds, and aerosols from both CAM5 and WRF-Chem by extracting simulated quantities in space and in time to match observations collected along aircraft flight tracks. We focus on comparing the model predictions with black carbon (BC) measurements from Single-Particle Soot Photometers (SP2) and aerosol composition measurements from Aerosol Mass Spectrometers (AMS) that were de-25 ployed on the DC-8 and P-3B aircraft (Brock et al., 2011;Cubison et al., 2011;Mc-Naughton et al., 2011;Spackman et al., 2010). Cloud liquid water content and ice water content measured by Convair-580 (Jackson et al., 2012)  In addition to the aircraft measurements, we also use other observational data to evaluate the model simulations. For aerosols, we use the surface BC measurements from a Particle Soot Absorption Photometer (PSAP) at Barrow, Alaska, and the aerosol optical thickness (AOT) measurements from Aerosol Robotic NETwork (AERONET) (Holben et al., 1998) at Barrow and Bonanza Creek, Alaska. Measurements at the 5 DOE's ARM North Slope of Alaska (NSA) site at Barrow, Alaska, are used to evaluate the simulated clouds. The Liu and Illingworth (2000) retrieval is used to evaluate the ice water content (IWC) and the Liao and Sassen (1994) data product is used to evaluate the liquid water content (LWC). For the column integrated ice water path (IWP) and liquid water path (LWP), we use the ARM Best Estimate (ARMBE) data products (Xie 10 et al., 2010), while acknowledging the fact that multiple retrievals show 20-30 % (for LWP) and 10-60 % (for IWP) uncertainty (Zhao et al., 2012). The ARMBE data product is also used to evaluate the total column precipitable water. To assess the fidelity of the simulated precipitation, we used the 1-degree Global Precipitation Climatology Project (GPCP) daily precipitation data product. 15

Resolution dependence of long-range transport of aerosols and aerosolcloud interactions
Aerosol-cloud interactions remain one of the largest uncertainties in climate projections. These interactions occur at subgrid scales in most atmospheric models except for large eddy simulations. Deficiencies in the description of these interactions are be-20 lieved to lead to high estimates of aerosol indirect forcing  and the underestimation of aerosol concentrations in remote regions such as the Arctic (Koch et al., 2009;Lee et al., 2013;Rasch et al., 2000;Shindell et al., 2008;Textor et al., 2006;Wang et al., 2013). The aerosol indirect forcing can also be amplified due to data aggregation from grid-box averages (McComiskey and Feingold, 2012). Increas- 25 ing model resolution to reduce the subgrid inhomogeneity and to better resolve aerosol plumes (Weigum et al., 2012)  the resolution dependence of aerosol concentration and aerosol-cloud interactions in a high latitude region. The offline CAM5 simulation was started on 1 January 2008. Water vapor, condensate, and aerosol fields were allowed to spin up for 3 months, and simulation results from 1 April to 1 May 2008 were analyzed and used for initial and boundary conditions 5 for regional downscaling modeling. The regional modeling domain encompasses northeastern Asia, the northern Pacific Ocean, and northwestern North America (Fig. 1, right panel) so that the primary aerosol transport pathway from Asia to Alaska is included. For all our WRF-Chem simulations, the winds and temperature were nudged towards the offline CAM5 fields (which comes from the ERA-Interim reanalysis) using a time 10 scale of one hour to constrain those fields to be close to the large-scale analyzed fields since our focus was the lifecycle of aerosols and clouds. Simulations were performed using horizontal grid spacings of 160 km, 80 km, 40 km, 20 km, and 10 km (labeled as WRF_160 km, WRF_80 km, WRF_40 km, WRF_20 km, and WRF_10 km, respectively) to explore the behavior of the cloud and aerosol parameterizations as a function of 15 grid spacing. Note that since the lowest resolution of WRF-Chem used in this study (160 km) is still finer than the CAM5 grid (1.9 by 2.5 • ) in the mid-latitudes where most aerosol emissions and wet scavenging take place, there might be differences between these two model simulations. In addition, the differences in dynamical cores and time steps between the two models might also contribute to different results. 20 Figure 3 shows that the simulated meteorology agrees well with the ARCTAS DC-8 aircraft measurements for both the nudged winds and temperature and then freely evolving water vapor mixing ratio at all altitudes around Alaska. There is a 2degree temperature difference between model simulations and aircraft measurements at higher altitudes ( Fig. 3c) that appears to be due to the disagreement between the 25 average temperature within the model layer (typically about 600 m thick at 500 hPa) regridded from the ERA-Interim reanalysis and the instantaneous temperature measurements. As expected, the winds and temperatures from the WRF-Chem simulations are similar to those from CAM5 with no systematic bias associated with resolution. For GMDD 6,2013 Assessing the CAM5 physics suite in the WRF-Chem model clouds, it is challenging to simulate the condensate values exactly when and where they occur along an aircraft flight path given the large spatial and temporal variability of clouds as shown in Fig. 4a and b. Therefore, we have also summarized the performance of simulated in-cloud LWC and IWC along 15 Convair-580 flights in terms of the 95th, 75th, 50th, 25th, and 5th percentiles as a function of the model grid spacing 5 ( Fig. 4c and d). The simulated median instantaneous in-cloud LWC and IWC are generally within the range of the aircraft measurements, but the LWC for the CAM5 and 160 km WRF-Chem simulations is about one order of magnitude lower than observations. This bias is reduced with increasing resolution; the model mean values are within a factor of three for the 10 km simulation. While CAM5 and WRF_160km produce similar features, the remaining differences between these two simulations can be attributed to the differences in their resolutions, time steps, and dynamical cores. As previously mentioned, global models have difficulties simulating aerosols in the Arctic, producing a significant low bias of aerosol concentrations especially near the surface (Shindell et al., 2008). Using the same aircraft transect in Fig. 3, we show that 15 aerosol mass concentrations of black carbon, organic matter, ammonia, and sulfate are low in all model simulations by 1-3 orders of magnitude (Fig. 5). Increasing resolution reduces the bias, but even with the smallest grid spacing of 10 km a significant bias remains. Since these aircraft measurements were taken a few days before an episode with much higher aerosol concentrations that started around 20 April, this result sug-20 gests that the model has a much lower background aerosol concentration. Meanwhile, the variability of the simulated aerosols increases with resolution because the spatially inhomogeneous distribution of aerosols is better resolved in high-resolution simulations. Figure 6 summarizes profiles of BC concentration from a total of eight flights from surface, and about one to two orders of magnitude too low aloft. In general, the higher resolution model simulations show greater variability of BC concentration that results from narrower aerosol plumes with higher concentrations. The concentrations outside of the plume's centers are also lower since the low-resolution simulations tend to dilute the aerosols over a larger region. Decreasing grid spacing from the 160 km to 10 km 5 grid spacing reduces the low bias of BC at the surface by about a factor of 5. Figure 7 compares simulated surface BC concentrations at Barrow with PSAP measurements of equivalent BC (EBC) (Sharma et al., 2006) at the surface with the specific attenuation of 10 m 2 g −1 (S. Sharma, personal communication, 2012). The monthly mean EBC concentration is about two to three orders of magnitude higher than model simula-10 tions (Fig. 7a), and the model consistently under-predict the surface BC concentration by about three to five orders of magnitude before the episode of high BC concentration that starts around 20 April, and about one to three orders of magnitude when the episode starts (Fig. 7b). This model bias reduces monotonically with increasing resolution. These results show that the high-resolution simulations are able to deliver BC 15 from the source regions over Asia to Barrow when the high BC concentration episode took place, even though all resolutions fail to produce the background BC concentration at the observed level. The low background BC might be attributed to either model deficiencies in transporting aerosols into the Arctic (Ma et al., 2013b;Wang et al., 2013), the underestimation or omission of aerosol sources (Stohl et al., 2013;Wang et al., 20 2011b), or a combination of both. The uncertainty associated with the emission inventory is considered to be a factor of 2 or more (Bond et al., 2004;Ramanathan and Carmichael, 2008), and its effect on the simulation requires further investigations. Figure 8 shows snapshots of the simulated distribution of cloud fraction and BC at 00:00 UTC on 20 April. It illustrates the effect of increasing horizontal resolution. Al-25 though the domain averaged cloud fraction from simulations of different resolutions is about the same (∼ 0.2) (Fig. 8a-f), the surface BC over Barrow increases by over one order of magnitude with resolution, from 3.7 ng kg −1 for the 160 km grid-spacing simulation to 51.6 ng kg −1 for the 10 km grid-spacing simulation ( Fig. 8h and l) the 160 km grid cell over Barrow, the instantaneous maximum BC concentration increases by about a factor of 30 (from 3.7 ng kg −1 to 95.6 ng kg −1 ), suggesting that the spatial variability of BC also increases with resolution, consistent with previous studies (Gustafson et al., 2011;Qian et al., 2010). These results suggest that the unrealistic assumption of homogeneous aerosol distribution within each grid cell in CAM5 might 5 be at least partly responsible for the underestimation of aerosol transport into the Arctic, and some of the bias that remains in the high resolution simulations are likely due to the low values of BC present in the boundary conditions for the regional simulations delivered to the region by the low resolution CAM simulations. Figure 8g-l shows that the large difference in surface BC concentration at Barrow is related to the filamentary structure of aerosol plumes, which can only be resolved in the high-resolution simulations, and the transport of BC from northeast Asia to Barrow is sensitive to this flow feature. In high-resolution simulations, a concentrated aerosol plume located in the gap between the clouds of two frontal systems is evident ( Fig. 8f and l). Since clouds and the associated precipitation play a critical role remov- 15 ing aerosols in CAM5 (Liu et al., 2012;Wang et al., 2013) and the poleward transport of aerosols in these latitudes depends on features driven by these eddies (Ma et al., 2013b), this cloudless transport pathway between two mesoscale eddies facilitates the transport of aerosols into the Arctic. In low-resolution simulations, the coarse grids are unable to resolve this filamentary features associated with mesoscale eddies; hence, 20 aerosols are subject to wet removal within clouds along the path. All six aerosol species in MAM3 show behavior similar to BC ( Fig. 9) with more aerosols delivered to the Arctic (north of 66.5 • N) with increasing resolution and the effect gradually accumulates over time. The transport of sea salt aerosols into the Arctic shows the largest sensitivity to resolution for this domain. This is likely due to the fact 25 that both clouds and the resolution-dependent cloud-free pathways are more frequent over ocean than land in the simulations, so the impact of the cloud-free pathways is greatest for aerosols produced over ocean. Another possible explanation is that the sea salt aerosols have a source that is not remote in the regional simulations, so it is not GMDD Introduction

Tables Figures
Back Close

Full Screen / Esc
Printer-friendly Version

Interactive Discussion
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | dependent on the boundary conditions. Although sea salt emission would generally be sensitive to surface wind speeds, and thus could have a resolution dependence (higher wind speeds being resolved with smaller grid spacings), this effect is not present in our simulation because the sea salt emissions have been prescribed from CAM5 for all the WRF-Chem simulations. 5 As aerosol and cloud distributions and loadings evolve, the aerosol-cloud interactions may change accordingly. One convenient means of characterizing the relationship between aerosols and clouds is "cloud susceptibility to aerosols", defined as the fractional increase of cloud liquid water path to fractional increase of aerosol optical thickness (e.g., Quaas et al., 2009;Wang et al., 2011a), often expressed as the slope of a measure of fractional change in a cloud property related to its radiative effect (in this case using liquid water column burden) to fractional change in aerosol amount (in this case using the aerosol optical thickness). CAM5 has a very high value of susceptibility , and part of this behavior can be attributed to parametric uncertainty (Liu et al., 2013). Figure 10 displays the susceptibility as a function of hor- 15 izontal resolution. The cloud susceptibility to aerosols decreases by about 45 % in the highest-resolution simulation. The reduction can be largely explained by the fact that aerosol plumes and clouds are less collocated in high-resolution simulations as shown in Fig. 8. Clouds and precipitation occur in narrower spatial regions, and the occurrence of clouds and precipitation in any given model column is less frequent in many columns, 20 producing weaker aerosol-cloud interactions. Other physical mechanisms may also be sensitive to resolution (e.g. autoconversion and accretion) and affect cloud susceptibility and aerosol indirect forcing. These issues are being explored and will form the basis of a later study. 25 The previous section demonstrated the behavior of aerosols and clouds in the CAM5 physics suite by tightly constraining the meteorological component of the simulations to observations through nudging. In this section, we use "freely evolving" simulations of the WRF-Chem model using the CAM5 physics suite (labeled as "WRF_CAM5") and a commonly used parameterization suite within the WRF-Chem mesoscale modeling community (labeled as "WRF_MESO") to expose differences in the two different physics packages. These simulations indicate where the parameterizations used by the 5 global and mesoscale modeling communities are significantly different, and provide an opportunity to quantify the performance of alternative process modules before they are incorporated and used by global climate models. Such comparison could be expanded to provide insights into the model's structural uncertainty (deficiencies associated with specific treatments of various physical and chemical processes in the model parame-10 terizations) and help to identify the components that produce the largest differences in simulations by changing process representations one at a time. For brevity, we restrict our analysis to a demonstration of capability rather than quantifying effects of each process in the two different parameterization suites. The WRF_MESO configuration includes a double moment cloud microphysics 15 scheme , the Grell and Devenyi (2002) cumulus ensemble parameterization, an operational boundary layer scheme for the Eta model (Janjic, 2002) which employs a turbulent kinetic energy formulation (Mellor and Yamada, 1982), and the 8 size bin version of Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) aerosol model (Zaveri et al., 2008). The WRF_MESO configu-20 ration does not treat aerosol wet removal and vertical transport by subgrid convective clouds, but does represent wet removal when the cell averaged values indicate clouds are present, and aerosols undergo vertical transport by resolved winds and turbulent processes. Our analysis later in this section shows that these differences do not contribute to large differences in the simulations over this high-latitude domain. domain that encompasses Alaska and its vicinity (area enclosed by purple lines in Fig. 12f), approximately 3000 km by 3000 km in total area. All WRF simulations use the global CAM5 results for the initial and boundary conditions for model state (meteorology, cloud condensate, aerosols, and trace gases). The same emission inventory is used to drive all model simulations. The MAM3 aerosol concentrations and emissions 5 from CAM5 are mapped to the corresponding MOSAIC aerosol species in the emission interface subroutine as well as the initial and boundary conditions of WRF_MESO simulations. While the WRF-Chem model with the CAM5 physics suite uses prescribed sea salt and dust emissions as well as surface moisture and heat fluxes archived from the offline CAM5 simulation, the WRF_MESO simulations compute these fields as a func-10 tion of surface wind speed. Figure 11 shows that when driven by the same initial and boundary conditions from the CAM5 simulation, all the WRF_CAM5 and WRF_MESO simulations under-predict AOT by one order of magnitude compared to AERONET observations at Barrow and Bonanza Creek. Increasing resolution only makes very modest improvement, and the 15 WRF_CAM5 and WRF_MESO simulations produce similar results. Vertical profiles are also similar between model configurations, and all simulations are biased-low by one order of magnitude near the surface compared with the observations calculated from the samples obtained from the ARCPAC field campaign (Fig. 12, and see Koch et al., 2009 for details on data processing). 20 Local emissions have been suggested to have a very large contribution efficiency to the aerosol concentrations in the Arctic (Ma et al., 2013a;Stohl et al., 2013). However, they are underestimated in the POLMIP emission inventory, and most aerosols in the Arctic originate from remote sources (Ma et al., 2013a) and enter this WRF domain through the lateral boundary conditions. Since aerosol concentrations in the 25 coarse-resolution CAM5 simulation are significantly underestimated (Ma et al., 2013b;Wang et al., 2013) due to a combination of strong wet scavenging over the ocean and weaker eddy transport in the low resolution CAM5 simulation, and underestimation of emissions of aerosols and their precursors, the regional WRF simulations receives low aerosol inflows from the model's western boundary over Bering Sea and southern boundary over the Pacific Ocean, leading to a low bias for aerosol concentrations and AOT. In addition, since the aerosol plumes reaching this region are already aged, different treatments of aerosol chemistry between MAM3 and MOSAIC make only marginal effects on the simulation. For regional domains containing a major source region where 5 local emissions dominate and aerosols are less aged, differences between MAM3 and MOSAIC can lead to significantly different aerosol simulations (not shown).

GMDD
In spite of the strong influence from the lateral boundaries, there is still one noticeable difference between WRF_CAM5 and WRF_MESO:BC concentration near the surface. It is about a factor of two higher in the WRF_CAM5 simulations (Fig. 12). 10 One possible explanation contributing to this difference is that the CAM5 physics suite takes into account the resuspension of aerosol particles from evaporated raindrops, which is omitted in the WRF_MESO physics suite. Note that both the WRF_CAM5 and the WRF_MESO physics suite consider the aerosol resuspension process in decaying clouds (when cloud droplets evaporate). Further investigation is needed to understand 15 all differences between the two parameterizations and their effects on the simulation, and this modeling framework can be useful for such process-level studies.
The simulated total column water vapor and precipitation rates in all model simulations are generally in good agreements with observations with high correlations and small biases (Figs. 13 and 14). All model simulations appear to have a bias toward too 20 little light precipitation. The correlation between model simulations and GPCP precipitation rate is very good for rain rates larger than 0.4 mm day −1 and poor below that. Increasing resolution generally has little effect on the daily precipitation in this highlatitude domain that is largely over land. However, resolution dependence is expected for instantaneous precipitation rates and precipitation over regions where synoptic and 25 mesoscale features of moisture convergence are better resolved at higher resolution. Figure 15 shows that the WRF_CAM5 simulations produce higher in-cloud liquid water condensate amount and much more realistic (compared with ARM observations at Barrow) frequency of occurrence of liquid cloud than the WRF_MESO simulations,

Tables Figures
Back Close

Full Screen / Esc
Printer-friendly Version

Interactive Discussion
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | despite the fact that the cloud statistics within cloudy cells in all model simulations are qualitatively similar and in good agreement with observations (Fig. 16). One possible explanation for this difference is that the CAM5 physics suite considers subgrid clouds, which allows clouds to form within a fraction of the grid box even if the whole grid box is not saturated. Gustafson et al. (2013) reached a similar conclusion that omitting the treatment of subgrid cloud in WRF running with CAM5 physics suite for typical mesoscale resolutions (grid spacing ranging from 4 to 32 km) results in a reduction of liquid and ice condensate amount. In addition, Fig. 15 also shows that the WRF_CAM5 simulations produce low liquid water condensate above 1.5 km compared with observations, and the bias is improved with increasing resolution. In contrast, the production can be shifted from autoconversion to accretion as the model changes its rainwater treatment from diagnostic rain (as used in the WRF_CAM5 cloud microphysics) to prognostic rain (as used in the WRF_MESO cloud microphysics), which reduces the aerosol effect on clouds (Posselt and Lohmann, 2009). Quantification of the effect of each process in the two parameterizations is beyond the scope of this 20 study and should be addressed by an independent study. Figure 17 shows that all model simulations produce lower ice water condensate amount and higher frequency of occurrence of ice clouds compared to observations, and the WRF_CAM5 simulations agree better with observations than the WRF_MESO simulations. The in-cloud ice water path distributions between model simulations and 25 observations are qualitatively similar, but the WRF_CAM5 simulations show more events of high ice water path whereas the WRF_MESO simulations have more low ice water path events (Fig. 18). In addition, the observations show a mean liquid-tototal water ratio of about 29.8 % in mixed phase clouds, higher than the WRF_CAM5 by computing both the evaporation of liquid condensates and the deposition of water vapor in mixed phase clouds using cell averaged quantities, the Morrison and Gettelman (2008) parameterization computes the deposition rate by converting the subgrid in-cloud liquid to ice condensates directly. This difference can lead to different ice deposition rates. Insights into the improvement on the ice cloud parameterization could be gained from further investigation. In summary, this modeling framework can be useful to probe effects of different parameterizations. We demonstrated that the CAM5 physics suite in WRF-Chem produces a realistic regional meteorology and aerosol distributions similar to simulations using a common set of the WRF parameterizations. Some improvements in the BC con-15 centration near the surface, the frequency of occurrence of liquid water content, and the ice water condensate amount were produced by the CAM5 physics suite. However, the underestimation of aerosols in the global model remains in the RCM. The low bias is insensitive to model resolutions and physics suites for this particular domain and the simulated period of time in this study, suggesting that much of the underestimation is 20 caused by the low estimates of local emissions in this emission inventory and the low concentrations of BC entering the region from the boundary condition values provided by CAM5.

Concluding remarks
The CAM5 physics parameterization suite has been recently ported to the WRF-Chem 25 model. A downscaling modeling framework has been developed that minimizes inconsistencies between the global and regional models. This allows us to run the CAM5 physics suite over a range of scales. We use this downscaling modeling framework to evaluate the CAM5 physics at high resolution, and to directly compare model predictions with high-resolution field campaign data. The suite was released as part of WRF (and WRF-Chem) version 3.5 in April 2013 for broader use by both the WRF and CAM communities for various research objectives. This manuscript describes how the CAM5 5 parameterizations were implemented in WRF-Chem and provides an initial evaluation. We ported the CAM5 physics suite in WRF using interface routines to minimize changes to the CAM5 codes (both formulation and programming). This approach (1) eases implementation of future updates; (2) minimizes the likelihood of changing the behavior of the parameterizations because of implementation issues; and (3) makes the CAM5 physics suite behave similarly in WRF and CAM when running at similar resolutions. Necessary minor modifications appear to have only marginal effects on the simulation. We extended the capabilities of WRF by adding parameterizations that are designed for larger space and time scales, and provided a means to compare simulations using different parameterizations within the same modeling framework to 15 better assess the effects of different model treatments of aerosols, clouds, convections, and aerosol-cloud interactions. The new parameterization suite supports a consistent treatment of subgrid scale clouds, a feature that was previously neglected in WRF and WRF-Chem.
We demonstrated the use of this downscaling modeling framework by exploring the 20 resolution dependence of the CAM5 physics suite for one testbed case. Some model biases in CAM5 (e.g., low aerosol concentration near the surface in the Arctic and high cloud susceptibility to aerosols) were found to reduce with increasing horizontal resolution without any modification to the model physics. While cloud distributions become more realistic at higher resolution as expected, the domain averages of some cloud 25 properties (liquid and ice water path, precipitation rate, etc.) were not very sensitive to varying horizontal grid spacing using the CAM5 physics suite. Our analysis shows that the collocation of aerosols and clouds is very different at different model resolutions, resulting in differences in long-range aerosol transport and cloud susceptibility to aerosols. At higher-resolutions, filamentary aerosol transport pathways evolve according to the circulation associated with the resolved mesoscale eddies that are not present at low resolution. This feature leads to a significant increase of aerosol concentration over the Arctic due to greater eddy transport and weaker wet scavenging (resulting from changes in the frequency of occurrence of clouds, precipitation, and the collocation of clouds and aerosols) during transport. We also show that the WRF-Chem model running with the CAM5 physics suite produces realistic meteorological conditions over the high-latitude domain, and yields higher liquid and ice water condensate simulations and higher BC concentrations near the surface than the example WRF physics suite we examined. We believe this framework can be used to guide future 10 parameterization development that improves the representation of aerosol, cloud, and aerosol-cloud interaction processes and their subgrid variability. However, the present study focuses on only one case and there are many parameterization choices in WRF; other parameterization combinations are likely to perform differently. Additional analyses and more case studies are needed to understand the overall performance of the 15 CAM5 physics compared to other model treatments. Sharma, Hailong Wang, Kai Zhang, Yang Zhang, and Atmos. Chem. Phys., 11, 12049-12064, doi:10.5194/acp-11-12049-2011 Healy Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Ghan, S. J., Liu, X., Easter, R. C., Zaveri, R., Rasch, P. J., Yoon, J. H., and Eaton, B.: Toward a minimal representation of aerosols in climate models: comparative decomposition of aerosol direct, semidirect, and indirect radiative forcing, J. Climate, 25, 6461-6476, doi:10.1175/Jcli-D-11-00650.1, 2012. Ginoux, P., Chin, M., Tegen, I., Prospero, J. M., Holben, B., Dubovik, O., andLin, S. J.: Sources 5 and distributions of dust aerosols simulated with the GOCART model, J. Geophys. Res.-Atmos., 106, 20255-20273, doi:10.1029/2000jd000053, 2001. Giorgi, F. and Bi, X. Q.: A study of internal variability of a regional climate model, J. Geophys.