Donato Stilla
email: donystilla@hotmail.it

Mehrez Zribi
email: mehrez.zribi@ird.fr

Nazzareno Pierdicca
email: nazzareno.pierdicca@uniroma1.it

Nicolas Baghdadi
email: nicolas.baghdadi@teledetection.fr

Mireille Huc
email: mireille.huc@cesbio.cnes.fr

remote sensing

Keywords: CYGNSS, GNSS-R, ALOS-2, roughness, aerodynamic roughness, desert

The aim of this paper is to assess the potential use of data recorded by the Global Navigation Satellite System Reflectometry (GNSS-R) Cyclone Global Navigation Satellite System (CYGNSS) constellation to characterize desert surface roughness. The study is applied over the Sahara, the largest non-polar desert in the world. This is based on a spatio-temporal analysis of variations in Cyclone Global Navigation Satellite System (CYGNSS) data, expressed as changes in reflectivity (Γ). In general, the reflectivity of each type of land surface (reliefs, dunes, etc.) encountered at the studied site is found to have a high temporal stability. A grid of CYGNSS Γ measurements has been developed, at the relatively fine resolution of 0.03 • × 0.03 • , and the resulting map of average reflectivity, computed over a 2.5-year period, illustrates the potential of CYGNSS data for the characterization of the main types of desert land surface (dunes, reliefs, etc.). A discussion of the relationship between aerodynamic or geometric roughness and CYGNSS reflectivity is proposed. A high correlation is observed between these roughness parameters and reflectivity. The behaviors of the GNSS-R reflectivity and the Advanced Land Observing Satellite-2 (ALOS-2) Synthetic Aperture Radar (SAR) backscattering coefficient are compared and found to be strongly correlated. An aerodynamic roughness (Z 0 ) map of the Sahara is proposed, using four distinct classes of terrain roughness.

Introduction

Desert areas play an essential role in regional climates, due primarily to the transfer of aerosols [START_REF] Gillette | Modeling dust emission caused by wind erosion[END_REF][START_REF] Callot | Geomorphological approach for modeling the surface features over arid environments in a model of dust emission: Application to the Sahara desert[END_REF][START_REF] Marticorena | Surface and aerodynamic roughness in arid and semiarid areas and their relation to radar backscatter coefficient[END_REF]. In this context, numerous studies have mapped the properties of desert surfaces, in particular parameters such as their aerodynamic roughness. [START_REF] Prigent | Estimation of the aerodynamic roughness length in arid and semi-arid regions over the globe with the ERS scatterometer[END_REF] [START_REF] Prigent | Estimation of the aerodynamic roughness length in arid and semi-arid regions over the globe with the ERS scatterometer[END_REF] proposed the mapping of the aerodynamic roughness length (Z0) using data from the European Remote-Sensing Satellite (ERS) scatterometer. Z0 is defined as the height above the displacement plane at which the mean wind becomes zero when extrapolating the logarithmic wind speed profile downward through the surface layer. These authors have shown the ability of this instrument to distinguish between different classes of aerodynamic roughness Z0, which can then be used as input data for aerosol transfer models [START_REF] Marticorena | Modeling the atmospheric dust cycle: 2. Simulation of Saharan dust sources[END_REF]. Simultaneous, high spatial resolution low frequency SAR radar data can also be used to characterize the geometric roughness of surfaces, generally defined by two statistical parameters calculated from the height autocorrelation function: i.e., the Root Mean Square Height (rmsh) (vertical statistical parameter) and correlation length (horizontal statistical parameter) [START_REF] Zribi | Evaluation of a rough soil surface description with ASAR-ENVISAT radar data[END_REF][START_REF] Hegarat-Mascle | Retrieval of elevation by radarclinometry in arid or semi-arid regions[END_REF][START_REF] Baghdadi | A New Empirical Model for Radar Scattering from Bare Soil Surfaces[END_REF][START_REF] Gorrab | Retrieval of Both Soil Moisture and Texture Using TerraSAR-X Images[END_REF]. Roughness information can be exploited to gain insight into sub-surface properties through the strong penetration of electromagnetic waves [START_REF] Baghdadi | Apport de l'imagerie satellitaire radar pour l'exploration géologique en zones arides[END_REF][START_REF] Grandjean | A Volume Scattering Model for Coupled Interpretation of Ground-Penetrating Radar (GPR) and Synthetic Aperture Radar (SAR)[END_REF][START_REF] Paillou | Subsurface imaging in south-central egypt using low-frequency radar: bir safsaf revisited[END_REF].

Remote Sens. 2020, 12, 743 2 of [START_REF] Rius | Feasibility of GNSS-R Ice Sheet Altimetry in Greenland Using TDS-1[END_REF] In this context, over the last 15 years, we have witnessed the arrival of GNSS-R (Global Navigation Satellite System Reflectometry) technology, which initially demonstrated its considerable potential for the characterization of ocean surface properties (waves, etc.) [START_REF] Martin-Neira | A Passive Reflectometry and Interferometry System(PARIS)-Application to ocean altimetry[END_REF][START_REF] Zavorotny | Tutorial on Remote Sensing Using GNSS Bistatic Radar of Opportunity[END_REF][START_REF] Li | First spaceborne phase altimetry over sea ice using TechDemoSat-1 GNSS-R signals[END_REF][START_REF] Rius | Feasibility of GNSS-R Ice Sheet Altimetry in Greenland Using TDS-1[END_REF][START_REF] Li | Measuring Greenland Ice Sheet Melt Using Spaceborne GNSS Reflectometry from TechDemoSat-1[END_REF]. It was then tested over continental surfaces, where it revealed its strong potential for the monitoring of various surface parameters such as the water content of the soil or vegetation biomass [START_REF] Masters | Initial results of land-reflected GPS bistatic radar measurements in SMEX02[END_REF][START_REF] Katzberg | Utilizing calibrated GPS reflected signals to estimate soil reflectivity and dielectric constant: Results from SMEX02[END_REF][START_REF] Rodriguez-Alvarez | Soil Moisture Retrieval Using GNSS-R Techniques: Experimental Results Over a Bare Soil Field[END_REF][START_REF] Egido | Global Navigation Satellite Systems Reflectometry as a Remote Sensing Tool for Agriculture[END_REF][START_REF] Egido | Airborne GNSS-R Polarimetric Measurements for Soil Moisture and Above-Ground Biomass Estimation[END_REF][START_REF] Pei | A Fully Software GNSS-R receiver for Soil Monitoring[END_REF][START_REF] Sanchez | On the Synergy of Airborne GNSS-R and Landsat 8 for Soil Moisture Estimation[END_REF][START_REF] Jia | Estimation of Surface Characteristics Using GNSS LH-Reflected Signals: Land Versus Water[END_REF]. Following various recent in situ and airborne campaigns collecting GNSS-R measurements [START_REF] Egido | Airborne GNSS-R Polarimetric Measurements for Soil Moisture and Above-Ground Biomass Estimation[END_REF][START_REF] Zribi | CAROLS: A New Airborne L-Band Radiometer for Ocean Surface and Land Observations[END_REF][START_REF] Motte | A GNSS-R Dual Polarization Airborne Instrument for Land Surface Monitoring[END_REF], several studies have illustrated their potential for the development of GNSS-R based applications [START_REF] Camps | Sensitivity of GNSS-R Spaceborne Observations to Soil Moisture and Vegetation[END_REF][START_REF] Chew | Demonstrating soil moisture remote sensing with observations from the UK TechDemoSat-1 satellite mission[END_REF][START_REF] Jia | Sensing soil moisture and vegetation using GNSS-R polarimetric measurement[END_REF][START_REF] Carreno-Luengo | Spaceborne GNSS-R from the SMAP Mission: First Assessment of Polarimetric Scatterometry over Land and Cryosphere[END_REF][START_REF] Zribi | Potential Applications of GNSS-R Observations over Agricultural Areas: Results from the GLORI Airborne Campaign[END_REF][START_REF] Camps | Sensitivity of TDS-1 GNSS-R Reflectivity to Soil Moisture: Global and Regional Differences and Impact of Different Spatial Scales[END_REF][START_REF] Chew | Soil Moisture Sensing Using Spaceborne GNSS Reflections: Comparison of CYGNSS Reflectivity to SMAP Soil Moisture[END_REF][START_REF] Clarizia | Analysis of CYGNSS Data for Soil Moisture Retrieval[END_REF][START_REF] Wan | Using CYGNSS Data to Monitor China's Flood Inundation during Typhoon and Extreme Precipitation Events in 2017[END_REF][START_REF] Al-Khaldi | Time-Series Retrieval of Soil Moisture Using CYGNSS[END_REF][START_REF] Calabia | Soil Moisture Content from GNSS Reflectometry Using Dielectric Permittivity from Fresnel Reflection Coefficients[END_REF]. [START_REF] Camps | Sensitivity of GNSS-R Spaceborne Observations to Soil Moisture and Vegetation[END_REF][START_REF] Carreno-Luengo | Spaceborne GNSS-R from the SMAP Mission: First Assessment of Polarimetric Scatterometry over Land and Cryosphere[END_REF][START_REF] Camps | Sensitivity of TDS-1 GNSS-R Reflectivity to Soil Moisture: Global and Regional Differences and Impact of Different Spatial Scales[END_REF] have shown the strong potential of GNSS-R TechDemoSat-1 receiver and reflections collected by the SMAP (Soil Moisture Active Passive) receiver for the mapping of surface conditions, such as soil moisture or vegetation properties. In this context, it has been shown that surface scattering is predominantly coherent over flat land surfaces, if we may neglect the effect of topography [START_REF] Pierdicca | A Simulator of GNSS Reflections from Bare and Vegetated Soils[END_REF].

The CYGNSS mission, developed by the University of Michigan and the Southwest Research Institute and sponsored by National Aeronautics and Space Administration (NASA), was launched in 2017 [START_REF] Ruf | New Ocean Winds Satellite Mission to Probe Hurricanes and Tropical Convection[END_REF]. It was designed to improve hurricane forecasting with measurements of wind speed and to study interactions between air and water surfaces during the tropical cyclone season, including in high precipitation conditions. This satellite system comprises eight micro-satellites, fitted with the Global Positioning System (GPS) receiver payloads. These are referred to as Delay Doppler Map Instruments (DDMI), and their operating principle relies on the reception of the GPS signal powers reflected by the Earth's surface, given as a function of the delay and Doppler shift. The CYGNSS constellation offers high space-time sampling of several geophysical variables, thanks to the simultaneous collection, by each CYGNSS satellite, of the signals transmitted by four GPS satellites and then scattered by the ocean or terrain. Various studies have been proposed to estimate land surface properties, such as soil moisture and vegetation biomass [START_REF] Chew | Soil Moisture Sensing Using Spaceborne GNSS Reflections: Comparison of CYGNSS Reflectivity to SMAP Soil Moisture[END_REF][START_REF] Clarizia | Analysis of CYGNSS Data for Soil Moisture Retrieval[END_REF][START_REF] Wan | Using CYGNSS Data to Monitor China's Flood Inundation during Typhoon and Extreme Precipitation Events in 2017[END_REF][START_REF] Al-Khaldi | Time-Series Retrieval of Soil Moisture Using CYGNSS[END_REF][START_REF] Calabia | Soil Moisture Content from GNSS Reflectometry Using Dielectric Permittivity from Fresnel Reflection Coefficients[END_REF].

The aim of the present study is to analyze, for the first time, the potential of CYGNSS GNSS-R data for desert surface applications. This paper is organized as follows: Section 2 provides a general description of the studied site and the datasets considered in this study; Section 3 describes the methodology used to interpret and map the CYGNSS data recorded over the Sahara; Section 4 focuses on the results, and discusses the relationship between reflectivity and surface roughness (aerodynamic or geometric); finally, our conclusions are provided in Section 5.

Studied Site and Database

Studied Site

The results described and discussed in this paper refer to data recorded over the Sahara Desert in North Africa (Figure 1); the considered areas are situated between 15 • 0'N and 37 • 30'N and between 17 • 0'W and 39 • 30'E, with some specific areas being selected for further analysis. The Sahara contains a large number of heterogeneous environments: sand dunes, flat rocky surfaces, high mountains and volcanic rocks. Furthermore, the Sahara encompasses various distinctly different climate zones, including a Mediterranean climate along its northern extremity, an extremely arid central zone with rare rainfall events, and a southern zone, which receives some rain during the monsoon season.

Database

Satellite Data a) CYGNSS data

The main dataset used for this analysis is the Level 1 (L1) CYGNSS 0 data, version 2.1, obtainable from the Physical Oceanography Distributed Active Archive Center (PODAAC) (https://podaac.jpl.nasa.gov). The L1 data are stored as daily Network Common Data Form (NETCDF) files, in which each day is represented by eight files (one for each CYGNSS platform), with each of these containing a large number of useful variables that can be used for the subsequent processing steps. The main pieces of information are the Delay/Doppler Maps (DDMs) of the analogue scattered power. For each data point, the specular point properties (i.e., the DDM sample corresponding to the point on the surface minimum electromagnetic path length) are taken into account, together with several geometrical parameters relevant to each acquisition (e.g., specular point incidence angle, transmitter and receiver ranges, etc.), and metadata related to the quality of the incoming and reflected GPS signals (e.g., direct Signal to Noise Ratio (SNR), receiver gain, SNR for each DDM etc.). In this research, almost 2.5 years of CYGNSS acquisitions, recorded between March 2017 and July 2019, have been processed and analyzed. In order to filter the data correctly, bitwise quality flags are provided to select the best data, and in the present study, all available quality flags were considered and set to 1 [START_REF] Gleason | Level 1B DDM Calibration Algorithm Theoretical Basis Document[END_REF][START_REF] Gleason | The CYGNSS Level 1 Calibration Algorithm and Error Analysis Based on On-Orbit Measurements[END_REF]. The theoretical footprint of a reflected GNSS signal is ~0.5 × 0.5 km, for the case of a smooth surface and a receiver at the altitude of CYGNSS. The receivers integrate signals over 1 s. During this period of time, the receiver travels ~7 km. Thus, the smallest area of ground from which reflections are received is ~7 × 0.5 km.

The signal over land includes both coherent and incoherent scattering, with a relative weight that is subject to discussion [START_REF] Balakhder | On the Coherency of Ocean and Land Surface Specular Scattering for GNSS-R and Signals of Opportunity Systems[END_REF][START_REF] Camps | Spatial Resolution in GNSS-R Under Coherent Scattering[END_REF]. The main observables related to the two mechanisms are the reflectivity and the Normalized Radar Cross Section, respectively. Here, we refer to the reflectivity ( Γ ), which should then be considered as an "equivalent reflectivity" that is derived from the CYGNSS data as suggested in [START_REF] Clarizia | Analysis of CYGNSS Data for Soil Moisture Retrieval[END_REF] using the following expression:

Γ(𝜃) = (4𝜋) (𝑃 -𝑁)(𝑅 + 𝑅 ) 𝜆 𝐺 𝐺 𝑃 ( 1
)
where 𝑃 is the maximum value of the analog power in the Delay/Doppler Maps (DDM), N is the noise floor related to the DDM, 𝑅 is the receiver-specular point (SP) distance, 𝑅 is the transmitter-

Database

Satellite Data (a) CYGNSS data

The main dataset used for this analysis is the Level 1 (L1) CYGNSS 0 data, version 2.1, obtainable from the Physical Oceanography Distributed Active Archive Center (PODAAC) (https://podaac.jpl. nasa.gov). The L1 data are stored as daily Network Common Data Form (NETCDF) files, in which each day is represented by eight files (one for each CYGNSS platform), with each of these containing a large number of useful variables that can be used for the subsequent processing steps. The main pieces of information are the Delay/Doppler Maps (DDMs) of the analogue scattered power. For each data point, the specular point properties (i.e., the DDM sample corresponding to the point on the surface minimum electromagnetic path length) are taken into account, together with several geometrical parameters relevant to each acquisition (e.g., specular point incidence angle, transmitter and receiver ranges, etc.), and metadata related to the quality of the incoming and reflected GPS signals (e.g., direct Signal to Noise Ratio (SNR), receiver gain, SNR for each DDM etc.). In this research, almost 2.5 years of CYGNSS acquisitions, recorded between March 2017 and July 2019, have been processed and analyzed. In order to filter the data correctly, bitwise quality flags are provided to select the best data, and in the present study, all available quality flags were considered and set to 1 [START_REF] Gleason | Level 1B DDM Calibration Algorithm Theoretical Basis Document[END_REF][START_REF] Gleason | The CYGNSS Level 1 Calibration Algorithm and Error Analysis Based on On-Orbit Measurements[END_REF]. The theoretical footprint of a reflected GNSS signal is ~0.5 × 0.5 km, for the case of a smooth surface and a receiver at the altitude of CYGNSS. The receivers integrate signals over 1 s. During this period of time, the receiver travels ~7 km. Thus, the smallest area of ground from which reflections are received is ~7 × 0.5 km.

The signal over land includes both coherent and incoherent scattering, with a relative weight that is subject to discussion [START_REF] Balakhder | On the Coherency of Ocean and Land Surface Specular Scattering for GNSS-R and Signals of Opportunity Systems[END_REF][START_REF] Camps | Spatial Resolution in GNSS-R Under Coherent Scattering[END_REF]. The main observables related to the two mechanisms are the reflectivity and the Normalized Radar Cross Section, respectively. Here, we refer to the reflectivity (Γ), which should then be considered as an "equivalent reflectivity" that is derived from the CYGNSS data as suggested in [START_REF] Clarizia | Analysis of CYGNSS Data for Soil Moisture Retrieval[END_REF] using the following expression:

Γ(θ) = (4π) 2 (P DDM -N)(R r + R t ) 2 λ 2 G r G t P t (1)
where P DDM is the maximum value of the analog power in the Delay/Doppler Maps (DDM), N is the noise floor related to the DDM, R r is the receiver-specular point (SP) distance, R t is the transmitter-specular point (SP) distance, λ is the wave length, θ is the incidence angle, G r is the receiver antenna gain in the direction of the specular point (SP) and G t P t is the transmitter equivalent isotropically radiated power (EIRP). The noise floor is computed as the mean value of the DDM subset, when the signal is absent (located above the characteristic horseshoe shape of the DDMs) as suggested, for example, in [START_REF] Clarizia | Analysis of CYGNSS Data for Soil Moisture Retrieval[END_REF].

The data are collected in daily data frames, taking into account various empirical filters used to remove noisy data, which could otherwise be expected to degrade the quality of the imagery and subsequent analysis.

The filtering criteria considered for this purpose were the following:

(1) CYGNSS reflectivity values should range between -35 dB and -5 dB;

(2) Incidence angles should range between 0 • and 30 • ;

(3) DDM signal-to-noise ratio (SNR) should be greater than 3 dB; (4) G rx (receiver antenna gain towards SP) should be greater than 5 dB;

(5) Data corresponding to terrain at altitudes greater than 650 m above the WGS84 ellipsoid model of the Earth should be removed.

Altitudes greater than 650 m are not taken into account, in order to ensure that the detected reflections occur within the Delay/Doppler acquisition window, which is reduced and compressed for DDMs released as L1 CYGNSS products [START_REF] Ruf | DDM Data Compression and Decimation Algorithm[END_REF].

(b) ALOS-2 data

The open dataset provided by the ALOS-2 mission, an advanced L-band platform developed by Japan Aerospace Exploration Agency (JAXA), is considered. The data are recorded at a global resolution of 25 m, created by mosaicking SAR images of backscattering coefficients, which were measured by the Phased Array type L-band Synthetic Aperture Radar 2 (PALSAR-2), freely available at the JAXA Earth Observation Research Center (EORC/JAXA) (https://www.eorc.jaxa.jp/ALOS/en/ palsar_fnf/data/index.htm). Each file contains pre-processed information related to the local incidence angle, acquisition date, and the backscattering coefficient for the horizontal-horizontal (HH) and horizontal-vertical (HV) polarizations. The pre-processing includes ortho-rectification (correction of geometric distortion) and slope correction (topographic effects on image intensity). In the present paper, HH backscattering maps from the year 2017 are used, whereas the HV data are not analyzed due to their greater sensitivity to volume rather than surface scattering. As described in the PALSAR-2 handbook [START_REF] Eorc | Alos Data Users Handbook[END_REF], the backscattering maps are not calibrated, and the results are provided in the form of digital counts. In practice, a calibration coefficient is needed to transform the product into a more convenient decimal format, according to the following expression:

σ o | dB = 10 log 10 DN 2 + CF (2)
where σ 0 is the output backscattering coefficient in decibels, DN is the input digital number contained in the dataset, and CF is the calibration factor, which is set to -83 dB [START_REF] Shimada | PALSAR initial calibration and validation results[END_REF].

(c) SRTM

The altitude values, along with quantities derived from the height of the terrain, were retrieved from the Shuttle Radar Topography Mission (SRTM) data, which provides a three arcsec (90 m) gridded, near-global scale digital elevation model (DEM) between 56 • S and 60 • N, produced by NASA-Jet Propulsion Laboratory (JPL) and available online at: http://srtm.csi.cgiar.org/srtmdata/.

Ground Measurements of Aerodynamic Roughness

In this study, we made use of ground measurements of aerodynamic roughness, measured over the Sahara desert, provided in the study of [START_REF] Marticorena | Modeling the atmospheric dust cycle: 2. Simulation of Saharan dust sources[END_REF] [START_REF] Marticorena | Modeling the atmospheric dust cycle: 2. Simulation of Saharan dust sources[END_REF]. As listed in Table 1, 23 areas are considered, corresponding to a wide range of surface conditions (dunes, smooth surfaces, reliefs etc.). The same database was also used in the past to develop aerodynamic roughness maps, based on data recorded by the ERS scatterometer [START_REF] Prigent | Estimation of the aerodynamic roughness length in arid and semi-arid regions over the globe with the ERS scatterometer[END_REF].

Methodology

This section proposes a methodology to process CYGNSS data over desert surfaces.

CYGNSS Data Normalisation

Several studies have highlighted the influence of incidence angle on reflectivity, and in the present study, the role of this effect needs to be assessed for different types of surface topography. In Figure 2, the computed reflectivity derived from Equation ( 1) is plotted as a function of the incidence angle for two different areas with the whole temporal dataset. These test areas are described in Table 2.

present study, the role of this effect needs to be assessed for different types of surface topography. In Figure 2, the computed reflectivity derived from Equation ( 1) is plotted as a function of the incidence angle for two different areas with the whole temporal dataset. These test areas are described in Table 2.

It can be seen that the incidence angle has a very limited effect for incidence angle values between 2° and 30°.

As shown in Figure 2, a simple linear normalization is thus proposed, referring to an incidence angle of 20°, in order to account for the small influence of the incidence angle on the reflectivity. We use the following expression:

Γ(20°) = Γ(θ) -A . (θ -20°) (3)
where the parameter A is the slope of the retrieved least squares linear relationship between reflectivity and incidence angle for a given type of surface.

In our context, with a limited time series of CYGNSS data, the azimuth angle effect is not considered in the proposed data analysis. This effect is negligible for land surfaces, particularly after data temporal averaging.

Regions Coordinates

Great Eastern Erg-Algeria (Algeria dunes) 31

• 51 N• • • 7 • 33 E 31 • 45 N• • • 7 • 39 E Great Libyan Desert-Egypt (Egypt dunes) 28 • N• • • 25 • 36 E 27 • 54 N• • • 25 • 42 E Kufrah Basin-Libya 23 • 54 N• • • 24 • 9 E 23 • 48 N• • • 24 • 15 E Tamanarasset (Mali) 23 • 49 N• • • 4 • 46 W 23 • 43 N• • • 4 • 40 W
It can be seen that the incidence angle has a very limited effect for incidence angle values between 2 • and 30 • .

As shown in Figure 2, a simple linear normalization is thus proposed, referring to an incidence angle of 20 • , in order to account for the small influence of the incidence angle on the reflectivity. We use the following expression:

Γ(20 • ) = Γ(θ) -A . (θ -20 • ) (3)
where the parameter A is the slope of the retrieved least squares linear relationship between reflectivity and incidence angle for a given type of surface.

In our context, with a limited time series of CYGNSS data, the azimuth angle effect is not considered in the proposed data analysis. This effect is negligible for land surfaces, particularly after data temporal averaging.

Temporal Stability of CYGNSS Data

In the case of desert areas, remotely sensed signals, and in particular microwave measurements, are in general found to be very stable [START_REF] Hegarat-Mascle | Retrieval of elevation by radarclinometry in arid or semi-arid regions[END_REF][START_REF] Baghdadi | A New Empirical Model for Radar Scattering from Bare Soil Surfaces[END_REF]. In this context, in order to ensure that the retrieved values of reflectivity are accurate, it is essential to evaluate their temporal stability. Figure 3 plots the signals retrieved in the four areas cited in Table 2, characterizing different types of surface (dunes, smooth surfaces, relief). A 2.5-year time series is shown for each zone, with the plotted values corresponding to the mean reflectivity within a 3-month period and a relatively small (0.1 • × 0.1 • ) zone. For some months of the year, no data were available since no specular points corresponding to such a small zone were detected. The reflectivities plotted in Figure 3 confirm the high temporal stability of the reflectivity data, with a negligible standard deviation (std dev), less than 0.7 dB for three of these sites.

of reflectivity are accurate, it is essential to evaluate their temporal stability. Figure 3 plots the signals retrieved in the four areas cited in Table 2, characterizing different types of surface (dunes, smooth surfaces, relief). A 2.5-year time series is shown for each zone, with the plotted values corresponding to the mean reflectivity within a 3-month period and a relatively small (0.1° × 0.1°) zone. For some months of the year, no data were available since no specular points corresponding to such a small zone were detected. The reflectivities plotted in Figure 3 confirm the high temporal stability of the reflectivity data, with a negligible standard deviation (std dev), less than 0.7 dB for three of these sites.

In the lower right graph (Figure 3d), corresponding to the fourth site, Tamanrasset-Mali, small seasonal effects are revealed by variations in reflectivity; in particular, towards the end of the summer of 2018. A clear correlation can be seen between the peak of the monsoon season in the Sahel region, which falls between August and September, and the peak in reflectivity at Tamanrasset in 2018 (-15.5dB for Q3). The reflectivity observed for July-September is approximately +2.26 dB greater than that observed in January. This is certainly due to the monsoon rainfall [START_REF] Cook | Contemporary Climate Change of the African Monsoon Systems[END_REF], which increases the dielectric constant of the soil and leads to higher values of reflectivity. In the lower right graph (Figure 3d), corresponding to the fourth site, Tamanrasset-Mali, small seasonal effects are revealed by variations in reflectivity; in particular, towards the end of the summer of 2018. A clear correlation can be seen between the peak of the monsoon season in the Sahel region, which falls between August and September, and the peak in reflectivity at Tamanrasset in 2018 (-15.5 dB for Q3). The reflectivity observed for July-September is approximately +2.26 dB greater than that observed in January. This is certainly due to the monsoon rainfall [START_REF] Cook | Contemporary Climate Change of the African Monsoon Systems[END_REF], which increases the dielectric constant of the soil and leads to higher values of reflectivity.

Reflectivity Mapping over Sahara

Following the normalization of the CYGNSS measurements and validation of the stability of the Sahara's reflectivity during the studied period, the data were arranged in a grid format, each pixel of which corresponds to the mean value over the full period of observation; i.e., nearly 2.5 years. The challenge of this process is to retain the finest possible grid resolution, with the constraint of using a sufficiently large number of observations to ensure the accuracy of the resulting average value.

For this purpose, three different resolutions were tested. Figure 4 shows the reflectivity histograms obtained with resolutions of 0.1 • , 0.05 • and 0.03 • (close to 3.5 km), and simple Gaussian distributions fitting the experimental data. It can be seen that this approximation is better adapted to the Sudan-Nile site (Figure 4b) than to the Kufrah site (Figure 4a). This could be explained by the presence of strong heterogeneities at the former site, which includes both flat and mountainous areas. Once the 0.03° resolution was selected, a minimum number of measurements (four CYGNSS data points) was used to achieve a reliable time average for each pixel. This threshold was determined empirically, using qualitative observations of retrieved maps, and our analysis of the sensitivity of the estimated reflectivity to the number of measurements.

Grid points that were not covered by four measurements were not included in the final grid. These points correspond mainly to high altitude zones (>650 m), for which the data were not analyzed (see above regarding data filtering). Figure 5 shows the resulting reflectivity map at a resolution of 0.03° × 0.03°, which covers the latitudes between 15°N and 37°30'N and the longitudes between 17°W and 39°30'E. Mountainous regions higher than 650 m-for example, the Hoggar Mountains in the South-East of Algeria, the Tibesti Mountains in the southern Libya, and the Atlas Mountains in the North-East of Morocco (Maghreb)-are not mapped. Regions with sand dunes, such as the Western and Eastern Erg in the North-West and North-East of Algeria respectively, the Great Libyan Desert in the East of Libya bordering with Egypt, and the Aoukar Erg in the South-East of Mauritania, are recognizable by their lower received signals due probably to a high wave penetration in soil. In Figure 5, mountains and areas with low received signals are masked in white. At a resolution of 0.1 • , the data are smoothed by relatively high spatial averaging, such that the dynamic range (very low or very high reflectivity) is very much reduced. This corresponds to a stronger degradation (smoothing) of fine geometric structures when compared to the results obtained at higher resolutions; i.e., 0.05 • and 0.03 • . At these finer resolutions, very similar histograms are observed. We finally chose the grid at a resolution of 0.03 • , thus allowing the finest possible geometrical structures to be resolved from the processed data. It is interesting to confirm that this corresponds to a CYGNSS spatial resolution of about 3.5 km. It is important to note that the purely coherent signal has a spatial resolution corresponding to the first Fresnel zone [START_REF] Camps | Spatial Resolution in GNSS-R Under Coherent Scattering[END_REF], whereas the spatial resolution (~25 km) of the purely incoherent signal, represented by an iso-delay ellipse projected on the ground, depends on the chip range resolution (~300 m) of the GPS L1 C/A signal. This empirical result shows that we can distinguish details, which are larger than the Fresnel zone, but also much smaller than the resolution of the incoherent signal.

Once the 0.03 • resolution was selected, a minimum number of measurements (four CYGNSS data points) was used to achieve a reliable time average for each pixel. This threshold was determined empirically, using qualitative observations of retrieved maps, and our analysis of the sensitivity of the estimated reflectivity to the number of measurements.

Grid points that were not covered by four measurements were not included in the final grid. These points correspond mainly to high altitude zones (>650 m), for which the data were not analyzed (see above regarding data filtering). Figure 5 shows the resulting reflectivity map at a resolution of 0.03 • × 0.03 • , which covers the latitudes between 15 • N and 37

Results and Discussions

In this section, we focus on the relationship between the CYGNSS reflectivity and the large-scale roughness observed over different types of the Sahara environment. We discuss the influence of geometric roughness and the aerodynamic roughness length Z 0 on the signals received by CYGNSS. As SAR data are often used to map roughness in deserts and other regions of the globe, it is worthwhile to carefully compare the CYGNSS reflectivity results with SAR roughness measurements derived from ALOS-2 backscattering coefficients.

Comparison between CYGNSS and ALOS-2 Data

Prior to our analysis of the influence of roughness on CYGNSS Γ data, the latter was compared with the ALOS-2 backscattering coefficient (or sigma nought: σ 0 ) data recorded over the regions described in Section 2.2. As shown in Figure 6, a linear relationship can be established between the ALOS-2 backscattering coefficient and the CYGNSS reflectivity. This comparison was carried out for the Kufrah area, where each point corresponds to a single 0.03 • × 0.03 • pixel for both types of data (Γ and σ 0 ). A trend towards decreasing reflectivity, as a function of the increasing backscattering coefficient, can be seen in this figure. As it is known that in desert areas, variations of SAR data depend mainly on soil roughness, this relationship between the CYGNSS Γ and SAR σ 0 data could potentially be used to retrieve surface roughness from CYGNSS observations. The relatively noisy correlation between CYGNSS Γ and SAR σ 0 data is attributed mainly to the fact that incoherent SAR radar signals (ALOS-2 data) are compared with reflectivities that are associated to the coherent components of CYGNSS data. This means that different surface roughness properties are contributing to each of these remotely sensed parameters.

Remote Sens. 2020, 12, x; doi: FOR PEER REVIEW www.mdpi.com/journal/remotesensing

Results and Discussions

In this section, we focus on the relationship between the CYGNSS reflectivity and the large-scale roughness observed over different types of the Sahara environment. We discuss the influence of geometric roughness and the aerodynamic roughness length Z0 on the signals received by CYGNSS. As SAR data are often used to map roughness in deserts and other regions of the globe, it is worthwhile to carefully compare the CYGNSS reflectivity results with SAR roughness measurements derived from ALOS-2 backscattering coefficients.

Comparison between CYGNSS and ALOS-2 Data

Prior to our analysis of the influence of roughness on CYGNSS Γ data, the latter was compared with the ALOS-2 backscattering coefficient (or sigma nought: σ0) data recorded over the regions described in Section 2.2. As shown in Figure 6, a linear relationship can be established between the ALOS-2 backscattering coefficient and the CYGNSS reflectivity. This comparison was carried out for the Kufrah area, where each point corresponds to a single 0.03° × 0.03° pixel for both types of data (Γ and σ0). A trend towards decreasing reflectivity, as a function of the increasing backscattering coefficient, can be seen in this figure. As it is known that in desert areas, variations of SAR data depend mainly on soil roughness, this relationship between the CYGNSS Γ and SAR σ0 data could potentially be used to retrieve surface roughness from CYGNSS observations. The relatively noisy correlation between CYGNSS Γ and SAR σ0 data is attributed mainly to the fact that incoherent SAR radar signals (ALOS-2 data) are compared with reflectivities that are associated to the coherent components of CYGNSS data. This means that different surface roughness properties are contributing to each of these remotely sensed parameters.

Relationship between Geometric Roughness and CYGNSS Γ/ALOS-2 Sigma Naught Data

In order to define large-scale roughness, mean slope values were considered in the same areas as those for which measurements of aerodynamic roughness of the terrain had been carried out (Table 1). Slope values were computed starting from the digital elevation model (DEM) of the SRTM mission, and the Horn algorithm was then used to compute the mean slope [START_REF] Horn | Hill shading and the reflectance map[END_REF] from consecutive

Relationship between Geometric Roughness and CYGNSS Γ/ALOS-2 Sigma Naught Data

In order to define large-scale roughness, mean slope values were considered in the same areas as those for which measurements of aerodynamic roughness of the terrain had been carried out (Table 1). Slope values were computed starting from the digital elevation model (DEM) of the SRTM mission, and the Horn algorithm was then used to compute the mean slope [START_REF] Horn | Hill shading and the reflectance map[END_REF] from consecutive blocks of eight pixels. The mean slope values were then gridded onto maps with 0.03 • × 0.03 • pixels in order to obtain the same resolution as that used for the CYGNSS and ALOS-2 data.

Figure 7 shows the mean values of CYGNSS Γ and ALOS-2 σ 0 data corresponding to the 0.03 • pixels as a function of the mean slopes. An exponentially decreasing trend can be seen for the CYGNSS Γ data as a function of increasing values of mean slope. Conversely, an empirical logarithmic increasing trend can be seen for the ALOS-2 backscattering coefficients as a function of increasing values of mean slope. At low values of slope, the corresponding values of CYGNSS Γ and ALOS-2 data are strongly dispersed. This is probably due to the presence of different physical, micro-topographical contexts (flat rocky surfaces, sandy soils, etc.), which can significantly affect the coherent component of the CYGNSS signals and backscattered ALOS-2 data [START_REF] Zribi | A new soil roughness parameter for the modelling of radar backscattering over bare soil[END_REF]. An increase in the strength of the incoherent component, as a function of roughness, could also lead to a certain degree of saturation of the CYGNSS reflectivity [START_REF] Al-Khaldi | Time-Series Retrieval of Soil Moisture Using CYGNSS[END_REF]. The correlation between roughness and the ALOS-2 data reveals a stronger degree of correlation at high levels of roughness, corresponding to areas with relief.

Remote Sens. 2020, 12, x FOR PEER REVIEW 2 of 18 blocks of eight pixels. The mean slope values were then gridded onto maps with 0.03° × 0.03° pixels in order to obtain the same resolution as that used for the CYGNSS and ALOS-2 data. Figure 7 shows the mean values of CYGNSS Γ and ALOS-2 σ0 data corresponding to the 0.03° pixels as a function of the mean slopes. An exponentially decreasing trend can be seen for the CYGNSS Γ data as a function of increasing values of mean slope. Conversely, an empirical logarithmic increasing trend can be seen for the ALOS-2 backscattering coefficients as a function of increasing values of mean slope. At low values of slope, the corresponding values of CYGNSS Γ and ALOS-2 data are strongly dispersed. This is probably due to the presence of different physical, microtopographical contexts (flat rocky surfaces, sandy soils, etc.), which can significantly affect the coherent component of the CYGNSS signals and backscattered ALOS-2 data [START_REF] Zribi | A new soil roughness parameter for the modelling of radar backscattering over bare soil[END_REF]. An increase in the strength of the incoherent component, as a function of roughness, could also lead to a certain degree of saturation of the CYGNSS reflectivity [START_REF] Al-Khaldi | Time-Series Retrieval of Soil Moisture Using CYGNSS[END_REF]. The correlation between roughness and the ALOS-2 data reveals a stronger degree of correlation at high levels of roughness, corresponding to areas with relief.

(a) (b)

Aerodynamic Roughness Length vs. CYGNSS/ALOS-2 Data

Following the above analysis of the influence of geometric roughness, expressed in terms of the mean slope, this section discusses the influence of aerodynamic roughness. The analysis is based on the ground measurements illustrated in Table1.

For each area with available ground measurements, the mean value of CYGNSS Γ is computed. Figure 8 shows the empirical relationship between reflectivity Γ and Z0 for the ground measurement areas considered in Table 1. An approximately linear trend is found, with a correlation R 2 equal to 0.81. The standard deviation of the reflectivity during the 2.5 years covered by this study is also shown, in order to appreciate the stability of the relationship.

Aerodynamic Roughness Length vs. CYGNSS/ALOS-2 Data

Following the above analysis of the influence of geometric roughness, expressed in terms of the mean slope, this section discusses the influence of aerodynamic roughness. The analysis is based on the ground measurements illustrated in Table 1.

For each area with available ground measurements, the mean value of CYGNSS Γ is computed. Figure 8 shows the empirical relationship between reflectivity Γ and Z 0 for the ground measurement areas considered in Table 1. An approximately linear trend is found, with a correlation R 2 equal to 0.81. The standard deviation of the reflectivity during the 2.5 years covered by this study is also shown, in order to appreciate the stability of the relationship.

The same procedure was applied to the ALOS-2 backscattering coefficients, taking into account the mean value of the backscattering coefficients inside the areas of interest. A linear trend is also found for this relationship (R 2 = 0.63), which is very similar to that retrieved from ERS measurements [START_REF] Prigent | Estimation of the aerodynamic roughness length in arid and semi-arid regions over the globe with the ERS scatterometer[END_REF]. In conclusion, the results of this study illustrate the relevance of CYGNSSS reflectivity measurements for the retrieval of aerodynamic roughness lengths in desert locations. The same procedure was applied to the ALOS-2 backscattering coefficients, taking into account the mean value of the backscattering coefficients inside the areas of interest. A linear trend is also found for this relationship (R 2 = 0.63), which is very similar to that retrieved from ERS measurements [START_REF] Prigent | Estimation of the aerodynamic roughness length in arid and semi-arid regions over the globe with the ERS scatterometer[END_REF]. In conclusion, the results of this study illustrate the relevance of CYGNSSS reflectivity measurements for the retrieval of aerodynamic roughness lengths in desert locations.

Aerodynamic Roughness Length Mapping

On the basis of the empirical relationship developed in Section 4.3 the reflectivity is inverted in order to generate a map of aerodynamic roughness over the entire Sahara. As the accuracy of this relationship is moderate, four aerodynamic roughness classes are considered, as shown in Figure 9: the first class indicates comparatively smooth areas (Z0 < 0.1), corresponding to flat rocky surfaces; the second class, with Z0 lying between 0.1 and 0.2, indicates medium rough surfaces; the third class, with Z0 lying between 0.2 and 0.3, corresponds to rough rocks and hills; and finally, the highest roughness class, with Z0 > 0.3, corresponds to non-masked mountainous areas. The map shown in Figure 9a was produced using an inversion of the empirical relationship in Figure 8a over a 10 × 10 pixel window to reduce local heterogeneities and retrieve the main structural features. It does not include the areas in the CYGNSS data that were filtered out, due either to their high altitudes or to the presence of dunes characterized by scattered signals with low SNR that were probably highly sensitive to subsurface features. The elimination of dunes was also considered in studies using data from the ERS scatterometer (Figure 9b). Despite the differences in the definition of roughness ranges for each class, maps derived from CYGNSS and from the ERS scatterometer [START_REF] Prigent | Estimation of the aerodynamic roughness length in arid and semi-arid regions over the globe with the ERS scatterometer[END_REF] present significant similarities. With the ERS scatterometer, classes are particularly concentrated in very low Z0 values.

Aerodynamic Roughness Length Mapping

On the basis of the empirical relationship developed in Section 4.3 the reflectivity is inverted in order to generate a map of aerodynamic roughness over the entire Sahara. As the accuracy of this relationship is moderate, four aerodynamic roughness classes are considered, as shown in Figure 9: the first class indicates comparatively smooth areas (Z 0 < 0.1), corresponding to flat rocky surfaces; the second class, with Z 0 lying between 0.1 and 0.2, indicates medium rough surfaces; the third class, with Z 0 lying between 0.2 and 0.3, corresponds to rough rocks and hills; and finally, the highest roughness class, with Z 0 > 0.3, corresponds to non-masked mountainous areas. The map shown in Figure 9a was produced using an inversion of the empirical relationship in Figure 8a over a 10 × 10 pixel window to reduce local heterogeneities and retrieve the main structural features. It does not include the areas in the CYGNSS data that were filtered out, due either to their high altitudes or to the presence of dunes characterized by scattered signals with low SNR that were probably highly sensitive to subsurface features. The elimination of dunes was also considered in studies using data from the ERS scatterometer (Figure 9b). Despite the differences in the definition of roughness ranges for each class, maps derived from CYGNSS and from the ERS scatterometer [START_REF] Prigent | Estimation of the aerodynamic roughness length in arid and semi-arid regions over the globe with the ERS scatterometer[END_REF]

Conclusions

GNSS-R data derived from CYGNSS observations have been used for the first time to characterize surface properties in a major desert area. This study demonstrates the strong stability of GNSS-R reflectivities (Γ) in the Sahara desert, meaning that the long-term average of the reflectivity can be used to characterize surface properties. Thanks to this high stability, a grid could be produced to map the average reflectivity, with a fine resolution of 0.03° × 0.03°. The main structures of the dunes and reliefs are thus qualitatively observable. Our comparison of CYGNSS Γ data and SAR ALOS-2/PALSAR σ0 measurements reveals a robust correlation. The SAR radar signal is found to decrease with increasing reflectivity. This is directly related to roughness, with the main parameter affecting radar signal strength, in the absence of variations in soil water content. Using ground measurements of aerodynamic roughness or geometric roughness computed from the SRTM product, a relatively strong correlation is observed between the surface parameters and the reflectivity, with the correlation R 2 equal to 0.81 for the aerodynamic roughness. From the resulting empirical relationship, a map of the Sahara showing four classes of aerodynamic roughness is proposed. These classes range from very smooth surfaces to mountainous surfaces corresponding to the highest values of roughness. This map strongly resembles the roughness products generated by scatterometer data. In the future, it would be of significant importance to gain a deeper understanding of the signals observed on the dunes, which remain difficult to interpret. This is probably due to the influence of sub-surface features and is not yet understood. It would also be important to validate the experimental observations with physical models, in order to reproduce the DDMs observed on each type of surface.

Conclusions

GNSS-R data derived from CYGNSS observations have been used for the first time to characterize surface properties in a major desert area. This study demonstrates the strong stability of GNSS-R reflectivities (Γ) in the Sahara desert, meaning that the long-term average of the reflectivity can be used to characterize surface properties. Thanks to this high stability, a grid could be produced to map the average reflectivity, with a fine resolution of 0.03 • × 0.03 • . The main structures of the dunes and reliefs are thus qualitatively observable. Our comparison of CYGNSS Γ data and SAR ALOS-2/PALSAR σ 0 measurements reveals a robust correlation. The SAR radar signal is found to decrease with increasing reflectivity. This is directly related to roughness, with the main parameter affecting radar signal strength, in the absence of variations in soil water content. Using ground measurements of aerodynamic roughness or geometric roughness computed from the SRTM product, a relatively strong correlation is observed between the surface parameters and the reflectivity, with the correlation R 2 equal to 0.81 for the aerodynamic roughness. From the resulting empirical relationship, a map of the Sahara showing four classes of aerodynamic roughness is proposed. These classes range from very smooth surfaces to mountainous surfaces corresponding to the highest values of roughness. This map strongly resembles the roughness products generated by scatterometer data. In the future, it would be of significant importance to gain a deeper understanding of the signals observed on the dunes, which remain difficult to interpret. This is probably due to the influence of sub-surface features and is not yet understood. It would also be important to validate the experimental observations with physical models, in order to reproduce the DDMs observed on each type of surface.

Figure 1 .

 1 Figure 1. Illustration of the studied site: the Sahara desert.

Figure 1 .

 1 Figure 1. Illustration of the studied site: the Sahara desert.

Figure 2 .

 2 Figure 2. Plot of variations in Cyclone Global Navigation Satellite System (CYGNSS) reflectivity as a function of incidence angles: (a) Algeria dunes, (b) Kufrah area.

Figure 2 .

 2 Figure 2. Plot of variations in Cyclone Global Navigation Satellite System (CYGNSS) reflectivity as a function of incidence angles: (a) Algeria dunes, (b) Kufrah area.

Figure 3 .

 3 Figure 3. Temporal variations of CYGNSS reflectivity over a period of 2.5 years: (a) Algerian dunes, (b) Egypt dunes, (c) Kufrah (Libya), (d) Tamanrasset-Mali.

Figure 4 .

 4 Figure 4. Reflectivity histograms for three spatial resolutions: 0.1°, 0.05° and 0.03°, at two different locations: (a) Kufrah area (Libya), (b) Sudan-Nile area.

Figure 4 .

 4 Figure 4. Reflectivity histograms for three spatial resolutions: 0.1 • , 0.05 • and 0.03 • , at two different locations: (a) Kufrah area (Libya), (b) Sudan-Nile area.

18 Figure 5 .

 185 Figure 5. Average reflectivity (dB) map (mosaic) of the Sahara region, with a gridded-pixel size of 0.03° × 0.03°. Mountains and areas with low received signals are masked (white color).

Figure 5 .

 5 Figure 5. Average reflectivity (dB) map (mosaic) of the Sahara region, with a gridded-pixel size of 0.03 • × 0.03 • . Mountains and areas with low received signals are masked (white color).

Figure 6 .

 6 Figure 6. Plot of Advanced Land Observing Satellite-2 (ALOS-2) backscattering coefficient vs. CYGNSS reflectivity data collected over the Kufrah site, where each point corresponds to a pixel with a spatial resolution of 0.03° × 0.03°.

Figure 6 .

 6 Figure 6. Plot of Advanced Land Observing Satellite-2 (ALOS-2) backscattering coefficient vs. CYGNSS reflectivity data collected over the Kufrah site, where each point corresponds to a pixel with a spatial resolution of 0.03 • × 0.03 • .

Figure 7 .

 7 Figure 7. Plot of roughness mean slope and satellite data: (a) CYGNSS reflectivity, (b) ALOS-2 backscattering coefficient.

Figure 7 .

 7 Figure 7. Plot of roughness mean slope and satellite data: (a) CYGNSS reflectivity, (b) ALOS-2 backscattering coefficient.

Figure 8 .

 8 Figure 8. Plot of the aerodynamic roughness length as a function of satellite data: (a) CYGNSS reflectivity, (b) ALOS-2 backscattering coefficient.

Figure 8 .

 8 Figure 8. Plot of the aerodynamic roughness length as a function of satellite data: (a) CYGNSS reflectivity, (b) ALOS-2 backscattering coefficient.

Figure 9 .

 9 Figure 9. Cont.

Figure 9 .

 9 Figure 9. Maps of aerodynamic roughness classes over the Sahara desert: (a) derived from a CYGNSS reflectivity database, (b) derived from the European Remote-Sensing Satellite (ERS) scatterometer (Prigent et al. [4]).

Figure 9 .

 9 Figure 9. Maps of aerodynamic roughness classes over the Sahara desert: (a) derived from a CYGNSS reflectivity database, (b) derived from the European Remote-Sensing Satellite (ERS) scatterometer (Prigent et al. [4]).

Table 1 .

 1 Aerodynamic roughness ground measurements over the Sahara site[START_REF] Marticorena | Modeling the atmospheric dust cycle: 2. Simulation of Saharan dust sources[END_REF].

	Number	Regions	Type of Terrain	Coordinates	Z 0 [cm]
	1	Adrar Desert-Mauritania	Dunes	20 • 44 N 9 • 44 W	20 • 18 N 9 • 18 W	0.0023
	2	Great Eastern Erg-Algeria	Dunes	31 • 36 N 7 • 00 E	31 • 14 N 7 • 25 E	0.0023
	3	Adrar Desert-Mauritania	Dunes	21 • 48 N 7 • 42 W	21 • 20 N 7 • 15 W	0.0023
	4	Great Eastern Erg-Algeria	Dunes	31 • 39 N 8 • 34 E	31 • 13 N 9 • 03 E	0.0023
	5	Bordj Badji Mokhtar-Algeria	Sandy flat surface	23 • 42 N 0 • 59 E	23 • 18 N 1 • 22 E	0.01
	6	Tiris Zemmour Plateau-Mauritania	Rocky surface	26 • 12 N 7 • 48 W	25 • 48 N 7 • 18 W	0.025
	7	Tamanrasset-Algeria	Rocky surface	21 • 42 N 5 • 20 E	21 • 15 N 5 • 48 E	0.05
	8	Tamanrasset-Algeria	Rocky surface	21 • 40 N 4 • 20 E	21 • 21 N 4 • 48 E	0.05
	9	Nafusah Plateau, Tripolitania-Libya	Rocky surface	30 • 42 N 11 • 09 E	30 • 15 N 11 • 12 E	0.15
	10	Ubari Desert-Libya	Rocky surface (hammada)	30 • 10 N 10 • 45 E	29 • 45 N 11 • 12 E	0.15
	11	Tassili N'Ajjer Mountains-Algeria	Mountains	25 • 42 N 8 • 10 E	25 • 18 N 8 • 32 E	0.5
	12	Atlas Mountains-Algeria	Mountains	33 • 36 N 1 • 18 W	33 • 12 N 0 • 48 W	0.873
	13	Tanezrouft-Algeria	Rocky surface	26 • 42 N 4 • 54 W	26 • 18 N 4 • 27 W	0.05
	14	Tassili N'Ajjer Mountains-Algeria	Mountains	26 • 12 N 8 • 17 E	25 • 46 N 8 • 32 E	0.5
	15	Atlas Mountains-Algeria	Mountains	33 • 28 N 2 • 53 E	33 • 12 N 3 • 30 E	0.347
	16	Ouled Nail Mountains-Algeria	Mountains	33 • 37 N 3 • 31 E	33 • 12 N 4 • 01 E	0.347
	17	El Eglab Massif-Algeria	Mountains	26 • 59 N 3 • 44 W	26 • 46 N 3 • 09 W	0.131
	18	Hammada du Draa-Morocco	Rocky surface (hammada)	27 • 47 N 8 • 43 W	27 • 28 N 8 • 17 W	0.131
	19	Tidikelt Plain-Algeria	Sandy flat surface	26 • 41 N 2 • 00 E	27 • 09 N 2 • 24 E	0.0873
	20	Tademait Plateau-Algeria	Sandy flat surface	28 • 41 N 2 • 40 E	28 • 09 N 2 • 55 E	0.0873
	21	Hassi Gara, Great Western Erg-Algeria	Sandy flat surface	29 • 46 N 2 • 58 E	29 • 14 N 3 • 35 E	0.0873
	22	Tiris Zemmour Plateau-Mauritania	Sandy flat surface	23 • 50 N 7 • 52 W	23 • 27 N 7 • 37 W	0.0174
	23	Tanerzrouft-n-Ahnet-Algeria	Sandy flat surface	22 • 30 N 0 • 53 E	22 • 17 N 0 • 37 E	0.01

Table 2 .

 2 Coordinates of test areas used for CYGNSS data processing.

Remote Sens. 2020, 12, 743; doi:10.3390/rs12040743 www.mdpi.com/journal/remotesensing

Remote Sens. 2020, 12, x FOR PEER REVIEW

Funding: This research was funded by CNES/TOSCA program.

Conflicts of Interest:

The authors declare no conflict of interest.

Author Contributions: M.Z. supervised and designed the methodologies; D.S. developed the algorithms; M.H. participated in the data processing, N.P. provided advice on CYGNSS data processing and methodologies; N.B. provided advice on results interpretation; D.S. and M.Z. wrote the paper. All authors have read and agreed to the published version of the manuscript.