The European Space Agency's (ESA) BIOMASS mission is a pioneering Earth observation satellite mission launched on April 29, 2025. Utilizing a P-band synthetic aperture radar (SAR), the objective of BIOMASS is to deliver estimates of above-ground forest biomass, forest height (FH), and forest disturbance (FD), with unprecedented accuracy. The mission's primary scientific goal is to quantify the distribution and changes in forest biomass, thereby reducing uncertainties in carbon flux estimates and informing climate models. The satellite's advanced instrumentation and innovative approach allow it to penetrate dense forest canopies, capturing data even in challenging environments. The mission will operate in two distinct phases: the tomographic phase and the interferometric phase, which will support polarimetric interferometric SAR (Pol-InSAR) and tomographic SAR (TomoSAR) processing. Additionally, BIOMASS will provide valuable observational data for ice sheets, deserts, the ionosphere, below canopy topography, and other domains.

In the framework of the comparison of Synthetic Aperture Radar (SAR) imagery from the Magellan space mission and the VISAR and VenSAR radar instruments which will be onboard the forthcoming VERITAS and EnVision missions to Venus, the problem of the disparity between the resolutions of the images arises when attempting to define a test statistic with which to detect changes. Reliable change detection requires equivalent spatial resolutions which, for the two different images, inevitably involves different equivalent number of looks after speckle-reduction processing. This study presents a method to address this scenario using a Generalized Beta Prime Distribution as a probability density function (PDF) which is fit to the histogram of the ratio between the two intensity images. The work demonstrates and verifies the properties of the function, highlights its most useful traits, and elaborates on the mathematical procedure required to achieve a meaningful change detection in line with the classic theory of equal number of looks. The results show that the method accurately describes the ratio histogram of two SAR intensity images with different number of looks. Furthermore, they demonstrate the adaptability of the method to the presence of high pixel correlation between the images, and validate its robustness in the presence of textural complexity when the texture patterns of the images are similar.

This paper presents an electromagnetic propagation and scattering model that has been developed to address the feasibility of mapping the basal topography of the Antarctic ice shelves with space-based SAR tomography. Modeling results suggest that volume clutter is a major challenge that can only partly be alleviated by increasing the number of TomoSAR orbits.

Information on ice sheet crystal orientation fabrics (COF) is valuable, as it tells about past deformation and current flow properties, which can potentially improve ice flow modelling. Due to the birefringence of ice sheets, ice sounding radars can provide the desired COF information if multiple polarizations are employed. Until now such data have mostly been acquired with ground-based single-polarization radars by rotating the antenna between multiple data takes. In this paper, a novel method based on polarimetric radar data, polarization synthesis, and complex coherence is proposed and tested with data from an airborne, polarimetric ice sounding radar. The results are consistent with ice core analyses, ice sheet physics, and existing radar results.

Ice crystal orientation fabric (COF) records information about past ice-sheet deformation and influences the present-day flow of ice. Polarimetric radar sounding provides a means to infer anisotropic COF patterns due to the associated birefringence of polar ice. Here, we develop a polarimetric coherence (phase-based) method to determine horizontal properties of the COF. The method utilizes the azimuth and depth dependence of the vertical gradient of the hhvv coherence phase to infer the dielectric principal axes and birefringence, which are then related to the second-order fabric orientation tensor. Specifically, under the assumption that one of the orientational eigenvectors is vertical, we can determine the horizontal eigenvectors and the difference between the horizontal eigenvalues (a measure of horizontal fabric asymmetry). The method exploits single-polarized data acquired with varying antenna orientation. It applies to ground-based “multi-polarization” surveys and is demonstrated using data acquired by Center for Remote Sensing of Ice Sheets (CReSIS) using Multi-Channel Coherent Radar Depth Sounder (MCRDS) from the North Greenland Eemian Ice Drilling (NEEM) ice core region in Greenland. The analysis is validated using a combination of polarimetric matrix backscatter simulations and comparison with COF data from the NEEM ice core. The results are consistent with a conventional model of ice deformation at an ice divide where a lateral tension component is present, with minor horizontal COF asymmetry and the greatest horizontal concentration of crystallographic axes orientated near parallel to the ice divide.

Space-based radio echo sounding of the continental ice sheets can potentially provide full coverage with uniform sampling and data quality as well as detection of change in environmentally sensitive areas. This paper addresses the feasibility of sounding the Antarctic ice sheets with a space-based P-band radar. The assessment firstly makes use of an electromagnetic model of the ice sheets where the key model parameters are determined from data that have been acquired in Antarctica with an airborne P-band ice sounding radar. The performance of a space-based radar with a nadir-looking geometry but otherwise similar to ESA's Biomass SAR, is then simulated for a set of Antarctic scenarios that are defined based on the statistics of key ice regions. It is found that in about 2/3 of the simulation scenarios clutter and/or thermal noise will obscure the echo from the ice bed.

The BIOMASS mission gives completely new challenges in external calibration arising from the orbital pattern needed for the tomographic and Pol-InSAR phases of the mission, the strong effects of the ionosphere at P-band, and the lack of pre-existing P-band data except over very limited parts of the globe. Together these create problems that can only be solved by combining infrequent visits to instrumented calibration sites with systematic exploitation of the properties of distributed targets and targets of opportunity. Proposed approaches to performing radiometric and polarimetric calibration are described, together with meeting geolocation accuracy requirements.

The Earth Explorer Biomass mission will provide the scientific community with accurate maps of tropical, temperate and boreal forest biomass, including height and disturbance patterns. This information is urgently needed to improve our understanding of the global carbon cycle and to reduce uncertainties in the calculation of carbon stocks and fluxes associated to the terrestrial biosphere. It is also crucial for approaches to managing climate, such as the UNFCCC initiative known as Reducing Emissions through Degradation and Deforestation (REDD+), aimed at climate change mitigation through conservation and better management of tropical forests The required measurements are forest biomass and forest height at resolution of 200 m, and detection of deforestation at 50 m. Global maps of biomass are required with accuracy of 20% (or l0 t ha-1when above-ground biomass are less than 50 t ha-1). To achieve this Biomass will be implemented as a P-band SAR mission. It will exploit the unique sensitivity of P-band SAR together with advanced retrieval methods including polarimetric interferometry (Pol-InSAR) and SAR tomography to measure biomass, height and disturbances across the entire biomass range every 6 months. The mission will also support important secondary objectives, including sub-surface imaging in arid zones, production of a bare-earth DTM and ice applications.

Data-driven classification algorithms have proved to do well for automatic target recognition (ATR) in synthetic aperture radar (SAR) data. Collecting data sets suitable for these algorithms is a challenge in itself as it is difficult and expensive. Due to the lack of labeled data sets with real SAR images of sufficient size, simulated data play a big role in SAR ATR development, but the transferability of knowledge learned on simulated data to real data remains to be studied further. In this letter, we show the first study of Transfer Learning between a simulated data set and a set of real SAR images. The simulated data set is obtained by adding a simulated object radar reflectivity to a terrain model of individual point scatters, prior to focusing. Our results show that a Convolutional Neural Network (Convnet) pretrained on simulated data has a great advantage over a Convnet trained only on real data, especially when real data are sparse. The advantages of pretraining the models on simulated data show both in terms of faster convergence during the training phase and on the end accuracy when benchmarked on the Moving and Stationary Target Acquisition and Recognition data set. These results encourage SAR ATR development to continue the improvement of simulated data sets of greater size and complex scenarios in order to build robust algorithms for real life SAR ATR applications.

In this paper, we analyze the direction-of-arrival (DOA) of the ice-sheet data collected over Jakobshavn Glacier with the airborne Multichannel Radar Depth Sounder (MRDS) during the 2006 field season. We extracted weak ice-bed echoes buried in signals scattered by the rough surface of the fast-flowing Jakobshavn Glacier by analyzing the DOA of signals received with a five-element receive-antenna array. This allowed us to obtain ice thickness information, which is a key parameter when generating bed topography of glaciers. We also estimated ice-bed roughness and bed slope from the combined analysis of the DOA and radar waveforms. The bed slope is about 8° and the roughness in terms of rms slope is about 16°.

In this paper the electromagnetic properties of the continental ice sheets are modeled with the aim of assessing the feasibility of space-based radar ice sounding at P-band. The paper focuses on the model itself, as distinct from the estimation of model parameters and the exploitation of the model. The scattering and propagation properties of the ice sheets are characterized using an empirical approach. The model comprises surface scattering from the air/ice interface and the ice/bed interface as well as volume scattering from the firn and the ice. Also specular reflection from the internal layers is modeled. In combination with simple models of the ice sounding radar and the geometry, the electromagnetic model estimates the return waveform and its constituents, such that sensitivity and clutter masking issues can be assessed.

In this paper, first results concerning the characterization of the subsurface of ice sheets and glaciers through single and multipolarization synthetic aperture radar (SAR) tomography (TomoSAR) are illustrated. To this aim, the processing of data acquired in the framework of the European Space Agency IceSAR 2012 campaign is discussed. IceSAR 2012 was conceived so as to support the secondary objectives of the future Earth Explorer mission BIOMASS, which will be a SAR instrument with media penetration capabilities due to the use of the P-band frequency. In this regard, a tomographic study of ice was motivated by the fact that cryospheric remote sensing is of fundamental importance in order to understand more in depth the morphology and the dynamic processes regulating ice sheets. The main objective of the tomographic experiment of the campaign herein discussed was indeed to assess the capability of P-band SAR to retrieve any information about ice subsurface structure. Imaging has been achieved through TomoSAR techniques, applied to airborne multibaseline data acquired in the southwest of Greenland. Different imaging approaches are compared, and the main results achieved are presented: It is found that scattering in the upper layers of glacial subsurface can be achieved up to an extent of about 20-60 m, conditional on the different types of glaciological zone observed. Moreover, clear morphological structures have been found beneath the ice surface at one of the investigated sites.

Ice sounding radars are able to measure ice sheets by profiling their glaciological features from the surface to the bedrock. The current airborne and, in particular, future space-based systems are suffering from off-nadir surface clutter, which can mask the depth signal of interest. The most recent surface clutter suppression techniques are based on multi-phase-center systems combined with sophisticated coherent postprocessing. The performance of the techniques can be improved by accurate direction-of-arrival (DOA) estimates of the surface clutter. This paper deals with data-driven DOA estimation for surface clutter signals, which includes a formulation of the mathematical foundation of spatial aliasing. DOA estimation is applied to data acquired with the P-band POLarimetric Airborne Radar Ice Sounder at the Jutulstraumen Glacier, Antarctica. The effects of spatial aliasing related to a large phase center spacing are analyzed, and an unwrapping procedure is presented and applied to the data. Finally, DOA estimation of full-scene data is analyzed and used to show an along-track and incidence (off-nadir) angle dependent variation of the effective scattering center of the surface return, which is caused by a varying penetration depth.

European Space Agency's Cryosat-2 comes with the first ever SAR (Synthetic Aperture Radar) altimeter onboard a satellite. In this work precise sea surface heights and gravity fields are determined using Cryosat-2 SAR data. These determinations through satellite altimetry are difficult in the Arctic because of the presence of sea ice and coastal areas. Traditional retrackers generate erroneous results due to the superposition of echoes from the sea surface and sea ice. This work establishes an adaptation of the SAMOSA3 retracker as a suitable candidate; this model is improved and customized for the Arctic. Through this research it has been demonstrated that the SAMOSA3 retracker has a better performance as compared to other SAR retrackers when sea surface height and gravity field determination needs to be done. The performance evaluation of the SAMOSA3 retracker as compared to other retrackers has been done using sea surface height anomaly method and gravity field anomaly method.

Due to the increased melting of ice sheets over the last decades, monitoring of ice dynamics and structure with remote sensing instruments is of extreme importance to achieve a deeper insight on related environmental issues. The study presented in this paper documents an attempt of mapping ice structure with P-band SAR tomography. First results from ESA IceSAR 2012 campaign carried out in south-west Greenland are presented. It is found that significant penetration in the upper layers of glacial subsurface can be achieved up to an extent of about 20-60 m, conditional on the different type of glaciological zone observed.

Glacier and ice sheet dynamics are currently mapped with X-, C-, and L-band SAR. With the prospect of a P-band SAR, Biomass, to be launched within the next decade it is interesting to look into the potential of P-band for ice velocity mapping. In this paper first results are presented. Airborne P-band SAR data have been acquired in Greenland, and both offset tracking and DInSAR have been applied to the full resolution data as well as to data degraded to the resolution of Biomass. Generally, ice velocity maps are successfully generated, but in the ablation zone, DInSAR fails in the melt season at both resolutions, while feature tracking fails at the coarser resolution.

In February 2011, the Polarimetric Airborne Radar Ice Sounder (POLARIS) was flown in Antarctica in order to assess the feasibility of a potential space-based radar ice sounding mission. The campaign has demonstrated that the basal return is detectable in areas with up to 3 km thick cold ice, in areas with up to several hundred meters thick warm shelf ice, and in areas with up to 700 m thick crevassed glacier ice. However, major gaps in the basal return are observed, presumably due to excessive absorption, scattering from ice inclusions in the firn, low basal reflectivity, and the masking effect of the surface clutter. Internal layers are observed down to depths exceeding 2 km. The polarimetric data show that the internal layers are strongly anisotropic at a ridge, where the ice flow is supposed to be highly unidirectional. In case of space-based ice sounding, the antenna pattern cannot offer sufficient surface clutter suppression, but improved clutter suppression has been demonstrated with novel multi-phase-center techniques.

The BIOMASS mission is designed to provide unique information on the biomass in the world's forests at spatial and temporal resolutions suitable for characterizing their dynamics and their contribution to carbon cycle estimates. To achieve this it combines biomass estimates from direct inversion of polarimetric backscattering coefficients with Pol-InSAR forest height estimates. The mission will also support important secondary objectives, including sub-surface imaging in arid zones, production of a bare-earth DTM and ice applications, and is optimized to be robust against environmental and ionospheric disturbances.

For polar ice sheets, valuable stress and strain information can be deduced from crystal orientation fabrics (COF) and their prevailing c-axis alignment. Polarimetric radio echo sounding is a promising technique to measure the anisotropic electromagnetic propagation and reflection properties associated with COFs. In this paper, fully polarimetric P-band data acquired with the airborne POLARIS system near the ice divide of the Greenland ice sheet are analyzed. Based on a simple electromagnetic model, these data are interpreted, and a pronounced birefringence is found.

SAR processing can be applied to ice sounder data to improve along-track resolution and clutter suppression. This paper presents a time-domain back-projection technique for SAR focusing of ice sounder data. With this technique, variations in flight track and ice surface slope can be accurately accommodated at the expense of computation time. The back-projection algorithm can be easily parallelized however, and can advantageously be implemented on a graphics processing unit (GPU). Results from using the back-projection algorithm on POLARIS ice sounder data from North Greenland shows that the quality of data is improved by the processing, and the performance of the GPU implementation allows for very fast focusing.

The European Space Agency (ESA) released a Call for Proposals for the next Earth Explorer Core Mission in March 2005, with the aim to select the 7thEarth Explorer (EE-7) mission for launch in the next decade. Twenty-four proposals were received and subject to scientific and technical assessment. Six candidate missions were selected and further investigated in the preliminary feasibility studies (Phase 0). One of these missions is BIOMASS, which has recently been selected to proceed to Phase-A. BIOMASS is a response to the urgent need for greatly improved mapping of global biomass and the lack of any current space systems capable of addressing this need.

For polar ice sheets valuable stress and strain information can be deduced from the crystal orientation fabric (COF) and its prevailing c-axis alignment. Polarimetric radio echo sounding is a promising technique to measure the anisotropic electromagnetic propagation and reflection properties associated with the COF. In this paper, dual-polarized P-band data acquired with the airborne POLARIS system near the ice divide of the Greenland ice sheet are analyzed. The internal layers in the uppermost few hundred meters of the ice sheet look the same at HH and VV polarizations, whereas the layering differs further down. Accordingly, the magnitude of the complex HH-VV correlation coefficient decreases with depth and, interestingly, the phase gives evidence of polarization dependent reflection and birefringence effects.

ESA has assigned the Technical University of Denmark to develop an airborne P-band ice sounding radar demonstrator. The intention is to obtain a better understanding of the electromagnetic properties of the Antarctic ice sheet at P-band and to test novel ice sounding techniques in preparation for a potential spaceborne ice sounding radar. The airborne system is a coherent, high-resolution and fully polarimetric radar. Aperture synthesis is applied in the along track direction and an experimental surface clutter suppression technique based on a multi-aperture antenna can be applied in the across track direction. In May 2008, a proof-of-concept campaign was organized in Greenland, where data were acquired over the ice sheet. The system proved capable of detecting the bedrock under 3 km thick ice and of mapping the internal ice layers down to a depth of at least 1.3 km. In this paper, the system concept is outlined and first results are presented.

C- and L-band airborne synthetic aperture radar (SAR) imagery acquired at like- and cross-polarizations over sea ice under winter conditions is examined with the objective to study the discrimination between level ice and ice deformation features. High-resolution low-noise data were analyzed in the first paper. In this second paper, the main topics are the effects of spatial resolution and signal-to-noise ratio. Airborne high-resolution SAR scenes are used to generate a sequence of images with increasingly coarser spatial resolution from 5 to 25 m, keeping the number of looks constant. The signal-to-noise ratio is varied between typical noise levels for airborne imagery and satellite data. Areal fraction of deformed ice and average deformation distance are determined for each image product. At L-band, the retrieved values of the areal fraction get larger as the image resolution is degraded. The areal fraction at C-band remains constant. The retrieved average distance between deformation features increases both at C- and L-bands as the image resolution gets coarser. The influence of noise becomes noticeable if its level is equal or larger than the average intensity backscattered from the level ice. The retrieval of deformation parameters using simulated images that resemble ERS-2 SAR, Envisat ASAR, and ALOS PALSAR data products is discussed. Basic differences between real and simulated ERS-2 SAR images are analyzed.

This paper presents the top-level design of an airborne, P-band ice sounding radar under development at the Technical University of Denmark. The ice sounder is intended to provide more information on the electromagnetic properties of the Antarctic ice sheet at P-band. A secondary objective is to test new ice sounding techniques, e.g. polarimetry, synthetic aperture processing, and coherent clutter suppression. A system analysis involving ice scattering models confirms that it is feasible to detect the bedrock through 4 km of ice and to detect deep ice layers. The ice sounder design features a digital signal generator, a microstrip antenna array, a conventional RF-architecture with a central transmitter, four receivers, and internal calibration loops. In 2008 the first data acquisition campaign will take place in Greenland.

Natural media like cold-land ice, vegetation, and dry sand are subject to a substantial penetration at microwave frequencies. For such media, the synthetic aperture radar (SAR) phase center is located below the surface, and consequently, the surface elevation determined with SAR interferometry (InSAR) is biased downward. For infinitely deep uniform volumes, the elevation bias is often equated with the penetration depth, but in this paper, it is shown that the two quantities generally differ. The interferometric bias is approximately equal to the two-way power-penetration depth if the latter is small compared to the ambiguity height, but for increasing penetration depth, the bias approaches one quarter of the ambiguity height. Consequently, no phase wrapping results even if the penetration depth exceeds the ambiguity height. The ratio of the InSAR elevation bias to the ambiguity height depends only on the ratio of the penetration depth to the ambiguity height, and the bias can be expressed in terms of the InSAR coherence magnitude, which makes it possible to correct the InSAR surface elevation for the bias. The volume depth can be considered infinite if it exceeds the penetration depth by a factor of two to five and if the surface scattering from the top and the bottom of the volume is negligible.

In this paper, we present a quantitative comparison of L- and C-band airborne synthetic aperture radar imagery acquired at like- and cross-polarizations over deformed sea ice under winter conditions. The parameters characterizing the deformation state of the ice are determined at both radar bands and at different polarizations. The separation of deformed and level ice is based on a target detection technique. The threshold is set such that image pixels with intensities equal to or larger than the highest 2% of the level-ice intensity distribution are classified as deformed ice, independent of the radar configuration and ice conditions. Optical imagery of sufficient quality for comparison is available only in a very few cases. To characterize the deformation state, the areal fraction of deformation features and the average distance between these features are evaluated. The values obtained for both parameters are very sensitive to the radar frequency. Aeral fractions are larger, and average distances are smaller at L-band than at C-band because of the much higher intensity contrast between the deformed and level ice at L-band. The differences between polarizations at one radar band are smaller but not always negligible. In comparison to optical images, it was observed that deformed-ice areas can be distinguished from level ice over their whole length and extension at L-band, whereas at C-band, often, only prominent parts are visible.

A novel ambiguity suppression technique is proposed. Range ambiguities in synthetic aperture radar (SAR) images are eliminated with an azimuth filter after having applied an azimuth phase modulation to the transmitted pulses and a corresponding demodulation to the received pulses. The technique excels by actually eliminating the ambiguities rather than just defocusing them as most other techniques do. This makes the proposed technique applicable to distributed targets. The range ambiguity suppression permits the pulse repetition frequency (PRF) to exceed the upper limit otherwise defined by the antenna elevation dimension. The fundamental antenna area constraint still applies, but the PRF can be chosen with more freedom. In addition to ambiguity suppression, potential applications include nadir return elimination and signal-to-noise ratio improvement

This paper presents first results obtained from the analysis of fully polarimetric / interferometric land ice data. It is demonstrated that L-band data from an ice cap located in the percolation zone of Greenland show a strong polarization de- pendent interferometric behavior. No model-based inversion has yet been attempted, but a radar reflector deployed on the ice surface suggests an effective penetration depth of about 15 m.

EMISAR is a fully polarimetric, dual frequency (L- and C-band) SAR system designed for remote sensing applications. The data are usually processed to 2/spl times/2 m resolution. The system has the capability of C-band cross-track single-pass interferometry and fully polarimetric repeat-pass interferometry. The SAR is operated at high altitudes on a Gulfstream G-3 jet aircraft. The system is very well calibrated and has low sidelobes and low cross-polar contamination. Digital technology has been utilized to realize a flexible and highly stable radar with variable resolution, swath width, and imaging geometry. Thermal control and several calibration loops have been built into the system to ensure system stability and absolute calibration. Accurately measured antenna gains and radiation patterns are included in the calibration. The processing system is developed to support data calibration, which is the key to most of the current applications.

Cross-calibration is a fully automated algorithm for calibration of interferometric synthetic aperture radar (IFSAR) data. It has been developed for single-pass interferometry, but the principles may be applicable to multi-pass interferometry, too. The algorithm is based on natural distributed targets and it excels by neither requiring surveyed ground control points nor dedicated calibration scenes. However, the parameters to be calibrated must be stable during mapping. The algorithm has been applied to data from the Danish airborne SAR, EMISAR, and the performance has been assessed. The algorithm appears to be fairly robust with respect to the terrain type. However, the result of the calibration may deteriorate if the terrain elevation, as measured with the SAR, changes systematically with the incidence angle or the aspect angle.

EMISAR is a high-resolution (2/spl times/2 m), fully polarimetric, dual-frequency (L- and C-band) synthetic aperture radar (SAR) system designed for remote-sensing applications. The SAR is operated at high altitudes on a Gulfstream G-3 jet aircraft. The system is very well calibrated and has low sidelobes and low cross-polar contamination. Digital technology has been utilized to realize a flexible and highly stable radar with variable resolution, swath width, and imaging geometry. Thermal control and several calibration loops have been built into the system to ensure system stability and absolute calibration. Accurately measured antenna gains and radiation patterns are included in the calibration. The processing system is developed to support data calibration, which is the key to most of the current applications. Recent interferometric enhancements are important for many scientific applications.

High-resolution synthetic aperture radar (SAR) imaging requires the motion of the radar platform to be known very accurately. Otherwise, phase errors are induced in the processing of the raw SAR data, and bad focusing results. In particular, a constant error in the measured along-track velocity or the cross-track acceleration leads to a phase error that varies quadratically over the synthetic aperture. The process of estimating this quadratic phase error directly from the radar data is termed autofocus. A novel autofocus algorithm with a computational complexity which is at least an order of magnitude lower than that of other algorithms providing comparable accuracies is presented. The algorithm has been tested on data from the Danish Airborne SAR, and the performance is compared with that of the traditional map drift algorithm.>

In January 1986, the design of a high-resolution airborne C-band synthetic aperture radar (SAR) started at the Electromagnetics Institute of the Technical University of Denmark. The initial system test flights took place in November and December 1989. The authors describe the design of the system, its implementation, and its performance. They show how digital technology has been utilized to realize a very flexible radar with variable resolution, swath-width, and imaging geometry. The motion-compensation algorithms implemented to obtain the high resolution and the special features built into the system to ensure proper internal calibration are outlined. The data processing system, developed for image generation and quality assurance, is sketched, with special emphasis on the flexibility of the system.>