Nothing Special   »   [go: up one dir, main page]

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

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (839)

Search Parameters:
Keywords = geostationary

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
35 pages, 27811 KiB  
Article
Machine Learning to Retrieve Gap-Free Land Surface Temperature from Infrared Atmospheric Sounding Interferometer Observations
by Fabio Della Rocca, Pamela Pasquariello, Guido Masiello, Carmine Serio and Italia De Feis
Remote Sens. 2025, 17(4), 694; https://doi.org/10.3390/rs17040694 - 18 Feb 2025
Abstract
Retrieving LST from infrared spectral observations is challenging because it needs separation from emissivity in surface radiation emission, which is feasible only when the state of the surface–atmosphere system is known. Thanks to its high spectral resolution, the Infrared Atmospheric Sounding Interferometer (IASI) [...] Read more.
Retrieving LST from infrared spectral observations is challenging because it needs separation from emissivity in surface radiation emission, which is feasible only when the state of the surface–atmosphere system is known. Thanks to its high spectral resolution, the Infrared Atmospheric Sounding Interferometer (IASI) instrument onboard Metop polar-orbiting satellites is the only sensor that can simultaneously retrieve LST, the emissivity spectrum, and atmospheric composition. Still, it cannot penetrate thick cloud layers, making observations blind to surface emissions under cloudy conditions, with surface and atmospheric parameters being flagged as voids. The present paper aims to discuss a downscaling–fusion methodology to retrieve LST missing values on a spatial field retrieved from spatially scattered IASI observations to yield level 3, regularly gridded data, using as proxy data LST from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) flying on Meteosat Second Generation (MSG) platform, a geostationary instrument, and from the Advanced Very High-Resolution Radiometer (AVHRR) onboard Metop polar-orbiting satellites. We address this problem by using machine learning techniques, i.e., Gradient Boosting, Random Forest, Gaussian Process Regression, Neural Network, and Stacked Regression. We applied the methodology over the Po Valley region, a very heterogeneous area that allows addressing the trained models’ robustness. Overall, the methods significantly enhanced spatial sampling, keeping errors in terms of Root Mean Square Error (RMSE) and bias (Mean Absolute Error, MAE) very low. Although we demonstrate and assess the results primarily using IASI data, the paper is also intended for applications to the IASI follow-on, that is, IASI Next Generation (IASI-NG), and much more to the Infrared Sounder (IRS), which is planned to fly this year, 2025, on the Meteosat Third Generation platform (MTG). Full article
(This article belongs to the Special Issue Remote Sensing in Geomatics (Second Edition))
Show Figures

Figure 1

Figure 1
<p>The red box indicates the Po Valley target region with the CLC 2018 as shapefile.</p>
Full article ">Figure 2
<p>Flowchart of the proposed framework. (<b>a</b>) Retrieval of LST; (<b>b</b>) Training; (<b>c</b>) L3 LST.</p>
Full article ">Figure 3
<p>Comparison of L2 IASI observations and the derived prediction mask for August 2022. The left panel shows the spatial distribution of L2 observations across the 9 years of data, while the right panel shows the spatial domain used for prediction.</p>
Full article ">Figure 4
<p>Comparison of MAE cross validated errors of the tested ML algorithms: Random Forest (blue), Boosting (orange), Neural Network (yellow), Gaussian Process Regression (purple), and Stacked Regression (green). The numbers in the legend represent the average MAE for all methods calculated across all months.</p>
Full article ">Figure 5
<p>Comparison of RMSE cross validated errors of the tested ML algorithms: Random Forest (blue), Boosting (orange), Neural Network (yellow), Gaussian Process Regression (purple), and Stacked Regression (green). The numbers in the legend represent the average RMSE for all methods calculated across all months.</p>
Full article ">Figure 6
<p>Example of the L3 LST for the months January–June. The first column represents the IASI LST L2 observations, while the second column shows the LST L3 predicted with Stacked Regression; each row represents a month.</p>
Full article ">Figure 7
<p>Example of the L3 LST for the months July–December. The first column represents the IASI LST L2 observations, while the second column shows the LST L3 predicted with Stacked Regression; each row represents a month.</p>
Full article ">Figure 8
<p>Comparison between the predicted LST for August 2022 and the mean of nine years of AVHRR and SEVIRI data for the same month. The top-left panel shows the IASI L2 observations, while the right panel displays the difference maps with SEVIRI (<b>top</b>) and AVHRR (<b>bottom</b>) including also the L2 observations represented by the small black dots. The bottom-left panel presents the KDE plot of these differences, including the mean and standard deviation of the errors.</p>
Full article ">Figure 9
<p>Comparison between the predicted LST for March 2022 and the mean of nine years of AVHRR and SEVIRI data for the same month. The top-left panel shows the IASI L2 observations, while the right panel displays the difference maps with SEVIRI (<b>top</b>) and AVHRR (<b>bottom</b>) including also the L2 observations represented by the small black dots. The bottom-left panel presents the KDE plot of these differences, including the mean and standard deviation of the errors.</p>
Full article ">Figure 10
<p>L2/L3 Differences IASI - MODIS for the months January–June. The first column displays the L2/L3 differences using KDE plots with the mean and standard deviation: the red curves display the L2 errors, and the blue curves display the L3 errors. The second column shows the scatterplots between the predicted IASI L3 LST values and the MODIS L3 LST values, the linear fits, and the <math display="inline"><semantics> <mrow> <msup> <mi mathvariant="normal">R</mi> <mn>2</mn> </msup> </mrow> </semantics></math> indexes.</p>
Full article ">Figure 11
<p>L2/L3 Differences IASI - MODIS for the months July–December. The first column displays the L2/L3 differences using KDE plots with the mean and standard deviation: the red curves display the L2 errors, and the blue curves display the L3 errors. The second column shows the scatterplots between the predicted IASI L3 LST values and the MODIS L3 LST values, the linear fits, and the <math display="inline"><semantics> <mrow> <msup> <mi mathvariant="normal">R</mi> <mn>2</mn> </msup> </mrow> </semantics></math> indexes.</p>
Full article ">Figure 12
<p>Comparison LSTs from IASI, AVHRR, SEVIRI and MODIS for the months January–June. The first column displays the L2 differences using KDE plots with the mean and standard deviation: the blue curves represent the L2 differences between IASI and AVHRR, the red curves represent the L2 differences between IASI and SEVIRI, and the yellow curves represent the L2 differences between IASI and MODIS. The second column shows the same differences using boxplots, with the <math display="inline"><semantics> <mrow> <msup> <mi mathvariant="normal">R</mi> <mn>2</mn> </msup> </mrow> </semantics></math> indexes included on the <span class="html-italic">x</span>-axis.</p>
Full article ">Figure 13
<p>Comparison LSTs from IASI, AVHRR, SEVIRI and MODIS for the months July–December. The first column displays the L2 differences using KDE plots with the mean and standard deviation: the blue curves represent the L2 differences between IASI and AVHRR, the red curves represent the L2 differences between IASI and SEVIRI, and the yellow curves represent the L2 differences between IASI and MODIS. The second column shows the same differences using boxplots, with the <math display="inline"><semantics> <mrow> <msup> <mi mathvariant="normal">R</mi> <mn>2</mn> </msup> </mrow> </semantics></math> indexes included on the <span class="html-italic">x</span>-axis.</p>
Full article ">
23 pages, 10437 KiB  
Article
A Continuous B2b-PPP Model Considering Interruptions in BDS-3 B2b Orbits and Clock Corrections as Well as Signal-in-Space Range Error Residuals
by Rui Shang, Zhenhao Xu, Chengfa Gao, Xiaolin Meng, Wang Gao and Qi Liu
Remote Sens. 2025, 17(4), 618; https://doi.org/10.3390/rs17040618 - 11 Feb 2025
Viewed by 256
Abstract
In 2020, BDS-3 began broadcasting high-precision positioning correction products through B2b signals, effectively addressing the limitations of ground-based augmentation. However, challenges such as the “south wall effect” from geostationary orbit (GEO) satellites, issues of data (IOD) mismatch, and signal priority conflicts often result [...] Read more.
In 2020, BDS-3 began broadcasting high-precision positioning correction products through B2b signals, effectively addressing the limitations of ground-based augmentation. However, challenges such as the “south wall effect” from geostationary orbit (GEO) satellites, issues of data (IOD) mismatch, and signal priority conflicts often result in interruptions and anomalies during real-time positioning with the B2b service. This paper proposes a continuous B2b-PPP (B2b signal-based Precise Point Positioning) model that incorporates signal-in-space range error (SISRE) residuals and predictions for B2b orbits and clock corrections to achieve seamless, high-precision continuous positioning. In our experiments, we first analyze the characteristics of B2b SISRE for both BDS-3 and GPS. We then evaluate the positioning accuracy of several models, B2b-PPP, EB2b-PPP, PB2b-PPP, EB2bS-PPP, and PB2bS-PPP, through simulated and real dynamic experiments. Here, ‘E’ indicates the direct utilization of the previous observation corrections from B2b before the signal interruption, ‘P’ represents B2b prediction products, and ‘S’ signifies the incorporation of the SISRE residuals. The results show that EB2b-PPP exhibits significant deviations as early as 10 min into a B2b signal interruption. Both PB2b-PPP and EB2bS-PPP demonstrate comparable performances, with PB2bS-PPP emerging as the most effective method. Notably, in real dynamic experiments, PB2bS-PPP maintains positioning accuracy in the E/N directions like B2b-PPP, even after 40 min of signal interruption, ensuring continuous and stable positioning upon signal restoration. This achievement significantly enhances the capability for high-precision continuous positioning based on B2b signals. Full article
(This article belongs to the Special Issue Advanced Multi-GNSS Positioning and Its Applications in Geoscience)
Show Figures

Figure 1

Figure 1
<p>B2b data collection experimental setup. The <b>left panel</b> is the antenna we used, while the <b>right panel</b> displays the ComNav K823 receiver board.</p>
Full article ">Figure 2
<p>SISRE of BDS-3 for different interruption scenarios (spanning from 10 min to 40 min).</p>
Full article ">Figure 3
<p>SISRE of GPS for different interruption scenarios (spanning from 10 min to 40 min).</p>
Full article ">Figure 4
<p>STD of SISRE for BDS-3 satellites under different interruption durations in non-Interrupted, Ext, and Pred Scenarios (Each label contains three bars representing the three methods: the first column represents B2b, the second column represents the Ext-B2b, and the third column represents the Pred-B2b).</p>
Full article ">Figure 5
<p>STD of SISRE for GPS satellites under different interruption durations in non-Interrupted, Ext, and Pred Scenarios (Each label contains three bars representing the three methods: the first column represents B2b, the second column represents the Ext-B2b, and the third column represents the Pred-B2b).</p>
Full article ">Figure 6
<p>The ZHD0 station used for collecting static experiment data. The <b>left panel</b> displays the antenna used, while the <b>right panel</b> shows the receiver.</p>
Full article ">Figure 7
<p>The effective satellites and PDOP for the ZHD0 on DOY 120 (<b>left</b>) and 121 (<b>right</b>).</p>
Full article ">Figure 8
<p>The positioning errors in the E/N/U directions for the B2b-PPP, EB2b-PPP, PB2b-PPP, EB2bS-PPP, and PB2bS-PPP under signal interruptions of 10 min, 20 min, 30 min, and 40 min on DOY 120 and 121.</p>
Full article ">Figure 9
<p>The maximum positioning errors in the E/N/U directions for the B2b-PPP, EB2b-PPP, PB2b-PPP, EB2bS-PPP, and PB2bS-PPP under signal interruptions of 10 min, 20 min, 30 min, and 40 min on DOY 120 and 121.</p>
Full article ">Figure 10
<p>Dynamic vehicle-mounted positioning experimental setup (<b>left</b>) and trajectory of the measurement campaign (<b>right</b>).</p>
Full article ">Figure 11
<p>The number of effective satellites and PDOP for the kinematic experiment.</p>
Full article ">Figure 12
<p>Positioning errors in the E/N/U directions for various PPP models (B2b-PPP, EB2b-PPP, PB2b-PPP, EB2bS-PPP, and PB2bS-PPP) during a 40 min B2b signal interruption in the kinematic experiment.</p>
Full article ">
22 pages, 24659 KiB  
Article
A Multi-Scale Fusion Deep Learning Approach for Wind Field Retrieval Based on Geostationary Satellite Imagery
by Wei Zhang, Yapeng Wu, Kunkun Fan, Xiaojiang Song, Renbo Pang and Boyu Guoan
Remote Sens. 2025, 17(4), 610; https://doi.org/10.3390/rs17040610 - 11 Feb 2025
Viewed by 324
Abstract
Wind field retrieval, a crucial component of weather forecasting, has been significantly enhanced by recent advances in deep learning. However, existing approaches that are primarily focused on wind speed retrieval are limited by their inability to achieve real-time, full-coverage retrievals at large scales. [...] Read more.
Wind field retrieval, a crucial component of weather forecasting, has been significantly enhanced by recent advances in deep learning. However, existing approaches that are primarily focused on wind speed retrieval are limited by their inability to achieve real-time, full-coverage retrievals at large scales. To address this problem, we propose a novel multi-scale fusion retrieval (MFR) method, leveraging geostationary observation satellites. At the mesoscale, MFR incorporates a cloud-to-wind transformer model, which employs local self-attention mechanisms to extract detailed wind field features. At large scales, MFR incorporates a multi-encoder coordinate U-net model, which incorporates multiple encoders and utilises coordinate information to fuse meso- to large-scale features, enabling accurate and regionally complete wind field retrievals, while reducing the computational resources required. The MFR method was validated using Level 1 data from the Himawari-8 satellite, covering a geographic range of 0–60°N and 100–160°E, at a resolution of 0.25°. Wind field retrieval was accomplished within seconds using a single graphics processing unit. The mean absolute error of wind speed obtained by the MFR was 0.97 m/s, surpassing the accuracy of the CFOSAT and HY-2B Level 2B wind field products. The mean absolute error for wind direction achieved by the MFR was 23.31°, outperforming CFOSAT Level 2B products and aligning closely with HY-2B Level 2B products. The MFR represents a pioneering approach for generating initial fields for large-scale grid forecasting models. Full article
(This article belongs to the Special Issue Image Processing from Aerial and Satellite Imagery)
Show Figures

Figure 1

Figure 1
<p>Vector diagrams for (<b>a</b>) normal conditions, (<b>b</b>) Typhoon Hinnamnor, and (<b>c</b>) Typhoon Nanmadol in the study area. Arrow directions and lengths indicate wind direction and speed, respectively.</p>
Full article ">Figure 2
<p>Multi-scale fusion retrieval architecture. Res, resolution.</p>
Full article ">Figure 3
<p>Sliding window sampling method.</p>
Full article ">Figure 4
<p>Structure of the C2W-Former model.</p>
Full article ">Figure 5
<p>(<b>a</b>) The Swin Transformer block architecture; (<b>b</b>) W-MSA and SW-MSA, which are multi-head self-attention modules with regular and shifted windowing configurations, respectively. Blue and red boxes represent patches and windows, respectively.</p>
Full article ">Figure 6
<p>Discontinuous and blurred boundaries in the preliminary UV.</p>
Full article ">Figure 7
<p>Architecture of the Multi-encoder Coordinate U-net (M-CoordUnet) model. (<b>a</b>) Overall architecture. (<b>b</b>–<b>d</b>) Structural details of the encoder, centre, and decoder blocks, respectively.</p>
Full article ">Figure 8
<p>Analysis of MFR and OSR in comparison to ERA5 for land and sea regions.</p>
Full article ">Figure 9
<p>Scatter plots of each model versus ERA5 at 00:00 on 19 April 2021 (Super Typhoon Surigae). Closer proximity to the x = y line indicates better agreement between the two models. Warmer colours indicate higher frequency.</p>
Full article ">Figure 10
<p>Analysis of ERA5, IFS, and MFR results in comparison to weather station data.</p>
Full article ">Figure 11
<p>UV MAE statistics of MFR wind fields for land (green), sea (orange), and the total study area (blue) across different months and regions in the test data. Solid line indicates average error; shaded area indicates the 95% confidence interval.</p>
Full article ">Figure 12
<p>Comparison of wind field characteristics among different models and data products during Super Typhoon Surigae (19 April 2021, 00:00 UTC).</p>
Full article ">Figure 13
<p>Comparison of wind field characteristics among different models and data products during Super Typhoon Mindulle (28 September 2021, 12:00 UTC).</p>
Full article ">
26 pages, 2493 KiB  
Article
Resource Allocation and Interference Coordination Strategies in Heterogeneous Dual-Layer Satellite Networks
by Jinhong Li, Rong Chai, Tianyi Zhou and Chengchao Liang
Sensors 2025, 25(4), 1005; https://doi.org/10.3390/s25041005 - 8 Feb 2025
Viewed by 263
Abstract
In the face of rapidly evolving communication technologies and increasing user demands, traditional terrestrial networks are challenged by the need for high-quality, high-speed, and reliable communication. This paper explores the integration of heterogeneous satellite networks (HSN) with emerging technologies such as Mobile Edge [...] Read more.
In the face of rapidly evolving communication technologies and increasing user demands, traditional terrestrial networks are challenged by the need for high-quality, high-speed, and reliable communication. This paper explores the integration of heterogeneous satellite networks (HSN) with emerging technologies such as Mobile Edge Computing (MEC), in-network caching, and Software-Defined Networking (SDN) to enhance service efficiency. By leveraging dual-layer satellite networks combining Low Earth Orbit (LEO) and Geostationary Earth Orbit (GEO) satellites, the study addresses resource allocation and interference coordination challenges. This paper proposes a novel resource allocation and interference coordination strategy for dual-layer satellite networks integrating LEO and GEO satellites. We formulate a mathematical optimization problem to optimize resource allocation while minimizing co-channel interference and develop an ADMM-based distributed algorithm for efficient problem-solving. The proposed scheme enhances service efficiency by incorporating MEC, in-network caching, and SDN technologies into the satellite network. Simulation results demonstrate that our proposed algorithm significantly improves network performance by effectively managing resources and reducing interference. Full article
(This article belongs to the Topic Advances in Wireless and Mobile Networking)
Show Figures

Figure 1

Figure 1
<p>Heterogeneous Satellite Network Scene Model.</p>
Full article ">Figure 2
<p>Schematic diagram of antenna off-axis angle and interference distance.</p>
Full article ">Figure 3
<p>Problem Decomposition.</p>
Full article ">Figure 4
<p>Average user vMOS for different computing capacities.</p>
Full article ">Figure 5
<p>Average user vMOS for different numbers of access points.</p>
Full article ">Figure 6
<p>Average MEC server load for different numbers of users.</p>
Full article ">Figure 7
<p>Average MEC server load for different cache capacities.</p>
Full article ">Figure 8
<p>Average MEC server load for different computing capacities.</p>
Full article ">
26 pages, 17033 KiB  
Article
Cost-Effective Satellite Ground Stations in Real-World Development for Space Classrooms
by Pirada Techavijit and Polkit Sukchalerm
Aerospace 2025, 12(2), 105; https://doi.org/10.3390/aerospace12020105 - 30 Jan 2025
Viewed by 731
Abstract
This paper presents the development and outcomes of a cost-effective satellite ground station designed as a learning tool for satellite communication and wireless communication education. The study investigates accessible satellites and the methods for accessing them. The developed ground station has the capability [...] Read more.
This paper presents the development and outcomes of a cost-effective satellite ground station designed as a learning tool for satellite communication and wireless communication education. The study investigates accessible satellites and the methods for accessing them. The developed ground station has the capability to access satellites in the V, U, and L frequency bands, allowing it to receive a variety of satellite data. This includes full-disk meteorological images, high-resolution multispectral images, and scientific data from payloads of satellites in both low Earth orbit (LEO) and geostationary orbit (GEO). The ground station demonstrates capabilities similar to those of large organizations but at a significantly lower cost. This is achieved through a process of identifying educational requirements and optimizing the system for cost-efficiency. This paper presents the design demonstration, actual construction of the ground station, and results. Additionally, it compiles characteristics from real signal reception experiences from various satellites. Full article
(This article belongs to the Section Astronautics & Space Science)
Show Figures

Figure 1

Figure 1
<p>Ground station implement according to this research paper. (<b>a</b>) L-band station, (<b>b</b>) terrestrial and LEO satellite station, (<b>c</b>) ground control system.</p>
Full article ">Figure 2
<p>The ground stations and space classrooms: (<b>a</b>) space classroom, (<b>b</b>) control station in the space classroom, (<b>c</b>) ground stations for learning components.</p>
Full article ">Figure 3
<p>Examples of satellites that the satellite ground station in this research can receive and decode signals from.</p>
Full article ">Figure 4
<p>Circuit of ground station receiver.</p>
Full article ">Figure 5
<p>Conceptual design of satellite ground station.</p>
Full article ">Figure 6
<p>System diagram from antennas to SDR server (numbers in figure refer to BOM in <a href="#aerospace-12-00105-t004" class="html-table">Table 4</a>).</p>
Full article ">Figure 7
<p>L-band antenna system: (<b>a</b>) helical feed 2.5 turns with LHCP and LNA with SAW filter and (<b>b</b>) 1.9 m dish antenna setup.</p>
Full article ">Figure 8
<p>Terrestrial and LEO satellite antenna system. (<b>a</b>) Yagi-Uda antenna system, (<b>b</b>) Grid antenna with LPDA feed.</p>
Full article ">Figure 9
<p>Antenna rotator system: (<b>a</b>) Yaesu G-5500 rotator, (<b>b</b>) the diagram of the rotator controller, (<b>c</b>) the error values from the azimuth axis, and (<b>d</b>) the error values from the elevation axis.</p>
Full article ">Figure 10
<p>SDR server: (<b>a</b>) SDR server of satellite ground station and (<b>b</b>) software diagram on computer.</p>
Full article ">Figure 11
<p>Directional antenna patterns: (<b>a</b>) Yagi–Uda VHF test on 145 MHz, (<b>b</b>) Yagi–Uda UHF antenna tested on 437 MHz, (<b>c</b>) grid antenna of 60 cm × 90 cm tested on 1600 MHz, and (<b>d</b>) dish antenna 1.9 m tested on 1700 MHz.</p>
Full article ">Figure 12
<p>Satellite images from FY-2H after the software process.</p>
Full article ">Figure 13
<p>Raw images downloaded from FY-2H.</p>
Full article ">Figure 14
<p>Image from Electro-L N3 after the software process.</p>
Full article ">Figure 15
<p>Raw images from Electro-L N3.</p>
Full article ">Figure 16
<p>Images from Electro-L N3 after the software process.</p>
Full article ">Figure 17
<p>Image from Electro-L N4 after the software process.</p>
Full article ">Figure 18
<p>Images from GK-2A after the SatDump process.</p>
Full article ">Figure 19
<p>Images form GK-2A: (<b>a</b>) Raw image download, and (<b>b</b>) Meteorological data of surface pressure analysis map.</p>
Full article ">Figure 20
<p>Satellite images from NOAA-19 after the process from WXtoImg (<b>a</b>) Natural color composite image showing cloud cover and land features. (<b>b</b>) False-color composite image highlighting cloud intensity and convection areas.</p>
Full article ">Figure 21
<p>STD−C data from Inmarsat4−F2.</p>
Full article ">Figure 22
<p>Inmarsat Aero from Inmarsat4−F2. The waterfall graph represents signal intensity over time. Brighter colors represent stronger signals, while darker colors indicate weaker or no signals.</p>
Full article ">Figure 23
<p>ADS-B data.</p>
Full article ">
14 pages, 18673 KiB  
Technical Note
A Machine Learning Algorithm to Convert Geostationary Satellite LST to Air Temperature Using In Situ Measurements: A Case Study in Rome and High-Resolution Spatio-Temporal UHI Analysis
by Andrea Cecilia, Giampietro Casasanta, Igor Petenko and Stefania Argentini
Remote Sens. 2025, 17(3), 468; https://doi.org/10.3390/rs17030468 - 29 Jan 2025
Viewed by 391
Abstract
Air temperature (Ta) measurements are crucial for characterizing phenomena like the urban heat island (UHI), which can create critical conditions in cities during summer. This study aims to develop a machine learning-based model, namely gradient boosting, to estimate Ta [...] Read more.
Air temperature (Ta) measurements are crucial for characterizing phenomena like the urban heat island (UHI), which can create critical conditions in cities during summer. This study aims to develop a machine learning-based model, namely gradient boosting, to estimate Ta from geostationary satellite LST data and to apply these estimates to investigate UHI dynamics. Using Rome, Italy, as a case study, the model was trained with Ta data from 15 weather stations, taking multi-temporal LST values (instantaneous and lagged up to 4 h) and additional predictors. The model achieved an overall RMSE of 0.9 °C. The resulting Ta fields, with a 3 km spatial and hourly temporal resolution, enabled a detailed analysis of UHI intensity and dynamics during the summers of 2019–2020, significantly enhancing the spatial and temporal detail compared to previous studies based solely on in situ data. The results also revealed a slightly higher nocturnal UHI intensity than previously reported, attributed to the inclusion of rural areas with near-zero imperviousness, thanks to the complete mapping of Ta across the domain now accessible. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Spatial work domain, highlighting the edges of the LST measurement cells (pixels). In light blue are the cells where in situ meteorological stations are present, indicated by orange triangles. (<b>b</b>) Zoom on the domain to better shows the urbanized areas and the parks, the latter highlighted in green with the name in white. The blue lines indicate the main roads, while the light blue color represents the main watercourses. Finally, the green line indicates the A90 Grande Raccordo Anulare motorway, within which we conventionally define the city of Rome.</p>
Full article ">Figure 2
<p>Some of the dataset predictors, resampled to LST resolution. (<b>a</b>) Altitude a.s.l. (<b>b</b>) Imperviousness. (<b>c</b>) Land Cover. (<b>d</b>) Tree Cover. (<b>e</b>) Grassland. (<b>f</b>) NDVI (July 2019).</p>
Full article ">Figure 3
<p>(<b>a</b>) Diurnal cycle of RMSE; (<b>b</b>) spatial distribution of RMSE, the mean over the entire time period.</p>
Full article ">Figure 4
<p>(<b>a</b>) Importance of predictors in percentage. (<b>b</b>) Correlation matrix between predictors. The abbreviations used are as follows: cell, cell ID in the domain; lon, longitude; lat, latitude; lst_C, synchronous LST; dtm, elevation; imp, imperviousness; lst_DX, LST with a lag of X hours backwards.</p>
Full article ">Figure 5
<p>Spatial distribution of daily extrema of LST and air temperature, including for the latter maps based only on in situ measurements and on model estimates. (<b>a</b>) Minimum <math display="inline"><semantics> <msub> <mi>T</mi> <mi>a</mi> </msub> </semantics></math> (in situ data). (<b>b</b>) Minimum <math display="inline"><semantics> <msub> <mi>T</mi> <mi>a</mi> </msub> </semantics></math> (model output). (<b>c</b>) Minimum LST. (<b>d</b>) Maximum <math display="inline"><semantics> <msub> <mi>T</mi> <mi>a</mi> </msub> </semantics></math> (in situ data). (<b>e</b>) Maximum <math display="inline"><semantics> <msub> <mi>T</mi> <mi>a</mi> </msub> </semantics></math> (model output). (<b>f</b>) Maximum LST.</p>
Full article ">Figure 6
<p>Scatter plots of daily extremes of LST and air temperature (model output data) against IMP (imperviousness), with respective Pearson correlation indices R. (<b>a</b>) Minimum LST. (<b>b</b>) Minimum temperature. (<b>c</b>) Maximum LST. (<b>d</b>) Maximum temperature.</p>
Full article ">Figure 7
<p>(<b>a</b>) Diurnal cycle of UHI intensity estimated using air temperature data; (<b>b</b>) diurnal cycle of SUHI intensity estimated using LST and imperviousness data, along with the diurnal cycle of the correlation coefficient between LST and imperviousness.</p>
Full article ">Figure 8
<p>(<b>a</b>) Diurnal cycle of UHI intensity estimated using air temperature data output from the machine learning model, and SUHI intensity measured using LST data, compared, (<b>b</b>) and the diurnal cycle of the SUHI and UHI gradients.</p>
Full article ">
22 pages, 6054 KiB  
Article
Evaluation and Adjustment of Precipitable Water Vapor Products from FY-4A Using Radiosonde and GNSS Data from China
by Xiangping Chen, Yifei Yang, Wen Liu, Changzeng Tang, Congcong Ling, Liangke Huang, Shaofeng Xie and Lilong Liu
Atmosphere 2025, 16(1), 99; https://doi.org/10.3390/atmos16010099 - 17 Jan 2025
Viewed by 458
Abstract
The geostationary meteorological satellite Fengyun-4A (FY-4A) has rapidly advanced, generating abundant high spatiotemporal resolution atmospheric precipitable water vapor (PWV) products. However, remote sensing satellites are vulnerable to weather conditions, and these latest operational PWV products still require systematic validation. This study presents a [...] Read more.
The geostationary meteorological satellite Fengyun-4A (FY-4A) has rapidly advanced, generating abundant high spatiotemporal resolution atmospheric precipitable water vapor (PWV) products. However, remote sensing satellites are vulnerable to weather conditions, and these latest operational PWV products still require systematic validation. This study presents a comprehensive evaluation of FY-4A PWV products by separately using PWV data retrieved from radiosondes (RS) and the Global Navigation Satellite System (GNSS) from 2019 to 2022 in China and the surrounding regions. The overall results indicate a significant consistency between FY-4A PWV and RS PWV as well as GNSS PWV, with mean biases of 7.21 mm and −8.85 mm, and root mean square errors (RMSEs) of 7.03 mm and 3.76 mm, respectively. In terms of spatial variability, the significant differences in mean bias and RMSE were 6.50 mm and 2.60 mm between FY-4A PWV and RS PWV in the northern and southern subregions, respectively, and 5.36 mm and 1.73 mm between FY-4A PWV and GNSS PWV in the northwestern and southern subregions, respectively. The RMSE of FY-4A PWV generally increases with decreasing latitude, and the bias is predominantly negative, indicating an underestimation of water vapor. Regarding temporal differences, both the monthly and daily biases and RMSEs of FY-4A PWV are significantly higher in summer than in winter, with daily precision metrics in summer displaying pronounced peaks and irregular fluctuations. The classic seasonal, regional adjustment model effectively reduced FY-4A PWV deviations across all regions, especially in the NWC subregion with low water vapor distribution. In summary, the accuracy metrics of FY-4A PWV show distinct spatiotemporal variations compared to RS PWV and GNSS PWV, and these variations should be considered to fully realize the potential of multi-source water vapor applications. Full article
(This article belongs to the Special Issue GNSS Remote Sensing in Atmosphere and Environment (2nd Edition))
Show Figures

Figure 1

Figure 1
<p>Distribution of RS sites and GNSS sites from 2019–2022 in the research area.</p>
Full article ">Figure 2
<p>Observation mode of the AGRI on FY-4A satellite. The vertical axis represents UTC time in hours, while the horizontal axis represents the minutes within each hour.</p>
Full article ">Figure 3
<p>Fitting plots between RS PWV and FY-4A PWV from 2019 to 2022 for different regions, with correlation, annual mean bias, and RMSE values.</p>
Full article ">Figure 4
<p>Site distribution maps of the mean bias and mean RMSE between FY-4A PWV and RS PWV from 2019 to 2022.</p>
Full article ">Figure 5
<p>Histograms of annual mean bias and RMSE between FY-4A PWV and RS PWV in different regions.</p>
Full article ">Figure 6
<p>Seasonal average distribution of FY-4A PWV and GNSS PWV for 2022.</p>
Full article ">Figure 7
<p>Fitting plots between FY-4A PWV and GNSS PWV from 2019 to 2022 for different regions, with correlation, annual mean bias, and RMSE values.</p>
Full article ">Figure 8
<p>Site distribution maps of the mean bias and mean RMSE between FY-4A PWV and GNSS PWV from 2019 to 2022.</p>
Full article ">Figure 9
<p>Bar charts of monthly mean bias and RMSE for four seasons between FY-4A PWV and GNSS PWV from 2019–2022.</p>
Full article ">Figure 10
<p>Box plots of monthly mean bias and RMSE between FY-4A PWV and GNSS PWV from 2019–2022 in different regions. Q1 and Q3 of the box are the first and third quartiles, respectively. The distance between Q1 and Q3 reflects the degree of fluctuation of the data; Q2 is the median value, which reflects the average level of the data; Q4 is the outlier.</p>
Full article ">Figure 11
<p>Time series of daily mean bias and RMSE between FY-4A PWV and GNSS PWV in different regions from 2019 to 2022.</p>
Full article ">Figure 12
<p>Bar charts of the mean MAE and RMSE between FY-4A PWV and GNSS PWV before and after adjustment in different regions and seasons for 2022. The length of the arrows represents the degree of improvement in mean MAE and RMSE.</p>
Full article ">Figure 13
<p>Site-level distribution of seasonal average improvements in MAE and RMSE between corrected and uncorrected FY-4A PWV and GNSS PWV for 2022. IMAE and IRMSE represent the improved MAE and RMSE values, respectively.</p>
Full article ">
20 pages, 17962 KiB  
Article
Conversion of 10 min Rain Rate Time Series into 1 min Time Series: Theory, Experimental Results, and Application in Satellite Communications
by Emilio Matricciani and Carlo Riva
Appl. Sci. 2025, 15(2), 743; https://doi.org/10.3390/app15020743 - 13 Jan 2025
Viewed by 811
Abstract
We propose a semi-empirical method—based on a filtered Markov process—to convert 10 min rain rate time series into 1 min time series, i.e., quasi-instantaneous rainfall—the latter to be used as input to the synthetic storm technique, which is a very reliable tool for [...] Read more.
We propose a semi-empirical method—based on a filtered Markov process—to convert 10 min rain rate time series into 1 min time series, i.e., quasi-instantaneous rainfall—the latter to be used as input to the synthetic storm technique, which is a very reliable tool for calculating rain attenuation time series in satellite communication systems or for estimating runoff, erosion, pollutant transport, and other applications in hydrology. To develop the method, we used a very large data bank of 1 min rain rate time series collected in several sites with different climatic conditions. The experimental and simulated 1 min rain rate time series agree very well. Afterward, we used them to simulate rain attenuation time series at 20.7 GHz, in 35.5° slant paths to geostationary satellites. The two simulated annual rain attenuation probability distributions show very small differences. We conclude that the rain rate conversion method is very reliable. Full article
(This article belongs to the Special Issue Advanced Technologies in Optical and Microwave Transmission)
Show Figures

Figure 1

Figure 1
<p><math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> (cyan) and corresponding <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>10</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> (magenta). Both rain rates are expressed in mm/h. Spino d’Adda, 20 October 2000; the event starts at 10:32.</p>
Full article ">Figure 2
<p>Histograms of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> in the ranges of <math display="inline"><semantics> <mrow> <mn>0</mn> <mo>–</mo> <mn>2</mn> </mrow> </semantics></math> mm/h of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>10</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> (<b>left panel</b>) and <math display="inline"><semantics> <mrow> <mn>2</mn> <mo>–</mo> <mn>4</mn> </mrow> </semantics></math> mm/h (<b>right panel</b>) of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>10</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Histograms of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> in the ranges of <math display="inline"><semantics> <mrow> <mn>4</mn> <mo>–</mo> <mn>6</mn> </mrow> </semantics></math> mm/h of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>10</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> (<b>left panel</b>) and <math display="inline"><semantics> <mrow> <mn>6</mn> <mo>–</mo> <mn>8</mn> </mrow> </semantics></math> mm/h (<b>right panel</b>) of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>10</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Histograms of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> in the ranges of <math display="inline"><semantics> <mrow> <mn>8</mn> <mo>–</mo> <mn>10</mn> </mrow> </semantics></math> mm/h of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>10</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> (<b>left panel</b>) and <math display="inline"><semantics> <mrow> <mn>10</mn> <mo>–</mo> <mn>15</mn> </mrow> </semantics></math> mm/h (<b>right panel</b>) of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>10</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>Histograms of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> in the ranges of <math display="inline"><semantics> <mrow> <mn>15</mn> <mo>–</mo> <mn>20</mn> </mrow> </semantics></math> mm/h of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>10</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> (<b>left panel</b>) and <math display="inline"><semantics> <mrow> <mn>20</mn> <mo>–</mo> <mn>30</mn> </mrow> </semantics></math> mm/h (<b>right panel</b>) of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>10</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>Histograms of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> in the ranges of <math display="inline"><semantics> <mrow> <mn>30</mn> <mo>–</mo> <mn>40</mn> </mrow> </semantics></math> mm/h of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>10</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> (<b>left panel</b>) and <math display="inline"><semantics> <mrow> <mo>&gt;</mo> <mn>40</mn> </mrow> </semantics></math> mm/h (<b>right panel</b>) of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>10</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 7
<p>Experimental <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> (mm/h) (blue, original) and simulated <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> <mo>,</mo> <mi>s</mi> <mi>i</mi> <mi>m</mi> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> (mm/h) time series (black, simul). Left: low-intensity rain rate event. Right panel: high-intensity rain rate event. The 10 min quantity of water is conserved.</p>
Full article ">Figure 8
<p>Mean value (<b>left panel</b>, mm/h) and standard deviation (<b>right panel</b>, mm/h) of <math display="inline"><semantics> <mrow> <mi>ε</mi> </mrow> </semantics></math> versus <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>f</mi> </mrow> <mrow> <mi>c</mi> </mrow> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 9
<p>Example of 1 min rain rate time series, measured (blue line, original) and simulated (red line, gener), after filtering and water conservation. (<b>Left panel</b>): a low rain rate event. (<b>Right panel</b>): a high-intensity rain rate event (see also <a href="#applsci-15-00743-f007" class="html-fig">Figure 7</a>).</p>
Full article ">Figure 10
<p>Probability distribution (PD) that the 1 min rain rate in abscissa is exceeded in the experimental data <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </semantics></math>; blue line (original), and in the simulated 1 min data, <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math>, black line (simul).</p>
Full article ">Figure 11
<p>Scatterplots of mean values (<b>left panel</b>), standard deviations (<b>central panel</b>), and correlation coefficients (<b>right panel</b>) between the values of the sites in <a href="#applsci-15-00743-t001" class="html-table">Table 1</a> and Spino d’Adda. Gera Lario: green; Fucino: blue; Madrid: cyan; Prague: yellow; Tampa: red; White Sands: magenta; Vancouver: black.</p>
Full article ">Figure 12
<p><b>Gera Lario.</b> Probability distribution that the 1 min rain rate in abscissa is exceeded in the experimental data, <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </semantics></math>, blue line, and in the simulated 1 min data, <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math>, black line. <b>Left panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using local values of the conditional PDFs. <b>Right panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using Spino d’Adda conditional PDFs (<a href="#applsci-15-00743-t002" class="html-table">Table 2</a>).</p>
Full article ">Figure 13
<p><b>Fucino.</b> Probability distribution that the 1 min rain rate in abscissa is exceeded in the experimental data, <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </semantics></math>, blue line, and in the simulated 1 min data, <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math>, black line. <b>Left panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using local values of the conditional PDFs. <b>Right panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using Spino d’Adda conditional PDFs (<a href="#applsci-15-00743-t002" class="html-table">Table 2</a>).</p>
Full article ">Figure 14
<p><b>Madrid.</b> Probability distribution that the 1 min rain rate in abscissa is exceeded in the experimental data, <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </semantics></math>, blue line, and in the simulated 1 min data, <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math>, black line. <b>Left panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using local values of the conditional PDFs. <b>Right panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using Spino d’Adda conditional PDFs (<a href="#applsci-15-00743-t002" class="html-table">Table 2</a>).</p>
Full article ">Figure 15
<p><b>Prague.</b> Probability distribution that the 1 min rain rate in abscissa is exceeded in the experimental data, <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </semantics></math>, blue line, and in the simulated 1 min data, <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math>, black line. <b>Left panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using local values of the conditional PDFs. <b>Right panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using Spino d’Adda conditional PDFs (<a href="#applsci-15-00743-t002" class="html-table">Table 2</a>).</p>
Full article ">Figure 16
<p><b>Tampa.</b> Probability distribution that the 1 min rain rate in abscissa is exceeded in the experimental data, <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </semantics></math>, blue line, and in the simulated 1 min data, <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math>, black line. <b>Left panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using local values of the conditional PDFs. <b>Right panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using Spino d’Adda conditional PDFs (<a href="#applsci-15-00743-t002" class="html-table">Table 2</a>).</p>
Full article ">Figure 17
<p><b>White Sands.</b> Probability distribution that the 1 min rain rate in abscissa is exceeded in the experimental data, <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </semantics></math>, blue line, and in the simulated 1 min data, <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math>, black line. <b>Left panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using local values of the conditional PDFs. <b>Right panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using Spino d’Adda conditional PDFs (<a href="#applsci-15-00743-t002" class="html-table">Table 2</a>).</p>
Full article ">Figure 18
<p><b>Vancouver.</b> Probability distribution that the 1 min rain rate in abscissa is exceeded in the experimental data, <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </semantics></math>, blue line, and in the simulated 1 min data, <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math>, black line. <b>Left panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using local values of the conditional PDFs. <b>Right panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>R</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using Spino d’Adda conditional PDFs (<a href="#applsci-15-00743-t002" class="html-table">Table 2</a>).</p>
Full article ">Figure 19
<p>Average annual probability distribution <math display="inline"><semantics> <mrow> <mi>P</mi> <mfenced separators="|"> <mrow> <mi>A</mi> </mrow> </mfenced> </mrow> </semantics></math>—namely, the fraction of time of a year—that the rain attenuation <math display="inline"><semantics> <mrow> <mi>A</mi> </mrow> </semantics></math> (dB) in abscissa is exceeded, estimated with the SST. Cyan line: experimental <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>; magenta line: simulated <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 20
<p><b>Gera Lario.</b> Average annual probability distribution <math display="inline"><semantics> <mrow> <mi>P</mi> <mfenced separators="|"> <mrow> <mi>A</mi> </mrow> </mfenced> </mrow> </semantics></math> that the rain attenuation <math display="inline"><semantics> <mrow> <mi>A</mi> </mrow> </semantics></math> (dB) in abscissa is exceeded, estimated with the SST; blue line: experimental <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>, black line: simulated <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>. <b>Left panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>A</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using local values of the conditional rain rate PDFs. <b>Right panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>A</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using Spino d’Adda conditional PDFs (<a href="#applsci-15-00743-t002" class="html-table">Table 2</a>).</p>
Full article ">Figure 21
<p><b>Fucino.</b> Average annual probability distribution <math display="inline"><semantics> <mrow> <mi>P</mi> <mfenced separators="|"> <mrow> <mi>A</mi> </mrow> </mfenced> </mrow> </semantics></math> that the rain attenuation <math display="inline"><semantics> <mrow> <mi>A</mi> </mrow> </semantics></math> (dB) in abscissa is exceeded, estimated with the SST; blue line: experimental <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>, black line: simulated <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>. <b>Left panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>A</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using local values of the conditional rain rate PDFs. <b>Right panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>A</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using Spino d’Adda conditional PDFs (<a href="#applsci-15-00743-t002" class="html-table">Table 2</a>).</p>
Full article ">Figure 22
<p><b>Madrid.</b> Average annual probability distribution <math display="inline"><semantics> <mrow> <mi>P</mi> <mfenced separators="|"> <mrow> <mi>A</mi> </mrow> </mfenced> </mrow> </semantics></math> that the rain attenuation <math display="inline"><semantics> <mrow> <mi>A</mi> </mrow> </semantics></math> (dB) in abscissa is exceeded, estimated with the SST; blue line: experimental <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>, black line: simulated <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>. <b>Left panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>A</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using local values of the conditional rain rate PDFs. <b>Right panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>A</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using Spino d’Adda conditional PDFs (<a href="#applsci-15-00743-t002" class="html-table">Table 2</a>).</p>
Full article ">Figure 23
<p><b>Prague.</b> Average annual probability distribution <math display="inline"><semantics> <mrow> <mi>P</mi> <mfenced separators="|"> <mrow> <mi>A</mi> </mrow> </mfenced> </mrow> </semantics></math> that the rain attenuation <math display="inline"><semantics> <mrow> <mi>A</mi> </mrow> </semantics></math> (dB) in abscissa is exceeded, estimated with the SST; blue line: experimental <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>, black line: simulated <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>. <b>Left panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>A</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using local values of the conditional rain rate PDFs. <b>Right panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>A</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using Spino d’Adda conditional PDFs (<a href="#applsci-15-00743-t002" class="html-table">Table 2</a>).</p>
Full article ">Figure 24
<p><b>Tampa.</b> Average annual probability distribution <math display="inline"><semantics> <mrow> <mi>P</mi> <mfenced separators="|"> <mrow> <mi>A</mi> </mrow> </mfenced> </mrow> </semantics></math> that the rain attenuation <math display="inline"><semantics> <mrow> <mi>A</mi> </mrow> </semantics></math> (dB) in abscissa is exceeded, estimated with the SST; blue line: experimental <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>, black line: simulated <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>. <b>Left panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>A</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using local values of the conditional rain rate PDFs. <b>Right panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>A</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using Spino d’Adda conditional PDFs (<a href="#applsci-15-00743-t002" class="html-table">Table 2</a>).</p>
Full article ">Figure 25
<p><b>White Sands.</b> Average annual probability distribution <math display="inline"><semantics> <mrow> <mi>P</mi> <mfenced separators="|"> <mrow> <mi>A</mi> </mrow> </mfenced> </mrow> </semantics></math> that the rain attenuation <math display="inline"><semantics> <mrow> <mi>A</mi> </mrow> </semantics></math> (dB) in abscissa is exceeded, estimated with the SST; blue line: experimental <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>, black line: simulated <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>. <b>Left panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>A</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using local values of the conditional rain rate PDFs. <b>Right panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>A</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using Spino d’Adda conditional PDFs (<a href="#applsci-15-00743-t002" class="html-table">Table 2</a>).</p>
Full article ">Figure 26
<p><b>Vancouver.</b> Average annual probability distribution <math display="inline"><semantics> <mrow> <mi>P</mi> <mfenced separators="|"> <mrow> <mi>A</mi> </mrow> </mfenced> </mrow> </semantics></math> that the rain attenuation <math display="inline"><semantics> <mrow> <mi>A</mi> </mrow> </semantics></math> (dB) in abscissa is exceeded, estimated with the SST; blue line: experimental <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>, black line: simulated <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>. <b>Left panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>A</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using local values of the conditional rain rate PDFs. <b>Right panel</b>: <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>A</mi> <mo>)</mo> </mrow> </semantics></math> is obtained by using Spino d’Adda conditional PDFs (<a href="#applsci-15-00743-t002" class="html-table">Table 2</a>).</p>
Full article ">
21 pages, 13194 KiB  
Article
A Multi-Layer Perceptron Approach to Downscaling Geostationary Land Surface Temperature in Urban Areas
by Alexandra Hurduc, Sofia L. Ermida and Carlos C. DaCamara
Remote Sens. 2025, 17(1), 45; https://doi.org/10.3390/rs17010045 - 27 Dec 2024
Viewed by 499
Abstract
Remote sensing of land surface temperature (LST) is a fundamental variable in analyzing temperature variability in urban areas. Geostationary sensors provide sufficient observations throughout the day for a diurnal analysis of temperature, however, lack the spatial resolution needed for highly heterogeneous areas such [...] Read more.
Remote sensing of land surface temperature (LST) is a fundamental variable in analyzing temperature variability in urban areas. Geostationary sensors provide sufficient observations throughout the day for a diurnal analysis of temperature, however, lack the spatial resolution needed for highly heterogeneous areas such as cities. Polar orbiting sensors have the advantage of a higher spatial resolution, enabling a better characterization of the surface while only providing one to two observations per day. This work aims at using a multi-layer perceptron-based method to downscale geostationary-derived LST based on a polar-orbit-derived one. The model is trained on a pixel-by-pixel basis, which reduces the complexity of the model while requiring fewer auxiliary data to characterize the surface conditions. Results show that the model is able to successfully downscale LST for the city of Madrid, from approximately 4.5 km to 750 m. Performance metrics between training and validation datasets show no overfitting. The model was applied to a different time period and compared to data derived from three additional sensors, which were not used in any stage of the training process, yielding a R2 of 0.99, root mean square errors between 1.45 and 1.58 and mean absolute errors ranging from 1.07 to 1.15. The downscaled LST is shown to improve the representation of both the temporal variability and spatial heterogeneity of temperature, when compared to geostationary- and polar-orbit-derived LST individually. The resulting downscaled data take advantage of the high observation frequency of geostationary data, combined with the spatial resolution of polar orbiting sensors and may be of added value for the study of diurnal and seasonal patterns of LST in urban environments. Full article
(This article belongs to the Special Issue Advances in Thermal Infrared Remote Sensing II)
Show Figures

Figure 1

Figure 1
<p>Urban/rural distinction for the city of Madrid. Boxes represent the sub-regions used to analyze the diurnal cycle of the downscaled LST. Red—urban, blue—rural, yellow—mixed urban and rural.</p>
Full article ">Figure 2
<p>The structure of the MLP model used to downscale LST. The first layer corresponds to the input layer with three neurons, the second one is the hidden layer with five neurons (N<sub>1…5</sub>) and the final layer represents the output layer, with one neuron.</p>
Full article ">Figure 3
<p>Diagram of the workflow employed for the downscaling of geostationary LST to a regular 750 m.</p>
Full article ">Figure 4
<p>(<b>a</b>) Loss associated with each model/pixel and (<b>b</b>) number of iterations until convergence.</p>
Full article ">Figure 5
<p>Estimated LST versus target LST: (<b>a</b>) training dataset, (<b>b</b>) test dataset.</p>
Full article ">Figure 6
<p>Estimated LST versus observed LST by sensor.</p>
Full article ">Figure 7
<p>Distribution of LST difference between estimated and observed LST by sensor.</p>
Full article ">Figure 8
<p>Maps of bias (<b>a</b>–<b>d</b>) and RMSE (<b>e</b>–<b>h</b>) for each sensor.</p>
Full article ">Figure 9
<p>LST difference per class of estimated LST. The number of pixels corresponding to each boxplot is shown above each one.</p>
Full article ">Figure 10
<p>Difference between coarse and downscaled LST by land cover fraction.</p>
Full article ">Figure 11
<p>Seasonal diurnal cycle of SEVIRI observations, polar orbiting sensor observations and estimated LST for the three sub-regions identified In <a href="#remotesensing-17-00045-f001" class="html-fig">Figure 1</a>. (<b>a</b>–<b>d</b>) Rural, (<b>e</b>–<b>h</b>) mixed rural and urban, (<b>i</b>–<b>l</b>) urban.</p>
Full article ">Figure 12
<p>LST during 23 September 2023, during four times of day: (<b>a</b>–<b>d</b>) 0200UTC; (<b>e</b>–<b>h</b>) 1030UTC; (<b>i</b>–<b>l</b>) 1300UTC; (<b>m</b>–<b>p</b>) 2145UTC.</p>
Full article ">
34 pages, 10549 KiB  
Review
Multi-Sensor Precipitation Estimation from Space: Data Sources, Methods and Validation
by Ruifang Guo, Xingwang Fan, Han Zhou and Yuanbo Liu
Remote Sens. 2024, 16(24), 4753; https://doi.org/10.3390/rs16244753 - 20 Dec 2024
Viewed by 774
Abstract
Satellite remote sensing complements rain gauges and ground radars as the primary sources of precipitation data. While significant advancements have been made in spaceborne precipitation estimation since the 1960s, the emergence of multi-sensor precipitation estimation (MPE) in the early 1990s revolutionized global precipitation [...] Read more.
Satellite remote sensing complements rain gauges and ground radars as the primary sources of precipitation data. While significant advancements have been made in spaceborne precipitation estimation since the 1960s, the emergence of multi-sensor precipitation estimation (MPE) in the early 1990s revolutionized global precipitation data generation by integrating infrared and microwave observations. Among others, Global Precipitation Measurement (GPM) plays a crucial role in providing invaluable data sources for MPE by utilizing passive microwave sensors and geostationary infrared sensors. MPE represents the current state-of-the-art approach for generating high-quality, high-resolution global satellite precipitation products (SPPs), employing various methods such as cloud motion analysis, probability matching, adjustment ratios, regression techniques, neural networks, and weighted averaging. International collaborations, such as the International Precipitation Working Group and the Precipitation Virtual Constellation, have significantly contributed to enhancing our understanding of the uncertainties associated with MPEs and their corresponding SPPs. It has been observed that SPPs exhibit higher reliability over tropical oceans compared to mid- and high-latitudes, particularly during cold seasons or in regions with complex terrains. To further advance MPE research, future efforts should focus on improving accuracy for extremely low- and high-precipitation events, solid precipitation measurements, as well as orographic precipitation estimation. Full article
(This article belongs to the Special Issue Synergetic Remote Sensing of Clouds and Precipitation II)
Show Figures

Figure 1

Figure 1
<p>A brief history of precipitation-observing techniques, experiments, and products.</p>
Full article ">Figure 2
<p>GPM constellation. The left figure was obtained from <a href="https://gpm.nasa.gov/image-gallery/gpm" target="_blank">https://gpm.nasa.gov/image-gallery/gpm</a> (accessed on 1 December 2024).</p>
Full article ">Figure 3
<p>Summary of major global satellite precipitation products currently available.</p>
Full article ">Figure 4
<p>Number of SPP validation studies published over the last two decades (covered in Web of Science Core Collection). We used the keywords “validation” or “evaluation” or “assessment” for the topic and “IMERG”, “PERSIANN”, “CMORPH”, “GSMaP”, “CMAP and Merged Analysis of Precipitation”, “GPCP” or “TMPA or 3B42” for the abstract, focusing on the period between 2020 and 2024, the period between 2015 and 2019, the period between 2010 and 2014, and the period between 2000 and 2009.</p>
Full article ">Figure 5
<p>Schematic diagram showing the SPE validation process.</p>
Full article ">
21 pages, 13566 KiB  
Article
Assimilation of Fengyun-4A Atmospheric Motion Vectors and Its Impact on China Meteorological Administration—Beijing System Forecasts
by Yanhui Xie, Shuting Zhang, Xin Sun, Min Chen, Jiancheng Shi, Yu Xia and Ruixia Liu
Remote Sens. 2024, 16(23), 4561; https://doi.org/10.3390/rs16234561 - 5 Dec 2024
Viewed by 578
Abstract
The ever-increasing capacity of numerical weather prediction (NWP) models requires accurate flow information at higher spatial and temporal resolutions. The atmospheric motion vectors (AMVs) extracted from the Advanced Geostationary Radiation Imager (AGRI) mounted on the Fengyun-4A (FY-4A) satellite can provide information about atmospheric [...] Read more.
The ever-increasing capacity of numerical weather prediction (NWP) models requires accurate flow information at higher spatial and temporal resolutions. The atmospheric motion vectors (AMVs) extracted from the Advanced Geostationary Radiation Imager (AGRI) mounted on the Fengyun-4A (FY-4A) satellite can provide information about atmospheric flow fields on small scales. This study focused on the assimilation of FY-4A AMVs and its impact on forecasts in the regional NWP system of the China Meteorological Administration—Beijing (CAM-BJ). The statistical characterization of FY-4A AMVs was firstly analyzed, and an optimal observation error in each vertical level was obtained. Three groups of retrospective runs over a one-month period were conducted, and the impact of assimilating the AMVs with different strategies on the forecasts of the CMA-BJ system were compared and evaluated. The results suggested that the optimal observation errors reduced the standard deviation of the background departures for U and V wind, leading to an improvement in the standard deviation in the corresponding analysis departures of about 8.3% for U wind and 7.3% for V wind. Assimilating FY-4A AMV data with a quality indicator (QI) above 80 and the optimal observation errors reduced the error of upper wind forecast in the CMA-BJ system. A benefit was also obtained in the error of surface wind forecast after 6 h of the forecasts, although it was not significant. For rainfall forecast with different thresholds, the score skills increased slightly after 6 h of the forecasts. There was an overall improvement for the overprediction of 24 h accumulated precipitation forecast including the AMVs, even when conventional observations were relatively rich. The application of FY-4A AMVs with a QI > 80 and adjustment to observation errors has a positive impact on the upper wind forecast in the CMA-BJ system, improving the score skill of rainfall forecasting. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The coverage of the two domains in the CMA-BJ system and the distribution of observation data used for assimilation.</p>
Full article ">Figure 2
<p>FY-4A AMV spatial patterns at 00 UTC 1 July 2021.</p>
Full article ">Figure 3
<p>Observation numbers and error characterization of AMVs with different QI values in vertical levels against GFS reanalysis data over a period of one month from 1 to 31 July 2021.</p>
Full article ">Figure 4
<p>The bias and RMSE of all FY-4A AMVs and the AMVs with a QI &gt; 80 against the reanalysis data from the NCEP over a period of one month from 1 to 31 July 2021.</p>
Full article ">Figure 5
<p>The observation errors and the corresponding sample data numbers by vertical level for the AMVs with a QI &gt; 80 over one month from 1 to 31 July 2021.</p>
Full article ">Figure 6
<p>First guess (blue) and analysis (red) of U (the first line) and V (the second line) wind versus their corresponding values of the AMV data derived from the FY-4A satellite with (<b>a</b>) the default and (<b>b</b>) the optimized observation errors.</p>
Full article ">Figure 7
<p>Statistics of the background departures for U and V wind in amv_qi80 and amv_qi80uperr.</p>
Full article ">Figure 8
<p>Time series of the observation numbers and the departures of first guess and analysis for U and V wind components versus the corresponding values of the AMVs over one month from 1 to 31 July 2021.</p>
Full article ">Figure 9
<p>Mean biases of wind forecast in CTRL, amv_qi80, and amv_qi80uperr averaged over the 9 km domain in the CMA-BJ system at 12 h and 24 h forecasts. The first line is for U (<b>a</b>) and V (<b>b</b>) winds for 12 h forecasts, and the second line is for U (<b>c</b>) and V (<b>d</b>) for 24 h.</p>
Full article ">Figure 10
<p>RMSEs of forecast wind in the CTRL, amv_qi80, and amv_qi80uperr averaged over the 9 km domain in the CMA-BJ system for 12 h and 24 h forecasts. The first line is for U (<b>a</b>) and V (<b>b</b>) winds for 12 h forecasts, and the second line is for U (<b>c</b>) and V (<b>d</b>) for 24 h.</p>
Full article ">Figure 11
<p>(<b>a</b>) Mean biases and (<b>b</b>) RMSEs over forecast time for 10 m wind forecasts from CTRL, amv_qi80 and amv_qi80uperr against observations averaged over 9 km domain in CMA-BJ system.</p>
Full article ">Figure 12
<p><span class="html-italic">TS</span> scores of 6 h accumulated rainfall forecast from three retrospective runs of CTRL, amv_qi80, and amv_qi80uperr over 9 km domain of CMA-BJ system.</p>
Full article ">Figure 13
<p>Bias scores for 6 h accumulated rainfall forecast from the three retrospective runs of CTRL, amv_qi80, and amv_qi80uperr over the 9 km domain of the CMA-BJ system.</p>
Full article ">Figure 14
<p>Performance diagram of 24 h accumulated rainfall forecast for the three retrospective runs over the 62 forecasts starting at 00 UTC and 12 UTC every day.</p>
Full article ">Figure 15
<p>The comprehensive scorecard for humidity, temperature, and wind forecasts from amv_qi80uperr against the CTRL over 62 forecasts starting at 00 UTC and 12 UTC every day.</p>
Full article ">
24 pages, 13737 KiB  
Article
Generating a 30 m Hourly Land Surface Temperatures Based on Spatial Fusion Model and Machine Learning Algorithm
by Qin Su, Yuan Yao, Cheng Chen and Bo Chen
Sensors 2024, 24(23), 7424; https://doi.org/10.3390/s24237424 - 21 Nov 2024
Viewed by 944
Abstract
Land surface temperature (LST) is a critical parameter for understanding climate change and maintaining hydrological balance across local and global scales. However, existing satellite LST products face trade-offs between spatial and temporal resolutions, making it challenging to provide all-weather LST with high spatiotemporal [...] Read more.
Land surface temperature (LST) is a critical parameter for understanding climate change and maintaining hydrological balance across local and global scales. However, existing satellite LST products face trade-offs between spatial and temporal resolutions, making it challenging to provide all-weather LST with high spatiotemporal resolution. In this study, focusing on Chengdu city, a framework combining a spatiotemporal fusion model and machine learning algorithm was proposed and applied to retrieve hourly high spatial resolution LST data from Chinese geostationary weather satellite data and multi-scale polar-orbiting satellite observations. The predicted 30 m hourly LST values were evaluated against in situ LST measurements and Sentinel-3 SLSTR data on 11 August 2019 and 21 April 2022, respectively. The results demonstrate that validation based on the in situ LST, the root mean squared error (RMSE) of the predicted LST using the proposed framework are around 0.89 °C to 1.23 °C. The predicted LST is highly consistent with the Sentinel-3 SLSTR data, and the RMSE varies from 0.95 °C to 1.25 °C. In addition, the proposed framework was applied to Xi’an City, and the final validation results indicate that the method is accurate to within about 1.33 °C. The generated 30 m hourly LST can provide important data with fine spatial resolution for urban thermal environment monitoring. Full article
(This article belongs to the Section Environmental Sensing)
Show Figures

Figure 1

Figure 1
<p>Location and land-cover maps of the study area.</p>
Full article ">Figure 2
<p>Sample points from the Google Earth image used to validate the accuracy of land cover classification.</p>
Full article ">Figure 3
<p>Flowchart of the proposed framework with FY-4A, MOD11A1, and downscaled LST data.</p>
Full article ">Figure 4
<p>Flowchart of LST downscaling procedure using the machine learning methods.</p>
Full article ">Figure 5
<p>Comparison between the observed LST and predicted LST. (<b>a</b>) observed FY-4A LST, 11:00 local time, 11 August 2019. (<b>b</b>) observed MOD11A1, 11:00 local time, 11 August 2019. (<b>c</b>) observed FY-4A LST, 14:00 local time, 11 August 2019. (<b>d</b>) predicted LST, 14:00 local time, 11 August 2019. (<b>e</b>) observed MYD11A1, 14:00 local time, 11 August 2019. (<b>f</b>) observed FY-4A LST, 11:00 local time, 21 April 2022. (<b>g</b>) observed MOD11A1, 11:00 local time, 21 April 2022. (<b>h</b>) observed FY-4A LST, 14:00 local time, 21 April 2022. (<b>i</b>) predicted LST, 14:00 local time, 21 April 2022. (<b>j</b>) observed MYD11A1, 14:00 local time, 21 April 2022.</p>
Full article ">Figure 6
<p>Scatter plot comparison between predicted LST by CFSDAF and observed MYD11A1 LST for (<b>a</b>) 14:00 on 11 August 2019 and (<b>b</b>) 14:00 on 21 April 2022.</p>
Full article ">Figure 7
<p>The generated 1 km hourly LST dataset: (<b>a</b>) 11 August 2019, and (<b>b</b>) 21 April 2022.</p>
Full article ">Figure 8
<p>Comparison between observed 100 m Landsat 8 LST and downscaled LST with 100 m spatial resolution using machine learning algorithms.</p>
Full article ">Figure 9
<p>The generated 30 m hourly LST dataset: (<b>a</b>) 11 August 2019, and (<b>b</b>) 21 April 2022.</p>
Full article ">Figure 10
<p>Scatter plots of the relationship between the predicted LST results and observed in situ LST on 11 August 2019 and 21 April 2022: (<b>a</b>) 00:00; (<b>b</b>) 01:00; (<b>c</b>) 02:00; (<b>d</b>) 03:00; (<b>e</b>) 04:00; (<b>f</b>) 05:00; (<b>g</b>) 06:00; (<b>h</b>) 07:00; (<b>i</b>) 08:00; (<b>j</b>) 09:00; (<b>k</b>) 10:00; (<b>l</b>) 11:00; (<b>m</b>) 12:00; (<b>n</b>) 13:00; (<b>o</b>) 14:00; (<b>p</b>) 15:00; (<b>q</b>) 16:00; (<b>r</b>) 17:00; (<b>s</b>) 18:00; (<b>t</b>) 19:00; (<b>u</b>) 20:00; (<b>v</b>) 21:00; (<b>w</b>) 22:00; (<b>x</b>) 23:00.</p>
Full article ">Figure 10 Cont.
<p>Scatter plots of the relationship between the predicted LST results and observed in situ LST on 11 August 2019 and 21 April 2022: (<b>a</b>) 00:00; (<b>b</b>) 01:00; (<b>c</b>) 02:00; (<b>d</b>) 03:00; (<b>e</b>) 04:00; (<b>f</b>) 05:00; (<b>g</b>) 06:00; (<b>h</b>) 07:00; (<b>i</b>) 08:00; (<b>j</b>) 09:00; (<b>k</b>) 10:00; (<b>l</b>) 11:00; (<b>m</b>) 12:00; (<b>n</b>) 13:00; (<b>o</b>) 14:00; (<b>p</b>) 15:00; (<b>q</b>) 16:00; (<b>r</b>) 17:00; (<b>s</b>) 18:00; (<b>t</b>) 19:00; (<b>u</b>) 20:00; (<b>v</b>) 21:00; (<b>w</b>) 22:00; (<b>x</b>) 23:00.</p>
Full article ">Figure 11
<p>Scatter plots of predicted LST results against observed Sentinel-3 SLSTR LST product: (<b>a</b>) 11:00 local time on 11 August 2019, (<b>b</b>) 23:00 local time on 11 August 2019, (<b>c</b>) 11:00 local time on 21 April 2022, and (<b>d</b>) 23:00 local time on 21 April 2022.</p>
Full article ">Figure 11 Cont.
<p>Scatter plots of predicted LST results against observed Sentinel-3 SLSTR LST product: (<b>a</b>) 11:00 local time on 11 August 2019, (<b>b</b>) 23:00 local time on 11 August 2019, (<b>c</b>) 11:00 local time on 21 April 2022, and (<b>d</b>) 23:00 local time on 21 April 2022.</p>
Full article ">Figure 12
<p>Distribution of LST errors between the predicted LST results and observed Sentinel-3 SLSTR LST product: (<b>a</b>) 11:00 local time on 11 August 2019, (<b>b</b>) 23:00 local time on 11 August 2019, (<b>c</b>) 11:00 local time on 21 April 2022, and (<b>d</b>) 23:00 local time on 21 April 2022.</p>
Full article ">Figure 13
<p>Comparison of the predicted hourly LSTs using the proposed framework and the predicted 100 m hourly LSTs using STARFM, FSDAF, and CFSDAF, respectively, with observed in situ LST using infrared thermometer on 11 August 2019. (<b>a</b>) subarea 1 (vegetation). (<b>b</b>) subarea 2 (water body). (<b>c</b>) subarea 3 (bare soil). (<b>d</b>) subarea 4 (ISA).</p>
Full article ">Figure 14
<p>(<b>a</b>) <span class="html-italic">R</span><sup>2</sup> and (<b>b</b>) RMSE between the predicted hourly LSTs using the proposed framework and the observed in situ LST using infrared thermometer on 11 August 2019.</p>
Full article ">Figure 15
<p>The generated 30 m hourly LST dataset using the proposed framework in Xi’an on 7 April 2022.</p>
Full article ">
16 pages, 6426 KiB  
Article
Unveiling Illumination Variations During a Lunar Eclipse: Multi-Wavelength Spaceborne Observations of the January 21, 2019 Event
by Min Shu, Tianyi Xu, Wei Cai, Shibo Wen, Hengyue Jiao and Yunzhao Wu
Remote Sens. 2024, 16(22), 4181; https://doi.org/10.3390/rs16224181 - 9 Nov 2024
Viewed by 821
Abstract
Space-based observations of the total lunar eclipse on 21 January 2019 were conducted using the geostationary Earth-orbiting satellite Gaofen-4 (GF-4). This study represents a pioneering effort to address the observational gap in full-disk lunar eclipse photometry from space. With its high resolution and [...] Read more.
Space-based observations of the total lunar eclipse on 21 January 2019 were conducted using the geostationary Earth-orbiting satellite Gaofen-4 (GF-4). This study represents a pioneering effort to address the observational gap in full-disk lunar eclipse photometry from space. With its high resolution and ability to capture the entire lunar disk, GF-4 enabled both quantitative and qualitative analyses of the variations in lunar brightness, as well as spectra and color changes, across two spatial dimensions, from the whole lunar disk to resolved regions. Our results indicate that before the totality phase of the lunar eclipse, the irradiance of the Moon diminishes to below approximately 0.19% of that of the uneclipsed Moon. Additionally, we observed an increase in lunar brightness at the initial entry into the penumbra. This phenomenon is attributed to the opposition effect, providing scientific evidence for this unexpected behavior. To investigate detailed spectral variations, specific calibration sites, including the Chang’E-3 landing site, MS-2 in Mare Serenitatis, and the Apollo 16 highlands, were analyzed. Notably, the red-to-blue ratio dropped below 1 near the umbra, contradicting the common perception that the Moon appears red during lunar eclipses. The red/blue ratio images reveal that as the Moon enters Earth’s umbra, it does not simply turn red; instead, a blue-banded ring appears at the boundary due to ozone absorption and the lunar surface composition. These findings significantly enhance our understanding of atmospheric effects on lunar eclipses and provide crucial reference information for the future modeling of lunar eclipse radiation, promoting the integration of remote sensing science with astronomy. Full article
(This article belongs to the Special Issue Laser and Optical Remote Sensing for Planetary Exploration)
Show Figures

Figure 1

Figure 1
<p>The effects of removing bad pixels and bad columns for GF-4 B2. (<b>a</b>) Before bad pixels removal; (<b>b</b>) After bad pixels removal; (<b>c</b>) before bad columns removal; (<b>d</b>) after bad columns removal.</p>
Full article ">Figure 2
<p>GF-4 B4 image mosaic (<b>Top</b>) and true color image mosaic (red: B4; green: B3; and blue: B2) (<b>Bottom</b>) before and after flat-field correction ((<b>Left</b>): before; (<b>Right</b>): after). The non-uniformity problems between the two stripe areas are significantly resolved.</p>
Full article ">Figure 3
<p>An overview of lunar radiation images obtained with a 30 ms exposure time during the lunar eclipse on 21 January 2019, presented in true color (red: B4; green: B3; and blue: B2). A 2% linear stretch was applied to these images for display enhancement to improve visibility.</p>
Full article ">Figure 4
<p>Disk-integrated irradiance at the standard distances during the lunar eclipse on 21 January 2019, measured by GF-4 across spectral bands B2–B5. Six sets of double-dotted lines depict each stage of the eclipse, denoted as P1–P4.</p>
Full article ">Figure 5
<p>Three sites in GF-4 color mosaic images captured at 02:30 UTC. (1) CE-3, (2) MS-2, and (3) Apollo-16 highlands. Due to the influence of observational geometry and fact that Site (3) is located in highlands, the brightness observed at site (3) is significantly higher than that of other sites. Consequently, a 2% linear stretch was specifically applied to Site (3) to enhance image contrast.</p>
Full article ">Figure 6
<p>The radiance spectra variation of CE-3 (<b>Top</b>), MS-2 (<b>Middle</b>) and Apollo 16 highlands (<b>Bottom</b>).</p>
Full article ">Figure 7
<p>Ratio of eclipsed irradiance to uneclipsed irradiance at corresponding phase angles over time, utilizing the lunar photometric model for GF-4 B2.</p>
Full article ">Figure 8
<p>Ratio images (654 nm/491 nm) from GF-4 data captured at 03:30 UTC, 03:40 UTC, 03:50 UTC, and 04:10 UTC on 21 January 2019.</p>
Full article ">
17 pages, 533 KiB  
Article
Statistical Analysis of LEO and GEO Satellite Anomalies and Space Radiation
by Jeimmy Nataly Buitrago-Leiva, Mohamed El Khayati Ramouz, Adriano Camps and Joan A. Ruiz-de-Azua
Aerospace 2024, 11(11), 924; https://doi.org/10.3390/aerospace11110924 - 8 Nov 2024
Viewed by 1108
Abstract
Exposure to space radiation substantially degrades satellite systems, provoking severe partial or, in some extreme cases, total failures. Electrostatic discharges (ESD), single event latch-up (SEL), and single event upsets (SEU) are among the most frequent causes of those reported satellite anomalies. The impact [...] Read more.
Exposure to space radiation substantially degrades satellite systems, provoking severe partial or, in some extreme cases, total failures. Electrostatic discharges (ESD), single event latch-up (SEL), and single event upsets (SEU) are among the most frequent causes of those reported satellite anomalies. The impact of space radiation dose on satellite equipment has been studied in-depth. This study conducts a statistical analysis to explore the relationships between low-Earth orbit (LEO) and geostationary orbit (GEO) satellite anomalies and particle concentrations, solar and geomagnetic activity in the period 2010–2022. Through a monthly and daily timescale analysis, the present work explores the temporal response of space disturbances on satellite systems and the periods when satellites are vulnerable to those disturbances. Full article
(This article belongs to the Section Astronautics & Space Science)
Show Figures

Figure 1

Figure 1
<p>Radiation particles and their effects on satellite systems.</p>
Full article ">Figure 2
<p>Anomaly selection criteria and their categorization per subsystem. F1, F2, F3 F4, and F5 corresponds to the filters described in <a href="#sec3dot1-aerospace-11-00924" class="html-sec">Section 3.1</a>.</p>
Full article ">Figure 3
<p>LEO/GEO satellites’ anomalies and solar and geomagnetic activity correlation (2010–2022). Sunspot number and CME speed index are used to quantify solar activity, whereas Kp and Dst indices are used to quantify geomagnetic activity. Despite using a monthly timescale for analysis, this figure is plotted by grouping four months for visual simplicity.</p>
Full article ">Figure 4
<p>Relationship between LEO anomalies and the average SSN, Kp, CME speed, and Dst indices of the month the anomaly occurred and the previous three months, respectively (<b>a</b>–<b>d</b>). Relationship between GEO anomalies and the average SSN, Kp, CME speed, and Dst indices of the month the anomaly occurred and the previous three months, respectively (<b>e</b>–<b>h</b>). Note: M-1,2,3 indicate the months before the anomalous event, respectively.</p>
Full article ">Figure 5
<p>Relationship between LEO (<b>a</b>–<b>c</b>) and GEO (<b>d</b>–<b>f</b>) anomalies and number of days per month with Kp-index ≥ 5 and Dst index ≤−16 of the month the anomaly occurred and the prior three months, respectively.</p>
Full article ">Figure 6
<p>Monthly LEO/GEO anomaly rate correlation with the average SSN, Kp, CME speed, and Dst indices, computed for every month by averaging the 13 samples throughout the study period (2010–2022).</p>
Full article ">Figure 7
<p>Assessment of solar and geomagnetic indicators during seven days (anomaly day + 6 prior days), for those days (anomaly day, D) with LEO and GEO anomalies ≥ 3.</p>
Full article ">Figure 8
<p>Proton, electron, and X-ray flux time series comparison between days with no reported anomalies (see <a href="#aerospace-11-00924-t002" class="html-table">Table 2</a>) and those selected days with <math display="inline"><semantics> <msub> <mi>A</mi> <mi>D</mi> </msub> </semantics></math> ≥ 3 during a seven-day window (see <a href="#aerospace-11-00924-t001" class="html-table">Table 1</a>). Given that one or a few days can have significantly higher particle concentrations than others, logarithmic charts are employed to respond to the large value ranges. Before the log() application, proton (P), electron (E), and X-ray flux (X) were in protons/(cm·day·sr), electrons/(cm·day·sr), and W/m<sup>2</sup>, respectively.</p>
Full article ">Figure 9
<p>Orbital inclination classification for LEO and GEO satellites whose anomalies were selected for the analysis. The satellites failures percentage in each orbit inclination range is normalized by dividing the total number of satellites failed by active satellites in that range from 2010 to 2022 according to Seradata [<a href="#B38-aerospace-11-00924" class="html-bibr">38</a>].</p>
Full article ">
19 pages, 6344 KiB  
Article
Evaluation of Fengyun-4B Satellite Temperature Profile Products Using Radiosonde Observations and ERA5 Reanalysis over Eastern Tibetan Plateau
by Yuhao Wang, Xiaofei Wu, Haoxin Zhang, Hong-Li Ren and Kaiqing Yang
Remote Sens. 2024, 16(22), 4155; https://doi.org/10.3390/rs16224155 - 7 Nov 2024
Viewed by 1010
Abstract
The latest-generation geostationary meteorological satellite, Fengyun-4B (FY-4B), equipped with the Geostationary Interferometric Infrared Sounder (GIIRS), offers high-spatiotemporal-resolution three-dimensional temperature structures. Its deployment serves as a critical complement to atmospheric temperature profile (ATP) observation in the Tibetan Plateau (TP). Based on radiosonde observation (RAOB) [...] Read more.
The latest-generation geostationary meteorological satellite, Fengyun-4B (FY-4B), equipped with the Geostationary Interferometric Infrared Sounder (GIIRS), offers high-spatiotemporal-resolution three-dimensional temperature structures. Its deployment serves as a critical complement to atmospheric temperature profile (ATP) observation in the Tibetan Plateau (TP). Based on radiosonde observation (RAOB) and the fifth-generation ECMWF global climate atmospheric reanalysis (ERA5), this study validates the availability and representativeness of FY-4B/GIIRS ATP products in the eastern TP region. Due to the issue of satellite zenith, this study focuses solely on examining the eastern TP region. Under a clear sky, FY-4B/GIIRS ATP exhibits good consistency with RAOB compared to cloudy conditions, with an average root mean square error (RMSE) of 2.57 K. FY-4B/GIIRS tends to underestimate temperatures in the lower layers while overestimating temperatures in the upper layers. The bias varies across seasons. Except for summer, the horizontal and vertical bias distribution patterns are similar, though there are slight differences in values. Despite the presence of bias, FY-4B/GIIRS ATP maintains a good consistency with observations and reanalysis data, indicating commendable product quality. These results demonstrate that it can play a vital role in augmenting the ATP observation network limited by sparse radiosonde stations in the eastern TP, offering crucial data support for numerical weather prediction, weather monitoring, and related meteorological research in this region. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Distribution map of the nine RAOB stations (red triangles) over the TP. (<b>b</b>) The FY-4B/GIIRS observation pixels (blue dots) for the Garze station in the MW method at 12 UTC on 17 January 2023. The color shading represents the elevation (units, m), and the red line in (<b>a</b>) indicates the border of the TP.</p>
Full article ">Figure 2
<p>RMSE (green bars) and the number of effective data (orange bars) for the IDW and the MW method at nine RAOB stations.</p>
Full article ">Figure 3
<p>The percentages of the FY-4B/GIIRS ATP products quality flags during clear sky (green bars) and cloudy sky (orange bars).</p>
Full article ">Figure 4
<p>The average ATP observed by FY-4B/GIIRS (blue line) and RAOB (orange line) and the average bias of FY-4B/GIIRS referring to RAOB (cyan line with triangles) for (<b>a</b>–<b>i</b>) 00 UTC and (<b>j</b>–<b>r</b>) 12 UTC. The light cyan shading accompanied with the bias line indicates one standard variation of the bias.</p>
Full article ">Figure 5
<p>Scatter plot of FY-4B/GIIRS ATP versus the RAOB ATP (black dashed line represents the 1:1 line, red line represents regression line). (<b>a</b>–<b>i</b>) represent nine RAOB stations arranged in order of elevation from lowest to highest.</p>
Full article ">Figure 6
<p>Same as <a href="#remotesensing-16-04155-f005" class="html-fig">Figure 5</a>, but for ERA5 ATP versus the RAOB ATP. (<b>a</b>–<b>i</b>) represent nine RAOB stations arranged in order of elevation from lowest to highest.</p>
Full article ">Figure 7
<p>(<b>a</b>) Satellite zenith angle (shaded, degree) of FY-4B/GIIRS at 11:00 UTC on 17 January 2024 and the annual mean troposphere temperature (<b>b</b>) before and (<b>c</b>) after filtering based on the satellite zenith angle of 60° as the red line shown in (<b>b</b>). The black line in (<b>a</b>–<b>c</b>) indicates the TP region. Points A and B in (<b>b</b>) are the intersection points of the contour line of 60° and the borderline of the TP region in (<b>a</b>).</p>
Full article ">Figure 8
<p>The spatial distribution of annual mean temperature bias between FY-4B/GIIRS and ERA5 ATP: (<b>a</b>) horizontal distribution of troposphere (600–100 hPa) averaged bias and (<b>b</b>) vertical distribution of regional averaged bias for the blue box in (<b>a</b>). The shading in (<b>b</b>) indicates one STD range of the bias.</p>
Full article ">Figure 9
<p>Scatter plot of seasonal average FY-4B/GIIRS ATP versus ERA5 ATP for each of the four seasons among the eastern TP, the black dashed line represents the 1:1 line, and the red line represents the regression line. (<b>a</b>–<b>d</b>) correspond to winter, spring, summer, and autumn, respectively.</p>
Full article ">Figure 10
<p>The horizontal distribution of annual mean troposphere (600–100 hPa) averaged temperature bias between FY-4B/GIIRS and ERA5 ATP. (<b>a</b>–<b>d</b>) correspond to winter, spring, summer, and autumn, respectively.</p>
Full article ">Figure 11
<p>The vertical distribution of annual mean regional averaged temperature bias between FY-4B/GIIRS and ERA5 ATP for the blue box in <a href="#remotesensing-16-04155-f008" class="html-fig">Figure 8</a>a. The shading indicates one STD range of the bias. (<b>a</b>–<b>d</b>) correspond to winter, spring, summer, and autumn, respectively.</p>
Full article ">
Back to TopTop