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

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (833)

Search Parameters:
Keywords = geostationary

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
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 309
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 713
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 421
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 617
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 502
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 815
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 735
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 865
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 906
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 ">
20 pages, 2978 KiB  
Article
Considerations for a Micromirror Array Optimized for Compressive Sensing (VIS to MIR) in Space Applications
by Ulrike Dauderstädt, Peter Dürr, Detlef Kunze, Sara Francés González, Donato Borrelli, Lorenzo Palombi, Valentina Raimondi and Michael Wagner
J. Imaging 2024, 10(11), 282; https://doi.org/10.3390/jimaging10110282 - 5 Nov 2024
Viewed by 878
Abstract
Earth observation (EO) is crucial for addressing environmental and societal challenges, but it struggles with revisit times and spatial resolution. The EU-funded SURPRISE project aims to improve EO capabilities by studying space instrumentation using compressive sensing (CS) implemented through spatial light modulators (SLMs) [...] Read more.
Earth observation (EO) is crucial for addressing environmental and societal challenges, but it struggles with revisit times and spatial resolution. The EU-funded SURPRISE project aims to improve EO capabilities by studying space instrumentation using compressive sensing (CS) implemented through spatial light modulators (SLMs) based on micromirror arrays (MMAs) to improve the ground sampling distance. In the SURPRISE project, we studied the development of an MMA that meets the requirements of a CS-based geostationary instrument working in the visible (VIS) and mid-infrared (MIR) spectral ranges. This paper describes the optical simulation procedure and the results obtained for analyzing the performance of such an MMA with the goal of identifying a mirror design that would allow the device to meet the optical requirements of this specific application. Full article
Show Figures

Figure 1

Figure 1
<p>Compressive sensing principle.</p>
Full article ">Figure 2
<p>Schematic illustration of the optical principle.</p>
Full article ">Figure 3
<p>Mirror actuation modes.</p>
Full article ">Figure 4
<p>Mirror deflection vs. addressing voltage.</p>
Full article ">Figure 5
<p>On-axis illumination–angles.</p>
Full article ">Figure 6
<p>Diagonal vs. orthogonal mirror.</p>
Full article ">Figure 7
<p>Airy disks for different wavelengths (amplitude, <math display="inline"><semantics> <msub> <mi>E</mi> <mi>Airy</mi> </msub> </semantics></math>). The mirror size in this example is <math display="inline"><semantics> <mrow> <mn>20</mn> <mo> </mo> <mo>μ</mo> <mi mathvariant="normal">m</mi> </mrow> </semantics></math>, as indicated by the gridlines.</p>
Full article ">Figure 8
<p>Optical signal flow.</p>
Full article ">Figure 9
<p>Intensity in the Fourier plane, <math display="inline"><semantics> <mrow> <msub> <mi>I</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mrow> <mo>|</mo> <msub> <mi>E</mi> <mn>2</mn> </msub> <mo>|</mo> </mrow> <mn>2</mn> </msup> </mrow> </semantics></math>, for one pixel in the ‘ON’ state surrounded by pixels in the ‘OFF’ state with the configurations from <a href="#jimaging-10-00282-t003" class="html-table">Table 3</a>. The red circles indicate the area seen via the collimator.</p>
Full article ">Figure 10
<p>Efficiency for ‘ON’ and ‘OFF’ pixels, minimum and maximum mirror size (<a href="#jimaging-10-00282-t003" class="html-table">Table 3</a>), all parameters in <a href="#jimaging-10-00282-t004" class="html-table">Table 4</a>.</p>
Full article ">Figure 11
<p>Intensity, <math display="inline"><semantics> <mrow> <msub> <mi>I</mi> <mn>2</mn> </msub> <mo>=</mo> <msup> <mrow> <mo>|</mo> <msub> <mi>E</mi> <mn>2</mn> </msub> <mo>|</mo> </mrow> <mn>2</mn> </msup> </mrow> </semantics></math> in collimator plane for ‘OFF’ pixel with <math display="inline"><semantics> <mrow> <mi>p</mi> <mo>=</mo> <msub> <mi>p</mi> <mi>min</mi> </msub> </mrow> </semantics></math>, diagonal mirror.</p>
Full article ">Figure 12
<p>Dependency on deflection angle, <math display="inline"><semantics> <mrow> <mi>p</mi> <mo>=</mo> <mn>30</mn> <mo> </mo> <mrow> <mo>μ</mo> <mi mathvariant="normal">m</mi> </mrow> </mrow> </semantics></math>, diagonal mirrors, all parameters in <a href="#jimaging-10-00282-t004" class="html-table">Table 4</a>.</p>
Full article ">Figure 13
<p>Dependency on mirror size, <math display="inline"><semantics> <msub> <mi>N</mi> <mi>collect</mi> </msub> </semantics></math>, <math display="inline"><semantics> <mi>δ</mi> </semantics></math> adjusted for <span class="html-italic">p</span>, diagonal mirrors, all parameters in <a href="#jimaging-10-00282-t004" class="html-table">Table 4</a>.</p>
Full article ">Figure 14
<p>Dependency on the mirror size, <math display="inline"><semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>10.35</mn> <mo>°</mo> </mrow> </semantics></math>, <math display="inline"><semantics> <msub> <mi>N</mi> <mi>collect</mi> </msub> </semantics></math> adjusted for <span class="html-italic">p</span>, diagonal mirrors, all parameters in <a href="#jimaging-10-00282-t004" class="html-table">Table 4</a>.</p>
Full article ">Figure 15
<p>Dependency on <math display="inline"><semantics> <msub> <mi>N</mi> <mi>collect</mi> </msub> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>6</mn> <mo>°</mo> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>p</mi> <mo>=</mo> <mn>30</mn> <mo> </mo> <mrow> <mo>μ</mo> <mi mathvariant="normal">m</mi> </mrow> </mrow> </semantics></math>, diagonal mirrors, all parameters in <a href="#jimaging-10-00282-t004" class="html-table">Table 4</a>.</p>
Full article ">Figure 16
<p>Dependency on mirror size, <math display="inline"><semantics> <mrow> <msub> <mi>f</mi> <mi>collect</mi> </msub> <mo>=</mo> <mn>6.17</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>6</mn> <mo>°</mo> </mrow> </semantics></math>, diagonal mirrors, all parameters in <a href="#jimaging-10-00282-t004" class="html-table">Table 4</a>.</p>
Full article ">Figure 17
<p>Diffraction efficiency for different micropixel sizes, <math display="inline"><semantics> <mrow> <msub> <mi>N</mi> <mi>collect</mi> </msub> <mo>=</mo> <mrow> <mn>2.35</mn> </mrow> </mrow> </semantics></math>.</p>
Full article ">Figure 18
<p>Efficiencies for modified system specifications (<math display="inline"><semantics> <mrow> <mi>G</mi> <mi>S</mi> <mi>D</mi> </mrow> </semantics></math>) as in <a href="#jimaging-10-00282-t005" class="html-table">Table 5</a>.</p>
Full article ">
19 pages, 1475 KiB  
Review
Overview of Space-Based Laser Communication Missions and Payloads: Insights from the Autonomous Laser Inter-Satellite Gigabit Network (ALIGN)
by Othman I. Younus, Amna Riaz, Richard Binns, Eamon Scullion, Robert Wicks, Jethro Vernon, Chris Graham, David Bramall, Jurgen Schmoll and Cyril Bourgenot
Aerospace 2024, 11(11), 907; https://doi.org/10.3390/aerospace11110907 - 5 Nov 2024
Viewed by 2382
Abstract
This paper examines the growing adoption of laser communication (lasercom) in space missions and payloads for identifying emerging trends and key technology drivers of future optical communications satellite systems. It also presents a comprehensive overview of commercially available and custom-designed lasercom terminals, outlining [...] Read more.
This paper examines the growing adoption of laser communication (lasercom) in space missions and payloads for identifying emerging trends and key technology drivers of future optical communications satellite systems. It also presents a comprehensive overview of commercially available and custom-designed lasercom terminals, outlining their characteristics and specifications to meet the evolving demands of global satellite networks. The analysis explores the technical considerations and challenges associated with integrating lasercom terminals into LEO constellations and the Inter-satellite communications service provision in LEO due to their power, size, and weight constraints. By analyzing advancements in CubeSat lasercom technology designed to cater for the emergence of future mega constellations of interacting small satellites, the paper underscores its promising role in establishing high-performance satellite communication networks for future space exploration and data transmission. In addition, a brief overview of our ALIGN planned mission is provided, which highlights the main key operational features in terms of PAT and link budget analysis. Full article
(This article belongs to the Special Issue Space Telescopes & Payloads)
Show Figures

Figure 1

Figure 1
<p>Architecture of satellite communication systems via different communication link types.</p>
Full article ">Figure 2
<p>Distribution of link types in satellite communication systems by country.</p>
Full article ">Figure 3
<p>Current design of FOCUS.</p>
Full article ">Figure 4
<p>Mission concept in space, showing positional uncertainties (yellow circles) along the dashed path, with Tx and Rx alignment in green and blue cones.</p>
Full article ">
18 pages, 10136 KiB  
Article
The Combination Application of FY-4 Satellite Products on Typhoon Saola Forecast on the Sea
by Chun Yang, Bingying Shi and Jinzhong Min
Remote Sens. 2024, 16(21), 4105; https://doi.org/10.3390/rs16214105 - 2 Nov 2024
Viewed by 967
Abstract
Satellite data play an irreplaceable role in global observation data systems. Effective comprehensive application of satellite products will inevitably improve numerical weather prediction. FengYun-4 (FY-4) series satellites can provide not only radiance data but also retrieval data with high temporal and spatial resolutions. [...] Read more.
Satellite data play an irreplaceable role in global observation data systems. Effective comprehensive application of satellite products will inevitably improve numerical weather prediction. FengYun-4 (FY-4) series satellites can provide not only radiance data but also retrieval data with high temporal and spatial resolutions. To evaluate the potential benefits of the combination application of FY-4 Advanced Geostationary Radiance Imager (AGRI) products on Typhoon Saola analysis and forecast, two group of experiments are set up with the Weather Research and Forecasting model (WRF). Compared with the benchmark experiment, whose sea surface temperature (SST) is from the National Centers for Environmental Prediction (NCEP) reanalysis data, the SST replacement experiments with FY-4 A/B SST products significantly improve the track and precipitation forecast, especially with the FY-4B SST product. Based on the above results, AGRI clear-sky and all-sky assimilations with FY-4B SST are implemented with a self-constructed AGRI assimilation module. The results show that the AGRI all-sky assimilation experiment can obtain better analyses and forecasts. Furthermore, it is proven that the combination application of AGRI radiance and SST products is beneficial for typhoon prediction. Full article
Show Figures

Figure 1

Figure 1
<p>The weighting function of channels 9-14 of FY-4A AGRI with RTTOV and the U.S. standard atmospheric profile.</p>
Full article ">Figure 2
<p>(<b>a</b>) The evolution of the best track, (<b>b</b>) the central sea level pressure (units: hPa) and maximum wind (units: knot) for Typhoon Saola from 0000 UTC 22 August to 1200 UTC 3 September 2023.</p>
Full article ">Figure 3
<p>Initial SST (units: K) from (<b>a</b>) <span class="html-italic">CON</span>, (<b>b</b>) SSTA, and (<b>c</b>) SSTB.</p>
Full article ">Figure 4
<p>The predicted (<b>a</b>) track, (<b>b</b>) track errors (units: km), (<b>c</b>) CSLP errors (units: hPa), and (<b>d</b>) MW errors (units: knot) in <span class="html-italic">CON</span> (light blue lines), SSTA (red lines), and SSTB (light green lines) are compared to the JMA best track estimates (blue lines) from 0600 UTC 30 August to 1200 UTC 2 September 2023.</p>
Full article ">Figure 5
<p>Time series of the U and V components of average steering flow (units: m/s) from 0600 UTC 30 August to 1200 UTC 2 September 2023.</p>
Full article ">Figure 6
<p>The 24 h accumulated precipitation (units: mm) from 1200 UTC 1 September to 1200 UTC 2 September 2023 of (<b>a</b>) the Micaps observation; (<b>b</b>) the interpolated Micaps observation with a horizontal resolution of 0.5° × 0.5°; (<b>c</b>) <span class="html-italic">CON</span>; (<b>d</b>) SSTA; and (<b>e</b>) SSTB. The dots with different colors in (<b>a</b>) represent different accumulated precipitation, as shown in the color bar.</p>
Full article ">Figure 7
<p>Performance diagram for the 24 h accumulated precipitation for the <span class="html-italic">CON</span> (light blue), SSTA (red), and SSTB (light green) with a threshold of (<b>a</b>) 0.01 mm; (<b>b</b>) 10 mm; (<b>c</b>) 25 mm; (<b>d</b>) 50 mm; and (<b>e</b>) 75 mm from 1200 UTC 1 September to 1200 UTC 2 September 2023.</p>
Full article ">Figure 8
<p>(<b>a</b>,<b>b</b>) The AGRI observed brightness temperature (units: K) distributions at channel 9 after QC in (<b>a</b>) CLR and (<b>b</b>) ALL valid at 1500 UTC 30 August 2023. (<b>c</b>) The counts of assimilated AGRI observations at channel 9 in ALL and CLR with different cloud mask types every 3 hr from 0900 UTC 30 August to 1500 UTC 30 August 2023.</p>
Full article ">Figure 9
<p>The IPs (units: %) over <span class="html-italic">CON</span> of individual experiments every 3 h from 0900 UTC 30 August to 1500 UTC 30 August 2023 in (<b>a</b>) <span class="html-italic">CTTs</span> and (<b>b</b>) agreements on sky conditions.</p>
Full article ">Figure 10
<p>The predicted (<b>a</b>) track, (<b>b</b>) track errors (units: km), (<b>c</b>) CSLP errors (units: hPa), and (<b>d</b>) MW errors (units: knot) in <span class="html-italic">CON</span> (light blue lines), SSTB (light green lines), CLR (light yellow lines), ALL (orange lines), CLR + SSTB (light red lines), and ALL + SSTB (brown lines) are compared to the JMA best track estimates (blue lines) from 1800 UTC 30 August to 1200 UTC 2 September 2023.</p>
Full article ">Figure 11
<p>The (<b>a</b>) U and (<b>b</b>) V components of steering flows (units: m/s) from 700 to 200 hPa with an interval of 50 hPa in individual experiments at 0600 UTC 1 September 2023.</p>
Full article ">Figure 12
<p>The same as <a href="#remotesensing-16-04105-f007" class="html-fig">Figure 7</a> but for (<b>a</b>) the interpolated Micaps observation with a horizontal resolution of 0.5° × 0.5°; (<b>b</b>) <span class="html-italic">CON</span>; (<b>c</b>) CLR; (<b>d</b>) ALL; (<b>e</b>) SSTB; (<b>f</b>) CLR + SSTB; and (<b>g</b>) ALL + SSTB.</p>
Full article ">Figure 13
<p>The same as <a href="#remotesensing-16-04105-f006" class="html-fig">Figure 6</a> but for <span class="html-italic">CON</span> (light blue), SSTB (light green), CLR (light yellow), ALL (orange), CLR + SSTB (light red), and ALL + SSTB (brown) with a threshold of (<b>a</b>) 0.01 mm; (<b>b</b>) 10 mm; (<b>c</b>) 25 mm; (<b>d</b>) 50 mm; and (<b>e</b>) 75 mm.</p>
Full article ">
24 pages, 7524 KiB  
Article
A Study on Typhoon Center Localization Based on an Improved Spatio-Temporally Consistent Scale-Invariant Feature Transform and Brightness Temperature Perturbations
by Chaoyu Yan, Jie Guang, Zhengqiang Li, Gerrit de Leeuw and Zhenting Chen
Remote Sens. 2024, 16(21), 4070; https://doi.org/10.3390/rs16214070 - 31 Oct 2024
Viewed by 910
Abstract
Extreme weather events like typhoons have become more frequent due to global climate change. Current typhoon monitoring methods include manual monitoring, mathematical morphological methods, and artificial intelligence. Manual monitoring is accurate but labor-intensive, while AI offers convenience but requires accuracy improvements. Mathematical morphology [...] Read more.
Extreme weather events like typhoons have become more frequent due to global climate change. Current typhoon monitoring methods include manual monitoring, mathematical morphological methods, and artificial intelligence. Manual monitoring is accurate but labor-intensive, while AI offers convenience but requires accuracy improvements. Mathematical morphology methods, such as brightness temperature perturbation (BTP) and a spatio-temporally consistent (STC) Scale-Invariant Feature Transform (SIFT), remain mainstream for typhoon positioning. This paper enhances BTP and STC SIFT methods for application to Fengyun 4A (FY-4A) Advanced Geosynchronous Radiation Imager (AGRI) L1 data, incorporating parallax correction for more accurate surface longitude and latitude positioning. The applicability of these methods for different typhoon intensities and monitoring time resolutions is analyzed. Automated monitoring with one-hour observation intervals in the northwest Pacific region demonstrates high positioning accuracy, reaching 25 km or better when compared to best path data from the China Meteorological Administration (CMA). For 1 h remote sensing observations, BTP is more accurate for typhoons at or above typhoon intensity, while STC SIFT is more accurate for weaker typhoons. In the current era of a high temporal resolution of typhoon monitoring using geostationary satellites, the method presented in this paper can serve the national meteorological industry for typhoon monitoring, which is beneficial to national pre-disaster prevention work as well as global meteorological research. Full article
Show Figures

Figure 1

Figure 1
<p>The yellow rectangular box represents the study area.</p>
Full article ">Figure 2
<p>Typhoon automatic center positioning process.</p>
Full article ">Figure 3
<p>Parallax correction geometry relationship model [<a href="#B32-remotesensing-16-04070" class="html-bibr">32</a>].</p>
Full article ">Figure 4
<p>The application of the typhoon detection algorithm to identify the typhoon cloud system over the northwest Pacific using FY-4A AGRI level 1 data on 6 November 2019, at 02:00 Beijing time (typhoon HL). (<b>a</b>) The initial data after projection conversion pre-processing, (<b>b</b>) the spatial distribution of the BT over the study area (unit: K), (<b>c</b>) the binarized image of the BT spatial distribution, and (<b>d</b>) further processing that shows the locations of candidate target cloud systems.</p>
Full article ">Figure 4 Cont.
<p>The application of the typhoon detection algorithm to identify the typhoon cloud system over the northwest Pacific using FY-4A AGRI level 1 data on 6 November 2019, at 02:00 Beijing time (typhoon HL). (<b>a</b>) The initial data after projection conversion pre-processing, (<b>b</b>) the spatial distribution of the BT over the study area (unit: K), (<b>c</b>) the binarized image of the BT spatial distribution, and (<b>d</b>) further processing that shows the locations of candidate target cloud systems.</p>
Full article ">Figure 5
<p>The BT distribution of the typhoon cloud system using FY-4A AGRI level 1 data on 6 November 2019, at 02:00 Beijing time: the BT distribution of the target cloud system before (<b>a</b>) and after (<b>b</b>) parallax correction. (<b>c</b>) and (<b>d</b>) show details of the typhoon eye area before and after parallax correction, respectively.</p>
Full article ">Figure 6
<p>Results from the BTP typhoon localization algorithm applied to FY-4A AGRI level 1 data over the northwest Pacific on 6 November 2019, at 02:00 Beijing time (typhoon HL): (<b>a</b>) the longitudinal BT gradient near the typhoon eye area (unit: km); (<b>b</b>) latitudinal BT gradient near the typhoon eye area; (<b>c</b>) spatial distribution of BT divergence (unit: <math display="inline"><semantics> <mrow> <mi mathvariant="normal">N</mi> <mo>·</mo> <msup> <mrow> <mi mathvariant="normal">m</mi> </mrow> <mrow> <mo>−</mo> <mn>3</mn> </mrow> </msup> <mo>×</mo> <msup> <mrow> <mn>10</mn> </mrow> <mrow> <mn>7</mn> </mrow> </msup> </mrow> </semantics></math>); (<b>d</b>) spatial distribution of the BT curl (unit: <math display="inline"><semantics> <mrow> <mi mathvariant="normal">N</mi> <mo>·</mo> <msup> <mrow> <mi mathvariant="normal">m</mi> </mrow> <mrow> <mo>−</mo> <mn>3</mn> </mrow> </msup> <mo>×</mo> <msup> <mrow> <mn>10</mn> </mrow> <mrow> <mn>7</mn> </mrow> </msup> </mrow> </semantics></math>); (<b>e</b>) spatial distribution of BTP. The location of the typhoon as determined from the BTP distribution is indicated with a red + and the location of the typhoon center provided by CMA is indicated with a blue star.</p>
Full article ">Figure 6 Cont.
<p>Results from the BTP typhoon localization algorithm applied to FY-4A AGRI level 1 data over the northwest Pacific on 6 November 2019, at 02:00 Beijing time (typhoon HL): (<b>a</b>) the longitudinal BT gradient near the typhoon eye area (unit: km); (<b>b</b>) latitudinal BT gradient near the typhoon eye area; (<b>c</b>) spatial distribution of BT divergence (unit: <math display="inline"><semantics> <mrow> <mi mathvariant="normal">N</mi> <mo>·</mo> <msup> <mrow> <mi mathvariant="normal">m</mi> </mrow> <mrow> <mo>−</mo> <mn>3</mn> </mrow> </msup> <mo>×</mo> <msup> <mrow> <mn>10</mn> </mrow> <mrow> <mn>7</mn> </mrow> </msup> </mrow> </semantics></math>); (<b>d</b>) spatial distribution of the BT curl (unit: <math display="inline"><semantics> <mrow> <mi mathvariant="normal">N</mi> <mo>·</mo> <msup> <mrow> <mi mathvariant="normal">m</mi> </mrow> <mrow> <mo>−</mo> <mn>3</mn> </mrow> </msup> <mo>×</mo> <msup> <mrow> <mn>10</mn> </mrow> <mrow> <mn>7</mn> </mrow> </msup> </mrow> </semantics></math>); (<b>e</b>) spatial distribution of BTP. The location of the typhoon as determined from the BTP distribution is indicated with a red + and the location of the typhoon center provided by CMA is indicated with a blue star.</p>
Full article ">Figure 7
<p>Typhoon center localization with the STC SIFT feature method using FY-4A AGRI level 1 data on 6 November 2019, at 02:00 Beijing time (typhoon HL): (<b>a</b>) feature point distributions in the extracted historical image (<b>left</b>) and the current image (<b>right</b>); (<b>b</b>) results of matching the remaining feature points after STC filtering and rotation uniform distribution filtering; (<b>c</b>) the comparison of the final positioning result (red +) with the typhoon center location provided by CMA (blue star).</p>
Full article ">Figure 8
<p>The comparison of the optimal paths of SD, HMN, HL, and HK determined using two typhoon localization methods, with the best path provided by CMA. The left column shows results from the BTP typhoon localization method and the right column shows results from the STC SIFT feature typhoon localization method for typhoons SD (<b>a</b>,<b>b</b>), HNM (<b>c</b>,<b>d</b>), HL (<b>e</b>,<b>f</b>), and HK (<b>g</b>,<b>h</b>). The average accuracy, from comparison with the CMA path, is indicated in each figure. The region outlined by blue lines delineates the portion of the typhoon characterized by lower intensity (before developing into typhoon intensity, TY). <a href="#remotesensing-16-04070-t004" class="html-table">Table 4</a> presents the typhoon localization accuracy for the blue-outlined area.</p>
Full article ">Figure 9
<p>The error analysis of different typhoon localization methods under various typhoon intensities. Panel (<b>a</b>) illustrates the error analysis for STC SIFT at different typhoon intensities; panel (<b>b</b>) shows the error analysis for BTP typhoon localization under varying intensities; and panel (<b>c</b>) depicts the error analysis for BTP localization without parallax correction across different typhoon intensities. The orange line represents the median error, while the green dashed line indicates the mean error. The top and bottom edges of each box correspond to the upper and lower quartiles of the error distribution, and the whiskers denote the maximum and minimum error values.</p>
Full article ">
20 pages, 13089 KiB  
Article
Development of Vertical Radar Reflectivity Profiles Based on Lightning Density Using the Geostationary Lightning Mapper Dataset in the Subtropical Region of Brazil
by Tiago Bentes Mandú, Laurizio Emanuel Ribeiro Alves, Éder Paulo Vendrasco and Thiago Souza Biscaro
Remote Sens. 2024, 16(20), 3767; https://doi.org/10.3390/rs16203767 - 11 Oct 2024
Viewed by 812
Abstract
The study aims to develop vertical radar reflectivity profiles based on lightning density data from the Geostationary Lightning Mapper (GLM) on the GOES-16 satellite in the subtropical region of Brazil. The primary objective is to improve the assimilation of lightning data in numerical [...] Read more.
The study aims to develop vertical radar reflectivity profiles based on lightning density data from the Geostationary Lightning Mapper (GLM) on the GOES-16 satellite in the subtropical region of Brazil. The primary objective is to improve the assimilation of lightning data in numerical weather prediction models. The methodology involves the analysis of polarimetric radar data from Chapecó-SC and Jaraguari-MS, spanning from January 2019 to December 2023, and their correlation with lightning data from the GLM. Radar reflectivity profiles were created for different lightning density classes, categorized into six classes based on geometric progression. Results show a significant relationship between lightning activity and radar reflectivity, with distinct profiles for convective and stratiform events. These findings demonstrate the potential of using GLM data to enhance short-term weather forecasting, particularly for severe weather events. The study concludes that the integration of GLM data into weather models can lead to more accurate predictions of intense precipitation events, contributing to better preparedness and response strategies. Full article
(This article belongs to the Special Issue Remote Sensing of Extreme Weather Events: Monitoring and Modeling)
Show Figures

Figure 1

Figure 1
<p>Study area.</p>
Full article ">Figure 2
<p>Diagram illustrating the relationship between a single GLM pixel and the corresponding 9 × 9 radar grid points used for data matching.</p>
Full article ">Figure 3
<p>Decision tree schematic illustrating the logic used to determine if a reflectivity bin is classified as convective or stratiform. Source: <a href="https://vlab.noaa.gov/web/wdtd/-/convective-stratiform-precipitation-separation-csps-algorithm" target="_blank">https://vlab.noaa.gov/web/wdtd/-/convective-stratiform-precipitation-separation-csps-algorithm</a> (accessed on 26 August 2024).</p>
Full article ">Figure 4
<p>Vertical radar reflectivity profiles based on lightning density classes.</p>
Full article ">Figure 5
<p>Percentage of data used for composing the average profile by lightning density class.</p>
Full article ">Figure 6
<p>Number of radar profiles by height level and lightning density class.</p>
Full article ">Figure 7
<p>Sum of stratiform vertical profiles by lightning density class.</p>
Full article ">Figure 8
<p>Sum of convective vertical profiles by lightning density class.</p>
Full article ">Figure 9
<p>Vertically Integrated Liquid (VIL) for stratiform and convective profiles.</p>
Full article ">Figure 10
<p>Map of streamlines at different atmospheric levels (1000, 850, 500, and 250 hPa), Mean Sea Level Pressure (MSLP), and Geopotential Height (GH) on 29 November 2020, at 12:00 UTC over South America.</p>
Full article ">Figure 11
<p>Map of the brightness temperature from Band 13 of the GOES-16 satellite on 29 November 2020, at 12 UTC over South America.</p>
Full article ">Figure 12
<p>Map of the thermodynamic indices CAPE and LI on 29 November 2020, at 12 UTC over South America.</p>
Full article ">Figure 13
<p>Map of the CAPPI composition (dBZ) at 3 km, 5 km, and 7 km from the Chapecó-SC radar and Lightning Density (UNIT) on 29 November 2020, at 12 UTC over South America.</p>
Full article ">
18 pages, 12939 KiB  
Article
Dust Monitoring and Three-Dimensional Transport Characteristics of Dust Aerosol in Beijing, Tianjin, and Hebei
by Siqin Zhang, Jianjun Wu, Jiaqi Yao, Xuefeng Quan, Haoran Zhai, Qingkai Lu, Haobin Xia, Mengran Wang and Jinquan Guo
Atmosphere 2024, 15(10), 1212; https://doi.org/10.3390/atmos15101212 - 10 Oct 2024
Viewed by 726
Abstract
Global dust events have become more frequent due to climate change and increased human activity, significantly impacting air quality and human health. Previous studies have mainly focused on determining atmospheric dust pollution levels through atmospheric parameter simulations or AOD values obtained from satellite [...] Read more.
Global dust events have become more frequent due to climate change and increased human activity, significantly impacting air quality and human health. Previous studies have mainly focused on determining atmospheric dust pollution levels through atmospheric parameter simulations or AOD values obtained from satellite remote sensing. However, research on the quantitative description of dust intensity and its cross-regional transport characteristics still faces numerous challenges. Therefore, this study utilized Fengyun-4A (FY-4A) satellite Advanced Geostationary Radiation Imager (AGRI) imagery, Cloud-Aerosol Lidar, and Infrared Pathfinder Satellite Observation (CALIPSO) lidar, and other auxiliary data, to conduct three-dimensional spatiotemporal monitoring and a cross-regional transport analysis of two typical dust events in the Beijing–Tianjin–Hebei (BTH) region of China using four dust intensity indices Infrared Channel Shortwave Dust (Icsd), Dust Detection Index (DDI), dust value (DV), and Dust Strength Index (DSI)) and the HYSPLIT model. We found that among the four indices, DDI was the most suitable for studying dust in the BTH region, with a detection accuracy (POCD) of >88% at all times and reaching a maximum of 96.14%. Both the 2021 and 2023 dust events originated from large-scale deforestation in southern Mongolia and the border area of Inner Mongolia, with dust plumes distributed between 2 and 12 km being transported across regions to the BTH area. Further, when dust aerosols are primarily concentrated below 4 km and PM10 concentrations consistently exceed 600 µg/m3, large dust storms are more likely to occur in the BTH region. The findings of this study provide valuable insights into the sources, transport pathways, and environmental impacts of dust aerosols. Full article
(This article belongs to the Section Aerosols)
Show Figures

Figure 1

Figure 1
<p>Administrative map. (<b>a</b>) National 1 km DEM elevation map. (<b>b</b>) PM<sub>10</sub> monitoring station distribution in Beijing-Tianjin-Hebei Region. (<b>c</b>) Bar chart of dust source management project construction in Beijing-Tianjin-Hebei Region (2015–2019).</p>
Full article ">Figure 2
<p>Technical flowchart.</p>
Full article ">Figure 3
<p>Histogram of frequency distribution for thin clouds, thick clouds, and dust under four dust intensity indices.</p>
Full article ">Figure 4
<p>Dust identification results in the Beijing–Tianjin–Hebei Region. AGRI true-color images for 15 March 2021, UTC 03:00–06:00 (<b>a<sub>1</sub></b>–<b>a<sub>4</sub></b>), and DDI distribution maps (<b>b<sub>1</sub></b>–<b>b<sub>4</sub></b>); AGRI true-color images for 22 March 2023, UTC 03:00–06:00 (<b>c<sub>1</sub></b>–<b>c<sub>4</sub></b>), and DDI distribution maps (<b>d<sub>1</sub></b>–<b>d<sub>4</sub></b>); DDI violin and boxplot statistics for 15 March 2021, and 22 March 2023, UTC 03:00–06:00 (<b>e</b>).</p>
Full article ">Figure 5
<p>HYSPLIT backward trajectory simulations and FY-4A true-color images for the two dust events: (<b>a</b>,<b>b</b>) Beijing backward trajectory simulation for 15 March 2021; (<b>d</b>,<b>e</b>) Beijing backward trajectory simulation for 22 March 2023; (<b>c</b>) an FY-4A true-color image for 15 March 2021, at UTC 04:00; (<b>f</b>) an FY-4A true-color image for 22 March 2023, at UTC 04:00.</p>
Full article ">Figure 6
<p>Vertical distribution characteristics of aerosols and hourly changes in PM<sub>10</sub> concentration in the BTH and Inner Mongolia regions: 15 March 2021, BTH and Inner Mongolia regions (<b>a<sub>1</sub></b>–<b>a<sub>3</sub></b>, <b>b<sub>1</sub></b>–<b>b<sub>3</sub></b>); 21 March 2023, BTH and Inner Mongolia regions (<b>c<sub>1</sub></b>–<b>c<sub>3</sub></b>, <b>d<sub>1</sub></b>–<b>d<sub>3</sub></b>).</p>
Full article ">
Back to TopTop