Fast observations of rapid surface large-changes are demanded in disaster evaluations and scientific studies. Digital elevation model (DEM) differencing before and after the events is an effective way to retrieve the changes. Owing to a short repeat cycle, geosynchronous synthetic aperture radar (GEO SAR) systems can quickly obtain repeat-pass data and generate postevent DEMs by interferometry. However, interferometric baselines under its quick revisit cases are short, resulting in generating low-accuracy postevent DEMs. Moreover, surface large-changes can bring height ambiguity problems under the single-baseline interferometric processing. In this letter, we address the problem through a multibaseline (MB) processing. Since GEO SAR MB data can derive from the repeat-pass interferometric data of different subapertures and revisits, a subaperture-decomposition-based temporal and spatial MB method is proposed. The simulation results verify the effectiveness of the proposed method, where the quickly generated postevent DEM can help to realize the rapid large-elevation change observations.

Tropospheric delays are one of the main contributors to the interferometric phase in synthetic aperture radar (SAR) interferometry. When the phase contributions from surface deformation, topography, and ionospheric delays are negligible or known, the interferogram can be used to estimate the differential tropospheric delay (DTD), which can help to improve tropospheric delay predictions from weather models and in situ measurements. In conventional repeat-pass interferometric SAR (InSAR), however, the estimation of the DTD can still be significantly hindered by baseline errors. In addition, a single interferogram provides only relative DTDs, as the delays can be retrieved up to an unknown offset. To address such issues, this article presents a method for the estimation of DTDs on large scales by using repeat-pass simultaneous multi-angle SAR systems. Complementary simultaneous observations of the correlated troposphere from multiple angles are used to retrieve estimates of the absolute DTD and, at the same time, to mitigate the effect of baseline knowledge errors. Finally, a performance evaluation is presented for the Harmony Earth Explorer 10 candidate mission. A centimeter-level absolute accuracy and a submillimeter-level relative accuracy of the DTD estimation are achieved under the multistatic Harmony case when at least one companion satellite has an inter-satellite distance longer than 300 km to provide enough sensitivity.

This article describes the observation techniques and suggests processing methods to estimate dynamical sea-ice parameters from data of the Earth Explorer 10 candidate Harmony. The two Harmony satellites will fly in a reconfigurable formation with Sentinel-1D. Both will be equipped with a multi-angle thermal infrared sensor and a passive radar receiver, which receives the reflected Sentinel-1D signals using two antennas. During the lifetime of the mission, two different formations will be flown. In the stereo formation, the Harmony satellites will fly approximately 300km in front and behind Sentinel-1, which allows for the estimation of instantaneous sea-ice drift vectors. We demonstrate that the addition of instantaneous sea-ice drift estimates on top of the daily integrated values from feature tracking have benefits in terms of interpretation, sampling and resolution. The wide-swath instantaneous drift observations of Harmony also help to put high-temporal-resolution instantaneous buoy observations into a spatial context. Additionally, it allows for the extraction of deformation parameters, such as shear and divergence. As a result, Harmony's data will help to improve sea-ice statistics and parametrizations to constrain sea-ice models. In the cross-track interferometry (XTI) mode, Harmony's satellites will fly in close formation with an XTI baseline to be able to estimate surface elevations. This will allow for improved estimates of sea-ice volume and also enables the retrieval of full, two-dimensional swell-wave spectra in sea-ice-covered regions without any gaps. In stereo formation, the line-of-sight diversity allows the inference of swell properties in both directions using traditional velocity bunching approaches. In XTI mode, Harmony's phase differences are only sensitive to the ground-range direction swell. To fully recover two-dimensional swell-wave spectra, a synergy between XTI height spectra and intensity spectra is required. If selected, the Harmony mission will be launched in 2028.

Based on its ability to obtain two-dimensional (2D) high-resolution images in all-time and all-weather conditions, spaceborne synthetic aperture radar (SAR) has become an important remote sensing technique and the study of such systems has entered a period of vigorous development. Advanced imaging modes such as radar interferometry, tomography, and multi-static imaging, have been demonstrated. However, current in-orbit spaceborne SARs, which all operate in low Earth orbits, have relatively long revisit times ranging from several days to dozens of days, restricting their temporal sampling rate. Geosynchronous SAR (GEO SAR) is an active research area because it provides significant new capability, especially its much-improved temporal sampling. This paper reviews the research progress of GEO SAR technologies in detail. Two typical orbit schemes are presented, followed by the corresponding key issues, including system design, echo focusing, main disturbance factors, repeat-track interferometry, etc, inherent to these schemes. Both analysis and solution research of the above key issues are described. GEO SAR concepts involving multiple platforms are described, including the GEO SAR constellation, GEO-LEO/airborne/unmanned aerial vehicle bistatic SAR, and formation flying GEO SAR (FF-GEO SAR). Due to the high potential of FF-GEO SAR for three-dimensional (3D) deformation retrieval and coherence-based SAR tomography (TomoSAR), we have recently carried out some research related to FF-GEO SAR. This research, which is also discussed in this paper, includes developing a formation design method and an improved TomoSAR processing algorithm. It is found that GEO SAR will continue to be an active topic in the aspect of data processing and multi-platform concept in the near future.

Multiple CubeSat altimeters can work independently or corporately to form altimeter constellations. Different configurations of the constellations can acquire distinguished advantages: improved spatial/temporal sampling and high cross-track resolution, which will be helpful for observations of oceanic small-scale structures and weather forecasting. Compared to single conventional altimeters, CubeSat altimeter constellations may achieve better performances with lower costs. To fully understand these systems, this article focuses on the performance analysis and methodology for CubeSat altimeter constellations. Besides the typical analyses of the resolution, revisit, and absolute sea surface height (SSH) accuracy, the performance analysis was conducted by considering the characteristics of multiple measurements provided by CubeSat altimeter constellations. Local and global spatial sampling performances are investigated for various constellations and compared by sampling density and swath size. Moreover, relative SSH accuracy is introduced and evaluated based on the spatial structure functions of errors to effectively evaluate the measurement performance. Related system requirements on power, delta-v, etc., to achieve the performance are also discussed, which ensures that the analysis fits the boundary conditions of implementation. Finally, different concepts of the CubeSat altimeter constellations are compared, where their limitations and possible solutions are also discussed.

Coherence-based synthetic aperture radar (SAR) tomography (TomoSAR) exploits the complex coherences of SAR images to achieve 3-D imaging. Utilizing two-sensor spaceborne SAR formation flying to realize coherence-based TomoSAR has attracted increasing attention because temporal decorrelation-free interferograms can be constructed; therefore, TomoSAR has excellent potential for inverting the vertical structures of natural scenes such as forests and glaciers. However, low earth orbit (LEO) TomoSAR is disadvantaged by limited data and nonuniform sampling in the elevation direction. Geosynchronous (GEO) TomoSAR can overcome these limitations owing to its short revisit time of no more than 24 h. For the first time, this article discusses coherence-based TomoSAR exploiting GEO SAR formation flying. The benefits of GEO-formation coherence-based TomoSAR, including the low cost of slave satellites, rich data sets, and uniform sampling, are noted. The key problems of system design, including the formation design and data acquisition, are discussed. A formation design method based on the minimum along-track baseline is proposed that can realize uniform elevation sampling. The geometric correlation of a general SAR observation geometry is derived; on this basis, an optimal data acquisition method based on the optimal height measurement Cramer-Rao lower bound (CRLB) is proposed. Finally, the performance of GEO-formation coherence-based TomoSAR is analyzed; in particular, the ambiguity height in the altitude direction, the Rayleigh resolution in the altitude direction, and the theoretical optimal geometric correlation are evaluated. Finally, computer simulations validate the proposed formation design method, data acquisition scheme, and performance analysis formula.

Geosynchronous spaceborne-airborne bistatic synthetic aperture radar (GEO SA-BSAR), consisting of GEO transmitter and airborne receiver, has stable coverage for a long time and benefits moving target detection. However, the performance of GEO SA-BSAR moving target indication (MTI) system varies widely between bistatic configurations. The traditional configuration design for GEO SA-BSAR system only considers the imaging performance, which may cause the poor MTI performance. In this paper, we propose a bistatic configuration design method to jointly optimize the MTI and SAR imaging performance for GEO SA-BSAR MTI system. The relationship between the MTI performance and bistatic configuration parameters is derived analytically and analyzed based on the maximum output signal to clutter and noise ratio (SCNR) criterion. Then, the MTI performance and SAR imaging performance are jointly considered to model the configuration design problem as a multi-objective optimization problem under the constrained condition. Finally, the optimal configuration for GEO SA-BSAR MTI system is given.

Despite the increasing interest in geosynchronous synthetic aperture radar (GEO SAR) systems, the ground moving target indication and motion parameters estimation aspect have never been addressed in GEO SAR scenarios. In this article, we tackle the issue of multiple moving target velocity estimation in GEO SAR. We develop new closed-form expressions that relate both Doppler centroid and Doppler rate to the target motion parameters in GEO SAR by considering the specific features of a geosynchronous orbit. Furthermore, we propose a new velocity estimation algorithm that combines the nonuniform cubic phase function (NU-CPF) algorithm with the newly developed models to estimate the moving target's two velocity components. Moreover, based on the abovementioned technique, we further propose a solution to address the multiple moving targets' problem. Simulation results along with an estimation accuracy analysis are provided to demonstrate the effectiveness of the proposed multitarget GEO SAR velocity estimation technique.

The single-pass geosynchronous synthetic aperture radar interferometry (GEO InSAR) adopts the formation of a slave satellite accompanying the master satellite, which can reduce the temporal decorrelation caused by atmospheric disturbance and observation time gap between repeated tracks. Current formation design methods for spaceborne SAR are based on the Relative Motion Equation (RME) in the Earth-Centered-Inertial (ECI) coordinate system (referred to as ECI-RME). Since the Earth rotation is not taken into account, the methods will lead to a significant error for the baseline calculation while applied to formation design for GEO InSAR. In this paper, a formation design method for single-pass GEO InSAR based on Coordinate Rotational Transformation (CRT) is proposed. Through CRT, the RME in Earth-Centered-Earth-Fixed (ECEF) coordinate system (referred to as ECEF-RME) is derived. The ECEF-RME can be used to describe the accurate baseline of close-flying satellites for different orbital altitudes, but not limited to geosynchronous orbit. Aiming at the problem that ECEF-RME does not have a regular geometry as ECI-RME does, a numerical formation design method based on the minimum baseline error criterion is proposed. Then, an analytical formation design method is proposed for GEO InSAR, based on the Minimum Along-track Baseline Criterion (MABC) subject to a fixed root mean square of the perpendicular baseline. Simulation results verify the validity of the ECEF-RME and the analytical formation design method. The simulation results also show that the proposed method can help alleviate the atmospheric phase impacts and improve the retrieval accuracy of the digital elevation model (DEM) compared with the ECI-RME-based approach.

The work investigates staggered and random PRF (Pulse Repetition Frequency) strategies for a close formation of small Synthetic Aperture Radar (SAR) satellites operating in a multistatic configuration. The satellites are positioned within a fraction of the along-track critical baseline, hence allowing for the application of Displaced Phase Center image formation approaches. The performance of regular and random pulse sampling schemes is in particular assessed for a single-input multiple-output (SIMO) S-Band constellation, whose feasibility is further analyzed in relation to the number of satellites and their antenna size.

Ambiguities in short-baseline ATI interferometry need to be treated not as noise that lowers the coherence, but as a source of bias. A mathematical formulation of the interferometric ambiguity model is given, and an approach to correct ambiguities is proposed and illustrated with simulation results.

The sub-satellite track of geosynchronous synthetic aperture radar (GEO SAR) presents the figure "8" or "O", which causes the great changes of platform motion direction and the different projection of anisotropic irregularities along the line-of-sight (LOS) direction. Due to the almost equal angle velocity to that of Earth, the GEO SAR has smaller ionospheric penetration point (IPP) scanning velocity which is much smaller to the counterpart of the low earth orbit SAR (LEO SAR) while is comparable to the drifting velocity of irregularities, which will affect the effective azimuthal velocity. These facts lead to the consequence that the satellite signals from the GEO SAR would become more vulnerable when they are transmitted in the environment where the ionospheric scintillation occurs. However, few works are focused on these mentioned issues towards the GEO SAR system. In this paper, the impacts of ionospheric scintillation on GEO SAR imaging will be analyzed considering the anisotropy and drifting velocity of irregularities. The anisotropy and drifting velocity effects can both originate from the effect on power spectral density (PSD) of phase screen which is used to model the ionospheric scintillation effects. Based on the data from international geomagnetic reference field (IGRF) and satellite tool kit (STK), the GEO SAR imaging simulations for different GEO SAR orbital configurations and positions are carried out. The simulation results demonstrate that the anisotropy and the drifting velocity of irregularities will cause the changes of stripe direction and affect the quality of GEO SAR images.

This paper debates the performance of the HARMONY (STEREOID) ESA EE-10 candidate mission in measuring the two-dimensional (2D) terrain deformation. Thanks to its Stereo configuration, where the two passive spacecrafts span a large along-track baseline centered on the Sentinel-1 satellite, a large observation angle diversity in azimuth can be achieved. This theoretically leads to promising deformation performance in the north-south direction component, which will play an extremely important role for the surface displacement analysis in the future.