[go: up one dir, main page]

 
 
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

Search Results (261)

Search Parameters:
Keywords = Himawari

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
30 pages, 11305 KiB  
Article
A Case Study on the Integration of Remote Sensing for Predicting Complicated Forest Fire Spread
by Pingbo Liu and Gui Zhang
Remote Sens. 2024, 16(21), 3969; https://doi.org/10.3390/rs16213969 - 25 Oct 2024
Viewed by 733
Abstract
Forest fires can occur suddenly and have significant environmental, economic, and social consequences. The timely and accurate evaluation and prediction of their progression, particularly the spread speed in difficult-to-access areas, are essential for emergency management departments to proactively implement prevention strategies and extinguish [...] Read more.
Forest fires can occur suddenly and have significant environmental, economic, and social consequences. The timely and accurate evaluation and prediction of their progression, particularly the spread speed in difficult-to-access areas, are essential for emergency management departments to proactively implement prevention strategies and extinguish fires using scientific methods. This paper provides an analysis of models for predicting forest fire spread in China and globally. Incorporating remote sensing (RS) technology and forest fire science as the theoretical foundation, and utilizing the Wang Zhengfei forest fire spread model (1983), which is noted for its broad adaptability in China as the technical framework, this study constructs a forest fire spread model based on remote sensing interpretation. The model improves the existing model by adding elevation an factor and optimizes the method for acquiring certain parameters. By considering regional landforms (ridge lines, valley lines, and slopes) and vegetation coverage, this paper establishes three-dimensional visual interpretation markers for identifying hotspots; the orientation of the hotspots can be identified to simulate the spread of the fire uphill, downhill, in the direction of the wind, left-level slope, and right-level slope. Then, the data of Sentinel-2 and DEM were used to invert the fuel humidity and slope of pixels in the fire line areas. The statistical inversion data from pixels, which replaced fixed-point values in traditional models, were utilized for predicting forest fire spread speed. In this paper, the model was applied to the case of a forest fire in Mianning County, Sichuan Province, China, and verified using high-time-resolution Himawari-8 data, Gaofen-4 data, and historical data. The results demonstrate that the direction and maximum speed of fire spread for the fire lines in Baifen Mountai, Jiaguer Villageand, Muchanggou, Xujiabaozi, and Zhaizigou are uphill, 16.5 m/min; wind direction, 17.32 m/min; wind direction, 1.59 m/min; and wind direction, 5.67 m/min. The differences are mainly due to the locations of the fire lines, moisture content of combustibles, and maximum slopes being different. Across the entire fire line area, the average rate of increase in the area of open flames within one hour was 3.257 hm2/10 min (square hectares per 10 min), closely matching the average increase rate (3.297 hm2/10 min) monitored by the Himawari-8 satellite in 10 min intervals. In contrast, conventional fixed-point fire spread models predicted an average rate of increase of 3.5637 hm2/10 min, which shows a larger discrepancy compared to the Himawari-8 satellite monitoring results. Moreover, when compared to the fire spot monitoring results from the Gaofen-4 satellite taken 54 min after the initial location of the fire line, the predictions from the RS-enabled fire spread model, which integrates remote sensing interpretations, closely matched the actual observed fire boundaries. Although the predictions from the RS-enabled fire spread model and the traditional model both align with historical data in terms of the overall fire development trends, the RS-enabled model exhibits higher reliability and can provide more accurate information for forest fire emergency departments, enabling effective pre-emptive measures and scientific firefighting strategies. Full article
Show Figures

Figure 1

Figure 1
<p>Forest fire area in Mianning County, Sichuan, in April 2021.</p>
Full article ">Figure 2
<p>Technical route of forest fire spread simulation based on RS interpretation.</p>
Full article ">Figure 3
<p>Relationship of temperature and fire spread speed.</p>
Full article ">Figure 4
<p>The angle <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="sans-serif">α</mi> </mrow> <mrow> <mi mathvariant="normal">N</mi> <mi mathvariant="normal">I</mi> <mi mathvariant="normal">R</mi> </mrow> </msub> </mrow> </semantics></math> formed at the NIR by the reflectance at bands red, NIR, and SWIR1. An additive offset was applied to make spectral values equal at the NIR band (adapted from Shruti Khanna et al. [<a href="#B12-remotesensing-16-03969" class="html-bibr">12</a>]).</p>
Full article ">Figure 5
<p>True-color composite image of the Mianning 4.24 forest fire.</p>
Full article ">Figure 6
<p>False-color composite image of the Mianning 4.24 forest fire ( numbers 1–6 are the range of fire lines).</p>
Full article ">Figure 7
<p>Three-dimensional visualization of the fire scene ( numbers 1–6 are the range of fire lines).</p>
Full article ">Figure 8
<p>Normalized vegetation index after error correction.</p>
Full article ">Figure 9
<p>The estimated moisture content in combustibles.</p>
Full article ">Figure 10
<p>Result of slope extraction from the DEM.</p>
Full article ">Figure 11
<p>Prediction of fire line spread distance one hour later.</p>
Full article ">Figure 12
<p>Fire area land cover-type map.</p>
Full article ">Figure 13
<p>Map of the fire line and Jingkun Expressway locations.</p>
Full article ">Figure 14
<p>Dynamic map of the increase in the visible fire area at the Mianning 4.24 forest fire site monitored by the Himawai-8 satellite.</p>
Full article ">Figure 15
<p>Himawari-8 remote sensing imagery brightness temperature anomaly monitoring results from 24 April 2021 T03:40 to 24 April 2021 T05:00.</p>
Full article ">Figure 15 Cont.
<p>Himawari-8 remote sensing imagery brightness temperature anomaly monitoring results from 24 April 2021 T03:40 to 24 April 2021 T05:00.</p>
Full article ">Figure 16
<p>Comparison of forest fire spread model simulation results and remote sensing monitoring. (<b>a</b>) Forest fire spread simulation results; (<b>b</b>) Himawari-8 remote sensing imagery brightness temperature anomaly monitoring at T03:40; (<b>c</b>) Himawari-8 remote sensing imagery brightness temperature anomaly monitoring at T05:00. (The green line is the Jingkun Expressway; the yellow line is the comparison region).</p>
Full article ">Figure 17
<p>Gaofen-4 mid-infrared brightness temperature map.</p>
Full article ">Figure 18
<p>Gaofen-4 mid-infrared brightness temperature map (the yellow indicates the boundary of the abnormal brightness temperature).</p>
Full article ">Figure 19
<p>Overlaid results of the abnormal brightness temperature boundary and the predicted forest fire spread area (the green line indicates the abnormal brightness temperature boundary; the yellow line indicates the forest fire spread area).</p>
Full article ">
18 pages, 6507 KiB  
Article
Estimation of PM2.5 Using Multi-Angle Polarized TOA Reflectance Data from the GF-5B Satellite
by Ruijie Zhang, Hui Chen, Ruizhi Chen, Chunyan Zhou, Qing Li, Huizhen Xie and Zhongting Wang
Remote Sens. 2024, 16(21), 3944; https://doi.org/10.3390/rs16213944 - 23 Oct 2024
Viewed by 559
Abstract
The use of satellite data to estimate PM2.5 is an appropriate approach for long-term, substantial monitoring and assessment. To estimate PM2.5, the majority of the algorithms now in use utilize the top-of-atmosphere (TOA) reflectance or aerosol optical depth (AOD) derived [...] Read more.
The use of satellite data to estimate PM2.5 is an appropriate approach for long-term, substantial monitoring and assessment. To estimate PM2.5, the majority of the algorithms now in use utilize the top-of-atmosphere (TOA) reflectance or aerosol optical depth (AOD) derived from scalar satellite data. However, there is relatively little research on the retrieval of PM2.5 using multi-angle polarized data. With its directional polarimetric camera (DPC), the Chinese new-generation satellite Gaofen 5B (henceforth referred to as GF-5B) offers a unique opportunity to close this gap in multi-angle polarized observation data. In this research, we utilized TOA data from the DPC payload and applied the gradient boosting machine method to simulate the impact of the observation angle, wavelength, and polarization information on the accuracy of PM2.5 retrieval. We identified the optimal conditions for the effective estimation of PM2.5. The quantitative results indicated that, under these optimal conditions, the PM2.5 concentrations retrieved by GF-5B showed a strong correlation with the ground-based data, achieving an R2 of 0.9272 and an RMSE of 7.38 µg·m−3. By contrast, Himawari-8’s retrieval accuracy under similar data conditions consisted of an R2 of 0.9099 and RMSE of 7.42 µg·m−3, indicating that GF-5B offers higher accuracy. Furthermore, the retrieval results in this study demonstrated an R2 of 0.81 when compared to the CHAP dataset, confirming the feasibility and effectiveness of the use of GF-5B for PM2.5 retrieval and providing support for PM2.5 estimation through multi-angle polarized data. Full article
Show Figures

Figure 1

Figure 1
<p>Geographical location and administrative division of the study area.</p>
Full article ">Figure 2
<p>Process flow diagram for estimation of PM<sub>2.5</sub> concentrations proposed in this study.</p>
Full article ">Figure 3
<p>Impact of model parameters on inversion ((<b>Left</b>): n_estimators, (<b>Right</b>): n_depth).</p>
Full article ">Figure 4
<p>Relative importance ranking of the variables in the model.</p>
Full article ">Figure 5
<p>Impact of different angles and bands on inversion accuracy (only TOA remote sensing data).</p>
Full article ">Figure 6
<p>Inversion accuracy of PM<sub>2.5</sub> under auxiliary data.</p>
Full article ">Figure 7
<p>Impact of different angles and bands on inversion accuracy (TOA data combined with meteorological and auxiliary data).</p>
Full article ">Figure 8
<p>Comparison of inversion accuracy between Scheme 1 (<b>Left</b>) and Scheme 2 (<b>Right</b>).</p>
Full article ">Figure 9
<p>Monthly results of PM<sub>2.5</sub> estimation by GF-5B ((<b>a</b>–<b>l</b>) representing January to December).</p>
Full article ">Figure 10
<p>Comparison of annual average values between GF-5B (<b>a</b>) results and the CHAP dataset (<b>b</b>).</p>
Full article ">Figure 11
<p>Comparison of estimation results between GF-5B (<b>left</b>) and Himawari-8 (<b>right</b>).</p>
Full article ">Figure 12
<p>The maps of the R<sup>2</sup>, RMSE, and MAE for various stations in the Beijing–Tianjin–Hebei region derived from GF-5B (<b>a</b>–<b>c</b>) and Himawari-8 (<b>d</b>–<b>f</b>).</p>
Full article ">Figure 13
<p>Scatter plot comparing GF-5B results with CHAP data..</p>
Full article ">
11 pages, 3674 KiB  
Communication
Characterizing the Supercooled Cloud over the TP Eastern Slope in 2016 via Himawari-8 Products
by Qiuyu Wu, Jinghua Chen and Yan Yin
Remote Sens. 2024, 16(19), 3643; https://doi.org/10.3390/rs16193643 - 29 Sep 2024
Viewed by 492
Abstract
Supercooled liquid water (SLW) refers to droplets in clouds that remain unfrozen at temperatures below 0 °C. SLW is an important intermediate hydrometeor in the processes of snowfall and rainfall that can modulate the radiation budget. This study investigates the distribution of supercooled [...] Read more.
Supercooled liquid water (SLW) refers to droplets in clouds that remain unfrozen at temperatures below 0 °C. SLW is an important intermediate hydrometeor in the processes of snowfall and rainfall that can modulate the radiation budget. This study investigates the distribution of supercooled cloud water over mainland China using the East Asia–Pacific cloud macro- and microphysical properties dataset (2016), derived from Himawari-8 observations. The results show that the highest frequency of SLW in liquid-phase stratus clouds occur at the eastern slope of the Tibetan Plateau, the western side of the Sichuan Basin. Additional SLW is mostly found in liquid-phase clouds over the Sichuan Basin and its adjacent areas in southern China. In the region with the highest frequency of SLW, the mechanical forcing of the Tibetan Plateau causes the convergence of low-level airflow within the basin, which also carries moisture that is forced to ascend stably, creating a favorable condition for the formation of supercooled clouds. As the airflow continues to ascend, it encounters the mid-to-upper-level westerlies and temperature inversion. At the mid-to-upper level, the westerlies exhibit stronger wind speeds, directing flow towards the basin. Concurrently, the temperature inversion stabilizes the atmospheric stratification, limiting the further ascent of airflow. This inversion can also restrain convection and upward motion within the clouds, allowing for SLW to exist and persist for an extended period. Full article
Show Figures

Figure 1

Figure 1
<p>Typical area geographical location (the area within the red box is a typical region).</p>
Full article ">Figure 2
<p>The distribution of supercooled cloud water in 2016: (<b>a</b>) the frequency of supercooled cloud water occurrences; (<b>b</b>) the probability of supercooled cloud water occurring within clouds.</p>
Full article ">Figure 3
<p>Distribution of supercooled water in different cloud types and cloud phases. (<b>a</b>) Distribution of supercooled water in liquid cloud. (<b>b</b>) Distribution of supercooled water in mixed-phase clouds. (<b>c</b>) Distribution of supercooled water in altostratus clouds. (<b>d</b>) Distribution of supercooled water in nimbostratus clouds.</p>
Full article ">Figure 4
<p>Cloud top pressure and cloud top temperature statistics. (<b>a</b>) Statistics of cloud top pressure in January, February, November, and December. (<b>b</b>) Statistics of cloud top temperature in January, February, November, and December (the values on the Y-axis are normalized and then aggregated).</p>
Full article ">Figure 5
<p>Average vertical divergence.</p>
Full article ">Figure 6
<p>The average water vapor flux and the average winds at 850 hPa in (<b>a</b>) January, (<b>b</b>) February, (<b>c</b>) November, and (<b>d</b>) December.</p>
Full article ">Figure 7
<p>The average IVT (integrated water vapor transport) and the average winds in (<b>a</b>) January, (<b>b</b>) February, (<b>c</b>) November, and (<b>d</b>) December (the values &gt; 8 or &lt;−8, treated as the same value for interpolation, and positive values show that water vapor is transported in the direction of the wind, while negative values indicate the opposite).</p>
Full article ">
24 pages, 6198 KiB  
Article
The China Coastal Front from Himawari-8 AHI SST Data—Part 2: South China Sea
by Igor M. Belkin, Shang-Shang Lou, Yi-Tao Zang and Wen-Bin Yin
Remote Sens. 2024, 16(18), 3415; https://doi.org/10.3390/rs16183415 - 14 Sep 2024
Viewed by 441
Abstract
High-resolution (2 km) high-frequency (hourly) SST data from 2015 to 2021 provided by the Advanced Himawari Imager (AHI) onboard the Japanese Himawari-8 geostationary satellite were used to study spatial and temporal variability of the China Coastal Front (CCF) in the South China Sea. [...] Read more.
High-resolution (2 km) high-frequency (hourly) SST data from 2015 to 2021 provided by the Advanced Himawari Imager (AHI) onboard the Japanese Himawari-8 geostationary satellite were used to study spatial and temporal variability of the China Coastal Front (CCF) in the South China Sea. The SST data were processed with the Belkin and O’Reilly (2009) algorithm to generate monthly maps of the CCF’s intensity (defined as SST gradient magnitude GM) and frontal frequency (FF). The horizontal structure of the CCF was investigated from cross-frontal distributions of SST along 11 fixed lines that allowed us to determine inshore and offshore boundaries of the CCF and calculate the CCF’s strength (defined as total cross-frontal step of SST). Combined with the results of Part 1 of this study, where the CCF was documented in the East China Sea, the new results reported in this paper allowed the CCF to be traced from the Yangtze Bank to Hainan Island. The CCF is continuous in winter, when its intensity peaks at 0.15 °C/km (based on monthly data). In summer, when the Guangdong Coastal Current reverses and flows eastward, the CCF’s intensity is reduced to 0.05 °C/km or less, especially off western Guangdong, where the CCF vanishes almost completely. Owing to its breadth (50–100 km, up to 200 km in the Taiwan Strait), the CCF is a very strong front, especially in winter, when the total SST step across the CCF peaks at 9 °C in the Taiwan Strait. The CCF’s strength decreases westward to 6 °C off eastern Guangdong, 5 °C off western Guangdong, and 2 °C off Hainan Island, all in mid-winter. Full article
(This article belongs to the Section Ocean Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Long-term (2015–2021) mean monthly SST (°C) in the northern South China Sea.</p>
Full article ">Figure 2
<p>Histograms of long-term (2015–2021) mean monthly SST gradient magnitude GM.</p>
Full article ">Figure 3
<p>Long-term (2015–2021) mean monthly gradient magnitude GM of SST. Color scales of GM are adjusted monthly using the respective monthly histograms of GM (<a href="#remotesensing-16-03415-f002" class="html-fig">Figure 2</a>).</p>
Full article ">Figure 4
<p>Long-term (2015–2021) mean monthly frontal frequency FF at GM ≥ 0.1 °C/km.</p>
Full article ">Figure 5
<p>Bathymetry of the northern South China Sea and locations of 11 fixed lines.</p>
Full article ">Figure 6
<p>Long-term (2015–2021) mean monthly distributions of SST along 11 meridional and zonal lines across the northern South China Sea in January–June. The SST curve numbers in the plot legends correspond to the fixed line numbers in <a href="#remotesensing-16-03415-f005" class="html-fig">Figure 5</a>.</p>
Full article ">Figure 7
<p>Long-term (2015–2021) mean monthly distributions of SST along 11 meridional and zonal lines across the northern South China Sea in July–December. The SST curve numbers in the plot legends correspond to the fixed line numbers in <a href="#remotesensing-16-03415-f005" class="html-fig">Figure 5</a>.</p>
Full article ">
27 pages, 4362 KiB  
Article
Himawari-8 Sea Surface Temperature Products from the Australian Bureau of Meteorology
by Pallavi Govekar, Christopher Griffin, Owen Embury, Jonathan Mittaz, Helen Mary Beggs and Christopher J. Merchant
Remote Sens. 2024, 16(18), 3381; https://doi.org/10.3390/rs16183381 - 11 Sep 2024
Viewed by 845
Abstract
As a contribution to the Integrated Marine Observing System (IMOS), the Bureau of Meteorology introduces new reprocessed Himawari-8 satellite-derived Sea Surface Temperature (SST) products. The Radiative Transfer Model and a Bayesian cloud clearing method is used to retrieve SSTs every 10 min from [...] Read more.
As a contribution to the Integrated Marine Observing System (IMOS), the Bureau of Meteorology introduces new reprocessed Himawari-8 satellite-derived Sea Surface Temperature (SST) products. The Radiative Transfer Model and a Bayesian cloud clearing method is used to retrieve SSTs every 10 min from the geostationary satellite Himawari-8. An empirical Sensor Specific Error Statistics (SSES) model, introduced herein, is applied to calculate bias and standard deviation for the retrieved SSTs. The SST retrieval and compositing method, along with validation results, are discussed. The monthly statistics for comparisons of Himawari-8 Level 2 Product (L2P) skin SST against in situ SST quality monitoring (iQuam) in situ SST datasets, adjusted for thermal stratification, showed a mean bias of −0.2/−0.1 K and a standard deviation of 0.4–0.7 K for daytime/night-time after bias correction, where satellite zenith angles were less than 60° and the quality level was greater than 2. For ease of use, these native resolution SST data have been composited using a method introduced herein that retains retrieved measurements, to hourly, 4-hourly and daily SST products, and projected onto the rectangular IMOS 0.02 degree grid. On average, 4-hourly products cover ≈10% more of the IMOS domain, while one-night composites cover ≈25% more of the IMOS domain than a typical 1 h composite. All available Himawari-8 data have been reprocessed for the September 2015–December 2022 period. The 10 min temporal resolution of the newly developed Himawari-8 SST data enables a daily composite with enhanced spatial coverage, effectively filling in SST gaps caused by transient clouds occlusion. Anticipated benefits of the new Himawari-8 products include enhanced data quality for applications like IMOS OceanCurrent and investigations into marine thermal stress, marine heatwaves, and ocean upwelling in near-coastal regions. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) SST, (<b>b</b>) probability of a pixel being clear, (<b>c</b>) sensitivity, and (<b>d</b>) assigned quality levels for one random L2P on 15 December 2020, 20:00:00 UTC, for all quality levels.</p>
Full article ">Figure 2
<p>(<b>a</b>) <math display="inline"><semantics> <mrow> <mi mathvariant="monospace">sses</mi> <mo>_</mo> <mi mathvariant="monospace">bias</mi> </mrow> </semantics></math>, (<b>b</b>) <math display="inline"><semantics> <mrow> <mi mathvariant="monospace">sses</mi> <mo>_</mo> <mi mathvariant="monospace">standard</mi> <mo>_</mo> <mi mathvariant="monospace">deviation</mi> </mrow> </semantics></math>, and (<b>c</b>) <math display="inline"><semantics> <mrow> <mi mathvariant="monospace">sses</mi> <mo>_</mo> <mi mathvariant="monospace">count</mi> </mrow> </semantics></math> for one random L2P on 15th December 2020, 20:00:00 UTC.</p>
Full article ">Figure 3
<p>Full disk spatial coverage Himawari-8 L2P validation against drifting buoys and tropical moorings, September 2015–December 2022, showing the impact of bias correction on day (<b>top</b>) and night (<b>bottom</b>) retrievals, <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Δ</mi> <mi>SST</mi> </mrow> </semantics></math>. Left—hand panels show variables before bias correction, and right—hand panels after bias correction have been applied to the SST values. The grey region indicates pixels with no data.</p>
Full article ">Figure 4
<p>Same as <a href="#remotesensing-16-03381-f003" class="html-fig">Figure 3</a>, for <math display="inline"><semantics> <mrow> <mi>σ</mi> <mi mathvariant="sans-serif">Δ</mi> <mi>SST</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>Same as <a href="#remotesensing-16-03381-f003" class="html-fig">Figure 3</a>, for <math display="inline"><semantics> <mrow> <mi>σ</mi> <mi>zSST</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>Annual and diurnal performance of full disk Himawari-8 L2P validation against drifting buoys and tropical moorings for September 2015–December 2022 for (<b>a</b>) northern and (<b>b</b>) southern part of the disk. The colour indicates mean bias, when bias-corrected SST is compared with drifting buoys and tropical moorings.</p>
Full article ">Figure 7
<p>Same as <a href="#remotesensing-16-03381-f006" class="html-fig">Figure 6</a>. Here, the colour indicates standard deviation, when bias-corrected SST compared with drifting buoys and tropical moorings for (<b>a</b>) northern and (<b>b</b>) southern part of the disk.</p>
Full article ">Figure 8
<p>Same as <a href="#remotesensing-16-03381-f006" class="html-fig">Figure 6</a>. Here, the colour indicates <math display="inline"><semantics> <msub> <mi>z</mi> <mrow> <mi>s</mi> <mi>s</mi> <mi>e</mi> <mi>s</mi> </mrow> </msub> </semantics></math>, when bias-corrected SST is compared with drifting buoys and tropical moorings for (<b>a</b>) northern and (<b>b</b>) southern part of the disk.</p>
Full article ">Figure 9
<p>Full disk Himawari-8 L2P skin SST validation against drifting buoys and tropical moorings, 30-day running statistics, September 2015–December 2022, (<b>a</b>) uncorrected mean, (<b>b</b>) bias-corrected mean, (<b>c</b>) uncorrected standard deviation, and (<b>d</b>) bias-corrected standard deviation, for QL ≥ 3.</p>
Full article ">Figure 9 Cont.
<p>Full disk Himawari-8 L2P skin SST validation against drifting buoys and tropical moorings, 30-day running statistics, September 2015–December 2022, (<b>a</b>) uncorrected mean, (<b>b</b>) bias-corrected mean, (<b>c</b>) uncorrected standard deviation, and (<b>d</b>) bias-corrected standard deviation, for QL ≥ 3.</p>
Full article ">Figure 10
<p>Composite SSTs for (<b>a</b>) 1 h, (<b>b</b>) 4 h and (<b>c</b>) 1 night on the IMOS domain for 15 December 2020.</p>
Full article ">Figure 11
<p>Monthly statistics for validation of 1-hour L3C skin SST against drifting buoys and tropical moorings for September 2015–December 2022 on the IMOS domain, uncorrected (<b>a</b>) mean and (<b>c</b>) standard deviation, bias-corrected (<b>b</b>) mean and (<b>d</b>) standard deviation, for QL ≥ 3. Daytime validations are shown in blue, and night-time validations in orange.</p>
Full article ">Figure 11 Cont.
<p>Monthly statistics for validation of 1-hour L3C skin SST against drifting buoys and tropical moorings for September 2015–December 2022 on the IMOS domain, uncorrected (<b>a</b>) mean and (<b>c</b>) standard deviation, bias-corrected (<b>b</b>) mean and (<b>d</b>) standard deviation, for QL ≥ 3. Daytime validations are shown in blue, and night-time validations in orange.</p>
Full article ">Figure 12
<p>Same as <a href="#remotesensing-16-03381-f011" class="html-fig">Figure 11</a>, for L3C-4hour SSTs.</p>
Full article ">Figure 12 Cont.
<p>Same as <a href="#remotesensing-16-03381-f011" class="html-fig">Figure 11</a>, for L3C-4hour SSTs.</p>
Full article ">Figure 13
<p>Monthly statistics for 1-day Night L3C skin SST validation against drifting buoys and tropical moorings for September 2015–December 2022 for the IMOS domain (<b>a</b>) mean, (<b>b</b>) standard deviation, for QL ≥ 3. The brown line denotes uncorrected data, whereas the cyan line corresponds to bias−corrected data.</p>
Full article ">Figure 14
<p>SST data coverage from (<b>a</b>) MultiSensor and (<b>b</b>) GeoPolar MultiSensor L3S SST product on the IMOS domain for 15th December 2020.</p>
Full article ">
19 pages, 6287 KiB  
Article
Research on Multiscale Atmospheric Chaos Based on Infrared Remote-Sensing and Reanalysis Data
by Zhong Wang, Shengli Sun, Wenjun Xu, Rui Chen, Yijun Ma and Gaorui Liu
Remote Sens. 2024, 16(18), 3376; https://doi.org/10.3390/rs16183376 - 11 Sep 2024
Viewed by 638
Abstract
The atmosphere is a complex nonlinear system, with the information of its temperature, water vapor, pressure, and cloud being crucial aspects of remote-sensing data analysis. There exist intricate interactions among these internal components, such as convection, radiation, and humidity exchange. Atmospheric phenomena span [...] Read more.
The atmosphere is a complex nonlinear system, with the information of its temperature, water vapor, pressure, and cloud being crucial aspects of remote-sensing data analysis. There exist intricate interactions among these internal components, such as convection, radiation, and humidity exchange. Atmospheric phenomena span multiple spatial and temporal scales, from small-scale thunderstorms to large-scale events like El Niño. The dynamic interactions across different scales, along with external disturbances to the atmospheric system, such as variations in solar radiation and Earth surface conditions, contribute to the chaotic nature of the atmosphere, making long-term predictions challenging. Grasping the intrinsic chaotic dynamics is essential for advancing atmospheric analysis, which holds profound implications for enhancing meteorological forecasts, mitigating disaster risks, and safeguarding ecological systems. To validate the chaotic nature of the atmosphere, this paper reviewed the definitions and main features of chaotic systems, elucidated the method of phase space reconstruction centered on Takens’ theorem, and categorized the qualitative and quantitative methods for determining the chaotic nature of time series data. Among quantitative methods, the Wolf method is used to calculate the Largest Lyapunov Exponents, while the G–P method is used to calculate the correlation dimensions. A new method named Improved Saturated Correlation Dimension method was proposed to address the subjectivity and noise sensitivity inherent in the traditional G–P method. Subsequently, the Largest Lyapunov Exponents and saturated correlation dimensions were utilized to conduct a quantitative analysis of FY-4A and Himawari-8 remote-sensing infrared observation data, and ERA5 reanalysis data. For both short-term remote-sensing data and long-term reanalysis data, the results showed that more than 99.91% of the regional points have corresponding sequences with positive Largest Lyapunov exponents and all the regional points have correlation dimensions that tended to saturate at values greater than 1 with increasing embedding dimensions, thereby proving that the atmospheric system exhibits chaotic properties on both short and long temporal scales, with extreme sensitivity to initial conditions. This conclusion provided a theoretical foundation for the short-term prediction of atmospheric infrared radiation field variables and the detection of weak, time-sensitive signals in complex atmospheric environments. Full article
(This article belongs to the Topic Atmospheric Chemistry, Aging, and Dynamics)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Analysis process of chaotic nature.</p>
Full article ">Figure 2
<p>(<b>a</b>) 60 × 0 LLE diagram of FY-4A AGRI China Area 4KM L1 data for Channel 1; (<b>b</b>) 60 × 60 LLE diagram of FY-4A AGRI China Area 4KM L1 data for Channel 4; (<b>c</b>) 60 × 60 LLE diagram of FY-4A AGRI China Area 4KM L1 data for Channel 10; (<b>d</b>) 60 × 60 LLE diagram of FY-4A AGRI China Area 4KM L1 data for Channel 13.</p>
Full article ">Figure 3
<p>(<b>a</b>) 40 × 40 LLE diagram of Himawari-8 AHI Full Disk 2KM L1 data for Channel 2; (<b>b</b>) 40 × 40 LLE diagram of Himawari-8 AHI Full Disk 2KM L1 data for Channel 12.</p>
Full article ">Figure 4
<p>(<b>a</b>) 41 × 41 LLE diagram of ERA5 hourly data for Z500; (<b>b</b>) 41 × 41 LLE diagram of ERA5 hourly data for T850.</p>
Full article ">Figure 5
<p>(<b>a</b>) Log–log plots of the correlation integral <math display="inline"><semantics> <mrow> <mi>C</mi> <mfenced> <mi>r</mi> </mfenced> </mrow> </semantics></math> versus the distance <math display="inline"><semantics> <mi>r</mi> </semantics></math> for different embedding dimensions <span class="html-italic">m</span> for channel 1 of FY-4A AGRI China area 4KM L1 data; (<b>b</b>) Log–log plots of the correlation integral <math display="inline"><semantics> <mrow> <mi>C</mi> <mfenced> <mi>r</mi> </mfenced> </mrow> </semantics></math> versus the distance <math display="inline"><semantics> <mi>r</mi> </semantics></math> for different embedding dimensions <span class="html-italic">m</span> for channel 4 of FY-4A AGRI China area 4KM L1 data; (<b>c</b>) Log–log plots of the correlation integral <math display="inline"><semantics> <mrow> <mi>C</mi> <mfenced> <mi>r</mi> </mfenced> </mrow> </semantics></math> versus the distance <math display="inline"><semantics> <mi>r</mi> </semantics></math> for different embedding dimensions <span class="html-italic">m</span> for channel 10 of FY-4A AGRI China area 4KM L1 data; (<b>d</b>) Log–log plots of the correlation integral <math display="inline"><semantics> <mrow> <mi>C</mi> <mfenced> <mi>r</mi> </mfenced> </mrow> </semantics></math> versus the distance <math display="inline"><semantics> <mi>r</mi> </semantics></math> for different embedding dimensions <span class="html-italic">m</span> for channel 13 of FY-4A AGRI China area 4KM L1 data of FY-4A AGRI China area 4KM L1 data.</p>
Full article ">Figure 6
<p>Correlation Dimension vs. Embedding Dimension Curve for channel 1, channel 4, channel 10 and channel 13 of FY-4A AGRI China area 4KM L1 data.</p>
Full article ">Figure 7
<p>(<b>a</b>) Log–log plots of the correlation integral <math display="inline"><semantics> <mrow> <mi>C</mi> <mfenced> <mi>r</mi> </mfenced> </mrow> </semantics></math> versus the distance <math display="inline"><semantics> <mi>r</mi> </semantics></math> for different embedding dimensions <span class="html-italic">m</span> for channel 2 of Himawari-8 AHI full-disk 2KM L1 data; (<b>b</b>) Log–log plots of the correlation integral <math display="inline"><semantics> <mrow> <mi>C</mi> <mfenced> <mi>r</mi> </mfenced> </mrow> </semantics></math> versus the distance <math display="inline"><semantics> <mi>r</mi> </semantics></math> for different embedding dimensions <span class="html-italic">m</span> for channel 12 of Himawari-8 AHI full-disk 2KM L1 data.</p>
Full article ">Figure 8
<p>Correlation Dimension vs. Embedding Dimension Curve for channel 2 and channel 12 of Himawari-8 AHI Full disk 2KM L1 data at location (4050, 5250).</p>
Full article ">Figure 9
<p>(<b>a</b>) Log–log plots of the correlation integral <math display="inline"><semantics> <mrow> <mi>C</mi> <mfenced> <mi>r</mi> </mfenced> </mrow> </semantics></math> versus the distance <math display="inline"><semantics> <mi>r</mi> </semantics></math> for different embedding dimensions <span class="html-italic">m</span> for Z500 hourly data of ERA5 at 20.00°S, 175.00°W; (<b>b</b>) Log–log plots of the correlation integral <math display="inline"><semantics> <mrow> <mi>C</mi> <mfenced> <mi>r</mi> </mfenced> </mrow> </semantics></math> versus the distance <math display="inline"><semantics> <mi>r</mi> </semantics></math> for different embedding dimensions <span class="html-italic">m</span> for T850 hourly data of ERA5 at 20.00°S, 175.00°W.</p>
Full article ">Figure 10
<p>(<b>a</b>) Correlation Dimension vs. Embedding Dimension Curve for Z500 hourly data of ERA5 at 20.00°S, 175.00°W; (<b>b</b>) Correlation Dimension vs. Embedding Dimension Curve for T850 hourly data of ERA5 at 20.00°S, 175.00°W.</p>
Full article ">
11 pages, 3215 KiB  
Article
Heat and Drought Have Exacerbated the Midday Depression Observed in a Subtropical Fir Forest by a Geostationary Satellite
by Qianqian Xie, Kexin Chen, Tong Li, Jia Liu, Yuqiu Wang and Xiaolu Zhou
Forests 2024, 15(9), 1572; https://doi.org/10.3390/f15091572 - 7 Sep 2024
Viewed by 615
Abstract
Recently, increasing heat and drought events have threatened the resilience of Chinese fir forests. Trees primarily respond to these threats by downregulating photosynthesis including through stomatal limitation that causes a drop in productivity at noon (known as the midday depression). However, the effects [...] Read more.
Recently, increasing heat and drought events have threatened the resilience of Chinese fir forests. Trees primarily respond to these threats by downregulating photosynthesis including through stomatal limitation that causes a drop in productivity at noon (known as the midday depression). However, the effects of these events on midday and afternoon GPP inhibition are rarely analyzed on a fine timescale. This may result in negligence of critical responses. Here, we investigated the impact of climatic events on the midday depression of photosynthesis at a subtropical fir forest in Huitong from 2016 to 2022 using data from the Himawari 8 meteorological satellite and flux tower. Our results indicated that the highest number of midday depression occurred in 2022 (126 times) with the highest average temperature (29.1 °C). A higher incidence of midday depression occurred in summer and autumn, with 48 and 34 occurrences, respectively. Compound drought, heat, and drought events induced increases in midday depression at 74.3%, 66.0%, and 47.5%. Thus, trees are more likely to adopt midday depression as an adaptive strategy during compound drought and heat events. This study can inform forest management and lead to improvements in Earth system models. Full article
(This article belongs to the Section Forest Inventory, Modeling and Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>The site location of the flux tower in the Huitong fir forest site.</p>
Full article ">Figure 2
<p>Comparison of GPP estimates with measured values of the flux tower between 2016 and 2022. The number of available samples was 14633. Each sample represents a half-hourly observation. (<b>a</b>) Prediction using observed climate data from the flux tower. The black dashed line is the data fit line. The red solid line is the 45° tangent. The color intensity of the right bar represents the density of data points. Unit for vertical coordinates is g C m<sup>−2</sup> h<sup>−1</sup>. (<b>b</b>) Residual distribution, RMSE, and MAE corresponding to (<b>a</b>). (<b>c</b>) Prediction using the satellite data. (<b>d</b>) Residual distribution, RMSE, and MAE corresponding to (<b>c</b>).</p>
Full article ">Figure 3
<p>Observed variables at the Huitong fir forest site from 2016 to 2022. (<b>a</b>) GPP value, (<b>b</b>) air temperature, (<b>c</b>) precipitation. (<b>d</b>) Average daily temperature and (<b>e</b>) precipitation during climatic events. In (<b>d</b>,<b>e</b>), 1 represents average daily temperature or precipitation during heat events, 2 indicates drought events, 3 indicates compound events, and Ave. indicates overall climatic event averages.</p>
Full article ">Figure 4
<p>The number of midday depression occurrences in different years (<b>a</b>), different months (<b>b</b>), and different seasons (<b>c</b>) at the Huitong fir forest site from 2016 to 2022.</p>
Full article ">Figure 5
<p>Regression between the number of midday depression occurrences and the number of climatic event occurrences from 2016 to 2022. (<b>a</b>) Regression between the number of midday depression occurrences and the number of heat event occurrences. (<b>b</b>) Regression between the number of midday depression occurrences and the number of drought event occurrences. (<b>c</b>) Regression between the number of midday depression occurrences and the number of compositeevents. (<b>d</b>) Mean annual number of climatic event occurrences from 2016 to 2022. The C event, D event, and HT event represent composite events (compound drought and heat), drought events, and heat events, respectively.</p>
Full article ">Figure 6
<p>Probability of midday depression being concurrently triggered by climatic events between 2016 and 2022. HT event, D event, and C event represent heat events, drought events, and composite events (compound drought and heat), respectively.</p>
Full article ">
23 pages, 24947 KiB  
Article
Quality Assessment and Application Scenario Analysis of AGRI Land Aerosol Product from the Geostationary Satellite Fengyun-4B in China
by Nan Wang, Bingqian Li, Zhili Jin and Wei Wang
Sensors 2024, 24(16), 5309; https://doi.org/10.3390/s24165309 - 16 Aug 2024
Viewed by 567
Abstract
The Advanced Geostationary Radiation Imager (AGRI) sensor on board the geostationary satellite Fengyun-4B (FY-4B) is capable of capturing particles in different phases in the atmospheric environment and acquiring aerosol observation data with high spatial and temporal resolution. To understand the quality of the [...] Read more.
The Advanced Geostationary Radiation Imager (AGRI) sensor on board the geostationary satellite Fengyun-4B (FY-4B) is capable of capturing particles in different phases in the atmospheric environment and acquiring aerosol observation data with high spatial and temporal resolution. To understand the quality of the Land Aerosol (LDA) product of AGRI and its application prospects, we conducted a comprehensive evaluation of the AGRI LDA AOD. Using the 550 nm AGRI LDA AOD (550 nm) of nearly 1 year (1 October 2022 to 30 September 2023) to compare with the Aerosol Robotic Network (AERONET), MODIS MAIAC, and Himawari-9/AHI AODs. Results show the erratic algorithmic performance of AGRI LDA AOD, the correlation coefficient (R), mean error (Bias), root mean square error (RMSE), and the percentage of data with errors falling within the expected error envelope of ±(0.05+0.15×AODAERONET) (within EE15) of the LDA AOD dataset are 0.55, 0.328, 0.533, and 34%, respectively. The LDA AOD appears to be overestimated easily in the southern and western regions of China and performs poorly in the offshore areas, with an R of 0.43, a Bias of 0.334, a larger RMSE of 0.597, and a global climate observing system fraction (GCOSF) percentage of 15% compared to the inland areas (R = 0.60, Bias = 0.163, RMSE = 0.509, GCOSF = 17%). Future improvements should focus on surface reflectance calculation, water vapor attenuation, and more suitable aerosol model selection to improve the algorithm’s accuracy. Full article
(This article belongs to the Special Issue Recent Trends in Air Quality Sensing)
Show Figures

Figure 1

Figure 1
<p>The terrain, land type, and distribution of AERONET sites in the study area: (<b>a</b>) topographic maps; (<b>b</b>) land type; and AERONET site distribution.</p>
Full article ">Figure 2
<p>Overall distribution and validation scatter plots: (<b>a</b>) LDA AOD, (<b>b</b>) AHI AOD, (<b>c</b>) MAIAC AOD, and (<b>d</b>–<b>f</b>) Scatter plots of AGRI LDA AOD, AHI AOD, and MAIAC AOD compared with AERONET measurement data, respectively. The orange dashed line represents the linear fit between the ground-based measurements and the AOD dataset need be validated. The black dashed line indicates the envelope of the expected error (EE15).</p>
Full article ">Figure 3
<p>Validation results of LDA AOD and AERONET observations in inland and offshore areas. (<b>a</b>) Scatter plot for inland region validation; (<b>b</b>) Box plot for inland region error distribution; (<b>c</b>) Distribution of average percentage error for inland region; (<b>d</b>) Scatter plot for offshore region validation; (<b>e</b>) Box plot for offshore region error distribution; and (<b>f</b>) Distribution of average percentage error for offshore region.</p>
Full article ">Figure 4
<p>Time series of 10 min (solid line) and monthly average AOD (solid line with marked points) for FY4B and AERONET, October 2022–September 2023. (<b>a</b>) Inland areas and (<b>b</b>) offshore areas.</p>
Full article ">Figure 5
<p>Scatter plots of FY4B LDA AOD 550 versus AERONET observations in the inland region for different time periods. (<b>a</b>) 8:00~9:00, (<b>b</b>) 9:00~10:00, (<b>c</b>) 10:00~11:00, (<b>d</b>) 11:00~12:00, (<b>e</b>) 12:00~13:00, (<b>f</b>) 13:00~14:00, (<b>g</b>) 14:00~15:00, (<b>h</b>) 15:00~16:00, and (<b>i</b>) 16:00~17:00.</p>
Full article ">Figure 6
<p>Scatter plots of FY4B LDA AOD 550 versus AERONET observations in the coastal region for different time periods. (<b>a</b>) 8:00~9:00, (<b>b</b>) 9:00~10:00, (<b>c</b>) 10:00~11:00, (<b>d</b>) 11:00~12:00, (<b>e</b>) 12:00~13:00, (<b>f</b>) 13:00~14:00, (<b>g</b>) 14:00~15:00, (<b>h</b>) 15:00~16:00, and (<b>i</b>) 16:00~17:00.</p>
Full article ">Figure 7
<p>Angular dependence of the AOD 550 bias (<math display="inline"><semantics> <mrow> <mi>B</mi> <mi>i</mi> <mi>a</mi> <mi>s</mi> <mo>=</mo> <msub> <mi>τ</mi> <mrow> <mi>L</mi> <mi>D</mi> <mi>A</mi> </mrow> </msub> <mo>−</mo> <msub> <mi>τ</mi> <mrow> <mi>aeronet</mi> </mrow> </msub> </mrow> </semantics></math>) on the number of match points (blue bars) between FY4B and AERONET: (<b>a</b>) zenith angle of satellite observations in inland and (<b>b</b>) offshore regions; (<b>c</b>) AERONET solar zenith angle of observations in inland and (<b>d</b>) offshore stations as well as AERONET; and (<b>e</b>) scattering angle of inland and (<b>f</b>) offshore observations.</p>
Full article ">Figure 8
<p>Scatter plots of FY4B LDA AOD 550 and AERONET observations in different seasons. (<b>a</b>–<b>d</b>) Inland areas; (<b>e</b>–<b>h</b>) offshore areas.</p>
Full article ">Figure 9
<p>Distribution of FY4B LDA AOD and AERONET observations in inland and offshore regions under different AE conditions. (<b>a</b>–<b>d</b>) inland areas and (<b>e</b>–<b>h</b>) offshore areas.</p>
Full article ">Figure 10
<p>Scatterplot of FY4B LDA AOD 550 with AERONET observations and different aerosol types. (<b>a</b>) Dust type; (<b>b</b>) Mixture type; (<b>c</b>) Non-absorbing type; (<b>d</b>) Slightly-absorbing type; (<b>e</b>) Moderately-absorbing type; and (<b>f</b>) Highly-absorbing type.</p>
Full article ">Figure 11
<p>Daily AOD time series of FY4B and AERONET measurements at three AERONET stations for different major aerosol types; the temporal variation AOD (<b>a</b>) and aerosol type distribution (<b>b</b>) at XiangHe Site, the temporal variation AOD (<b>c</b>) and aerosol type distribution (<b>d</b>) at Dhaka University Site, and the temporal variation AOD (<b>e</b>) and aerosol type distribution (<b>f</b>) at Kaohsiung Site.</p>
Full article ">Figure 12
<p>Scatterplot of LDA AOD 550 nm versus AERONET AOD measurements under different land use types obtained by the MCD12Q1 product. (<b>a</b>) broadleaf evergreen forest; (<b>b</b>) broadleaf deciduous forest; (<b>c</b>) mixed forest; (<b>d</b>) woody savanna; (<b>e</b>) savanna; (<b>f</b>) grassland; (<b>g</b>) agricultural land; (<b>h</b>) urban and built-up areas; (<b>i</b>) agricultural land/natural vegetation.</p>
Full article ">Figure 13
<p>Time series (solid line with markers) of daily mean AOD 550 nm bias (shown by the left <span class="html-italic">y</span>-axis, <math display="inline"><semantics> <mrow> <mi>B</mi> <mi>i</mi> <mi>a</mi> <mi>s</mi> <mo>=</mo> <msub> <mi>τ</mi> <mrow> <mi>L</mi> <mi>D</mi> <mi>A</mi> </mrow> </msub> <mo>−</mo> <msub> <mi>τ</mi> <mrow> <mi>A</mi> <mi>E</mi> <mi>R</mi> <mi>O</mi> <mi>N</mi> <mi>E</mi> <mi>T</mi> </mrow> </msub> </mrow> </semantics></math>), and MODIS NDVI (right <span class="html-italic">y</span>-axis) for the inland region. (<b>a</b>) April, (<b>b</b>) May.</p>
Full article ">Figure 14
<p>Retrieval error bias (<math display="inline"><semantics> <mrow> <mi>B</mi> <mi>i</mi> <mi>a</mi> <mi>s</mi> <mo>=</mo> <msub> <mi>τ</mi> <mrow> <mi>L</mi> <mi>D</mi> <mi>A</mi> </mrow> </msub> <mo>−</mo> <msub> <mi>τ</mi> <mrow> <mi>A</mi> <mi>E</mi> <mi>R</mi> <mi>O</mi> <mi>N</mi> <mi>E</mi> <mi>T</mi> </mrow> </msub> </mrow> </semantics></math>) dependence analysis of LDA AOD, AHI AOD, and MODIS AOD. (<b>a</b>) MODIS NDVI, (<b>b</b>) AE index, (<b>c</b>) solar zenith angle, (<b>d</b>) AERONET water vapor content, (<b>e</b>) ozone concentration, and (<b>f</b>) NO<sub>2</sub> concentration.</p>
Full article ">Figure 15
<p>Validation results of FY4B LDA 550 nm AOD and MODIS 550 nm AOD with AERONET observations in the inland region (<b>a</b>,<b>c</b>) and offshore region (<b>b</b>,<b>d</b>), respectively.</p>
Full article ">Figure 16
<p>Validated scatterplots of LDA AOD, MODIS AOD, and AERONET data for the inland region after removing angular effects (removing SOZ_AERONET &lt; 45° and SAZ_FY4B &lt; 45°). (<b>a</b>) LDA AOD, (<b>b</b>) MODIS AOD.</p>
Full article ">
22 pages, 21022 KiB  
Article
Forest Fire Detection Based on Spatial Characteristics of Surface Temperature
by Houzhi Yao, Zhigao Yang, Gui Zhang and Feng Liu
Remote Sens. 2024, 16(16), 2945; https://doi.org/10.3390/rs16162945 - 12 Aug 2024
Viewed by 1539
Abstract
Amidst the escalating threat of global warming, which manifests in more frequent forest fires, the prompt and accurate detection of forest fires has ascended to paramount importance. The current surveillance algorithms employed for forest fire monitoring—including, but not limited to, fixed threshold algorithms, [...] Read more.
Amidst the escalating threat of global warming, which manifests in more frequent forest fires, the prompt and accurate detection of forest fires has ascended to paramount importance. The current surveillance algorithms employed for forest fire monitoring—including, but not limited to, fixed threshold algorithms, multi-channel threshold algorithms, and contextual algorithms—rely primarily upon the degree of deviation between the pixel temperature and the background temperature to discern pyric events. Notwithstanding, these algorithms typically fail to account for the spatial heterogeneity of the background temperature, precipitating the consequential oversight of low-temperature fire point pixels, thus impeding the expedited detection of fires in their initial stages. For the amelioration of this deficiency, the present study introduces a spatial feature-based (STF) method for forest fire detection, leveraging Himawari-8/9 imagery as the main data source, complemented by the Shuttle Radar Topography Mission (SRTM) DEM data inputs. Our proposed modality reconstructs the surface temperature information via selecting the optimally designated machine learning model, subsequently identifying the fire point through utilizing the difference between the reconstructed surface temperatures and empirical observations, in tandem with the spatial contextual algorithm. The results confirm that the random forest model demonstrates superior efficacy in the reconstruction of the surface temperature. Benchmarking the STF method against both the fire point datasets disseminated by the China Forest and Grassland Fire Prevention and Suppression Network (CFGFPN) and the Wild Land Fire (WLF) fire point product validation datasets from Himawari-8/9 yielded a zero rate of omission errors and a comprehensive evaluative index, predominantly surpassing 0.74. These findings show that the STF method proposed herein significantly augments the identification of lower-temperature fire point pixels, thereby amplifying the sensitivity of forest surveillance. Full article
Show Figures

Figure 1

Figure 1
<p>Overview map of the study area.</p>
Full article ">Figure 2
<p>Vegetation area and DEM in Hunan Province.</p>
Full article ">Figure 3
<p>Histogram of the frequency distribution of the surface temperatures in vegetation areas in Hunan Province on different dates.</p>
Full article ">Figure 4
<p>Flowchart of fire point detection algorithm.</p>
Full article ">Figure 5
<p>Feature correlation heatmap at different moments during the daytime.</p>
Full article ">Figure 6
<p>Scatter density plot of validation data for RF at different moments of the day.</p>
Full article ">Figure 7
<p>Scatter density plot of reconstructed LST versus original LST.</p>
Full article ">Figure 8
<p>LST of original vs. reconstructed vegetation area during daytime.</p>
Full article ">Figure 9
<p>LST of original vs. reconstructed area at nighttime.</p>
Full article ">Figure 10
<p>The result of fire point identification at 15:30 on 18 October 2022.</p>
Full article ">Figure 11
<p>The result of fire point identification at 10:30 on 19 October 2022.</p>
Full article ">Figure 12
<p>The result of fire point identification at 15:20 on 23 October 2022.</p>
Full article ">Figure 13
<p>The results of fire detection.</p>
Full article ">Figure 14
<p>Identification results of fire point image elements in Xintian County, Hunan Province, at four moments on 18 and 19 October 2022. (<b>a</b>) Mid-infrared 7th band of Himawari-9 image and its bright temperature. (<b>b</b>) Identification results of the algorithm of this study. (<b>c</b>) Results of WLF fire point product.</p>
Full article ">
18 pages, 5958 KiB  
Article
Oceanic Precipitation Nowcasting Using a UNet-Based Residual and Attention Network and Real-Time Himawari-8 Images
by Xianpu Ji, Xiaojiang Song, Anboyu Guo, Kai Liu, Haijin Cao and Tao Feng
Remote Sens. 2024, 16(16), 2871; https://doi.org/10.3390/rs16162871 - 6 Aug 2024
Viewed by 1333
Abstract
Qualitative precipitation forecasting plays a vital role in marine operational services. However, predicting heavy precipitation over the open ocean presents a significant challenge due to the limited availability of ground-based radar observations far from coastal regions. Recent advancements in deep learning models offer [...] Read more.
Qualitative precipitation forecasting plays a vital role in marine operational services. However, predicting heavy precipitation over the open ocean presents a significant challenge due to the limited availability of ground-based radar observations far from coastal regions. Recent advancements in deep learning models offer potential for oceanic precipitation nowcasting using satellite images. This study implemented an enhanced UNet model with an attention mechanism and a residual architecture (RA-UNet) to predict the precipitation rate within a 90 min time frame. A comparative analysis with the standard UNet and UNet with an attention algorithm revealed that the RA-UNet method exhibited superior accuracy metrics, such as the critical ratio index and probability of detection, with fewer false alarms. Two typical cases demonstrated that RA-UNet had a better ability to forecast monsoon precipitation as well as intense precipitation in a tropical cyclone. These findings indicate the greater potential of the RA-UNet approach for nowcasting heavy precipitation over the ocean using satellite imagery. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>An example of the input data valid at 1800 UTC on 3 June 2017. (<b>a</b>–<b>g</b>) Brightness temperature (K) in each spectral band from the Himawari-8 satellite images; (<b>h</b>) precipitation rate from the GPM IMERG products (mm h<sup>−1</sup>).</p>
Full article ">Figure 2
<p>Flowchart of data preprocessing for the Himawari-8 and GPM datasets.</p>
Full article ">Figure 3
<p>Probability density function of precipitation samples from the GPM data (unit: mm h<sup>−1</sup>) before (<b>a</b>) and after (<b>b</b>) performing the <math display="inline"><semantics> <mrow> <mi>l</mi> <mi>o</mi> <mi>g</mi> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>x</mi> <mo>)</mo> </mrow> </semantics></math> transformation.</p>
Full article ">Figure 4
<p>A diagram of RA-UNet.</p>
Full article ">Figure 5
<p>Different variants of convolutional and recurrent convolutional units: (<b>a</b>) forward convolutional units; (<b>b</b>) residual convolutional units.</p>
Full article ">Figure 6
<p>The critical success index, probability of detection, and false alarm ratio for three different intensity thresholds (0.1, 1, and 5 mm h<sup>−1</sup>). The metrics are shown as a function of the lead time. The blue, orange, and red lines, respectively, represent the forecast performance of the UNet model, Att-UNet model, and RA-UNet model with the lead time.</p>
Full article ">Figure 7
<p>The critical success index (CSI), probability of detection (POD), and false alarm ratio (FAR) for three different intensity thresholds (0.1, 1, and 5 mm h<sup>−1</sup>) averaged over 30–90 min forecasts for each month of 2018. The metrics are shown as a function of the month. The blue, orange, and red lines, respectively, represent the forecast performance of the UNet model, Att-UNet model, and RA-UNet model in each month.</p>
Full article ">Figure 8
<p>An example of a forecast during a summer monsoon. The brightness temperature (band 16) of inputs (<b>a</b>–<b>c</b>), the ground truth (<b>d</b>–<b>f</b>), and nowcasts from the UNet (<b>g</b>–<b>i</b>), Att-UNet (<b>j</b>–<b>l</b>), and RA-UNet (<b>m</b>–<b>o</b>) models corresponding to 30 min (first column), 60 min (second column), and 90 min forecasts.</p>
Full article ">Figure 9
<p>A prediction example that occurred during a typhoon. Brightness temperature (band 16) of inputs ((<b>a</b>–<b>c</b>), 29 October 2018 at 11:30, 12:30, and 13:30, respectively), ground truth ((<b>d</b>–<b>f</b>), 29 October 2018 at 14:30, 15:00, and 15:30, respectively), and nowcasts from the UNet (<b>g</b>–<b>i</b>), Att-UNet (<b>j</b>–<b>l</b>), and RA-UNet (<b>m</b>–<b>o</b>) models corresponding to 30 min (first column), 60 min (second column), and 90 min forecasts.</p>
Full article ">
13 pages, 3316 KiB  
Article
Evaluation of the Forecasting Performance of Supercooled Clouds for the Weather Modification Model of the Cloud and Precipitation Explicit Forecasting System
by Jia Wang, Qin Mei, Haixia Mei, Jun Guo and Tongchang Liu
Atmosphere 2024, 15(8), 928; https://doi.org/10.3390/atmos15080928 - 3 Aug 2024
Viewed by 596
Abstract
Through the application of cloud top temperature data and the extraction of supercooled cloud information in cloud-type data from the next-generation Himawari-8 geostationary satellite with high spatial–temporal resolution, a quantitative evaluation of the forecasting performance of the weather modification model named the Cloud [...] Read more.
Through the application of cloud top temperature data and the extraction of supercooled cloud information in cloud-type data from the next-generation Himawari-8 geostationary satellite with high spatial–temporal resolution, a quantitative evaluation of the forecasting performance of the weather modification model named the Cloud and Precipitation Explicit Forecasting System (CPEFS) was conducted. The evaluation, based on selected forecast cases from 8 days in September and October 2018 initialized at 00 and 12 UTC every day, focused especially on the forecasting performance in supercooled clouds (vertical integrated supercooled liquid water, VISL > 0), including the comprehensive spatial distribution of cloud top temperature (CTT) and 3 h precipitation over 0.1 mm (R3 > 0.1). The results indicated that the forecasting performance for VISL > 0 was relatively good, with the Threat Score (TS) ranging from 0.46 to 0.67. The forecasts initialized at 12 UTC slightly outperformed the forecasts initialized at 00 UTC. Additionally, the corresponding spatial Anomaly Correlation Coefficient (ACC) of CTT between forecasts and observations was 0.23, and the TS for R3 > 0.1 reached as high as 0.87. For a mix of cold and warm cloud systems, there was a correlation between the forecasting performance of VISL > 0 and CTT. The trends in the TS for VISL > 0 and the ACC of CTT aligned with the forecast lead-time. Full article
(This article belongs to the Special Issue Numerical Weather Prediction Models and Ensemble Prediction Systems)
Show Figures

Figure 1

Figure 1
<p>The 48 h average spatial Anomaly Correlation Coefficient of cloud top temperature (CTT_ACC), Threat Score for 3 h precipitation exceeding 0.1 mm, and vertical integrated supercooled liquid water exceeding 0 (R3 &gt; 0.1_TS, VISL &gt; 0_TS) between observations and model forecast results (All represents the average for all 16 cases, while S00 and S12 represent the averages for cases initialized at 00 and 12 UTC, respectively).</p>
Full article ">Figure 2
<p>Performance diagram summarizing the average Threat Score (TS, labeled solid contours), False Alarm Rate (FAR, <span class="html-italic">X</span>-axis), and Missed Detection Rate (MDR, <span class="html-italic">Y</span>-axis) of supercooled clouds (VISL &gt; 0) in 48 h forecasts for different cases (different shapes, first two digits represent the month and last two digits represent the date of the cases) initialized at 00 (S00, dark blue) and 12 UTC (S12, light blue) over 8 days in the model (ALL is the average of S00 and S12, black).</p>
Full article ">Figure 3
<p>The scatter plot between the TS of supercooled clouds and spatial Anomaly Correlation Coefficient (ACC) of cloud top temperature (CTT) for 8 cases (0904, 0919, 0921, 0922, 0923, 0934, 1021, and 1021) initialized at 00 and 12 UTC (S00, S12, and All as shown in <a href="#atmosphere-15-00928-f002" class="html-fig">Figure 2</a>).</p>
Full article ">Figure 4
<p>Comparison of the 27th-hour forecasts of CTT (K, <b>right</b>) for cases 1021 (<b>top</b>) and 0924 (<b>bottom</b>) initialized at 12 UTC with observations (<b>left</b>).</p>
Full article ">Figure 5
<p>Comparison of the 30th-hour forecasts (<b>right</b>) of CTT (K, <b>top</b>) and supercooled water (pink areas, <b>bottom</b>) for case 0921 initialized at 12 UTC with observations (<b>left</b>).</p>
Full article ">Figure 6
<p>Sequence of variations in the spatial ACC of CTT and TS of supercooled clouds with forecast lead-time for cases 0924 (<b>top</b>), 0921 (<b>middle</b>), and 1021 (<b>bottom</b>) initialized at 12 UTC.</p>
Full article ">Figure 7
<p>Comparison of the 30th-hour forecasts (<b>right</b>) of CTT (K, <b>top</b>) and supercooled water (pink areas, <b>bottom</b>) for case 0904 initialized at 12 UTC with observations (<b>left</b>).</p>
Full article ">
22 pages, 3924 KiB  
Article
Diurnal Variation in Surface Incident Solar Radiation Retrieved by CERES and Himawari-8
by Lu Lu, Ying Li, Lingjun Liang and Qian Ma
Remote Sens. 2024, 16(14), 2670; https://doi.org/10.3390/rs16142670 - 22 Jul 2024
Viewed by 599
Abstract
The diurnal variation of surface incident solar radiation (Rs) has a significant impact on the Earth’s climate. Satellite-retrieved Rs datasets display good spatial and temporal continuity compared with ground-based observations and, more importantly, have higher accuracy than reanalysis datasets. Facilitated by these advantages, [...] Read more.
The diurnal variation of surface incident solar radiation (Rs) has a significant impact on the Earth’s climate. Satellite-retrieved Rs datasets display good spatial and temporal continuity compared with ground-based observations and, more importantly, have higher accuracy than reanalysis datasets. Facilitated by these advantages, many scholars have evaluated satellite-retrieved Rs, especially based on monthly and annual data. However, there is a lack of evaluation on an hourly scale, which has a profound impact on sea–air interactions, climate change, agriculture, and prognostic models. This study evaluates Himawari-8 and Clouds and the Earth’s Radiant Energy System Synoptic (CERES)-retrieved hourly Rs data covering 60°S–60°N and 80°E–160°W based on ground-based observations from the Baseline Surface Radiation Network (BSRN). Hourly Rs were first standardized to remove the diurnal and seasonal cycles. Furthermore, the sensitivities of satellite-retrieved Rs products to clouds, aerosols, and land cover types were explored. It was found that Himawari-8-retrieved Rs was better than CERES-retrieved Rs at 8:00–16:00 and worse at 7:00 and 17:00. Both satellites performed better at continental sites than at island/coastal sites. The diurnal variations of statistical parameters of Himawari-8 satellite-retrieved Rs were stronger than those of CERES. Relatively larger MABs in the case of stratus and stratocumulus were exhibited for both hourly products. Smaller MAB values were found for CERES covered by deep convection and cumulus clouds and for Himawari-8 covered by deep convection and nimbostratus clouds. Larger MAB values at evergreen broadleaf forest sites and smaller MAB values at open shrubland sites were found for both products. In addition, Rs retrieved by Himawari-8 was more sensitive to AOD at 10:00–16:00, while that retrieved by CERES was more sensitive to COD at 9:00–15:00. The CERES product showed larger sensitivity to COD (at 9:00–15:00) and AOD (at 7:00–10:00) than Himawari-8. This work helps data producers know how to improve their future products and helps data users be aware of the uncertainties that exist in hourly satellite-retrieved Rs data. Full article
Show Figures

Figure 1

Figure 1
<p>Geographical distribution of observational sites used for the evaluation of satellite-retrieved Rs data. The red and magenta circles indicate the site location on the continent and island/coast.</p>
Full article ">Figure 2
<p>The sample size of the observational data at different times. The red line indicates 40% of the total sample size.</p>
Full article ">Figure 3
<p>Scatter plots of the annual average of hourly satellite-retrieved and observed Rs from 2015 to 2021 at 7:00–17:00 (<b>a</b>–<b>k</b>) and all hours (<b>l</b>).</p>
Full article ">Figure 4
<p>Diurnal variations of statistical parameters between hourly satellite-retrieved Rs and observed Rs for different types of sites. (<b>a</b>) Bias %; (<b>b</b>) mean absolute bias (MAB) %; (<b>c</b>) root mean square error (RMSE) %; and (<b>d</b>) correlation coefficient (R).</p>
Full article ">Figure 5
<p>Taylor diagram describing the standard deviation and correlation coefficient between the hourly satellite-retrieved Rs and observed Rs at 15 selected stations. The circles and crosses denote Himawari-8-retrieved Rs and CERES-retrieved Rs. “REF” can be regarded as the point of perfection, where the value closer to the point indicates a better evaluation.</p>
Full article ">Figure 6
<p>MAB between satellites and BSRN hourly Rs under different cloud types from 7:00 to 17:00. (<b>a</b>) CERES and (<b>b</b>) Himawari-8.</p>
Full article ">Figure 7
<p>MAB between satellites and BSRN hourly Rs under different cloud optical depth (COD) categories from 7:00 to 17:00.</p>
Full article ">Figure 8
<p>MAB between satellites and BSRN hourly Rs under different aerosol optical depth (AOD) categories from 7:00 to 17:00.</p>
Full article ">Figure 9
<p>Ranges (largest value minus smallest value) and variation (the standard deviation of the values) in MAB (shown in <a href="#remotesensing-16-02670-f007" class="html-fig">Figure 7</a> and <a href="#remotesensing-16-02670-f008" class="html-fig">Figure 8</a>) from satellite-retrieved hourly Rs at all ground-based sites for each hour under different COD and AOD conditions.</p>
Full article ">Figure 10
<p>Diurnal variations of MAB between hourly satellite-retrieved Rs and observed Rs at nine sites covered by different land cover types for 2015–2021. Solid lines for CERES and dashed lines for Himawari-8.</p>
Full article ">Figure 11
<p>Diurnal variations of MAB between hourly CERES-retrieved Rs and observed Rs at 39 sites covered by different land cover types for 2000–2021.</p>
Full article ">
25 pages, 15972 KiB  
Article
CACM-Net: Daytime Cloud Mask for AGRI Onboard the FY-4A Satellite
by Jingyuan Yang, Zhongfeng Qiu, Dongzhi Zhao, Biao Song, Jiayu Liu, Yu Wang, Kuo Liao and Kailin Li
Remote Sens. 2024, 16(14), 2660; https://doi.org/10.3390/rs16142660 - 20 Jul 2024
Viewed by 602
Abstract
Accurate cloud detection is a crucial initial stage in optical satellite remote sensing. In this study, a daytime cloud mask model is proposed for the Advanced Geostationary Radiation Imager (AGRI) onboard the Fengyun 4A (FY-4A) satellite based on a deep learning approach. The [...] Read more.
Accurate cloud detection is a crucial initial stage in optical satellite remote sensing. In this study, a daytime cloud mask model is proposed for the Advanced Geostationary Radiation Imager (AGRI) onboard the Fengyun 4A (FY-4A) satellite based on a deep learning approach. The model, named “Convolutional and Attention-based Cloud Mask Net (CACM-Net)”, was trained using the 2021 dataset with CALIPSO data as the truth value. Two CACM-Net models were trained based on a satellite zenith angle (SZA) < 70° and >70°, respectively. The study evaluated the National Satellite Meteorological Center (NSMC) cloud mask product and compared it with the method established in this paper. The results indicate that CACM-Net outperforms the NSMC cloud mask product overall. Specifically, in the SZA < 70° subset, CACM-Net enhances accuracy, precision, and F1 score by 4.8%, 7.3%, and 3.6%, respectively, while reducing the false alarm rate (FAR) by approximately 7.3%. In the SZA > 70° section, improvements of 12.2%, 19.5%, and 8% in accuracy, precision, and F1 score, respectively, were observed, with a 19.5% reduction in FAR compared to NSMC. An independent validation dataset for January–June 2023 further validates the performance of CACM-Net. The results show improvements of 3.5%, 2.2%, and 2.8% in accuracy, precision, and F1 scores for SZA < 70° and 7.8%, 11.3%, and 4.8% for SZA > 70°, respectively, along with reductions in FAR. Cross-comparison with other satellite cloud mask products reveals high levels of agreement, with 88.6% and 86.3% matching results with the MODIS and Himawari-9 products, respectively. These results confirm the reliability of the CACM-Net cloud mask model, which can produce stable and high-quality FY-4A AGRI cloud mask results. Full article
Show Figures

Figure 1

Figure 1
<p>Coverage of all CALIPSO and AGRI matched points during daytime throughout 2021 and January–June 2023.</p>
Full article ">Figure 2
<p>Schematic of the daytime SZA &gt; 70° and SZA &lt; 70° portions of FY-4A AGRI, with the green line indicating SZA = 70° and the red line indicating SZA = 70°.</p>
Full article ">Figure 3
<p>Conceptual diagram of the structure of CACM-Net, which consists mainly of a training step and a prediction step, with the sizes of the input and output vectors shown at the bottom of the picture, representing the dimensional sizes of the channels, rows, and columns.</p>
Full article ">Figure 4
<p>CBAM block; the module has two consecutive submodules, namely the channel attention module and the spatial attention module. ⊗ denotes element-wise multiplication. The sizes of the input and output vectors are shown below the image, representing the dimensional sizes of the channels, rows, and columns.</p>
Full article ">Figure 5
<p>Confusion matrix schematic.</p>
Full article ">Figure 6
<p>Overall accuracy trend of CACM-Net cloud mask results and NSMC cloud mask product on the 2021 dataset. (<b>a</b>) CACM-Net cloud mask result. (<b>b</b>) NSMC cloud mask product.</p>
Full article ">Figure 7
<p>Box plots of per-batch accuracy for the last epoch after convergence for all models.</p>
Full article ">Figure 8
<p>CACM-Net training set data with rising and falling nodes showing cloud probability distributions. Vertical dashed lines indicate the probability thresholds that distinguish clear from probably clear (0.09), probably clear from probably cloudy (0.56), and probably cloudy from cloudy (0.87).</p>
Full article ">Figure 9
<p>CACM-Net (SZA &lt; 70°), CACM-Net (SZA &gt; 70°), and CACM-Net (Full) cloud mask evaluation metrics including accuracy, POD, precision, F1 score, and FAR referenced to the 2021 test dataset.</p>
Full article ">Figure 10
<p>Schematic comparison of the results of CACM-Net and NSMC on 21 April 2021, 05:00 UTC. (<b>a</b>) Reflectance of 0.65 μm; (<b>b</b>) difference in cloud mask between CACM-Net and NSMC, where red pixels are mostly judged as cloudy by NSMC, blue pixels are mostly judged as cloudy by CACM-Net, and white pixels represent agreement between the models; (<b>c</b>) cloud mask results for CACM-Net; (<b>d</b>) Cloud mask results for NSMC.</p>
Full article ">Figure 11
<p>CACM-Net cloud mask result, as well as NSMC cloud mask product evaluation metrics, including accuracy, POD, precision, and F1 score, (<b>a</b>) for the SZA &lt; 70° portion and (<b>b</b>) for the SZA &gt; 70° portion. (<b>c</b>) Comparison of the FAR metrics for the CACM-Net cloud mask result, as well as the NSMC cloud mask product, for the SZA &lt; 70° and SZA &gt; 70° portions using the 2023 test set data as a reference.</p>
Full article ">Figure 12
<p>Overall accuracy trends for CACM-Net cloud mask results and NSMC cloud mask products on the 2023 independently validated dataset.</p>
Full article ">Figure 13
<p>A diagram showing the formation and dissipation of Typhoon Nanmadol from 13 September 2022 to 20 September 2022, with the red dot representing the center of the typhoon.</p>
Full article ">Figure 14
<p>Schematic comparison of the results of CACM-Net and MODIS for 16 January 2023, 01:00 UTC (<b>a</b>–<b>d</b>); (<b>a</b>) reflectance of 0.65 μm; (<b>b</b>) difference in cloud mask between CACM-Net and MODIS; (<b>c</b>) cloud mask results for CACM-Net; (<b>d</b>) cloud mask results for MODIS. Schematic comparison of the results of CACM-Net and Himawari 9 for 1 January 2023, 01:00 UTC (<b>e</b>–<b>h</b>), (<b>e</b>) reflectance of 0.65 μm; (<b>f</b>) difference in cloud mask between CACM-Net and Himawari-9; (<b>g</b>) cloud mask results for CACM-Net; (<b>h</b>) cloud mask results for Himawari-9.</p>
Full article ">Figure 14 Cont.
<p>Schematic comparison of the results of CACM-Net and MODIS for 16 January 2023, 01:00 UTC (<b>a</b>–<b>d</b>); (<b>a</b>) reflectance of 0.65 μm; (<b>b</b>) difference in cloud mask between CACM-Net and MODIS; (<b>c</b>) cloud mask results for CACM-Net; (<b>d</b>) cloud mask results for MODIS. Schematic comparison of the results of CACM-Net and Himawari 9 for 1 January 2023, 01:00 UTC (<b>e</b>–<b>h</b>), (<b>e</b>) reflectance of 0.65 μm; (<b>f</b>) difference in cloud mask between CACM-Net and Himawari-9; (<b>g</b>) cloud mask results for CACM-Net; (<b>h</b>) cloud mask results for Himawari-9.</p>
Full article ">
15 pages, 5639 KiB  
Article
Identification Method for Spring Dust Intensity Levels Based on Multiple Remote Sensing Parameters
by Qi Jiang, Linchang An, Fei Wang, Guozhou Wu, Jianwei Wen, Bin Li, Yuchen Jin and Yapeng Wei
Remote Sens. 2024, 16(14), 2606; https://doi.org/10.3390/rs16142606 - 17 Jul 2024
Viewed by 658
Abstract
The advancement of more precise remote sensing inversion technology for dust aerosols has long been a hot topic in the field of the atmospheric environment. In 2023, China experienced 18 dust-related weather events, predominantly in spring. These high-intensity and frequent dust events have [...] Read more.
The advancement of more precise remote sensing inversion technology for dust aerosols has long been a hot topic in the field of the atmospheric environment. In 2023, China experienced 18 dust-related weather events, predominantly in spring. These high-intensity and frequent dust events have attracted considerable attention. However, gridded observation data of dust intensity levels are not collected in current dust monitoring and forecasting operations. Based on the Himawari 9 geostationary satellite data, this study establishes a new method to identify spring dust events. This method integrates the brightness temperature difference method and the multiple infrared dust index, taking into account the response discrepancies of the multiple infrared dust index under various underlying surfaces. Furthermore, by obtaining dynamic background brightness temperature values eight times a day, threshold statistics are applied to analyze the correlation between the infrared difference dust index and ground-observed dust level, so as to establish a satellite-based near-surface dust intensity level identification algorithm. This algorithm aims to improve dust detection accuracy, and to provide more effective gridded observation support for dust forecasting and monitoring operations. The test results indicate that the algorithm can effectively identify the presence or absence of dust, with a misjudgment rate of less than 3%. With regard to dust intensity, the identification of blowing sand and floating dust aligns relatively well with ground-based observations, but notable uncertainties exist in determining a dust intensity of sand-storm level or above. Among these uncertainties, the differences between ground-based observations and satellite identification caused by non-grounded dust in the upper air, and the selection of dust identification thresholds, are two important error sources in the dust identification results of this study. Full article
(This article belongs to the Special Issue Application of Satellite Aerosol Remote Sensing in Air Quality)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) The underlying surface type; (<b>b</b>) the frequency of dust events and the average PM10 concentration at 30 stations from March to May 2023.</p>
Full article ">Figure 2
<p>(<b>a</b>) Correlation analysis between brightness temperature difference (BTD) and multiple infrared dust index (MIDI). The gray circles represent all samples. The colored triangles represent the statistical distribution of different levels of surface dust events. (<b>b</b>) Correlation distribution for non-primary sand source, desert, and gobi. The highlighted in the red circled and yellow circled area is the misclassification of some non-dust data. The yellow dotted line represents different thresholds for MIDI and BTD.</p>
Full article ">Figure 3
<p>(<b>a</b>) The statistical distribution of MIDIs for different dust levels on different underlying surfaces. The frequency distribution of MIDIs under different underlying surface conditions during (<b>b</b>) the entire study period and (<b>c</b>) dust weather periods. (<b>d</b>) The statistical distribution of BTDs corresponding to different dust levels. (<b>e</b>,<b>f</b>) are the same as (<b>b</b>,<b>c</b>), but for BTD. In (<b>a</b>,<b>d</b>), the horizontal dashed lines denote the average statistical values of the corresponding underlying surfaces, the dots represent the mean values of the corresponding dust levels, the vertical lines represent the 10th percentile (lower) and 90th percentile (upper), and the horizontal lines from top to bottom represent the 25th, 50th, and 75th percentiles (same for Figure 5).</p>
Full article ">Figure 4
<p>The highest background brightness temperature distributions during 07–09 UTC on dust occurrence days of 21 March 2023, 10 April 2023, 19 April 2023, and 19 May 2023.</p>
Full article ">Figure 5
<p>The distribution of the infrared difference dust index (IDDI) statistical values under different PM<sub>10</sub> ranges in (<b>a</b>) non-primary sand source areas and (<b>b</b>) desert–gobi areas, and the statistics of (<b>c</b>) PM<sub>10</sub> and (<b>d</b>) IDDI for different dust intensity levels.</p>
Full article ">Figure 6
<p>Comparison of satellite identification results of dust levels with near-surface station observations at 12:00 UTC on 21 March, 2023, 09:00 UTC on 10 April, 10:00 UTC on 19 April, and 00:00 UTC on 19 May, 2023. The black-and-white backgrounds of the figures represent the cloud images at the corresponding time. The blue dots denote the distribution of ground observation stations.</p>
Full article ">Figure 7
<p>(<b>a</b>) Lidar data from Jiuquan station from 13:00 UTC on April 18, 2023, to 13:00 UTC on 20 April, 2023; (<b>b</b>) wind profiling radar data from Jiuquan station on 19 April. The direction of the arrows represents the wind direction, while the length of the arrows indicates the wind speed. The black box represents the wind field near the altitude of 2500m; (<b>c</b>) true-color composite image from H9 satellite at 12:00 UTC on 19 April and 18:00 UTC on 19 April. The blue circles represent the dust storm areas; (<b>d</b>) CALIOP aerosol subtype image from 8:06 UTC to 8:19 UTC on 19 April, 2023. The blue box indicates the area around Jiuquan at an altitude of 2–5 km.</p>
Full article ">Figure 7 Cont.
<p>(<b>a</b>) Lidar data from Jiuquan station from 13:00 UTC on April 18, 2023, to 13:00 UTC on 20 April, 2023; (<b>b</b>) wind profiling radar data from Jiuquan station on 19 April. The direction of the arrows represents the wind direction, while the length of the arrows indicates the wind speed. The black box represents the wind field near the altitude of 2500m; (<b>c</b>) true-color composite image from H9 satellite at 12:00 UTC on 19 April and 18:00 UTC on 19 April. The blue circles represent the dust storm areas; (<b>d</b>) CALIOP aerosol subtype image from 8:06 UTC to 8:19 UTC on 19 April, 2023. The blue box indicates the area around Jiuquan at an altitude of 2–5 km.</p>
Full article ">
21 pages, 9076 KiB  
Article
Satellite Remote Sensing False Forest Fire Hotspot Excavating Based on Time-Series Features
by Haifeng Wang, Gui Zhang, Zhigao Yang, Haizhou Xu, Feng Liu and Shaofeng Xie
Remote Sens. 2024, 16(13), 2488; https://doi.org/10.3390/rs16132488 - 7 Jul 2024
Viewed by 822
Abstract
Satellite remote sensing has become an important means of forest fire monitoring because it has the advantages of wide coverage, few ground constraints and high dynamics. When utilizing satellites for forest fire hotspot monitoring, two types of ground hotspots, agricultural and other fire [...] Read more.
Satellite remote sensing has become an important means of forest fire monitoring because it has the advantages of wide coverage, few ground constraints and high dynamics. When utilizing satellites for forest fire hotspot monitoring, two types of ground hotspots, agricultural and other fire hotspots can be ruled out through ground object features. False forest fire hotspots within forested areas must be excluded for a more accurate distinction between forest fires and non-forest fires. This study utilizes spatio-temporal data along with time-series classification to excavate false forest fire hotspots exhibiting temporal characteristics within forested areas and construct a dataset of such false forest fire hotspots, thereby achieving a more realistic forest fire dataset. Taking Hunan Province as the research object, this study takes the satellite ground hotspots in the forests of Hunan Province as the suspected forest fire hotspot dataset and excludes the satellite ground hotspots in the forests such as fixed heat sources, periodic heat sources and recurring heat sources which are excavated. The validity of these methods and results was then analyzed. False forest fire hotspots, from satellite ground hotspots extracted from 2019 to 2023 Himawari-8/9 satellite images, closely resemble the official release of actual forest fires data and the accuracy rate in the actual forest fire monitoring is 95.12%. This validates that the method employed in this study can improve the accuracy of satellite-based forest fire monitoring. Full article
Show Figures

Figure 1

Figure 1
<p>Overview map of the study area.</p>
Full article ">Figure 2
<p>Flow chart of the methodology.</p>
Full article ">Figure 3
<p>Sliding Window Working Principle Diagram. We used suspected forest fire hotspot data in Excel format to mine fixed heat sources according to window settings and step size.</p>
Full article ">Figure 4
<p>Schematic diagram of the curved path.</p>
Full article ">Figure 5
<p>Map of land classification results for Xinshao County.</p>
Full article ">Figure 6
<p>Map of land classification results for Xintian County.</p>
Full article ">Figure 7
<p>Map of the results of land classification of landforms in Hunan Province using random forests. The classification results of the six land classes in Hunan Province are shown in (<b>a</b>). The forest types are extracted from (<b>a</b>) to obtain (<b>b</b>). Figure (<b>b</b>) depicts the high forest coverage in Hunan Province predominantly concentrated in the western region, particularly in Huaihua City, Xiangxi Tujia and Miao Autonomous Prefecture and Zhangjiajie City. The northeastern region of the country is primarily characterized by extensive land use for watersheds and water management infrastructure, leading to a limited forested area.</p>
Full article ">Figure 8
<p>Fixed Heat Source Excavation Results Map.</p>
Full article ">Figure 9
<p>Periodic Heat Source Excavation Results Map.</p>
Full article ">Figure 10
<p>Seasonal Statistical Map of Periodic Heat Source Hotspots.</p>
Full article ">Figure 11
<p>Recurring Heat Source Excavation Results Map.</p>
Full article ">Figure 12
<p>Forest fire monitoring flowchart. When conducting forest fire monitoring, we are able to arbitrarily select a certain moment and first exclude satellite ground hotspots outside the forest based on the feature characteristics as the suspected forest fire hotspot dataset. The forest fire can be obtained by excluding the false forest fire hotspot dataset based on the false forest fire hotspot dataset constructed in this study. Emergency management forest fire verification data and comparison of land classes in high resolution imagery were also utilized in this study to verify the accuracy of forest fire monitoring.</p>
Full article ">Figure 13
<p>Map of forest fire monitoring results. Figure (<b>a</b>) shows the distribution of suspected forest fire hotspots on 12 January 2024 in Hunan Province. Figure (<b>b</b>) shows 41 forest fires and 22 false forest fire hotspots in Hunan Province on 12 January 2024 monitored in this study. Figure (<b>c</b>) shows a close-up screenshot of some of the correctly classified forest fire hotspots monitored in this study in (<b>b</b>) with small light blue boxes located in Huaihua, Zhuzhou and Chenzhou. Figure (<b>d</b>) shows the small red boxes in (<b>b</b>), located in Xiangtan and Changde, and is a close-up screenshot of 2 of the 41 forest fires monitored in this study that are false forest fire hotspots.</p>
Full article ">
Back to TopTop