Sentinel-1 backscatter time series for characterization of evapotranspiration dynamics over temperate coniferous forests

Forests’ ecosystems are an essential part of the global carbon cycle with vast carbon storage potential. These systems are currently under external pressures showing increasing change due to climate change. A better understanding of the biophysical properties of forests is, therefore, of paramount importance for research and monitoring purposes. While there are many biophysical properties, the focus of this study is on the in-depth analysis of the connection between the C-band Copernicus Sentinel-1 SAR backscatter and evapotranspiration (ET) estimates based on in situ meteorological data and the FAO-based Penman–Monteith equation as well as the well-established global terrestrial ET product from the Terra and Aqua MODIS sensors. The analysis was performed in the Free State of Thuringia, central Germany, over coniferous forests within an area of 2452 km2, considering a 5-year time series (June 2016–July 2021) of 6- to 12-day Sentinel-1 backscatter acquisitions/observations, daily in situ meteorological measurements of four weather stations as well as an 8-day composite of ET products of the MODIS sensors. Correlation analyses of the three datasets were implemented independently for each of the microwave sensor’s acquisition parameters, ascending and descending overpass direction and co- or cross-polarization, investigating different time series seasonality filters. The Sentinel-1 backscatter and both ET time series datasets show a similar multiannual seasonally fluctuating behavior with increasing values in the spring, peaks in the summer, decreases in the autumn and troughs in the winter months. The backscatter difference between summer and winter reaches over 1.5 dB, while the evapotranspiration difference reaches 8 mm/day for the in situ measurements and 300 kg/m2/8-day for the MODIS product. The best correlation between the Sentinel-1 backscatter and both ET products is achieved in the ascending overpass direction, with datasets acquired in the late afternoon, and reaches an R2-value of over 0.8. The correlation for the descending overpass direction reaches values of up to 0.6. These results suggest that the SAR backscatter signal of coniferous forests is sensitive to the biophysical property evapotranspiration under some scenarios

Review Article: Global Monitoring of Snow Water Equivalent Using High-Frequency Radar Remote Sensing

Seasonal snow cover is the largest single component of the cryosphere in areal extent, covering an average of 46 × 106 km2 of Earth\u27s surface (31 % of the land area) each year, and is thus an important expression and driver of the Earth\u27s climate. In recent years, Northern Hemisphere spring snow cover has been declining at about the same rate (∼ −13 % per decade) as Arctic summer sea ice. More than one-sixth of the world\u27s population relies on seasonal snowpack and glaciers for a water supply that is likely to decrease this century. Snow is also a critical component of Earth\u27s cold regions\u27 ecosystems, in which wildlife, vegetation, and snow are strongly interconnected. Snow water equivalent (SWE) describes the quantity of water stored as snow on the land surface and is of fundamental importance to water, energy, and geochemical cycles. Quality global SWE estimates are lacking. Given the vast seasonal extent combined with the spatially variable nature of snow distribution at regional and local scales, surface observations are not able to provide sufficient SWE information. Satellite observations presently cannot provide SWE information at the spatial and temporal resolutions required to address science and high-socio-economic-value applications such as water resource management and streamflow forecasting. In this paper, we review the potential contribution of X- and Ku-band synthetic aperture radar (SAR) for global monitoring of SWE. SAR can image the surface during both day and night regardless of cloud cover, allowing high-frequency revisit at high spatial resolution as demonstrated by missions such as Sentinel-1. The physical basis for estimating SWE from X- and Ku-band radar measurements at local scales is volume scattering by millimeter-scale snow grains. Inference of global snow properties from SAR requires an interdisciplinary approach based on field observations of snow microstructure, physical snow modeling, electromagnetic theory, and retrieval strategies over a range of scales. New field measurement capabilities have enabled significant advances in understanding snow microstructure such as grain size, density, and layering. We describe radar interactions with snow-covered landscapes, the small but rapidly growing number of field datasets used to evaluate retrieval algorithms, the characterization of snowpack properties using radar measurements, and the refinement of retrieval algorithms via synergy with other microwave remote sensing approaches. This review serves to inform the broader snow research, monitoring, and application communities on progress made in recent decades and sets the stage for a new era in SWE remote sensing from SAR measurements

Biomass Representation in Synthetic Aperture Radar Interferometry Data Sets

This work makes an attempt to explain the origin, features and potential applications of the elevation bias of the synthetic aperture radar interferometry (InSAR) datasets over areas covered by vegetation. The rapid development of radar-based remote sensing methods, such as synthetic aperture radar (SAR) and InSAR, has provided an alternative to the photogrammetry and LiDAR for determining the third dimension of topographic surfaces. The InSAR method has proved to be so effective and productive that it allowed, within eleven days of the space shuttle mission, for acquisition of data to develop a three-dimensional model of almost the entire land surface of our planet. This mission is known as the Shuttle Radar Topography Mission (SRTM). Scientists across the geosciences were able to access the great benefits of uniformity, high resolution and the most precise digital elevation model (DEM) of the Earth like never before for their a wide variety of scientific and practical inquiries. Unfortunately, InSAR elevations misrepresent the surface of the Earth in places where there is substantial vegetation cover. This is a systematic error of unknown, yet limited (by the vertical extension of vegetation) magnitude. Up to now, only a limited number of attempts to model this error source have been made. However, none offer a robust remedy, but rather partial or case-based solutions. More work in this area of research is needed as the number of airborne and space-based InSAR elevation models has been steadily increasing over the last few years, despite strong competition from LiDAR and optical methods. From another perspective, however, this elevation bias, termed here as the “biomass impenetrability”, creates a great opportunity to learn about the biomass. This may be achieved due to the fact that the impenetrability can be considered a collective response to a few factors originating in 3D space that encompass the outermost boundaries of vegetation. The biomass, presence in InSAR datasets or simply the biomass impenetrability, is the focus of this research. The report, presented in a sequence of sections, gradually introduces terminology, physical and mathematical fundamentals commonly used in describing the propagation of electromagnetic waves, including the Maxwell equations. The synthetic aperture radar (SAR) and InSAR as active remote sensing methods are summarised. In subsequent steps, the major InSAR data sources and data acquisition systems, past and present, are outlined. Various examples of the InSAR datasets, including the SRTM C- and X-band elevation products and INTERMAP Inc. IFSAR digital terrain/surface models (DTM/DSM), representing diverse test sites in the world are used to demonstrate the presence and/or magnitude of the biomass impenetrability in the context of different types of vegetation – usually forest. Also, results of investigations carried out by selected researchers on the elevation bias in InSAR datasets and their attempts at mathematical modelling are reviewed. In recent years, a few researchers have suggested that the magnitude of the biomass impenetrability is linked to gaps in the vegetation cover. Based on these hints, a mathematical model of the tree and the forest has been developed. Three types of gaps were identified; gaps in the landscape-scale forest areas (Type 1), e.g. forest fire scares and logging areas; a gap between three trees forming a triangle (Type 2), e.g. depending on the shape of tree crowns; and gaps within a tree itself (Type 3). Experiments have demonstrated that Type 1 gaps follow the power-law density distribution function. One of the most useful features of the power-law distributed phenomena is their scale-independent property. This property was also used to model Type 3 gaps (within the tree crown) by assuming that these gaps follow the same distribution as the Type 1 gaps. A hypothesis was formulated regarding the penetration depth of the radar waves within the canopy. It claims that the depth of penetration is simply related to the quantisation level of the radar backscattered signal. A higher level of bits per pixels allows for capturing weaker signals arriving from the lower levels of the tree crown. Assuming certain generic and simplified shapes of tree crowns including cone, paraboloid, sphere and spherical cap, it was possible to model analytically Type 2 gaps. The Monte Carlo simulation method was used to investigate relationships between the impenetrability and various configurations of a modelled forest. One of the most important findings is that impenetrability is largely explainable by the gaps between trees. A much less important role is played by the penetrability into the crown cover. Another important finding is that the impenetrability strongly correlates with the vegetation density. Using this feature, a method for vegetation density mapping called the mean maximum impenetrability (MMI) method is proposed. Unlike the traditional methods of forest inventories, the MMI method allows for a much more realistic inventory of vegetation cover, because it is able to capture an in situ or current situation on the ground, but not for areas that are nominally classified as a “forest-to-be”. The MMI method also allows for the mapping of landscape variation in the forest or vegetation density, which is a novel and exciting feature of the new 3D remote sensing (3DRS) technique. Besides the inventory-type applications, the MMI method can be used as a forest change detection method. For maximum effectiveness of the MMI method, an object-based change detection approach is preferred. A minimum requirement for the MMI method is a time-lapsed reference dataset in the form, for example, of an existing forest map of the area of interest, or a vegetation density map prepared using InSAR datasets. Preliminary tests aimed at finding a degree of correlation between the impenetrability and other types of passive and active remote sensing data sources, including TerraSAR-X, NDVI and PALSAR, proved that the method most sensitive to vegetation density was the Japanese PALSAR - L-band SAR system. Unfortunately, PALSAR backscattered signals become very noisy for impenetrability below 15 m. This means that PALSAR has severe limitations for low loadings of the biomass per unit area. The proposed applications of the InSAR data will remain indispensable wherever cloud cover obscures the sky in a persistent manner, which makes suitable optical data acquisition extremely time-consuming or nearly impossible. A limitation of the MMI method is due to the fact that the impenetrability is calculated using a reference DTM, which must be available beforehand. In many countries around the world, appropriate quality DTMs are still unavailable. A possible solution to this obstacle is to use a DEM that was derived using P-band InSAR elevations or LiDAR. It must be noted, however, that in many cases, two InSAR datasets separated by time of the same area are sufficient for forest change detection or similar applications

Physics-based Modeling for High-fidelity Radar Retrievals.

Knowledge of soil moisture on a global scale is crucial for understanding the Earth's water, energy, and carbon cycles. This dissertation is motivated by the need for accurate soil moisture estimates and focuses on the improvement of soil moisture retrieval based on active remote sensing over vegetated areas. It addresses important, but often neglected, aspects in radar imaging: effects related to the ionosphere, multispecies vegetation (heterogeneity at pixel level), and heterogeneity at landscape level. The first contribution is the development of a generalized radar scattering model as an advancement of current radar modeling techniques for vegetated areas at fine-scale pixel level. It consists of realistic representations of multispecies and subsurface soil layer modeling, and includes terrain topography. This modeling improvement allows greater applicability to different land cover types and higher soil moisture retrieval accuracy. Most coarse-scale satellite pixels (km-scale or coarser) contain highly heterogeneous scenes with fine-scale (100 m or finer) variability of soil moisture, soil texture, topography, and vegetation cover. The second contribution is the development of spatial scaling techniques to investigate effects of landscape-level heterogeneity on radar scattering signatures. Using the above radar forward scattering model, which assumes homogeneity over fine scales, tailor-made models are derived for the contribution of fine-scale heterogeneity to the coarse-scale satellite pixel for effective soil moisture retrieval. Finally, the third contribution is the development of a self-contained calibration technique based on an end-to-end radar system model. The model includes ionospheric effects allowing the use of spaceborne radar signals for accurate soil moisture retrieval from lower frequencies, such as L- and P-band. These combined contributions will greatly increase the usability of low-frequency spaceborne radar data for soil moisture retrieval: ionospheric effects are mitigated, landscape level heterogeneity is resolved, and fine-scale scenes are better modeled. These contributions ultimately allow improved fidelity in soil moisture retrieval and are immediately applicable in current missions such as the ongoing AirMOSS mission that observes root-zone soil moisture with a P-band radar at fine-scale resolution (100 m), and NASA's upcoming SMAP spaceborne mission, which will assess surface soil moisture with an L-band radar and radiometer at km-scale resolution (3 km).PHDElectrical EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/107290/1/mburgin_1.pd

Multi-source Remote Sensing for Forest Characterization and Monitoring

As a dominant terrestrial ecosystem of the Earth, forest environments play profound roles in ecology, biodiversity, resource utilization, and management, which highlights the significance of forest characterization and monitoring. Some forest parameters can help track climate change and quantify the global carbon cycle and therefore attract growing attention from various research communities. Compared with traditional in-situ methods with expensive and time-consuming field works involved, airborne and spaceborne remote sensors collect cost-efficient and consistent observations at global or regional scales and have been proven to be an effective way for forest monitoring. With the looming paradigm shift toward data-intensive science and the development of remote sensors, remote sensing data with higher resolution and diversity have been the mainstream in data analysis and processing. However, significant heterogeneities in the multi-source remote sensing data largely restrain its forest applications urging the research community to come up with effective synergistic strategies. The work presented in this thesis contributes to the field by exploring the potential of the Synthetic Aperture Radar (SAR), SAR Polarimetry (PolSAR), SAR Interferometry (InSAR), Polarimetric SAR Interferometry (PolInSAR), Light Detection and Ranging (LiDAR), and multispectral remote sensing in forest characterization and monitoring from three main aspects including forest height estimation, active fire detection, and burned area mapping. First, the forest height inversion is demonstrated using airborne L-band dual-baseline repeat-pass PolInSAR data based on modified versions of the Random Motion over Ground (RMoG) model, where the scattering attenuation and wind-derived random motion are described in conditions of homogeneous and heterogeneous volume layer, respectively. A boreal and a tropical forest test site are involved in the experiment to explore the flexibility of different models over different forest types and based on that, a leveraging strategy is proposed to boost the accuracy of forest height estimation. The accuracy of the model-based forest height inversion is limited by the discrepancy between the theoretical models and actual scenarios and exhibits a strong dependency on the system and scenario parameters. Hence, high vertical accuracy LiDAR samples are employed to assist the PolInSAR-based forest height estimation. This multi-source forest height estimation is reformulated as a pan-sharpening task aiming to generate forest heights with high spatial resolution and vertical accuracy based on the synergy of the sparse LiDAR-derived heights and the information embedded in the PolInSAR data. This process is realized by a specifically designed generative adversarial network (GAN) allowing high accuracy forest height estimation less limited by theoretical models and system parameters. Related experiments are carried out over a boreal and a tropical forest to validate the flexibility of the method. An automated active fire detection framework is proposed for the medium resolution multispectral remote sensing data. The basic part of this framework is a deep-learning-based semantic segmentation model specifically designed for active fire detection. A dataset is constructed with open-access Sentinel-2 imagery for the training and testing of the deep-learning model. The developed framework allows an automated Sentinel-2 data download, processing, and generation of the active fire detection results through time and location information provided by the user. Related performance is evaluated in terms of detection accuracy and processing efficiency. The last part of this thesis explored whether the coarse burned area products can be further improved through the synergy of multispectral, SAR, and InSAR features with higher spatial resolutions. A Siamese Self-Attention (SSA) classification is proposed for the multi-sensor burned area mapping and a multi-source dataset is constructed at the object level for the training and testing. Results are analyzed by different test sites, feature sources, and classification methods to assess the improvements achieved by the proposed method. All developed methods are validated with extensive processing of multi-source data acquired by Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR), Land, Vegetation, and Ice Sensor (LVIS), PolSARproSim+, Sentinel-1, and Sentinel-2. I hope these studies constitute a substantial contribution to the forest applications of multi-source remote sensing

Earth resources: A continuing bibliography with indexes (issue 62)

This bibliography lists 544 reports, articles, and other documents introduced into the NASA scientific and technical information system between April 1 and June 30, 1989. Emphasis is placed on the use of remote sensing and geophysical instrumentation in spacecraft and aircraft to survey and inventory natural resources and urban areas. Subject matter is grouped according to agriculture and forestry, environmental changes and cultural resources, geodesy and cartography, geology and mineral resources, hydrology and water management, data processing and distribution systems, instrumentation and sensors, and economic analysis

Relationship between Lidar-Derived Canopy Densities and the Scattering Phase Center of High-Resolution TanDEM-X Data

Abstract: The estimation of forestry parameters is essential to understanding the three-dimensional structure of forests. In this respect, the potential of X-band synthetic aperture radar (SAR) has been recognized for years. Many studies have been conducted on deriving tree heights with SAR data, but few have paid attention to the effects of the canopy structure. Canopy density plays an important role since it provides information about the vertical distribution of dominant scatterers in the forest. In this study, the position of the scattering phase center (SPC) of interferometric X-band SAR data is investigated with regard to the densest vegetation layer in a deciduous and coniferous forest in Germany by applying a canopy density index from high-resolution airborne laser scanning data. Two different methods defining the densest layer are introduced and compared with the position of the TanDEM-X SPC. The results indicate that the position of the SPC often coincides with the densest layer, with mean differences ranging from −1.6 m to +0.7 m in the deciduous forest and +1.9 m in the coniferous forest. Regarding relative tree heights, the SAR signal on average penetrates up to 15% (3.4 m) of the average tree height in the coniferous forest. In the deciduous forest, the difference increases to 18% (6.2 m) during summer and 24% (8.2 m) during winter. These findings highlight the importance of considering not only tree height but also canopy density when delineating SAR-based forest heights. The vertical structure of the canopy influences the position of the SPC, and incorporating canopy density can improve the accuracy of SAR-derived forest height estimations

Estimation of biophysical parameters in boreal forests from ERS and JERS SAR interferometry

The thesis describes investigations concerning the evaluation of ERS and JERS SAR images and repeat-pass interferometric SAR images for the retrieval of biophysical parameters in boreal forests. The availability of extensive data sets of images over several test sites located in Sweden, Finland and Siberia has allowed analysis of temporal dynamics of ERS and JERS backscatter and coherence, and of ERS interferometric phase. Modelling of backscatter, coherence and InSAR phase has been performed by means of the Water Cloud Model (WCM) and the Interferometric Water Cloud Model (IWCM); sensitivity analysis and implications for the retrieval of forest biophysical parameters have been thoroughly discussed. Model inversion has been carried out for stem volume retrieval using ERS coherence, ERS backscatter and JERS backscatter, whereas for tree height estimation the ERS interferometric phase has been used. Multi-temporal combination of ERS coherence images, and to a lesser extent of JERS backscatter images, can provide stem volume estimates comparable to stand-wise ground-based measurements. Since the information content of the interferometric phase is strongly degraded by phase noise and uncorrected atmospheric artefacts, the retrieved tree height shows large errors

Temporal changes in mediterranean pine forest biomass using synergy models of ALOS PALSAR-Sentinel 1-Landsat 8 Sensors

Currently, climate change requires the quantification of carbon stored in forest biomass. Synthetic aperture radar (SAR) data offers a significant advantage over other remote detection measurement methods in providing structural and biomass-related information about ecosystems. This study aimed to develop non-parametric Random Forest regression models to assess the changes in the aboveground forest biomass (AGB), basal area (G), and tree density (N) of Mediterranean pine forests by integrating ALOS-PALSAR, Sentinel 1, and Landsat 8 data. Variables selected from the Random Forest models were related to NDVI and optical textural variables. For 2015, the biomass models with the highest performance integrated ALS-ALOS2-Sentinel 1-Landsat 8 data (R2 = 0.59) by following the model using ALS data (R2 = 0.56), and ALOS2-Sentinel 1-Landsat 8 (R2 = 0.50). The validation set showed that R2 values vary from 0.55 (ALOS2-Sentinel 1-Landsat 8) to 0.60 (ALS-ALOS2-Sentinel 1-Landsat 8 model) with RMSE below 20 Mg ha−1. It is noteworthy that the individual Sentinel 1 (R2 = 0.49). and Landsat 8 (R2 = 0.47) models yielded equivalent results. For 2020, the AGB model ALOS2-Sentinel 1-Landsat 8 had a performance of R2 = 0.55 (validation R2 = 0.70) and a RMSE of 9.93 Mg ha−1. For the 2015 forest structural variables, Random Forest models, including ALOS PAL-SAR 2-Sentinel 1 Landsat 8 explained between 30% and 55% of the total variance, and for the 2020 models, they explained between 25% and 55%. Maps of the forests’ structural variables were generated for 2015 and 2020 to assess the changes during this period using the ALOS PALSAR 2-Sentinel 1-Landsat 8 model. Aboveground biomass (AGB), diameter at breast height (dbh), and dominant height (Ho) maps were consistent throughout the entire study area. However, the Random Forest models underestimated higher biomass levels (>100 Mg ha−1) and overestimated moderate biomass levels (30–45 Mg ha−1). The AGB change map showed values ranging from gains of 43.3 Mg ha−1 to losses of −68.8 Mg ha−1 during the study period. The integration of open-access satellite optical and SAR data can significantly enhance AGB estimates to achieve consistent and long-term monitoring of forest carbon dynamics

Comparing synthetic aperture radar and LiDAR for above-ground biomass estimation in Glen Affric, Scotland

Quantifying above-ground biomass (AGB) and carbon sequestration has been a significant focus of attention within the UNFCCC and Kyoto Protocol for improvement of national carbon accounting systems (IPCC, 2007; UNFCCC, 2011). A multitude of research has been carried out in relatively flat and homogeneous forests (Ranson & Sun, 1994; Beaudoin et al.,1994; Kurvonen et al., 1999; Austin et al., 2003; Dimitris et al., 2005), yet forests in the highlands, which generally form heterogeneous forest cover and sparse woodlands with mountainous terrain have been largely neglected in AGB studies (Cloude et al., 2001; 2008; Lumsdon et al., 2005; 2008; Erxue et al., 2009, Tan et al., 2010; 2011a; 2011b; 2011c; 2011d). Since mountain forests constitute approximately 28% of the total global forest area (Price and Butt, 2000), a better understanding of the slope effects is of primary importance in AGB estimation. The main objective of this research is to estimate AGB in the aforementioned forest in Glen Affric, Scotland using both SAR and LiDAR data. Two types of Synthetic Aperture Radar (SAR) data were used in this research: TerraSAR-X, operating at X-band and ALOS PALSAR, operating at L-band, both are fully polarimetric. The former data was acquired on 13 April 2010 and of the latter, two scenes were acquired on 17 April 2007 and 08 June 2009. Airborne LiDAR data were acquired on 09 June 2007. Two field measurement campaigns were carried out, one of which was done from winter 2006 to spring 2007 where physical parameters of trees in 170 circular plots were measured by the Forestry Commission team. Another intensive fieldwork was organised by myself with the help of my fellow colleagues and it comprised of tree measurement in two transects of 200m x 50m at a relatively flat and dense plantation forest and 400m x 50m at hilly and sparse semi-natural forest. AGB is estimated for both the transects to investigate the effectiveness of the proposed method at plot-level. This thesis evaluates the capability of polarimetric Synthetic Aperture Radar data for AGB estimation by investigating the relationship between the SAR backscattering coefficient and AGB and also the relationship between the decomposed scattering mechanisms and AGB. Due to the terrain and heterogeneous nature of the forests, the result from the backscatter-AGB analysis show that these forests present a challenge for simple AGB estimation. As an alternative, polarimetric techniques were applied to the problem by decomposing the backscattering information into scattering mechanisms based on the approach by Yamaguchi (2005; 2006), which are then regressed to the field measured AGB. Of the two data sets, ALOS PALSAR demonstrates a better estimation capacity for AGB estimation than TerraSAR-X. The AGB estimated results from SAR data are compared with AGB derived from LiDAR data. Since tree height is often correlated with AGB (Onge et al., 2008; Gang et al., 2010), the effectiveness of the tree height retrieval from LiDAR is evaluated as an indicator of AGB. Tree delineation was performed before AGB of individual trees were calculated allometrically. Results were validated by comparison to the fieldwork data. The amount of overestimation varies across the different canopy conditions. These results give some indication of when to use LiDAR or SAR to retrieve forest AGB. LiDAR is able to estimate AGB with good accuracy and the R2 value obtained is 0.97 with RMSE of 14.81 ton/ha. The R2 and RMSE obtained for TerraSAR-X are 0.41 and 28.5 ton/ha, respectively while for ALOS PALSAR data are 0.70 and 23.6 ton/ha, respectively. While airborne LiDAR data with very accurate height measurement and consequent three-dimensional (3D) stand profiles which allows investigation into the relationship between height, number density and AGB, it's limited to small coverage area, or large areas but at large cost. ALOS PALSAR, on the other hand, can cover big coverage area but it provide a lower resolution, hence, lower estimation accuracy