Time-series retrievals of soil moisture obtained from the Cyclone Global Navigation Satellite System (CYGNSS) constellation are presented. The retrieval approach assumes that vegetation and roughness changes occur on timescales longer than those associated with soil moisture changes to allow soil moisture sensing in the presence of vegetation and surface roughness contributions as well as the varying incidence angles associated with spaceborne Global Navigation Satellite System-Reflectometry (GNSS-R) systems. The approach is focused on incoherent scattering from land surfaces due to the expectation that coherent land surface returns arise primarily from inland water body contributions that are not directly representative of soil moisture. An approach for discarding coherent CYGNSS measurements is therefore developed and described. Because the approach requires the retrieval of N temporal soil moisture samples at a given location but uses only N-1 ratios of CYGNSS measured quantities, ancillary information is incorporated in the retrieval through the use of maximum and minimum monthly soil moisture maps obtained from the Soil Moisture Active Passive (SMAP) mission. Retrieved soil moistures are presented for the 6-month period December 2017-May 2018 and are compared against values reported by the SMAP mission. The comparisons suggest that there exists the potential for using spaceborne GNSS-R systems for global soil moisture retrievals with an rms error on the order of 0.04 cm ³ /cm ³ over varied terrain.

An algorithm for detecting coherence in Cyclone Global Navigation Satellite System (CYGNSS) mission delay-Doppler maps (DDMs) is presented. Because CYGNSS DDMs report only the observed power without phase information, the algorithm uses estimates of power ``spread'' within the DDM to flag coherency. Since the estimate used is a ratio of the powers in differing portions of the DDM, it is less sensitive to absolute power calibration and to the GPS C/A code type observed, and is applied to CYGNSS Level-1 uncalibrated DDMs. The basic detector formulation is described along with modifications to improve performance in lower signal-to-noise ratio (SNR) situations. The required detection thresholds are determined using matchups with CYGNSS ``Raw I/F'' mode measurements for which the DDM phase can be computed and used to identify coherence more precisely. Application of the final detector over a large CYGNSS data set suggests that approximately 8.9% of all inland returns are coherent. Inland regions persistently identified as coherent were found largely to be associated with the presence of water bodies. A smaller set of desert locations apparently having very low surface roughness were also found to be associated with persistent coherence. The detector was also applied to a set of ocean measurements, with the results showing that persistent coherence is limited to areas with sheltered waters. Ocean tests avoiding such regions indicate that the detector's false-alarm rate is approximately 0.0012% for the detection threshold used.

This work demonstrates the creation of dynamic inland water body masks at spatial resolutions ranging from 1 to 3 km through the use of a recently developed coherence detector for the delay-Doppler maps produced by the cyclone global navigation satellite system (CYGNSS) constellation. The use of the coherence of the observed measurements reduces many of the uncertainties associated with previous signal-to-noise ratio-based water body detection approaches for CYGNSS. Using data from January 2018 to February 2020 and producing maps representing time intervals ranging from 3 months to 2 years, the water body masks created are found to be associated with a probability of detection that exceeds 80% as compared to the Pekel water mask developed from Landsat observations. The analysis presented in this work highlights the potential of using spaceborne Global Navigation Satellite Systems Reflectometry (GNSS-R) systems for dynamic inland water body mapping.

An analysis of the coherency properties of specular scattering from ocean and land surfaces as observed in global navigation satellite system-reflectometry (GNSS-R) and signals of opportunity systems is presented. This analysis applies existing approximate models for the coherent and incoherent contributions. Approximate expressions are developed for when one component of the return is likely to dominate as a function of surface and observing system properties. The model developed is then applied for sea surface returns, and the relative contribution of the coherent term is expressed as a function of the receiver height, frequency, incidence angle, and wind speed. For L-band spaceborne measurements, it is shown that coherence is expected only for wind speeds less than 2-3 m/s, while for P-band spaceborne measurements, coherence can dominate returns for wind speeds up to 5-7 m/s. For land surface measurements from space, it is shown that the surface rms height needs to be sufficiently low for coherent components to dominate returns. Coherence dominates for roughness values not exceeding a range of 5-7 cm for the L-band and 15-30 cm for the P-band. For the L-band, these conditions over land are likely to be created primarily by inland water bodies. A model for the specular scattering from a water body, including earth curvature effects, is then developed to highlight the strong dependence of the resulting coherent field on the shape of the water body and any offset in its location from the specular point. These results further clarify the significant variability that should be expected in coherent scattering from inland water bodies.

This paper presents the processing algorithms for geolocating and calibration of the Cyclone Global Navigation Satellite System (CYGNSS) level 1 land data products, as well as analysis of the spatial resolution of Global Navigation Satellite System Reflectometry (GNSS-R) coherent reflections. Accurate and robust geolocation and calibration of GNSS-R land observations are necessary first steps that enable subsequent geophysical parameter retrievals. The geolocation algorithm starts with an initial specular point location on the Earth’s surface, predicted by modeling the Earth as a smooth ellipsoid (the WGS84 representation) and using the known transmitting and receiving satellite locations. Information on terrain topography is then compiled from the Shuttle Radar Topography Mission (SRTM) generated Digital Elevation Map (DEM) to generate a grid of local surface points surrounding the initial specular point location. The delay and Doppler values for each point in the local grid are computed with respect to the empirically observed location of the Delay Doppler Map (DDM) signal peak. This is combined with local incident and reflection angles across the surface using SRTM estimated terrain heights. The final geolocation confidence is estimated by assessing the agreement of the three geolocation criteria at the estimated surface specular point on the local grid, including: the delay and Doppler values are in agreement with the CYGNSS observed signal peak and the incident and reflection angles are suitable for specular reflection. The resulting geolocation algorithm is first demonstrated using an example GNSS-R reflection track that passes over a variety of terrain conditions. It is then analyzed using a larger set of CYGNSS data to obtain an assessment of geolocation confidence over a wide range of land surface conditions. Following, an algorithm for calibrating land reflected signals is presented that considers the possibility of both coherent and incoherent scattering from land surfaces. Methods for computing both the bistatic radar cross section (BRCS, for incoherent returns) and the surface reflectivity (for coherent returns) are presented. a flag for classifying returns as coherent or incoherent developed in a related paper is recommended for use in selecting whether the BRCS or reflectivity should be used in further analyses for a specific DDM. Finally, a study of the achievable surface feature detection resolution when coherent reflections occur is performed by examining a series of CYGNSS coherent reflections across an example river. Ancillary information on river widths is compared to the observed CYGNSS coherent observations to evaluate the achievable surface feature detection resolution as a function of the DDM non-coherent integration interval.

This article presents a detailed examination of fluctuating input noise power levels on the analog-to-digital convertor (ADC) sampling hardware of the NASA Cyclone Global Navigation Satellite System (CYGNSS) instruments and the associated impacts on the level 1 normalized bistatic radar cross section (NBRCS) estimation performance. The impact of external noise variations on both the CYGNSS science and navigation channels is quantified with respect to how the fluctuating received power impacts the low-level ADC sampling distribution and subsequent NBRCS estimation. This work demonstrates that there are clear quantifiable geospatially dependent noise variations linked to Global Navigation Satellite System (GNSS) space-based augmentation systems [notably Japanese Quasi-Zenith Satellite System (QZSS) and U.S. Wide Area Augmentation System (WAAS)] and that these additional noise sources significantly alter the digital sampling distribution of the CYGNSS instruments, actively degrading the level 1 NBRCS estimation if not corrected. A derivation of the theoretical correction for both the science and navigation channel ADC sampling variations is presented which is later tuned based on empirical performance metrics in an effort to minimize the induced calibration errors. The impacts of the enhancements outlined in this work on CYGNSS level 1 calibration are evaluated using one year of observations before and after the digital sampling corrections, using model European Centre for Medium-Range Weather Forecasts (ECMWF) wind and Wavewatch III mean square slope (MSS) surface validation datasets.