The synthetic aperture radar interferometry (InSAR) technique has been applied in monitoring the deformation of infrastructures, such as bridges, highways, railways and subways. Persistent scatterer (PS)-InSAR is one of the InSAR techniques, which utilises persistent scatterers to derive long-term displacements. This study applied time-series methods to post-process the PS-InSAR-derived time-series displacements with the use of 86 Sentinel-1A acquisitions spanning from 6 January 2018 to 27 November 2020. Empirical mode decomposition (EMD) and seasonal and trend decomposition using loess (STL) were combined to estimate the seasonal component of the total time-series displacements. Then, a temperature correlation map was generated by correlating the seasonal component with the temperature variation. Results show that the thermal expansion phenomenon is pronounced on the buildings of the Zhuhai–Macao Passenger Terminal as well as the bridge and road connecting to the Hong Kong International Airport (HKIA), while it is less obviously observed at the main Hong Kong-Zhuhai-Macao Bridge (HZMB). In addition, sudden changes between subsidence and uplift can be detected through the p-values derived by applying the augmented Dickey-Fuller (ADF) test to the residual signals after removing the linear and seasonal components from the original ones.

The proper use of distributed scatterer (DS) can improve both the density and quality of synthetic aperture radar (SAR) interferometry (InSAR) measurements. A critical step in DS interferometry (DSI) is the restoration of a consistent phase series from SAR interferogram stacks. Most state-of-the-art algorithms adopt an approximate likelihood function to calculate the likelihood by replacing the true coherence matrix with its estimation, more specifically, the sample coherence matrix (SCM). However, this approximation has a drawback in that the coherence estimates are greatly biased when the coherence is low. In this study, we derive a new likelihood function without such an approximation. Accordingly, a DSI framework using this function for phase estimation and point selection is provided. In this framework, the new likelihood function serves as a cost function for phase estimation and a quality measure for DS selection. Its performance is investigated by experiments in a simulation study and a real-world case study using Sentinel-1 data over Shenzhen airport in China. The results reveal that the proposed DSI framework outperforms the existing state-of-the-art approaches in different scenarios, in terms of providing a more accurate estimation and improving DS density and coverage.

Englacial layering reflects ice dynamics within the ice bodies, which improves understanding of ice flow variation, past accumulation rates and vertical flows transferring between the surface and the underlying bedrock. The internal layers can be observed by using Radar Echo Sounding (RES), such as the Multi-channel Coherent Radar Depth Sounder (MCoRDS) used in NASA’s Operation IceBridge (OIB) mission. Since the 1960s, the accumulation of the RES data has prompted the development of automated methods to extract the englacial layers. In this study, we propose a new automated method that combines peak detection methods, namely the CWT-based peak detection or the Automatic Phase Picker (APP), with a Hough Transform (HT) to trace boundaries of englacial layers. For CWT-based peak detection, we test it using two different wavelets. The proposed method is tested with twelve MCoRDS radio echograms, which are acquired south of the Northern Greenland Eemian (NEEM) ice drilling site, where the folding of ice layers was observed. The method is evaluated in comparison to the isochrones that were extracted in an independent study. In comparison, the proposed new automated method can restore more than 70% of the englacial layers. This new automated layer-tracing method is publicly available on github.

The High-Resolution Imaging Science Experiment (HiRISE) onboard the Mars Reconnaissance Orbiter provides remotely sensed imagery at the highest spatial resolution at 25–50 cm/pixel of the surface of Mars. However, due to the spatial resolution being so high, the total area covered by HiRISE targeted stereo acquisitions is very limited. This results in a lack of the availability of high-resolution digital terrain models (DTMs) which are better than 1 m/pixel. Such high-resolution DTMs have always been considered desirable for the international community of planetary scientists to carry out fine-scale geological analysis of the Martian surface. Recently, new deep learning-based techniques that are able to retrieve DTMs from single optical orbital imagery have been developed and applied to single HiRISE observational data. In this paper, we improve upon a previously developed single-image DTM estimation system called MADNet (1.0). We propose optimisations which we collectively call MADNet 2.0, which is based on a supervised image-to-height estimation network, multi-scale DTM reconstruction, and 3D co-alignment processes. In particular, we employ optimised single-scale inference and multi-scale reconstruction (in MADNet 2.0), instead of multi-scale inference and single-scale reconstruction (in MADNet 1.0), to produce more accurate large-scale topographic retrieval with boosted fine-scale resolution. We demonstrate the improvements of the MADNet 2.0 DTMs produced using HiRISE images, in comparison to the MADNet 1.0 DTMs and the published Planetary Data System (PDS) DTMs over the ExoMars Rosalind Franklin rover’s landing site at Oxia Planum. Qualitative and quantitative assessments suggest the proposed MADNet 2.0 system is capable of producing pixel-scale DTM retrieval at the same spatial resolution (25 cm/pixel) of the input HiRISE images.

The lack of adequate stereo coverage and where available, lengthy processing time, various artefacts, and unsatisfactory quality and complexity of automating the selection of the best set of processing parameters, have long been big barriers for large-area planetary 3D mapping. In this paper, we propose a deep learning-based solution, called MADNet (Multi-scale generative Adversarial u-net with Dense convolutional and up-projection blocks), that avoids or resolves all of the above issues. We demonstrate the wide applicability of this technique with the ExoMars Trace Gas Orbiter Colour and Stereo Surface Imaging System (CaSSIS) 4.6 m/pixel images on Mars. Only a single input image and a coarse global 3D reference are required, without knowing any camera models or imaging parameters, to produce high-quality and high-resolution full-strip Digital Terrain Models (DTMs) in a few seconds. In this paper, we discuss technical details of the MADNet system and provide detailed comparisons and assessments of the results. The resultant MADNet 8 m/pixel CaSSIS DTMs are qualitatively very similar to the 1 m/pixel HiRISE DTMs. The resultant MADNet CaSSIS DTMs display excellent agreement with nested Mars Reconnaissance Orbiter Context Camera (CTX), Mars Express’s High-Resolution Stereo Camera (HRSC), and Mars Orbiter Laser Altimeter (MOLA) DTMs at large-scale, and meanwhile, show fairly good correlation with the High-Resolution Imaging Science Experiment (HiRISE) DTMs for fine-scale details. In addition, we show how MADNet outperforms traditional photogrammetric methods, both on speed and quality, for other datasets like HRSC, CTX, and HiRISE, without any parameter tuning or re-training of the model. We demonstrate the results for Oxia Planum (the landing site of the European Space Agency’s Rosalind Franklin ExoMars rover 2023) and a couple of sites of high scientific interest.

Spaceborne Synthetic Aperture Radar (SAR) interferometry has been utilised to acquire high-resolution Digital Elevation Models (DEMs) with wide coverage, particularly for persistently cloud-covered regions where stereophotogrammetry is hard to apply. Since the discovery of sand buried drainage systems by the Shuttle Imaging Radar-A (SIR-A) L-band mission in 1982, radar images have been exploited to map subsurface features beneath a sandy cover of extremely low loss and low bulk humidity in some hyper-arid regions such as from the Japanese Earth Resources Satellite 1 (JERS-1) and Advanced Land Observing Satellite/Phased Array type L-band Synthetic Aperture Radar (ALOS/PALSAR). Therefore, we hypothesise that a Digital Elevation Model (DEM) derived by InSAR in hyper-arid regions is likely to represent a subsurface elevation model, especially for lower frequency radar systems, such as the L-band system (1.25 GHz). In this paper, we compare the surface appearance of radar images (L-band and C-band) with that of optical images to demonstrate their different abilities to show subsurface features. Moreover, we present an application of L-band InSAR to measure penetration depths in the eastern Sahara Desert. We demonstrate how the retrieved L-band InSAR DEM appears to be of a consistently 1–2 m lower elevation than the C-band Shuttle Radar Topography Mission (SRTM) DEM over sandy covered areas, which indicates the occurrence of penetration and confirms previous studies.