Miracle-dz/UniFAR-FTUs huggingface dataset:
Recent Changes of Glacial Lakes in the High Mountain Asia and Its Potential Controlling Factors Analysis

The current glacial lake datasets in the High Mountain Asia (HMA) region still need to be improved because their boundary divisions in the land–water transition zone are not precisely delineate, and also some very small glacial lakes have been lost due to their mixed reflectance with backgrounds. In addition, most studies have only focused on the changes in the area of a glacial lake as a whole, but do not involve the actual changes of per pixel on its boundary and the potential controlling factors. In this research, we produced more accurate and complete maps of glacial lake extent in the HMA in 2008, 2012, and 2016 with consistent time intervals using Landsat satellite images and the Google Earth Engine (GEE) cloud computing platform, and further studied the formation, distribution, and dynamics of the glacial lakes. In total, 17,016 and 21,249 glacial lakes were detected in 2008 and 2016, respectively, covering an area of 1420.15 ± 232.76 km2 and 1577.38 ± 288.82 km2; the lakes were mainly located at altitudes between 4400 m and 5600 m. The annual areal expansion rate was approximately 1.38% from 2008 to 2016. To explore the cause of the rapid expansion of individual glacial lakes, we investigated their long-term expansion rates by measuring changes in shoreline positions. The results show that glacial lakes are expanding rapidly in areas close to glaciers and had a high expansion rate of larger than 20 m/yr from 2008 to 2016. Glacial lakes in the Himalayas showed the highest expansion rate of more than 2 m/yr, followed by the Karakoram Mountains (1.61 m/yr) and the Tianshan Mountains (1.52 m/yr). The accelerating rate of glacier ice and snow melting caused by global warming is the primary contributor to glacial lake growth. These results may provide information that will help in the understanding of detailed lake dynamics and the mechanism, and also facilitate the scientific recognition of the potential hazards associated with glacial lakes in this region.
Vertical Structure of Air Pollutant Transport Flux as Determined by Ground-Based Remote Sensing Observations in Fen-Wei Plain, China

Air pollutant transport plays an important role in local air quality, but field observations of transport fluxes, especially their vertical distributions, are very limited. We characterized the vertical structures of transport fluxes in central Luoyang, Fen-Wei Plain, China, in winter based on observations of vertical air pollutant and wind profiles using multi-axis differential optical absorption spectroscopy (MAX-DOAS) and Doppler wind lidar, respectively. The northwest and the northeast are the two privileged wind directions. The wind direction and total transport scenarios were dominantly the northwest during clear days, turning to the northeast during the polluted days. Increased transport flux intensities of aerosol were found at altitudes below 400 m on heavily polluted days from the northeast to the southwest over the city. Considering pollution dependence on wind directions and speeds, surface-dominated northeast transport may contribute to local haze events. Northwest winds transporting clean air masses were dominant during clean periods and flux profiles characterized by high altitudes between 200 and 600 m in Luoyang. During the COVID-19 lockdown period in late January and February, clear reductions in transport flux were found for NO2 from the northeast and for HCHO from the northwest, while the corresponding main transport altitude remained unchanged. Our findings provide better understandings of regional transport characteristics, especially at different altitudes.
Changes in Urban Mobility Related to the Public Bike System with Regard to Weather Conditions and Statutory Retail Restrictions

The main purpose of the paper is to determine changes in transport behaviour of users of the public bike-share (PBS) scheme in a large Polish city, Łódź. By tracking GPS signals for individual trips taken by PBS users, it was possible to analyse their changeability (time and spatial) for periods before the implementation of statutory Sunday retail restrictions (2017) and after their partial introduction (2018). The study also took into account weather conditions, namely maximum and minimum daily temperatures and daily totals of precipitation recorded by a weather station in Łodź. In order to determine the correlations between certain weather conditions and PBS trips, the authors applied regression analysis. The results of the study showed that weekend cycling is less susceptible to the impact of weather than cycling on weekdays. At the same time, a comparative analysis of trading and non-trading Sundays proved that, during Sundays with retail restrictions, public bikes were used for longer, farther, and slower trips. These observations were confirmed by analyses of maps of traffic structure.
Deep Learning for Chlorophyll-a Concentration Retrieval: A Case Study for the Pearl River Estuary

The abundance of phytoplankton is generally estimated by measuring the chlorophyll-a concentration (Cchla), which is an important factor in photosynthesis and can be used to analyze the density and biomass of phytoplankton in the ecosystem. The band-ratio-based empirical or semi-analytical algorithms are operationally applied to retrieve Cchla in global oceans, which generally experience difficulties from the diversity of optical properties and the complexity of the radiative transfer equations in analytical analyses, respectively. With an attempt to develop an accurate Cchla retrieval model for the optically complex coastal and estuarine waters, this study aimed to explore the deep learning (DL) methods in satellite retrieval of Cchla. A two-stage convolutional neural network (CNN), named Cchla-Net, was proposed, which utilized the spectral information of remote sensing reflectances at MODIS/Aqua’s visible bands. In the first-stage phase, the Cchla-Net was pretrained by a set of remote sensing patches, in which the Cchla was generated from an existing model (OC3M). The pretrained results were than used as the initial values to refine the network with the synthetic oversampled in-situ dataset in the second-stage training phase. Using in-situ samples for training with the new initial values has a higher probability to reach the global optimum. The quantitative analyses showed that the two-stage training was more likely to achieve a global optimum in the optimization than the one-stage training. Matchups of the in-situ Cchla measurements were used to evaluate the retrieval models. Results showed that the proposed Cchla-Net produced obvious better performance than the empirical and semi-analytical algorithms, implying the DL method was more effective for optically complex waters with extremely high Cchla. This study provided an applicable method for remote sensing retrieval of Cchla, which should be helpful for studying the spatial distribution and temporal variability in the productive Pearl River estuary (PRE) waters.
A Spectral Spatial Attention Fusion with Deformable Convolutional Residual Network for Hyperspectral Image Classification

Convolutional neural networks (CNNs) have exhibited excellent performance in hyperspectral image classification. However, due to the lack of labeled hyperspectral data, it is difficult to achieve high classification accuracy of hyperspectral images with fewer training samples. In addition, although some deep learning techniques have been used in hyperspectral image classification, due to the abundant information of hyperspectral images, the problem of insufficient spatial spectral feature extraction still exists. To address the aforementioned issues, a spectral–spatial attention fusion with a deformable convolution residual network (SSAF-DCR) is proposed for hyperspectral image classification. The proposed network is composed of three parts, and each part is connected sequentially to extract features. In the first part, a dense spectral block is utilized to reuse spectral features as much as possible, and a spectral attention block that can refine and optimize the spectral features follows. In the second part, spatial features are extracted and selected by a dense spatial block and attention block, respectively. Then, the results of the first two parts are fused and sent to the third part, and deep spatial features are extracted by the DCR block. The above three parts realize the effective extraction of spectral–spatial features, and the experimental results for four commonly used hyperspectral datasets demonstrate that the proposed SSAF-DCR method is superior to some state-of-the-art methods with very few training samples.
Remote Sensing of Turbidity in the Tennessee River Using Landsat 8 Satellite

The Tennessee River in the United States is one of the most ecologically distinct rivers in the world and serves as a great resource for local residents. However, it is also one of the most polluted rivers in the world, and a leading cause of this pollution is storm water runoff. Satellite remote sensing technology, which has been used successfully to study surface water quality parameters for many years, could be very useful to study and monitor the quality of water in the Tennessee River. This study developed a numerical turbidity estimation model for the Tennessee River and its tributaries in Southeast Tennessee using Landsat 8 satellite imagery coupled with near real-time in situ measurements. The obtained results suggest that a nonlinear regression-based numerical model can be developed using Band 4 (red) surface reflectance values of the Landsat 8 OLI sensor to estimate turbidity in these water bodies with the potential of high accuracy. The accuracy assessment of the estimated turbidity achieved a coefficient of determination (R2) value and root mean square error (RMSE) as high as 0.97 and 1.41 NTU, respectively. The model was also tested on imagery acquired on a different date to assess its potential for routine remote estimation of turbidity and produced encouraging results with R2 value of 0.94 and relatively high RMSE.
Response of Vegetation Photosynthetic Phenology to Urbanization in Dongting Lake Basin, China

Urbanization can induce environmental changes such as the urban heat island effect, which in turn influence the terrestrial ecosystem. However, the effect of urbanization on the phenology of subtropical vegetation remains relatively unexplored. This study analyzed the changing trend of vegetation photosynthetic phenology in Dongting Lake basin, China, and its response to urbanization using nighttime light and chlorophyll fluorescence datasets. Our results indicated the start of the growing season (SOS) of vegetation in the study area was significantly advanced by 0.70 days per year, whereas the end of the growing season (EOS) was delayed by 0.24 days per year during 2000–2017. We found that urbanization promoted the SOS advance and EOS delay. With increasing urbanization intensity, the sensitivity of SOS to urbanization firstly increased then decreased, while the sensitivity of EOS to urbanization decreased with urbanization intensity. The climate sensitivity of vegetation phenology varied with urbanization intensity; urbanization induced an earlier SOS by increasing preseason minimum temperatures and a later EOS by increasing preseason precipitation. These findings improve our understanding of the vegetation phenology response to urbanization in subtropical regions and highlight the need to integrate human activities into future vegetation phenology models.
Liquid Hould-up in Packed Towers:The Generalized Correlation of Loading Velocities

ガスを流さぬ場合の灌液充填塔のOperating hold-upと灌液の物性及び操作條件との関係については前報にて報告した。 ガスを向流させた場合でもその流速の小さい範囲ではhold-upはガス速度に関係しないことがElginら及びPiretらによつて報告されているから,この関係式はそのまゝ適用できるわけである。 液量を一定にして,ガス流速を増して行つた場合の充填塔の性状についてはElginらによつて詳細に論じられており,これによると或る流速に至るとhold-upはガス流速と共に増加し始め,さらに或る限界速度に達すると再びhold-upは急激に増加してくる。 第1の屈折点がLoading pointであり,第2の屈折点がFlooding pointに相当する。 本報文では既往文献のhold-upの測定値よりLoading velocityを各種充填物について求め,Loading velocityと操作條件との一般的関係を導いたものである。 なおガス流速とhold-upとの関係についても言及した。
Testing the Height Variation Hypothesis with the R rasterdiv Package for Tree Species Diversity Estimation

Forest biodiversity is a key element to support ecosystem functions. Measuring biodiversity is a necessary step to identify critical issues and to choose interventions to be applied in order to protect it. Remote sensing provides consistent quality and standardized data, which can be used to estimate different aspects of biodiversity. The Height Variation Hypothesis (HVH) represents an indirect method for estimating species diversity in forest ecosystems from the LiDAR data, and it assumes that the higher the variation in tree height (height heterogeneity, HH), calculated through the ’Canopy Height Model’ (CHM) metric, the more complex the overall structure of the forest and the higher the tree species diversity. To date, the HVH has been tested exclusively with CHM data, assessing the HH only with a single heterogeneity index (the Rao’s Q index) without making use of any moving windows (MW) approach. In this study, the HVH has been tested in an alpine coniferous forest situated in the municipality of San Genesio/Jenesien (eastern Italian Alps) at 1100 m, characterized by the presence of 11 different tree species (mainly Pinus sylvestris, Larix decidua, Picea abies followed by Betula alba and Corylus avellana). The HH has been estimated through different heterogeneity measures described in the new R rasterdiv package using, besides the CHM, also other LiDAR metrics (as the percentile or the standard deviation of the height distribution) at various spatial resolutions and MWs (ALS LiDAR data with mean point cloud density of 2.9 point/m2). For each combination of parameters, and for all the considered plots, linear regressions between the Shannon’s H′ (used as tree species diversity index based on field data) and the HH have been derived. The results showed that the Rao’s Q index (singularly and through a multidimensional approach) performed generally better than the other heterogeneity indices in the assessment of the HH. The CHM and the LiDAR metrics related to the upper quantile point cloud distribution at fine resolution (2.5 m, 5 m) have shown the most important results for the assessment of the HH. The size of the used MW did not influence the general outcomes but instead, it increased when compared to the results found in the literature, where the HVH was tested without MW approach. The outcomes of this study underline that the HVH, calculated with certain heterogeneity indices and LiDAR data, can be considered a useful tool for assessing tree species diversity in considered forest ecosystems. The general results highlight the strength and importance of LiDAR data in assessing the height heterogeneity and the related biodiversity in forest ecosystems.
Digital Elevation Models: Terminology and Definitions

Digital elevation models (DEMs) provide fundamental depictions of the three-dimensional shape of the Earth’s surface and are useful to a wide range of disciplines. Ideally, DEMs record the interface between the atmosphere and the lithosphere using a discrete two-dimensional grid, with complexities introduced by the intervening hydrosphere, cryosphere, biosphere, and anthroposphere. The treatment of DEM surfaces, affected by these intervening spheres, depends on their intended use, and the characteristics of the sensors that were used to create them. DEM is a general term, and more specific terms such as digital surface model (DSM) or digital terrain model (DTM) record the treatment of the intermediate surfaces. Several global DEMs generated with optical (visible and near-infrared) sensors and synthetic aperture radar (SAR), as well as single/multi-beam sonars and products of satellite altimetry, share the common characteristic of a georectified, gridded storage structure. Nevertheless, not all DEMs share the same vertical datum, not all use the same convention for the area on the ground represented by each pixel in the DEM, and some of them have variable data spacings depending on the latitude. This paper highlights the importance of knowing, understanding and reflecting on the sensor and DEM characteristics and consolidates terminology and definitions of key concepts to facilitate a common understanding among the growing community of DEM users, who do not necessarily share the same background.
Remote Tracking of Distributed Dynamic Sources Over a Random Access Channel With One-Bit Updates

In this work, we consider a network, where distributed information sources whose states evolve according to a random process transmit their time-varying states to a remote estimator over a shared wireless channel. Each source generates packets in a decentralized manner and employs a slotted random access mechanism to transmit the packets. In particular, we are interested in networks with a large number of low-complexity devices that share low-capacity random access channels. Accordingly, we investigate update strategies for remote tracking of source states that require each update to constitute as few bits as possible. To that end, we develop update strategies requiring only one-bit of information per update that employ a local cancellation strategy. We further analytically compare the performance of the cancellation-enabled update policy to the optimal policy that does not restrict the number of bits for each update, which show that an asymptotic upper bound of the optimality ratio is $\frac{13\sqrt{2}}{12}$. Through simulations, we compare the proposed cancellation-enabled one-bit update policy with zero-wait sampling and threshold-based sampling policies that require more than one-bit of information per update. The comparisons show that the cancellation-enabled update policy at its optimal threshold level outperforms the multi-bit update policies.
Retrieval of Daily Mean Top-of-Atmosphere Reflected Solar Flux Using the Advanced Very High Resolution Radiometer (AVHRR) Instruments

The records of the Advanced Very High Resolution Radiometer (AVHRR) instrument observations can resolve the current lack of a long global climate data record of Reflected Solar Flux (RSF), by transforming these measurements into broadband flux at the top-of-atmosphere. This paper presents a methodology for obtaining daily mean RSF (Wm−2) from AVHRR. First, the narrowband reflectances are converted to broadband reflectance using empirical regressions with the Clouds and the Earth’s Radiant Energy System (CERES) observations. Second, the anisotropy is corrected by applying Angular Distribution Models (ADMs), which convert directional reflectance into a hemispherical albedo. Third, the instantaneous albedos are temporally interpolated by a flexible diurnal cycle model, capable of ingesting any number of observations at any time of day, making it suitable for any orbital configuration of NOAA and MetOp satellites. Finally, the twilight conditions prevailing near sunrise and sunset are simulated with an empirical model. The entire day is then integrated into a single daily mean RSF. This paper furthermore demonstrates the methodology by validating a full year (2008) of RSF daily means with the CERES SYN1deg data record, both on daily and subdaily scale. Several configurations are tested, each excluding particular satellites from the constellation in order to mimic orbital changes (e.g., orbital drift), and to assess their relative importance to the daily mean RSF. The best performance is obtained by the combination of at least one mid-morning (NOAA-17 or MetOp-A) and one early afternoon (NOAA-18) orbit. In this case, the RMS difference with CERES is about 7 Wm−2. Removing NOAA-18 degrades the performance to an RMS difference of 12 Wm−2, thereby providing an estimate of the impact of NOAA-19’s orbital drift between 2016 and 2020. Very early or late observations (NOAA-15, NOAA-16) provide little added value, and both mid-morning orbits turn out to be almost interchangeable given their close temporal proximity.
Approaching Global Instantaneous Precise Positioning with the Dual- and Triple-Frequency Multi-GNSS Decoupled Clock Model

Precise Point Positioning (PPP), as a global precise positioning technique, suffers from relatively long convergence times, hindering its ability to be the default precise positioning technique. Reducing the PPP convergence time is a must to reach global precise positions, and doing so in a few minutes to seconds can be achieved thanks to the additional frequencies that are being broadcast by the modernized GNSS constellations. Due to discrepancies in the number of signals broadcast by each satellite/constellation, it is necessary to have a model that can process a mix of signals, depending on availability, and perform ambiguity resolution (AR), a technique that proved necessary for rapid convergence. This manuscript does so by expanding the uncombined Decoupled Clock Model to process and fix ambiguities on up to three frequencies depending on availability for GPS, Galileo, and BeiDou. GLONASS is included as well, without carrier-phase ambiguity fixing. Results show the possibility of consistent quasi-instantaneous global precise positioning through an assessment of the algorithm on a network of global stations, as the 67th percentile solution converges below 10 cm horizontal error within 2 min, compared to 8 min with a triple-frequency solution, showing the importance of having a flexible PPP-AR model frequency-wise. In terms of individual datasets, 14% of datasets converge instantaneously when mixing dual- and triple-frequency measurements, compared to just 0.1% in that of dual-frequency mode without ambiguity resolution. Two kinematic car datasets were also processed, and it was shown that instantaneous centimetre-level positioning with a moving receiver is possible. These results are promising as they only rely on ultra-rapid global satellite products, allowing for instantaneous real-time precise positioning without the need for any local infrastructure or prior knowledge of the receiver’s environment.
Simulating the Leaf Area Index of Rice from Multispectral Images

Accurate estimation of the leaf area index (LAI) is essential for crop growth simulations and agricultural management. This study conducted a field experiment with rice and measured the LAI in different rice growth periods. The multispectral bands (B) including red edge (RE, 730 nm ± 16 nm), near-infrared (NIR, 840 nm ± 26 nm), green (560 nm ± 16 nm), red (650 nm ± 16 nm), blue (450 nm ± 16 nm), and visible light (RGB) were also obtained by an unmanned aerial vehicle (UAV) with multispectral sensors (DJI-P4M, SZ DJI Technology Co., Ltd.). Based on the bands, five vegetation indexes (VI) including Green Normalized Difference Vegetation Index (GNDVI), Leaf Chlorophyll Index (LCI), Normalized Difference Red Edge Index (NDRE), Normalized Difference Vegetation Index (NDVI), and Optimization Soil-Adjusted Vegetation Index (OSAVI) were calculated. The semi-empirical model (SEM), the random forest model (RF), and the Extreme Gradient Boosting model (XGBoost) were used to estimate rice LAI based on multispectral bands, VIs, and their combinations, respectively. The results indicated that the GNDVI had the highest accuracy in the SEM (R2 = 0.78, RMSE = 0.77). For the single band, NIR had the highest accuracy in both RF (R2 = 0.73, RMSE = 0.98) and XGBoost (R2 = 0.77, RMSE = 0.88). Band combination of NIR + red improved the estimation accuracy in both RF (R2 = 0.87, RMSE = 0.65) and XGBoost (R2 = 0.88, RMSE = 0.63). NDRE and LCI were the first two single VIs for LAI estimation using both RF and XGBoost. However, putting more than one VI together could only increase the LAI estimation accuracy slightly. Meanwhile, the bands + VIs combinations could improve the accuracy in both RF and XGBoost. Our study recommended estimating rice LAI by a combination of red + NIR + OSAVI + NDVI + GNDVI + LCI + NDRE (2B + 5V) with XGBoost to obtain high accuracy and overcome the potential over-fitting issue (R2 = 0.91, RMSE = 0.54).