Offshore Wind Energy Resource in the Kingdom of Morocco: Assessment of the Seasonal Potential Variability Based on Satellite Data

: This study provides a ﬁrst estimate of the offshore wind power potential along the Mo-roccan Atlantic shelf based on remotely sensed data. An in-depth knowledge of wind potential characteristics allows assessment of the offshore wind energy project. Based on consistent daily satellite data retrieved from the Advanced Scatterometer (ASCAT) spanning the period from 2008 to 2017, the seasonal wind characteristics were statistically analyzed using the climatological Weibull distribution functions and an assessment of the Moroccan potential coastal wind energy resources was qualitatively analyzed across a range of sites likely to be suitable for possible exploitation. Also, an atlas of wind power density (WPD) at a height of 80 m was provided for the whole Moroccan coast. An examination of the bathymetrical conditions of the study area was carried out since bathymetry is among the primary factors that need to be examined with the wind potential during offshore wind project planning. The results were presented based on the average wind intensity and the prevailing direction, and also the wind power density was shown at monthly, seasonal and interannual time scale. The analysis indicated that the coastal wind regime of the southern area of Morocco has the greatest energy potential, with an average power density which can reach in some places a value around 450 W/m 2 at heights of 10 m and 80 m above sea level (a.s.l) (wind turbine hub height) more particularly in the south of the country.


Introduction
We are currently experiencing a new era of climate change caused potentially by the use of traditional energy power-generating methods such as fossil fuels, which have a high environmental cost, both locally and globally [1][2][3].
A global trend is currently observed in terms of demand and acceptability of efficient and renewable energy.These are already considered as a reliable alternative with regard to the negative impact.
Worldwide, several countries are already resolutely turned towards wind energy, which is the case of China, world leader with installed capacity of 237 GW, the USA by 227 GW, Germany 78.9 GW, India 42.8 GW and Spain with 48.1 GW installed capacity [4].
Offshore wind has emerged as one of the most dynamic technologies in the energy system.In 2018, a total of 4.3 GW of new offshore wind capacity was completed.From 3 GW of offshore wind in operation in 2010, installed capacity expanded to 23 GW in 2018.Annual deployment has increased by nearly 30% per year [5].
So far, in Morocco, no project has been installed or planned in offshore waters and to date no study has been launched on this segment.We hope that this paper contributes to arousing interest in the offshore wind power potential.
The wind resource is generally much greater offshore than over land for Morocco, with estimated potential of 250 GW, equivalent of 10 times the national wind energy potential onshore, of which 6 GW could be installed by 2030 [6], and offshore wind farms can be built with very low negative environmental impact as stipulated by Haddouche [7].Indeed, offshore wind farms are mostly located far from the main electricity consumption areas potentially limiting the benefits of their network services.Nevertheless, many large Moroccan cities are located close to the shoreline and already benefit from wind farm projects carried out or planned inland (Figure 1), which is advantageous for energy storage and can be attractive for deploying offshore power farms and reducing transmission lines.[8].Also, the energy map emphasizes the location of present and future wind farms projects by 2022.
Wind energy resources along the coastline are not known in much detail in Morocco despite the fact that several wind farms are being planned and are mostly erected on land [9,10] with the exception of some works conducted to estimate the wind potential close to the shoreline such as Allouhi et al. [11] who tried to assess the wind energy potential at six coastal locations in the Kingdom of Morocco, namely Al Hoceima, Tetouane, Assila, Essouira, Laayoune and Dakhla based on hourly in situ measurements from 2011 to 2015 at heights of 50 m and 75 m inland but not at the body of the sea.Also, Elkhchine et al. [10] used the Weibull probability distribution function to evaluate the wind energy potential in two regions, Taza and Dakhla, over three years duration.The monthly mean wind speed variation, wind power density and energy output at different heights was also determined.
In this study, we used remote-sensing tools mainly derived from the Advanced Scatterometer (ASCAT) for retrieving time series and determining the spatial and temporal variability of wind speed frequency distribution along the Moroccan Atlantic Ocean along the continental shelf.
The scatterometer provides observations of the ocean surface wind variability over large areas at higher spatial and temporal resolutions from 2008 to 2017 since in situ measurements were not available.Furthermore, the satellite information is not yet demystified or fully exploited [12].Such a long time series is valuable to assess the wind regimes and the corresponding energy quantities over coastal areas witch the most appropriate locations for wind project installation.
The analysis was carried out along the entire Atlantic shoreline with a focus on nine particular sites mostly located in the vicinity of the Cape.First, the data measured over a period of 10 years (2008-2017) was extracted and analyzed.The characteristics of wind such as mean wind speed, wind rose, extreme wind speed of each site were deduced.However, the deduced information was still exploratory for estimating the whole offshore wind resources around the Moroccan coastline until validated by in situ measurement.Despite low fourth degree resolution (0.25 • × 0.25 • ), satellite products alone can be used to reliably assess offshore wind resource [13].
Based on ASCAT products, mathematical modelling of the wind field was carried out to determine the Weibull parameters giving the best fit with measured data to estimate the adequate wind distribution based on the probability density and cumulative distribution functions.Also, an atlas of the spatio-temporal wind power climatology describing the main monthly, seasonal and interannual variability pattern was provided.Results of wind resources assessment around along areas of Morocco were summarized in this paper.
This paper is structured as follows: first we present the study area including renewable resources in Morocco, particularly wind resources.The strategy was to develop the share of these resources in electricity production and the status of wind projects.In Section 3, materials and methodology to assess wind potential in selected sites.The results and discussion are in Section 4 and the conclusions in Section 5.

Study Area
Located in the extreme North-West of the African continent, along the Atlantic and Mediterranean sea between latitudes 20 • N to 36 • N and longitudes 1 • W to 17 • W, Morocco is just 14 km from Europe across the strait of Gibraltar.Thanks to its large latitudinal extension, Morocco benefits from an important Atlantic shore (2934 km) and a very wide continental shelf (Figure 1) which is under the direct forcing of the subtropical anticyclones that surround the North Atlantic Ocean.
Morocco, a non-oil producing country, is over 96% dependent on the outside world for its domestic energy supplies [14], and the energy bill weighs heavily on its economic and financial balances.This dependence left Morocco rather exposed to the volatility of global energy prices.
The Kingdom of Morocco since 2008 has refined its energy strategy to better meet the challenges in terms of supply security, access to energy and environmental protection that aim to simultaneously decrease the dependence on foreign markets and increase the role of clean energy sources [15].We can already see progress after recent years transitioning towards a decarbonisation of the electricity sector.This transition reduced the energy dependency ratio from 98% in 2008 to 93% in 2018.The target is to increase the share of renewable energy to 42% (2 GW solar and 2 GW wind) of installed capacity by 2020 and to 52% by 2030 [15].As cited by Kouskou et al. [16], at the end of 2017, the kingdom had already reached 34% of its capacity production [16].
Morocco is endowed with real advantages, notably its strategic position and the prospects for development of the renewable energy sector are very favorable thanks in particular to: -Huge wind energy potential due to it 3500 km coast line with annual average wind speed exceeding 9 m s −1 at 40 m elevation.Therefore, the estimated total theoretical potential of wind power in Morocco is 25 GW [17].-Solar potential illustrated by 3000 h of sunshine per year and 5 KWh/m 2 /day of irradiation [11,18].
Despite these advantages, at the end of 2018 the installed capacity of renewable sources was 3685 MW and equivalent to 33.8% of total production, including 1215 MW of wind energy, representing 11.3% of the country's electricity production [19].In 2019, the objectives registered a shift, and no new installed power was recorded (Table 1), which will impact the ambitious target of meeting 42% of its energy requirements using renewable resources by 2020.The integrated wind program with 850 MW has been delayed because of the recent reorganisation of energy sector.However, several projects are in progress, and by the end of 2022, around 2000 MW will be installed [16,20].

Remote-Sensing Data Collection
The remotely sensed wind velocity components at 10 m above the sea surface are extracted from the Advanced Scatterometer Satellite (ASCAT) database available at http://apdrc.soest.hawaii.edu/datadoc/ascat.php.The ASCAT data were processed by the National Oceanic and Atmospheric Administration (NOAA) utilizing measurements from ASCAT aboard the EUMETSAT (European Organisation for the Exploitation of Mete-orological Satellites)/METOP (polar orbiting meteorological satellites) satellite at C-band.The data is used over 10 years (2008-2017) with daily chronological period which estimates near-surface wind speed at 10 ma.s.l by measuring the ocean's surface roughness.The data are mapped into 0.25 • × 0.25 • grid latitude/longitude along the Moroccan domain extended from 20-36 • N of latitude and 30-5 • W of longitude.The wind climatology is constructed in order to highlight of the offshore wind along the national coastline.
Although, in order to ensure the reliability of the ASCAT product for the estimation of the wind potential, an attempt was made to compare at daily timescale the in situ data initially recorded three-hourly along a full year (2018) at the Dakhla weather station with the satellite data.

Methods
Statistical analysis was carried out to determine the wind energy potential along the entire Moroccan area with special attention to some given sites (Table 2) and estimate the wind energy output.To describe the statistical distribution of wind speed, various probability functions exist throughout literature for wind regimes.According to [21], Weibull distribution has a very wide range of applicability, with an acceptable accuracy level.This function has the advantage of making it possible to quickly determine the average of annual production of a given wind turbine [22].However, the main limitation of the Weibull probability density function (PDF) is that it does not exhibit exactly the probabilities of low wind speeds [23].The detailed geographic information of the nine sites along the Moroccan coastline are shown in Table 2 and Figure 1.

Weibull Distribution
Weibull modelling is general, It's encompasses exponential distributions or Rayleigh that are only specific cases of this function [23].The function of Weibull can be described by two or three parameters.A Weibull model with three settings is more general and more flexible but the parameters are much more complexes to be determined [24].In order to comply with the standards of the wind industry, we use the two-parameters Weibull probability density function (PDF), and the mathematical expression is as follows [25,26]: where k is the Weibull shape parameter which describes the dispersion of the data, (it indicates a more or less great dissymmetry of the distribution of the wind speeds), c is the scale parameter which is related to the average speed by the function gamma (it characterizes, like the average, the intensity of the wind), and PDF (V) is the probability of observing wind speed V (m s −1 ).
The corresponding function to the PDF is the CDF which stands for Weibull cumulative density function (CDF) of observing wind speed V obtained by integration the Equation (1) [27] given as: The fitted probability density function is characterized by a scale parameter (c) and a shape parameter (k), which may be used to determine the mean wind speed and the wind power density.This was performed for each pixel along the coastline of the study area.
Wind speed is the key indicator for assessing wind potential [21].The expression of the averaged daily climatological wind speed (V avg ) is obtained as follows: and the associated standard of deviation is: For 1 ≤ k ≤ 10 approximation of the following equation is used [28]: where σv is the standard deviation.Scale factor 'c' can be found from the following empirical approximation:

Evaluation of Wind Power Density (WPD)
The evaluation of the wind power density (WPD) is of fundamental importance in assessing wind power project.In our study, WPD was calculated from the probability density and the wind speed at 10 ma.s.l for the nine sites.
Wind power density (E in W/m 2 ) is directly proportional to air density and can be defined from the instantaneous values of V as in [29]: At sea level and under standard conditions (temperature of 25 • C and pressure of 1 atmosphere), the air density (ρ) is about 1.225 kg/m 3 .

Wind Speed Vertical Extrapolation
The ASCAT wind speed is given at the height of 10 m a.s.l, while the wind speed increases with height.The energy extracted by the wind turbine is estimated using wind speed at the turbine hub height.However, there are theoretical expressions used for determining the wind speed profile.The extrapolation is expressed by the logarithmic wind profile [30] that correlates the wind speed at two different heights, for neutral conditions of atmospheric stability in which thermal effects have been neglected by the following equation: in which, v is the speed that must be calculated at hub height h, v 0 is the wind speed measured at height 10 m above sea level, and z 0 is the surface roughness parameter length which is frequently assumed as a value of 0.0002 m [31,32].
The power law was used to extrapolate the Moroccan offshore mean wind speed derived from satellite data to the hub power height of 80 m a.s.l.

Comparison of Scatterometer Winds and In Situ Measurement
In order to ensure the reliability of the ASCAT product for the estimation of wind potential along the Moroccan shoreline, the Scatterometer winds from ASCAT were roughly compared with in situ data recorded over entire year 2018 at the Dakhla meteorological station by using the density plot.
Globally, the wind speeds from Scatterometer potentially matched the in situ wind measurements where the number of collocated Scatterometer and in situ wind measurements was higher (Figure 2).However, some discrepancies could be observed between both products about the regression line.We noticed that the ASCAT product tended to underestimate the wind potential mainly during summer months when the wind is strong.In fact, the goodness of fit strongly depended on the datasets which contained too few samples in our case (data is over a single year with a daily time step), which made it difficult to compare results from both wind products.Second, we compared a geographically localized in-situ data at land with ASCAT data of coarse spatial resolution (0.25 • ) in the body of the sea.Regardless of that, it is very hard to comment on the statistical uncertainties while the adequacy between both products is acceptable (R = 0.82) and followed a similar pattern distribution (Figure 3).

Seasonal Variability
In order to highlight the general pattern of the climatological seasonal wind speed changes at 10 m a.s.l and the behavior of daily-averaged variability along the coastal band, the spatio-temporal plot drawn as Hovmoller diagram was computed (Figure 4) [33].
The processing was carried out based on the average wind velocity intensity over the first three pixels furthest from the coastline where the wind retrieval within 25-75 km (1-3 pixel widths) from the coast line [33] is not land contaminated.Afterwards, the daily zonal (u) and meridional (v) components were separately extracted and used for computing the wind intensity and mapping the wind direction.
The spatio-temporal changes of the offshore wind speed (Figure 4) showed a general southward increase of the wind intensity thus causing a very marked seasonality between north and south along the coast.The resulting seasonal wind pattern showed a windy region with a slight seasonal maximum variability between North and South.Large variability is seen in spring and winter, whereas weak variation are recorded in summer and autumn seasons.As depicted, the northern part (30-36 • N) of the studied area (20-36 • N) shows a more irregular wind field with moderate mean wind speeds occurs mainly during the summer and autumn (2-4 m s −1 ) and low mean intensity during the winter and late spring (3-2 m s −1 ), while the central part (26-33 • N) shows a strongly seasonal wind intensity during summer (~9 m s −1 ).However, the southern part (20-26 • N) has a quite regular, highly intense (more than 10 m s −1 ) and permanent wind which is blowing throughout the year.
It is clear that the summer months have a high wind speed compared with other months.

Wind Roses
The assessment of a region's wind energy consists of determining the distribution of wind speed and direction.The result is a wind rose which indicates the frequency of wind blowing in different directions and suitable locations for the installation of wind turbines.
Wind rose plots are shown in this study, and may contain an additional information from that obtained by the Hovmoller diagram.The climatological daily frequency distributions of wind speed and direction at 10 m a.s.l show the intensity of prevailing wind for the nine selected sites (Figure 5).
On a climatological time scale, wind roses reveals steadily variations in the annual wind behavior.The mainly prevailing wind directions of the Moroccan atmospheric circulation at sea surface over the selected sites are the North (N), North-north-east (NNE) and North-east (NE).Those three directions comprise about 80% of all climatological daily directions.On all the selected sites, it is rare that the wind blows from the south.
The most intense wind is recorded in sites located much further south especially in the Dakhla, Boujdour and Laayoune followed by the sites of Tarfaya, Tan Tan Cap Ghir and Essaouira.
The diagrams indicate also that the most frequent and highest wind direction originates from the North-East reaching a speed in the range of 9-12 m s −1 .Southerly and Westerlies winds are mainly observed at Cap spartel and Safi, and are least frequent for the rest of the sites and give rise to low and moderate speeds.However, it can be estimated that the wind circulation is much more intense where days of calm or very low wind (<3 m s −1 ) [34][35][36] represent only a very small percentage which does not exceed 20% of the cases.In all directions, high speeds are very common.Fairly strong winds (9-12 m s −1 ) is the dominant pattern of North-east winds.Regarding direction, the analysis shows that there is a predominance of wind in the sector from North-Northeast (NNE) to East-Northeast (ENE).The higher wind is recorded in the North-East direction (NE).In addition, this interval is the largest contributor to total energy.Despite a frequency of significant speeds in the North-west (NW) to South-west (SW) sector, the speed mean wind remains low.

Interannual Variability
The space-time diagram of the offshore wind, as computed from daily ASCAT wind (Figure 6), shows the seasonal and interannual variations of the wind regime for the study period at 10 m a.s.l.It highlights the latitudinal variations in wind intensity oscillation and seasonality, as well as in terms of spatio-temporal extent.
The offshore wind speeds show a very marked oscillation from year-to-year over continental Morocco as whole.We have noted also that, usually, the wind pattern repeats itself year-on-year.The interannual variability exhibit evident regional differences.It's more pronounced in the northern area (North of Safi) with lower intensity and in the central part (from Tan Tan to Safi) with peaked intensity occurring mainly at the middle of each year.The smaller interannual variation occurring in southern Morocco (21 • N-26 • N) with the greatest intensity where the wind is blowing permanently.
The whole time period was characterized by a considerable year to year variability and oscillations between periods of strong (2008, 2010, 2014, 2015 and 2016) and weaker wind intensity (2009,2013,2017).

Weibull Distribution
Using the Weibull distribution presented in Section 3.2.1 for the variation of the standard deviation values, the shape (K) and scale (C) parameters for each grid point of the region along the coastal band were performed.This allowed the calculation of average speed, average power density and the cumulative distribution function (CDF) for the nine sites as shown in Figure 4.
The Hovmollers plots (Figures 6 and 7) and histograms (Figure 8) show, respectively, the monthly and seasonal variation of the Weibull parameters and the probability density function for various values of k and c at height 10 m a.s.l for the whole region.It can be depicted that when the value of k increase, the wind speed variation is less marked and the standard deviation is reduced (Figure 7) giving useful information for industry providing a stability index of the wind over the years.The relatively higher scale parameter for the months of March to August with a peak value of 8.97 m s −1 in July at the Dakhla site, and a slight decrease in September (4.27 m s −1 ) and October (4.61 m s −1 ) was observed.We can also notice that it is maximal in summer months (~7 m s −1 ) while the shape factor ranges from 2 January to 5 September.In addition, the shape factor increases considerably in the middle of the year between May and October.Indeed, these maximum values are 4.3 and 3.7 when calculated, respectively, in spring and summer months.It explains that during the warm months, even if the wind is relatively high, it is uniform, constant and stable during this period.
In order to complete our analysis of the wind behavior, the mean wind frequency classification for the nine sites is expressed and displayed (Figure 8).
Frequency analysis of wind speed being expressed directly (as a percentage) or cumulatively highlights the frequency of occurrence, which allows us to know the probability that a speed value is not exceeded.
It can also be observed from Figure 8 that the forms of the PDFs are almost symmetric where the highest wind speed probabilities are in the range of 6-8 m s −1 .
A predominance of the 6-8 m s −1 range is shown which represents 80% of the cases with a minimum for the range of 1-3 m/s (5%).Almost more than half of the values are greater than or equal to 6-7 m s −1 (high wind).
Moderate winds (less than 5 m s −1 ) represent 15% of the cases and fairly strong and strong winds (6-8) 80%.Finally, strong and very strong winds (greater than or equal 10 m s −1 ) represent only 5% of cases.We noted that for the Cap Spartel, Safi and Essaouira sites, 85% of the data are accumulated in the range of 4 m s −1 to 8 m s −1 and the range of the velocity distribution is up to 9 m s −1 .For the Tarfaya site, 85% of the velocity data is is higher than 6 m s −1 and only 10% is in the range of 1 m s −1 to 4 m s −1 .As for the Tan Tan site, 85% of the data are recorded in the range of 5 m s −1 to 8 m s −1 and the speed range extends up to 10 m s −1 .
The most windy sites are the Boujdour, Laayoune and Dakhla sites with 7 m s −1 annual average speed.They have also the most interesting available energy output and thus represent the most favorable sites for using wind energy for electricity production.For more precision and to give a complementary snapshot on what is happening on the Moroccan continental shelf, the seasonal variability of the power density is described by the mean wind power densities via the Hovmoller space-time diagram at a 25 km/1-day spatio-temporal resolution, from 20 • N to 36 • N (Figure 11).The original data are smoothed with an 11-term moving average in time (11 days) and a 3-term moving average in space (75 km).The result shows a high spatio-temporal definition of the power density pattern, with well defined seasonal transitions from north to south.
Regarding the power density parameter presented in Figure 12, more important energetic peaks are noticed at Dakhla and Essaouira (>450 W/m 2 ) followed by Boujdour, Laayoune and Tarfaya with values in the range of 300-350 W/m 2 .For the northern site mainly at Cap Spartel, the values do not exceed 200 W/m 2 .Overall, this seasonality is similar to that found in other studies conducted by Allouhi [9], Elkhchine [19] and Daoudi [37] mainly focused at the Dakhla and Laayoune sites.Despite the discrepancies in terms of the maximum, this was defined in [19] as equal to 565.26 W/m 2 .
As can be seen, the one exception to this pattern is in the Southern area, where WPD prevailed throughout the year as a consequence of the persistent wind.The northern area exhibited pronounced seasonal change, with WPD reaching values slightly below ~250 W/m 2 in summer, declining to 200 W/m 2 through autumn, winter and spring.
Based on the power density classification approach adopted by the National Renewable Energy Laboratory (NREL) to consider whether a site has favorable power potential or not [38], we can conclude that southern Dakhla, Laayoune and Essaouira, with average annual wind power densities of more than 300 W/m 2 are economically viable for the grid-connected wind power systems.However, Cap Ghir and Tan Tan with annual power density average values that approach 200 W/m 2 (172 W/m 2 and 154 W/m 2 , respectively) were only suitable for the stand-alone wind power application.

Considerations about Wind Technology and Bathymetry
Offshore wind power is maturing, and the implementation of projects needs a comprehension of the wind potential and geological challenges posed by a given physical location.After ensuring the wind potential in the study area, we will focus on two parameters that are relevant factors that may affect and need to be considered when investigating the prospect of an offshore wind project: bathymetry and the distance from the shoreline [39].
Indeed, understand the variations in the geology of the sea is a huge determinant of the type of turbines that can operate above.A detailed assessment of ground conditions is given bellow showing that shallow water depth is highly variable from North to South (Figure 13).[37,38] was used to distinguish the possible floating technology resource from the likely reasonable upper economic limit of fixedbottom resource.
For physical characteristics of the study area, It was observed that shallow water extends the continental shelf from about 15 km at Cap Spartel to 75 km at Dakhla.The average width of the continental shelf is some 40 km to 50 km [40].The head of the continental slope is situated at a depth of 200 m making the Moroccan Atlantic coast a potential area for housing all devices of wind energy.
Otherwise, knowing that the more offshore the farms, the greater the energy potential recovered, we have an interest in pushing into the deep limit offshore and adapting our technology to take advantage of powerful offshore wind fields.
Most offshore wind farms are located at distances typically 8 km (>3 nautical miles) to 30 km (~12 nautical miles) from the shoreline and placed in between 8 to 30 m. Mainly, the offshore wind turbines are in general classified as: (i) shallow water foundation (5 m-30 m) water depth, (ii) transitional water foundation (30 m-60 m), and (iii) deep water foundation (more than 60 m depth) [41].
Both these shallow water and transitional offshore wind turbines have fixed foundations whereas the deep water offshore wind turbines are generally floating structures.
Based on that, we can propose a synthetic classification scheme for the different technologies to be deployed for each site according to topographic features of the continental shelf and to the seaward distance at which the shallow water is located (Table 3).The shelf break of the sites of this group (Northern Essaouira and Safi, Laayoune and Boujdour) is up to 12 nautical miles, and the distance offshore of the shallower water is closer to the shoreline.This region can be considered as transitional region between the closer coast area and offshore area.
-Group 3: Very wide continental shelf The continental shelf of this group mostly in the area including Tan Tan and Dakhla is very wide with a very gentle slope and the shallower waters are extended beyond 24 nautical mile up to a distance of 75 km offshore.These sites present a good energetic potential with mean wind speed >8 m s −1 and WPED > 300 W/m 2 at 80 m a.s.l.
Fixed foundation offshore turbines could potentially deployed in up to 60 m of water for the sites of groups 2 and 3.However, for the sites of group 1, floating support structures offer the possibility for deployment in much deeper waters in the sites of this group.

Conclusions
The objective of the study presented is mainly concerned with the evaluation of the available wind energy potential along the Moroccan continental shelf with a focus on particular sites over one decade from 2008 to 2017.The use of ASCAT data has led to a very fine synoptic analysis of energy resources both spatially and temporally where in situ measurements are missing.Based on these data, we presented a spatial continuity evaluation of offshore wind resources with a spatial resolution of 0.25 degrees and daily time step.
The choice of the climatological cycle has proven to be useful to map and analyze the variations over the monthly and seasonal time scale of the wind speed and for investigating the wind potential density.
A detailed statistical study of wind at 10 m and 80 m a.s.l along the Kingdom of Morocco shoreline is presented.The most important outcomes of the study can be summarized as follows: 1.
Wind characteristic analysis for nine sites of Moroccan shoreline was conducted and the seasonal maps and wind roses are presented.

2.
The wind roses and Hovmoller analysis showed the prevailing wind directions for 10 m height over sea level.

3.
Wind distribution was modeled using the Weibull probability function and the Weibull parameters were calculated.The results show the lowest variability of wind mainly on the South.4.
ASCAT data can be used to usefully map wind resources in countries with large extensional areas.

5.
The offshore winds that blow over the Atlantic shelves are a strong asset for the country.Wind resources are suitable for economically attractive electric generation, with three regions of higher productivity, the northern area at Cap Spartel, central area near Essaouira and the southern area in front and south of Dakhla.6.
Regarding the power density parameter, the wind potential is on average 200 W/m 2 for the entire region.This value varies depending on the season and month.The maximum energy density at height of 80 m a.s.l is more than 500 W/m 2 at the Dakhla site followed by Boujdour, Laayoune, Tarfaya and Essaouira with values in the range of 350-450 W/m 2 .For the northern stations, the average values do not exceed 200 W/m 2 .On the other hand, a production peak is observed at the Essaouira site during the summer, although the average annual production does not exceed 300 W/m 2 .7.
The large energy resource in the southern area is attributed to a large wind intensity.These southern areas tend to have a high quantity of resources at wind speeds between 7 m s −1 and 8 m s −1 and wide continental shelves.8.
These regions can potentially benefit from their proximity to harbor structures and planned projects which make them promising regions for possible exploitation of wind energy.
The processed data set provides informative potential with precise spatial and temporal remotely sensed information.Comparative studies should be conducted in order to improve and complete the assessment of the wind energy potential.A combination of different Scatterometer data, and the use of in situ measurements derived from offshore buoys, are highly desirable.

Figure 1 .
Figure 1.Study area showing the major sampled sites superimposed on the bathymetry of the Moroccan continental margin derived from the gridded chart of the European Marine Observation and Data Network (EMODNET) with the dashed isobath line of 200 m highlighting the width of the continental shelf (https://www.emodnet-bathymetry.eu/data-products) [8].Also, the energy map emphasizes the location of present and future wind farms projects by 2022.

Figure 2 .
Figure 2. Density plot and regression line for Scatterometer winds with respect to in situ measurements winds for 2018 at the Dakhla site; the colour bar indicates the number of collocated Scatterometer winds.

Figure 3 .
Figure 3. Temporal wind variability based on daily in-situ measurement and Advanced Scatterometer (ASCAT) product during the year of 2018.

Figure 4 .
Figure 4. Hovmoller diagram of the climatological ASCAT wind intensity with a spatial resolution of 0.25 degrees and daily time step along the entire Moroccan area at 10 m height a.s.l; the dashed lines denote the sites sampled and considered for an assessment of the energy potential.

Figure 5 .
Figure 5. Wind roses diagrams giving measured wind speed frequency and its corresponding wind directions for the nine selected sites; different scales of [0-15%] are used for (g) to (i).

Figure 6 .
Figure 6.Space-time Hovmoller plot, illustrating the seasonal and interannual variability of the wind intensity from January 2008 to December 2017.

Figure 7 .
Figure 7. Weibull parameters, wind speed and the corresponding standard deviation.

Figure 8 .
Figure 8. Climatological comparison of the observed and the predicted wind speed frequencies associated with the cumulative distribution function for the nine sites.

Figure 10 .
Figure 10.Annual average of WPD (W/m 2 ) computed from 3653 daily images for the period 2008-2017; the likely potential areas of future offshore farms are highlighted according to the WPD intensity.

Figure 11 .
Figure 11.Seasonal mean of the power wind energy density.

Figure 12 .
Figure 12.Climatological monthly variation of the WPD along each selected sites.

Figure 13 .
Figure13.Depth profiles of the different sites from the shoreline to the offshore 200 m isobath limit.The red arrows point at the shallowest waters and the distance at which they are found.The depth delimitation at 60 m[37,38] was used to distinguish the possible floating technology resource from the likely reasonable upper economic limit of fixedbottom resource.

Table 1 .
Status of Morocco wind projects by the end of 2020.

Table 2 .
Geographic information of the selected sites.

Table 3 .
Synthetic table of the likely technologies to be deployed for each site off Morocco according to the topographic features of the continental shelf.This group includes sites having a very narrow continental shelf, mainly Cap Spartel and Cap Ghir, where the seaward distance from the shoreline of the shallower waters is very short (<3 NM).