Last updated: 3 Sep 2023
March 2020 to February 2022, with phase 2 extended to September 2024
Total FTEs used:
initial project (2020)
During the last year, a new version of RRTM has been introduced into ICON, which is the new ECMWF radiation scheme for IFS, ecRAD (Hogan and Bozzo, 2018). Based on this preparation, a window for further development and testing is now open. In this project plan, we propose the implementation of considerable new developments regarding cloud optics and related aerosols-radiation or aerosols-cloud interaction into ICON, so as to be active in ecRAD. The main part of this development has already been done during PP T2(RC)2 and is now a part the COSMO model. Although related extensions have even been transferred to the current version of RRTM in ICON (henceforth referred to as RRTM), ecRAD has not yet been extended accordingly. Further, new sources of associated aerosol data are not yet accessible for ICON. Moreover, we want to investigate to what extend the statistical information of a stochastic shallow convection scheme, which is already present in a sub-development branch of ICON, can be integrated into the formulation of cloud-radiation interaction. Finally, we also intend to perform massive testing of the new code by comparison of model results with observations.
project extension (phase2, 2022)
The project started in March 2020 with a collaboration of IMS, RHM and DWD colleagues. We proposed the implementation of considerable new developments regarding cloud optics and related aerosols-radiation or aerosols-cloud interaction into ICON, and specifically into the new ecRAD radiation scheme recently introduced in ICON. Three main subjects were covered in the original proposal. The first is the introduction of new cloud droplets and ice particles optical properties (Hu & Stamness 1993, Fu 1996, 1998). The second is the introduction of new aerosols inputs to the radiation scheme: CAMS climatology, MacV2 climatology, CAMS forecasted aerosols and the 2D AOD advection scheme. The third is an investigation of the stochastic shallow convection scheme.
On top of all of these, massive testing of the new code by comparison of model results with default setups and observations was planned. Some of the mentioned tasks were finished or significantly promoted, while few have not yet started or ought to be altered for this second phase.
The major changes from the original project plan are:
Radiation transfer models are one of the main elements in NWP models, which can drastically contribute to heating and cooling rates in all levels of the atmosphere. This can critically change the dynamics in the atmosphere and dramatically affect the model skills. On the other hand, radiation schemes are extremely complicated due to the large amount of gases, hydrometeors and aerosol species being involved in radiation-transfer processes. For each of these air components, the governing optical properties need to be parametrized in each considered spectral interval of the model as a function of proper bulk properties of these species valid for a grid box, which is primarily the effective radius.
During the years 2015-2019, the radiation scheme in the COSMO model was updated as part of a Priority Project (PP) Testing and Tuning of the Revised Cloud Radiation Coupling T2(RC)2. The revision included updating the optical properties of water droplets Hu & Stamnes (1993) and ice particles Fu (1996, 1998).
The new parametrizations expanded the range of validity, namely towards larger dimensions of the included hydrometeors, up to the size of rain and snow particles. These parametrizations have been implemented in the official radiation schemes for both COSMO (the Ritter & Geleyn, 1992) and ICON (RRTM). Other related improvement, up to now only implemented in COSMO, includes the utilization of MACv2 aerosols climatology (Kinne 2013), prognostic ICON-ART dust and the prognostic CAMS aerosols as a direct effect in the radiation scheme.
This aerosol information has also been employed as an input for the process of water and ice nucleation, representing an indirect aerosol-effect on radiation transfer. The project also included a new water droplets effective radius parametrization and a revised formulation of the impact of sub-grid scale (SGS) cloud generation on radiation processes. Finally, computational optimizations, like the Monte-Carlo Spectral Integration (MCSI) scheme, have been implemented into COSMO in order to reduce the run-time of radiation calculations. The respective code and its new features have thoroughly been tested and tuned against ground based and satellite measurements, performed in Israel, Germany and Russia under different weather conditions. Through this, a positive impact of the new development on net radiation fluxes and near-surface temperatures could be demonstrated.
During the second phase of T2(RC)2 the new cloud optics have also been implemented into the currently operational radiation scheme of ICON (RRTM), including the new formulation of the SGS cloud-impact. Like in COSMO, all the mentioned improvement together has proven to cause better radiation fluxes also in ICON. However, the new ecRAD code, which is now available in ICON, has not yet been adapted to these promising developments. Moreover, for liquid clouds, the effective radius calculation in the ICON code needs to be updated, as it has already been done in COSMO.
The first topic of this project concerns with cloud optics. As mentioned earlier, new ice optical properties parametrizations following Fu (1996, 1998) have already been introduced and tested in ICON-RRTM. In ecRAD the ice optics are also based on the same author (Hogan and Bozzo, 2018), but they are implemented with a different spectral averaging technique. Now we wish to analyse the different schemes and to define the best-suited solution for ICON-ecRAD. Regarding the water droplets optical properties, ecRAD uses the SOCRATES package (Hogan and Bozzo, 2018 and references therein) while, in ICON, we lately introduced the Hu & Stamnes (1993) parametrization.
We wish to implement in the ecRAD framework due to several advantages such as a single simplistic formulation for the whole spectrum (long and short wave intervals) and an advanced spectral integration suitable for NWP models. The new proposed optical properties are now parametrized using formulas that uses effective radius Reff or aspect ratio (AR) as a parameter. However the future plan for ecRAD is using look-up table (LUT) approach proposed by Robin Hogan (ECMWF). Therefore a transformation from formulas to LUT is needed.
The cloud droplets bulk optical properties are defined by the effective radius Reff. For the use in ecRAD, Reff is treated either as a constant or, using the ECHAM5 method, as a function of liquid water content and cloud number concentration. We plan on implementing a prognostic approach, which was tested in COSMO, where the Segal & Khain (2006) activation scheme is used to define the cloud number concentration from prognostic CAMS aerosols number concentration and effective updraft speed. This prognostic cloud number concentration is later used to define a realistic Reff using the Khain et al. (2019) scheme.
The second goal of this project is to introduce new aerosols inputs for radiation, suited for the operational ICON mode. Currently the Tegen (1997) climatology is employed in ICON operational setup. These monthly averaged optical depth values of aerosols are used to define a number concentration profile, which, consequently, is not affected by the model dynamics. It is a well-known fact that aerosols are a major element in both radiative transfer and for cloud formation. The direct effect on radiation can be easily seen in dust outbreaks and volcanic ash eruptions, but even on regular conditions aerosols concentration can have large day-to-day variability imposing large errors when climatological averages are used. This is the reason that in COSMO model we added three new aerosols inputs in addition to the old climatologies of Tanre (1983) and Tegen (1997).
First is the CAMS prognostic aerosols provide a 5-day lead-time global forecast for the 3D mixing ratios of five types of aerosols, which are divided into eleven aerosols tracers with an hourly temporal resolution. As the CAMS aerosols optical properties are already a part of ecRAD, their implementation can be done straightforwardly.
The second is the global ICON-ART dust forecast, which is now running semi-operationally twice a day at DWD. These data have been made available for COSMO as well.
The third is the most up to date aerosols climatology available, namely the MACv2 by Kinne (2013). In the same way, these data should be integrated into ICON NWP-runs and can then be attached to ecRAD as well. Lately, a simplified and fast 2D prognostic scheme for mineral dust has been developed by Günther Zängl (DWD), using the Tegen (1997) climatology as initial conditions. In the course of our project, we aim to test this method and further develop it for the application to other aerosols species and different initial conditions. One or more of the mentioned aerosols inputs can also be used as condensation nuclei for the microphysical scheme using Segal & Khain method, which can improve precipitation rates and can improve the effective radius calculation for the radiation-indirect effect. Notice, that the prognostic ART component usually can’t be active in NWP runs with ICON. Hence, ICON-ART data, which so far are restricted to mineral dust only, need to be provided as offline products as well. However, DWD plans to include also sea salt and dust from vegetation burn in the next future (project PERMASTROM).
Another target of this project deals with the shallow convection parametrization and its effect on the radiation. Shallow convection (SC) parametrization plays a crucial role in the evolution of simulated boundary layer structure. This effect and the interaction with the grid scale dynamics have a significant impact on cloudiness and precipitation. For example, over the eastern Mediterranean, COSMO SC parametrization reasonably well mixes the boundary layer in fair weather, but dries out the lower levels and thus prevents grid scale precipitation in unstable situations. As part of this project we plan to analyze the behavior of the recently developed stochastic shallow convection parametrization (Sakradzija et al. 2016, 2018), focusing on the boundary layer structure. This parametrization couples a stochastic model based on a canonical ensemble of shallow cumuli to the Eddy-Diffusivity Mass-Flux scheme. The moist-convective area fraction is perturbed by subsampling the distribution of sub-grid convective states, yielding fluctuations around the quasi-equilibrium state of shallow cumulus ensemble. Unlike the deterministic approach, this parameterization reproduces the average and higher order statistics of the shallow cumulus case adequately.
Using this parametrization, we aim to implement and verify a detrainment-based cloud cover estimation, since cloud cover plays an important role in the radiation transfer. We further plan to utilize the sub-grid variability of the stochastic parametrization to better estimate the sub-grid microphysical rates. This may have a positive impact on the estimation of precipitation from shallow convection and on the dynamics of the boundary layer. We plan to test this parametrization on real cases, defining its weaknesses and the advantages over the default Tiedke-Bechtold scheme. This product will be verified versus observations in real cases, and it will be compared to idealized LES simulations of cloud cover in BOMEX case. Precipitation from shallow convection can have an important effect on the boundary layer dynamics. In order to estimate the microphysical rates (such as autoconversion, accretion and evaporation rate) a sub-grid PDF of cloud and rainwater content is needed (Larson et al. 2015, Cogan et al. 2014). We plan to utilize the sub-grid variability of the stochastic parametrization to estimate this PDF. If successful, this may have a positive impact on the estimation of precipitation from shallow convection, compared to the current method used in Tiedke-Bechtold scheme.
The last and not less important part of our suggested project deals with testing. Preliminary results (Poliukhov et al., 2018) show that the Tegen aerosol climatology (1997) provides too high aerosol loading. At the same time, the sensitivity studies revealed high impact of aerosol on temperature near surface (Chubarova et al., 2018). The new aerosol dataset proposed by Kinne et al. (2013) is considered to better describe real aerosol loading (Mueller and Träger-Chatterjee, 2014). We plan to verify both the computationally cheap new MACV2 aerosol climatology (Kinne et al. (2013), the new CAMS climatology as well as the prognostic CAMS aerosol against the data of Moscow State University Meteorological Observatory. This dataset contains different aerosol characteristics from its AERONET sun-photometer (Chubarova et al., 2011), hourly visual cloud observations, longwave and shortwave components of net radiation according to the CNR4 KIPP@Zonnen radiometer, etc. (www.momsu.ru). We also wish to test the ICON algorithm in different clear sky conditions (various solar elevation, surface albedo, aerosol loading, etc.). If necessary, another dataset is possible to apply from the AERONET and WRDC archives. We propose to provide verification in clear sky conditions for checking both the radiative code and the different aerosol inputs. We also propose to test the ICON shortwave code with different aerosol loading against other accurate reference radiation transfer (RT) codes.
For studies in cloudy conditions, we plan to use different types of radiative and cloud /aerosol microphysical measurements - performed at the Lindenberg Observatory (Germany) as well as data from the Moscow State University (MSU) observatory and possibly from some other supersites - for checking the new cloud algorithm in the ecRAD radiation scheme. A similar analysis had been applied before with respect to the radiation code of COSMO, which still reveals differences between simulated and measured water vapour profiles at high altitudes of the troposphere (Chubarova et al, 2018, SPIE). Further, it is planned to estimate the sensitivity of the most important meteorological forecast quantities to the new aerosol/cloud input in ecRAD and to assess the change in forecast quality due to all the related measures.
Massive and systematic testing and tuning of the new scheme and its multiple aerosols and cloud options should be performed under different weather conditions for the future decision making for the operational setup in ICON.
Estimated resources: 0.1 FTE per year
In T2(RC)2 new ice particles optical properties were introduced both in COSMO and ICON radiation schemes. Single scattering albedo, asymmetry factor, extinction coefficient and forward scattered fraction based on Fu (1996, 1998, 2007) were calculated. The calculation was based on the single particle data by Fu but instead of using the in-situ particle size distributions (PSD) used by Fu, we generated 7500 PSDs based on the generalized gamma distribution function allowing the extension of the effective size validity range from 70 to 300 micrometres. Using these functions and a wise spectral averaging technique we calculated the bulk optical properties as a function of the effective size of the needle shape particles or its aspect ratio in the case of the asymmetry factor parametrization. The parametrization was based on an ansatz of the form: f(x)= Σ i=0 n aixi / Σ i=0 m bixi where f can be each of the optical properties and x is either the effective size or the aspect ratio and it is asymptotically stable even for particles as large as a few centimeters. The mentioned spectral averaging was adapted for both COSMO and ICON spectral intervals and was separately tested and verified extensively. ecRAD scheme also uses Fu's ice parametrization but with different spectral averaging and bulk optical properties parametrizations.
Therefore, we propose the following subtasks:
Implementation of modified ice particles bulk properties parametrizations based of Fu (1996, 1998, 2007) into ecRAD by transferring the formulas into LUT. Analyzing and comparing the results for both parametrizations and their effects on cloud optical depth. Testing the effects of choosing aspect ratio diagnostics vs. effective radius on the asymmetry factor.
Testing both ICON-RRTM and ecRAD Fu based parametrizations using real case studies/Idealized runs including cirrus clouds. Statistical analysis of the radiative sensitivity to the new parametrization using global ICON runs especially focusing on polar areas.
(02.2021, 0.25 FTE, Harel 0.2, Uli 0.05) Implementation of the new bulk parametrizations based on Fu in ecRAD
(02.2021, 0.2 FTE, Pavel 0.15, Alon 0.1) Case studies and sensitivity analysis
FTEs altogether: Harel 0.2, Uli 0.05, Pavel 0.15, Alon 0.1
Estimated resources: 0.5 FTE
In 1993, Hu and Stamnes (HS93) published their work on liquid clouds optics based on Mie calculations. They claimed that cloud optical properties depend mainly on equivalent effective radius and liquid water content, while the details of the cloud droplets size distribution can be neglected. Based on these assumptions, the authors had computed extinction coefficient, asymmetry factor and single scattering- albedo for 74 narrow wavelengths between 0.29 µm – 150 µm; and they parameterized these optical properties for three droplets size bins in the range of 2.5 µm – 60 µm.
Using the HS93-parametrizations, we have already recalculated the mentioned optical properties for COSMO and ICON spectral intervals as part of PP T2(RC)2. With the new parametrizations, we could extend the validity range of their application, including now also larger liquid particles. The current ecRAD implementation, on the other hand, has two alternatives. The first is the parametrization based on Slingo (1989) for shortwave, and the one of Lindner and Li (2000) for longwave (SLL). The second is based on the SOCRATES Suite of Community Radiation Codes (Edwards and Slingo, 1996).
Therefore, we propose the following subtasks:
Implementation of HS93 parametrization in the ecRAD framework by transferring the formulas into LUT. The spectral intervals in ICON and ecRAD should be similar but in case of differences the optical properties should be calculated to match ecRAD spectral bands.
Radiation and surface temperature sensitivity analysis for the use of each of the mentioned methods (SLL, SOCRATES and HS93) under different cloud types conditions. Real case studies performed comparing the methods and evaluating the model performance in different weather situations against ground based measurements. The domain will be discussed at later time.
(02.2021, 0.25 FTE, Harel 0.2, Uli 0.05) Implementation of HS93 parametrization in the ecRAD framework
(02.2021, 0.3 FTE, Alon 0.1, Pavel 0.15, Yoav 0.05) Sensitivity analysis and case studies performed
FTEs altogether: Harel 0.2 FTE, Uli 0.05 FTE, Pavel 0.15 FTE, Alon 0.1 FTE, Yoav 0.05 FTE
Estimated resources: 0.55 FTE
In Segal & Khain, 2006 (SK), cloud droplets number concentration (QNC) at cloud base is determined by aerosols particle size distribution, number concentration, mode radius, geometric standard deviation (σg) and updraft speed at cloud base. The values can be extracted by parametrization or from 4D look-up tables. In the course of PP T2(RC)2 the look-up table version was applied for COSMO. The aerosols input can be either the Tegen AOD climatology using a fixed mode radius and σg (allowing the 4D table to be reduced to a 2D table) and using the AOD value to define an exponential decay for the aerosols number concentration. Alternatively, the CAMS 3D prognostic aerosols mixing ratios are used to define the number concentration of the different aerosols at each grid point.
In this option some averaged "effective" aerosol mode radii and σg are calculated at each point to account for all hydrophilic tracers available in the CAMS forecast. The number concentration of droplets is used for the evaluation of the effective radius of water clouds, which effects the radiation fluxes (the so-called "indirect radiative effect of aerosols"). There are two ways of calculating the droplet effective radius from the droplet number concentration. First is using the formula of liquid water content (LWC) to QNC power law R‚eff= c1(LWC/QNC)c2 where c1,c2 are determined from the particle size distribution used in the microphysical scheme. The second one is valid for sub-grid scale shallow convection clouds where an analytical adiabatic LWC is used to calculate Reff using the formula R‚eff=1.15(LWCad/QNC)1/3 The actual LWC is calculated from Reff using the same formula with the calculated R‚eff and QNC (Khain et al. 2019). The new cloud number concentration is also used in the 1-mom microphysical scheme for auto-conversion of cloud water to precipitation. Using this formula can dramatically effect precipitation compared to the default fixed cloud number concentration. We propose using the same approach in ICON model.
In addition, for phase 2, yet another subtask (3.4) was proposed:
Through a second way to benefit from CAMS forecasted aerosols, we want to evaluate ice nuclei (IN) from forecasted dust fields according to the DeMott (2015) (DM2015) heterogeneous nucleation scheme. By this approach, IN number-concentration is derived from dust number- concentration and temperature in a computationally very cheap and straightforward way. The results are strikingly realistic, as was shown in COSMO PP T2(RC)2s final report.
Therefore, we propose the following subtasks:
Implementation of SK method for cloud droplets activation in ICON 1-mom microphysics scheme using both Tegen/CAMS climatology and also CAMS prognostic aerosols to produce a realistic cloud number concentration compared with the default fixed number concentration.
Implementation of Reff calculation methods from the prognostic cloud number concentration as described above and applying the result in ecRAD radiation scheme.
Sensitivity analysis and documentation of the effects/case studies.
(02.2022, 0.5 FTE, Harel 0.4, Uli 0.1) Implementation of SK method in ICON 1-mom scheme with different aerosols inputs
(02.2022, 0.2 FTE, Pavel 0.2) Implementation of Reff parametrization methods using the calculated prognostic cloud droplets number concentration
(02.2022, 0.35 FTE, Alon 0.1, Pavel 0.1, Harel 0.1, Yoav 0.05) Sensitivity analysis and case studies performed
FTEs altogether: Harel 0.5, Pavel 0.3, Alon 0.1, Uli 0.1, Yoav 0.05
Estimated resources: 1.05 FTE
Subtask 3.4 (in phase 2)
Coupling of CAMS forecasted aerosols to DM2015 ice-nucleation scheme to produce a realistic IN number-concentration compared to the default fixed value.
(08.2023, 0.2 FTE, Harel 0.2) Coupling of CAMS forecasted mixing ratios with SK2006 for QNC calculation.
(08.2023, 0.2 FTE, Pavel 0.2) Implementation of Reff parametrization methods
(08.2023, 0.325 FTE, Pavel 0.1, Harel 0.2, Yoav 0.025) Sensitivity analysis and case studies performed
(08.2023, 0.2 FTE, Harel 0.2) Coupling of CAMS dust forecasted mixing ratios with DM2015 for IN number concentration calculation.
FTEs altogether: Harel 0.6, Pavel 0.3, Yoav 0.025
Estimated resources: 0.925 FTE
The Copernicus Atmosphere Monitoring Service (CAMS) is an ECMWF project, which provides air quality analysis and forecast including aerosols and other tracers. The model is based on IFS but with a coarser spatial resolution. The output forecast lead time is 5 days with an hourly temporal resolution and includes the 3D mixing ratios of 11 aerosols tracers: 3 size bins of mineral dust, 3 size bins of sea salt, sulphate, hydrophobic and hydrophilic tracers of both black carbon and organic matter. These aerosols tracers are used as both direct and indirect effects on radiation in the COSMO model in an experimental version running twice a day in IMS showing positive results.
The same can be achieved in ICON and more easily in ecRAD framework since the model is already coupled to CAMS climatology input so no optical properties adaptation to this version is needed. At first step, the CAMS data will be used only for cases where the used boundary data from IFS has the same 137 model levels as in CAMS.
Therefore, we propose the following subtasks:
Enabling the download of all CAMS aerosols tracers as additional fields or files via MARS archive as part of the mars4icon_smi pre-processing package. Implementing the CAMS 4D fields in the iconremap interpolation tool to allow the interpolation of the aerosols tracers in both space and time. As was done in COSMO, the aerosols can be interpolated as "boundary conditions" but for the entire domain which can be either global or LAM.
Implementing the CAMS 4D fields in ecRAD framework as a direct effect for the radiation as it was done with CAMS climatology. These fields will also be used for the indirect effect via cloud formation as part of task 3.
Sensitivity analysis and documentation of the effects on case studies.
(02.2022, 0.1 FTE, Harel 0.05, Daniel R. 0.05) Enabling the download of all CAMS fields from MARS (mars4icon_smi package), implementing the CAMS 4D fields in the iconremap interpolation tool.
(02.2022, 0.2 FTE, Harel 0.2) Implementation of the prognostic CAMS aerosols in ecRAD
(02.2022, 0.2 FTE, Alon 0.1, Pavel 0.1) Sensitivity analysis and case studies performed
FTEs altogether: Harel 0.25, Pavel 0.1, Alon 0.1, Daniel R. 0.05
Estimated resources: 0.5 FTE
Based on the combination of short-range forecast and reanalysis (CAMSiRA) of the CAMS model, a new monthly aerosol climatology for the period 2003-2014 was developed at ECMWF (Bozzo et al., ECMWF 801 Technical report, 2017). The new climatology suggests a different AOD 2D distribution of tracers mentioned in the previous sections and also different vertical distributions compared with previous climatologies. The new application showed improvements in both SW and LW spectral intervals mostly for local biases compared to the previous operational Tegen (1997) climatology that was used in IFS model.
Therefore, we propose the following subtasks:
Implementation of CAMS-based aerosol climatology as an external file input for ICON model. Any further work on radiation-aerosols interaction is unnecessary, since this interaction is already described in ecRad.
Testing the new radiation code using the new aerosols CAMS-climatology against satellite observations and ground-based measurements such as BSRN stations.
(02.2021, 0.2 FTE: Poliukhov 0.05, Chubarova 0.05, Shatunova 0.1) Implementing the CAMS-based aerosol climatology fields as an external parameters file for ICON model.
(02.2022, 0.35 FTE: Poliukhov 0.1, Chubarova 0.1, Shatunova 0.1, Khlestova 0.05) Testing and documentation of the effects compared with satellite observations and ground based easements.
FTEs altogether: Poliukhov 0.15, Chubarova 0.15, Shatunova 0.2, Khlestova 0.05
Estimated resources: 0.55 FTE
ICON has an option for a 2-dimensional prognostic aerosol scheme. For this, the horizontal wind components are vertically averaged (v‚H, j) R‚eff and used to calculate advection of the 2-dimensional aerosol optical depth ψj(x,y).
In addition, the scheme contains also simple source functions Se,j (emissions) and sink terms Sw,j (wet deposition). For the aerosol optical depth a relaxation to its climatological ψdim,j value on a certain time scale τdim,j is performed. The following equation is solved for each of the five aerosol optical depth species j:
∂ψj(x,y) ∂t = v‚H, j ∇ψj(x,y)+Se, j+Sw, j+ (fdim, jψdim,j(x,y)-ψj(x,y)) τdim, j
With the tuning factor fdim,j, the target value of the climatological AOD can be modified. The aim of this task is to improve the representation of the source functions Se,j. Consequently, the subtasks can be structured according to the different aerosol species.
Subtask 6.1 Mineral dust optical depth source function
Based on brittle fragmentation theory, Kok (2011a) predicts that the size distribution of emitted mineral dust particles does not depend on the wind speed (Kok, 2011b) which is in good agreement with field measurements. This is in contrast to commonly used mineral dust emission parameterizations that predict such a dependency based on wind tunnel measurements. The prediction of this dependency makes commonly used emission parameterizations like the one used for ICON-ART (Rieger et al., 2017) computationally expensive. Following their previous studies, Kok et al. (2014) present a computationally efficient parameterization for mineral dust emission fluxes. This parameterization will be implemented as a source function for mineral dust. Mie calculations will be performed to derive the dust optical depth source function Se,dist.
Subtask 6.2 Sea salt optical depth source function
In their review publication, Grythe et al. (2014) summarize the state of the sea spray aerosol (SSA) source functions. They compare 21 commonly used parameterizations to different measurement datasets representing the various conditions under which sea salt aerosol particles are emitted globally. In general, the source functions depend on 10-m wind speed, sea surface temperature and, in some cases, on the salinity of the ocean. In this subtask, one of these source functions will be chosen to be implemented into ICON. The decision which function will be implemented will be based on computational efficiency and the quality of the results. Using Mie calculations, a SSA optical depth source function Se,ssa can be derived.
Subtask 6.3 Anthropogenic optical depth source function
Finding an adequate source function for aerosol optical depth from anthropogenic sources will be a major challenge due to the multitude of complex non-linear processes involved in the formation of anthropogenic aerosol particles. These processes include the emission of primary particles and gaseous precursors, atmospheric chemistry, secondary aerosol formation and the interaction with organic compounds. The aim is to find an anthropogenic aerosol optical depth source function Se,anth based on parameters which are available within ICON like land use fraction (e.g., urban, agricultural, etc.) in combination with a daily cycle. Alternatively, the usage of an emission climatology will be explored.
(02.2022, 0.2 FTE, Rieger 0.2) Mineral dust optical depth source function
(02.2022, 0.2 FTE, Rieger 0.2) Sea salt optical depth source function
(02.2022, 0.5 FTE, Rieger 0.5) Anthropogenic optical depth source function
FTEs altogether: Daniel Rieger 0.9 FTE (0.45 per year)
Estimated resources: 0.9 FTE
Subtask 7.1 (cancelled in phase 2)
Since aerosol properties have a noticeable effect on air temperature, their adequate description will improve the quality of weather forecast. The cheapest way to do is to use the state of the art MACv2 aerosol climatology (Kinne et al., 2013) which has been developed recently at Max Planck Institute. This modern aerosol climatology will be implemented in ICON radiation scheme as an alternative to the current default scheme which uses the Tegen (1997) climatology.
We will test its efficiency by comparisons with the prognostic CAMS aerosol as well as with the new CAMS aerosol climatology. For clear sky conditions, the accurate comparison of ICON radiative simulations will be made in the off-line regime for the same statistics with benchmark Monte-Carlo simulations.
(02.2021, 0.3 FTE: Poliukhov 0.15, Chubarova 0.05, Rivin 0.05, Khlestova 0.05) The implementation of Kinne MACv2 aerosol climatology in the ICON model.
(02.2022, 0.25 FTE: Poliukhov 0.1, Chubarova 0.1, Rivin 0.05) The results of intercomparisons of different aerosol ICON simulations with the accurate off-line model simulations in clear sky conditions
FTEs altogether: Poliukhov 0.25, Chubarova 0.15, Rivin 0.1, Khlestova 0.05
Estimated resources: 0.55 FTE
The verification of the output results (aerosols, radiation, main meteorological parameters) after application of the updated Kinne (MACv2) climatology and CAMS aerosol fields from the ECMWF as part of the Copernicus Program in clear sky conditions against different aerosol sites in different meteorological conditions. The baseline site for comparisons is the dataset of Meteorological Observatory of Moscow State University (37.5E, 55.7N) for different geophysical conditions (solar elevation, aerosol properties, surface albedo, and trace gases amount (water vapor, ozone). The main instruments, which will be used, are AERONET CE-318 sun photometer (aerosol properties, column water vapor), Kipp@Zonen net radiometer CNR-4 (upward and downward longwave, shortwave fluxes, and surface albedo), standard meteorological observations (temperature, ground water vapor pressure, and relative humidity), and satellite retrievals (total ozone). The period for the comparisons will cover the previous years since 2018 and the future observations. The total number of cases will be about 15-25 depending on the availability of stable (at least 4-5 hours) clear sky conditions both in model output and real observations.
For cloudy conditions, we will test the new ecRAD radiation scheme (including the account for non-direct aerosol effects) against the cloud and radiation data of the Lindenberg and Moscow Observatories. For this purpose, we plan to use different, accurate measurements of shortwave and longwave radiation where the selected cases will cover preferably the range of one full year.
Further, we want to evaluate the sensitivity of important forecasted meteorological parameters (air temperature, for example, and some others) to the aerosol/cloud characteristics applied in the radiative scheme. Finally, we try to assess the influence of the implemented modifications (including the consideration of various aerosol input) on the absolute and relative forecast error.
(02.2021, 0.2 FTE: Poliukhov 0.05, Chubarova 0.05, Rivin 0.05, Khlestova 0.05) The results of intercomparisons of different aerosol ICON simulations with the accurate experimental measurements in clear sky conditions.
(02.2022, 0.4 FTE: Khlestova 0.15, Chubarova 0.05, Shatunova 0.15, Rivin 0.05) The assessment of the absolute and relative deviations between the forecasted and observed meteorological parameters due to new cloud scheme and different aerosol inputs.
(02.2022, 0.5 FTE: Poliukhov 0.1, Khlestova 0.1, Chubarova 0.1, Shatunova 0.15, Rivin 0.05) The results of intercomparison of different aerosol ICON simulations with the accurate experimental measurements in cloudy conditions. The assessment of the accuracy of implementation of new aerosol climatology to radiation fields and several meteorological parameters in the ICON model.
FTEs altogether: Poliukhov 0.15, Chubarova 0.2, Rivin 0.15, Shatunova 0.3, Khlestova 0.3
Estimated resources: 1.1 FTE
This task deals with the shallow convection parametrization and its effect on the radiation. This parametrization effects the cloud cover estimation especially for partial cloudiness conditions and therefore dramatically effects the atmosphere opacity and the total energy budget. Recently, a stochastic parametrization of shallow cumulus was developed and implemented in a separate branch of ICON (Sakradzija et al. 2016, 2018). At the moment, the scheme has not shown benefits and had some stability issues compared to the default scheme.
Nevertheless, further development is still ongoing in DWD and we plan to test this parametrization on real cases. The detrainment approach will be implemented in this parametrization to deduce the cloud cover of shallow convection. This product will be verified versus observations and compared to idealized LES simulations. We plan to utilize the sub-grid variability of the stochastic parametrization to estimate the microphysical rates, which may have a positive impact on the estimation of precipitation from shallow convection.
Therefore, we propose the following subtasks:
Subtask 9.1 Testing the stochastic shallow convection parametrization
Stochastic shallow convection parametrization will be verified against observations.
We plan to focus on profiles verification against soundings. The exact verification domains and periods will be defined depending on the available computer resources. This will allow us to define the weaknesses and the advantages of the stochastic parametrization over the default Tiedke-Bechtold scheme.
Subtask 9.2 Estimating the cloud cover from the stochastic shallow convection parametrization
Similarly to the Tiedke-Bechtold scheme, the detrainment approach will be implemented in the stochastic shallow convection parametrization to deduce the cloud cover of shallow convection. This product will be verified versus available observations (still not defined) in real cases. Again, the exact verification domains and periods will be defined depending on the available computer resources. In addition, the scheme will be simulated in idealized (preferably BOMEX or RICO) fair weather situations, and compared to idealized LES simulations of cloud cover. Analyze the two cloud cover parameterization in ICON, tune and to determine and the best scheme for different geographical domains.
(02.2021, FTE: Pavel 0.2) Testing the stochastic shallow convection parametrization
(02.2021, FTE: Pavel 0.2) Estimating the cloud cover from the stochastic shallow convection parametrization
FTEs altogether: Pavel 0.4
Estimated resources: 0.4 FTE
Following the project, the cloud-radiation scheme will include many dependencies whose contributions will be described by tens of new parameters. Besides, different options will be activated using several logical switches. First, we will perform sensitivity experiments to identify the parameters and logical switches, which have particularly high influence on the radiative fluxes in the model. Second, several most important logical switches combinations (cloud-radiation scheme configurations) will be selected, and the relevant parameters for each of the combinations will be tuned against global radiation data from surface stations and satellite data. The exact verification domains and periods will be defined depending on the available computer resources. For example, a full parameter tuning can be performed as a first step on the ICON-LAM scale (i.e. ICON-D2 domain) and then can be tested on a global scale when only the aerosols climatology is modified.
Therefore, we propose the following subtasks:
Subtask 10.1 Parameters sensitivity
At first stage, we will perform sensitivity experiments to identify the parameters and logical switches which have particularly high influence on the radiative fluxes in the model. The exact simulations domains and periods for these sensitivity tests will be defined depending on the available computer resources.
Subtask 10.2 Tuning the cloud-radiation scheme
At the second stage, several most important logical switches combinations of the cloud-radiation scheme will be selected, and the relevant parameters for each of the combinations will be tuned against global radiation data from surface stations and satellite data. The exact verification domains and periods will be defined depending on the available computer resources.
(02.2022, 0.15 FTE, Pavel 0.15) Parameters sensitivity
(02.2022, 0.15 FTE, Pavel 0.15) Tuning the cloud-radiation scheme
FTEs altogether: Pavel 0.3
Estimated resources: 0.3 FTE
In this second phase, we would like to propose a new task. The applied microphysical scheme significantly influences the radiation scheme. First, it determines the development and the lifetime of a cloud. Second, it determines its optical properties via the liquid water content and the effective radius.
Spectral Bin Microphysics (SBM) is a state-of-the-art microphysical scheme. It is included in several models worldwide, such as WRF, SAM, HUCM, JMA-NHM and the Goddard Cumulus Ensemble model. Due to the significant computer resources needed to run the scheme, it is usually applied for research purposes and the analysis of test cases. The implementation of SBM will provide similar possibilities in ICON in order to identify the pros and cons of the existing microphysical schemes with special regard to their impact on the simulation of radiation transport. We propose the following subtasks:
Subtask 11.1 Inclusion of the warm-phase SBM in ICON
At the first stage, we will include the warm-phase SBM. The schemes performance and quality will be analysed and compared with the already present 2M-scheme using idealized test cases.
Subtask 11.2 Inclusion of the mixed-phase SBM in ICON
At the second stage, we will include the mixed-phase SBM in the same manner as for the warm phase SBM before.
(02.2022, 0.2 FTE, Pavel 0.2) Inclusion of the warm-phase SBM in ICON
(02.2023, 0.3 FTE, Pavel 0.3) Inclusion of the mixed-phase SBM in ICON
FTEs altogether: Pavel 0.5
Estimated resources: 0.5 FTE
The work on ecRAD radiation scheme should be communicated with ECMWF and should be open for their use in the future. An open discussion with ECMWF scientists should be done during developments.
|Ulrich Blahak||0.1 FTEs||0.1 FTEs|
|Daniel Rieger||0.45 FTEs||0.5 FTEs|
|Pavel Khain||0.7 FTEs||0.7 FTEs|
|Harel Muskatel||0.7 FTEs||0.65 FTEs|
|Alon Shtivelman||0.2 FTEs||0.2 FTEs|
|Yoav Levi||0.05 FTEs||0.05 FTEs|
|Poliukhov||0.25 FTEs||0.3 FTEs|
|Chubarova||0.15 FTEs||0.35 FTEs|
|Shatunova||0.1 FTEs||0.4 FTEs|
|Khlestova||0.1 FTE||0.3 FTEs|
|Rivin||0.1 FTEs||0.15 FTEs|
|Total||2.9 FTEs||3.7 FTE|
|Daniel Rieger||0.25 FTEs||0.25 FTEs|
|Pavel Khain||0.30 FTEs||0.5 FTEs|
|Harel Muskatel||.35 FTEs||0.7 FTEs|
|Alon Shtivelman||.1 FTEs||0.1 FTEs|
|Yoav Levi||.025 FTEs||0.05 FTEs|
|Poliukhov||.2 FTEs||0.2 FTEs|
|Chubarova||.1 FTEs||0.25 FTEs|
|Shatunova||.15 FTEs||0.35 FTEs|
|Khlestova||.1 FTEs||0.25 FTEs|
|Rivin||.05 FTEs||0.15 FTEs|
|Total||1.625 FTEs||2.8 FTEs|