[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (3,806)

Search Parameters:
Keywords = land surface model

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
28 pages, 8072 KiB  
Article
Quantifying Evapotranspiration and Environmental Factors in the Abandoned Saline Farmland Using Landsat Archives
by Liya Zhao, Jingwei Wu, Qi Yang, Hang Zhao, Jun Mao, Ziyang Yu, Yanqi Liu and Anne Gobin
Land 2025, 14(2), 283; https://doi.org/10.3390/land14020283 - 30 Jan 2025
Viewed by 295
Abstract
This study investigates the complex interaction of biophysical and meteorological factors that drive evapotranspiration (ET) in saline environments. Leveraging a total of 182 cloud-free Landsat 5/8 time-series data from 1988 to 2019, we employed the Surface Energy Balance System (SEBS) model to quantify [...] Read more.
This study investigates the complex interaction of biophysical and meteorological factors that drive evapotranspiration (ET) in saline environments. Leveraging a total of 182 cloud-free Landsat 5/8 time-series data from 1988 to 2019, we employed the Surface Energy Balance System (SEBS) model to quantify ET and investigate its relationships with soil salinity, vegetation cover, groundwater depth, and landscape metrics. We validated the predicted ET at two experimental sites using ET observation calculated by a water balance model. The result shows an R2 of 0.78 and RMSE of 0.91 mm for the SEBS predicted ET, indicating high accuracy of the ET estimation. We detected abandoned saline farmland patches across Hetao and extracted the normalized difference vegetation index (NDVI), salinization index (SI), and the predicted ET for analysis. The results indicate that ET is negatively correlated with SI with a Pearson correlation coefficient (r) up to −0.7, while ET is positively correlated with NDVI (r = 0.4). In addition, we designed a control-variable experiment in the Yichang subdistrict to investigate the effects of groundwater depth, land aggregation index, soil salinity index, and the area of abandoned saline farmland patches on ET. The results indicate that increased NDVI could significantly enhance ET, while smaller saline farmland patches exhibited greater sensitivity to groundwater recharge, with higher averaged ET than larger patches. Moreover, we analyzed factor importance using Lasso regression and Random Forest (RF) regression. The result shows that the ranking of the importance of the features is consistent for both methods and for all the features, with NDVI being the most important (with an RF importance score of 0.4), followed by groundwater table depth (GWTD), and the influence of the surface area of abandoned saline farmland being the weakest. We found that smaller patches of abandoned saline farmland were more sensitive to changes in groundwater levels induced by nearby irrigation, affecting their averaged ET more dynamically than larger patches. Decreasing patch size over time indicates ongoing changes in land management and ecological conditions. This study, through a multifactor analysis of ET in abandoned saline farmland and its intrinsic factors, provides a reference for evaluating the dry drainage efficiency of abandoned saline farmland in a dry drainage system. Full article
(This article belongs to the Special Issue Salinity Monitoring and Modelling at Different Scales: 2nd Edition)
Show Figures

Figure 1

Figure 1
<p>Abandoned saline farmland (<b>a</b>) and its surface (<b>b</b>), photographed in the Hetao Irrigation District of Inner Mongolia, China.</p>
Full article ">Figure 2
<p>Schematic diagram of a dry drainage system, where the abandoned farmland acts as an evaporation sink with a lower groundwater table. The groundwater and dissolved salts move from the surrounding crop field to the abandoned farmland, which provides a natural-based solution for sustainable farming.</p>
Full article ">Figure 3
<p>Location of land cover cases and five sub-irrigation districts of the Hetao Irrigation District (<b>a</b>); Location of the experimental site and well locations of the Yichang subdistrict (<b>b</b>).</p>
Full article ">Figure 4
<p>Land cover of five cases with a surface area of 5765 ha. The location of the cases is shown in <a href="#land-14-00283-f003" class="html-fig">Figure 3</a>a, where (<b>a</b>) Case in Yichang sub-irrigation district, (<b>b</b>) Case in Yichang sub-irrigation district, (<b>c</b>) Case in Urat sub-irrigation district, (<b>d</b>) Case in Jiefangzha sub-irrigation district, and (<b>e</b>) Case in Wulanbuhe sub-irrigation district.</p>
Full article ">Figure 5
<p>Distribution of abandoned saline farmland in the Hetao Irrigation District, (<b>a</b>) 1988, (<b>b</b>) 1998, (<b>c</b>) 2008, (<b>d</b>) 2018.</p>
Full article ">Figure 6
<p>Daily ET map for 2019: (<b>a</b>) 21 April, (<b>b</b>) 8 June, (<b>c</b>) 10 July, (<b>d</b>) 27 August, (<b>e</b>) 28 September, (<b>f</b>) 30 October.</p>
Full article ">Figure 7
<p>Non-freezing period ET inversion values from 1988 to 2019: (<b>a</b>) 1988, (<b>b</b>) 1993, (<b>c</b>) 1998, (<b>d</b>) 2006, (<b>e</b>) 2011, (<b>f</b>) 2016, (<b>g</b>) 2018, (<b>h</b>) 2019.</p>
Full article ">Figure 8
<p>Scatterplot of ET estimated from field-scale water balance and remote sensing- based SEBS.</p>
Full article ">Figure 9
<p>Scatter density plots of ET versus SI for abandoned saline farmland patches in Hetao in (<b>a</b>) 1990, (<b>b</b>) 1999, (<b>c</b>) 2010, and (<b>d</b>) 2019.</p>
Full article ">Figure 10
<p>Scatter density plots of ET versus NDVI for abandoned saline farmland patches in Hetao in (<b>a</b>) 1990, (<b>b</b>) 1999, (<b>c</b>) 2010, and (<b>d</b>) 2019.</p>
Full article ">Figure 11
<p>Scatter density plot of annual average groundwater depth (m) and ET (mm) for abandoned saline farmland patches in 2019.</p>
Full article ">Figure 12
<p>Number of patches in 1988–2019 (<b>a</b>), number distribution of abandoned saline farmland (patches) in 2019 (<b>b</b>).</p>
Full article ">Figure 13
<p>Distribution of averaged ET (mm) over multiple years (1988–2011, 2013–2019) for different abandoned saline area patches.</p>
Full article ">Figure 14
<p>Non-freezing season ET (mm) for different area ranges of abandoned saline farmland in (<b>a</b>) 1990, (<b>b</b>) 1999, (<b>c</b>) 2010, and (<b>d</b>) 2019. The values of the independent variable use left-closed and right-open intervals.</p>
Full article ">Figure 15
<p>The relationship between AI and ET is based on five landscape cases in the Hetao Irrigation District. The blue dashed line shows the quadratic regression fit.</p>
Full article ">Figure 16
<p>Feature importance analysis based on (<b>a</b>) Lasso multivariate linear regression and (<b>b</b>) Random Forest regression.</p>
Full article ">Figure 17
<p>SHAP feature importance for the Lasso multivariate regression model.</p>
Full article ">Figure 18
<p>SHAP feature importance for the Random Forest (RF) model.</p>
Full article ">
21 pages, 4807 KiB  
Article
Spatial Analysis of Carbon Metabolism in Different Economic Divisions Based on Land Use and Cover Change (LUCC) in China
by Cui Yuan, Yaju Liu, Jingzhao Lu, Chengyi Guo, Tingting Quan and Wei Su
Atmosphere 2025, 16(2), 148; https://doi.org/10.3390/atmos16020148 - 29 Jan 2025
Viewed by 285
Abstract
Urbanization has greatly altered Earth’s surface form, and land use changes can lead to significant changes in carbon emissions. However, how these changes affect ecosystems remains unclear. Therefore, this study calculated the carbon absorption and emissions in 31 Chinese provinces using high-resolution (300 [...] Read more.
Urbanization has greatly altered Earth’s surface form, and land use changes can lead to significant changes in carbon emissions. However, how these changes affect ecosystems remains unclear. Therefore, this study calculated the carbon absorption and emissions in 31 Chinese provinces using high-resolution (300 m) land use data. Subsequently, a carbon flow model was used to evaluate the carbon transfer that occurred from the changes in land use in every province between 2000 and 2020. The standard deviation ellipse analytic techniques were also employed to research the spatiotemporal evolution features of carbon flow in various economic zones. Furthermore, the flux and utility analysis approaches in ecological network analysis were used to quantitatively examine the interaction relationship between two carbon metabolism land uses. The results revealed that the continuous expansion of China’s construction land has reduced the area of agricultural land, resulting in industrial land (53.14%) and urban land (39.38%) being the main contributors to the total carbon emissions. Among them, the five eastern provinces of Hebei, Jiangsu, Zhejiang, Shandong, and Guangdong had carbon emissions of more than 100 million tons. From 2000 to 2020, the center of gravity of the carbon flow in construction land had shifted significantly from Henan Province to Gansu Province. The ecological relationship of exploitation and control dominated the two land use types. It is mostly found in Xinjiang, Qinghai, Gansu, Inner Mongolia, and Ningxia provinces. The findings could provide relevant policy implications for the Chinese government to mitigate carbon metabolism on land. Full article
Show Figures

Figure 1

Figure 1
<p>Study area map.</p>
Full article ">Figure 2
<p>Land use spatial distribution (<b>a</b>) and land use transfer area (<b>b</b>) in China from 2000 to 2020. (A is agricultural land, F is forest land, G is grassland, W is water area and wetland, C is construction land, and B is bare ground).</p>
Full article ">Figure 3
<p>Carbon absorption and carbon emissions by provinces in China.</p>
Full article ">Figure 4
<p>The shift in gravity of different land types in China from 2000 to 2020 (<b>a</b>–<b>g</b>). (A is agricultural land, F is forest land, G is grassland, W is water area and wetland, C is construction land, and B is bare ground. The dashed box in (<b>d,f,g</b>) indicates the province range).</p>
Full article ">Figure 5
<p>The standard deviation ellipse of carbon flow from 2000 to 2020 (<b>a</b>–<b>e</b>). (α is the longitude direction and β is the latitude direction).</p>
Full article ">Figure 6
<p>Ecological relationships of land use types in China. (Yellow is mutualism, blue is exploitation, green is control, red is competition. “−, +” represents the positive and negative values of the elements in the complete utility matrix U. a is 2000–2005, b is 2005–2010, c is 2010–2015, d is 2015–2020).</p>
Full article ">Figure 7
<p>Spatial distribution of ecological relations of land use types in China. (The four parts of the gray area from deep to shallow represent northeast region, eastern region, central region and western region respectively).</p>
Full article ">
22 pages, 11317 KiB  
Article
Planning for a Warmer Future: Heat Risk Assessment and Mitigation in Lahti, Finland
by Ankur Negi, Rohinton Emmanuel and Eeva Aarrevaara
Atmosphere 2025, 16(2), 146; https://doi.org/10.3390/atmos16020146 - 29 Jan 2025
Viewed by 340
Abstract
With global climate change causing temperature increases, even cooler regions like Finland are facing increasing heat risks. The city of Lahti is expected to experience a higher-than-average temperature increase, making heat risk mitigation essential. This study aims to assess present and future heat [...] Read more.
With global climate change causing temperature increases, even cooler regions like Finland are facing increasing heat risks. The city of Lahti is expected to experience a higher-than-average temperature increase, making heat risk mitigation essential. This study aims to assess present and future heat risks in Lahti using exposure and social vulnerability indicators to identify heat risk hotspots and provide strategies for mitigation within the city’s urban planning framework. The method utilizes a combination of Land Surface Temperature (LST) data (2014–2024), climate projections, and microclimate analysis to identify heat risk in the city. Geographic Information Systems (GIS) and ENVI-met modeling were employed to assess the relationship between land surface temperatures (LST), urban structure, and green infrastructure. Risk assessments were conducted using social and environmental vulnerability indicators, and future projections were based on a combined SSP2-4.5 scenario. The results show a significant increase in high-risk areas by 2040, rising from 9.79% to 23.65% of Lahti’s core urban area. Although the current urban planning framework of the city (Masterplan 2035) is effective in terms of maintaining exposure levels, the continued increase in projected air temperatures, as modeled based on outputs of the EC-Earth3-veg GCM, remains a concern. Microclimate modeling confirmed that urban greenery significantly reduces heat stress and improves thermal comfort. To address future heat risks, Lahti must integrate more green infrastructure into its urban design and identify seasonal heat mitigation methodologies. Additionally, the findings emphasize the need for adaptive planning strategies to mitigate rising temperatures and ensure urban resilience. Full article
(This article belongs to the Section Biometeorology and Bioclimatology)
Show Figures

Figure 1

Figure 1
<p>Location of Lahti and Study Area in the urban core.</p>
Full article ">Figure 2
<p>UTCI model and equivalent temperature scale [<a href="#B39-atmosphere-16-00146" class="html-bibr">39</a>].</p>
Full article ">Figure 3
<p>Changing temperatures in Finland and Lahti (Data Source (<b>c</b>): Turku Urban Climate Research Group (TURCLIM), Geography Division, University of Turku, Finland).</p>
Full article ">Figure 4
<p>Land Surface Temperature profile of the Study area.</p>
Full article ">Figure 5
<p>(<b>a</b>) Present Risk Profile of the Study Area; (<b>b</b>) Future Heat Risk Profile of the Study Area.</p>
Full article ">Figure 5 Cont.
<p>(<b>a</b>) Present Risk Profile of the Study Area; (<b>b</b>) Future Heat Risk Profile of the Study Area.</p>
Full article ">Figure 6
<p>Bivariate Analysis for Hazard and Exposure.</p>
Full article ">Figure 7
<p>Bivariate Analysis for Hazard and Vulnerability.</p>
Full article ">Figure 8
<p>Bivariate Analysis for Vulnerability and Exposure.</p>
Full article ">Figure 9
<p>Potential air temperature and Mean Radiant temperature for sample points.</p>
Full article ">Figure 10
<p>UTCI comparison for sample points for Present and Future Scenarios.</p>
Full article ">Figure 11
<p>MRT and UTCI for high-albedo street simulation.</p>
Full article ">Figure 12
<p>Scatterplot for Vulnerability index and modeled future risk.</p>
Full article ">Figure 13
<p>Scatterplot for Exposure index and modeled future risk.</p>
Full article ">Figure 14
<p>Scatterplot for Hazard index and modeled future risk.</p>
Full article ">
21 pages, 3449 KiB  
Article
Indian Land Carbon Sink Estimated from Surface and GOSAT Observations
by Lorna Nayagam, Shamil Maksyutov, Rajesh Janardanan, Tomohiro Oda, Yogesh K. Tiwari, Gaddamidi Sreenivas, Amey Datye, Chaithanya D. Jain, Madineni Venkat Ratnam, Vinayak Sinha, Haseeb Hakkim, Yukio Terao, Manish Naja, Md. Kawser Ahmed, Hitoshi Mukai, Jiye Zeng, Johannes W. Kaiser, Yu Someya, Yukio Yoshida and Tsuneo Matsunaga
Remote Sens. 2025, 17(3), 450; https://doi.org/10.3390/rs17030450 - 28 Jan 2025
Viewed by 390
Abstract
The carbon sink over land plays a key role in the mitigation of climate change by removing carbon dioxide (CO2) from the atmosphere. Accurately assessing the land sink capacity across regions should contribute to better future climate projections and help guide [...] Read more.
The carbon sink over land plays a key role in the mitigation of climate change by removing carbon dioxide (CO2) from the atmosphere. Accurately assessing the land sink capacity across regions should contribute to better future climate projections and help guide the mitigation of global emissions towards the Paris Agreement. This study estimates terrestrial CO2 fluxes over India using a high-resolution global inverse model that assimilates surface observations from the global observation network and the Indian subcontinent, airborne sampling from Brazil, and data from the Greenhouse gas Observing SATellite (GOSAT) satellite. The inverse model optimizes terrestrial biosphere fluxes and ocean-atmosphere CO2 exchanges independently, and it obtains CO2 fluxes over large land and ocean regions that are comparable to a multi-model estimate from a previous model intercomparison study. The sensitivity of optimized fluxes to the weights of the GOSAT satellite data and regional surface station data in the inverse calculations is also examined. It was found that the carbon sink over the South Asian region is reduced when the weight of the GOSAT data is reduced along with a stricter data filtering. Over India, our result shows a carbon sink of 0.040 ± 0.133 PgC yr−1 using both GOSAT and global surface data, while the sink increases to 0.147 ± 0.094 PgC yr−1 by adding data from the Indian subcontinent. This demonstrates that surface observations from the Indian subcontinent provide a significant additional constraint on the flux estimates, suggesting an increased sink over the region. Thus, this study highlights the importance of Indian sub-continental measurements in estimating the terrestrial CO2 fluxes over India. Additionally, the findings suggest that obtaining robust estimates solely using the GOSAT satellite data could be challenging since the GOSAT satellite data yield significantly varies over seasons, particularly with increased rain and cloud frequency. Full article
(This article belongs to the Special Issue Remote Sensing of Carbon Fluxes and Stocks II)
Show Figures

Figure 1

Figure 1
<p>The locations of the Indian observation sites (orange squares), global surface observations from ObsPack (blue dots), and GOSAT (grey dots) observations. The data site SNG, used for independent validation, is also shown (red star).</p>
Full article ">Figure 2
<p>The South Asian flux estimates obtained by a subset of the OCO-2 MIP models (see [<a href="#B15-remotesensing-17-00450" class="html-bibr">15</a>] for individual model details) and our model. CT, TM5-4DVAR, CAMS, OU, and Baker are the selected OCO-2 MIP models. ISG, ISGH, and I4SG represent the different inversions carried out with the NTFVAR model.</p>
Full article ">Figure 3
<p>Time series of observation (green), forward (red), and optimized (blue) simulations for the selected background sites Syowa (SYO), Pallas (PAL), Hyytiälä (SMR) and Lampedusa (LMP). The smoothed values are the weekly averages of daily measures.</p>
Full article ">Figure 4
<p>Time series of observed (obs), forward (fwd, with prior fluxes), and inversion-optimized (opt) CO<sub>2</sub> concentration for the four Indian sites and Comilla, Bangladesh. The smoothed values are the weekly averages of daily measures.</p>
Full article ">Figure 5
<p>Monthly BIAS averaged for all Indian sites and Comilla (unit is ppm).</p>
Full article ">Figure 6
<p>Monthly bias for each Indian site and Comilla from inversions ISG, ISGH, and I4SG (unit ppm).</p>
Full article ">Figure 7
<p>The monthly posterior bias of GOSAT measurements averaged for each 10° latitude bin obtained from inversion ISG (unit ppm). The latitude given in the x-axis is the lower limit of each 10° bin.</p>
Full article ">Figure 8
<p>Same plots as <a href="#remotesensing-17-00450-f004" class="html-fig">Figure 4</a> but for SNG validation site.</p>
Full article ">
23 pages, 9644 KiB  
Article
Modeling Urban Microclimates for High-Resolution Prediction of Land Surface Temperature Using Statistical Models and Surface Characteristics
by Md Golam Rabbani Fahad, Maryam Karimi, Rouzbeh Nazari and Mohammad Reza Nikoo
Urban Sci. 2025, 9(2), 28; https://doi.org/10.3390/urbansci9020028 - 28 Jan 2025
Viewed by 478
Abstract
Surface properties in complex urban environments can significantly impact local-level temperature gradients and distribution on several scales. Studying temperature anomalies and identifying heat pockets in urban settings is challenging. Limited high-resolution datasets are available that do not translate into an accurate assessment of [...] Read more.
Surface properties in complex urban environments can significantly impact local-level temperature gradients and distribution on several scales. Studying temperature anomalies and identifying heat pockets in urban settings is challenging. Limited high-resolution datasets are available that do not translate into an accurate assessment of near-surface temperature. This study developed a model to predict land surface temperature (LST) at a high spatial–temporal resolution in urban areas using Landsat data and meteorological inputs from NLDAS. This study developed an urban microclimate (UC) model to predict air temperature at high spatial–temporal resolution for inner urban areas through a land surface and build-up scheme. The innovative aspect of the model is the inclusion of micro-features in land use characteristics, which incorporate surface types, urban vegetation, building density and heights, short wave radiation, and relative humidity. Statistical models, including the Generalized Additive Model (GAM) and spatial autoregression (SAR), were developed to predict land surface temperature (LST) based on surface characteristics and weather parameters. The model was applied to urban microclimates in densely populated regions, focusing on Manhattan and New York City. The results indicated that the SAR model performed better (R2 = 0.85, RMSE = 0.736) in predicting micro-scale LST variations compared to the GAM (R2 = 0.39, RMSE = 1.203) and validated the accuracy of the LST prediction model with R2 ranging from 0.79 to 0.95. Full article
Show Figures

Figure 1

Figure 1
<p>Administrative boundary of the borough of Manhattan in NYC and dominant land use types based on the latest National Land Cover Database (NLCD 2019).</p>
Full article ">Figure 2
<p>Key steps, required data processing, and methods for the proposed urban meteorological model.</p>
Full article ">Figure 3
<p>Snapshot of converted land surface temperature (°C) of New York City mapped within the zip codes.</p>
Full article ">Figure 4
<p>Comparison of downscaled NLDAS temperature with observed weather station data within three buffer zones (i.e., 1 km, 300 m, and 100 m).</p>
Full article ">Figure 5
<p>Results from Moran’s I index for spatial autocorrelation.</p>
Full article ">Figure 6
<p>GAM predicted temperature map, UC model predicted map results, and actual calculated temperature map using Landsat images.</p>
Full article ">Figure 7
<p>Scatter plots with R2 and RMSE values for observed vs. predicted values for the Landsat images.</p>
Full article ">
50 pages, 68393 KiB  
Article
Improved Stereophotogrammetric and Multi-View Shape-from-Shading DTMs of Occator Crater and Its Interior Cryovolcanism-Related Bright Spots
by Alicia Neesemann, Stephan van Gasselt, Ralf Jaumann, Julie C. Castillo-Rogez, Carol A. Raymond, Sebastian H. G. Walter and Frank Postberg
Remote Sens. 2025, 17(3), 437; https://doi.org/10.3390/rs17030437 - 27 Jan 2025
Viewed by 305
Abstract
Over the course of NASA’s Dawn Discovery mission, the onboard framing camera mapped Ceres across a wide wavelength spectrum at varying polar science orbits and altitudes. With increasing resolution, the uniqueness of the 92 km wide, young Occator crater became evident. Its central [...] Read more.
Over the course of NASA’s Dawn Discovery mission, the onboard framing camera mapped Ceres across a wide wavelength spectrum at varying polar science orbits and altitudes. With increasing resolution, the uniqueness of the 92 km wide, young Occator crater became evident. Its central cryovolcanic dome, Cerealia Tholus, and especially the associated bright carbonate and ammonium chloride deposits—named Cerealia Facula and the thinner, more dispersed Vinalia Faculae—are the surface expressions of a deep brine reservoir beneath Occator. Understandably, this made this crater the target for future sample return mission studies. The planning and preparation for this kind of mission require the characterization of potential landing sites based on the most accurate topography and orthorectified image data. In this work, we demonstrate the capabilities of the freely available and open-source USGS Integrated Software for Imagers and Spectrometers (ISIS 3) and Ames Stereo Pipeline (ASP 2.7) in creating high-quality image data products as well as stereophotogrammetric (SPG) and multi-view shape-from-shading (SfS) digital terrain models (DTMs) of the aforementioned spectroscopically challenging features. The main data products of our work are four new DTMs, including one SPG and one SfS DTM based on High-Altitude Mapping Orbit (HAMO) (CSH/CXJ) and one SPG and one SfS DTM based on Low-Altitude Mapping Orbit (LAMO) (CSL/CXL), along with selected Extended Mission Orbit 7 (XMO7) framing camera (FC) data. The SPG and SfS DTMs were calculated to a GSD of 1 and 0.5 px, corresponding to 136 m (HAMO SPG), 68 m (HAMO SfS), 34 m (LAMO SPG), and 17 m (LAMO SfS). Finally, we show that the SPG and SfS approaches we used yield consistent results even in the presence of high albedo differences and highlight how our new DTMs differ from those previously created and published by the German Aerospace Center (DLR) and the Jet Propulsion Laboratory (JPL). Full article
20 pages, 11525 KiB  
Article
Investigation of Spatio-Temporal Simulation of Mining Subsidence and Its Determinants Utilizing the RF-CA Model
by Jikun Xu, Chaode Yan, Baowei Zhang, Xuanchi Chen, Xu Yan, Rongxing Wang, Binhang Yu and Muhammad Waseem Boota
Land 2025, 14(2), 268; https://doi.org/10.3390/land14020268 - 27 Jan 2025
Viewed by 336
Abstract
It is important to carry out timely scientific assessments of surface subsidence in coal resource cities for ecological environmental protection. Traditional subsidence simulation methods cannot quantitatively describe the driving factors that contribute to or ignore the dynamic connections of subsidence across time and [...] Read more.
It is important to carry out timely scientific assessments of surface subsidence in coal resource cities for ecological environmental protection. Traditional subsidence simulation methods cannot quantitatively describe the driving factors that contribute to or ignore the dynamic connections of subsidence across time and space. Thus, a novel spatio-temporal subsidence simulation model is proposed that couples random forest (RF) and cellular automaton (CA) models, which are used to quantify the contributions of driving factors and simulate the spatio-temporal dynamic changes in subsidence. The RF algorithm is first utilized to clarify the contributions of the driving factors to subsidence and to formulate transformation rules for simulation. Then, a spatio-temporal simulation of subsidence is accomplished by combining it with the CA model. Finally, the method is validated based on the Yongcheng coalfield. The results show that the depth–thickness ratio (0.242), distance to the working face (0.159), distance to buildings (0.150), and lithology (0.147) play main roles in the development of subsidence. Meanwhile, the model can effectively simulate the spatio-temporal changes in mining subsidence. The simulation results were evaluated using 2021 subsidence data as the basis data; the simulation’s overall accuracy (OA) was 0.83, and the Kappa coefficient (KC) was 0.71. This method can obtain a more realistic representation of the spatio-temporal distribution of subsidence while considering the driving factors, which provides technological support for land-use planning and ecological and environmental protection in coal resource cities. Full article
Show Figures

Figure 1

Figure 1
<p>Flow chart of spatio-temporal subsidence simulation.</p>
Full article ">Figure 2
<p>Data processing procedure for InSAR technology.</p>
Full article ">Figure 3
<p>Results after reclassification of foundation subsidence data. (<b>a</b>–<b>c</b>) represent the corresponding subsidence images for 2019–2021, respectively.</p>
Full article ">Figure 4
<p>Geographic location map of the study area. (<b>a</b>) China; (<b>b</b>) Henan Province, China; (<b>c</b>) Yongcheng coalfield.</p>
Full article ">Figure 5
<p>Heat map of Pearson correlation coefficients of driving factors.</p>
Full article ">Figure 6
<p>Contribution of driving factors.</p>
Full article ">Figure 7
<p>Verification of simulation result accuracy. Numbers 1–5 represent no subsidence, slight subsidence, moderate subsidence, more serious subsidence, and serious subsidence, respectively. I, II, III, and IV represent four cross-sections.</p>
Full article ">Figure 8
<p>Analysis of local simulation accuracy in this study area. (<b>a</b>–<b>d</b>) represent the four cross-sections (I, II, III, and IV).</p>
Full article ">Figure 9
<p>Density distribution map of mining subsidence area.</p>
Full article ">Figure 10
<p>Trajectory of change in area of different types of subsidence for the years 2019–2026. Level-1 to Level-5 indicate the different levels of subsidence, respectively.</p>
Full article ">Figure 11
<p>Maps of center of gravity migration in subsidence areas.</p>
Full article ">Figure 12
<p>Comparison of simulation results obtained using different methods. (<b>a1</b>,<b>b1</b>,<b>c1</b>) are the RF-CA, CA, and FLUS model simulation results, respectively, and 1–5 represent no subsidence, slight subsidence, moderate subsidence, more serious subsidence, and serious subsidence, respectively; (<b>a2</b>,<b>b2</b>,<b>c2</b>) are the mis-simulated regions in the simulation results.</p>
Full article ">
23 pages, 24213 KiB  
Article
Optical Image Generation Through Digital Terrain Models for Autonomous Lunar Navigation
by Michele Ceresoli, Stefano Silvestrini and Michèle Lavagna
Aerospace 2025, 12(2), 92; https://doi.org/10.3390/aerospace12020092 - 27 Jan 2025
Viewed by 383
Abstract
In recent years, Vision-Based Navigation (VBN) techniques have emerged as a fundamental component to enable autonomous spacecraft operations, particularly in challenging environments such as planetary landings, where ground control may be limited or unavailable. Developing and testing VBN algorithms requires the availability of [...] Read more.
In recent years, Vision-Based Navigation (VBN) techniques have emerged as a fundamental component to enable autonomous spacecraft operations, particularly in challenging environments such as planetary landings, where ground control may be limited or unavailable. Developing and testing VBN algorithms requires the availability of a large number of realistic images of the application scenario; however, these are rarely available. This paper presents a novel rendering software tool to generate accurate synthetic optical images of the lunar surface by leveraging high-resolution Digital Terrain Models (DTMs). Unlike traditional ray-tracing algorithms, the method iteratively propagates camera rays to determine their intersection with the terrain surface defined by a Digital Elevation Model (DEM). The color information is then retrieved from the corresponding Digital Orthophoto Model (DOM) through the knowledge of the ray impact points, bypassing the need for the costly computation of shadows, reflections, and refractions effects. The rendering performance is demonstrated through a comprehensive selection of images of the lunar surface under different illumination conditions and camera orientations. Full article
(This article belongs to the Special Issue Space Navigation and Control Technologies)
Show Figures

Figure 1

Figure 1
<p>Local region comparison of (<b>a</b>) the 512 pix/degree LOLA-only DEM, (<b>b</b>) SLDEM2015 at 512 pix/degree, and (<b>c</b>) the Chang’e-2 20 m resolution DEM. The images are centered at around 104° W and 31° S and were rendered by using Telespazio’s Interactive Mission Modeling, Visualization &amp; Validation Tool (IMMV<sup>2</sup>) [<a href="#B16-aerospace-12-00092" class="html-bibr">16</a>,<a href="#B33-aerospace-12-00092" class="html-bibr">33</a>].</p>
Full article ">Figure 2
<p>(<b>a</b>) Northern and (<b>b</b>) southern polar stereographic NAC mosaics. Both images manifest local lighting inconsistencies, in the form of unrealistic shadows due to regions of near-permanent darkness or as discontinuities along the extended NAC radial strips.</p>
Full article ">Figure 3
<p>(<b>a</b>,<b>b</b>) Original and (<b>c</b>,<b>d</b>) corrected Chang’e-2 north pole DOM sections. Craters from the Robbins catalog [<a href="#B41-aerospace-12-00092" class="html-bibr">41</a>] with dimensions between 2.5 and 50 km have been projected onto the images as a reference. A comparison with the same projection on a portion of the global WAC 100 m/pixel mosaic is shown in (<b>e</b>,<b>f</b>). Only craters that are not completely shadowed are highlighted.</p>
Full article ">Figure 4
<p>Pinhole camera model. The image highlights the camera coordinate system (<span class="html-italic">left</span>) and the Field-of-View (FoV) influence on the image plane coordinates (<span class="html-italic">right</span>).</p>
Full article ">Figure 5
<p>Pixel coordinate conversions. The coordinates of a point are first expressed in raster space, then converted to NDC space (i.e., mapped to a range between 0 and 1), and finally converted to screen space.</p>
Full article ">Figure 6
<p>Comparison between a wrong (<span class="html-italic">left</span>) and a correct (<span class="html-italic">right</span>) ray propagation order. Black and red arrows indicate correct and miscalculated impact locations, respectively. In the left image, as the propagation is started from one of the distant impact points, the successive rays begin the propagation at a <span class="html-italic">t</span>-value that causes the algorithm to miss closer obstacles. The error is removed by starting the propagation from the ray closest to the surface.</p>
Full article ">Figure 7
<p>Illustration of the basic principles behind SSAA.</p>
Full article ">Figure 8
<p>Scene rendering without (<span class="html-italic">top</span>) and with (<span class="html-italic">bottom</span>) SSAA. The bottom scene is rendered by using 8 subsamples per pixel.</p>
Full article ">Figure 9
<p>Rendering examples of a local region under different levels of defocus blur, decreasing from left to right.</p>
Full article ">Figure 10
<p>Rendering examples at different resolution levels.</p>
Full article ">Figure 11
<p>Rendering examples of large impact craters from the northern lunar hemisphere.</p>
Full article ">Figure 12
<p>Projections of crater rims from the Robbins catalogue with diameters between 7 and 500 km superimposed on renderings of the Schomberger-C crater.</p>
Full article ">Figure 13
<p>Off-nadir camera pointings highlighting the central structures of the (<b>a</b>) Lovelace and (<b>b</b>) Pythagoras craters.</p>
Full article ">
21 pages, 1216 KiB  
Article
Long-Term Effects of Crop Treatments and Fertilization on Soil Stability and Nutrient Dynamics in the Loess Plateau: Implications for Soil Health and Productivity
by Farhat Ullah Khan, Faisal Zaman, Yuanyuan Qu, Junfeng Wang, Ojimamdov Habib Darmorakhtievich, Qinxuan Wu, Shah Fahad, Feng Du and Xuexuan Xu
Sustainability 2025, 17(3), 1014; https://doi.org/10.3390/su17031014 - 26 Jan 2025
Viewed by 485
Abstract
Soil degradation and erosion pose significant threats to agricultural sustainability in fragile ecosystems, such as the Loess Plateau in northern China. This study examines the long-term impacts of fertilization regimes and land-use systems on soil health, focusing on soil aggregate stability, fertility, and [...] Read more.
Soil degradation and erosion pose significant threats to agricultural sustainability in fragile ecosystems, such as the Loess Plateau in northern China. This study examines the long-term impacts of fertilization regimes and land-use systems on soil health, focusing on soil aggregate stability, fertility, and crop productivity. Six treatment combinations were evaluated in our study, including three continuous alfalfa fields (AL-CK, AL-P, and AL-NPM) and three continuous wheat fields (WH-NPM, WH-NP, and WH-P), each representing a combination of land use and three fertilization treatments: (1) no fertilization (CK), (2) inorganic fertilization (120 kg ha−1 N, 60 kg ha−1 P-NP), and (3) a combination of organic and inorganic fertilization (75 t ha−1 cow manure-NPM). Soil samples were collected from three depths (0–10 cm, 10–20 cm, and 20–30 cm) to assess physical and chemical properties. We evaluated the long-term effects of different fertilization treatments on soil stability, fertility, and crop yield to explore the interactions among soil’s physical and chemical properties under two land-use types and to assess the effectiveness of combined organic and inorganic fertilization strategies in improving soil health and mitigating erosion in vulnerable landscapes. The study revealed significant depth-specific variations with surface layers (0–10 cm) showing the greatest improvement under NPM treatments, particularly in continuous alfalfa fields, which exhibited higher soil fertility, improved soil structure, and crop yield. In contrast, continuous wheat fields with minimal fertilization demonstrated significantly lower soil quality and productivity. Using the combination of mineral fertilizers and organic amendments, such as cow manure, proved to be the most effective strategy for significantly enhancing nutrient availability and overall soil health. Partial Least Squares Modeling (PLS-M) and Mantel analysis highlighted the critical role of fertilization management in maintaining soil quality, boosting crop productivity, and mitigating erosion in high-risk areas. This study emphasizes the importance of integrated nutrient management for sustainable land use and soil conservation in erosion-prone regions. Full article
19 pages, 2367 KiB  
Article
Determining Sources of Air Pollution Exposure Inequity in New York City Through Land-Use Regression Modeling of PM2.5 Constituents
by Masha Pitiranggon, Sarah Johnson, Ariel Spira-Cohen, Holger Eisl and Kazuhiko Ito
Pollutants 2025, 5(1), 2; https://doi.org/10.3390/pollutants5010002 - 26 Jan 2025
Viewed by 275
Abstract
Differences in exposures and resources to manage personal health contribute to persistent inequities in air pollution burden despite vast air quality improvements over the past 2–3 decades in the United States. These factors are, partly, linked to historic racist practices, such as redlining, [...] Read more.
Differences in exposures and resources to manage personal health contribute to persistent inequities in air pollution burden despite vast air quality improvements over the past 2–3 decades in the United States. These factors are, partly, linked to historic racist practices, such as redlining, a discriminatory housing policy that was practiced legally between 1935 and 1968. Using 100 m × 100 m resolution land-use regression predicted surfaces of PM2.5 constituents (black carbon, nickel, vanadium, and copper) as pollution source indicators, we fit Bayesian generalized linear mixed-effects models to examine differences in source exposures over two study periods, 2008–2015 and 2016–2019, comparing (1) redlined to not redlined and (2) high-asthma to low-asthma neighborhoods. We examine redlining as an indicator of historical, and structural racism and asthma rates as an indicator of present-day community burden. Redlined areas saw near elimination of disparities in exposure to residual oil boilers and marine residual oil but persistent disparities in traffic. High-asthma neighborhoods continue to have disproportionately high exposures to both residual oil boilers and traffic, with no discernable disparities related to marine residual oil emissions. Overall exposure disparities are small, with PM2.5 disparities by both asthma morbidity and redlining amounting to less than 1 µg/m3 and NO2 disparities by asthma and redlining amounting to less than 2 ppb in the post-2016 period. For context, 2019 NYC average PM2.5 and NO2 were 8.5 µg/m3 and 20 ppb, respectively. Our findings suggest that local pollution policy should focus on reducing traffic and building boiler emissions in high-asthma neighborhoods to reduce exacerbations. Full article
(This article belongs to the Section Air Pollution)
25 pages, 3692 KiB  
Article
Physical Parameterization Sensitivity of Noah-MP for Hydrothermal Simulation Within the Active Layer on the Qinghai–Tibet Plateau
by Yongliang Jiao, Ren Li, Tonghua Wu, Xiaodong Wu, Shenning Wang, Jimin Yao, Guojie Hu, Xiaofan Zhu, Jianzong Shi, Yao Xiao, Erji Du and Yongping Qiao
Land 2025, 14(2), 247; https://doi.org/10.3390/land14020247 - 24 Jan 2025
Viewed by 297
Abstract
The accurate modeling of complex freeze–thaw processes and hydrothermal dynamics within the active layer is challenging. Due to the uncertainty in hydrothermal simulation, it is necessary to thoroughly investigate the parameterization schemes in land surface models. The Noah-MP was utilized in this study [...] Read more.
The accurate modeling of complex freeze–thaw processes and hydrothermal dynamics within the active layer is challenging. Due to the uncertainty in hydrothermal simulation, it is necessary to thoroughly investigate the parameterization schemes in land surface models. The Noah-MP was utilized in this study to conduct 23,040 ensemble experiments based on 11 physical processes, which were aimed at improving the understanding of parameterization schemes and reducing model uncertainty. Next, the impacts of uncertainty of physical processes on land surface modeling were evaluated via Natural Selection and Tukey’s test. Finally, Technique for Order of Preference by Similarity to Ideal Solution (TOPSIS) was used to identify the optimal combination of parameterization schemes for improving hydrothermal simulation. The results of Tukey’s test agreed well with those of Natural Selection for most soil layers. More importantly, Tukey’s test identified more parameterization schemes with consistent model performance for both soil temperature and moisture. Results from TOPSIS showed that the determination of optimal schemes was consistent for the simulation of soil temperature and moisture in each physical process except for frozen soil permeability (INF). Further analysis showed that scheme 2 of INF yielded better simulation results than scheme 1. The improvement of the optimal scheme combination during the frozen period was more significant than that during the thawed period. Full article
(This article belongs to the Section Land – Observation and Monitoring)
18 pages, 6072 KiB  
Article
Application of UAV Photogrammetry and Multispectral Image Analysis for Identifying Land Use and Vegetation Cover Succession in Former Mining Areas
by Volker Reinprecht and Daniel Scott Kieffer
Remote Sens. 2025, 17(3), 405; https://doi.org/10.3390/rs17030405 - 24 Jan 2025
Viewed by 440
Abstract
Variations in vegetation indices derived from multispectral images and digital terrain models from satellite imagery have been successfully used for reclamation and hazard management in former mining areas. However, low spatial resolution and the lack of sufficiently detailed information on surface morphology have [...] Read more.
Variations in vegetation indices derived from multispectral images and digital terrain models from satellite imagery have been successfully used for reclamation and hazard management in former mining areas. However, low spatial resolution and the lack of sufficiently detailed information on surface morphology have restricted such studies to large sites. This study investigates the application of small, unmanned aerial vehicles (UAVs) equipped with multispectral sensors for land cover classification and vegetation monitoring. The application of UAVs bridges the gap between large-scale satellite remote sensing techniques and terrestrial surveys. Photogrammetric terrain models and orthoimages (RGB and multispectral) obtained from repeated mapping flights between November 2023 and May 2024 were combined with an ALS-based reference terrain model for object-based image classification. The collected data enabled differentiation between natural forests and areas affected by former mining activities, as well as the identification of variations in vegetation density and growth rates on former mining areas. The results confirm that small UAVs provide a versatile and efficient platform for classifying and monitoring mining areas and forested landslides. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>A</b>) Overview of the study site (“Trassbruch Gossendorf”) based on the digital elevation model; (<b>B</b>) oblique photograph. Former mining and mine dump areas, access roads and the landslide area are highlighted in (<b>A</b>).</p>
Full article ">Figure 2
<p>(<b>A</b>) Study site with the boundaries of former mining, mine dump and landslide affected areas. (<b>B</b>) Subset at the southern slope, visualizing the segmentation and the effect of the 0.5 m buffer around the sampling points and the typical tree crown dimension (diameter ~2–3 m).</p>
Full article ">Figure 3
<p>Python-based OBIA workflow, including a summary of each processing step.</p>
Full article ">Figure 4
<p>Classified map datasets for all four classification periods. (<b>A</b>) November 2023 (sunny, oblique flight); (<b>B</b>) December 2023 (overcast, nadir flight); (<b>C</b>) April 2024 (overcast, nadir flight); (<b>D</b>) May 2024 (sunny, nadir flight). [X] = area prone to misclassification (Zone A2), [Y] = old mine dump (Zone B1), that was only partially cleared for operation.</p>
Full article ">Figure 5
<p>(<b>A</b>) Parameter variation during the cross-validation process (global performance metrics and class performance metrics). (<b>B</b>) Classification metrics for all flight epochs including combined confusion matrices. (<b>C</b>) Confusion matrices derived from holdout dataset (holdout confusion matrix). The confusion matrices were standardized in horizontal direction and the corresponding sample number is given in square brackets.</p>
Full article ">Figure 6
<p>Time series for the mean NDVI, NDRE, height above rDTM (dDTM), height above rDSM and (dDSM) extracted from the former mining zones (mine dump, mine), the landslide area and the natural forest.</p>
Full article ">
24 pages, 7033 KiB  
Article
geeSSEBI: Evaluating Actual Evapotranspiration Estimated with a Google Earth Engine Implementation of S-SEBI
by Jerzy Piotr Kabala, Jose Antonio Sobrino, Virginia Crisafulli, Dražen Skoković and Giovanna Battipaglia
Remote Sens. 2025, 17(3), 395; https://doi.org/10.3390/rs17030395 - 24 Jan 2025
Viewed by 385
Abstract
Quantifying and mapping evapotranspiration (ET) from land surfaces is crucial in the context of climate change. For decades, remote sensing data have been utilized to estimate ET, leading to the development of numerous algorithms. However, their application is still non-trivial, mainly due to [...] Read more.
Quantifying and mapping evapotranspiration (ET) from land surfaces is crucial in the context of climate change. For decades, remote sensing data have been utilized to estimate ET, leading to the development of numerous algorithms. However, their application is still non-trivial, mainly due to practical constraints. This paper introduces geeSSEBI, a Google Earth Engine implementation of the S-SEBI (Simplified Surface Energy Balance Index) model, for deriving ET from Landsat data and ERA5-land radiation. The source code and a graphical user interface implemented as a Google Earth Engine application are provided. The model ran on 871 images, and the estimates were evaluated against multiyear data of four eddy covariance towers belonging to the ICOS network, representative of both forests and agricultural landscapes. The model showed an RMSE of approximately 1 mm/day, and a significant correlation with the observed values, at all the sites. A procedure to upscale the data to monthly is proposed and tested as well, and its accuracy evaluated. Overall, the model showed acceptable accuracy, while performing better on forest ecosystems than on agricultural ones, especially at daily and monthly timescales. This implementation is particularly valuable for mapping evapotranspiration in data-scarce environments by utilizing Landsat archives and ERA5-land radiation estimates. Full article
(This article belongs to the Special Issue Remote Sensing and Modelling of Terrestrial Ecosystems Functioning)
Show Figures

Figure 1

Figure 1
<p>Example of the S-SEBI diagram, showing the relationship between LST and albedo, and the wet and dry limits, employed in the evaporative fraction computation, represented by the red and yellow regression lines, respectively. Blue dots represent several sample pixels of the area of interest.</p>
Full article ">Figure 2
<p>Graphical explanation of the monthly scaling procedure.</p>
Full article ">Figure 3
<p>Workflow of the Google Earth Engine implementation of the S-SEBI evapotranspiration algorithm presented in this paper.</p>
Full article ">Figure 4
<p>On the left, RGB composites of Landsat 8 images, depicting the areas where the model was applied. The red circle corresponds to the location of the eddy covariance tower. On the right, climatograms for each study area based on the ERA5-land data. The blue bars represent monthly precipitation, the red line the average monthly temperature, with the red shade identifying the range between minimum and maximum temperatures. Climatic data of the last 30 years are represented.</p>
Full article ">Figure 5
<p>Modelled LE values, against eddy covariance, corrected, half-hourly LE values. (<b>A</b>) SanRossore2 site, (<b>B</b>) Bosco Fontana site, (<b>C</b>) Lison site, and (<b>D</b>) Lamasquere site. In blue, the linear regression of the modelled against measured values is shown, with grey shadow representing its standard error. The line equation and correlation coefficient are reported in blue as well. The red dashed line is the identity line.</p>
Full article ">Figure 6
<p>Scatterplot of the daily estimates obtained with the S-SEBI model, against the daily average evapotranspiration values recorded at the eddy covariance towers. In blue, the regression line between the modelled and observed values is shown, along with the equation of the regression line, and Pearson correlation coefficient (R). The grey shadows represent the standard errors of the regression line. The red line is the identity line. (<b>A</b>) SanRossore2 site, (<b>B</b>) Bosco Fontana site, (<b>C</b>) Lison site, (<b>D</b>) Lamasquere site.</p>
Full article ">Figure 7
<p>Scatterplot comparing the modelled values against the observed ones, for non-overpass dates. Site names abbreviated as: BF (Bosco Fontana, magenta circles), LM (Lamasquere, blue circles), LS (Lison, yellow circles), and SR (San Rossore2, dark green circles). The identity line is shown as dashed black line.</p>
Full article ">Figure 8
<p>Comparison of the monthly scaled ET values obtained with S-SEBI against the eddy covariance data. (<b>A</b>) San Rossore2 site, (<b>B</b>) Bosco Fontana site, (<b>C</b>) Lison site, (<b>D</b>) Lamasquere site. Blue lines represent the linear regression of modelled against observed data, along with the standard error (grey shadow), line equation and correlation coefficient (written in blue). Black dots represent the observations, while the yellow dashed line is the identity line.</p>
Full article ">Figure 9
<p>Evaluation of LE and ET estimated with S-SEBI for the Lamasquere site on Landsat 5 images acquired between 2005 and 2011. The blue line represents the linear regression between the S-SEBI estimates and the eddy covariance measurements, along with the standard error (grey shadow). Black dots represent actual observations, while the red lines (solid and dashed) represent the identity line. (<b>A</b>) Instantaneous latent heat of evapotranspiration and (<b>B</b>) daily evapotranspiration.</p>
Full article ">Figure 10
<p>Comparison of the accuracy metrics of S-SEBI and MOD16A2 evapotranspiration product at the four eddy covariance sites encompassed in this study. Comparison was performed at an 8-day timescale, thus the unit of RMSE is mm/8 days, % bias is expressed as %, and of all the other metrics (R, R squared, NSE, and KGE) are dimensionless. Site names abbreviated as BF (Bosco Fontana), SR (San Rossore2), LS (Lison), and LM (Lamasquere).</p>
Full article ">
25 pages, 12496 KiB  
Article
Impact of Future Climate Change on Groundwater Storage in China’s Large Granary: A Study Based on LSTM and CMIP6 Models
by Haiqing Wang, Peng Qi, Moran Xu, Yao Wu and Guangxin Zhang
Water 2025, 17(3), 315; https://doi.org/10.3390/w17030315 - 23 Jan 2025
Viewed by 396
Abstract
Northeast China, as a primary grain-producing region, has long drawn attention for its intensive groundwater extraction for irrigation. However, previous studies on the future spatiotemporal changes of groundwater storage (GWS) are lacking. Utilizing the Global Land Data Assimilation System Version 2.2 (GLDAS-2.2), which [...] Read more.
Northeast China, as a primary grain-producing region, has long drawn attention for its intensive groundwater extraction for irrigation. However, previous studies on the future spatiotemporal changes of groundwater storage (GWS) are lacking. Utilizing the Global Land Data Assimilation System Version 2.2 (GLDAS-2.2), which simulates groundwater storage (as Equivalent Water Height) using the Catchment Land Surface Model (CLSM-F2.5) and calibrates it with terrestrial water storage data from the GRACE satellite, we analyzed the spatiotemporal variations of GWS in northeast China and employed a Long Short-Term Memory (LSTM) neural network model to quantify the responses of GWS to future climate change. Maintaining current socio–economic factors and combining climate factors from four scenarios (SSP126, SSP245, SSP370, and SSP585) under the CMIP6 model, we predicted GWS from 2022 to 2100. The results indicate that historically, groundwater storage exhibits a decreasing trend in the south and an increasing trend in the north, with a 44° N latitude boundary. Under the four scenarios, the predicted GWS increments in northeast China are 0.08 ± 0.09 mm/yr in SSP126, 0.11 ± 0.08 mm/yr in SSP245, 0.12 ± 0.09 mm/yr in SSP370, and 0.20 ± 0.07 mm/yr in SSP585. Although overall groundwater storage has slightly increased and the model projections indicate a continued increase, the southern part of the region may not return to past levels and faces water stress risks. This study provides an important reference for the development of sustainable groundwater management strategies. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Spatial distribution of mean GWS (as Equivalent Water Height, representing the spatial pattern of groundwater storage) from February 2003 to December 2022 based on GLDAS-2.2 [<a href="#B29-water-17-00315" class="html-bibr">29</a>].</p>
Full article ">Figure 2
<p>Comparison of GLDAS and CSR GRACE GWSA Data for northeast China, with a NSE of 0.758.</p>
Full article ">Figure 3
<p>Plot of correlation analysis between current and corrected data for climate influences from 1950 to 2014, from left to right, for precipitation, potential evapotranspiration, and temperature.</p>
Full article ">Figure 4
<p>LSTM network structure.</p>
Full article ">Figure 5
<p>Spatial distributions of GWS interannual spacing values in the northeast China, 2004–2022.</p>
Full article ">Figure 6
<p>Spatial trend of the GWS in northeast China, 2003.02–2022.12. (<b>a</b>) The spatial trend of the GWS in northeast China, (<b>b</b>) spatial distribution of <span class="html-italic">p</span>-value corresponding to the trend.</p>
Full article ">Figure 7
<p>Spatial distribution of monthly-scale groundwater storage anomaly (GWSA) in northeast China.</p>
Full article ">Figure 8
<p>Spatial distributions of monthly trends in GWS in northeast China.</p>
Full article ">Figure 9
<p>Spatial distribution of correlation coefficients between GWS and corresponding influencing factors in northeast China (inter-annual data), with no arable land in the Daxinganling region.</p>
Full article ">Figure 10
<p>Spatial distribution of LSTM model fit scores. The data from 2003.02 to 2016.09 were used as the training set, and the data from 2016.10 to 2019.12 were used as the test set. (<b>a</b>,<b>b</b>) were the NSE and RMSE spatial distribution plots of the model test set and the status quo data, respectively.</p>
Full article ">Figure 11
<p>Spatial distribution map of the future GWS. (<b>a</b>) Spatial trend of GWS in different contexts over time and its plotting. (<b>b</b>) Spatial variation of the mean GWS in different contexts over time and its plotting against the mean GWS in the historical period.</p>
Full article ">Figure 12
<p>Temporal changes in GWS in the northeast, with shading indicating the range of fluctuations in GWS under different scenarios. 2004–2019 is a plot of interannual changes in GWS for the historical period, and 2002–2100 is a plot of simulated interannual changes in GWS for the future.</p>
Full article ">Figure 13
<p>Box plots of monthly trends in mean values of future scenarios of groundwater storage in the northeast and changes in GWS. (<b>a</b>,<b>b</b>) Changes in inter-monthly trends for different periods. (<b>c</b>,<b>d</b>) Comparisons of inter-monthly mean values of GWS for different periods with historical GWS. (<b>a</b>,<b>c</b>) The period in 2022–2060; (<b>b</b>,<b>d</b>) the period in 2061–2100.</p>
Full article ">Figure 14
<p>Fitted plots for the four sampled data points. In the figure, the yellow line represents the historical data, the blue line represents the prediction result of the training set, and the green line represents the prediction result of the testing set. (<b>a</b>,<b>b</b>) LSTM model; (<b>c</b>,<b>d</b>) Conv-LSTM model.</p>
Full article ">Figure 15
<p>Schematic diagram of LSTM model.</p>
Full article ">Figure 16
<p>Graph of temporal trends in future climate factors.</p>
Full article ">
25 pages, 6944 KiB  
Article
Representation Learning of Multi-Spectral Earth Observation Time Series and Evaluation for Crop Type Classification
by Andrea González-Ramírez, Clement Atzberger, Deni Torres-Roman and Josué López
Remote Sens. 2025, 17(3), 378; https://doi.org/10.3390/rs17030378 - 23 Jan 2025
Viewed by 410
Abstract
Remote sensing (RS) spectral time series provide a substantial source of information for the regular and cost-efficient monitoring of the Earth’s surface. Important monitoring tasks include land use and land cover classification, change detection, forest monitoring and crop type identification, among others. To [...] Read more.
Remote sensing (RS) spectral time series provide a substantial source of information for the regular and cost-efficient monitoring of the Earth’s surface. Important monitoring tasks include land use and land cover classification, change detection, forest monitoring and crop type identification, among others. To develop accurate solutions for RS-based applications, often supervised shallow/deep learning algorithms are used. However, such approaches usually require fixed-length inputs and large labeled datasets. Unfortunately, RS images acquired by optical sensors are frequently degraded by aerosol contamination, clouds and cloud shadows, resulting in missing observations and irregular observation patterns. To address these issues, efforts have been made to implement frameworks that generate meaningful representations from the irregularly sampled data streams and alleviate the deficiencies of the data sources and supervised algorithms. Here, we propose a conceptually and computationally simple representation learning (RL) approach based on autoencoders (AEs) to generate discriminative features for crop type classification. The proposed methodology includes a set of single-layer AEs with a very limited number of neurons, each one trained with the mono-temporal spectral features of a small set of samples belonging to a class, resulting in a model capable of processing very large areas in a short computational time. Importantly, the developed approach remains flexible with respect to the availability of clear temporal observations. The signal derived from the ensemble of AEs is the reconstruction difference vector between input samples and their corresponding estimations, which are averaged over all cloud-/shadow-free temporal observations of a pixel location. This averaged reconstruction difference vector is the base for the representations and the subsequent classification. Experimental results show that the proposed extremely light-weight architecture indeed generates separable features for competitive performances in crop type classification, as distance metrics scores achieved with the derived representations significantly outperform those obtained with the initial data. Conventional classification models were trained and tested with representations generated from a widely used Sentinel-2 multi-spectral multi-temporal dataset, BreizhCrops. Our method achieved 77.06% overall accuracy, which is 6% higher than that achieved using original Sentinel-2 data within conventional classifiers and even 4% better than complex deep models such as OmnisCNN. Compared to extremely complex and time-consuming models such as Transformer and long short-term memory (LSTM), only a 3% reduction in overall accuracy was noted. Our method uses only 6.8k parameters, i.e., 400x fewer than OmnicsCNN and 27x fewer than Transformer. The results prove that our method is competitive in terms of classification performance compared with state-of-the-art methods while substantially reducing the computational load. Full article
(This article belongs to the Collection Sentinel-2: Science and Applications)
Show Figures

Figure 1

Figure 1
<p>Illustration of representation learning (RL) as a function <span class="html-italic">f</span>, mapping vectors from a dimensional space to a representation space.</p>
Full article ">Figure 2
<p>Example of an autoencoder architecture with mathematical definition as a function. In the present work, the reconstruction difference between the input and output is used as a representation and not the code itself.</p>
Full article ">Figure 3
<p>First level of the proposed workflow. A scene classification product provided by the European Space Agency (ESA) is used to mask out cloudy samples from a geographic point (pixel) shaped as a <math display="inline"><semantics> <mrow> <mi>T</mi> <mo>×</mo> <mi>B</mi> </mrow> </semantics></math> array.</p>
Full article ">Figure 4
<p>Proposed framework block diagram. The full methodology is composed of four main blocks: data preprocessing, model training, representation generation and evaluation.</p>
Full article ">Figure 5
<p>Dataset downloading process using the Google Earth Engine (GEE) database.</p>
Full article ">Figure 6
<p>Example of the expected output for positive and negative samples. The difference from the ensemble of autoencoders (AEs) constitutes the representations for the downstream task.</p>
Full article ">Figure 7
<p>Autoencoder (AE) training. Each autoencoder is trained with a finite set of individual spectral curves belonging to one of the crop types. The reconstructions from the <span class="html-italic">C</span> classes are used to calculate the difference vector across the ensemble that is the final set of representations.</p>
Full article ">Figure 8
<p>Inference workflow of the proposed framework. For each temporal set of cloud-free reflectance spectra, the average reconstruction difference vector is calculated for each of the <span class="html-italic">C</span> autoencoders (AEs) and concatenated to define the representations of this pixel.</p>
Full article ">Figure 9
<p>3D scatterplot of (<b>a</b>) S2 fixed-length time series (45 observations) and (<b>b</b>) representation over three principal components obtained by t-distributed Stochastic Neighbor Embedding (TSNE) only for visual interpretation.</p>
Full article ">Figure 10
<p>Overall accuracy (OA) of the random forest (RF), support vector machine (SVM), extreme gradient boosting (XGBoost) and fully connected network (FCN) trained with a variable percentage of training samples and using (i) representations (solid line) and (ii) original Sentinel-2 data (broken line).</p>
Full article ">Figure 11
<p>(<b>a</b>) True color image of the study area in 2017 and composites images generated by combining three random representations per map: (<b>b</b>) 9-64-30, (<b>c</b>) 59-84-81, (<b>d</b>) 30-11-141, (<b>e</b>) 45-66-57, (<b>f</b>) 20-10-32, (<b>g</b>) 5-142-83 and (<b>h</b>) 24-79-133.</p>
Full article ">Figure 12
<p>(<b>a</b>) Study area ground truth at field level (polygons), (<b>b</b>) representations-based fully connected network (FCN) pixel-wise classification (raster), (<b>c</b>) representations-based FCN field-based classification (polygons) and (<b>d</b>) map of correctly classified fields in green and misclassified fields in red.</p>
Full article ">Figure 12 Cont.
<p>(<b>a</b>) Study area ground truth at field level (polygons), (<b>b</b>) representations-based fully connected network (FCN) pixel-wise classification (raster), (<b>c</b>) representations-based FCN field-based classification (polygons) and (<b>d</b>) map of correctly classified fields in green and misclassified fields in red.</p>
Full article ">Figure A1
<p>Hyperparameters and quality indicators correlation matrix.</p>
Full article ">
Back to TopTop