Preliminary Investigation of the Potentialities of a Mesoscale Meteorological Model to Reproduce Experimental Statistics of Rain Attenuation on Earth-Space Links

Current spatial resolutions achieved by mesoscale weather forecast models allow them to be used to generate the state of the lowest layers of the atmosphere over areas as small as a few square kilometers which corresponds to the typical size of the tropospheric area crossed by Earth-space links. Furthermore, they allow the evolution of the troposphere to be predicted with a time stamp of five minutes instead of every hour with large-scale weather forecast models which makes them attractive for radio propagation predictions for satellite communication applications. This paper aims at studying the capability ofthe Weather Research andForecast (WRF) model coupledwith an electromagnetic physical model to reproduce rain attenuation statistics for Earth-space paths at Ka-band. To this purpose, one year of propagation measurements collected at 20GHz in different places at midlatitudes in Toulouse and Salon de Provence (France), Spino d’Adda (Italy), Aveiro (Portugal), and Madrid (Spain), at high latitudes in Svalbard (Norway) and at low latitudes in Kourou are used to make comparisons between simulations and measurements. Comparisons between the simulated and the experimental annual statistics considered in this paper provide encouraging results, with a similar accuracy as Recommendation ITU-R P.618–13 for midlatitude European locations and with better accuracy for a high latitude area in Svalbard and for an equatorial location in French Guiana.


Introduction
Nowadays, xed satellite communication systems deliver services everywhere in the world to provide access to the Internet with higher and higher data rates for the consumer to communicate or to use any kinds of multimedia applications at any time. Experiencing increasing success, these satellite systems have been widely developed using available standard technologies at Ku-band for big satellites placed in geostationary orbit. As a consequence of this success, before being deployed new satellite systems shall demonstrate that they will not cause interference to existing ones, through sharing studies controlled by the regulator: the International Telecommunication Union. ese sharing studies impose a lot of constraints on antenna patterns and on power ampli ers to avoid harmful interference, so much so that the complexity of the systems will increase a lot and it will then be more convenient to use a higher frequency band where more spectrum is available, even if the technology is more expensive. erefore, the combined needs for the end users to bene t from higher data rates and the regulatory rules are now resulting in the observed tendency to go to higher frequency bands such as Ka band (20)(21)(22)(23)(24)(25)(26)(27)(28)(29)(30) [1], Q/V band (40-50 GHz) [2] or probably W-band (70-80 GHz) in the future [3,4].
On the other hand, the maturity of the technology allows designers to propose new systems using non-GSO orbits and several projects of satellite constellations have appeared recently, constituted of tens of medium size satellites in MEO orbits [5] or of hundreds or thousands of small satellites in LEO orbits [6]. In addition, high frequency datalinks are needed not only for telecommunications but also to transfer data from Earth observation satellites or from scientific missions in near Earth or in deep space. However, with increasing frequencies, the impact of the Earth's atmosphere will degrade the propagation of transmitted signals due to attenuation (by gases, clouds, and rain), scintillation due to tropospheric turbulence or depolarization due to nonspherical particles on the path [7]. To close link budgets, that is to specify the appropriate margin to be considered to comply with service availability, it is necessary to know the statistics of these propagation effects. For more than 50 years, a lot of prediction methods of propagation impairments have been proposed in the open literature and the reference method to be followed to predict total impairments for Earth-space paths is given in Recommendation ITU-R P.618 [8]. More recently in the last 20 years, channel models have been developed allowing time series of propagation impairments to be synthesized in order to optimize by preliminary simulating the system design and operation through the use of fade mitigation techniques [9] such as power control, adaptive waveform, or site diversity. e reference channel model is described in Recommendation ITU-R P.1853 [10].
Today, statistical prediction methods and channel models have been tested for links established with GEO satellites thanks to propagation campaigns performed with GEO satellites, such as ACTS [11] in the US, Olympus [12], Italsat [13], and Alphasat [14] in Europe, or the ETS series in Japan [15]. However, to upgrade them for non-GSO configurations, only very few propagation experiments have been performed so far and their duration was not sufficient to obtain reliable statistics. Indeed, for a LEO configuration, a single satellite can be seen from a reference Earth station only around 10 to 15 min per day, which represents a number of samples 100 times lower than with a GEO configuration. erefore, if 5 years are needed to obtain reliable propagation statistics with a GEO experiment, the number of years needed to obtain meaningful annual statistics with a LEO satellite would be significantly higher. So other solutions have to be found to develop and test propagation models for non-GSO configurations.
A new approach to replace propagation campaigns could be to use a physical description of the atmosphere and to convert it into propagation attenuation. In a first attempt, the use of the UK Met'Office Unified Model was proposed [16] to generate the atmosphere and then to apply physical models to produce rain attenuation time series. e main drawback of this approach was the poor resolution of the rain attenuation time series: 5 km in space and 1 hour in time. Considering the progress of computational power and storage capacity, a promising solution would be in a first step to use high resolution mesoscale meteorological weather forecast models to generate the spatial-temporal behavior of the troposphere. en, in a second step this evolving volume would be converted into propagation parameters thanks to physical models of the interaction between an electromagnetic wave and the atmospheric medium, such as radiative transfer equation (RTE), turbulence theory, or Mie theory [17,18].
In a recent past, the use of mesoscale models with RTE proved to be appropriate to model clear sky propagation conditions with the MM5 model to perform link budget for deep space links [19] or for radio links with GEO satellites. In parallel, [18,[20][21][22]Outeiral et al., Jeannin et al., and Fayon et al. propose to use this type of approach to model rain attenuation whose variability in space and in time is significantly faster than oxygen and water vapor attenuation. However, if these preliminary results are quite promising, the validation of the method describing the behavior of rain attenuation on the long-term (that is in terms of local Complementary Cumulative Distribution Functions (CCDFs) of rain attenuation) has not been proved yet in the open literature. erefore, the objective of this paper is to show how a mesoscale weather forecast model used in reanalysis mode, coupled with a physical model relying on the Mie theory of scattering can be successfully used to predict long-term CCDFs of rain attenuation. In Section 2, the use of the Weather Research and Forecast (WRF) model coupled with an electromagnetic module (EMM) is described. e propagation data used to test the representativeness of the model are presented in Section 3. Finally, Section 4 shows the comparisons between the simulations and the statistics generated from propagation measurements collected in several types of climates.

Coupling of a Mesoscale Meteorological
Model with an Electromagnetic Module 2.1. Overview of WRF. WRF is a mesoscale numerical weather prediction system which can simulate atmospheric states using real data or idealized conditions [23]. In this paper, the high-resolution numerical weather prediction model WRF-ARW (v3.4) is used in reanalysis mode to generate a 4D description of the atmosphere at high temporal (5 min) and spatial (2 km) resolutions, thus computing data cubes in a horizontal square of approximately 160 × 160 km 2 and using 36 vertical levels of atmosphere up to 20 km height. WRF is initialized using the ECMWF ERA-Interim reanalysis database which provides surface and pressure levels parameters every 6 hours with a horizontal resolution of 0.75 × 0.75°(around 85 × 85 km 2 ) [24]. ERA-Interim data are also used by WRF-ARW as boundary conditions every 6 hours. Atmospheric data such as temperature, pressure, humidity, and wind profiles are downloaded from ERA-Interim for the desired time period and for the desired region around a central point defined by its geographic coordinates (longitude and latitude), as shown in Figure 1. is information is used to provide sufficiently detailed boundary conditions to WRF in order to solve the Navier-Stokes equations that describe the dynamics of the atmosphere and the related thermodynamic equations.
In order to synthesize realistic past states of the atmosphere, simulations are limited to one day with a previous 12-hour spin-up period (the results of which are deleted later on) giving the opportunity for the model to provide stable outputs. In order to reach the 2 × 2 km 2 horizontal 2 International Journal of Antennas and Propagation resolution, three nested domains are used: a coarser domain with a 30 × 30 km 2 resolution on a polar stereographic projection and static fields initialized directly through an interpolation from ERA-Interim, and two successive nests at 6 × 6 km 2 and 2 × 2 km 2 resolution (see Table 1 and Figure 2). As the solving of such equations requires detailed modelling of the various energy fluxes and of the phase changes, the WRF model includes different models to represent energy transfers between atmospheric layers. Among the parameters computed by WRF, the following ones are of interest to calculate rain attenuation: the pressure level P expressed in Pa, the temperature T in K, and the precipitating liquid water specific content Q r in kg/kg (mass of precipitating water per kg of air mass).

Electromagnetic Module (EMM).
At the output of WRF-ARW simulation, a 4D volume of tropospheric parameters is produced that are converted into propagation parameters and more particularly into rain attenuation using an electromagnetic module (EMM) developed at ONERA (see Figure 3). e specific attenuation due to rain, c r in dB/km, is derived from the following equation (25): where σ ext (D) in cm 2 is the extinction cross-section, N(D) in cm −4 is the rain drop size distribution per class of diameter, and D max in cm is the diameter of the largest raindrop (generally around 0.6 or 0.8 cm because raindrops larger than this size tend to break up when colliding with other large raindrops). Assuming rain drops as spherical particles, the extinction cross section σ ext (D) in cm 2 can be computed using Mie's theory [26,27], which inherently assumes spherical rain drops. As larger drops tend to be oblate spheroids, alternative method such as point matching technique would be better suited, in particular when depolarization needs to be studied. e rain drop size distribution N(D) is derived from the expression of the modified gamma distribution [28]: where Λ is the slope parameter and µ is the shape parameter of the modified gamma distribution. e slope parameter Λ depends on the precipitating water content M r in kg/m 3 as where ρ w is the density of water in kg/m 3 (ρ w � 997 kg/m 3 ) and Γ is the gamma distribution function. An example of the impact of µ parameter on rain attenuation is presented in for an equatorial area. e number of input parameters to be fixed can be reduced assuming a Marshall-Palmer raindrop size distribution [29], so that erefore, the only missing information is the precipitating water content M r . It is related to the precipitating water mixing ratio Q r given by WRF for each pixel cube (in kg/kg) [30] through the following formula: where ρ d is the dry air density (in kg/m 3 ) expressed from the dry air pressure P d in Pa and the temperature T in K as where the specific gas constant for dry air Once the specific attenuation due to rain c r is calculated at the considered frequency for each pixel of the d03 domain (horizontal resolution 2 × 2 km 2 ), the c r values are integrated along the path length to compute rain attenuation time series for the considered (azimuth and elevation angles) satellite radio-link until a specific attenuation value of 0 dB/ km is reached.
Azimuth φ and elevation θ angles are calculated from the following equations: where L es , l es , and h are the longitude, the latitude, and the altitude, respectively, amsl (above mean sea level) of the Earth station, L s is the satellite longitude on the geostationary orbit, R E is the mean radius of the Earth, and R s is the satellite altitude. In order to avoid boundary effects, the d03 domain is truncated from 11 pixels at each side of the grid, and the actual size of the d03 domain is therefore 102 × 102 km 2 . For all locations except SvalSat, the station is placed at the centre of the d01 domain; for the extreme case of SvalSat where a International Journal of Antennas and Propagation very low elevation is considered (2.65°), the station is shifted to have the middle of the path length (up to the 0°C isotherm height) placed at the centre of the d01 domain.
For the considered simulation duration, the computation of rain attenuation is performed every five minutes.
e first database is constituted of propagation measurements collected in Southern France by ONERA from the 20.2 GHz beacon of ASTRA-3B satellite located at 23.5°East longitude on the geostationary orbit in two locations: Toulouse and Salon de Provence in 2014 [32]. Toulouse is situated in South-Western France. It is subject to both Atlantic and Mediterranean climatic influences. Salon de Provence is located at 294 km from Toulouse in South-Eastern France and is subject to the Mediterranean climatic influence. e data processing methodology is detailed in the study of Boulanger et al. [31]. e second dataset considered in this analysis is the statistics of excess attenuation collected in Aveiro (Western Portugal) in 2014 by Universidade de Aveiro from the  19.7 GHz beacon of the EUTELSAT Ka-Sat satellite located at 9°E in GEO orbit. As reported in the studies of Amaya et al. and Garcia et al. [33,34], Aveiro is affected by an oceanic climate. Data processing methodology is detailed in the study of Amaya et al. [33]. e third dataset considered in this analysis is the statistics of excess attenuation collected in Madrid (central Spain) from September 2013 to August 2014 by Universidad Politécnica de Madrid from the 19.7 GHz beacon of Ka-Sat. Madrid is surrounded by mountains. It is located at an altitude of 630 m and is affected by a rather dry continental climate [34]. Data processing methodology is detailed in the study of Garcia et al. [34]. e fourth dataset considered in this analysis is the statistics of excess attenuation collected in Spino d'Adda (Northern Italy) in 1994 by Politecnico di Milano from the 18.7 GHz beacon of the ITALSAT F1 satellite located at 13.2°E in GEO orbit. As reported in the study of Riva [13], being near Milano, Spino d'Adda is affected both by the Mediterranean influence and the proximity of the Alps Mountain Range. Data processing methodology is detailed in the study of Riva [13]. e fifth dataset considered in this paper is constituted of statistics of total attenuation collected in a polar region at the satellite teleport (SvalSat) operated by KSAT in the Spitzbergen Island of the Svalbard archipelago in Norway. e measurements were collected by NASA from April 2016 to March 2017, from the 20.2 GHz beacon of the THOR 7 satellite located at 1°W in the GEO orbit and have been processed and analyzed by ONERA. As reported in the study of Houts et al. [36], two main characteristics have to be considered: on the one hand the climate is very dry because Svalbard lies on the polar circle, and on the other hand, the path length is very long because of the elevation angle lower than 3°. Data processing methodology is detailed in the study of Queyrel et al. [37]. e final dataset considered in this paper is the statistics of excess attenuation collected in an equatorial area at the Guiana Space Centre in Kourou (French Guiana) in 2018 by CNES from the 20.2 GHz beacon of the Amazonas 3 GEO satellite located at 61°W, as reported in the study of Boulanger et al. [38], Kourou is in a coastal area, subject to an equatorial climate in a relatively flat area without any tropical forest, with the main influence of the Atlantic Ocean. Data processing methodology is detailed in the study of Boulanger et al. [38].

Comparisons with Statistics Computed from Propagation Measurements
Comparisons with annual statistics are carried out firstly from propagation measurements collected in Toulouse in 2014, secondly from data measured in different locations in Southern Europe in Aveiro, Madrid, Salon de Provence, and Spino d'Adda, and thirdly for more specific climates at high latitude in Svalbard and at low latitude in Kourou.
Since the final temporal resolution obtained with WRF-EMM is 5 minutes, all experimental time series considered in this study have been undersampled to 5 minutes too. As a consequence, some differences may appear between the experimental results presented here and those published in the studies of Riva et al., Boulanger et al., Lacoste et al., and Garcia et al. [13,31,32,34]. All annual statistics of rain attenuation presented in this section are plotted for time percentages between 0.01% and 5%. e lowest limit is chosen to get a sufficient number of samples to achieve meaningful results, i.e., ∼10 points for the 5 min time step. e highest limit is chosen because it corresponds more or less to the probability of occurrence of rain in midlatitude European climates. For Svalbard, as total attenuation is considered, the plots are provided up to 100% of an average year.  Figure 4.

Comparisons with Propagation Measurements
In both examples presented in Figure 4, the model is not able to reproduce accurately the instantaneous temporal variation of rain attenuation observed on the measurements, even if the rainy period is well identified especially in Figure 4(b). is can be explained either by the lower spatial resolution of the model (2 × 2 km 2 ) or in time (one sample every 5 min). Another possible explanation is that the model is able to generate the expected rain events over a large area but not exactly where they are in reality within a few kilometres (as for any forecast model). erefore, the following analysis will focus on the statistical representativeness of the model instead of the instantaneous temporal behavior.
To appreciate the ability of WRF-EMM to reproduce the statistical behavior of the atmospheric propagation channel, one full year of data simulated at 20.2 GHz in the direction of ASTRA-3B is generated with WRF-EMM for the year 2014.
e configuration of the model is the one described in Section 2. e statistics computed from the time series generated with WRF-EMM and obtained from the measurements are presented in Figure 5.
ey are compared to the CCDF predicted with Recommendation ITU-R P.618-13, taking as     International Journal of Antennas and Propagation input the rainfall rate exceeded for 0.01% of an average year predicted by Recommendation ITU-R P.837-7. Figure 5 shows that for this particular year the annual CCDF derived from one year of WRF-EMM simulations is higher than the measured and predicted CCDF for time percentages lower than 0.5%. For the highest time percentages above 0.5%, better results are obtained with WRF-EMM which is able to predict very well the probability to have rain attenuation on the link where Recommendation ITU-R P.618-13 fails due to the shape of the theoretical distribution in Step 10 of this model. Nevertheless, over the range of probability as a whole, the differences between the two approaches are small.

Comparisons with Annual Attenuation Statistics in Aveiro, Madrid, Salon de Provence, and Spino d'Adda.
To broaden the discussion, in this section, the WRF-EMM simulator is tested against the other measured datasets listed in Section 3 in terms of annual statistics.
In Figure 6, the results obtained with the simulations are worse than the predicted ones using Recommendation ITU-R P.618-13 for almost all the considered time percentages, even if the shape of the statistics is better respected. A small bias seems to be present on the simulation that may be due to the coastal situation of this experiment. Indeed, the pixel considered on the one hand from ERA Interim to initialize WRF, and on the other hand as output of the WRF-EMM simulation is located both over the Atlantic Ocean and over lands.
In the case of Madrid which is affected by a continental climate, both the prediction and the simulation fail to reproduce the shape of the CCDF, even if the differences are not so high in absolute values in compliance with Figure 7. For this particular very dry year (between September 2013 and August 2014), WRF-EMM provides better results for percentages lower than 0.1%, whereas Recommendation ITU-R P.618-13 performs better above this time percentage.
For Salon de Provence, differently from Madrid, the measured CCDF is one of the worst cases ever measured at 20 GHz in Europe, which is confirmed by Météo France outlining very strong convective events called "Cevennes episodes" in this area during autumn 2014. Even if the estimate of the probability of rain attenuation and even if the shape of the statistics coming from the simulations performed with WRF-EMM are similar to the observations, further analyses need to be performed to better reproduce the particular strong convective events observed in 2014 in this area. Here, a similar trend as in Salon de Provence is observed, with better results obtained with Recommendation ITU-R P.618-13. However, here again, the simulated results obtained with WRF-EMM are not too far from the predicted ones.

Comparisons with Annual Attenuation Statistics in Specific Climates at Low and High Latitudes.
In this section, specific climatic conditions are considered to reinforce the interest of using high resolution weather forecast models: on the one hand, a low latitude area with an equatorial climate in French Guiana, and on the other hand, a high latitude area with a polar climate in Svalbard.
In a first step, the results of a similar analysis carried out for an equatorial area [39] from one year of propagation measurements collected in Kourou (French Guiana) [38] are presented in Figure 10.
For this equatorial area, quite good results are obtained with WRF-EMM for time percentages higher than 0.05%. In this case, the atmospheric simulator performs much better than Recommendation ITU-R P.618-13 fed with R 0.01 predicted with Rec. ITU-R P.837-7.
In a second step, results are presented for the specific configuration of a very low elevation link (2.65°) at high latitude (78°N) in Svalbard. Considering this very specific geometry, rain attenuation is not the dominant effect to be taken into account, so unlike the results presented in this paper, the analysis has been carried out considering here total attenuation. In this specific configuration, the ray bending effect is taken into account using Recommendation ITU-R P.676-11. e statistics of total attenuation in Svalbard computed from measurements, from the simulations with WRF-EMM and from Recommendations ITU-R P.618-13 for rain attenuation, P.840-8 for cloud attenuation and P.676-11 for gaseous attenuation are presented in Figure 11 for the period April 2016-March 2017 [37].
From the study of Houts et al. [36], here, both WRF-EMM and Recommendation ITU-R P.618-13 predict a CCDF a bit far from the measured one. However, for WRF-EMM simulations, there is an underprediction with respect to the measurements that could be explained by two missing effects: on the one hand by the absence of the impairments due to other hydrometeors than rain in this version of WRF-EMM, and on the other hand, by the specific configuration of this experiment with a very low elevation less than 3°, for which other contributions such as refraction and scintillation may play a more significant role than for usual elevation angles.
Eventually, the results presented in this section are an illustration of the potential of physical simulators based on meteorological models for link configurations for which the prediction method of rain attenuation recommended by ITU-R was not tested during its development (due to lack of corresponding data). For such cases, it can be expected that an atmospheric simulator relying on a high-resolution weather forecast model can provide fair results because of its physical nature.

Quantitative Comparisons with the ITU-R Metric.
In this section, the annual rain attenuation statistics obtained by WRF-EMM model as well as the prediction method of Rec. ITU-R P.618-13 are compared to the experimental data from a quantitative point of view. For this purpose, the ITU-R metric described in Rec. ITU-R P.311-15 is used [35]. is metric is defined by the following step-by-step methodology: -If A m,i > 10 dB: with M, a weighing factor representing the number of years upon which the measured data is based on.
(3) Calculate the final error ϵ ITUE,i : with and In this study, the probability vector p is created with only certain probabilities considered as valid. For each location, the valid probabilities are illustrated on Figures 5-11 by blue vertical lines.
For the purpose of this analysis, both the mean error ϵ ITUE,m and the RMSE ϵ ITUE,s are calculated and are gathered in Table 4.

Conclusion
is paper has investigated the interest of using a highresolution mesoscale weather forecast model coupled with a physical deterministic model to simulate the statistical behavior of rain attenuation affecting satellite communications at Ka-band. More specifically, Version 3.4 of the WRF-ARF weather forecast model initialized with the ECMWF ERA-Interim reanalysis database has been run to generate every five minutes the state of the troposphere crossed by the considered Earth-space links.
e physical deterministic method used to convert the volume of the troposphere into propagation parameters and more particularly into rain attenuation relies on the Mie scattering theory considering a Marshall-Palmer rain drop size distribution.
Unlike other papers published on the use of high-resolution weather forecast models for propagation purposes that provide comparisons about specific rain events in   International Journal of Antennas and Propagation 9 specific locations, this article focuses both on long-term simulations over a full period of one year as well as on several locations at mid latitudes and at high latitudes in Europe and in an equatorial area.
One year of simulations has been carried out to simulate time series of rain attenuation and then to compute the corresponding statistics. Several locations in different climates have been considered where rain attenuation measurements were available over one full year: in Toulouse and Salon de Provence (France), Aveiro (Portugal), Madrid (Spain), Spino d'Adda (Italy), SvalSat (Svalbard, Norway), and Kourou (French Guiana). e comparisons between the simulated and the measured annual statistics provide encouraging results, the accuracy obtained with the simulator being on the same order of magnitude as the ones predicted with the reference prediction method of Recommendation ITU-R P.618-13 for the five midlatitude European locations considered in this paper. e observed significant differences are probably due to attenuation values generated by the model averaged over an entire pixel, which should be lower than the point measurements over the slant path. For the two specific climates at a high latitude in Svalbard and at a low latitude in French Guiana, significantly better results are obtained with the atmospheric numerical simulator based on the WRF high-resolution model, probably because the conventional prediction methods have not been developed and tested for these specific climates and configurations. erefore, at this stage, this paper demonstrated that this type of mesoscale weather forecast models is a promising solution to generate the behavior of the propagation channel, not only in clear sky conditions but also for precipitating systems to assess rain attenuation annual statistics, for a wide variety of climates. Time series generation using WRF models might be useful for predicting rain attenuation in order to develop fade mitigation techniques, so any future endeavour in this direction would be useful. In the future after more intense validation activities, the goal will be to exploit this kind of atmospheric numerical simulator to generate the statistical behavior of the propagation channel where measurements are not available. More precisely, this could be particularly interesting for specific geometrical configurations (such as for MEO or LEO satellite paths), low elevations, specific areas where orography plays a significant role, higher frequencies (such as V or W bands), or specific climatic areas for instance in the tropics.
However, work still needs to be done in order to improve the accuracy of the results. Different axes can be followed to reach this objective. e primary step is to consider a more recent version of the WRF weather forecast model (version 4) and to initialize it with more recent reanalysis databases such as ERA5 from the ECMWF which is characterized by higher spatial and temporal resolutions than ERA Interim (0.25°× 0.25°instead of 0.75°× 0.75°in the space domain and 1 hour instead of 6 hours in the time domain). e analysis presented in this paper was performed by propagation researchers that are not experts of weather forecast models. erefore, the WRF model has been used with the standard configuration recommended by its designers. In the future, the activity on the configuration of the model itself needs to be undertaken, considering for instance the number of domains to be considered in the simulation, the radiative transfer, and cumulus schemes or the microphysics of precipitation systems in terms of rain drop granulometry. More precisely, if the raindrop size distribution of Marshall-Palmer seems to be appropriate at midlatitudes, its applicability to tropical and equatorial regions as well as to high latitudes is more questionable.

Data Availability
Statistical distributions of the data generated with the WRF-EMM atmospheric simulator of ONERA can be accessible on request to laurent.castanet@onera.fr.

Conflicts of Interest
e authors declare that they have no conflicts of interest.