Basic Facts about Numerical Simulations of Atmospheric Composition in the City of Sofia

The atmospheric composition in urban areas is one of the primary tasks in air pollution studies. The research aims to provide a statistically reliable assessment of the atmospheric composition climate of the city of Sofia—typical and extreme features of the special/temporal behavior, annual means, seasonal and diurnal variations. For that purpose, extensive numerical simulations of the atmospheric composition fields in Sofia city have been performed. Three models were chosen as modeling tools. We used WRF as a meteorological pre-processor, CMAQ as a chemical transport model, and SMOKE as the emission pre-processor of Models-3 system. We developed the following conclusions. The daily concentration changes of the two essential air pollution species—nitrogen dioxide (NO2) and fine particle matters (FPRM, particulate matter (PM2.5), which has a diameter between 0 and 2.5 micrometers)—have different magnitudes. Second, the emissions relative contributions to the concentration of different species could be different, varying from 0% to above 100%. The contributions of different emission categories to other species surface concentrations have various diurnal courses. Last, the total concentration change (ΔC) is different for each pollutant. The sign of the contributions of some processes is evident. Still, some may have different signs depending on the type of emissions, weather conditions, or topography.


Introduction
The atmospheric composition in urban areas is one of the primary tasks in air pollution studies. The air pollution climate in urban areas has not been systematically studied yet in Bulgaria, but some air pollution modeling for the city of Sofia had been performed, and air pollution forecast for the city is operationally occurring [1][2][3][4][5]. Recently, extensive studies for long enough simulation periods and reasonable resolution of the atmospheric composition status in Bulgaria have been carried out using up-to-date modeling tools and detailed and reliable input data [6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23][24][25]. The next step in studying the atmospheric composition climate is performing simulations on an urban scale. Research works dealing with this topic include domains with different sizes and scales from meters to kilometers [26][27][28][29][30][31][32][33][34]. The models' simulations in the smallest scales usually deal with sub-urban-sized areas as streets and neighborhoods. The simulations aim to construct an ensemble comprehensive enough to provide a statistically reliable assessment of the atmospheric composition climate of the city of Sofia-typical and extreme features of the special/temporal behavior, annual means, seasonal and diurnal variations, etc. Some evaluations of the contribution of different pollution sources and dynamic and chemical processes to the atmospheric composition of the city of Sofia will be presented in the paper.
Methodology: Extensive numerical simulations of the atmospheric composition fields in Sofia city have been recently performed, and an ensemble, comprehensive enough as to provide a statistically reliable assessment of the atmospheric composition climate, was The TNO inventory with resolution 0.25 • × 0.125 • in 10 SNAP categories for 2010 [39][40][41] is exploited for the territories outside Bulgaria in the mother CMAQs domain. For the Bulgarian domains, the National inventory for 2010 as provided by Bulgarian Executive Environmental Agency is used. All simulations were performed for a 7-year period from 2008 to 2014 with two-way nesting mod on. Special pre-processing procedures were created for introducing temporal profiles and speciation of the emissions [42].
The study is based on a large number of numerical simulations carried out daily for all periods. Different characteristics of the numerically obtained concentration fields and for five emission scenarios, with all the emissions included and reduced of factor 0.8, including emissions from energetics, non-industrial and industrial combustions, and road transport, will be demonstrated in the present paper. Results concerning the contribution of the different emission categories are demonstrated. The air pollution pattern is also formed as a result of the interaction of different processes; thus, knowing the contribution of each one of these processes for different meteorological conditions and given emission spatial configuration and temporal profiles is useful for understanding the atmospheric composition and air pollutants origin and behavior. To analyze the contribution of different dynamic and chemical processes, the CMAQ "integrated process rate analysis" option was applied. The procedure allows the concentration change for each compound to be presented as a sum of the contribution of each one of the processes, which determines the air pollution concentration.

Materials and Methods
The sample size of the 7-year simulation is large and comprehensive enough to allow processing a variety of statistical quantities. We can calculate different statistics for the entire Sofia city domain, as well as for chosen locations-dispersion, the absolute minimum and maximum concentrations, a percent of the averages in the given interval of values, the probability density function, skewness, kurtosis, etc.
We present plots for the main statistical characteristics of the ensemble characteristics of the air pollution, namely seasonal and annual ones for some pollutants NO 2 (nitrogen dioxide) and FPRM (fine particulate matter), averaged for the Sofia city region and in two locations in the city-Orlov most (point situated at the city center) and Bistrica (situated at semi-mountain outskirts). The plots give graphics of mean concentrations, maximum ones for entire simulation (period/ensemble), and curves noted as 10%, 25%, 75%, and 90%, for the 10th, 25th, 75th, and 90th percentile, respectively. The resulting curves suggest that 50% of the cases fall into the 25-75% interval, and 80% of the cases fall into the 10-90% interval. The plotted curves give a good enough idea for the statistical characteristics of the ensemble-dispersion, asymmetry, excess kurtosis, without being shown explicitly. The axes are given in logarithmic scale due to the high absolute minimal concentrations and the impossibility to distinguish the intervals of cases with different concentrations from the absolute ones on one plot. The absolute minimum concentrations are close to zero in most of the cases; thus, they are not shown on the logarithmic plots.
As already stated, the simulated fields ensemble is large enough to allow statistical treatment. In particular, the probability density functions for each of the atmospheric compounds can be calculated, with the respective seasonal and diurnal variations, for each of the points of the simulation grid or averaged over the territory of the city. Knowing the probability density functions means knowing everything about the ensemble. An example of spatial and diurnal variations of the annual ensembles of surface NO 2 and FPRM behaviors ate two locations, the typically urban site "Orlov most" and "Bistritsa", are shown in the paper.
The emission inventory, used in the simulations, includes 10 emission categories (SNAP categories) and allows the evaluation of the contribution of various anthropogenic activities to the overall picture of air pollution in the city of Sofia: 1.
Extraction and distribution of fossil fuels; 6.
Solvent and other product use; 7.
Other mobile sources and machinery; 9.
The used SNAP categorization of emissions, reduced by some factor, allows evaluation of the contribution of road transport, energetic, industrial, and non-industrial combustions to the atmospheric composition in the city. The concentrations for each scenario of reduced SNAP's were calculated for each day of the period, and the relative contribution of the emissions for each of the scenarios was calculated in the following way.
If an arbitrary pollution characteristic (concentration, deposition, process contribution, etc.) for a given grid point, or averaged over chosen domain, obtained with all the emissions accounted for, is denoted by φ, then φ m is the respective characteristic obtained when the emissions form source category m is reduced by a factor of α. In such a case, the quantity ϕ m can be interpreted as the relative (in %) contributions of emission category m to the formation of the characteristic ϕ: More than one selected nomenclature for sources of air pollution (SNAP) category emissions can be reduced by a factor of α, and thus the joint contribution of several or all SNAP categories to the formation of the pollution characteristic φ can be evaluated. Obtained relative source contributions can also be averaged for the entire ensemble, thus providing the "climate" of the emission contributions, in particular the "typical" annual and seasonal contributions.
The reason why the emissions from a given category are reduced by a chosen factor and not simply at zero is that by completely removing the emissions from a given category, we can obtain much smaller concentrations, which may change the rate of some nonlinear chemical reactions. Moreover, the significant reduction of the concentrations may change the compound diffusion through the domain boundaries, which is why it is a general practice for the simulations to evaluate the contribution of emissions from a given category to be performed as applied and not removed, but reduced emissions from this category, which is normally performed in such studies. The reduction of 20% is preferred by many authors in atmospheric pollution studies.
Five emission scenarios will be considered here: Simulations with all the emissions included and with the emissions from all the SNAP categories (SNALL), SNAP categories 1 (energetics-SN1), SNAP categories 2 (non-industrial combustions-SN2), SNAP categories 3 (industrial combustions-SN3) and SNAP categories 7 (road transport-SN7) for Sofia reduced by a factor of 0.8. This makes it possible to evaluate the contribution of all the emissions, as well as the emissions from road transport, energetics, industrial and non-industrial combustions to the atmospheric composition in the city. The relative contribution of the emissions for each scenario was calculated for each day of this 7-year period, and then, by averaging over the ensemble, the typical fields of relative contributions of these emissions to the surface concentrations of each of the compounds were calculated for the four seasons and annually. For all the emission categories, the pattern of the contribution fields is complex, which reflects the emission source configuration, the heterogeneity of topography, land use, and meteorological conditions. In order to demonstrate the emission contribution behavior in a simpler and easy to comprehend way, the respective fields can be averaged over some domain, which makes it possible to follow and compare the diurnal behavior of the respective contributions for different species. Graphics of the diurnal evolution of the "typical" relative contribution annually and seasonal emissions of SNAP categories 1, 2, 3, 7 and all the emissions to the surface concentrations of NO 2 , FPRM, averaged for the territory of Sofia city and for Orlov most and Bistritsa, are shown in figures.
Atmospheric pollution is a result of the interaction of different dynamic and chemical processes. The consideration of the interaction and contribution of these processes provides a possibility for an explanation of the entire picture of the air pollution in Sofia city. The processes that influence the formation of the air pollution patterns are HADV-horizontal advection, ZADV-vertical advection, HDIF-horizontal diffusion, VDIF-vertical diffusion, EMIS-emissions, DDEP-dry deposition, CLDS-cloud processes, CHEM-chemical processes, and AERO-aerosol processes. In the current section, we present mainly results from the high-performance computing simulations, which evaluate the contribution of different dynamic transportation and transformation processes of air pollutants, which form the air pollution climate in Sofia city. The function of the CMAQ model-"integrated process rate analysis" was used for this task-a specific option that gives an opportunity to estimate the role of each one of the former processes in air pollution formation. In that way, the concentration change ∆C i of the i-th pollutant for a given time interval from t to [t+ ∆t], can be present as a sum of the contributions of the different processes that determine the concentration change, i.e., the equation for transport and transformation of pollutants can be written in the form: The figures show the annual and seasonal averaged contributions of the HADV, ZADV, HDIF, VDIF, EMIS, DDEP, CLDS, CHEM, AERO leading to the different air pollutants (NO 2 , FPRM) formations, averaged for some locations of the domain, as well as for the entire one-Sofia city. All the plots demonstrated are for the first model layer. The total change of the concentration (∆C) is also plotted, and we can see its sign-positive or negative and that it has a well pronounced diurnal course. Different contributions with different values dominate the formation of the air pollutants, which is traced on the graphics and which contribution of given processes dominate, in a given time and sign.

Results
The Section 3 is separated into three parts and presents characteristics of the numerically obtained concentration fields, the contribution of different emission sources, and process analysis of the atmospheric composition for the entire domain Sofia city and two locations-Orlov most and Bistritsa.
Graphics of the diurnal evolution of the "typical" relative contribution-annual and seasonal emissions-of SNAP categories 1, 2, 3, 7 and the emissions of the surface concentrations of NO 2 and FPRM, averaged for the territory of Sofia city and for Orlov most and Bistritsa, are shown.
The annual and seasonal averaged contributions of the processes leading to the formation of pollutants (NO 2 , FPRM), averaged for some locations of the domain as well as for the entire Sofia city, are presented. The total change of the concentration (∆C) is also plotted, and we can see its sign-positive or negative-has a well pronounced diurnal course. Different contributions with different values dominate the formation of the air pollutants, which is traced on the graphics-where a given process dominates, at a given time and sign. The values of ∆C differ in different seasons and are defined from the superposition of the contributions of the different processes, which are also with a different signs.

Characteristics of the Numerically Obtained Concentration Fields, Contribution of Different Emission Sources and Different Dynamic and Chemical Processes to the Atmospheric Composition in Sofia
The graphics in Figure 2 show that the average NO 2 concentrations have a wellexpressed diurnal course with a maximum in the early hours and a minimum in the afternoon. The average concentrations are mostly in the 80% interval of cases (between the two green curves 10-90%). The average concentrations are asymmetrically located in the different pieces of a number of case intervals during the day as well as in the different seasons. The seasonal course shows that the absolute maximum of NO 2 concentrations are highest in the autumn and the winter, which is probably due to the bigger frequency of the stable atmospheric stratification cases and the impeded vertical turbulent transport of NO 2 . The averaged ensemble annual and seasonal contributions of the different sources ( Figure 3), leading to the formation of surface NO 2 in Sofia, sheds light from another perspective. The results suggest that the contribution of the different sources varies in each season but with a similar diurnal course. Everywhere during the night, the biggest natural contribution is from emissions of all SNAP categories, followed by one of the road transport (SNALL and SN7), about 30% around midnight during all seasons. The contribution of the other sources in the morning hours has a small peak of around 10%. Figure 4 shows the annually and seasonally averaged contributions of the processes of NO 2 formation averaged for the entire domain. The main contribution for the NO 2 formation has chemical processes with a positive sign. The horizontal advection is positive in the morning and afternoon and becomes negative at noon. The vertical advection is opposite to the horizontal one. The dry deposition and vertical diffusion have a negative contribution, although with smaller magnitudes. We can outline the positive contribution of the emissions at all hours of the day. The contribution of the other processes is almost zero.    The plots in Figure 5 suggest that the average FPRM concentrations have a wellexpressed diurnal and seasonal course and asymmetrical distribution during the day and seasons for different intervals. The average concentrations fall into the interval containing 80% of cases during the warm months and above that for the cold ones, as during the winter, coinciding with the 90% curve. The absolute maximum and average concentrations are highest in the autumn and winter (stable stratification). There is a well-expressed maximum during the early hours for NO 2 and a minimum around noon. The possible reasons for that maximum are the stable atmosphere and intensive road traffic in the early morning. The results for the averaged by ensemble annual and seasonal contributions of any source leading to the formation of surface FPRM for Sofia city are different (Figure 6), suggesting that the contribution of the different sources varies seasonally but has an equal diurnal course. The dominating natural contribution is one of the sources of all SNAP categories (SNALL), about 40%, followed by one of the road transport (SN7), at about 30%. The contribution of the other sources is about 10% in all seasons. The diurnal distribution shows that the different SNAP categories have maximum contribution in the morning hours and afternoon, and minimum around noon and during the night, which correspond to the concentrations graphics.
The annually and seasonally averaged contributions of processes for FPRM formation averaged for the entire Sofia city domain are shown in Figure 7. The main positive contribution has vertical diffusion for the entire day and vertical advection around noon. The contribution of the vertical advection in the morning and afternoon is negative. The horizontal advection is opposite to the vertical one. The dry deposition has a negative contribution with a maximum around noon. The aerosol processes have a negative sign in all seasons except for winter. The processes are more active during the winter in comparison with other seasons. The winter is outlined with a large positive contribution of emissions and decreases to a negative contribution of vertical diffusion. The other processes have almost zero contributions.     Figure 8 shows the annually and seasonally averaged NO 2 concentrations for the Orlov most station. The results on the plots show that the average concentrations have a wellexpressed diurnal course with a maximum in the early morning hours and a minimum in the afternoon. The average concentrations are mostly in the 80% interval of cases (between the two green curves, at 10-90%). The average concentrations are asymmetrically located in the different intervals of a number of cases during the day as well as the seasons. The seasonal course suggests that the absolute maximum NO 2 concentrations are highest in autumn and winter. The last is probably due to the more cases with stable atmospheric stratification and the impeded turbulent exchange of NO 2 in the vertical direction. The averaged contributions from the different sources for NO 2 formation in each season are different but with an equal diurnal course (Figure 9). The main contribution during the night is from the emissions of all SNAP sources (SNALL)-about 80%, followed by traffic (SN7)-about 50%, in spring, summer, and autumn. The second largest contribution in winter is from the industry (SN3). The contribution of all other sources in the morning hours has a peak of about 20%. The contribution of the sources from non-industrial burning (SN2) is the least.

Characteristics of the Numerically Obtained Concentration Fields, Contribution of Different Emission Sources and Different Dynamic and Chemical Processes to the Atmospheric Composition in Orlov Most
The average annual and seasonal contributions of the processes of NO 2 for Orlov most are shown in Figure 10. The main contribution for NO 2 formation is chemical processes, which have a positive sign, as well as a high positive contribution of emissions in winter. Horizontal advection is positive during the entire day, and vertical advection is opposite. Vertical diffusion also has a negative contribution, although with smaller values and a maximum around noon. The contribution of the other processes is at almost zero magnitude.    The annually and seasonally averaged concentrations of the FPRM for the Orlov most station are shown in Figure 11. The curves show that the average concentrations have well-expressed diurnal and seasonal courses and asymmetric location in different intervals of numbers of cases during the day and seasons. The average concentrations fall into the 80% cases interval during the warm months and the above interval for the cold months. The winter average concentrations coincide with the 90% curve. The absolute maximum and average concentrations are highest for autumn and winter (stable stratification). The early morning hours are characterized by well-expressed maximum and minimum around noon. Probable reasons for this maximum are the stable atmosphere and intensive traffic early morning, which probably play a role in the results for the ensemble annual and seasonal average contributions of different sources for the formation of FPRM (Figure 12). The contribution of the different types of sources varies among the seasons but with the same diurnal course. The main contribution is from sources of all SNAPs (SNALL)-about 70%-followed by the one of the road transport (SN7)-about 60%. The contribution from the other sources is about 20% in all seasons. The distribution during the course of the day has maximum contributions in the morning and afternoon and a minimum around noon and during the night. Figure 13 shows the same characteristics as in Figure 10, but for the FPRM formation. The main positive contribution to the FPRM formation is vertical in the afternoon with vertical advection for all hours of the day. The phase of horizontal advection is opposite the vertical, which shows a negative contribution. The dry deposition is negative with a maximum around noon. The aerosol processes during all seasons except for winter have a negative contribution. The processes are more intense during the winter in comparison with the other seasons, which is outlined with a large positive contribution of emissions and which decreases to negative for vertical diffusion. The contribution of the other processes is almost zero.     Figure 14 shows the annually and seasonally averaged NO 2 concentrations for the Bistritsa station. The graphics show that the average concentrations have well-expressed diurnal course with a maximum in the early morning hours and a minimum in the afternoon. The average concentration curves follow the 75% and fall into the 80% interval of cases. The locations of the average concentrations are asymmetric for different pieces of the number of cases interval during the day and for the seasons. The seasonal course suggests that the maximum nitrogen dioxide concentrations are highest in autumn and winter. The ensemble-averaged contributions for surface NO 2 formation at Bistritsa point ( Figure 15) differ for each type of source and in each season, but their diurnal course is almost equal, although not for Orlov most. The contributions are positive, and the largest one in the afternoon and during the night is SNALL, followed by road transport (SN7). The contribution of the other sources in the morning hours peaks at about 20%. The contributions of the sources from energetics (SN1) and non-industrial burning (SN2) reach 15% in the afternoon during the spring and the winter.

Characteristics of the Numerically Obtained Concentration Fields, Contribution of Different Emission Sources and Different Dynamic and Chemical Processes to the Atmospheric Composition in Bistritsa
The averaged annual and seasonal contributions of the processes of NO 2 averaged for Bistritsa are shown in Figure 16. The main contribution for NO 2 formation is chemical processes, which have a positive sign. Horizontal advection is positive in morning hours and negative at noon, with a peak for all seasons, and vertical advection is opposite. Vertical diffusion also has a negative contribution, although with smaller values and a maximum around noon. The contribution of the other processes is of almost zero magnitude.    The annually and seasonally averaged concentrations of the FPRM for the Bistritsa station are shown in Figure 17. The plots show that the average concentrations have a well-expressed diurnal and seasonal course and are asymmetrically located in different intervals of numbers of cases during the day and seasons. The average concentrations fall into the 80% cases interval during the warm months and above it for the cold months. The winter average concentrations coincide with the 90% curve. The absolute maximum and average concentrations are highest for autumn and winter. The early morning hours are characterized by a well-expressed maximum for NO 2 , with a minimum around noon, which may be due to the stable atmosphere and intensive traffic in the early morning, which probably influence the patterns of the contributions of different sources for surface FPRM formation around Bistritsa ( Figure 18). The averaged contributions are different for each season. The main contribution is naturally from sources of all SNAP categories (SNALL)-about 40%-followed by road transport (SN7)-about 30%. In the afternoon, the sources from combustion in the production and transformation of energy (SN1) have a maximum of about 20%. The contribution from other sources is about 10% for all seasons. The daily distribution suggests that SNAPs have a maximum in the morning and afternoon and a minimum around noon and night. Figure 19 shows the same characteristics as in Figure 16, but for FPRM formation. The main positive contribution to the FPRM formation is vertical diffusion for all hours of the day and vertical advection around noon. The phase of horizontal advection is opposite to the vertical one, which has a negative contribution. The dry deposition is negative with a maximum around noon. The processes are more intense during the winter in comparison with other seasons. The contribution of the other processes is almost zero.

Conclusions
Statistical processing was performed, and the probability density function was calculated for each of the atmospheric compounds with the corresponding seasonal and diurnal fluctuations for each of the points of the grid or the territory of the city. Two points were selected: a mountain type "Bistritsa" and a typical city "Orlov most." Spatial and diurnal variations of the annual ensemble of NO 2 and FPRM are considered. The minimum and maximum curves and the 25%, 75%, 10%, and 90% probability curves are shown. The curves show the imaginary concentrations for which the probability of each simulation is less than 25%, 75%, 10%, and 90%, respectively.
The concentrations of NO 2 , FPRM of Bistritsa are of the same order as Orlov Most. The NO 2 and FPRM concentrations around noon are at a local minimum and are larger during morning hours due to the combination of factors, including traffic and atmospheric stability. The NO 2 concentrations are at a minimum around noon, probably because of the more intense turbulent mixing and the slope effect. The FPRM concentrations do not have significant diurnal variations. The ensemble behavior of NO 2 and FPRM is significantly asymmetric for both selected sites.
For all emission categories, the pattern of the contribution fields is complex, which reflects the emission source configuration and the heterogeneity of topography, land use, and meteorological conditions. Plots of this kind can give a good qualitative impression of the spatial complexity of emission contributions. In order to demonstrate emission contribution behavior in a simpler and easier way, the respective fields can be averaged over some domain, which makes it possible to follow and compare the diurnal behavior of the respective contributions of different species. The results presented in the paper are a first glance at the atmospheric composition status in urban areas; thus, few decisive conclusions can be made at this stage of the study. Different emissions, relative to contribution of the concentration of different species, can be different, varying from 80% to above 100%. The contributions of different emission categories to different species surface concentrations have different diurnal courses. For all of the pollutants, the contribution of SNALL is dominant, but this contribution of emissions is less than 100%, which means that part of the concentrations is formed from sources outside the Sofia city, due to transport into the domain. The contribution of all SN7 (road transport) to NO 2 surface concentrations is positive and reaches about 50% around the busiest traffic roads. The SN7 emissions have dominant contributions to the NO 2 and FPRM surface concentrations. The relative contribution of SNALL and SN7 to the formation of FPRM, as well as for NO 2 at the busiest traffic place, have two maximums in the diurnal course.
The estimation of the contribution of emissions from different source categories is valuable information that can be useful for the definition of measures for improving the air quality in Sofia by reducing emissions. Moreover, knowing the diurnal course of the contributions for a specific time period can suggest an optimal emission temporal regime in order to mitigate air pollution for the given episode.
The results produced by the CMAQ "integrated process rate analysis" demonstrate the complex behavior and interaction of the different processes. Further analysis of these processes, their spatial, diurnal, and seasonal variability, and interaction can be helpful for an explanation of the overall picture and origin of the pollution in the considered region. For the entire domain of Sofia, and for each of the selected items, the total concentration change (∆C), leading to a change in a concentration, is determined mainly by a small number of dominating processes that have large values and may have opposite signs and phases. The total concentration change (∆C) is different for each pollutant. The sign of the contributions of some of the processes is obvious, but some may have different signs, depending on the type of emissions as well as weather conditions and topography. In general, it can be concluded that the contributions of different processes have different behaviors and interact in a complex way.
The "integrated process rate analysis" is not often applied in urban air quality simulations; thus, in the case of the present study, it enriches the entire picture of the atmospheric composition climate of the city of Sofia. This ensemble treatment of the contribution of different processes, however, does not give an easy answer to how the processes interact and how exactly they form the air composition. It will probably be more fruitful to consider the process contribution for a given episode together with the specific meteorological conditions for that episode.
The obtained results and the corresponding conclusions made in this paper are in good agreement with the general idea of how urban atmospheric composition is formed; thus, they do not add general knowledge on the subject. The research and applied contribution of the paper are that it presents quantitative estimations specific to the city of Sofia. Such an extensive and comprehensive study for the city of Sofia had not been made before.
The models and the entire procedure used in the present paper can, of course, be applied to other cities. Such studies have been made for many cities. A crucial point, however, is the emission inventory, and the activities for preparing it are specific for each city, depending on the raw data available.

Data Availability Statement:
The input data-large scale meteorological background and the emission inventories are described in the paper. The output from the computer simulations, which, after processing and generalization is the basis of the results, reported in the paper, is, unfortunately, still not publicly available. This is due to the requirements of the projects, which financially support the present study.