This article shows the first spectral analysis of fully‐focused Synthetic Aperture Radar (FFSAR) altimetry data with the objective of studying backscatter modulations caused by swells. Swell waves distort the backscatter in altimetry radargrams by means of velocity and range bunching. These swell signatures are visible in the tail of the waveform. By locally normalizing the backscatter and projecting the waveforms on an along‐/cross‐track grid, satellite altimetry can be exploited to retrieve swell information. The analysis of FFSAR spectra is supported by buoy‐derived swell‐wave spectra of the National Oceanic and Atmospheric Administration network. Using cases with varying wave characteristics, we discuss the altimetry‐derived spectra and relate them to what is known from side‐looking SAR imaging systems. Besides having a vast amount of additional data for swell‐wave analysis, altimeter data can also help us to better understand the side‐looking SAR spectra.
13 November 2024 in Journal of Geophysical Research: Oceans
This study investigates the influence of ocean currents on wave modeling and satellite observations using in situ wave measurements from the One Ocean Expedition 2021–2023. In January 2023, six OpenMetBuoy drifters were deployed in the Agulhas Current region. Their high immersion ratio minimized wind effects, allowing them to follow the current and return to the Indian Ocean by the Agulhas retroflection, collecting data for about 2 months. Comparing surface current velocities from both the Mercator model and Globcurrent product with drifter data reveals underestimation for velocities over with Mercator showing greater variability. Significant wave height and Stokes drift parameters from MFWAM and ERA5 were also evaluated against drifters. Both models tend to overestimate Stokes drift more noticeable in ERA5, indicating sensitivity to wind seas. For significant wave height, both models agree well with drifter measurements with correlations of 0.90 for MFWAM and 0.83 for ERA5. However, ERA5's lack of surface current data combined with its coarse resolution (0.5) lead to underestimation of wave heights exceeding 2.5 m. MFWAM products including and excluding currents exhibit root mean square errors of 0.39 and 0.45 m, respectively, when compared to drifter measurements. This confirms that neglecting currents introduces additional errors particularly in areas with sharp current gradients. Analyzing MFWAM wave spectra, including and excluding currents, reveals wave energy transfer attributed to wave‐current interactions. The spatial extent of these interactions is captured by satellite altimeters, revealing wave modulations with considerable wave height variations when waves cross eddies and the current core.
This study presents the first azimuth cutoff analysis in Synthetic Aperture Radar (SAR) altimetry, aiming to assess its applicability in characterizing sea-state dynamics. In SAR imaging, the azimuth cutoff serves as a proxy for the shortest waves, in terms of wavelength, that can be detected by the satellite under certain wind and wave conditions. The magnitude of this parameter is closely related to the wave orbital velocity variance, a key parameter for characterizing wind-wave systems. We exploit wave modulations exhibited in the tail of fully-focused SAR waveforms and extract the azimuth cutoff from the radar signal through the analysis of its along-track autocorrelation function. We showcase the capability of Sentinel-6A in deriving these two parameters based on analyses in the spatial and wavenumber domains, accompanied by a discussion of the limitations. We use Level-1A high-resolution Sentinel-6A data from one repeat cycle (10 days) globally to verify our findings against wave modeled data. In the spatial domain analysis, the estimation of azimuth cutoff involves fitting a Gaussian function to the along-track autocorrelation function. Results reveal pronounced dependencies on wind speed and significant wave height, factors primarily determining the magnitude of the velocity variance. In extreme sea states, the parameters are underestimated by the altimeter, while in relatively calm sea states and in the presence of swells, a substantial overestimation trend is observed. We introduce an alternative approach to extract the azimuth cutoff by identifying the fall-off wavenumber in the wavenumber domain. Results indicate effective mitigation of swell-induced errors, with some additional sensitivity to extreme sea states compared to the spatial domain approach.
17 May 2024 in IEEE Transactions on Geoscience and Remote Sensing
Until recently, intensity modulations in synthetic aperture radar (SAR) altimetry waveform tails have been considered a nuisance for geophysical-parameter retrieval. These modulations are actually predictable and might be exploited using a spectral analysis of the waveform tails. After Altiparmaki et al. (2022), a more elaborated analysis is performed to improve the interpretation of these SAR altimeter spectra. A fast numerical model is developed to explain the modulation mechanisms in focused SAR altimetry waveform tails. Using numerical solutions, standard analytical closed-form solutions, are demonstrated to be invalid to retrieve ocean-wave-spectra retrievals from nadir altimeters. Although not valid, a closed-form derivation provides intuitive insights about the information contained in a SAR altimetry cross-spectrum. Under moderate environmental conditions (significant wave heights of ~2 m), a closed-form solution might still be useful to infer swell-wave spectra from swath-altimetry SAR spectra at incident angles of ~4°. Comparable to side-looking SAR ocean processing, the cross-spectral analysis for nadir signals reduces noise and might remove the 180-degree ambiguity of the wave direction. Since the synthetic aperture length of nadir altimeters is larger than sidelooking imaging SARs (e.g. Sentinel-1, RadarSat, Gaofen-3), sublook processing can be performed to compute multiple cross-spectra for the same scene. With a slightly changing observation geometry, the cross-spectra reveal slightly different parts of the ocean-wave spectrum. The resulting stack of cross-spectra can thus be used to improve the retrieval of ocean-wave parameters. Retrieved ocean-wave parameters shall then enhance the sampling of the global wave field, but also serve to advance more consistent sea-state-bias corrections.
In this study, the full-focusing (FF) algorithm is reviewed with the objective of optimizing it for processing data from different types of surfaces probed in altimetry. In particular, this work aims to provide a set of optimal FF processing parameters for the Sentinel-6 Michael Freilich (S6-MF) mission. The S6-MF satellite carries an advanced radar altimeter offering a wide range of potential FF-based applications which are just beginning to be explored and require prior optimization of this processing. In S6-MF, the Synthetic Aperture Radar (SAR) altimeter acquisitions are known to be aliased in the along-track direction. Depending on the target, aliasing can be tolerated or may be a severe impairment to provide the level of performance expected from FF processing. Another key aspect to consider in this optimization study is the unprecedented resolution of the FF processing, which results in a higher posting rate than the standard SAR processing. This work investigates the relationship between posting rate and noise levels and provides recommendations for optimal algorithm configurations in various scenarios, including transponder, open ocean, and specular targets like sea-ice and inland water scenes. The Omega–Kappa (WK) algorithm, which has demonstrated superior CPU efficiency compared to the back-projection (BP) algorithm, is considered for this study. But, unlike BP, it operates in the Doppler frequency domain, necessitating further precise spectral and time domain settings. Based on the results of this work, real case studies using S6-MF acquisitions are presented. We first compare S6-MF FF radargrams with Sentinel-1 (S1) images to showcase the potential of optimally configured FF processing. For highly specular surfaces such as sea-ice, distinct techniques are employed for lead signature identification. S1 relies on image-based lineic reconstruction, while S6-MF utilizes phase coherency of focalized pulses for lead detection. The study also delves into two-dimensional wave spectra derived from the amplitude modulation of image/radargrams, with a focus on a coastal example. This case is especially intriguing, as it vividly illustrates different sea states characterized by varying spectral peak positions over time.