European Remote Sensing Research Articles

A global long-term, high-resolution satellite radar backscatter data record (1992–2022+): merging C-band ERS/ASCAT and Ku-band QSCAT

Abstract. Satellite radar backscatter contains unique information on land surface moisture, vegetation features, and surface roughness and has thus been used in a range of Earth science disciplines. However, there is no single global radar data set that has a relatively long wavelength and a decades-long time span. We here provide the first long-term (since 1992), high-resolution (∼8.9 km instead of the commonly used ∼25 km resolution) monthly satellite radar backscatter data set over global land areas, called the long-term, high-resolution scatterometer (LHScat) data set, by fusing signals from the European Remote Sensing satellite (ERS; 1992–2001; C-band; 5.3 GHz), Quick Scatterometer (QSCAT, 1999–2009; Ku-band; 13.4 GHz), and the Advanced SCATterometer (ASCAT; since 2007; C-band; 5.255 GHz). The 6-year data gap between C-band ERS and ASCAT was filled by modelling a substitute C-band signal during 1999–2009 from Ku-band QSCAT signals and climatic information. To this end, we first rescaled the signals from different sensors, pixel by pixel. We then corrected the monthly signal differences between the C-band and the scaled Ku-band signals by modelling the signal differences from climatic variables (i.e. monthly precipitation, skin temperature, and snow depth) using decision tree regression. The quality of the merged radar signal was assessed by computing the Pearson r, root mean square error (RMSE), and relative RMSE (rRMSE) between the C-band and the corrected Ku-band signals in the overlapping years (1999–2001 and 2007–2009). We obtained high Pearson r values and low RMSE values at both the regional (r≥0.92, RMSE ≤ 0.11 dB, and rRMSE ≤ 0.38) and pixel levels (median r across pixels ≥ 0.64, median RMSE ≤ 0.34 dB, and median rRMSE ≤ 0.88), suggesting high accuracy for the data-merging procedure. The merged radar signals were then validated against the European Space Agency (ESA) ERS-2 data, which provide observations for a subset of global pixels until 2011, even after the failure of on-board gyroscopes in 2001. We found highly concordant monthly dynamics between the merged radar signals and the ESA ERS-2 signals, with regional Pearson r values ranging from 0.79 to 0.98. These results showed that our merged radar data have a consistent C-band signal dynamic. The LHScat data set (https://doi.org/10.6084/m9.figshare.20407857; Tao et al., 2023) is expected to advance our understanding of the long-term changes in, e.g., global vegetation and soil moisture with a high spatial resolution. The data set will be updated on a regular basis to include the latest images acquired by ASCAT and to include even higher spatial and temporal resolutions.

Gangotri glacier dynamics from multi-sensor SAR and optical data

The present study has analyzed dynamics of Gangotri glacier using multiple remote sensing (RS) datasets and ground based observations. Interferometric Synthetic Aperture Radar (InSAR) data pairs from European Remote Sensing satellite (ERS 1/2) tandem pair for spring of 1996, Sentinel-1 SAR pairs and Japanese’s Advance Land Observation System (ALOS) PALSAR-2 SAR data for Spring of 2015 were used to derive glacier-surface velocity at seasonal time scale using Differential InSAR (DInSAR) techniques. Bi-static TanDEM-X (Experimental) data was used for the 1st time to estimate glacier surface elevation changes for a period of 22, 44, 88 days during summer of 2012 using InSAR techniques in this study. Annual glacier velocity was also estimated using temporal panchromatic data of LANDSAT-5 (30 m), LANDSAT-7/8 (15 m), Sentinel-2 (10 m) and Indian Remote Sensing Satellite IRS-1C/1D panchromatic (5 m) data during 1998–2019 with feature tracking approach. This study has estimated glacier surface velocity and surface elevation changes for the major parts of Gangotri glacier and its tributary glaciers using medium to high resolution optical and SAR datasets, at annual and seasonal time scale, which is an improvement over earlier studies, wherein snout based glacier recession or only main glacier velocities were reported. The velocity and slope were used to assess glacier-ice thickness distribution using Glabtop-2, slope dependent and laminar flow based methods over the Gangotri group of glaciers. The estimated ice thickness was estimated in the range of 58–550 m for the complete glacier while few small areas in middle & upper regions carry higher thickness of about 607 m. The estimated glacier-ice thickness was found in the range of 58–67 m at the snout region. The estimation was validated using 2014 field measurements from Terrestrial Laser Scanner (TLS) for the first time and correlation was found to be 0.799 at snout of the glacier.

Global land surface 250 m 8 d fraction of absorbed photosynthetically active radiation (FAPAR) product from 2000 to 2021

Abstract. The fraction of absorbed photosynthetically active radiation (FAPAR) is a critical land surface variable for carbon cycle modeling and ecological monitoring. Several global FAPAR products have been released and have become widely used; however, spatiotemporal inconsistency remains a large issue for the current products, and their spatial resolutions and accuracies can hardly meet the user requirements. An effective solution to improve the spatiotemporal continuity and accuracy of FAPAR products is to take better advantage of the temporal information in the satellite data using deep learning approaches. In this study, the latest version (V6) of the FAPAR product with a 250 m resolution was generated from Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance data and other information, as part of the Global LAnd Surface Satellite (GLASS) product suite. In addition, it was aggregated to multiple coarser resolutions (up to 0.25∘ and monthly). Three existing global FAPAR products (MODIS Collection 6; GLASS V5; and PRoject for On-Board Autonomy–Vegetation, PROBA-V, V1) were used to generate the time-series training samples, which were used to develop a bidirectional long short-term memory (Bi-LSTM) model. Direct validation using high-resolution FAPAR maps from the Validation of Land European Remote sensing Instrument (VALERI) and ImagineS networks revealed that the GLASS V6 FAPAR product has a higher accuracy than PROBA-V, MODIS, and GLASS V5, with an R2 value of 0.80 and root-mean-square errors (RMSEs) of 0.10–0.11 at the 250 m, 500 m, and 3 km scales, and a higher percentage (72 %) of retrievals for meeting the accuracy requirement of 0.1. Global spatial evaluation and temporal comparison at the AmeriFlux and National Ecological Observatory Network (NEON) sites revealed that the GLASS V6 FAPAR has a greater spatiotemporal continuity and reflects the variations in the vegetation better than the GLASS V5 FAPAR. The higher quality of the GLASS V6 FAPAR is attributed to the ability of the Bi-LSTM model, which involves high-quality training samples and combines the strengths of the existing FAPAR products, as well as the temporal and spectral information from the MODIS surface reflectance data and other information. The 250 m 8 d GLASS V6 FAPAR product for 2020 is freely available at https://doi.org/10.5281/zenodo.6405564 and https://doi.org/10.5281/zenodo.6430925 (Ma, 2022a, b) as well as at the University of Maryland for 2000–2021 (http://glass.umd.edu/FAPAR/MODIS/250m, last access 1 November 2022).

Mapping Samudra Tapu glacier: A holistic approach utilizing radar and optical remote sensing data for glacier radar facies mapping and velocity estimation

Himalayan glaciers have shown more sensitivity and visible changes to the climate change and global warming in the last 150 years. The highly rugged topography and inaccessible remote areas makes satellite images as the most appropriate source of information retrieval. We performed remote sensing based glacier change study for Samudra Tapu glacier, located in the Chandra basin of North-West Himalaya. In the present study, the capabilities of both optical and microwave remote sensing data was analysed in glacier change study in terms of its coverage, shift in equilibrium line altitude (ELA) and surface velocity over a period from 2000 to 2021. Multi Sensor (RISAT-1, Sentinel-1) time series of C-band SAR data along with a object oriented classification technique were used to identify different glacier facies such as percolation facies, icefalls, bare ice facies, refreeze snow and supraglacial debris. These classified maps were also used to detect the snow line and firn line along with ELA, aided with elevation information from digital elevation model (DEM). It was identified that more than 50 % of the total glacier area still lies into accumulation region. Further, we estimated the glacier surface velocity using Differential Interferometric Synthetic Aperture Radar (DInSAR) technique using European Remote Sensing Satellite (ERS-1/2) tandem data of 1996. High value of coherence was observed from the SAR return signal for one-day temporal difference. A mean velocity of 17–24 cm/day was found for the months of March and May 1996, highest flow rates were seen in the high accumulation area located in the Eastern and Southern Aspect of glacier. Spatial analysis of velocity patterns with respect to slope and aspect show that high rates of flow was found in southern slopes and movement rates generally increase with increase in slope. Feature tracking approach was used to estimate the glacier flow for long term and seasonal basis using optical and SAR datasets (IRS-1C, 1D PAN, Landsat-7, 8 PAN, and TANDEM-x) during 1999–2020 period. The results suggest that glacier flow varies with season, i.e., high velocity during spring-summer season, as compared to late summer or winter and, the rate of ice flow changes over the years. The mean glacier velocity reduced to 49.5 m/year during 2013–2020 time, as compared to 67.67 m/year during 1999–2003 time. These results of reducing glacier velocity and changing snow line altitude indicates enhanced glacier’s melt rate and overall negative mass balance for Smudra tapu glacier.

Evaluation of Global Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) Products at 500 m Spatial Resolution

The fraction of absorbed photosynthetically active radiation (FAPAR) is a key biophysical variable directly associated with the photosynthetic activity of plants. Several global FAPAR products with different spatial resolutions have been generated from remote sensing data, and much work has focused on validating them. However, those studies have primarily evaluated global FAPAR products at a spatial resolution of 1 km or more, whereas few studies have evaluated the global 500 m resolution FAPAR product distributed in recent years. Furthermore, there are a few FAPAR products, including black-sky, white-sky and blue-sky FAPAR datasets, and almost no studies have evaluated these products. In this article, three global FAPAR products at 500 m resolution, namely MODIS (only black-sky FAPAR), MUSES and EBR (black-sky, white-sky and blue-sky FAPAR) were compared to evaluate their temporal and spatial discrepancies and direct validation was conducted to compare these FAPAR products with the FAPAR values derived from the high-resolution reference maps from the Validation of Land European Remote Sensing Instrument (VALERI) and Implementing Multi-Scale Agricultural Indicators Exploiting Sentinels (IMAGINES) projects. The results showed that the MUSES FAPAR product exhibited the best spatial integrity, whereas the MODIS and EBR FAPAR products had many missing pixels in the equatorial rainforest regions and at high latitudes in the Northern Hemisphere. The MODIS, MUSES and EBR FAPAR products were generally consistent in their spatial patterns. However, a relatively large discrepancy among these FAPAR products was present in the equatorial rainforest regions and the middle and high latitude regions where the main vegetation type was forest. The differences between the black-sky and white-sky FAPAR datasets at the global scale were significant. In January, the MUSES and EBR black-sky FAPAR values were larger than their white-sky FAPAR values in the region north of 30° north latitude but they were smaller than their white-sky FAPAR values in the region south of 30° north latitude. In July, the MUSES and EBR black-sky FAPAR values were lower than their white-sky FAPAR values in the region north of 30° south latitude and they were larger than their white-sky FAPAR values in the region south of 30° south latitude. The temporal profiles of the MUSES FAPAR product were continuous and smooth, whereas those of the MODIS and EBR FAPAR products showed many fluctuations, particularly during the growing seasons. Direct validation indicated that the MUSES FAPAR product had the best accuracy (R2 = 0.6932, RMSE = 0.1495) compared to the MODIS FAPAR product (R2 = 0.6202, RMSE = 0.1710) and the EBR FAPAR product (R2 = 0.5746, RMSE = 0.1912).

A combined use of in situ and satellite-derived observations to characterize surface hydrology and its variability in the Congo River basin

Abstract. The Congo River basin (CRB) is the second largest river system in the world, but its hydroclimatic characteristics remain relatively poorly known. Here, we jointly analyse a large record of in situ and satellite-derived observations, including a long-term time series of surface water height (SWH) from radar altimetry (a total of 2311 virtual stations) and surface water extent (SWE) from a multi-satellite technique, to characterize the CRB surface hydrology and its variability. First, we show that SWH from altimetry multi-missions agrees well with in situ water stage at various locations, with the root mean square deviation varying from 10 cm (with Sentinel-3A) to 75 cm (with European Remote Sensing satellite-2). SWE variability from multi-satellite observations also shows a plausible behaviour over a ∼25-year period when evaluated against in situ observations from the subbasin to basin scale. Both datasets help to better characterize the large spatial and temporal variability in hydrological patterns across the basin, with SWH exhibiting an annual amplitude of more than 5 m in the northern subbasins, while the Congo River main stream and Cuvette Centrale tributaries vary in smaller proportions (1.5 to 4.5 m). Furthermore, SWH and SWE help illustrate the spatial distribution and different timings of the CRB annual flood dynamic and how each subbasin and tributary contribute to the hydrological regime at the outlet of the basin (the Brazzaville/Kinshasa station), including its peculiar bimodal pattern. Across the basin, we estimate the time lag and water travel time to reach the Brazzaville/Kinshasa station to range from 0–1 month in its vicinity in downstream parts of the basin and up to 3 months in remote areas and small tributaries. Northern subbasins and the central Congo region contribute highly to the large peak in December–January, while the southern part of the basin supplies water to both hydrological peaks, in particular to the moderate one in April–May. The results are supported using in situ observations at several locations in the basin. Our results contribute to a better characterization of the hydrological variability in the CRB and represent an unprecedented source of information for hydrological modelling and to study hydrological processes over the region.

A 30-year monthly 5 km gridded surface elevation time series for the Greenland Ice Sheet from multiple satellite radar altimeters

Abstract. A long-term time series of ice sheet surface elevation change (SEC) is an essential parameter to assess the impact of climate change. In this study, we used an updated plane-fitting least-squares regression strategy to generate a 30-year surface elevation time series for the Greenland Ice Sheet (GrIS) at monthly temporal resolution and 5×5 km grid spatial resolution using ERS-1 (European Remote Sensing), ERS-2, Envisat, and CryoSat-2 satellite radar altimeter observations obtained between August 1991 and December 2020. The ingenious corrections for intermission bias were applied using an updated plane-fitting least-squares regression strategy. Empirical orthogonal function (EOF) reconstruction was used to supplement the sparse monthly gridded data attributable to poor observations in the early years. Validation using both airborne laser altimeter observations and the European Space Agency GrIS Climate Change Initiative (CCI) product indicated that our merged surface elevation time series is reliable. The accuracy and dispersion of errors of SECs of our results were 19.3 % and 8.9 % higher, respectively, than those of CCI SECs and even 30.9 % and 19.0 % higher, respectively, in periods from 2006–2010 to 2010–2014. Further analysis showed that our merged time series could provide detailed insight into GrIS SEC on multiple temporal (up to 30 years) and spatial scales, thereby providing an opportunity to explore potential associations between ice sheet change and climatic forcing. The merged surface elevation time series data are available at https://doi.org/10.11888/Glacio.tpdc.271658 (Zhang et al., 2021).

A 250 m resolution global leaf area index product derived from MODIS surface reflectance data

ABSTRACT There are several global leaf area index (LAI) products currently available. The spatial resolution of these products is 500 m and above, which is unsuitable for many applications requiring higher spatial resolution. In the past several years, we developed a method to estimate the LAI from time series satellite remote sensing data using general regression neural networks. The method has been used to generate global LAI products at 500 m and 1000 m from Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance data, and a global LAI product at 0.05° from Advanced Very High Resolution Radiometer (AVHRR) surface reflectance data. In this study, the method was extended to generate a global LAI product at 250 m (one of the MUltiscale Satellite remotE Sensing (MUSES) product suite) from MODIS surface reflectance data in the red and near-infrared (NIR) bands. As far as we know, it is the first global LAI product at 250 m spatial resolution and is the highest spatial resolution global LAI product available. The spatial and temporal consistency of the MUSES LAI product was evaluated by comparing it with the MODIS LAI product, and the MUSES LAI product was validated by high-resolution reference maps at the Validation of Land European Remote Sensing Instruments (VALERI) and Implementing Multi-Scale Agricultural Indicators Exploiting Sentinels (IMAGINES) sites representative of different biomes. The root mean square error (RMSE) of the MUSES LAI product versus the LAI values derived from the high-resolution reference maps over the VALERI and IMAGINES sites was 0.9984, and the bias of the MUSES LAI product was −0.2005.

Evaluating the Reliability of Air Temperature From ERA5 Reanalysis Data

The reliability of European Remote Sensing 5 (ERA5) satellite-based air temperature data is under investigation in this letter. To evaluate this, the ERA5 data will be compared with land-based data obtained from weather stations on the global historical climatology network (GHCN). Two climate regions are taken into consideration, temperate and tropical. Five years worth of data is collected and compared through box plots, regression models, and statistical metrics. The results show that the satellite temperature performs better in the temperate region than the tropical region. This suggests that the time of year and climate region have an impact on the accuracy of the satellite data as milder temperatures produce better approximations.

Neural Network-Based Urban Change Monitoring with Deep-Temporal Multispectral and SAR Remote Sensing Data

Remote-sensing-driven urban change detection has been studied in many ways for decades for a wide field of applications, such as understanding socio-economic impacts, identifying new settlements, or analyzing trends of urban sprawl. Such kinds of analyses are usually carried out manually by selecting high-quality samples that binds them to small-scale scenarios, either temporarily limited or with low spatial or temporal resolution. We propose a fully automated method that uses a large amount of available remote sensing observations for a selected period without the need to manually select samples. This enables continuous urban monitoring in a fully automated process. Furthermore, we combine multispectral optical and synthetic aperture radar (SAR) data from two eras as two mission pairs with synthetic labeling to train a neural network for detecting urban changes and activities. As pairs, we consider European Remote Sensing (ERS-1/2) and Landsat 5 Thematic Mapper (TM) for 1991–2011 and Sentinel 1 and 2 for 2017–2021. For every era, we use three different urban sites—Limassol, Rotterdam, and Liège—with at least 500km2 each, and deep observation time series with hundreds and up to over a thousand of samples. These sites were selected to represent different challenges in training a common neural network due to atmospheric effects, different geographies, and observation coverage. We train one model for each of the two eras using synthetic but noisy labels, which are created automatically by combining state-of-the-art methods, without the availability of existing ground truth data. To combine the benefit of both remote sensing types, the network models are ensembles of optical- and SAR-specialized sub-networks. We study the sensitivity of urban and impervious changes and the contribution of optical and SAR data to the overall solution. Our implementation and trained models are available publicly to enable others to utilize fully automated continuous urban monitoring.

Evaluation of the version 5.0 global land surface satellite (GLASS) leaf area index product derived from MODIS data

ABSTRACT The Global Land Surface Satellite (GLASS) leaf area index (LAI) product is one of the most widely used global LAI products in the scientific community. The latest (version 5) GLASS LAI product has been generated from Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance data. The purpose of this paper is to evaluate the quality of the version 5 GLASS LAI product. The GLASS LAI product was compared with the latest MODIS LAI product (MCD15A2 H, Collection 6) and the second version of Geoland2 (GEOV2) LAI product to evaluate their temporal and spatial discrepancies. A direct validation was conducted to compare these LAI products to the LAI values derived from the high-resolution reference maps from the Validation of Land European Remote Sensing Instruments (VALERI) and Implementing Multi-Scale Agricultural Indicators Exploiting Sentinels (IMAGINES) sites. The results show that the GLASS and GEOV2 LAI products have great spatial integrity. However, the MODIS LAI product contains many missing pixels in tropical areas. These LAI products follow fairly consistent seasonal characteristics. The spatial discrepancies of these LAI products mainly exist in forest areas, especially evergreen broadleaf forests where the GLASS LAI values are generally lower than the GEOV2 LAI values by approximately 1.0 LAI units and lower than the MODIS LAI values by 0.5 to 1.0 LAI units. The spatial distribution of these LAI products has slight discrepancies in savannahs, broadleaf crops, grasses/cereal crops and shrubs. The GLASS and GEOV2 LAI products capture a complete and reasonable temporal profile, contrasting with the MODIS LAI product, which shows dramatic fluctuations, particularly during the growing seasons. These LAI products show similar temporal trajectories and interannual variations for all biome types except evergreen broadleaf forests. The direct validation shows that the accuracy of the GLASS LAI product is better than the accuracy of the MODIS and GEOV2 LAI products. The coefficient of determination (R2) of the GLASS, MODIS and GEOV2 LAI products versus the LAI values derived from the high-resolution reference maps are 0.68, 0.47 and 0.55, respectively, and the root mean square error (RMSE) of these products are 0.86, 1.22 and 1.21, respectively.

ERS-1/2 and Sentinel-1 SAR Data Mining for Flood Hazard and Risk Assessment in Lima, Peru

The coastline environment and urban areas of Peru overlooking the Pacific Ocean are among the most affected by El Niño-Southern Oscillation (ENSO) events, and its cascading hazards such as floods, landslides and avalanches. In this work, the complete archives of the European Space Agency (ESA)’s European Remote-Sensing (ERS-1/2) missions and European Commission’s Copernicus Sentinel-1 constellation were screened to select synthetic aperture radar (SAR) images covering the most severe and recent ENSO-related flooding events that affected Lima, the capital and largest city of Peru, in 1997–1998 and 2017–2018. Based on SAR backscatter color composites and ratio maps retrieved from a series of pre-, cross- and post-event SAR pairs, flooded areas were delineated within the Rímac River watershed. These are mostly concentrated along the riverbanks and plain, where low-lying topography and gentle slopes (≤5°), together with the presence of alluvial deposits, also indicate greater susceptibility to flooding. A total of 409 areas (58.50 km2) revealing change were mapped, including 197 changes (32.10 km2) due to flooding-related backscatter variations (flooded areas, increased water flow in the riverbed, and riverbank collapses and damage), and 212 (26.40 km2) due to other processes (e.g., new urban developments, construction of river embankments, other engineering works, vegetation changes). Urban and landscape changes potentially contributing, either detrimentally or beneficially, to flooding susceptibility were identified and considered in the overall assessment of risk. The extent of built-up areas within the basin was mapped by combining information from the 2011 Global Urban Footprint (GUF) produced by the German Aerospace Center (DLR), the Open Street Map (OSM) accessed from the Quantum GIS (QGIS) service, and 2011–2019 very high-resolution optical imagery from Google Earth. The resulting flooding risk map highlights the sectors of potential concern along the Rímac River, should flooding events of equal severity as those captured by SAR images occur in the future.

Cloud<i>_</i>cci ATSR-2 and AATSR data set version 3: a 17-year climatology of global cloud and radiation properties

Abstract. We present version 3 (V3) of the Cloud_cci Along-Track Scanning Radiometer (ATSR) and Advanced ATSR (AATSR) data set. The data set was created for the European Space Agency (ESA) Cloud_cci (Climate Change Initiative) programme. The cloud properties were retrieved from the second ATSR (ATSR-2) on board the second European Remote Sensing Satellite (ERS-2) spanning 1995–2003 and the AATSR on board Envisat, which spanned 2002–2012. The data are comprised of a comprehensive set of cloud properties: cloud top height, temperature, pressure, spectral albedo, cloud effective emissivity, effective radius, and optical thickness, alongside derived liquid and ice water path. Each retrieval is provided with its associated uncertainty. The cloud property retrievals are accompanied by high-resolution top- and bottom-of-atmosphere shortwave and longwave fluxes that have been derived from the retrieved cloud properties using a radiative transfer model. The fluxes were generated for all-sky and clear-sky conditions. V3 differs from the previous version 2 (V2) through development of the retrieval algorithm and attention to the consistency between the ATSR-2 and AATSR instruments. The cloud properties show improved accuracy in validation and better consistency between the two instruments, as demonstrated by a comparison of cloud mask and cloud height with co-located CALIPSO data. The cloud masking has improved significantly, particularly in its ability to detect clear pixels. The Kuiper Skill score has increased from 0.49 to 0.66. The cloud top height accuracy is relatively unchanged. The AATSR liquid water path was compared with the Multisensor Advanced Climatology of Liquid Water Path (MAC-LWP) in regions of stratocumulus cloud and shown to have very good agreement and improved consistency between ATSR-2 and AATSR instruments. The correlation with MAC-LWP increased from 0.4 to over 0.8 for these cloud regions. The flux products are compared with NASA Clouds and the Earth's Radiant Energy System (CERES) data, showing good agreement within the uncertainty. The new data set is well suited to a wide range of climate applications, such as comparison with climate models, investigation of trends in cloud properties, understanding aerosol–cloud interactions, and providing contextual information for co-located ATSR-2/AATSR surface temperature and aerosol products. The following new digital identifier has been issued for the Cloud_cci ATSR-2/AATSRv3 data set: https://doi.org/10.5676/DWD/ESA_Cloud_cci/ATSR2-AATSR/V003 (Poulsen et al., 2019).

A Scaling-Based Method for the Rapid Retrieval of FPAR From Fine-Resolution Satellite Data in the Remote-Sensing Trend-Surface Framework

Accurate estimation of the fine-resolution fraction of absorbed photosynthetically active radiation (FPAR) across broad spatial extents and long time periods requires efficient and applicable methods. The existing methods can hardly provide a balance between accuracy, simplicity, and transferability through space and time. Within the remote-sensing trend-surface conceptual framework, this article proposes a scaling-based method to efficiently retrieve FPAR from fine-resolution satellite data using coarse-resolution FPAR products as a reference. The method was particularly developed and applied to Moderate Resolution Imaging Spectroradiometer (MODIS) FPAR product and Landsat imagery. First, necessary prior knowledge related to FPAR retrieval and scaling theories was used to explicitly linearize the complex relationship between MODIS FPAR and Landsat surface reflectance. Second, the explicit linear model for FPAR estimation was trained through one-pair image learning for each date to estimate FPAR from Landsat imagery in real time. Both homogeneous and heterogeneous cases were considered. The method was validated at ten selected worldwide sites from the Validation of Land European Remote Sensing Instruments (VALERI) program and derived an overall root mean squared error (RMSE) of 0.133. A long time series of FPAR data set at the 30-m resolution was generated at the regional scale (approximately 2000 km <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sup> ) for 13 years (2000–2012). The results were accurate (RMSE = 0.072) and MODIS-consistent, which were significantly better than those of the normalized difference vegetation index (NDVI) downscaling-based and regression tree methods. The scaling-based method provides accurate, MODIS-consistent and spatially consistent FPAR estimates in real time, is highly transferrable through space and time, and allows for future extension of FPAR estimates to the era of the Landsat series satellites.

Gravity Anomaly in Kelud, Kasinan-Songgoriti, and Arjuno-Welirang Volcano Hosted Geothermal Area, East Java, Indonesia

This study aims to determine the relationship of heat reservoirs in the Kelud, Kasinan-Songgoriti, and Arjuno-Welirang geothermal systems based on gravity data analysis. Gravity data are obtained from Geodetic Satellite (GEOSAT) and European Remote Sensing-1 (ERS-1) Satellite which have been corrected to free air correction. The result of gravity data analysis is in the form of a complete Bouguer anomaly which represents the gravity anomaly below the surface. The results of the complete Bouguer anomaly value obtained were -15,238 mGal to 86,087 mGal. Based on these results, regional anomalies and residual anomalies will be separated to determine the depth of the two anomalies. 3D modeling was carried out based on the complete Bouguer anomaly data to determine the reservoir relationships in the Kelud, Kasinan-Songgoriti, and Arjuno-Welirang geothermal systems.

Transformative Urban Changes of Beijing in the Decade of the 2000s

The rapid economic growth, the exodus from rural to urban areas, and the associated extreme urban development that occurred in China in the decade of the 2000s have severely impacted the environment in Beijing, its vicinity, and beyond. This article presents an innovative approach for assessing mega-urban changes and their impact on the environment based on the use of decadal QuikSCAT (QSCAT) satellite data, acquired globally by the SeaWinds scatterometer over that period. The Dense Sampling Method (DSM) is applied to QSCAT data to obtain reliable annual infrastructure-based urban observations at a posting of ~1 km. The DSM-QSCAT data, along with different DSM-based change indices, were used to delineate the extent of the Beijing infrastructure-based urban area in each year between 2000 and 2009, and assess its development over time, enabling a physical quantification of its urbanization which reflects the implementation of various development policies during the same time period. Eventually, as a proxy for the impact of Beijing urbanization on the environment, the decadal trend of its infrastructure-based urbanization is compared with that of the corresponding tropospheric nitrogen dioxide (NO2) column densities as observed from the Global Ozone Monitoring Experiment (GOME) instrument aboard the second European Remote Sensing satellite (ERS-2) between 2000 and 2002, and from the SCanning Imaging Absorption SpectroMeter for Atmospheric CHartographY aboard of the ESA’s ENVIronmental SATellite (SCIAMACHY /ENVISAT) between 2003 and 2009. Results reveal a threefold increase of the yearly tropospheric NO2 column density within the Beijing infrastructure-based urban area extent in 2009, which had quadrupled since 2000.

Impact of high-resolution sea surface temperature representation on the forecast of small Mediterranean catchments' hydrological responses to heavy precipitation

Abstract. Operational meteo-hydrological forecasting chains are affected by many sources of uncertainty. In coastal areas characterized by complex topography, with several medium-to-small size catchments, quantitative precipitation forecast becomes even more challenging due to the interaction of intense air–sea exchanges with coastal orography. For such areas, which are quite common in the Mediterranean Basin, improved representation of sea surface temperature (SST) space–time patterns can be particularly important. The paper focuses on the relative impact of different resolutions of SST representation on regional operational forecasting chains (up to river discharge estimates) over coastal Mediterranean catchments, with respect to two other fundamental options while setting up the system, i.e. the choice of the forcing general circulation model (GCM) and the possible use of a three-dimensional variational assimilation (3D-Var) scheme. Two different kinds of severe hydro-meteorological events that affected the Calabria region (southern Italy) in 2015 are analysed using the WRF-Hydro atmosphere–hydrology modelling system in its uncoupled version. Both of the events are modelled using the 0.25∘ resolution global forecasting system (GFS) and the 16 km resolution integrated forecasting system (IFS) initial and lateral atmospheric boundary conditions, which are from the European Centre for Medium-Range Weather Forecasts (ECMWF), applying the WRF mesoscale model for the dynamical downscaling. For the IFS-driven forecasts, the effects of the 3D-Var scheme are also analysed. Finally, native initial and lower boundary SST data are replaced with data from the Medspiration project by Institut Français de Recherche pour L'Exploitation de la Mer (IFREMER)/Centre European Remote Sensing d'Archivage et de Traitement (CERSAT), which have a 24 h time resolution and a 2.2 km spatial resolution. Precipitation estimates are compared with both ground-based and radar data, as well as discharge estimates with stream gauging stations' data. Overall, the experiments highlight that the added value of high-resolution SST representation can be hidden by other more relevant sources of uncertainty, especially the choice of the general circulation model providing the boundary conditions. Nevertheless, in most cases, high-resolution SST fields show a non-negligible impact on the simulation of the atmospheric boundary layer processes, modifying flow dynamics and/or the amount of precipitated water; thus, this emphasizes the fact that uncertainty in SST representation should be duly taken into account in operational forecasting in coastal areas.

Decadal‐Scale Aquifer Dynamics and Structural Complexities at a Municipal Wellfield Revealed by 25 Years of InSAR and Recent Groundwater Temperature Observations

AbstractOver the past 35 years the Buckman wellfield near Santa Fe, New Mexico, experienced production well drawdowns in excess of 180 m, resulting in ground subsidence and surface cracks. Increased reliance on surface water diversions since 2011 has reduced pumping and yielded water level recovery. To characterize the impact of wellfield management decisions on the aquifer system, we reconstruct the surface deformation history through the European Remote Sensing Satellite, Advanced Land Observing Satellite, and Sentinel‐1 Interferometric Synthetic Aperture Radar (InSAR) time series analysis during episodes of drawdown (1993–2000), recovery (2007–2010), and modern management (2015–2018) in discontinuous observations over a 25‐year period. The observed deformation generally reflects changes in hydraulic head. However, at times during the wellfield recovery, the deformation signal is complex, with patterns of uplift and subsidence suggesting a compartmentalized aquifer system. Recent records of locally high geothermal gradients and an overall warming of the system (~0.5°C during the water level recovery) obtained from repeat temperature measurements between 2013 and 2018 constrain a conceptual model of convective heat transfer that requires a vertical permeable zone near an observed fault. To reproduce observed temperature patterns at monitoring wells, high basal heat flow and convective cooling associated with downward flow of water from cool shallow aquifers during the drawdown period is necessary. The fault, however, appears to die out southward or may be locally permeable, as conceptual cross‐sectional hydrologic modeling reproduces the surface deformation without such a structure. Our work demonstrates the importance of incorporating well‐constrained stratigraphy and structure when modeling near‐surface deformation induced by, for instance, groundwater production.

Spatiotemporal variability and decadal trends of snowmelt processes on Antarctic sea ice observed by satellite scatterometers

Abstract. The timing and intensity of snowmelt processes on sea ice are key drivers determining the seasonal sea-ice energy and mass budgets. In the Arctic, satellite passive microwave and radar observations have revealed a trend towards an earlier snowmelt onset during the last decades, which is an important aspect of Arctic amplification and sea-ice decline. Around Antarctica, snowmelt on perennial ice is weak and very different than in the Arctic, with most snow surviving the summer. Here we compile time series of snowmelt onset dates on seasonal and perennial Antarctic sea ice from 1992 to 2014/15 using active microwave observations from the European Space Agency's (ESA) European Remote Sensing (ERS) 1 and 2 missions (ERS-1 and ERS-2), Quick Scatterometer (QSCAT), and Advanced Scatterometer (ASCAT) radar scatterometers. We define two snowmelt transition stages: a weak backscatter rise, indicating the initial warming and destructive metamorphism of the snowpack (pre-melt), followed by a rapid backscatter rise, indicating the onset of thaw–freeze cycles (snowmelt). Results show large interannual variability, with an average pre-melt onset date of 29 November and melt onset of 10 December, respectively, on perennial ice, without any significant trends over the study period, consistent with the small trends of Antarctic sea-ice extent. There was a latitudinal gradient from early snowmelt onsets in mid-November in the northern Weddell Sea to late (end of December) or even absent snowmelt conditions in the southern Weddell Sea. We show that QSCAT Ku-band-derived (13.4 GHz signal frequency) pre-melt and snowmelt onset dates are earlier by 20 and 18 d, respectively, than ERS and ASCAT C-band-derived (5.6 GHz) dates. This offset has been considered when constructing the time series. Snowmelt onset dates from passive microwave observations (37 GHz) are later by 14 and 6 d than those from the scatterometers, respectively. Based on these characteristic differences between melt onset dates observed by different microwave wavelengths, we developed a conceptual model which illustrates how the seasonal evolution of snow temperature profiles may affect different microwave bands with different penetration depths. These suggest that future multi-frequency active and passive microwave satellite missions could be used to resolve melt processes throughout the vertical snow column of thick snow on perennial Antarctic sea ice.

Arctic Ocean Sea Level Record from the Complete Radar Altimetry Era: 1991–2018

In recent years, there has been a large focus on the Arctic due to the rapid changes of the region. Arctic sea level determination is challenging due to the seasonal to permanent sea-ice cover, lack of regional coverage of satellites, satellite instruments ability to measure ice, insufficient geophysical models, residual orbit errors, challenging retracking of satellite altimeter data. We present the European Space Agency (ESA) Climate Change Initiative (CCI) Technical University of Denmark (DTU)/Technischen Universität München (TUM) sea level anomaly (SLA) record based on radar satellite altimetry data in the Arctic Ocean from the European Remote Sensing satellite number 1 (ERS-1) (1991) to CryoSat-2 (2018). We use updated geophysical corrections and a combination of altimeter data: Reprocessing of Altimeter Product for ERS (REAPER) (ERS-1), ALES+ retracker (ERS-2, Envisat), combination of Radar Altimetry Database System (RADS) and DTUs in-house retracker LARS (CryoSat-2). Furthermore, this study focuses on the transition between conventional and Synthetic Aperture Radar (SAR) altimeter data to make a smooth time series regarding the measurement method. We find a sea level rise of 1.54 mm/year from September 1991 to September 2018 with a 95% confidence interval from 1.16 to 1.81 mm/year. ERS-1 data is troublesome and when ignoring this satellite the SLA trend becomes 2.22 mm/year with a 95% confidence interval within 1.67–2.54 mm/year. Evaluating the SLA trends in 5 year intervals show a clear steepening of the SLA trend around 2004. The sea level anomaly record is validated against tide gauges and show good results. Additionally, the time series is split and evaluated in space and time.