Weather Effects of Aerosols in the Global Forecast Model

: The weather e ﬀ ects of aerosol types were investigated using well-posed aerosol climatology through the aerosol sensitivity test of thermodynamic and hydrometeor ﬁelds, and the weather forecast performances in July of 2017. The largest aerosol direct radiative forcing (ADRF) in July was due to dust aerosols at the surface and atmosphere, and sulfate at the top of the atmosphere (TOA), respectively. The ADRF of total aerosols had unilateral tendencies in thermodynamic and hydrometeor ﬁelds. The contribution of individual aerosols was linearly additive to those of total aerosols in the heat ﬂuxes, heating rates, humidity, and convective precipitation. However, no such linearity existed in temperature, geopotential height, cloud liquid or ice contents, and large-scale precipitation. Dust was the most inﬂuential forcing agent in July among ﬁve aerosol types due to the largest light-absorption capacity. Such unilateral tendencies of total aerosols and a part of the linearity of individual aerosols were exerted on the weather systems. The veriﬁcation of medium-range forecasts showed that aerosols alleviated the overestimation of surface shortwave (SW) downward ﬂuxes, the negative biases of temperature and geopotential heights at TOA and surface, and the underestimation in light and moderate precipitation. In contrast, they enhanced warm biases at the mid-atmosphere and underestimation in heavy precipitations, particularly negative biases in the intertropical convergence zone (ITCZ). Weather forecast scores including current aerosol information were improved in geopotential height (GPH) of the northern hemisphere (NH); however, they got worse in the temperature and the upper atmosphere GPH of the southern hemisphere (SH), which was mostly due to black carbon (BC) aerosols in the tropical regions. The missing mechanisms such as aerosol–cloud interactions, better aerosol spectral optical properties including mixing states and aging, and the near-real-time (NRT) based aerosol loading data are worthwhile to be tried in the near future for ﬁxing the intrinsic underestimation of precipitation in ITCZ and surface radiative ﬂuxes in the desert and biomass burning area. Eighteen AERONET sites were selected as the location where both AOD and AAOD were available. MODIS is onboard Aqua that is at 1:30 p.m. ascending node, sun-synchronous, and near-polar orbit. The MODIS_AOD indicates global coverage monthly data derived from MODIS level 3 daily mean data of Combined Dark Target and Deep Blue AOD at 550 nm for land and ocean. MISR is onboard Terra. The MISR AOD indicates level 3 Global Aerosol product covering a one-month period such as the optical depth average. precipitation to aerosol types. The changes in convective precipitation (PRECC) and large-scale precipitation (PRECL) relative to aerosol-free conditions for a 10-day forecast starting from 1 July of 2017.


Introduction
Aerosols play essential roles in Earth's Weather and Climate by modifying radiative transfer in the atmosphere. Aerosols directly scatter and absorb the electromagnetic waves, which are called aerosol-radiation interactions and indirectly affect them by changing cloud microphysical properties such as cloud albedo or cloud lifetime, which are called aerosol-cloud interactions. The climate system perturbed by aerosols induces a radiative imbalance at the top of the atmosphere (stratosphere or tropopause), and it reaches re-equilibrium after decadal to more than a hundred years. The aerosol direct forcing is the energy available to impact the climate system with the change in global mean surface temperature. Another perturbation by aerosols is rapidly adjusted after days to weeks without affecting the global mean surface temperature. This instantaneous direct forcing affects atmospheric thermodynamic fields in the weather time scale, which renders the cloud and precipitation changes and has been referred to as the semidirect effects [1][2][3].
Many efforts have been made to study the climate effects of aerosols since the WCP-12 [4] suggested the insightful hypotheses on the climate changes. The IPCC has published the progress on climate changes due to aerosols such as the magnitudes of radiative forcing by total and individual aerosol types along with the distributions of aerosols emission and mass mixing ratios, and physicochemical and optical characteristics of aerosols [1]. The characteristics of aerosol properties have been unveiled not only in the emission source regions but also in transport, mixing, and aging processes. These efforts have made it possible to explore the roles of individual aerosol types in climate. The nonabsorbing aerosols such as sea salts and sulfate are paid attention because they decrease energy budget in the earth-atmospheric systems and have opposite effects on the greenhouse effects [5][6][7]. The absorbing aerosols such as dust and carbonaceous aerosols convert the solar radiation into thermal and moist static energy to affect regional and global dynamics and hydrological cycles [8]. The climate effects are not necessarily identical among aerosol emission sources, aerosol optical depth (AOD), and radiative forcing [1,9,10]. The interhemispheric asymmetry of aerosol emissions results in climate responses beyond aerosol emission regions [1,11]: (i) the shifts of ITCZ by sulfate aerosols and black carbon (BC) [1,10,[12][13][14], (ii) weakening of monsoon circulation by sulfate or SO 2 emissions [15][16][17], and (iii) increases in temperature and precipitation due to reduction in SO 2 emission across the hemisphere while BC emissions affect localized [10,18,19]. Various climate change scenarios are examined regarding the glacial era, industrial revolution, volcano eruption, and future projections. Due to relatively small limitation of computing times and the larger time intervals of integration (~one month), climate models include the various aspects of aerosol properties: aloft emission fields of biomass burning aerosols, the vertical location of aerosol layers associated with clouds, the BC in the snow in the high altitudes or arctic regions.
In contrast, the weather effects of aerosols have been relatively less studied because the weather forecast takes a priority in fast and accurate predictions. Whereas simulating prognostic aerosols requires considerable computing time and expenses and the challengeable assimilation of initial data for aerosols fields [20]. Therefore, aerosol climatology has been alternatively used in NWP models. The main findings of weather effects appear in the local, specific events with the high aerosol concentration area, and tropical regions or monsoon regions where the intense solar radiation is available [21][22][23]. For instance, biomass burning events [24,25] and monsoon circulation [21,[24][25][26][27]. However, the aerosol effects are not evident on a global scale, although some studies detected the large-scale mean flows related to the teleconnection between the tropical region and the extratropical regions [28]. On the other hand, Remy et al. [27] suggested interactions between aerosols and boundary layer meteorology due to feedback of aerosol radiative forcing to aerosol production mechanisms.
In terms of implementing aerosol data into NWP models, several experiments were made. Interactive aerosol models with prognostic aerosols [22,29] were tried, but the results were in noise level despite significant computational expenses. Consideration of direct and indirect radiative forcing of aerosols resulted in not explicit impacts of the indirect forcing because of high uncertainty in cloud microphysical processes [21,30]. Instead, progressing aerosol climatology is more or less promising to the improved weather forecast [26][27][28]31]. In particular, Bozzo et al. [23] constructed aerosol climatology ameliorated with recent observation data and data assimilation and showed significant local impacts of the regional aerosol radiative forcing. However, few studies treated the roles of the individual aerosol types in thermodynamic and hydrometeor fields and their weather effects.
Aerosol climatology was often used as an initial condition in the NWP models. It is not dynamic but static, and thus it should be logical, science-based, and as realistic as possible with the most recent information and technology. Notably, the spectral optical characteristics and atmospheric loading of aerosol types are fundamental to control the sign and magnitudes of their radiative forcing, respectively. Their products are decisive in determining the roles of aerosol types in the weather system. However, some studies used aerosol optical properties out of date, or calculated the products of spectral optical characteristics and loading data of each aerosol type without consistency [32]; one type was based on the chemical composition, and other types were based on the geographical regions.
Atmosphere 2020, 11, 850 3 of 29 In general, the aerosol types of geographic regions are composed of a few aerosol types of chemical composition [33,34].
To overcome the problems mentioned above with aerosol climatology in the global forecast model, we suggest well-posed aerosol climatology data. Well-posed aerosol climatology should be consistent between aerosol loading and single-particle optical properties, and among different fields such as thermodynamic and radiation fields. It leads to not only assess the aerosol-radiation interaction (ARI) but also extend to the various roles of aerosols such as aerosol-cloud interaction (ACI), surface albedo, and aerosol-chemistry interactions. We constructed a well-posed aerosol climatology by incorporating decent spectral optical properties of aerosols into the recent aerosol concentration data from observations or model outputs. Then, we implemented the well-posed aerosol climatology in the Korea Integrated Model (KIM) and evaluated aerosol radiative properties. To examine the behaviors of aerosol types in the thermodynamic and hydrometeor fields through ARI, the sensitivity to individual aerosol types was taken in the NWP model. Finally, we link the roles of aerosol species in thermodynamic and hydrometeor fields to the performance of weather forecasts. This paper will expound as follows. Section 2 describes the configuration of aerosol climatology and the calculation of aerosol optical properties and experimental setups. The radiative properties of modeled aerosols are analyzed in Section 3. Sections 4 and 5 present the roles of aerosols in the global NWP models through sensitivity test and forecast performance, while Section 6 provides conclusions.

Construction of Aerosol Climatology and Implementation of Aerosols into Radiation Module
Aerosol climatology is composed of microphysical, chemical, and optical properties and loading data of aerosol types. It was implemented in the KIM through initial data of aerosol loading and look-up tables of aerosol spectral optical characteristics. Kimaerclim used aerosol optical characteristics embedded in the rapid radiative transfer model for global circulation models (RRTMG) [32]. Aerosol loading data was manipulated with the inverse of monthly AOD at 550 nm from Monitoring Atmospheric Composition and Climate Reanalysis (MACCRA) [35] and the three-dimensional aerosol mass concentration fields of Model for OZone and Related Chemical Tracers (MOZART) [36]. The look-up table (LUT) of aerosol spectral optical characteristics contains the relative ratio of AOD at 550 nm, single scattering albedo (SSA), and asymmetry parameter (ASY) at 14 shortwaves (SW) wavelengths [37]. However, the aerosol climatology of KIM-v3.3 (Kimaerclim) has inconsistency in terms of several points. First of all, aerosol optical characteristics were not consistently matched with aerosol loading data. The aerosol loading data in Kimaerclim are classified by chemical composition such as sea salt, dust, organic matter (OM), BC, and sulfate, while aerosol types of optical properties are classified by geographic regions such as maritime, desert, continental, urban, and volcano. Second, some aerosol optical characteristics are out of date. Besides, the aerosol optical characteristics were treated without consistency. For instance, the extinction coefficients and asymmetry parameters of desert aerosols were calculated with the refractive index of Tanré et al. [37] and the size distribution of McClatchy et al. [38]. However, its single scattering albedo was recently switched to the values calculated with refractive indices of Sinyuk et al. [39] and the size distribution of Drury et al. [40] for the improvement of NWP performances. To keep consistencies of microphysical and chemical data among desert aerosol optical characteristics, this study used the preexisting values of desert aerosols as Kimaerclim. The optical characteristics data of the Kimaerclim were compared with those of the new aerosol climatology (Newaerclim) that this study constructed, in Appendix A ( Figures A1-A3), which was explained in Section 2.2.
Newaerclim adopted aerosol loading data from the monthly climatology of CAMSiRA (Corpernicus Atmospheric Monitoring Services Interim ReAnalysis) over the period 2003-2014 [41], which is aerosol column concentrations at 3 • by 3 • with a 60 model level [20]. It was provided as an initial condition and interpolated in time and space before the radiation modules of the KIM. The KIM is composed of the Atmosphere 2020, 11, 850 4 of 29 nonhydrostatic dynamic core using spectral elemental methods and the physics package in the cubed sphere [42]. The vertical profiles of aerosol loading data were made by interpolation of the fraction of the total aerosol column concentrations to the KIM model grid and then its normalization to the original total aerosol column mass. The products of those aerosol species column mass concentrations and the single-particle spectral optical characteristics turned to optical properties of individual aerosol types in the radiation modules. The summations of aerosol species optical properties were combined with the optical properties of the other atmospheric components such as clouds and gaseous species, which were fed into the calculation of spectral radiative fluxes. Newaerclim.mod indicates the data produced by modifying the aerosol loading of Newaerclim (more discussion will be in Section 2.3). The three aerosol climatology data used in this study are presented in Table 1.

Aerosol Optical Properties
The aerosol optical characteristics were calculated based on the Mie theory. Mie calculation solves the interactions between electromagnetic waves at a wavelength of λ and homogeneous spherical particles with a radius, r, and refractive index, m [43,44]. Mie theory gives the solution of the vector wave equation, including scattering amplitude functions (S) and scattering phase matrix (P) that relates the Stokes parameter of scattered radiation. S and P provide extinction and scattering factors, and asymmetry parameters. Extinction and scattering factors are obtained: where k = 2π λ is the wavenumber and x = 2πr λ is the size parameter. s g = πr 2 is particle cross-section and Re{z} means the real part of a complex number, z. θ is the scattering angle between incident radiation and object, and φ is the zenith angle. S(θ, φ) is scattering amplitude, which is expressed with Mie angular function (π n and τ n ) and Mie coefficient (a n and b n ). For more detail, information can be referred to in Bohren and Hoffman [45] and Boucher [3]. a n and b n are independent of angles but depend on x and m. The absorption factor, Q abs (λ) is obtained from the differences between extinction and scattering factors.
Single scattering albedo, 0 , that is the ratio of scattering to extinction is expressed as: Atmosphere 2020, 11, 850 5 of 29 Phase function, P(θ) is equal to the element at (1,1) of the scattering phase matrix, P 11 (θ). If we ignore the polarization of the radiation, the normalized phase function can be written as: Asymmetry parameter, g (λ), that is the first moment of phase function, is obtained by: Re a n a * n+1 + b n b * n+1 + (2n + 1) n(n + 1) Re(a n b * n ) (6) where z* means the conjugate of a complex number, z. Therefore, optical characteristics of single particles such as Q ext (λ), Q sca (λ), Q abs (λ), 0 (λ), P(θ, λ), and g(λ) are uniquely determined by given sizes (r) and refractive index (m) of aerosols. In the atmosphere, particle size distribution n(r), are assumed to be lognormal, where N is the total aerosol number. r 0 is the mean geometric radius, σ g = exp(σ o ) is the geometric standard deviation and σ o is the standard deviation. The mass extinction coefficients at λ, σ ext m (λ) is as follows.
where ρ p is particle density. rmax and rmin are the maximum radius and minimum radius, respectively. In the atmosphere, extinction aerosol optical depth and absorption aerosol optical depth at λ of aerosol type, iaer are calculated as follows. σ abs m, iaer (λ)ds (10) where σ ext m, iaer (λ) and σ abs m, iaer (λ) are mass extinction and absorption coefficients (m 2 /g) at λ of aerosol type, iaer, respectively. ds is the path length of aerosols, which is aerosol column mass in the unit of (kg/m 2 ). Then, AOD (absorbing aerosol optical depth (AAOD)) of total aerosols at λ, τ aer_ext (λ) (τ aer_abs (λ)) is obtained by the summation of AOD (AAOD) of individual aerosols at λ: τ aer_ext (λ) = iaer τ ext,iaer (λ) (11) τ aer_abs (λ) = iaer τ ext,iaer (λ) (12) Single scattering albedo of total aerosols at a wavelength of λ, ω aer (λ) is ω aer (λ) = iaer 0 iaer (λ)τ ext,iaer (λ) iaer τ ext,iaer (λ) (13) where, 0 iaer is a single scattering albedo of iaer. The asymmetry parameter of total aerosols at a wavelength of λ, g aer (λ) is g aer (λ) = iaer g iaer (λ) 0 iaer (λ)τ ext,iaer (λ) iaer 0 iaer (λ)τ ext,iaer (λ) (14) where, g iaer is the asymmetry parameter of iaer. τ ext,iaer (λ) and τ abs,iaer (λ) are linearly additive to τ aer_ext (λ) and τ aer_abs (λ). On the other hand, ω aer (λ) and g aer (λ) are weighted by sums of scattering intensity of individual aerosols. Sizes and compositions of five aerosol types were referred to as CAMS-based aerosol climatology [20], as shown in Table 2. Sea salt and dust have three size bins in nucleation, accumulation, and coarse modes. Sea salt has a bimodal distribution, while dust, organic matter (OM), black carbon (BC), and sulfate have monomodal distribution. OM is assumed as the external mixtures of water-soluble, water-insoluble, and soot aerosols, with the mixing ratio of 84:13:3 [20]. External mixing keeps the microphysical and chemical characteristics of the individual aerosol species in a mixture. However, internal mixing makes a mixture with new microphysical and chemical characteristics that differ from the original ones of individual aerosols. External mixtures were reported as less absorbing than internal mixtures to be assumed, such as core-shell or homogeneous mixing types [46][47][48][49]. Although the modeled internal mixtures were not exactly the same as the observed in the atmosphere, the absorption intensity of core-shell types was reported as that of between external mixtures and homogeneous mixing types [49]. The refractive indices of aerosol chemical composition, the sizes of aerosols, and mixing states are essential factors determining the light absorption of aerosols [46]. The aerosol absorption is also affected by the hygroscopicity that changes the sizes and compositions of dry aerosols. We assumed that all sea salt and sulfate aerosols are in 80% relative humidity (RH). Hydrophobic and hydrophilic parts of carbonaceous aerosols were supposed in 20% RH and 50% RH, respectively.  Figure A1 presents the Mie calculations of Newaerclim in blue lines. All aerosol optical properties were interpolated/extrapolated into 14 SW bands in RRTMG. Aerosol optical characteristics of Newaerclim were compared with those of Kimaerclim (black line) [4,47]. There were distinct differences between them. First, the classification standard of aerosol types is the chemical composition in Newaerclim but geographic regions in Kimaerclim. The inconsistency in aerosol types between aerosol loading and aerosol optical characteristics causes unreasonable optical properties in Kimaerclim. SSA of urban aerosols at 532 nm matching with BC loading is too high 0.78 than the well-known BC data (0.3-0.35) [34] ( Figure A2). The asymmetry parameter of urban aerosols has a smaller spectral dependency in the range of 0.6-0.8 than the asymmetry parameter of BC in Newaerclim in the range of 0.0-0.7 ( Figure A3). Second, extinction characteristics in Newaerclim indicate a mass extinction coefficient (m 2 /g), but in Kimaerclim represents the ratio relative to the inverse of AOD at 550 nm (unitless) ( Figure A1). Third, Newaerclim uses more realistic size spectra such as three size bins (fine, accumulate, coarse) for sea salts and dust aerosols, but Kimaerclim uses one size bin for maritime and desert aerosols. Fourth, some aerosol optical properties in Kimaerclim are out of date. SSA of dust-like at 532 nm is as low as 0.75, much lower than Newaerclim dust SSA (0.96, 0.90, and 0.84). It implies that dust-like of the Kimaerclim is more absorbing than the dust of the Newaerclim.
MODIS is onboard Aqua that is at 1:30 p.m. ascending node, sun-synchronous, and near-polar orbit. The MODIS_AOD indicates global coverage monthly data derived from MODIS level 3 daily mean data of Combined Dark Target and Deep Blue AOD at 550 nm for land and ocean. MISR is onboard Terra. The MISR AOD indicates level 3 Global Aerosol product covering a one-month period such as the optical depth average.
For the comparison of modeled AOD and AAOD with observations, interpolation/extrapolation of modeled AOD, MODIS_AOD, and MISR_AOD were made concerning 18 AERONET sites. At most AERONET sites, modeled AODs were in good agreement with monthly variations of three observed AODs. Modeled AAOD agreed well with AERONET_AAOD. Figure 1 shows the scatter plots of AOD and AAOD, and the model results mostly agree with AERONET observations within a factor of two. The linear regression analysis turned out that modeled AOD and AAOD was in good and moderate agreement with AERONET observations as correlation coefficients of 0.807 and 0.631, respectively. Monthly variations of modeled AOD and AAOD of individual aerosols were shown in Figure 2. Dust and sea salt were out of phase, and BC aerosols were the smallest fractions of total AOD. Total AOD peak appeared in May when dust emission was dominant, and the second peak appeared in September when biomass burning was dominant. As for the monthly variations of AAOD, the peak appeared in September due to BC absorption.
On the other hand, Newaerclim revealed a few shortcomings. In the site-comparisons with AERONET data, the modeled AODs overestimated in the summer months at such sites as Banizoumbou, Beijing, Carpentras, Cartel, Ispra, Mexico_City, Sede Boker, and Venise, while they underestimated in winter months. Also, the verification of SW downward fluxes with Clouds and Earth's Radiant Energy System (CERES) turned out to be a dominant underestimation in the desert area (shown in Section 5). It was partly related to the overloading of dust aerosols in CAMSiRA climatology [20]. To obtain the better values of AOD and improve the biases in radiative fluxes, we tried to modulate the dust loadings and produce Newaerclim.mod. We did a curve-fitting to the monthly dust loadings with maximum reduction by 29% in May, and maximum increment by 22% in December. As a result, the monthly variation of modeled AOD turned to double peaks in May and September. The correlation coefficients of AOD and AAOD between the modeled and AERONET improved as 0.863 and 0.693, respectively. When it came to underestimating SW downward radiative fluxes against the Clouds and Earth's Radiant Energy System (CERES), Newaerclim.mod was improved than Newaerclim in July (shown in Section 5).
AODs. Modeled AAOD agreed well with AERONET_AAOD. Figure 1 shows the scatter plots of AOD and AAOD, and the model results mostly agree with AERONET observations within a factor of two. The linear regression analysis turned out that modeled AOD and AAOD was in good and moderate agreement with AERONET observations as correlation coefficients of 0.807 and 0.631, respectively. Monthly variations of modeled AOD and AAOD of individual aerosols were shown in Figure 2. Dust and sea salt were out of phase, and BC aerosols were the smallest fractions of total AOD. Total AOD peak appeared in May when dust emission was dominant, and the second peak appeared in September when biomass burning was dominant. As for the monthly variations of AAOD, the peak appeared in September due to BC absorption.  The left axis is AOD and AAOD of individual aerosol types (itype), and the right axis is AOD and AAOD of total aerosols.
On the other hand, Newaerclim revealed a few shortcomings. In the site-comparisons with AERONET data, the modeled AODs overestimated in the summer months at such sites as Banizoumbou, Beijing, Carpentras, Cartel, Ispra, Mexico_City, Sede Boker, and Venise, while they underestimated in winter months. Also, the verification of SW downward fluxes with Clouds and Earth's Radiant Energy System (CERES) turned out to be a dominant underestimation in the desert area (shown in Section 5). It was partly related to the overloading of dust aerosols in CAMSiRA climatology [20]. To obtain the better values of AOD and improve the biases in radiative fluxes, we tried to modulate the dust loadings and produce Newaerclim.mod. We did a curve-fitting to the monthly dust loadings with maximum reduction by 29% in May, and maximum increment by 22% in December. As a result, the monthly variation of modeled AOD turned to double peaks in May and September. The correlation coefficients of AOD and AAOD between the modeled and AERONET improved as 0.863 and 0.693, respectively. When it came to underestimating SW downward radiative fluxes against the Clouds and Earth's Radiant Energy System (CERES), Newaerclim.mod was improved than Newaerclim in July (shown in Section 5).

Experimental Setup
The behaviors of individual aerosols in atmospheric thermodynamic and hydrometeor fields were examined through sensitivity tests. Each aerosol species was added to the aerosol-free The left axis is AOD and AAOD of individual aerosol types (itype), and the right axis is AOD and AAOD of total aerosols.

Experimental Setup
The behaviors of individual aerosols in atmospheric thermodynamic and hydrometeor fields were examined through sensitivity tests. Each aerosol species was added to the aerosol-free conditions in the KIM and its effects on the thermodynamic and hydrometeor fields were analyzed. The sensitivity results of a 10-day forecast simulation were provided to find the role of aerosol species in the medium-range weather forecasts. Then, we verified a medium-range forecast with aerosol climatology in July of 2017. July and January were conventionally used to verify the medium-range forecast in the KIM as representative months of boreal summer and winter. This study chose July when aerosol loading is more abundant, and solar radiation is more intense than January, in the northern hemisphere, for the distinct signals.
The experimental setups are presented in Tables 3 and 4. The experimental names for sensitivity tests are lower letters such as noaer and newaer, while those for the verification of medium-range forecasts are upper letters such as NOAer and NewAer. The NWP model is used with KIM-v3.3 that comprises the dynamic core and physics package. Physics module contains a rapid radiative transfer model for general circulation models (RRTMG)-based radiation [54], Noah-based land surface model [55], Yonsei University (YSU)-based planetary boundary layer [56], KIM-Arakawa-based orographic gravity wave [57], Simplified Arakawa-Schubert (SAS)-based deep-convection [58], SAS-based shallow convection [59], WRF Single Moment 5-class (WSM5)-based microphysics [60,61], and Tiedtke-based cloudiness [62]. Table 3. Experimental setups for sensitivity test to aerosol types (a 10-day forecast simulation).
x indicates a specific aerosol species such as sea salt, dust, organic matter (OM), black carbon (BC), and sulfate.

Exp. Name
Aerosol Climatology

noaer newaer_all
Newaerclim including all types of aerosols newaer_x Newaerclim including aerosol type, x For the sensitivity to aerosol types, we did experiment on the one 10-day forecast starting at 00:00 UTC on 1 July in 2017 with ne180np3 resolution (approximately 12 km grid spacing). The fifth generation of ECMWF atmospheric reanalysis of the global climate (ERA5) [63] were used as Initial and boundary conditions for atmosphere, and the National Center for Environmental Prediction (NCEP) Global Forecast System (GFS) analysis [64] for sea surface temperature (SST) and sea ice. noaer of the aerosol-free simulation was used as a control run. Then, we simulated newaer_all and newaer_x by adding column mass loading of each aerosol type of aerosol climatology, Newaerclim, to noaer, as shown in Table 3.
For medium-range forecasts, 31-member ensembles of 10-day forecasts were simulated from every 00:00 UTC in July 2017 with the same model resolution, and initial and boundary conditions as those in sensitivity experiments. The four medium-range forecast experiments were carried out with different aerosol climatology ( Table 4). The control run was NOAer with no aerosol climatology. KIMaer, NewAer, and NewAer.mod use Kimaerclim, Newaerclim, and Newaerclim.mod, respectively. For the verification of medium-range forecasts, the Clouds and Earth's Radiant Energy System (CERES), radiosonde, Automatic Weather Station (AWS) data in Korean Peninsula, and the Climate Prediction Center (CPC) ( [65]), and the Tropical Rainfall Measuring Mission (TRMM) Multisatellite Precipitation Analysis (TMPA) [66], were used. We assessed the biases and errors regarding meteorological factors and estimated the forecast score using the KIM self-analysis fields.

Radiative Effects of Aerosols
To analyze aerosol-radiation interactions, we produced aerosol loading in optics units and the radiative forcing at Earth-atmospheric systems. The radiative effects were calculated in the experiment of the 10-day forecast in Table 3, starting at 00:00 UTC on 1 July in 2017 and January 1 in 2017. Figure 3 shows the AOD, AAOD, and scattering AOD (SAOD). The global means of aerosol optical depths in July were larger than those in January by a factor of 1.3, 1.7, and 1.3 for AOD, AAOD, and SAOD, respectively. The fractions of AAOD to AOD were 0.06 in July and 0.04 in January. The distribution AOD in July was distinct in the desert area of Sahara, western India, and the Taklamakan. The primary sources of absorbing aerosols were deserts of the northern hemisphere (NH) and biomass burning regions in central Africa. The sources of AOD in January moved further southward than those in July, and the main AAOD regions are the central western coast of Africa, Bangladesh, and the east coast of China.
Atmosphere 2020, 11, x FOR PEER REVIEW 10 of 28 Radiative forcing is a direct measure of the amount that the Earth's energy budget is out of balance, and it is vital to give the idea how much various radiative forcing agents affect the weather and climate system since industrialization in 1750. Among radiative forcing agents, aerosols are difficult to assess the radiative forcing due to the most significant uncertainties that stem from the inhomogeneous distribution of an aerosol amount in space and time, composition and radiative properties, radiative transfer, the dependence of the forcing calculation on clouds or surface albedos, retrieval errors of remote sensing measurements, and non-aerosol-related errors ( [1], reference therein).
Aerosol radiative forcing due to aerosol-radiation interaction (RFari) is the changes in net irradiances at tropopause as stratospheric temperature profiles adjust to an equilibrium state, resulting in earth surface temperature changes [2]. However, it recently included the rapid adjustment processes in which forcing mechanism affects cloud cover or the other components of the earth-atmospheric Radiative forcing is a direct measure of the amount that the Earth's energy budget is out of balance, and it is vital to give the idea how much various radiative forcing agents affect the weather and climate system since industrialization in 1750. Among radiative forcing agents, aerosols are difficult to assess the radiative forcing due to the most significant uncertainties that stem from the inhomogeneous distribution of an aerosol amount in space and time, composition and radiative properties, radiative transfer, the dependence of the forcing calculation on clouds or surface albedos, retrieval errors of remote sensing measurements, and non-aerosol-related errors ( [1], reference therein).
Aerosol radiative forcing due to aerosol-radiation interaction (RFari) is the changes in net irradiances at tropopause as stratospheric temperature profiles adjust to an equilibrium state, resulting in earth surface temperature changes [2]. However, it recently included the rapid adjustment processes in which forcing mechanism affects cloud cover or the other components of the earth-atmospheric system during a short period (a few hours to a few weeks). The sum of RFari and the rapid adjustment processes is called effective radiative forcing from aerosol-radiation interactions (ERFari) [3]. RFari was estimated as −0.35 ± 0.5 Wm −2 [1] from reanalysis estimates [67], and observational studies [68], with uncertainties mentioned above. Rapid adjustment to aerosol-radiation interactions was assessed as a best estimate of −0.1 Wm −2 [69] in the uncertainty of −0.3 to +0.1 Wm −2 with a 5% to 95% range [70]. Therefore, the overall assessment for ERFari is −0.45 ± 0.5 Wm −2 [1].
The aerosol direct radiative forcing (ADRF) in this study is more likely applicable to the one due to ERFari. However, the ADRF in this study using 10-day forecast data may not be directly compared with those in climate models assessed as global annual mean values. We calculated ADRF as follows. Direct radiative forcing from atmospheric aerosols, ∆F, is defined as the difference in the energy levels between a situation where aerosols are present, F A , and a situation where these atmospheric particles are absent, Fc (Equation (15)). The ADRF was analyzed at the top of atmosphere (TOA), the atmosphere (ATM), and the surface (SFC).
where isky indicates the sky conditions such as "clr" for clear-sky and "all" for all-sky. ilevel indicates TOA or SFC. The arrows indicate the direction of the global fluxes: ↓ is downward flux and ↑ is upward flux. The direct radiative forcing at the atmosphere is the difference between TOA and SFC. The ADRFs at ATM in clear-sky and all-sky are calculated as ∆F clr,ATM and ∆F all,ATM , respectively: The sign of ∆F implies that negative values are associated with aerosol cooling, and positive values are associated with aerosol warming at each level.
The magnitudes of ADRF in the clear sky were larger than those in all-sky because the radiative effects of clouds were added to the modification of radiative transfer. In this study, the ADRFs in July were more significant than those in January, by a factor of 1.3, 2.1, and 1.6 at TOA, ATM, and SFC, respectively. The differences in the ADRFs between July and January were the smallest at TOA, but the differences were the largest at ATM. The signs of ADRF imply the gain of SW radiation at the atmosphere but the loss of SW radiation at the surface and the top of the atmosphere. Such differential vertical structures of ADRF were linked with the horizontal distributions of aerosol optical properties, particularly with aerosol sources. The most significant contributing species to ADRF in July at the surface and atmosphere were dust, and those at TOA were sulfate (Figure 4a). When it came to aerosol optical properties, the largest AODs and AAODs were attributed to dust aerosols in desert areas. In contrast, the most dominant contributors to ADRF in January at surface and TOA were sea salt, while those in the atmosphere were black carbon (Figure 4b). The largest AODs and AAODs were due to sea salts from the ocean and BC aerosols from anthropogenic sources, respectively.
surface and atmosphere were dust, and those at TOA were sulfate (Figure 4a). When it came to aerosol optical properties, the largest AODs and AAODs were attributed to dust aerosols in desert areas. In contrast, the most dominant contributors to ADRF in January at surface and TOA were sea salt, while those in the atmosphere were black carbon (Figure 4b). The largest AODs and AAODs were due to sea salts from the ocean and BC aerosols from anthropogenic sources, respectively.

Sensitivity in Thermodynamic and Hydrometeor Fields
The individual aerosols exert different radiative forcing in terms of its sign, magnitudes, and spatial distribution. The radiative forcing is not necessary to appear to the aerosol source regions [1,9,71], and weather effects of aerosols are not directly related to the aerosol source region, either. The individual aerosols radiative forcing is smaller than total aerosols forcing, and the impacts of individual aerosols on the weather forecast performance may not be distinguishable. For this reason, we tried to examine the responses of individual aerosols to the thermal and hydrometeor fields and to get useful information about the roles of specific aerosol types in the NWP models.
For the effects of aerosols on thermodynamic fields, the changes in sensible heat fluxes (SHflx), latent heat fluxes (LHflx), and PBL height (hpbl) relative to noaer were investigated ( Figure 5). Due to the reduced radiative fluxes, the surface temperature decreased, and stability in the lower atmosphere increased. Thus, the turbulent mixing of SHflx and LHflx decreased, resulting in boundary layer structures of wet surfaces and dry above. Overall, Earth's surface became cold and wet due to aerosols; however, it got too warm and wet due to BC aerosols. Among five types of aerosols, dust aerosols made the most contributions to the changes in heat fluxes and atmospheric stability. Sea salt aerosols made fewer contributions to the thermodynamic field than any other aerosol types in July.
Atmosphere 2020, 11, x FOR PEER REVIEW 12 of 28 The individual aerosols exert different radiative forcing in terms of its sign, magnitudes, and spatial distribution. The radiative forcing is not necessary to appear to the aerosol source regions [1,9,71], and weather effects of aerosols are not directly related to the aerosol source region, either. The individual aerosols radiative forcing is smaller than total aerosols forcing, and the impacts of individual aerosols on the weather forecast performance may not be distinguishable. For this reason, we tried to examine the responses of individual aerosols to the thermal and hydrometeor fields and to get useful information about the roles of specific aerosol types in the NWP models.
For the effects of aerosols on thermodynamic fields, the changes in sensible heat fluxes (SHflx), latent heat fluxes (LHflx), and PBL height (hpbl) relative to noaer were investigated ( Figure 5). Due to the reduced radiative fluxes, the surface temperature decreased, and stability in the lower atmosphere increased. Thus, the turbulent mixing of SHflx and LHflx decreased, resulting in boundary layer structures of wet surfaces and dry above. Overall, Earth's surface became cold and wet due to aerosols; however, it got too warm and wet due to BC aerosols. Among five types of aerosols, dust aerosols made the most contributions to the changes in heat fluxes and atmospheric stability. Sea salt aerosols made fewer contributions to the thermodynamic field than any other aerosol types in July. Unlike negative surface radiative forcing, positive atmospheric radiative forcing affected vertical profiles of atmospheric thermodynamic and hydrometeor fields ( Figure 6). Due to the less moisture mixing at the lower atmosphere and the evaporation in a warmer atmosphere, the differences in specific humidity (Q) had transition points about 800 hPa. Below 500 hPa, SW heating rates (SWhr) Unlike negative surface radiative forcing, positive atmospheric radiative forcing affected vertical profiles of atmospheric thermodynamic and hydrometeor fields ( Figure 6). Due to the less moisture mixing at the lower atmosphere and the evaporation in a warmer atmosphere, the differences in specific humidity (Q) had transition points about 800 hPa. Below 500 hPa, SW heating rates (SWhr) increased due to absorption by aerosols, and the longwave (LW) cooling rates increased due to increments of LW fluxes emitted from Earth's surface and clouds out to space. The cold and wet surface but the warm and dry atmosphere made the cloud water content (Qc) increase below 700 hPa and decrease above. In particular, the maximum of Qc appeared at 900 hPa because drying and warming atmosphere caused more evaporation of water vapor and stabilized the boundary layer to lower its top, while the abundance of water vapor caused to lower lifting condensation levels (LCL), so-called cloud-top thinning and cloud-bottom thickening (Figure 6b). Such a pattern of Qc was similar types of cloud cover profiles that Persad et al. [72] found when BC aerosols were assumed to exist in the boundary layer [72], Figure 4c of tropical regions. It was caused by the light-absorption by BC aerosols without consideration of the indirect effects of BC. The aerosol also reduced the cloud cover above the mid-level atmosphere and decreased cloud ice content (Qi) at the entire atmosphere. Thus, precipitation such as rain (Qr) and snow (Qs) mixing ratios decreased when compared to the aerosol-free conditions.
The changes in thermodynamic and hydrometeor fields of individual aerosols cases (newaer_x) were smaller than those of total aerosols cases (newaer_all) (Figures 7 and 8) and they depended on their concentrations and optical properties. Figure 7 shows that the differences in heating rates pretty linearly responded to the absorption capacities of individual aerosol types. Dust made the most contribution to the changes in SW heating rates and LW heating rates in July among aerosol species. In contrast, the differences in temperature of newaer_all increased with the double-positive peaks at the lower and upper atmospheres. Such temperature differences appeared in newaer_dust and newaer_BC. However, temperature differences became negative due to sulfate in the entire atmosphere and dust below 900 hPa. It turned out that the decreases in temperature at the lower atmosphere were due to scattering by dust and sulfate, whereas the increases in temperature above the mid-level atmosphere were due to absorption by dust and black carbon.
On the other hand, the differences in geopotential height (GPH) relative to aerosol-free conditions were determined by thermal information accumulated from the bottom of the atmosphere (Figure 7). All the individual aerosol cases (newaer_x) increased GPH except sulfate (newaer_sulfate). GPH differences were negative at the surface due to the cooling of dust and sulfate aerosols, and they were positive above the mid-atmosphere due to the heating of absorbing aerosols. The most significant contributor to GPH changes in July was BC, which had the most potent Thus, precipitation such as rain (Qr) and snow (Qs) mixing ratios decreased when compared to the aerosol-free conditions.
The changes in thermodynamic and hydrometeor fields of individual aerosols cases (newaer_x) were smaller than those of total aerosols cases (newaer_all) (Figures 7 and 8) and they depended on their concentrations and optical properties. Figure 7 shows that the differences in heating rates pretty linearly responded to the absorption capacities of individual aerosol types. Dust made the most contribution to the changes in SW heating rates and LW heating rates in July among aerosol species. In contrast, the differences in temperature of newaer_all increased with the double-positive peaks at the lower and upper atmospheres. Such temperature differences appeared in newaer_dust and newaer_BC. However, temperature differences became negative due to sulfate in the entire atmosphere and dust below 900 hPa. It turned out that the decreases in temperature at the lower atmosphere were due to scattering by dust and sulfate, whereas the increases in temperature above the mid-level atmosphere were due to absorption by dust and black carbon.   On the other hand, the differences in geopotential height (GPH) relative to aerosol-free conditions were determined by thermal information accumulated from the bottom of the atmosphere (Figure 7). All the individual aerosol cases (newaer_x) increased GPH except sulfate (newaer_sulfate). GPH differences were negative at the surface due to the cooling of dust and sulfate aerosols, and they were positive above the mid-atmosphere due to the heating of absorbing aerosols. The most significant contributor to GPH changes in July was BC, which had the most potent absorption characteristics among aerosol composition. As for the difference in specific humidity relative to noaer, most aerosol types had positive changes in the lower atmosphere and negative changes in the mid-atmosphere. Among five aerosol types, the dust seemed to be the most influential on the changes due to total aerosols (newaer_all), which were moistened at a lower atmosphere and dry at mid-level atmosphere. In contrast, sulfate made the atmosphere cool and dry as the same as its surface impacts (shown in Figure 5b) even though it made a smaller contribution to the changes in the specific humidity than any other aerosol species.
Atmosphere 2020, 11, x FOR PEER REVIEW 15 of 28 individual aerosols (newaer_x) in the lower atmosphere ( Figure 8b). However, Qi profiles of newaer_all at 550 hPa and 300 hPa were affected by dust and BC aerosols, respectively. The differences in rainwater content (Qr) in newaer_x were negative below 500 hPa except those in newaer_seasalt. Qr decreased in the lower atmosphere due to BC and OM and in the midatmosphere due to dust. As for the changes in snow content (Qs) profiles, all types of aerosols in newaer_x decreased Qs most substantially around 700 hPa as compared to noaer. The dust case showed considerable negative changes at 550 hPa, which was similar to the case of newaer_all. However, BC, OM, and sea salt increased Qs above 500 hPa due to the pleasant condition to form ice nuclei. Such features of hydrometeor fields due to aerosol species appeared in the global precipitation changes (Figure 9). Overall, the global changes in convective precipitation (PRECC) were induced by absorbing aerosols such as dust, BC, and OM in decreasing order. The negative changes in PRECC of newaer_x were linearly additive to those of newaer_all. In contrast, the global changes in large-scale precipitation (PRECL) were affected by most aerosol types.  Figure 8 shows how specific aerosols change the hydrometeor relative to aerosol-free conditions(noaer). Most aerosol species increased cloud liquid water contents (Qc) at the lower atmosphere and decreased Qc at the mid-atmosphere (Figure 8b). The Qc changes of individual aerosols were not linearly additive to those of the newaer_all. Instead, the vertical profiles of Qc differences in newaer_all were similar to those in absorbing aerosols such as dust, BC, and OM. The changes in cloud ice content (Qi) profiles of newaer_all were close to the mean values of those in individual aerosols (newaer_x) in the lower atmosphere (Figure 8b). However, Qi profiles of newaer_all at 550 hPa and 300 hPa were affected by dust and BC aerosols, respectively.
The differences in rainwater content (Qr) in newaer_x were negative below 500 hPa except those in newaer_seasalt. Qr decreased in the lower atmosphere due to BC and OM and in the mid-atmosphere due to dust. As for the changes in snow content (Qs) profiles, all types of aerosols in newaer_x decreased Qs most substantially around 700 hPa as compared to noaer. The dust case showed considerable negative changes at 550 hPa, which was similar to the case of newaer_all. However, BC, OM, and sea salt increased Qs above 500 hPa due to the pleasant condition to form ice nuclei. Such features of hydrometeor fields due to aerosol species appeared in the global precipitation changes (Figure 9). Overall, the global changes in convective precipitation (PRECC) were induced by absorbing aerosols such as dust, BC, and OM in decreasing order. The negative changes in PRECC of newaer_x were linearly additive to those of newaer_all. In contrast, the global changes in large-scale precipitation (PRECL) were affected by most aerosol types.
Atmosphere 2020, 11, x FOR PEER REVIEW 16 of 28 Figure 9. The sensitivity of precipitation to aerosol types. The changes in convective precipitation (PRECC) and large-scale precipitation (PRECL) relative to aerosol-free conditions for a 10-day forecast starting from 1 July of 2017.
The differential radiative fluxes with altitudes caused changes in thermodynamic fields. Through ARIs, aerosols tended to change the vertical gradients of temperature and humidity and to stabilize the atmosphere, leading to making the atmosphere with fewer clouds and less precipitation with respect to noaer. The differences due to individual aerosols were linearly additive to those due to total aerosols in heat fluxes, heating rates, humidity, and PRECC; however, they were not linear in temperature, geopotential height, and cloud liquid or ice contents, and PRECL. Dust was the most influential on the direct and semidirect effects of aerosols in July, among absorbing and nonabsorbing aerosols.

Weather Effects of Aerosols
To find the impacts of aerosols on weather forecasts, ensemble simulations of 10-day forecasts from 1 July to 31 July 2017, were verified. Figure 10 shows the differences in SW downward radiation in the 54-72 h forecast against CERES in NOAer, KIMAer, NewAer, and NewAer.mod. In clear-sky, NOAer showed the positive biases in the entire Earth's surface, particularly, desert areas. Adding aerosols to NOAer made such positive biases disappear, rather negative biases turned out in the west Indian ocean and downwind of the Saharan desert. The global mean biases in the experiment with aerosol climatology decreased by up to 99% as compared with those in NOAer, and the global mean root mean square errors (RMSEs) decreased by more than 50% (Table 5). In the KIMAer negative biases over the ocean and the Arctic regions appeared due to the scattering by the overloading of sea salt aerosols. On the other hand, most biases were reduced particularly over the ocean in the NewAer, while negative biases appeared over northern Africa, western India, and Taklamakan desert. In NewAer.mod, negative biases over the desert reduced, but the rest of the negative biases remained due to absorption by BC aerosols. However, the global mean of the biases was smallest in NewAer.mod.
In all-sky radiation, the biases against CERES included model biases caused by clouds and precipitation. NOAer had the positive biases in continents of NH and the west coast of South America, and the negative biases in the ocean at ITCZ and the Arctic sea. When the KIM in all-sky was compared with the KIM in clear-sky, the KIM seemed to underestimate clouds over the continents, but it overestimated over the tropical oceans. The global bias distribution patterns in all-sky seemed to be similar among the four cases. The ERFari in the KIM improved the biases in continents, but it did not drastically affect biases in the ocean and enhanced negative biases in the ITCZ. The global mean RMSEs in the experiment with aerosol climatology decreased only by 10% as compared with Figure 9. The sensitivity of precipitation to aerosol types. The changes in convective precipitation (PRECC) and large-scale precipitation (PRECL) relative to aerosol-free conditions for a 10-day forecast starting from 1 July of 2017.
The differential radiative fluxes with altitudes caused changes in thermodynamic fields. Through ARIs, aerosols tended to change the vertical gradients of temperature and humidity and to stabilize the atmosphere, leading to making the atmosphere with fewer clouds and less precipitation with respect to noaer. The differences due to individual aerosols were linearly additive to those due to total aerosols in heat fluxes, heating rates, humidity, and PRECC; however, they were not linear in temperature, geopotential height, and cloud liquid or ice contents, and PRECL. Dust was the most influential on the direct and semidirect effects of aerosols in July, among absorbing and nonabsorbing aerosols.

Weather Effects of Aerosols
To find the impacts of aerosols on weather forecasts, ensemble simulations of 10-day forecasts from 1 July to 31 July 2017, were verified. Figure 10 shows the differences in SW downward radiation in the 54-72 h forecast against CERES in NOAer, KIMAer, NewAer, and NewAer.mod. In clear-sky, NOAer showed the positive biases in the entire Earth's surface, particularly, desert areas. Adding aerosols to NOAer made such positive biases disappear, rather negative biases turned out in the west Indian ocean and downwind of the Saharan desert. The global mean biases in the experiment with aerosol climatology decreased by up to 99% as compared with those in NOAer, and the global mean root mean square errors (RMSEs) decreased by more than 50% (Table 5). In the KIMAer negative biases over the ocean and the Arctic regions appeared due to the scattering by the overloading of sea salt aerosols. On the other hand, most biases were reduced particularly over the ocean in the NewAer, while negative biases appeared over northern Africa, western India, and Taklamakan desert. In NewAer.mod, negative biases over the desert reduced, but the rest of the negative biases remained due to absorption by BC aerosols. However, the global mean of the biases was smallest in NewAer.mod. Atmosphere 2020, 11, x FOR PEER REVIEW 17 of 28    In all-sky radiation, the biases against CERES included model biases caused by clouds and precipitation. NOAer had the positive biases in continents of NH and the west coast of South America, and the negative biases in the ocean at ITCZ and the Arctic sea. When the KIM in all-sky was compared with the KIM in clear-sky, the KIM seemed to underestimate clouds over the continents, but it overestimated over the tropical oceans. The global bias distribution patterns in all-sky seemed to be similar among the four cases. The ERFari in the KIM improved the biases in continents, but it did not drastically affect biases in the ocean and enhanced negative biases in the ITCZ. The global mean RMSEs in the experiment with aerosol climatology decreased only by 10% as compared with those in NOAer while the global mean biases decreased by more than 75% (Table 4). The biases in the continents and the north Pacific improved most in NewAer.mod, and global mean biases were the smallest in KIMAer due to smaller positive biases in the high latitude than in the other two aerosol cases.
Under such radiative fluxes modified by aerosols, the 10-day forecasts were verified with radiosonde data of GPH, temperature, specific humidity, and wind speed ( Figure 11). In this study, the changes in specific humidity and wind speed due to aerosols were negligible. In NOAer, GPH had the negative biases at the upper and lower atmospheres after the fifth-day forecasts. Positive biases appeared in the mid-atmosphere in Asia in the early days of forecasts and the upper atmosphere in the tropics after fifth-day forecasts. The biases of temperature in NOAer had similar vertical stratification to those of GPH. However, the positive biases in temperature were located at lower altitudes than those in GPH were because temperature responses appeared closer to the aerosol source locations than GPH responses. On the reflection of the impacts of aerosol types on GPH and temperature in the sensitivity tests, BC and OM heated at the lower atmosphere, and dust and BC heated at the upper atmosphere ( Figure 7c). BC aerosols also had dominant effects on positive changes of GPH (Figure 7d). It resulted in decreasing negative biases of temperature and GPH at the upper and lower atmospheres relative to NOAer. However, the positive biases at mid-atmosphere did not improve because the roles of absorbing aerosols seemed predominant in GPH and temperature. The opposite (positive and negative) biases seemed not significantly to change the global and 10-day means of the biases because they tended to cancel each other. Among the three aerosol experiments, the biases in GPH and temperature improved the most in NewAer.mod. The positive biases in the mid-atmosphere, which tended to be enhanced due to the existence of absorbing aerosols, were the smallest in NewAer.mod. The negative biases at the lower and upper atmospheres in NewAer.mod were larger than KIMAer; however, they were improved as compared with NOAER.
global and 10-day means of the biases because they tended to cancel each other. Among the three aerosol experiments, the biases in GPH and temperature improved the most in NewAer.mod. The positive biases in the mid-atmosphere, which tended to be enhanced due to the existence of absorbing aerosols, were the smallest in NewAer.mod. The negative biases at the lower and upper atmospheres in NewAer.mod were larger than KIMAer; however, they were improved as compared with NOAER.

NOAer
KIMAer NewAer NewAer.mod The statistical skill scores for the precipitation were calculated against observations such as the CPC unified gauge-based precipitation analysis, the TMPA, and the AWS data. The CPC and TMPA were used for the global precipitation in continents and ocean, respectively, and the AWS was used for the regional precipitation over the Korean Peninsula. However, the verification results were very similar, and the AWS was the most distinct among the three observations. NOAer showed the overestimation in weak (~5 mm/d) or moderate (10-15 mm/d) precipitation but the underestimation in heavy (~20 mm/d) precipitation. It turned out that aerosols unilaterally reduced precipitation, and thus, they weakened positive biases in weak and moderate precipitation but enhanced negative biases in heavy precipitation. On the other hand, the Equitable Threat Score (ETS) that validates the accuracy of forecasts without contribution from the hit by chance in random forecast [73], improved in the three aerosol experiments as compared with the NOAer. Figure 12 showed the bias and ETS The statistical skill scores for the precipitation were calculated against observations such as the CPC unified gauge-based precipitation analysis, the TMPA, and the AWS data. The CPC and TMPA were used for the global precipitation in continents and ocean, respectively, and the AWS was used for the regional precipitation over the Korean Peninsula. However, the verification results were very similar, and the AWS was the most distinct among the three observations. NOAer showed the overestimation in weak (~5 mm/d) or moderate (10-15 mm/d) precipitation but the underestimation in heavy (~20 mm/d) precipitation. It turned out that aerosols unilaterally reduced precipitation, and thus, they weakened positive biases in weak and moderate precipitation but enhanced negative biases in heavy precipitation. On the other hand, the Equitable Threat Score (ETS) that validates the accuracy of forecasts without contribution from the hit by chance in random forecast [73], improved in the three aerosol experiments as compared with the NOAer. Figure 12 showed the bias and ETS in the verification of all the experiments with AWS in the Korean peninsula. Among the three experiments with aerosols, NewAer showed the best performance in weak and moderate precipitation, and KIMAer made better performance in heavy precipitation than the other experiments. In the changes in surface precipitation by aerosols, PRECL was reduced by all aerosol types, whereas PRECC was reduced by only absorbing aerosols (Figure 9). Since the changes in concentration of absorbing aerosols tended to affect both PRECC and PRECL, absorbing aerosols were expected to enhance biases in the heavy precipitation. In the global distribution of precipitation biases to CPC or TMPA, NOAer overestimated precipitation in the continent and at high latitudes and underestimated in the tropical regions. Aerosols weakened the overestimation at high latitude, central Africa, and northmost of South America due to BC and OM. However, the underestimation in Southeast Asia was enhanced by carbonaceous aerosols (not shown here). In order to assess the statistical scores in weather forecasts, we calculated anomaly correlation (AC), root mean square error (RMSE), and skill score (SS) against the KIM self-analysis fields in GPH, temperature, relative humidity, and wind speed. The highlight of statistical analyses in this study were shown in Figure 13. The forecast scores with aerosols were not very different from those without aerosols. Nevertheless, GPH improved in the RMSEs at the upper atmosphere of NH, the SS at the high latitude, and the AC at the 500 hPa. The AC of mean sea level pressure also improved. On the other hand, forecast scores in SH became worse, such as the ACs of temperature at lower and upper atmospheres and the AC of GPH at the upper atmosphere. As shown in the sensitivity test of thermodynamic fields to aerosol types (Section 4), the changes in GPH profiles resulted mainly from BC and dust aerosols. Double-peak temperature profiles by BC and dust aerosols (Figure 7c), caused overheating at the lower and upper atmospheres. It led to decreasing in the AC of temperature at tropics and SH. The accumulated information of such temperature anomalies made the AC of GPH decrease at the upper atmosphere of SH. On the other hand, sea salts from vast oceans in SH seemed to play an additional role in the current NWPs besides ARI. In order to assess the statistical scores in weather forecasts, we calculated anomaly correlation (AC), root mean square error (RMSE), and skill score (SS) against the KIM self-analysis fields in GPH, temperature, relative humidity, and wind speed. The highlight of statistical analyses in this study were shown in Figure 13. The forecast scores with aerosols were not very different from those without aerosols. Nevertheless, GPH improved in the RMSEs at the upper atmosphere of NH, the SS at the high latitude, and the AC at the 500 hPa. The AC of mean sea level pressure also improved. On the other hand, forecast scores in SH became worse, such as the ACs of temperature at lower and upper atmospheres and the AC of GPH at the upper atmosphere. As shown in the sensitivity test of thermodynamic fields to aerosol types (Section 4), the changes in GPH profiles resulted mainly from BC and dust aerosols. Double-peak temperature profiles by BC and dust aerosols (Figure 7c), caused overheating at the lower and upper atmospheres. It led to decreasing in the AC of temperature at tropics and SH. The accumulated information of such temperature anomalies made the AC of GPH decrease at the upper atmosphere of SH. On the other hand, sea salts from vast oceans in SH seemed to play an additional role in the current NWPs besides ARI.
Atmosphere 2020, 11, 850 21 of 29 BC and dust aerosols. Double-peak temperature profiles by BC and dust aerosols (Figure 7c), caused overheating at the lower and upper atmospheres. It led to decreasing in the AC of temperature at tropics and SH. The accumulated information of such temperature anomalies made the AC of GPH decrease at the upper atmosphere of SH. On the other hand, sea salts from vast oceans in SH seemed to play an additional role in the current NWPs besides ARI. Through ARI, aerosols worked unilaterally to decrease SW downward radiation at the surface, increase GPH and temperature, and decrease precipitation relative to aerosol-free atmosphere. The light-absorption capacity of aerosols played essential roles in weather forecasts. It led to weakening positive biases in SW downward radiation at the surface but occurring negative biases in dust and BC source regions; improving negative biases in GPH and temperature at the upper and lower Through ARI, aerosols worked unilaterally to decrease SW downward radiation at the surface, increase GPH and temperature, and decrease precipitation relative to aerosol-free atmosphere. The light-absorption capacity of aerosols played essential roles in weather forecasts. It led to weakening positive biases in SW downward radiation at the surface but occurring negative biases in dust and BC source regions; improving negative biases in GPH and temperature at the upper and lower atmospheres but worsening positive biases in the mid-atmosphere; reducing positive biases in weak and moderate precipitation but enhancing negative biases in heavy precipitation.
Those biases imply that there is room for the processes associated with the ERFari, such as aging and mixing of aerosol types and aerosol spectral optical properties in longwave radiations. Aerosol loading can be replaced with near-real-time (NRT)-base scenario in the specific or severe aerosol events. In the weather effects of ERFari, aerosols unilaterally changed the biases of the global forecast models, resulting in improving or exacerbating the biases of thermodynamic fields. Therefore, missing processes such as the aerosol-cloud interactions are likely to change the biases distribution patterns in the all-sky radiative fluxes and the tendency of reduction in cloud and precipitation.
This study suggests that Newaerclim.mod is the most desirable to use in the KIM among the three aerosol climatology data even though it does not drastically improve weather forecast scores as compared with the other two. Because it is a well-posed aerosol climatology, that is logical and consistent between aerosol loading and single-particle characteristics, and as realistic and recent as possible. Thus, Newaerclim.mod makes it easier to improve further bias corrections with the science-base and extend to the other aerosol roles in the hydrological cycles, surface albedo, and atmospheric chemistry, with consistency.

Conclusions
The roles of aerosols in NWP models were investigated by using well-posed aerosol data in Korean Integrated Model. The ARDF was calculated for a 10-day forecast on 1 July and 1 January 2017. In July, the ADRFs were larger than those in January, by a factor of 1.3, 2.1, and 1.6 at TOA, ATM, and SFC, respectively. To use a more distinct signal for the responses of thermodynamic and hydrometeor fields to individual aerosols, we carried out a sensitivity test for July. The medium-range weather forecasts with 31 ensemble simulations in July 2017 were verified with observations, and the model biases and the performances in weather forecasts were assessed in the cases of NOAer, KIMAer, NewAer, and NewAer.mod. We found the roles of individual aerosols in weather effects through the linkage between the sensitivity test and weather forecast performances. AOD were simulated well in terms of the peaks of individual aerosol types and relative ratios among aerosol species. The representative agents of July ADRF were dust at the surface and the atmosphere, and sulfate at the top of the atmosphere, respectively. • Negative ADRFs at the surface had an impact on decreasing in sensible and latent heat fluxes at the surface. Positive ADRFs in the atmosphere caused the differences in atmospheric thermodynamic and hydrometeor fields relative to the aerosol-free conditions. The sensitivity tests showed that double-positive maxima in temperature were due to dust and BC, and wet surface and dry above in moisture mainly due to dust. Also, decreasing in cloud liquid or ice contents but a sharp positive layer of cloud liquid contents at 900 hPa were due to dust, BC, and OM. Through ARIs, aerosols tended to change the vertical gradients of temperature and humidity and to stabilize the atmosphere, leading to making the atmosphere with fewer clouds and less precipitation with respect to noaer. The differences due to individual aerosols were linearly additive to those due to total aerosols in heat fluxes, heating rates, humidity, and PRECC; however, they have not additive linearity to the differences due to total aerosols in temperature, geopotential height, and cloud liquid or ice contents, and PRECL. Dust was the most influential among absorbing and nonabsorbing aerosols, on the direct and semidirect effects of aerosols in July.

•
Through ARI, aerosols worked unilaterally to decrease SW downward radiation at the surface, increase GPH and temperature, and decrease precipitation relative to aerosol-free atmosphere. The capacity of absorbing light due to aerosols played essential roles in weather forecasts. The verification of medium-range forecasts revealed that aerosols weakened positive biases in SW downward radiation but occurred negative biases in dust and BC source regions. Aerosols improved negative biases in GPH and temperature at the upper and lower atmospheres but worsened positive biases in the mid-atmosphere due to absorbing aerosols. Moreover, they reduced positive biases in weak and moderate precipitation but enhanced negative biases in heavy precipitation mainly due to absorbing aerosols. According to the sensitivity test of hydrometeor fields and the global distribution of precipitation biases, dust and OM weakened the overestimation at high latitude, central Africa, and northmost of South America. In contrast, BC and sea salt enhanced the underestimation in southeast Asia and the southern Indian Ocean, respectively. • It turned out that weather forecast scores, including current aerosol information, improved in GPH of NH, however, they got worse in temperature of SH and GPH at the upper atmosphere of SH.
The current well-posed aerosol climatology can extend to more realistic descriptions of aerosol optical properties such as aging processes, mixing states, and spectral optical properties in the longwave radiation. Among the three aerosol climatology data, Newaerclim.mod is the most desirable data for the KIM in terms of the well-posed data frame, the bias corrections in the KIM with the science-base, and potential in cooperating with various physical and chemical aspects in NWP systems. It is worthy that prognostic aerosols produce the near-real-time (NRT) aerosol distribution and much more accurate aerosol radiative forcing when severe aerosol events break out. Moreover, they make it possible to execute the radiative forcing from aerosol-cloud interactions (RFaci) that should be included in the NWP models.