Miracle-dz/UniFAR-FTUs huggingface dataset:
What Is the Spatial Resolution of grace Satellite Products for Hydrology?

The mass change information from the Gravity Recovery And Climate Experiment (grace) satellite mission is available in terms of noisy spherical harmonic coefficients truncated at a maximum degree (band-limited). Therefore, filtering is an inevitable step in post-processing of grace fields to extract meaningful information about mass redistribution in the Earth-system. It is well known from previous studies that a number can be allotted to the spatial resolution of a band-limited spherical harmonic spectrum and also to a filtered field. Furthermore, it is now a common practice to correct the filtered grace data for signal damage due to filtering (or convolution in the spatial domain). These correction methods resemble deconvolution, and, therefore, the spatial resolution of the corrected grace data have to be reconsidered. Therefore, the effective spatial resolution at which we can obtain mass changes from grace products is an area of debate. In this contribution, we assess the spatial resolution both theoretically and practically. We confirm that, theoretically, the smallest resolvable catchment is directly related to the band-limit of the spherical harmonic spectrum of the grace data. However, due to the approximate nature of the correction schemes and the noise present in grace data, practically, the complete band-limited signal cannot be retrieved. In this context, we perform a closed-loop simulation comparing four popular correction schemes over 255 catchments to demarcate the minimum size of the catchment whose signal can be efficiently recovered by the correction schemes. We show that the amount of closure error is inversely related to the size of the catchment area. We use this trade-off between the error and the catchment size for defining the potential spatial resolution of the grace product obtained from a correction method. The magnitude of the error and hence the spatial resolution are both dependent on the correction scheme. Currently, a catchment of the size ≈63,000 km 2 can be resolved at an error level of 2 cm in terms of equivalent water height.
Detection of Spatiotemporal Extreme Changes in Atmospheric CO2 Concentration Based on Satellite Observations

Atmospheric CO2 concentrations are sensitive to the effects of climate extremes on carbon sources and sinks of the land biosphere. Therefore, extreme changes of atmospheric CO2 can be used to identify anomalous sources and sinks of carbon. In this study, we develop a spatiotemporal extreme change detection method for atmospheric CO2 concentrations using column-averaged CO2 dry air mole fraction (XCO2) retrieved from the Greenhouse gases Observing SATellite (GOSAT) from 1 June 2009 to 31 May 2016. For extreme events identified, we attributed the main drivers using surface environmental parameters, including surface skin temperature, self-calibrating Palmer drought severity index, burned area, and gross primary production (GPP). We also tested the sensitivity of XCO2 response to changing surface CO2 fluxes using model simulations and Goddard Earth Observing System (GEOS)-Chem atmospheric transport. Several extreme high XCO2 events are detected around mid-2010 over Eurasia and in early 2016 in the tropics. The magnitudes of extreme XCO2 increases are around 1.5–1.8 ppm in the Northern Hemisphere and 1.2–1.4 ppm in Southern Hemisphere. The spatiotemporal pattern of detected high XCO2 events are similar to patterns of local surface environmental parameter extremes. The extreme high XCO2 events often occurred during periods of increased temperature, severe drought, increased wildfire or reduced GPP. Our sensitivity tests show that the magnitude of detectable anomalies varies with location, for example 25% or larger anomalies in local CO2 emission fluxes are detectable in tropical forest, whereas anomalies must be half again as large in mid-latitudes (~37.5%). In conclusion, we present a method for extreme high XCO2 detection, and large changes in land CO2 fluxes. This provides another tool to monitor large-scale changes in the land carbon sink and potential feedbacks on the climate system.
Groundwater Storage Changes in China from Satellite Gravity: An Overview

Groundwater plays a critical role in the global water cycle and is the drinking source for almost half of the world’s population. However, exact quantification of its storage change remains elusive due primarily to limited ground observations in space and time. The Gravity Recovery and Climate Experiment (GRACE) twin-satellite data have provided global observations of water storage variations at monthly sampling for over a decade and a half, and is enable to estimate changes in groundwater storage (GWS) after removing other water storage components using auxiliary datasets and models. In this paper, we present an overview of GWS changes in three main aquifers within China using GRACE data, and conduct a comprehensive accuracy assessment using in situ ground well observations and hydrological models. GRACE detects a significant GWS depletion rate of 7.2 ± 1.1 km3/yr in the North China Plain (NCP) during 2002–2014, consistent with ground well observations and model predictions. The Liaohe River Basin (LRB) experienced a pronounced GWS decline during 2005–2009, at a depletion rate of 5.0 ± 1.2 km3/yr. Since 2010, GRACE-based GWS reveal a slow recovery in the LRB, with excellent agreement with ground well observations. For the whole study period 2002–2014, no significant long-term GWS depletion is found in the LRB nor in the Tarim Basin. A case study in the Inner Tibetan Plateau highlights there still exist large uncertainties in GRACE-based GWS change estimates.
Characteristic and Driving Factors of Aerosol Optical Depth over Mainland China during 1980-2017

Since the reform and opening up of China, the increasing aerosol emissions have posted great challenges to the country’s climate change and human health. The aerosol optical depth (AOD) is one of the main physical indicators quantifying the atmospheric turbidity and air pollution. In this study, 38-years (1980–2017) of spatial and temporal variations of AOD in China were analyzed using AOD records derived from MODIS atmosphere products and the MERRA-2 dataset. The results showed that the annual mean AOD values throughout China have gone through an increasing, but fluctuating, trend, especially in 1982 and in 1992 due to two volcano eruptions; the AOD values experienced a dramatically increasing period during 2000–2007 with the rapid economic development and “population explosions” in China/after 2008, the AOD values gradually decreased from 0.297 (2008) to 0.257 (2017). The AOD values in China were generally higher in spring than that in other seasons. The Sichuan Basin has always been an area with high AOD values owing to the strong human activity and the basin topography (hindering aerosol diffusions in the air). In contrast, the Qinghai Tibet Plateau has always been an area with low AOD values due to low aerosol emissions and clear sky conditions there. The trend analysis of AOD values during 1980–2017 in China indicated that the significant increasing trend was mainly observed in Southeastern China. By contrast, the AOD values in the northernmost of China showed a significant decreasing trend. Then, the contributions (AODP) of the AOD for black carbon aerosol (BCAOD), dust aerosol (DUAOD), organic carbon aerosol (OCAOD), sea salt aerosol (SSAOD), and SO4 aerosol (SO4AOD) to the total AOD values were calculated. The results showed that DUAOD (25.43%) and SO4AOD (49.51%) were found to be the main driving factors for the spatial and temporal variations of AOD values. Finally, the effects of anthropogenic aerosol emissions, socioeconomic factors, and land-use and land coverage changes on AOD were analyzed. The GDP, population density, and passenger traffic volume were found to be the main socioeconomic drivers for AOD distributions. Relatively larger AOD values were mainly found in urban land and land covered by water, while lower AOD values were found in grassland and permanent glacier areas. Remote Sens. 2018, 10, 1064; doi:10.3390/rs10071064 www.mdpi. com/journal/remotesensing Remote Sens. 2018, 10, 1064 2 of 25
The Characteristics of the Aerosol Optical Depth within the Lowest Aerosol Layer over the Tibetan Plateau from 2007 to 2014

The characteristics of aerosol optical depth (AOD) over the Tibetan Plateau (TP) were analyzed using 8-year (from January 2007 to December 2014) Cloud-Aerosol LiDAR and Infrared Pathfinder Satellite Observation (CALIPSO) level 2 aerosol layer products. Firstly, the overall feature of AOD over the Tibetan Plateau was investigated, including the seasonal diversities of AODS (the sum of AODs from all aerosol layers), and A (the amounts of aerosol layers). Then we deeply studied the characteristics of AOD within the lowest aerosol layer over TP, including the seasonal variations of AOD1 (The AOD of the first aerosol layer), HB1 (the height of the first aerosol layer base), TL1 (the thickness of the first aerosol layer) and PAOD1 (The AOD proportion of the first aerosol layer). The AODS was generally low ( 0.9) indicated that the aerosols were mainly concentrated in the lowest layer in summer, fall, and winter in the main body of TP. In spring, the PAOD1 value was relatively low (~0.7–0.85) and the distribution exhibited obvious differences between the southern (~0.85) and the northern (~0.75) TP, which appeared to be consistent with A. Most of the aerosol loads in summer were concentrated in the lowest aerosol layer with high aerosol loads. Most of the aerosol loads in fall and winter were also concentrated in the lowest aerosol layer, but with low aerosol loads.
Estimating Biomass and Nitrogen Amount of Barley and Grass Using UAV and Aircraft Based Spectral and Photogrammetric 3D Features

The timely estimation of crop biomass and nitrogen content is a crucial step in various tasks in precision agriculture, for example in fertilization optimization. Remote sensing using drones and aircrafts offers a feasible tool to carry out this task. Our objective was to develop and assess a methodology for crop biomass and nitrogen estimation, integrating spectral and 3D features that can be extracted using airborne miniaturized multispectral, hyperspectral and colour (RGB) cameras. We used the Random Forest (RF) as the estimator, and in addition Simple Linear Regression (SLR) was used to validate the consistency of the RF results. The method was assessed with empirical datasets captured of a barley field and a grass silage trial site using a hyperspectral camera based on the Fabry-Pérot interferometer (FPI) and a regular RGB camera onboard a drone and an aircraft. Agricultural reference measurements included fresh yield (FY), dry matter yield (DMY) and amount of nitrogen. In DMY estimation of barley, the Pearson Correlation Coefficient (PCC) and the normalized Root Mean Square Error (RMSE%) were at best 0.95% and 33.2%, respectively; and in the grass DMY estimation, the best results were 0.79% and 1.9%, respectively. In the nitrogen amount estimations of barley, the PCC and RMSE% were at best 0.97% and 21.6%, respectively. In the biomass estimation, the best results were obtained when integrating hyperspectral and 3D features, but the integration of RGB images and 3D features also provided results that were almost as good. In nitrogen content estimation, the hyperspectral camera gave the best results. We concluded that the integration of spectral and high spatial resolution 3D features and radiometric calibration was necessary to optimize the accuracy.
Remote Sensing of Sub-Surface Suspended Sediment Concentration by Using the Range Bias of Green Surface Point of Airborne LiDAR Bathymetry

Suspended sediment concentrations (SSCs) have been retrieved accurately and effectively through waveform methods by using green-pulse waveforms of airborne LiDAR bathymetry (ALB). However, the waveform data are commonly difficult to analyze. Thus, this paper proposes a 3D point-cloud method for remote sensing of SSCs in calm waters by using the range biases of green surface points of ALB. The near water surface penetrations (NWSPs) of green lasers are calculated on the basis of the green and reference surface points. The range biases (ΔS) are calculated by using the corresponding NWSPs and beam-scanning angles. In situ measured SSCs (C) and range biases (ΔS) are used to establish an empirical C-ΔS model at SSC sampling stations. The SSCs in calm waters are retrieved by using the established C-ΔS model. The proposed method is applied to a practical ALB measurement performed by Optech Coastal Zone Mapping and Imaging LiDAR. The standard deviations of the SSCs retrieved by the 3D point-cloud method are less than 20 mg/L.
Application of a Continuous Terrestrial Photogrammetric Measurement System for Plot Monitoring in the Beijing Songshan National Nature Reserve

Monitoring sample areas is the basis of ecological management. Songshan National Nature Reserve is one of the most important components of the ecosystem of the central metropolitan area of Beijing, and has fallen behind in its monitoring technology and methods. So, updating the existing equipment and technology is necessary. The current system suffers from high equipment costs and is not convenient to carry, so the work efficiency is low. Furthermore, the data cannot be visualized in three dimensions (3D), and complex terrain conditions cannot be measured. Therefore, this study researched and developed a continuous terrestrial photogrammetric measurement system that is theoretically based on the principles of photogrammetry, image processing technology, and dendrometry. The system applies a self-developed personal digital assistant (PDA) photogrammetry-based dendrometer and software to continuously evaluate stand sampling areas. Through experimental verification, the relative root mean square error (RMSE) of the trunk diameter measurements was found to be 5.59%, and the relative RMSE of hypsometrical measurements was 3.93%, which are both higher than the accuracy required for traditional forestry surveys. Furthermore, the advantages of this system include its low cost, lightweight equipment, easy operation, high measurement efficiency, 3D visualization, and applicability under complex terrain conditions. Since it is highly accurate and efficient, the continuous terrestrial photogrammetric system can be easily applied to monitor stand sampling areas in Songshan National Nature Reserve. In addition, it can be applied to second-class forest surveys in China, thus guaranteeing the monitoring of big data for the ecological environment of China.
Single Channel Circular SAR Moving Target Detection Based on Logarithm Background Subtraction Algorithm

In this paper, a new moving target detection algorithm for single channel circular Synthetic Aperture Radar (SAR) is presented. Circular SAR can perform observation with circular trajectory and retrieve multiaspect information. Due to the circular trajectory, the moving target has long and complex apparent trace (smear trajectory), while the background (clutter) has relatively slow changing characteristic in limited azimuth integration. Thus, the background image can be filtered by using radiometric adjusted overlap subaperture images. Then, the Signal-to-Clutter Ratio (SCR) of the region where the strong backscattering clutter exists is enhanced with logarithm background subtraction technique. In the real data experiment, the maximum SCR improvement is 13 dB. Finally, the moving target can be detected using a Constant False Alarm (CFAR) detector. The algorithm is evaluated with one channel data of three channel GOTCHA-GMTI dataset.
Sentinel-1 InSAR Measurements of Elevation Changes over Yedoma Uplands on Sobo-Sise Island, Lena Delta

Yedoma—extremely ice-rich permafrost with massive ice wedges formed during the Late Pleistocene—is vulnerable to thawing and degradation under climate warming. Thawing of ice-rich Yedoma results in lowering of surface elevations. Quantitative knowledge about surface elevation changes helps us to understand the freeze-thaw processes of the active layer and the potential degradation of Yedoma deposits. In this study, we use C-band Sentinel-1 InSAR measurements to map the elevation changes over ice-rich Yedoma uplands on Sobo-Sise Island, Lena Delta with frequent revisit observations (as short as six or 12 days). We observe significant seasonal thaw subsidence during summer months and heterogeneous inter-annual elevation changes from 2016–17. We also observe interesting patterns of stronger seasonal thaw subsidence on elevated flat Yedoma uplands by comparing to the surrounding Yedoma slopes. Inter-annual analyses from 2016–17 suggest that our observed positive surface elevation changes are likely caused by the delayed progression of the thaw season in 2017, associated with mean annual air temperature fluctuations.
Modeling Solar Radiation in the Forest Using Remote Sensing Data: A Review of Approaches and Opportunities

Solar radiation, the radiant energy from the sun, is a driving variable for numerous ecological, physiological, and other life-sustaining processes in the environment. Traditional methods to quantify solar radiation are done either directly (e.g., quantum sensors), or indirectly (e.g., hemispherical photography). This study, however, evaluates literature which utilized remote sensing (RS) technologies to estimate various forms of solar radiation or components, thereof under or within forest canopies. Based on the review, light detection and ranging (LiDAR) has, so far, been preferably used for modeling light under tree canopies. Laser system’s capability of generating 3D canopy structure at high spatial resolution makes it a reasonable choice as a source of spatial information about light condition in various parts of forest ecosystem. The majority of those using airborne laser system (ALS) commonly adopted the volumetric-pixel (voxel) method or the laser penetration index (LPI) for modeling the radiation, while terrestrial laser system (TLS) is preferred for canopy reconstruction and simulation. Furthermore, most of the studies focused only on global radiation, and very few on the diffuse fraction. It was also found out that most of these analyses were performed in the temperate zone, with a smaller number of studies made in tropical areas. Nonetheless, with the continuous advancement of technology and the RS datasets becoming more accessible and less expensive, these shortcomings and other difficulties of estimating the spatial variation of light in the forest are expected to diminish.
Assessment of Classifiers and Remote Sensing Features of Hyperspectral Imagery and Stereo-Photogrammetric Point Clouds for Recognition of Tree Species in a Forest Area of High Species Diversity

Recognition of tree species and geospatial information on tree species composition is essential for forest management. In this study, tree species recognition was examined using hyperspectral imagery from visible to near-infrared (VNIR) and short-wave infrared (SWIR) camera sensors in combination with a 3D photogrammetric canopy surface model based on RGB camera stereo-imagery. An arboretum with a diverse selection of 26 tree species from 14 genera was used as a test area. Aerial hyperspectral imagery and high spatial resolution photogrammetric color imagery were acquired from the test area using unmanned aerial vehicle (UAV) borne sensors. Hyperspectral imagery was processed to calibrated reflectance mosaics and was tested along with the mosaics based on original image digital number values (DN). Two alternative classifiers, a k nearest neighbor method (k-nn), combined with a genetic algorithm and a random forest method, were tested for predicting the tree species and genus, as well as for selecting an optimal set of remote sensing features for this task. The combination of VNIR, SWIR, and 3D features performed better than any of the data sets individually. Furthermore, the calibrated reflectance values performed better compared to uncorrected DN values. These trends were similar with both tested classifiers. Of the classifiers, the k-nn combined with the genetic algorithm provided consistently better results than the random forest algorithm. The best result was thus achieved using calibrated reflectance features from VNIR and SWIR imagery together with 3D point cloud features; the proportion of correctly-classified trees was 0.823 for tree species and 0.869 for tree genus.
Monitoring and Characterizing Heterogeneous Mediterranean Landscapes with Continuous Textural Indices Based on VHSR Imagery

Remote sensing tools (RS) can contribute to a better understanding of the diversity of natural and semi-naturals habitats, their spatial distribution, and their conservation status. RS can also provide a generic set of derived indicators to support local to regional habitat monitoring. Here we propose a set of synthetic continuous textural indices computed from high spatial resolution airborne images for the characterization of vegetation structure in very heterogeneous landscape mosaics. These indices are based on Fourier-based textural ordination (FOTO) of very-high-resolution images. We investigate the relationship between textural indices and a set of common landscape metrics derived from vegetation maps, identifying four strata of interest: bare soil, herbs, low ligneous, and high ligneous. We identify two continuous textural indices, the first one being related to vegetation strata fragmentation and the second being related to the dominance of high ligneous. The combination of these two textural indices with the Normalized Difference Vegetation Index (NDVI) provides a synoptic and accurate overview of the spatial organization of the different vegetation strata. The methodological approach presented herein has a generic value in response to national conservation targets in the context of mapping relevant habitat indicators.
Retrieving Three-Dimensional Co-Seismic Deformation of the 2017 Mw7.3 Iraq Earthquake by Multi-Sensor SAR Images

The Mw7.3 Iraq earthquake on 12 November 2017 was the largest recorded earthquake in the Zagros Mountains since 1900. In order to quantitatively analyze the co-seismic deformation caused by this earthquake, both the ascending and descending SAR images from the Japan Aerospace Exploration Agency’s ALOS-2 and the European Space Agency’s Sentinel-1A satellites were collected to implement the conventional differential interferometric synthetic aperture radar (DInSAR), multiple aperture InSAR (MAI), and azimuth pixel offset (AZO) methods. Subsequently, the three-dimensional (3D) deformation field was reconstructed over an area of about 60 × 70 km2 by a combined use of the line-of-sight (LOS) motion (detected by the DInSAR method) and the along-track (AT) motion (detected by the MAI method) through the weighted least square method. The experiment indicates that the ALOS-2 satellite performs better than the Sentinel-1A sensor in larger-magnitude earthquake deformation monitoring. Furthermore, the MAI method based on phase differencing has a better performance than the AZO method based on SAR amplitude correlation, as long as the coherence of the interferograms is sufficient. The maximum co-seismic displacements in the up–down, north–south, and east–west directions are approximately 100 cm, 100 cm, and −50 cm, respectively. After comparative analysis between the obtained 3D deformation field and the simulated deformation field with the fault parameters published by the USGS (United States Geological Survey), both co-seismic deformation fields are highly coincident, and the residuals between both (in different directions/dimensional) are in the magnitude of centimeters. Considering the geological structure in the earthquake region and factors of the LOS and 3D co-seismic deformation—such as the trend and location of the deformation bound, the different sign of displacements in hanging wall and footwall, and the locations of mainshock and aftershock—the preliminary conclusion is that the Zagros Mountain Front fault is responsible for the earthquake.