Miracle-dz/UniFAR-FTUs huggingface dataset:
Remote Soil Moisture Measurement from Drone-Borne Reflectance Spectroscopy: Applications to Hydroperiod Measurement in Desert Playas

The extent, timing, and magnitude of soil moisture in wetlands (the hydropattern) is a primary physical control on biogeochemical processes in desert environments. However, determining playa hydropatterns is challenged by the remoteness of desert basin sites and by the difficulty in determining soil moisture from remotely sensed data at fine spatial and temporal scales (hundreds of meters to kilometers, and hours to days). Therefore, we developed a new, reflectance-based soil moisture index (continuum-removed water index, or CRWI) that can be determined via hyperspectral imaging from drone-borne platforms. We compared its efficacy at remotely determining soil moisture content to existing hyperspectral and multispectral soil moisture indices. CRWI varies linearly with in situ soil moisture content (R2 = 0.89, p < 0.001) and is comparatively insensitive to soil clay content (R2 = 0.4, p = 0.01), soil salinity (R2 = 0.82, p < 0.001), and soil grain size distribution (R2 = 0.67, p < 0.001). CRWI is negatively correlated with clay content, indicating it is not sensitive to hydrated mineral absorption features. CRWI has stronger correlation with surface soil moisture than other hyperspectral and multispectral indices (R2 = 0.69, p < 0.001 for WISOIL at this site). Drone-borne reflectance measurements allow monitoring of soil moisture conditions at the Alvord Desert playa test site over hectare-scale soil plots at measurement cadences of minutes to hours. CRWI measurements can be used to determine surface soil moisture at a range of desert sites to inform management decisions and to better reveal ecosystem processes in water-limited environments.
Vegetation Productivity Dynamics in Response to Climate Change and Human Activities under Different Topography and Land Cover in Northeast China

Net primary productivity (NPP) is the total amount of organic matter fixed by plants from the atmosphere through photosynthesis and is susceptible to the influences of climate change and human activities. In this study, we employed actual NPP (ANPP), potential NPP (PNPP), and human activity-induced NPP (HNPP) based on the Hurst exponent and statistical analysis to analyze the characteristics of vegetation productivity dynamics and to evaluate the effects of climate and human factors on vegetation productivity in Northeast China (NEC). The increasing trends in ANPP, PNPP, and HNPP accounted for 81.62%, 94.90%, and 89.63% of the total area, respectively, and ANPP in 68.64% of the total area will continue to increase in the future. Climate change played a leading role in vegetation productivity dynamics, which promoted an increase in ANPP in 71.55% of the area, and precipitation was the key climate factor affecting ANPP. The aggravation of human activities, such as increased livestock numbers and intensified agricultural activities, resulted in a decrease in ANPP in the western grasslands, northern Greater Khingan Mountains, and eastern Songnen Plain. In particular, human activities led to a decrease in ANPP in 53.84% of deciduous needleleaf forests. The impact of climate change and human activities varied significantly under different topography, and the percentage of the ANPP increase due to climate change decreased from 71.13% to 53.9% from plains to urgent slopes; however, the percentage of ANPP increase due to human activities increased from 3.44% to 21.74%, and the effect of human activities on the increase of ANPP was more obvious with increasing slope. At different altitudes, the difference in the effect of these two factors was not significant. The results are significant for understanding the factors influencing the vegetation productivity dynamics in NEC and can provide a reference for governments to implement projects to improve the ecosystem.
Usefulness of Global Root Zone Soil Moisture Product for Streamflow Prediction of Ungauged Basins

Using modelling approaches to predict stream flow from ungauged basins requires new model calibration strategies and evaluation methods that are different from the existing ones. Soil moisture information plays an important role in hydrological applications in basins. Increased availability of remote sensing data presents a significant opportunity to obtain the predictive performance of hydrological models (especially in ungauged basins), but there is still a limit to applying remote sensing soil moisture data directly to models. The Soil Moisture Active Passive (SMAP) satellite mission provides global soil moisture data estimated by assimilating remotely sensed brightness temperature to a land surface model. This study investigates the potential of a hydrological model calibrated using only global root zone soil moisture based on satellite observation when attempting to predict stream flow in ungauged basins. This approach’s advantage is that it is particularly useful for stream flow prediction in ungauged basins since it does not require observed stream flow data to calibrate a model. The modelling experiments were carried out on upstream watersheds of two dams in South Korea with high-quality stream flow data. The resulting model outputs when calibrated using soil moisture data without observed stream flow data are particularly impressive when simulating monthly stream flows upstream of the dams, and daily stream flows also showed a satisfactory level of predictive performance. In particular, the model calibrated using soil moisture data for dry years showed better predictive performance than for wet years. The performance of the model calibrated using soil moisture data was significantly improved under low flow conditions compared to the traditional regionalization approach. Additionally, the overall stream flow was also predicted better. In addition, the uncertainty of the model calibrated using soil moisture was not much different from that of the model calibrated using observed stream flow data, and showed more robust outputs when compared to the traditional regionalization approach. These results prove that the application of the global soil moisture product for predicting stream flows in ungauged basins is promising.
Deep Fully Convolutional Embedding Networks for Hyperspectral Images Dimensionality Reduction

Due to the superior spatial–spectral extraction capability of the convolutional neural network (CNN), CNN shows great potential in dimensionality reduction (DR) of hyperspectral images (HSIs). However, most CNN-based methods are supervised while the class labels of HSIs are limited and difficult to obtain. While a few unsupervised CNN-based methods have been proposed recently, they always focus on data reconstruction and are lacking in the exploration of discriminability which is usually the primary goal of DR. To address these issues, we propose a deep fully convolutional embedding network (DFCEN), which not only considers data reconstruction but also introduces the specific learning task of enhancing feature discriminability. DFCEN has an end-to-end symmetric network structure that is the key for unsupervised learning. Moreover, a novel objective function containing two terms—the reconstruction term and the embedding term of a specific task—is established to supervise the learning of DFCEN towards improving the completeness and discriminability of low-dimensional data. In particular, the specific task is designed to explore and preserve relationships among samples in HSIs. Besides, due to the limited training samples, inherent complexity and the presence of noise in HSIs, a preprocessing where a few noise spectral bands are removed is adopted to improve the effectiveness of unsupervised DFCEN. Experimental results on three well-known hyperspectral datasets and two classifiers illustrate that the low dimensional features of DFCEN are highly separable and DFCEN has promising classification performance compared with other DR methods.
Using A Sliding Window Phase Matching Method for Imaging of GNSS Radio Occultation Signals

Global Navigation Satellite System Radio Occultation (GNSS-RO) is a technique used to sound the atmosphere and derive vertical profiles of refractivity. Signals from GNSS satellites are received in a low-Earth orbit, and they are then processed to produce bending angle profiles, from which meteorological parameters can be retrieved. Generating two-dimensional images in the form of spectrograms from GNSS-RO signals is commonly done to, for instance, investigate reflections or estimate signal quality in the lower troposphere. This is typically implemented using, e.g., the Short-Time Fourier Transform (STFT) to produce a time-frequency representation that is subsequently transformed to bending angle (BA) and impact height (IH) coordinates by non-linear mapping. In this paper, we propose an alternative method based on a straightforward extension of the Phase Matching (PM) operator to produce two-dimensional spectral images in the BA-IH domain by applying a sliding window. This Sliding Window Phase Matching (SWPM) method generates the spectral amplitude on an arbitrary grid in BA and IH, e.g., along the coordinate axes. To illustrate, we show both SWPM and STFT methods applied to operational MetOp-A data. For SWPM we use a constant window in the BA-dimension, whereas for STFT we use a conventional constant time window. We show that the SWPM method produces the same result as STFT when the same window length is used for both methods. The sample points in impact parameter and bending angle are those generated by and the main advantage is that SWPM offers the user a convenient way to freely sample the BA-IH space. The cost for this is processing time that is somewhat longer than implementations based on the Fast Fourier Transform, such as the STFT method.
Meridional Changes in Satellite Chlorophyll and Fluorescence in Optically-Complex Coastal Waters of Northern Patagonia

Northern Patagonia is one of the largest estuarine systems worldwide. It is characterized by complex geography, including islands, peninsulas, channels, and fjords. Here, the Inner Sea of Chiloé (ISC) is the largest estuarine system extending about 230 km in the meridional direction. Phytoplankton’s long-term dynamics and the main physical drivers of their variability are not well understood yet. Time-space fluctuations of Chlorophyll-a (Chl-a) and Chlorophyll fluorescence (nFLH) within the ISC and their association with meteorological and oceanographic processes were analyzed using high resolution (1000 m) satellite data (2003–2019). Our results revealed a meridional Chl-a and nFLH gradient along the ISC, with higher concentrations north of the Desertores islands where the topography promotes a semi-closed system with estuarine characteristics yearlong. Satellite Chl-a and nFLH were characterized by asynchronous seasonal cycles (nFLH peaks in fall) that differed from the southern ISC where the maximum Chl-a and nFLH occurs in spring-summer. The adjacent coastal ocean influences the southern ISC, and thus, the Chl-a and nFLH variability correlated well with the seasonal variation of meridional winds. The northern ISC was clearly influenced by river discharges, which can bias the Chl-a retrievals, decoupling the annual cycles of Chl-a and nFLH. In situ data from a buoy in Seno Reloncaví reaffirmed this bias in satellite Chl-a and a higher correlation with nFLH, by which the construction of a local Chl-a algorithm for northern Patagonia is essential.
Variability in the Sea Surface Temperature Gradient and Its Impacts on Chlorophyll-a Concentration in the Kuroshio Extension

Sixteen years of satellite observational data in the Northwestern Pacific Ocean are used to describe the variability in the sea surface temperature (SST) gradient and its impact on chlorophyll-a concentrations (Chl-a). Spatially, a meridional dependence is identified in which the SST gradient increases to the north in association with elevated Chl-a. Temporally, the seasonal variability shows a large SST gradient and high Chl-a in winter and spring, while the SST gradient and Chl-a are much lower in summer. The seasonal variability in Chl-a leads the variability in the SST gradient by one month. A significant correlation between the SST gradient and Chl-a in the anomalous field is obtained only in the western section of the Kuroshio extension (KE) and the highest correlation is identified without any lags. An index for the section is defined as the proportion of the number of times that the SST gradient magnitude is anomalously large in each year, and the index is highly related to the stability of the KE and has a prominent influence on Chl-a in the region. An anomalously large positive (negative) SST gradient magnitude occurs when the KE is unstable (stable) and the corresponding Chl-a is high (low).
Simplified Priestley-Taylor Model to Estimate Land-Surface Latent Heat of Evapotranspiration from Incident Shortwave Radiation, Satellite Vegetation Index, and Air Relative Humidity

An operational and accurate model for estimating global or regional terrestrial latent heat of evapotranspiration (ET) across different land-cover types from satellite data is crucial. Here, a simplified Priestley–Taylor (SPT) model was developed without surface net radiation (Rn) by combining incident shortwave radiation (Rs), satellite vegetation index, and air relative humidity (RH). Ground-measured ET for 2000–2009 collected by 100 global FLUXNET eddy covariance (EC) sites was used to calibrate and evaluate the SPT model. A series of cross-validations demonstrated the reasonable performance of the SPT model to estimate seasonal and spatial ET variability. The coefficients of determination (R2) of the estimated versus observed daily (monthly) ET ranged from 0.42 (0.58) (p < 0.01) at shrubland (SHR) flux sites to 0.81 (0.86) (p < 0.01) at evergreen broadleaf forest (EBF) flux sites. The SPT model was applied to estimate agricultural ET at high spatial resolution (16 m) from Chinese Gaofen (GF)-1 data and monitor long-term (1982–2018) ET variations in the Three-River Headwaters Region (TRHR) of mainland China using the Global LAnd-Surface Satellite (GLASS) normalized difference vegetation index (NDVI) product. The proposed SPT model without Rn provides an alternative model for estimating regional terrestrial ET across different land-cover types.
Reconstruction of Snow Depth Data at Moderate Spatial Resolution (1 km) from Remotely Sensed Snow Data and Multiple Optimized Environmental Factors: A Case Study over the Qinghai-Tibetan Plateau

Snow depth distribution in the Qinghai-Tibetan plateau is important for atmospheric circulation and surface water resources. In-situ observations at meteorological stations and remote observation by passive microwave remote sensing technique are two main approaches for monitoring snow depth at regional or global levels. However, the meteorological stations are often scarce and unevenly distributed in mountainous regions because of inaccessibility, so are the in-situ snow depth measurements. Passive microwave remote sensing data can alleviate the unevenness issue, but accuracy and spatial (e.g., 25 km) and temporal resolutions are low; spatial heterogeneity in snow depth is thus hard to capture. On the other hand, optical sensors such as moderate resolution imaging spectroradiometer (MODIS) onboard Terra and Aqua satellites can monitor snow at moderate spatial resolution (1 km) and high temporal resolution (daily) but only snow area extent, not snow depth. Fusing passive microwave snow depth data with optical snow area extent data provides an unprecedented opportunity for generating snow depth data at moderate spatial resolution and high temporal resolution. In this article, a linear multivariate snow depth reconstruction (LMSDR) model was developed by fusing multisource snow depth data, optical snow area extent data, and environmental factors (e.g., spatial distribution, terrain features, and snow cover characteristics), to reconstruct daily snow depth data at moderate resolution (1 km) for 16 consecutive hydrological years, taking Qinghai-Tibetan Plateau (QTP) as a case study. We found that snow cover day (SCD) and environmental factors such as longitude, latitude, slope, surface roughness, and surface fluctuation have a significant impact on the variations of snow depth over the QTP. Relatively high accuracy (root mean square error (RMSE) = 2.26 cm) was observed in the reconstructed snow depth when compared with in-situ data. Compared with the passive microwave remote sensing snow depth product, constructing a nonlinear snow depletion curve product with an empirical formula and fusion snow depth product, the LMSDR model (RMSE = 2.28 cm, R2 = 0.63) demonstrated a significant improvement in accuracy of snow depth reconstruction. The overall spatial accuracy of the reconstructed snow depth was 92%. Compared with in-situ observations, the LMSDR product performed well regarding different snow depth intervals, land use, elevation intervals, slope intervals, and SCD and performed best, especially when the snow depth was less than 3 cm. At the same time, a long-time snow depth series reconstructed based on the LMSDR model reflected interannual variations of snow depth well over the QTP.
Anomaly Detection in Airborne Fourier Transform Thermal Infrared Spectrometer Images Based on Emissivity and a Segmented Low-Rank Prior

Although hyperspectral anomaly detection is commonly conducted in the visible, near-infrared, and shortwave infrared spectral regions, there has been less research on hyperspectral anomaly detection in the longwave infrared (LWIR) hyperspectral region. The radiance of thermal infrared hyperspectral imagery is determined by the temperature and emissivity. To avoid the detection uncertainty caused by the single factor of temperature, emissivity can be introduced to detect anomalies. However, in the emissivity domain, the spectral contrast and signal-to-noise ratio (SNR) are low, which makes it difficult to separate the anomalies from the background. In this paper, an anomaly detection method combining emissivity and a segmented low-rank prior (EaSLRP) is proposed for use with thermal infrared hyperspectral imagery. The EaSLRP method is divided into three parts—1) temperature/emissivity retrieval, 2) extraction of the thermal infrared hyperspectral background information, and 3) Mahalanobis distance detection. A homogeneous region generation method is also proposed to solve the problem of the complex global background leading to inaccurate background estimation. The GoDec method is used for matrix decomposition and background information extraction and to remove some of the noise. The proposed Mahalanobis distance detector then uses the background component and original image for anomaly detection, while highlighting the spectral difference between the anomalies and background. This method can also suppress the influence of noise, to some extent. The experimental results obtained with airborne Fourier transform thermal infrared spectrometer hyperspectral images demonstrate that the EaSLRP method is effective when compared with the Reed–Xiaoli detector (RXD), the segmented RX detector (SegRX), the low-rank and sparse representation-based detector (LRASR), the low-rank and sparse matrix decomposition (LRaSMD)-based Mahalanobis distance method (LSMAD), and the locally enhanced low-rank prior method (LELRP-AD).
Performance Evaluation of Different SAR-Based Techniques on the 2019 Ridgecrest Sequence

We evaluated the performances of different SAR-based techniques by analyzing the surface coseismic displacement related to the 2019 Ridgecrest seismic sequence (an Mw 6.4 foreshock on July 4th and an Mw 7.1 mainshock on July 6th) in the tectonic framework of the eastern California shear zone (Southern California, USA). To this end, we compared and validated the retrieved SAR-based coseismic displacement with the one estimated by a dense GNSS network, extensively covering the study area. All the SAR-based techniques constrained the surface fault rupture well; however, in comparison with the GNSS-based coseismic displacement, some significant differences were observed. InSAR data showed better performance than MAI and POT data by factors of about two and three, respectively, therefore confirming that InSAR is the most consolidated technique to map surface coseismic displacements. However, MAI and POT data made it possible to better constrain the azimuth displacement and to retrieve the surface rupture trace. Therefore, for cases of strike-slip earthquakes, all the techniques should be exploited to achieve a full synoptic view of the coseismic displacement field.
Comparison of Support Vector Machines and Random Forests for Corine Land Cover Mapping

Land cover information is essential in European Union spatial management, particularly that of invasive species, natural habitats, urbanization, and deforestation; therefore, the need for accurate and objective data and tools is critical. For this purpose, the European Union’s flagship program, the Corine Land Cover (CLC), was created. Intensive works are currently being carried out to prepare a new version of CLC+ by 2024. The geographical, climatic, and economic diversity of the European Union raises the challenge to verify various test areas’ methods and algorithms. Based on the Corine program’s precise guidelines, Sentinel-2 and Landsat 8 satellite images were tested to assess classification accuracy and regional and spatial development in three varied areas of Catalonia, Poland, and Romania. The method is dependent on two machine learning algorithms, Random Forest (RF) and Support Vector Machine (SVM). The bias of classifications was reduced using an iterative of randomized training, test, and verification pixels. The ease of the implementation of the used algorithms makes reproducing the results possible and comparable. The results show that an SVM with a radial kernel is the best classifier, followed by RF. The high accuracy classes that can be updated and classes that should be redefined are specified. The methodology’s potential can be used by developers of CLC+ products as a guideline for algorithms, sensors, and the possibilities and difficulties of classifying different CLC classes.
Best Matching Algorithm to Identify and Rank the Relevant Statutes

Automated identification of relevant documents or statutes for the given query is a time efficient process. Artificial Intelligence in Legal Assistance task (AILA) task is about using the artificial intelligence to identify the suitable prior documents or statute for the given query and semantic segmentation of the legal documents. We have identified the suitable statute for the given query using Best Matching (BM-25) ranking algorithm. We have calculated the scores using BM25 algorithm and sorted the statutes according to the scores and then identified the suitable statutes with better MAP score of 0.2975 than most of the other approaches on AILA@FIRE-2020 dataset.