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

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (172)

Search Parameters:
Keywords = SEVIRI

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
19 pages, 3947 KiB  
Article
Modeling of Biologically Effective Daily Radiant Exposures over Europe from Space Using SEVIRI Measurements and MERRA-2 Reanalysis
by Agnieszka Czerwińska and Janusz Krzyścin
Remote Sens. 2024, 16(20), 3797; https://doi.org/10.3390/rs16203797 - 12 Oct 2024
Viewed by 419
Abstract
Ultraviolet solar radiation at the Earth’s surface significantly impacts both human health and ecosystems. A biologically effective daily radiant exposure (BEDRE) model is proposed for various biological processes with an analytical formula for its action spectrum. The following processes are considered: erythema formation, [...] Read more.
Ultraviolet solar radiation at the Earth’s surface significantly impacts both human health and ecosystems. A biologically effective daily radiant exposure (BEDRE) model is proposed for various biological processes with an analytical formula for its action spectrum. The following processes are considered: erythema formation, previtamin D3 synthesis, psoriasis clearance, and inactivation of SARS-CoV-2 virions. The BEDRE model is constructed by multiplying the synthetic BEDRE value under cloudless conditions by a cloud modification factor (CMF) parameterizing the attenuation of radiation via clouds. The CMF is an empirical function of the solar zenith angle (SZA) at midday and the daily clearness index from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) measurements on board the second-generation Meteosat satellites. Total column ozone, from MERRA-2 reanalysis, is used in calculations of clear-sky BEDRE values. The proposed model was trained and validated using data from several European ground-based spectrophotometers and biometers for the periods 2014–2023 and 2004–2013, respectively. The model provides reliable estimates of BEDRE for all biological processes considered. Under snow-free conditions and SZA < 45° at midday, bias and standard deviation of observation-model differences are approximately ±5% and 15%, respectively. The BEDRE model can be used as an initial validation tool for ground-based UV data. Full article
(This article belongs to the Section Environmental Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>The location of the UV measuring stations shown in <a href="#remotesensing-16-03797-t001" class="html-table">Table 1</a> (created with Google My Maps: Map data 2024).</p>
Full article ">Figure 2
<p>Normalized action spectra for the specific biological effects: erythema appearance (black), photosynthesis of previtamin D<sub>3</sub> in human skin (blue), psoriasis clearance (green), and inactivation of SARS-CoV-2 virions (red).</p>
Full article ">Figure 3
<p>Scatter plot of UBE model against measured daily erythemal radiant exposure at Belsk for all-sky conditions for different ranges of noon SZA: (<b>a</b>) SZA &lt; 45°; (<b>b</b>) SZA ≥ 45° and SZA &lt; 60°; (<b>c</b>) SZA ≥ 60°. The dotted line is the 1–1 agreement line. The solid curve represents smoothed values from the LOWESS filter [<a href="#B43-remotesensing-16-03797" class="html-bibr">43</a>].</p>
Full article ">Figure 4
<p>Scatter plot RE<sub>BIOL</sub>(D) from the UBE model with <math display="inline"><semantics> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi mathvariant="sans-serif">α</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> <mrow> <mi mathvariant="normal">E</mi> <mi mathvariant="normal">R</mi> <mi mathvariant="normal">Y</mi> <mi mathvariant="normal">T</mi> </mrow> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi mathvariant="sans-serif">β</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> <mrow> <mi mathvariant="normal">E</mi> <mi mathvariant="normal">R</mi> <mi mathvariant="normal">Y</mi> <mi mathvariant="normal">T</mi> </mrow> </msub> </mrow> </semantics></math> when SZA<sub>N</sub> &lt; 45° versus corresponding values from spectral measurements at Belsk for the period 2011–2023: (<b>a</b>) for VITD, (<b>b</b>) for PSOR, and (<b>c</b>) for SARS.</p>
Full article ">Figure A1
<p>Scatter plot of the modeled (UBE model with <math display="inline"><semantics> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi mathvariant="sans-serif">α</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> <mrow> <mi mathvariant="normal">E</mi> <mi mathvariant="normal">R</mi> <mi mathvariant="normal">Y</mi> <mi mathvariant="normal">T</mi> </mrow> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi mathvariant="sans-serif">β</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> <mrow> <mi mathvariant="normal">E</mi> <mi mathvariant="normal">R</mi> <mi mathvariant="normal">Y</mi> <mi mathvariant="normal">T</mi> </mrow> </msub> </mrow> </semantics></math>) versus the measured daily radiant exposure for all-sky conditions and different ranges of SZA at noon: (<b>a</b>) Reading for SZA<sub>N</sub> &lt; 45°; (<b>b</b>) Reading for SZA<sub>N</sub> ≥ 45° and SZA<sub>N</sub> &lt; 60°; (<b>c</b>) Reading for SZA ≥ 60°; (<b>d</b>) Vienna for SZA<sub>N</sub> &lt;45°; (<b>e</b>) Vienna for SZA<sub>N</sub> ≥ 45° and SZA<sub>N</sub> &lt; 60°; (<b>f</b>) Vienna for SZA ≥ 60°. The dotted line is the 1–1 agreement line. The solid curve represents smoothed values from the Lowess filter [<a href="#B43-remotesensing-16-03797" class="html-bibr">43</a>].</p>
Full article ">Figure A2
<p>Scatter plot of the modeled (UBE model with <math display="inline"><semantics> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi mathvariant="sans-serif">α</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> <mrow> <mi mathvariant="normal">E</mi> <mi mathvariant="normal">R</mi> <mi mathvariant="normal">Y</mi> <mi mathvariant="normal">T</mi> </mrow> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi mathvariant="sans-serif">β</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> <mrow> <mi mathvariant="normal">E</mi> <mi mathvariant="normal">R</mi> <mi mathvariant="normal">Y</mi> <mi mathvariant="normal">T</mi> </mrow> </msub> </mrow> </semantics></math>) versus measured daily erythemal radiant exposure for different ranges of noon SZA: (SZA &lt; 45°; SZA ≥ 45° and SZA &lt; 60°; and SZA ≥ 60°: (<b>a</b>–<b>c</b>) Diekirch (Luxembourg); (<b>d</b>–<b>f</b>) Uccle (Belgium); (<b>g</b>–<b>i</b>) Davos (Switzerland); (<b>j</b>–<b>l</b>) Chisinau (Moldavia). As these stations were not used in UBE training, all available daily data in the period 2004–2023 have been used. The dotted line is the 1–1 perfect agreement line. The solid curve represents smoothed values from the Lowess filter [<a href="#B43-remotesensing-16-03797" class="html-bibr">43</a>].</p>
Full article ">Figure A2 Cont.
<p>Scatter plot of the modeled (UBE model with <math display="inline"><semantics> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi mathvariant="sans-serif">α</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> <mrow> <mi mathvariant="normal">E</mi> <mi mathvariant="normal">R</mi> <mi mathvariant="normal">Y</mi> <mi mathvariant="normal">T</mi> </mrow> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi mathvariant="sans-serif">β</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> <mrow> <mi mathvariant="normal">E</mi> <mi mathvariant="normal">R</mi> <mi mathvariant="normal">Y</mi> <mi mathvariant="normal">T</mi> </mrow> </msub> </mrow> </semantics></math>) versus measured daily erythemal radiant exposure for different ranges of noon SZA: (SZA &lt; 45°; SZA ≥ 45° and SZA &lt; 60°; and SZA ≥ 60°: (<b>a</b>–<b>c</b>) Diekirch (Luxembourg); (<b>d</b>–<b>f</b>) Uccle (Belgium); (<b>g</b>–<b>i</b>) Davos (Switzerland); (<b>j</b>–<b>l</b>) Chisinau (Moldavia). As these stations were not used in UBE training, all available daily data in the period 2004–2023 have been used. The dotted line is the 1–1 perfect agreement line. The solid curve represents smoothed values from the Lowess filter [<a href="#B43-remotesensing-16-03797" class="html-bibr">43</a>].</p>
Full article ">Figure A3
<p>Scatter plot of the modeled (UBE model with <math display="inline"><semantics> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi mathvariant="sans-serif">α</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> <mrow> <mi mathvariant="normal">E</mi> <mi mathvariant="normal">R</mi> <mi mathvariant="normal">Y</mi> <mi mathvariant="normal">T</mi> </mrow> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi mathvariant="sans-serif">β</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> <mrow> <mi mathvariant="normal">E</mi> <mi mathvariant="normal">R</mi> <mi mathvariant="normal">Y</mi> <mi mathvariant="normal">T</mi> </mrow> </msub> </mrow> </semantics></math>) versus the measured daily radiant exposure for all-sky conditions and SZA at noon less than 45°: (<b>a</b>) VITD for Reading; (<b>b</b>) PSOR for Reading; (<b>c</b>) SARS for Reading; (<b>d</b>) VITD for Uccle; (<b>e</b>) PSOR for Uccle; (<b>f</b>) SARS for Uccle. The dotted line is the 1–1 agreement line. The solid curve represents smoothed values from the Lowess filter [<a href="#B43-remotesensing-16-03797" class="html-bibr">43</a>].</p>
Full article ">Figure A3 Cont.
<p>Scatter plot of the modeled (UBE model with <math display="inline"><semantics> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi mathvariant="sans-serif">α</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> <mrow> <mi mathvariant="normal">E</mi> <mi mathvariant="normal">R</mi> <mi mathvariant="normal">Y</mi> <mi mathvariant="normal">T</mi> </mrow> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi mathvariant="sans-serif">β</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> <mrow> <mi mathvariant="normal">E</mi> <mi mathvariant="normal">R</mi> <mi mathvariant="normal">Y</mi> <mi mathvariant="normal">T</mi> </mrow> </msub> </mrow> </semantics></math>) versus the measured daily radiant exposure for all-sky conditions and SZA at noon less than 45°: (<b>a</b>) VITD for Reading; (<b>b</b>) PSOR for Reading; (<b>c</b>) SARS for Reading; (<b>d</b>) VITD for Uccle; (<b>e</b>) PSOR for Uccle; (<b>f</b>) SARS for Uccle. The dotted line is the 1–1 agreement line. The solid curve represents smoothed values from the Lowess filter [<a href="#B43-remotesensing-16-03797" class="html-bibr">43</a>].</p>
Full article ">
23 pages, 32897 KiB  
Article
On the Suitability of Different Satellite Land Surface Temperature Products to Study Surface Urban Heat Islands
by Alexandra Hurduc, Sofia L. Ermida and Carlos C. DaCamara
Remote Sens. 2024, 16(20), 3765; https://doi.org/10.3390/rs16203765 - 10 Oct 2024
Viewed by 959
Abstract
Remote sensing satellite data have been a crucial tool in understanding urban climates. The variety of sensors with different spatiotemporal characteristics and retrieval methodologies gave rise to a multitude of approaches when analyzing the surface urban heat island effect (SUHI). Although there are [...] Read more.
Remote sensing satellite data have been a crucial tool in understanding urban climates. The variety of sensors with different spatiotemporal characteristics and retrieval methodologies gave rise to a multitude of approaches when analyzing the surface urban heat island effect (SUHI). Although there are considerable advantages that arise from these different characteristics (spatiotemporal resolution, time of observation, etc.), it also means that there is a need for understanding the ability of sensors in capturing spatial and temporal SUHI patterns. For this, several land surface temperature products are compared for the cities of Madrid and Paris, retrieved from five sensors: the Spinning Enhanced Visible and InfraRed Imager onboard Meteosat Second Generation, the Advanced Very-High-Resolution Radiometer onboard Metop, the Moderate-resolution Imaging Spectroradiometer onboard both Aqua and Terra, and the Thermal Infrared Sensor onboard Landsat 8 and 9. These products span a wide range of LST algorithms, including split-window, single-channel, and temperature–emissivity separation methods. Results show that the diurnal amplitude of SUHI may not be well represented when considering daytime and nighttime polar orbiting platforms. Also, significant differences arise in SUHI intensity and spatial and temporal variability due to the different methods implemented for LST retrieval. Full article
(This article belongs to the Section AI Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Land cover resampled for the three projections of LST products (<b>a</b>–<b>f</b>) along with the percentage of urban pixels (<b>g</b>–<b>l</b>).</p>
Full article ">Figure 2
<p>Time of observation of each sensor: (<b>a</b>) for Madrid during daytime and time of minimum SUHI (SUHI<sub>min</sub>), (<b>b</b>) for Madrid during nighttime and time of maximum SUHI (SUHI<sub>max</sub>), (<b>c</b>) for Paris during daytime and SUHI<sub>max</sub>, (<b>d</b>) for Paris during nighttime and SUHI<sub>min</sub>. Colored bins are sampled every 15 min.</p>
Full article ">Figure 3
<p>Mean DJF (December, January, February) LST for all products considered. (<b>a1</b>–<b>a9</b>) The spatial pattern during daytime, and in the case of the MLST, the most frequent hour of the LST maximum and SUHI minimum are shown; (<b>b1</b>–<b>b9</b>) histograms of urban and rural LST shown in (<b>a1</b>–<b>a9</b>); (<b>c1</b>–<b>c7</b>) as in the first line but for nighttime and for the LST minimum and SUHI maximum; (<b>d1</b>–<b>d7</b>) as in (<b>b1</b>–<b>b9</b>) but for nighttime. Please note that color bars are different amongst the different products to allow a better visualization of patterns, but value ranges of the histograms are the same.</p>
Full article ">Figure 4
<p>As <a href="#remotesensing-16-03765-f003" class="html-fig">Figure 3</a> but for JJA (June, July, and August). (<b>a1</b>–<b>a9</b>) The spatial pattern during daytime, and in the case of the MLST, the most frequent hour of the LST maximum and SUHI minimum are shown; (<b>b1</b>–<b>b9</b>) histograms of urban and rural LST shown in (<b>a1</b>–<b>a9</b>); (<b>c1</b>–<b>c7</b>) as in the first line but for nighttime and for the LST minimum and SUHI maximum; (<b>d1</b>–<b>d7</b>) as in (<b>b1</b>–<b>b9</b>) but for nighttime. Please note that color bars are different amongst the different products to allow a better visualization of patterns, but value ranges of the histograms are the same.</p>
Full article ">Figure 5
<p>As <a href="#remotesensing-16-03765-f003" class="html-fig">Figure 3</a> but for Paris. (<b>a1</b>–<b>a9</b>) The spatial pattern during daytime, and in the case of the MLST, the most frequent hour of the LST and SUHI maximum; (<b>b1</b>–<b>b9</b>) histograms of urban and rural LST shown in (<b>a1</b>–<b>a9</b>); (<b>c1</b>–<b>c7</b>) as in the first line but for nighttime and for the LST and SUHI minimum; (<b>d1</b>–<b>d7</b>) as in (<b>b1</b>–<b>b9</b>) but for nighttime. Please note that color bars are different amongst the different products to allow a better visualization of patterns, but value ranges of the histograms are the same.</p>
Full article ">Figure 6
<p>As <a href="#remotesensing-16-03765-f003" class="html-fig">Figure 3</a> but for Paris and DJF; (<b>a1</b>–<b>a9</b>) The spatial pattern during daytime, and in the case of the MLST, the most frequent hour of the LST and SUHI maximum; (<b>b1</b>–<b>b9</b>) histograms of urban and rural LST shown in (<b>a1</b>–<b>a9</b>); (<b>c1</b>–<b>c7</b>) as in the first line but for nighttime and for the LST and SUHI minimum; (<b>d1</b>–<b>d7</b>) as in (<b>b1</b>–<b>b9</b>) but for nighttime, an extension of the histogram in (<b>d6</b>) is seen in (<b>d8</b>). Please note that color bars are different amongst the different products to allow a better visualization of patterns, but value ranges of the histograms are the same.</p>
Full article ">Figure 7
<p>Diurnal cycle of SUHI for Madrid: (<b>a</b>) DJF, (<b>b</b>) MAM, (<b>c</b>) JJA, (<b>d</b>) SON.</p>
Full article ">Figure 8
<p>As <a href="#remotesensing-16-03765-f007" class="html-fig">Figure 7</a> but for Paris.</p>
Full article ">Figure 9
<p>Correlation of monthly SUHI anomalies between all products considered: (<b>a</b>) daytime, (<b>b</b>) nighttime. Blank spaces correspond to pairs of products with no significant correlation (<span class="html-italic">p</span>-value &gt; 0.05).</p>
Full article ">Figure 10
<p>As <a href="#remotesensing-16-03765-f009" class="html-fig">Figure 9</a> but for Paris. (<b>a</b>) daytime, (<b>b</b>) nighttime. Blank spaces correspond to pairs of products with no significant correlation (<span class="html-italic">p</span>-value &gt; 0.05).</p>
Full article ">
27 pages, 20066 KiB  
Article
First Release of the Optimal Cloud Analysis Climate Data Record from the EUMETSAT SEVIRI Measurements 2004–2019
by Alessio Bozzo, Marie Doutriaux-Boucher, John Jackson, Loredana Spezzi, Alessio Lattanzio and Philip D. Watts
Remote Sens. 2024, 16(16), 2989; https://doi.org/10.3390/rs16162989 - 14 Aug 2024
Viewed by 703
Abstract
Clouds are key to understanding the atmosphere and climate, and a long series of satellite observations provide invaluable information to study their properties. EUMETSAT has published Release 1 of the Optimal Cloud Analysis (OCA) Climate Data Record (CDR), which provides a homogeneous time [...] Read more.
Clouds are key to understanding the atmosphere and climate, and a long series of satellite observations provide invaluable information to study their properties. EUMETSAT has published Release 1 of the Optimal Cloud Analysis (OCA) Climate Data Record (CDR), which provides a homogeneous time series of cloud properties of up to two overlapping layers, together with uncertainties. The OCA product is derived using the 15 min Spinning Enhanced Visible and Infrared Imager (SEVIRI) measurements onboard Meteosat Second Generation (MSG) in geostationary orbit and covers the period from 19 January 2004 until 31 August 2019. This paper presents the validation of the OCA cloud-top pressure (CTP) against independent lidar-based estimates and the quality assessment of the cloud optical thickness (COT) and cloud particle effective radius (CRE) against a combination of products from satellite-based active and passive instruments. The OCA CTP is in good agreement with the CTP sensed by lidar for low thick liquid clouds and substantially below in the case of high ice clouds, in agreement with previous studies. The retrievals of COT and CRE are more reliable when constrained by solar channels and are consistent with other retrievals from passive imagers. The resulting cloud properties are stable and homogeneous over the whole period when compared against similar CDRs from passive instruments. For CTP, the OCA CDR and the near-real-time OCA products are consistent, allowing for the use of OCA near-real time products to extend the CDR beyond August 2019. Full article
(This article belongs to the Special Issue Satellite-Based Cloud Climatologies)
Show Figures

Figure 1

Figure 1
<p>Spatial coverage of Meteosat SSP at 0° exploited for OCA up to 65 degrees. The area covered also by the backup platforms at 3.4°W and at 9.5°E are shown. The common retrieved area along the whole data record time coverage is shown by the grey hatched area.</p>
Full article ">Figure 2
<p>Illustrative vertical “section” view of the OCA product. Upper panel is a section of the SEVIRI image in “Natural colour” RGB and the pixels constituting the lower edge of this image are shown in the lower panel as the OCA retrieval: ice-phase clouds are shown blue, liquid-phase clouds in green. For purely visualization purposes, simple conversion factors were used to derive CTH (height in km) from the CTP retrievals and cloud geometrical thicknesses from the COT retrievals. Note that CRE, error estimates, and quality flags are not illustrated in this figure.</p>
Full article ">Figure 3
<p>Hovmöller diagram of zonal mean (+60° to −60°) monthly cloud-top pressure (hPa). The figure includes only filtered data using a threshold at 60% valid retrievals.</p>
Full article ">Figure 4
<p>Monthly Hovmöller diagrams for zonal mean (+60° to −60°) cloud optical thickness (<b>upper panel</b>) and cloud particle effective radius (<b>bottom panel</b>). No threshold on the available number of days per month is applied.</p>
Full article ">Figure 5
<p>(<b>a</b>) Monthly fraction of one-layer water clouds (upper left sub-panel), ice clouds (upper right sub-panel), two-layers clouds (lower left sub-panel), and valid percentage of daily values per month (lower right sub-panel). (<b>b</b>) Monthly averages of cloud-top pressure (upper left sub-panel), cloud effective radius (upper right sub-panel), cloud optical thickness (lower left sub-panel), and valid percentage of daily values per month (lower right sub-panel). The example is for January 2018 (Meteosat-10). A threshold of 60% valid retrievals is applied to all values.</p>
Full article ">Figure 6
<p>MODIS true color RGB image of the scene analyzed in <a href="#remotesensing-16-02989-f007" class="html-fig">Figure 7</a>. The blue line shows the track of AQUA ascending orbit. The image is rotated by 90° and the left side is at 57S latitude; the right side at 22S latitude.</p>
Full article ">Figure 7
<p>Example of identification of cloud tops from DARDAR and CALIPSO data for the granule section of <a href="#remotesensing-16-02989-f003" class="html-fig">Figure 3</a>. Red dots show the cloud top inferred from the ATrain data, green dots show the OCA cloud top where the identification of single- or two-layer agrees with the one from DARDAR data. The (<b>top panel</b>) shows the cloud top for single layers only, the (<b>middle panel</b>) the cloud top of the upper layer of two-layer profiles, and the (<b>lower panel</b>) the cloud top of the second layer of two-layer profiles. The orange dots in the first two panels show where the identification from OCA does not agree with DARDAR—OCA two-layer in the first panel and OCA single-layer in the middle panel.</p>
Full article ">Figure 8
<p>Time series of cloud-top height of single-layer clouds from OCA (continuous line) and DARDAR (dashed line) for daytime orbits. (<b>Top panel</b>): mean for collocated products within one granule crossing the SEVIRI disk. The colored shading shows one standard deviation of DARDAR data. (<b>Lower panel</b>): mean and standard deviation of the differences between OCA and DARDAR cloud-top height. Data are divided between ice (blue lines) and liquid (green lines) clouds.</p>
Full article ">Figure 9
<p>Same as <a href="#remotesensing-16-02989-f008" class="html-fig">Figure 8</a> but for two-layer pixels. Data are divided between upper layer (blue lines) and second layer (green lines).</p>
Full article ">Figure 10
<p>Scatter plots of cloud-top height (km) retrieved by OCA and DARDAR for all the collocated pixels in the daytime granules. (<b>Upper panels</b>): single-layer pixels ice (left) and liquid (right). (<b>Lower panels</b>): two-layer pixels top layer (left) and second layer (right).</p>
Full article ">Figure 11
<p>Scatter plots of ice cloud optical thickness and effective radius retrieved by OCA and DARDAR for all the collocated pixels in the daytime. (<b>Upper panels</b>): cloud optical thickness for pixels identified as single-layer (left) and two-layer (right). (<b>Lower panels</b>): cloud-top effective radius for pixels identified as single-layer (left) and two-layer (right). For the two-layer pixels, only the total COT from DARDAR where both upper and lower layers are of ice type are used.</p>
Full article ">Figure 12
<p>Comparison of cloud categorization in OCA and DARDAR (DD). The categories are single-layer ice (SL_i), single-layer liquid (SL_l), and multi-layer (ML). The numbers and colors refer to the percentage of cases with respect to the total in that category (blue = ice; green = liquid; red = multi-layer). The comparison is carried out for the daytime A-Train orbits collocated with SEVIRI with the quality control filter applied as explained in the text. In each column, the percentage of cases of each OCA category is reported with respect to the DARDAR “truth”, both in number and shades of color. For example, in the category identified by DARDAR as single-layer ice (DD_SL_i), OCA agrees in 63.3% of cases while in 22.1% of cases, OCA detects a multi-layer situation (OCA_ML).</p>
Full article ">Figure 13
<p>Winter (December–January–February) seasonal mean combined (ice–liquid water clouds) CTP from OCA, MODIS, CLAAS-3, and CALIPSO GEWEX L3 datasets. On the bottom the difference between OCA and respectively MODIS, CLAAS-3, and CALIPSO GEWEX L3 (“<span class="html-italic">passive CTP flavour</span>”). White areas in the OCA and CLAAS-3 figures indicate regions where not enough valid retrievals were available (see text for details).</p>
Full article ">Figure 14
<p>Weighted area average of retrieved cloud-top pressure (CTP) from OCA, MODIS, CLAAS-3, and CALIPSO L3 GEWEX (“passive CTP flavor”) datasets over a SEVIRI disk. Averages are determined for the areas within a maximum SEVIRI viewing angle of 50°. (<b>Top panel</b>): monthly means with a lowess (locally weighted scatter plot smooth) smoothing filter applied to each time series. (<b>Lower panel</b>): differences between OCA and, respectively, MODIS, CLAAS-3, and CALIPSO L3.</p>
Full article ">Figure 15
<p>Same as <a href="#remotesensing-16-02989-f014" class="html-fig">Figure 14</a> but using CALIPSO L3 GEWEX CTP unadjusted (“TopLayer flavor”).</p>
Full article ">Figure 16
<p>Winter (December–January–February) seasonal mean combined (ice–liquid water clouds) COT from OCA, MODIS, and CLAAS-3. On the bottom, the difference between OCA and, respectively, MODIS and CLAAS-3. White areas in the OCA and CLAAS-3 figures indicate regions where not enough valid retrievals were available (see text for details).</p>
Full article ">Figure 17
<p>Weighted area average of retrieved cloud optical thickness (COT) from OCA, MODIS, and CLAAS-3 datasets over a SEVIRI disk. Averages are determined for the areas within a maximum SEVIRI viewing angle of 50°. (<b>Top panel</b>): monthly means with a lowess (locally weighted scatter plot smooth) smoothing filter applied to each time series. (<b>Lower panel</b>): differences between OCA and, respectively, MODIS and CLAAS-3.</p>
Full article ">Figure 18
<p>Winter (December–January–February) seasonal mean cloud particle effective radius (CRE) (ice and liquid) from OCA (single-layer and two-layer), MODIS, and CLAAS-3 datasets. On the bottom, the difference between OCA and, respectively, MODIS and CLAAS-3. White areas in the OCA and CLAAS-3 figures indicate regions where not enough valid retrievals were available.</p>
Full article ">Figure 19
<p>Same as <a href="#remotesensing-16-02989-f018" class="html-fig">Figure 18</a> but for liquid clouds only. For the OCA dataset, areas with a fraction of liquid phase larger than 60% were selected. For CLAAS-3 and MODIS, the CRE for liquid clouds is shown.</p>
Full article ">Figure 20
<p>Weighted area average of retrieved ice water cloud effective radius (CRE) from OCA (single-layer and two-layer), MODIS, and CLAAS-3 datasets over a SEVIRI disk. Averages are determined for the areas within a maximum SEVIRI viewing angle of 50°. (<b>Top panel</b>): monthly means with a lowess (locally weighted scatter plot smooth) smoothing filter applied to each time series. (<b>Lower panel</b>): differences between OCA and, respectively, MODIS and CLAAS-3.</p>
Full article ">
24 pages, 6993 KiB  
Article
Advancing Volcanic Activity Monitoring: A Near-Real-Time Approach with Remote Sensing Data Fusion for Radiative Power Estimation
by Giovanni Salvatore Di Bella, Claudia Corradino, Simona Cariello, Federica Torrisi and Ciro Del Negro
Remote Sens. 2024, 16(16), 2879; https://doi.org/10.3390/rs16162879 - 7 Aug 2024
Viewed by 1498
Abstract
The global, near-real-time monitoring of volcano thermal activity has become feasible through thermal infrared sensors on various satellite platforms, which enable accurate estimations of volcanic emissions. Specifically, these sensors facilitate reliable estimation of Volcanic Radiative Power (VRP), representing the heat radiated during volcanic [...] Read more.
The global, near-real-time monitoring of volcano thermal activity has become feasible through thermal infrared sensors on various satellite platforms, which enable accurate estimations of volcanic emissions. Specifically, these sensors facilitate reliable estimation of Volcanic Radiative Power (VRP), representing the heat radiated during volcanic activity. A critical factor influencing VRP estimates is the identification of hotspots in satellite imagery, typically based on intensity. Different satellite sensors employ unique algorithms due to their distinct characteristics. Integrating data from multiple satellite sources, each with different spatial and spectral resolutions, offers a more comprehensive analysis than using individual data sources alone. We introduce an innovative Remote Sensing Data Fusion (RSDF) algorithm, developed within a Cloud Computing environment that provides scalable, on-demand computing resources and services via the internet, to monitor VRP locally using data from various multispectral satellite sensors: the polar-orbiting Moderate Resolution Imaging Spectroradiometer (MODIS), the Sea and Land Surface Temperature Radiometer (SLSTR), and the Visible Infrared Imaging Radiometer Suite (VIIRS), along with the geostationary Spinning Enhanced Visible and InfraRed Imager (SEVIRI). We describe and demonstrate the operation of this algorithm through the analysis of recent eruptive activities at the Etna and Stromboli volcanoes. The RSDF algorithm, leveraging both spatial and intensity features, demonstrates heightened sensitivity in detecting high-temperature volcanic features, thereby improving VRP monitoring compared to conventional pre-processed products available online. The overall accuracy increased significantly, with the omission rate dropping from 75.5% to 3.7% and the false detection rate decreasing from 11.0% to 4.3%. The proposed multi-sensor approach markedly enhances the ability to monitor and analyze volcanic activity. Full article
(This article belongs to the Special Issue Application of Remote Sensing Approaches in Geohazard Risk)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>(<b>1</b>) Derivation of the Normalized Thermal Index (NTI) obtained by combining the radiance of the MIR and the radiance of the TIR. (<b>2</b>) Application of the Spatial Standard Deviation (SSD) filter to each pixel in the image. (<b>3</b>) Definition of two statistical masks, Mask1 and Mask2, to identify “potential” and “true” hotspots, applied on the SSD and NTI of the volcanic area (VA). (<b>4</b>) Application of Gabor filter to extract the significant features of the image, resulting in a matrix called Gabor Weighted NTI (G-NTI). (<b>5</b>) Highlighting hotspots in the crater area and defining the Spatial Gabor Weighted NTI (SG-NTI). (<b>6</b>) Application of a statistical mask to the previously extracted matrix. (<b>7</b>) Calculation of the final VRP.</p>
Full article ">Figure 2
<p>Workflow image of the RSDF algorithm. Study cases: (<b>a</b>) Etna on 2 December 2023 at 01:10 UTC, MODIS sensor; (<b>b</b>) Etna on 15 January 2023 at 20:46 UTC, SLSTR sensor; (<b>c</b>) Stromboli on 3 October 2023 at 13:10 UTC, MODIS sensor; (<b>d</b>) Stromboli on 23 October 2023 at 09:08 UTC, SLSTR sensor.</p>
Full article ">Figure 3
<p>Time series of the Etna volcano. The panels show VRP calculated respectively from the RSDF Algorithm SLSTR (blue triangles) SLSTR Level 2 (red triangles), the RSDF Algorithm MODIS (blue triangles), MODIS Level 2 (red triangles). (<b>a</b>,<b>c</b>) shows data from January 2021 to April 2022, (<b>b</b>,<b>d</b>) from April 2022 to June 2023.</p>
Full article ">Figure 4
<p>Histograms (<b>a</b>,<b>c</b>) and probability plots (<b>b</b>,<b>d</b>) for Etna datasets. (<b>a</b>,<b>c</b>) Histograms display data distribution related to VRP (and FRP) in logarithmic scale; (<b>a</b>) blue bars represent the distribution of SLSTR-–RSDF algorithm processed data, and red bars represent SLSTR Level 2 product data; (<b>c</b>) blue bars represent the distribution of MODIS–RSDF algorithm processed data, and red bars represent MODIS active fire products. (<b>b</b>,<b>d</b>) Probability plots for normal distribution of RSDF algorithm processed data (blue), and Level 2 product data (red). The dashed grey lines represent the reference lines of the theoretical distributions, and the black dashed line in (<b>b</b>) corresponds to the slope change associated with the transition between regimes of background and high thermal activity.</p>
Full article ">Figure 5
<p>Stacked time series of VRPw (weekly mean) retrieved for the SLSTR–RSDF algorithm processed data (blue), SLSTR Level 2 product data (red), MODIS–RSDF algorithm processed data (green), and SLSTR Level 2 product data (black) at the Etna volcano, displayed on a logarithmic scale.</p>
Full article ">Figure 6
<p>VRP time series of the Stromboli volcano. The panels show VRP calculated respectively from the RSDF Algorithm SLSTR (blue triangles) SLSTR Level 2 (red triangles), the RSDF Algorithm MODIS (blue triangles), MODIS Level 2 (red triangles). (<b>a</b>,<b>c</b>) shows data from January 2021 to April 2022, (<b>b</b>,<b>d</b>) from April 2022 to June 2023.</p>
Full article ">Figure 7
<p>Histograms (<b>a</b>,<b>c</b>) and probability plots (<b>b</b>,<b>d</b>) for Stromboli datasets. (<b>a</b>,<b>c</b>) Histograms display data distribution related to VRP (and FRP) in logarithmic scale; (<b>a</b>) blue bars represent the distribution of SLSTR–RSDF algorithm processed data, and red bars represent SLSTR Level 2 product data; (<b>c</b>) blue bars represent the distribution of MODIS–RSDF algorithm processed data, and red bars represent MODIS active fire products. (<b>b</b>,<b>d</b>) Probability plots for normal distribution of RSDF algorithm processed data (blue), and Level 2 product data and MODIS active fire products (red). The dashed grey lines represent the reference lines of the theoretical distributions, and the black dashed line in (<b>b</b>) corresponds to the slope change associated with the transition between regimes of background and high thermal activity for Stromboli.</p>
Full article ">Figure 8
<p>Stacked time series of VRPw (weekly mean) retrieved for SLSTR–RSDF algorithm processed data (blue), SLSTR Level 2 product data (red), MODIS–RSDF algorithm processed data (green), and MODIS active fire products (black) at the Etna volcano, displayed on a logarithmic scale.</p>
Full article ">Figure 9
<p>Cumulative Volcanic Radiative Energy (VRE) calculated from VRP (and FRP) using the trapezoidal rule for integration. The blue line represents VRESLSTR, the red dashed line FREMODIS, the green dashed line VREMODIS, and the black dashed line FREMODIS. Panels (<b>a</b>,<b>c</b>) show data for Etna; panels (<b>b</b>,<b>d</b>) show data for Stromboli.</p>
Full article ">Figure 10
<p>Radiative power time series from SLSTR– and MODIS–RSDF algorithm data with intensity limits categorized as low, moderate, high, and extreme. (<b>a</b>) Etna, (<b>b</b>) Stromboli.</p>
Full article ">Figure 11
<p>Temporal trend of VRP values derived from the RSDF algorithm for SEVIRI, SLSTR, MODIS, and VIIRS over two periods at Mt. Etna: (<b>a</b>) 1 February 2021–30 April 2021, and (<b>b</b>) 27 September 2023–10 October 2023.</p>
Full article ">Figure 12
<p>TADR and lava flow volume flux during the effusive event at Etna from 14 May 2022 to 16 June 2022. TADR_max, TADR_mean, and TADR_min are represented by blue, red, and green points, respectively. The total volume_max, volume_mean, and volume_min are represented by blue, red, and green lines, respectively.</p>
Full article ">Figure 13
<p>TADR and lava flow volume flux during the effusive event at Stromboli from 27 September 2023 to 10 October 2023. TADR_max, TADR_mean, and TADR_min are represented by blue, red, and green points, respectively. The total volume_max, volume_mean, and volume_min are represented by blue, red, and green lines, respectively.</p>
Full article ">
18 pages, 2310 KiB  
Article
Data Assimilation of Satellite-Derived Rain Rates Estimated by Neural Network in Convective Environments: A Study over Italy
by Rosa Claudia Torcasio, Mario Papa, Fabio Del Frate, Alessandra Mascitelli, Stefano Dietrich, Giulia Panegrossi and Stefano Federico
Remote Sens. 2024, 16(10), 1769; https://doi.org/10.3390/rs16101769 - 16 May 2024
Viewed by 873
Abstract
The accurate prediction of heavy precipitation in convective environments is crucial because such events, often occurring in Italy during the summer and fall seasons, can be a threat for people and properties. In this paper, we analyse the impact of satellite-derived surface-rainfall-rate data [...] Read more.
The accurate prediction of heavy precipitation in convective environments is crucial because such events, often occurring in Italy during the summer and fall seasons, can be a threat for people and properties. In this paper, we analyse the impact of satellite-derived surface-rainfall-rate data assimilation on the Weather Research and Forecasting (WRF) model’s precipitation prediction, considering 15 days in summer 2022 and 17 days in fall 2022, where moderate to intense precipitation was observed over Italy. A 3DVar realised at CNR-ISAC (National Research Council of Italy, Institute of Atmospheric Sciences and Climate) is used to assimilate two different satellite-derived rain rate products, both exploiting geostationary (GEO), infrared (IR), and low-Earth-orbit (LEO) microwave (MW) measurements: One is based on an artificial neural network (NN), and the other one is the operational P-IN-SEVIRI-PMW product (H60), delivered in near-real time by the EUMETSAT HSAF (Satellite Application Facility in Support of Operational Hydrology and Water Management). The forecast is verified in two periods: the hours from 1 to 4 (1–4 h phase) and the hours from 3 to 6 (3–6 h phase) after the assimilation. The results show that the rain rate assimilation improves the precipitation forecast in both seasons and for both forecast phases, even if the improvement in the 3–6 h phase is found mainly in summer. The assimilation of H60 produces a high number of false alarms, which has a negative impact on the forecast, especially for intense events (30 mm/3 h). The assimilation of the NN rain rate gives more balanced predictions, improving the control forecast without significantly increasing false alarms. Full article
(This article belongs to the Special Issue Remote Sensing of Extreme Weather Events: Monitoring and Modeling)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>WRF model domain and rain gauge stations’ altitude (<b>a</b>); simulations’ scheme using the very short-term forecast approach (<b>b</b>).</p>
Full article ">Figure 2
<p>Rain rates estimated by the neural network (in millimetres per hour) (<b>a</b>), and the relative humidity innovation fields at 2094 m a.s.l. at 06 UTC on 18 August 2022 for the NN_3th (<b>b</b>) and the NN_1th (<b>c</b>) model configurations (in percentages).</p>
Full article ">Figure 3
<p>Rain rates estimated by H60 (<b>a</b>) and relative humidity field innovation (<b>b</b>) at 06 UTC on 18 August 2022 at 2094 m a.s.l. for the H60_3th configuration.</p>
Full article ">Figure 4
<p>Precipitation between 07 and 10 UTC on 18 August 2022, as reported by rain gauges (<b>a</b>) and MCM (<b>b</b>) and simulated by CTRL (<b>c</b>), H60_3th (<b>d</b>), NN_3th (<b>e</b>), and NN_1th (<b>f</b>) WRF model configurations.</p>
Full article ">Figure 5
<p>Performance diagrams for the CTRL, NN_3th, NN_1th, and H60_3th configurations and for the 1 mm/3 h, 10 mm/3 h, and 30 mm/3 h thresholds for the case studies in summer (<b>a</b>) and in fall (<b>b</b>) for the 1–4 h verification phase. Red hyperboles branches represent the threat score while cyan lines represent the frequency bias.</p>
Full article ">Figure 6
<p>Performance diagrams for the CTRL, NN_3th, NN_1th, and H60_3th configurations and for the 1 mm/3 h, 10 mm/3 h, and 30 mm/3 h thresholds for the case studies in summer (<b>a</b>) and in fall (<b>b</b>) for the 3–6 h verification phase. Red hyperboles branches represent the threat score while cyan lines represent the frequency bias.</p>
Full article ">
22 pages, 9756 KiB  
Article
Investigation of the Synoptic and Dynamical Characteristics of Cyclone Shaheen (2021) and Its Influence on the Omani Coastal Region
by Petros Katsafados, Pantelis-Manolis Saviolakis, George Varlas, Haifa Ben-Romdhane, Kosmas Pavlopoulos, Christos Spyrou and Sufian Farrah
Atmosphere 2024, 15(2), 222; https://doi.org/10.3390/atmos15020222 - 12 Feb 2024
Viewed by 1234
Abstract
Tropical Cyclone Shaheen (TCS), originating in the Arabian Sea on 30 September 2021, followed an east-to-west trajectory and made landfall as a category-1 cyclone in northern Oman on 3 October 2021, causing severe floods and damages before dissipating in the United Arab Emirates. [...] Read more.
Tropical Cyclone Shaheen (TCS), originating in the Arabian Sea on 30 September 2021, followed an east-to-west trajectory and made landfall as a category-1 cyclone in northern Oman on 3 October 2021, causing severe floods and damages before dissipating in the United Arab Emirates. This study aims to analyze the synoptic and dynamical conditions influencing Shaheen’s genesis and evolution. Utilizing ERA5 reanalysis data, SEVIRI-EUMETSAT imagery, and Sorbonne University Atmospheric Forecasting System (SUAFS) outputs, it was found that Shaheen manifested as a warm-core cyclone with moderate vertical wind shear within the eyewall. Distinctive features included a trajectory aligned with rising sea surface temperatures and increased specific humidity levels at 700 hPa in the Arabian Sea. As Shaheen approached the Gulf of Oman, a significant increase in rainfall rates occurred, correlated with variations in sea surface temperatures and vertical wind shear. Comparative analysis between SUAFS and ERA5 data revealed a slight northward shift in the SUAFS track and landfall. Advance warnings highlighted heavy rainfall, rough seas, and strong winds. This study provides valuable insights into the meteorological factors contributing to Shaheen’s formation and impact. Full article
(This article belongs to the Section Meteorology)
Show Figures

Figure 1

Figure 1
<p>The computational domain of SUAFS. The topography (m) of the area is also illustrated.</p>
Full article ">Figure 2
<p>Tracks of cyclones and depressions in the North Indian Ocean for the year 2021. (Source: Cyclone e-atlas RSMC India Meteorological Department).</p>
Full article ">Figure 3
<p>The development of Shaheen through satellite images in infrared. (<b>a</b>) is on 30/09/2021 at 12:00 UTC, (<b>b</b>) on 01/10/2021 at 00:00 UTC, (<b>c</b>) on 02/10/2021 at 00:00 UTC, (<b>d</b>) on 03/10/2021 at 00:00 UTC (Source: EUMETSAT Images, High Rate SEVIRI IR10.8μm).</p>
Full article ">Figure 4
<p>(<b>a</b>) The SST (K) in the region on 30/09/2021 at 12:00 UTC. (<b>b</b>) The SST (K) on 03/10/2021 at 12:00 UTC. (<b>c</b>) The total precipitation (mm/h) on 03/10/2021 12:00 UTC (Source Copernicus ERA5 data).</p>
Full article ">Figure 5
<p>Wind shear between 200−850 hPa (in m s<sup>−1</sup>) on 30/09/2021 at 12:00 UTC.</p>
Full article ">Figure 6
<p>Phase Space Hart Diagram for TCS (estimated using Copernicus ERA5 Data). The time interval between successive points is 6 h from 30/9/2021 12:00 UTC up to 03/10/2021 12:00 UTC.</p>
Full article ">Figure 7
<p>(<b>a</b>) Temperature at 850 hPa (K) and Geopotential height (m) at 500 hPa on 30/09/2021 12:00 UTC. (<b>b</b>) Temperature at 850 hPa (K) and Geopotential height (m) at 500 hPa on 01/10/2021 at 06:00 UTC. (<b>c</b>) Vertical velocity in a horizontal cross section of the cyclone and constant Lat 23.25° on 01/10/2021 at 06:00 UTC. (Source: Copernicus ERA5 data).</p>
Full article ">Figure 8
<p>(<b>a</b>) Specific humidity (kg/kg) and wind vectors at 700 hPa on 30/09/2021 12:00 UTC Based on Copernicus ERA5 data. (<b>b</b>) Potential vorticity at 850 hPa (1 PV = 10<sup>−6</sup> m<sup>2</sup> s<sup>−1</sup> K kg<sup>−1</sup>) and geopotential height (gpm) on 30/09/2021 12:00 UTC. (Source Copernicus ERA5 data).</p>
Full article ">Figure 9
<p>Wind shear between 200–850 hPa (in m s<sup>−1</sup>) from 30/09/2021 12:00 UTC to 03/10/2021 00:00 UTC. (<b>a</b>) 30/09/2021 12:00 UTC, (<b>b</b>) 01/10/2021 00:00 UTC, (<b>c</b>) 02/10/2021 00:00 UTC, (<b>d</b>) 03/10/2021 00:00 UTC. (Source Copernicus ERA5 data).</p>
Full article ">Figure 10
<p>IVT analysis in kg m<sup>−1</sup> s<sup>−1</sup>. (<b>a</b>) on 30/09/2021 12:00 UTC, (<b>b</b>) on 01/10/2021 00:00 UTC, (<b>c</b>) on 02/10/2021 00:00 UTC, (<b>d</b>) on 03/10/2021 00:00 UTC (Source Copernicus ERA5 data).</p>
Full article ">Figure 11
<p>(<b>a</b>) The locations and the WMO stations (showing the WMO ID number) used for the comparison against SUAFS model outputs, (<b>b</b>) Time plots of the modeled and observed temperature at 2 m, (<b>c</b>) Time plots of the modeled and observed wind speed at 10 m, (<b>d</b>) Time plots of the modeled and observed mean sea level pressure.</p>
Full article ">Figure 12
<p>Wind speed at 10 m (m s<sup>−1</sup>) from SUAFS (<b>a</b>) and from Copernicus ERA5 (<b>b</b>) on 03/10/2021 at 9:00 UTC. The vectors depict the wind direction at 10m.</p>
Full article ">Figure 13
<p>Simulated (<b>a</b>) latent heat flux and (<b>b</b>) sensible heat flux from SUAFS in W m<sup>−2</sup> on 03/10/2021 at 9:00 UTC. The vectors depict the wind direction at 10m.</p>
Full article ">Figure 14
<p>(<b>a</b>) Vertical velocity, (<b>b</b>) vertical profile of condensational heating rate, (<b>c</b>) vertical profile of potential vorticity on 03/10/2021 at 9:00 UTC, as simulated by SUAFS, (<b>d</b>) Potential Vorticity at 850 hPa from the SUAFS, on 03/10/2021 09:00 UTC in Potential Vorticity Units (1 PVU = 10<sup>−6</sup> m<sup>2</sup> s<sup>−1</sup>⋅K kg<sup>−1</sup>).</p>
Full article ">Figure 15
<p>The track of cyclone Shaheen as extracted from Copernicus ERA5 data (red) and the SUAFS model (blue). Dots correspond to 3 h time intervals.</p>
Full article ">Figure 16
<p>The drainage basin of the wadi Hawasnah and the wadi Bani (red line), where the fieldwork took place during November 2021.</p>
Full article ">Figure 17
<p>(<b>Top</b>): The flood plain of Wadi Hawasnah one month after the cyclone Shaheen passed as viewed from eastwards. (<b>Bottom</b>): Drainage basin as viewed from downstream/westwards. The wadi bedload deposits destroyed the main road.</p>
Full article ">
15 pages, 4077 KiB  
Article
Diurnal Variation in Urban Heat Island Intensity in Birmingham: The Relationship between Nocturnal Surface and Canopy Heat Islands
by Cong Wen, Ali Mamtimin, Jiali Feng, Yu Wang, Fan Yang, Wen Huo, Chenglong Zhou, Rui Li, Meiqi Song, Jiacheng Gao and Ailiyaer Aihaiti
Land 2023, 12(11), 2062; https://doi.org/10.3390/land12112062 - 13 Nov 2023
Cited by 5 | Viewed by 1308
Abstract
Urban heat islands have garnered significant attention due to their potential impact on human life. Previous studies on urban heat islands have focused on characterizing temporal and spatial variations over longer periods of time. In this study, we investigated the urban heat island [...] Read more.
Urban heat islands have garnered significant attention due to their potential impact on human life. Previous studies on urban heat islands have focused on characterizing temporal and spatial variations over longer periods of time. In this study, we investigated the urban heat island (UHI) in Birmingham from September 2013 to August 2014 using higher temporal resolution SEVIRI satellite surface temperature data along with data from the Birmingham Urban Climate Laboratory (BUCL) meteorological station and the UK Meteorological Office meteorological station. Our aim was to characterize the diurnal variations in the surface urban heat island intensity (SUHII) and canopy urban heat island intensity (CUHII) and to explore their relationship under the influence of three factors (day/nighttime, season, and wind speed) using regression analysis. Our findings reveal that SUHII and CUHII exhibit relatively stable patterns at night but vary significantly during the day with opposite diurnal trends. In addition, SUHII and CUHII were more variable in spring and summer but less variable in winter. During the nighttime, SUHII represents CUHII with high confidence, especially during spring and summer, but less so during the cold season. In addition, SUHII represents CUHII with greater confidence under low-wind conditions. This study deepens our understanding of the diurnal dynamics of urban heat islands and the influence of atmospheric conditions on the relationship between surface and canopy heat islands in urban areas. The results of this study can be used for heat island studies in cities that lack high-precision observation networks and to guide sustainable urban development. Full article
Show Figures

Figure 1

Figure 1
<p>Birmingham city weather station location.</p>
Full article ">Figure 2
<p>Diurnal variations in SEVIRI LST, MODIS LST, and air temperature at BUCL on 3 September 2013.</p>
Full article ">Figure 3
<p>Correlation analysis of SEVIRI LST data with BUCL observation data.</p>
Full article ">Figure 4
<p>Diurnal variation in SUHII throughout the year.</p>
Full article ">Figure 5
<p>Diurnal variation in CUHII throughout the year.</p>
Full article ">Figure 6
<p>Linear and elliptical trends of SUHII and CUHII based on the daytime and nighttime: (<b>a</b>) daytime; (<b>b</b>) nighttime.</p>
Full article ">Figure 7
<p>Linear and elliptical trends of SUHII and CUHII based on three wind speed conditions: (<b>a</b>) WS1, (<b>b</b>) WS2, and (<b>c</b>) WS3.</p>
Full article ">Figure 8
<p>Linear and elliptical trends of SUHII and CUHII based on the season: (<b>a</b>) spring, (<b>b</b>) summer, (<b>c</b>) autumn, and (<b>d</b>) winter.</p>
Full article ">
25 pages, 6080 KiB  
Article
Specific Features of the Land-Sea Contrast of Cloud Liquid Water Path in Northern Europe as Obtained from the Observations by the SEVIRI Instrument: Artefacts or Reality?
by Vladimir S. Kostsov and Dmitry V. Ionov
Meteorology 2023, 2(4), 464-488; https://doi.org/10.3390/meteorology2040027 - 11 Nov 2023
Viewed by 1202
Abstract
Liquid water path (LWP) is one of the most important cloud parameters and is crucial for global and regional climate modelling, weather forecasting, and modelling of the hydrological cycle and interactions between different components of the climate system: the atmosphere, the hydrosphere, and [...] Read more.
Liquid water path (LWP) is one of the most important cloud parameters and is crucial for global and regional climate modelling, weather forecasting, and modelling of the hydrological cycle and interactions between different components of the climate system: the atmosphere, the hydrosphere, and the land surface. Space-borne observations by the SEVIRI instrument have already provided evidence of the systematic difference between the cloud LWP values derived over the land surface in Northern Europe and those derived over the Baltic Sea and major lakes during both cold and warm seasons. In the present study, the analysis of this LWP land-sea contrast for the period 2011–2017 reveals specific temporal and spatial variations, which, in some cases, seem to be artefacts rather than of natural origin. The geographical objects of investigation are water bodies and water areas located in Northern Europe that differ in size and other geophysical characteristics: the Gulf of Finland and the Gulf of Riga in the Baltic Sea and large and small lakes in the neighbouring region. The analysis of intra-seasonal features has detected anomalous conditions in the Gulf of Riga and the Gulf of Finland, which show up as very low values of the LWP land-sea contrast in August with respect to the values in June and July every year within the considered time period. This anomaly is likely an artefact caused by the LWP retrieval algorithm since the transition from large LWP contrast to very low contrast occurs sharply, synchronically, and at a certain date every year at different places in the Baltic Sea. Full article
Show Figures

Figure 1

Figure 1
<p>The geographical domain and the water bodies under investigation: The Gulf of Riga (1), The Gulf of Finland (2), Lake Ladoga (3), Lake Onega (4), Lake Peipus (5), Lake Pihkva (6), Lake Ilmen (7), Lake Saimaa (8), and the Neva River Bay (9). This example map demonstrates the SEVIRI-derived mean LWP for July (2011–2017). Vector shoreline data: [<a href="#B27-meteorology-02-00027" class="html-bibr">27</a>]. The colour scale represents the LWP in kg m<sup>−2</sup>.</p>
Full article ">Figure 2
<p>The geographical domain and the specific locations of the LWP contrast measurements (shown as lines that connect ground pixels of the SEVIRI measurements, red—large distance LWP contrast, blue and green—small distance LWP contrast). Green colour is used to mark additional measurement locations in the Gulf of Finland, see text. The colour letters with numbers indicate the corresponding data sets (see <a href="#meteorology-02-00027-t001" class="html-table">Table 1</a> and <a href="#meteorology-02-00027-t002" class="html-table">Table 2</a>). Vector shoreline data: [<a href="#B27-meteorology-02-00027" class="html-bibr">27</a>].</p>
Full article ">Figure 3
<p>Typical statistical distributions (in terms of relative frequency of occurrence <span class="html-italic">R</span>) of the LWP land-sea contrast values for different measurement locations and seasons. Please note that for better visibility, the vertical axes are broken and have different scaling in the lower and upper parts.</p>
Full article ">Figure 4
<p>Inter-annual variability of the seasonal mean LWP land-sea contrast for measurement locations ML1–ML4 and different seasons (<b>a</b>–<b>d</b>). Green dots indicate the zero contrast.</p>
Full article ">Figure 5
<p>Inter-annual variations of the seasonal mean LWP land-sea contrast for measurement locations ML5–ML8 (<b>a</b>,<b>b</b>) and ML9 (<b>c</b>,<b>d</b>) and different seasons. Green dots indicate the zero contrast. Panels (<b>c</b>,<b>d</b>) demonstrate not only the seasonal mean values but also their standard errors.</p>
Full article ">Figure 6
<p>Intra-seasonal variability of the daily mean LWP land-sea contrast for measurement locations ML1-2, ML2-2, and ML9 (warm season, seven years of observations—see the legend).</p>
Full article ">Figure 7
<p>Intra-seasonal variability of the monthly mean LWP land-sea contrast (warm season, seven years of observations altogether).</p>
Full article ">Figure 8
<p>Same as <a href="#meteorology-02-00027-f007" class="html-fig">Figure 7</a> but for additional measurement locations.</p>
Full article ">Figure 9
<p>Intra-seasonal variability of the daily mean LWP land-sea contrast for all measurement locations in the Gulf of Riga and the Gulf of Finland. Averaging of daily mean values was done over all seven years of observations. Additionally, running averaging over 3 days was used for better visibility of the results.</p>
Full article ">Figure 10
<p>The map showing the geographical location of the ERA5 reanalysis grid points used for the calculations of the LWP land-sea contrast. Vector shoreline data: [<a href="#B27-meteorology-02-00027" class="html-bibr">27</a>].</p>
Full article ">Figure 11
<p>Typical statistical distributions (in terms of relative frequency of occurrence <span class="html-italic">R</span>) of the LWP land-sea contrast values for the RE2-1 location and two seasons as obtained from the ERA-Interim data (<b>a</b>) and the ERA5 data (<b>b</b>). (When comparing with <a href="#meteorology-02-00027-f003" class="html-fig">Figure 3</a>, please note that for better visibility, the vertical axes in <a href="#meteorology-02-00027-f003" class="html-fig">Figure 3</a> are broken and have different scaling in the lower and upper parts).</p>
Full article ">Figure 12
<p>(<b>a</b>–<b>h</b>) Inter-annual variability of the seasonal mean LWP land-sea contrast (large distance) for cold and warm seasons as derived from the SEVIRI measurements and the reanalysis data (ERA-Interim and ERA5, see the legends). Green dots indicate the zero-contrast line.</p>
Full article ">Figure 13
<p>Intra-seasonal variability of the monthly mean LWP land-sea contrast as obtained from the ERA-Interim and ERA5 reanalyses (warm season, 2011–2017) and the SEVIRI data at the locations, which correspond to the large-distance LWP differences (RE…/ML…). Symbols are connected by lines only for illustrative purposes.</p>
Full article ">Figure 14
<p>Intra-seasonal (June–August) variability of the daily mean LWP land-sea contrast obtained from the reanalyses ERA5 (<b>a</b>,<b>c</b>) and ERA-Interim (<b>b</b>,<b>d</b>) for locations RE1-1 (<b>a</b>,<b>b</b>) and RE4-1 (<b>c</b>,<b>d</b>).</p>
Full article ">
15 pages, 3421 KiB  
Technical Note
A Year of Volcanic Hot-Spot Detection over Mediterranean Europe Using SEVIRI/MSG
by Catarina Alonso, Rita Durão and Célia M. Gouveia
Remote Sens. 2023, 15(21), 5219; https://doi.org/10.3390/rs15215219 - 3 Nov 2023
Viewed by 1317
Abstract
Volcano eruption identification and watching is crucial to better understanding volcano dynamics, namely the near real-time identification of the eruption start, end, and duration. Eruption watching allows hazard assessment, eruption forecasting and warnings, and also risk mitigation during periods of unrest, to enhance [...] Read more.
Volcano eruption identification and watching is crucial to better understanding volcano dynamics, namely the near real-time identification of the eruption start, end, and duration. Eruption watching allows hazard assessment, eruption forecasting and warnings, and also risk mitigation during periods of unrest, to enhance public safety and reduce losses from volcanic events. The near real-time fire radiative power (FRP) product retrieved using information from the SEVIRI sensor onboard the Meteosat Second Generation (MSG) satellite are used to identify and follow up volcanic activity at the pan-European level, namely the Mount Etna and Cumbre Vieja eruptions which occurred during 2021. The FRP product is designed to record information on the location, timing, and fire radiative power output of wildfires. Measuring FRP from SEVIRI/MSG and integrating it over the lifetime of a fire provides an estimate of the total Fire Radiative Energy (FRE) released. Together with FRP data analysis, SO2 data from the Copernicus Atmosphere Monitoring Service (CAMS) is used to assess the relationship between daily emitted concentrations of SO2 and the radiative energy released during volcanic eruptions. Results show that the FRE data allows us to evaluate the amount of energy released and is related to the pollutant concentrations from volcanic emissions during the considered events. A good agreement between FRP detection and SO2 atmospheric concentrations was found for the considered eruption occurrences. The adopted methodology, due to its simplicity and near real-time availability, shows potential to be used as a management tool to help authorities monitor and manage resources during ongoing volcanic events. Full article
(This article belongs to the Special Issue Earth Observation Using Satellite Global Images of Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Study area inside the dashed box: Mount Etna (<b>right panel</b>) and the Cumbre Vieja (<b>left panel</b>). Elevation map is represented by gray colors (unit: meters).</p>
Full article ">Figure 2
<p>Daily FRE (GJ) accumulated over the dashed box for Mount Etna in <a href="#remotesensing-15-05219-f001" class="html-fig">Figure 1</a>, covering the total period (<b>top panel</b>) and the hourly FRE (GJ) for March (<b>left bottom panel</b>) and July (<b>right bottom panel</b>).</p>
Full article ">Figure 3
<p>Daily FRE (blue, GJ) and SO<sub>2</sub> (orange, µg/m<sup>3</sup> accumulated over the dashed box for Mount Etna in <a href="#remotesensing-15-05219-f001" class="html-fig">Figure 1</a>.</p>
Full article ">Figure 4
<p>FRE (GJ) and SO<sub>2</sub> (µg/m<sup>3</sup>) concentrations on 4 March (<b>left panel</b>) and FRE and SO<sub>2</sub> concentrations on 1 April (<b>right panel</b>). The dashed rectangle represents the area where the accumulated FRP, FRE, and SO<sub>2</sub> previously represented were calculated for Mount Etna.</p>
Full article ">Figure 5
<p>Daily FRE (GJ) from September to December (<b>top panel</b>) and hourly FRE (GJ) from 28 September to 12 October (<b>bottom panel</b>) for the Cumbre Vieja.</p>
Full article ">Figure 6
<p>Daily FRE (blue, GJ) and SO<sub>2</sub> (orange, µg/m<sup>3</sup>) for the same area and period for the Cumbre Vieja.</p>
Full article ">Figure 7
<p>As in <a href="#remotesensing-15-05219-f003" class="html-fig">Figure 3</a>, but for the Cumbre Vieja on (<b>left panel</b>) 2 October and on (<b>right panel</b>) 17 October 2021.</p>
Full article ">Figure 8
<p>Scater plot of FRE greater than 1 GJ and SO<sub>2</sub> greater than 200 µg/m<sup>3</sup>, values accumulated over 8 days between January and April for Mount Etna event (blue), and September and December for the Cumbre Vieja event (orange).</p>
Full article ">
23 pages, 9047 KiB  
Article
Flash Drought and Its Characteristics in Northeastern South America during 2004–2022 Using Satellite-Based Products
by Humberto Alves Barbosa
Atmosphere 2023, 14(11), 1629; https://doi.org/10.3390/atmos14111629 - 30 Oct 2023
Cited by 7 | Viewed by 4299
Abstract
The term flash drought describes a special category of drought with rapid onset and strong intensity over the course of days or weeks. To characterize the impact of flash droughts on vegetation coverage, this study assessed the influence of soil water deficits on [...] Read more.
The term flash drought describes a special category of drought with rapid onset and strong intensity over the course of days or weeks. To characterize the impact of flash droughts on vegetation coverage, this study assessed the influence of soil water deficits on vegetation dynamics in the northeastern South America region by combining time series of vegetation index, rainfall, and soil moisture based on satellite products at a daily time scale. An 18-year analysis, from 2004 to 2022, of the Normalized Difference Vegetation Index (NDVI), Standard Precipitation Index (SPI), and surface soil moisture (SSM) was performed based on three different satellite remote sensing estimates: the spinning enhanced visible and infrared imager (SEVIRI) and the integrated multi-satellite retrievals algorithm (IMERG), and the soil moisture and ocean salinity (SMOS). The results revealed that flash drought events exerted dramatic impacts on terrestrial ecosystems in the study region during the first two decades of the 2000s, with changes in seasonal and regional vegetation dynamics. Further, the fixed-threshold values to characterize flash drought events were suggested as the timing when the water deficit was less than −1.0 units and vegetation index reached the value equal to +0.3 during five consecutive weeks or more, coupled with soil moisture rates below 40% percentile, leading to a strong region-wide drought throughout the entire region. Additionally, the results of linear least squares trend analyses revealed a negative trend in the pentad-SEVIRI radiance for the solar channel 1 within the semiarid ecosystems of the study region (i.e., the Caatinga biome) that was suggested as a reduction in clouds in the 18 years of the study. Developing combined threshold measures of flash drought based on satellite remote sensing may lead to an accurate assessment of flash drought mitigation. Full article
(This article belongs to the Special Issue Climate Variability and Change in Brazil)
Show Figures

Figure 1

Figure 1
<p>Location map of the Caatinga biome and its geographic features (topography) within the northeastern South America region. It covers approximately 735,000 km<sup>2</sup> and comprises the following states: Alagoas (AL), Bahia (BA), Ceará (CE), Maranhão (MA), Paraíba (PB), Piauí (PI), Pernambuco (PE), Rio Grande do Norte (RN), and Sergipe (SE).</p>
Full article ">Figure 2
<p>Northeastern South American region: (<b>a</b>) its nine states, including the Caatinga biome, along with the spatial distribution of temporary and permanent crops in 2006 [<a href="#B29-atmosphere-14-01629" class="html-bibr">29</a>]; (<b>b</b>) Spatial distribution of average annual rainfall (mm) over the region based on the long-term means from 2004 to 2022 with data from IMERG.</p>
Full article ">Figure 3
<p>Near real-time SEVIRI data flow processing at Laboratório de Análise e Processamento de Imagens de Satélites (LAPIS: <a href="https://www.lapismet.com.br" target="_blank">https://www.lapismet.com.br</a>. Accessed on 25 September 2023).</p>
Full article ">Figure 4
<p>Procedure applied for reducing the high dimensionality of the NDVIa or SPI3 matrices.</p>
Full article ">Figure 5
<p>Spatially averaged NDVIa over the entire northeastern South America region grouped by the NDVI clusters for the period 1982–2012. The labels F1 to F8 located in top-right of each box indicates the decreasing order each NDVI cluster according to median calculated from spatially averaged NDVIa over the entire region. Each box shows the median and first and third quartiles, while the whiskers extend to the last values that are 1.5 times the inter-quartile range above or below the quartiles. The medians are equal to −0.115, −0.004, −0.554, 0.466, −0.301, −0.547, 0.149, and 0.140 by the clusters N1, N2, N3, N4, N5, N6, N7, and N8, respectively. Circles are outliers.</p>
Full article ">Figure 6
<p>As in <a href="#atmosphere-14-01629-f005" class="html-fig">Figure 5</a>, but here for the SPI clusters. The medians are equal to 0.400, −0.232, 0.781, −0.699, −0.132, −0.260, 0.155, and −0.169 by the clusters S1, S2, S3, S4, S5, S6, S7, and S8, respectively. Circles are outliers.</p>
Full article ">Figure 7
<p>Regional mean changes in annual mean rainfall (mm), annual mean air temperature (°C), and the impact of rainfall anomalies on the vegetation greenness response for northeastern South America region from 2004 to 2022. (<b>a</b>) Annual variability of mean rainfall amount (mm) over the 512 grid cells within study region. (<b>b</b>) Annual variability of mean air temperature (°C) over the 512 grid cells within study region. (<b>c</b>) Five-day averaged SEVIRI NDVI over the 512 grid cells within study region from 2004 to 2022. A drought pentad (5-day mean) is defined as when SPI-3 was less than −1.0, and a wet week is defined as when SPI-3 was greater than 1.0. Mean drought severity, defined as the mean SPI of drought period (SPI-3 &lt; −1.0). Each pentad SEVIRI NDVI is marked with a circle from 2004 to 2012. (<b>d</b>) Seasonal variations of SEVIRI NDVI profile during 2009 (wet year; the solid blue line) and 2012 (drought year; the solid yellow line) for study region for the entire period 2004–2022.</p>
Full article ">Figure 8
<p>Time series of 5-day averaged soil moisture in 2012 from SMOS-based SSM (m<sup>3</sup>/m<sup>3</sup>) over northeastern South America region. The pentad is the mean over 5 days of consecutive soil moisture. The orange and red dashed lines denote the 40th–20th percentile range of soil moisture values.</p>
Full article ">Figure 9
<p>Regional comparison of (<b>a</b>) averaged soil moisture from 4 June 2012 to 28 June 2012 and (<b>b</b>) averaged soil moisture from 9 July 2012 to 31 October 2012 for major flash drought events over the period 2004–2022. Blue line area within the northeastern South America region highlights its semi-arid domain (Caatinga biome).</p>
Full article ">Figure 10
<p>Regional comparison of (<b>a</b>) averaged SEVIRI NDVI from 4 June 2012 to 28 June 2012 and (<b>b</b>) averaged SEVIRI NDVI from 9 July 2012 to 31 October 2012 for major flash drought events over the period 2004–2022. Blue line area within the NE South America highlights (<b>a</b>) its semi-arid domain and (<b>b</b>) its Caatinga biome, respectively.</p>
Full article ">Figure 11
<p>Monthly averaged SEVIRI NDVI anomalies for each NDVIa pattern over the 2004–2022 period. The N1 to N8 label located in the top-center in each panel indicates the NDVI cluster. The F1 to F8 label located in the bottom-right in each panel indicates the decreasing order of each NDVI cluster according to the median calculated from spatially averaged NDVIa over the entire northeastern South American region.</p>
Full article ">Figure 12
<p>As in <a href="#atmosphere-14-01629-f011" class="html-fig">Figure 11</a>, but here for the SPI3 patterns. The S1 to S8 label located in the top-center in each panel indicates the SPI3 cluster.</p>
Full article ">Figure A1
<p>It shows the slope of the significant linear trends of (<b>a</b>) visible (0.64 µm) and (<b>b</b>) infrared thermal (10.8 µm) radiances from the pentad-SEVIRI spectral images over the 2004–2022 period. Green slopes show a significant positive trend in visible radiance (0.64 µm) coupled to an increase in infrared thermal radiance (10.8 µm). Red slopes show a coupling between a decrease in both radiances. For instance, the northeastern South American region shows a patchy result, with some areas showing a decrease in visible radiance, although the infrared thermal increases (orange areas in <a href="#atmosphere-14-01629-f0A1" class="html-fig">Figure A1</a>) are probably related to a reduction in clouds in the 18 years of the study.</p>
Full article ">
29 pages, 4450 KiB  
Article
Statistical Downscaling of SEVIRI Land Surface Temperature to WRF Near-Surface Air Temperature Using a Deep Learning Model
by Afshin Afshari, Julian Vogel and Ganesh Chockalingam
Remote Sens. 2023, 15(18), 4447; https://doi.org/10.3390/rs15184447 - 9 Sep 2023
Cited by 2 | Viewed by 2631
Abstract
The analysis of the near-surface air temperature is vital for many applications such as urban heat islands and climate change studies. In particular, extreme weather events are typically localized and so should the corresponding adaptation measures. However, climate scientists are often confronted with [...] Read more.
The analysis of the near-surface air temperature is vital for many applications such as urban heat islands and climate change studies. In particular, extreme weather events are typically localized and so should the corresponding adaptation measures. However, climate scientists are often confronted with the difficulty of providing reliable predictions at high spatial resolutions in the order of 1 km. We propose to train a convolutional neural network model to emulate the hourly high-resolution near-surface air temperature field simulated by the Weather Research and Forecasting (WRF) software over a period of 18 months. The model is driven by current and past lags of coarse SEVIRI land surface temperature fields. This nowcasting application implements a downscaling of the spatial resolution of the input by about a factor of four, while establishing a correlation between current and past land surface temperature maps and the current near-surface air temperature field. The U-Net variant that is proposed in this study uses regularization to prevent over-fitting and implements a novel space-time approach, where multiple time steps are fed into the model through 3D convolution layers. Besides LST, the model also uses urban density as additional static input to be able to predict the temperature more accurately in urban areas and to improve the generalizability of the trained model. The performance of our U-Net model is assessed via comparison to an MLR benchmark (ridge regularization). The model is shown to be superior on all performance metrics. It achieves a mean absolute error of 1.36 °C versus 1.49 °C for benchmark (a 9% relative improvement) and a root mean square error of 1.77 °C versus 1.91 °C for benchmark (a 7% relative improvement). The absolute error of the model is less than 2 °C for 77% of the prediction points versus 72% for the benchmark (a 7% relative improvement). The improvement over the benchmark is even more significant during extreme hot periods. We demonstrate the generalizability of the approach by testing the trained model on unseen spatial domains. Full article
(This article belongs to the Section Urban Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Land–use categories of the finest domain d03 of the WRF mesoscale model setup, where only occurring categories are shown. Also shown are positions and labels of DWD measurement stations used for model validation. The primary region of interest (ROI) was used for training the statistical downscaling models and the secondary ROI was used for additional validation of the deep learning model on mostly unseen data.</p>
Full article ">Figure 2
<p>Data splitting into training, validation, and testing sets. The data were split into 18 equal parts, where each part was approximately one month in length. The splits into different testing periods are denoted by index <span class="html-italic">j</span> and the split into different training/validation periods is denoted by index <span class="html-italic">i</span>. There are a total of 63 different split sequences (<math display="inline"><semantics> <mrow> <mi>J</mi> <mo>=</mo> <mn>9</mn> </mrow> </semantics></math> testing splits and <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>=</mo> <mn>7</mn> </mrow> </semantics></math> training/validation splits for each testing split).</p>
Full article ">Figure 3
<p>Deep learning downscaling model: U-Net variant “3D-space-time-residual”.</p>
Full article ">Figure 4
<p>Time-averaged model biases (98th percentile) for both the linear regression and the deep learning model.</p>
Full article ">Figure 5
<p>Time-averaged model inputs, outputs, and predictions for the deep learning model for different seasons.</p>
Full article ">Figure 6
<p>Box plots indicating the spatial variation of the error metrics within the 48 × 48 pixels of the core domain averaged for different seasons—full year, winter, spring, summer, and fall.</p>
Full article ">Figure 7
<p>Spatially-averaged time-series results for the deep learning model for different seasons. On the left-hand side, the daily average temperatures are shown for each day of the season, and on the right-hand side, the averaged temperature profile is shown for each hour of the day.</p>
Full article ">
6 pages, 2929 KiB  
Proceeding Paper
A Satellite-Based Evaluation of Upper-Level Aviation Turbulence Events over Europe during November 2009: A Case Study
by Vasileios T. Gerogiannis and Haralambos Feidas
Environ. Sci. Proc. 2023, 26(1), 61; https://doi.org/10.3390/environsciproc2023026061 - 25 Aug 2023
Viewed by 748
Abstract
Aviation turbulence is a major concern for flight safety. Detecting and nowcasting upper-level turbulence is usually associated with known sources of turbulence, such as convective clouds and transverse cirrus bands. However, in extended clear-air conditions where no optical indicators are present, this can [...] Read more.
Aviation turbulence is a major concern for flight safety. Detecting and nowcasting upper-level turbulence is usually associated with known sources of turbulence, such as convective clouds and transverse cirrus bands. However, in extended clear-air conditions where no optical indicators are present, this can be challenging for both aviation forecasters and pilots. This study aims to evaluate heavy–severe aviation scale turbulence events over 20.000 ft, by utilizing satellite data from MSG SEVIRI radiometer and in situ turbulence reports from en-route aircraft flights over Europe. We analyze 92 heavy–severe turbulence events during November 2009. The results could give an estimate of possible turbulence detection to pilots and aviation forecasters to identify and avoid upper-level turbulence, increasing flight safety. Full article
Show Figures

Figure 1

Figure 1
<p>The 92 ULT events over Europe in November 2009, as measured from aircraft over 20.000 ft.</p>
Full article ">Figure 2
<p>NCT H–S ULT event (circled red dot), recorded at 1430 UTC on 21 November 2009. The event was recorded over France, at FL226 (DEVG: 5.3 m/s) north of cumulonimbus clouds “C”, southeast of the baroclinic mature vortex “L”: (<b>a</b>) MSG RGB Natural Color, showing that the event is over clear-sky area but in the proximity of the convective clouds; (<b>b</b>) MSG RGB Day Microphysics, showing the isolated convective clouds over France.</p>
Full article ">Figure 3
<p>MSG water vapor image (WV6.2), at 2015 UTC on 24 November 2009. The H–S upper-level CAT event (circled red dot) was recorded over the Balkans, at FL232, (DEVG: 5.3 m/s) in the tropopause folding “F” at the outflow of a comma cloud.</p>
Full article ">Figure 4
<p>Transverse cirrus bands “T” extent radially away from the convection clouds “C”, almost perpendicular to the jet stream axis-associated H–S ULT event (circled red dot), recorded at 1900 UTC on 23 November 2009. The event was recorded over France, at FL230 (DEVG: 4.8 m/s): (<b>a</b>) MSG RGB Airmass; (<b>b</b>) MSG WV6.2 μm.</p>
Full article ">Figure 5
<p>H–S ULT related to embedded cumulonimbus clous “E” (circled red dot), over a baroclinic conveyor belt, recorded over France at FL226 at 1730 UTC on 25 November 2009 (DEVG: 4.8 m/s): (<b>a</b>) MSG RGB Airmass; (<b>b</b>) MSG WV6.2 μm.</p>
Full article ">
7 pages, 2281 KiB  
Proceeding Paper
Utilizing HSAF SE-E-SEVIRI (H10) Product to Track Seasonal and Monthly Variability in Snow Cover in Part of Southern Balkans
by Alexandros Paraskevas, Vasileios Skarmintzos, Theodoros Stamatopoulos, Antonios Polyzos, Ioannis Matsangouras, Konstantinos Kasapas and Panagiotis Nastos
Environ. Sci. Proc. 2023, 26(1), 55; https://doi.org/10.3390/environsciproc2023026055 - 25 Aug 2023
Viewed by 639
Abstract
Understanding the extent and location of snow cover is crucial for studying climate, hydrology, and ecosystems. Remote sensing, such as EUMETSAT’s H-SAF project, enables snow-cover monitoring using SEVIRI satellite data. This study analyzes the spatial variability of snow cover in the southern Balkans [...] Read more.
Understanding the extent and location of snow cover is crucial for studying climate, hydrology, and ecosystems. Remote sensing, such as EUMETSAT’s H-SAF project, enables snow-cover monitoring using SEVIRI satellite data. This study analyzes the spatial variability of snow cover in the southern Balkans over the past five years, focusing on the H10 product. Results reveal distinct seasonal dynamics, with limited snow cover during autumn and extensive cover during winter. January and February exhibit the highest snowfall probability and better snow cover persistence compared to autumn. Winter exhibits greater extent, duration, and quality of snow cover (>80% in December and January, ~70% in February). These findings contribute to regional climate understanding and water resource management. Full article
Show Figures

Figure 1

Figure 1
<p>H10 product consists of four different classes: snow (white), cloud (cyan), water (blue), bare ground (green), and no data (black) (source EUMETSAT HSAF).</p>
Full article ">Figure 2
<p>The area of study (black box) over the southern Balkans, 41.3° N to 39.80° N latitude and 19.50° E to 26.00° E longitude.</p>
Full article ">Figure 3
<p>Average number of days of snow cover for autumn of 2018–2022.</p>
Full article ">Figure 4
<p>Average number of days of snow cover for winter of 2018–2022.</p>
Full article ">
31 pages, 14168 KiB  
Article
Towards the Accurate Automatic Detection of Mesoscale Convective Systems in Remote Sensing Data: From Data Mining to Deep Learning Models and Their Applications
by Mikhail Krinitskiy, Alexander Sprygin, Svyatoslav Elizarov, Alexandra Narizhnaya, Andrei Shikhov and Alexander Chernokulsky
Remote Sens. 2023, 15(14), 3493; https://doi.org/10.3390/rs15143493 - 11 Jul 2023
Cited by 3 | Viewed by 2054
Abstract
Mesoscale convective systems (MCSs) and associated hazardous meteorological phenomena cause considerable economic damage and even loss of lives in the mid-latitudes. The mechanisms behind the formation and intensification of MCSs are still not well understood due to limited observational data and inaccurate climate [...] Read more.
Mesoscale convective systems (MCSs) and associated hazardous meteorological phenomena cause considerable economic damage and even loss of lives in the mid-latitudes. The mechanisms behind the formation and intensification of MCSs are still not well understood due to limited observational data and inaccurate climate models. Improving the prediction and understanding of MCSs is a high-priority area in hydrometeorology. One may study MCSs either employing high-resolution atmospheric modeling or through the analysis of remote sensing images which are known to reflect some of the characteristics of MCSs, including high temperature gradients of cloud-top, specific spatial shapes of temperature patterns, etc. However, research on MCSs using remote sensing data is limited by inadequate (in size) databases of satellite-identified MCSs and poorly equipped automated tools for MCS identification and tracking. In this study, we present (a) the GeoAnnotateAssisted tool for fast and convenient visual identification of MCSs in satellite imagery, which is capable of providing AI-generated suggestions of MCS labels; (b) the Dataset of Mesoscale Convective Systems over the European Territory of Russia (DaMesCoS-ETR), which we created using this tool, and (c) the Deep Convolutional Neural Network for the Identification of Mesoscale Convective Systems (MesCoSNet), constructed following the RetinaNet architecture, which is capable of identifying MCSs in Meteosat MSG/SEVIRI data. We demonstrate that our neural network, optimized in terms of its hyperparameters, provides high MCS identification quality (mAP=0.75, true positive rate TPR=0.61) and a well-specified detection uncertainty (false alarm ratio FAR=0.36). Additionally, we demonstrate potential applications of the GeoAnnotateAssisted labelling tool, the DaMesCoS-ETR dataset, and the MesCoSNet neural network in addressing MCS research challenges. Specifically, we present the climatology of axisymmetric MCSs over the European territory of Russia from 2014 to 2020 during summer seasons (May to September), obtained using MesCoSNet with Meteosat MSG/SEVIRI data. The automated identification of MCSs by the MesCoSNet artificial neural network opens up new avenues for previously unattainable MCS research topics. 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>An example of four mesoscale convective systems labeled by our expert in Meteosat imagery using our GeoAnnotateAssisted labeling and tracking tool. Here, the actual representations of transformed channels <math display="inline"><semantics><mrow><mi>c</mi><mi>h</mi><msub><mn>9</mn><mi>n</mi></msub><mo>,</mo><mo> </mo><mi>c</mi><mi>h</mi><msub><mn>5</mn><mi>n</mi></msub></mrow></semantics></math> and <math display="inline"><semantics><mover accent="true"><mi>b</mi><mo>˜</mo></mover></semantics></math> are shown (see <a href="#sec3dot2dot2-remotesensing-15-03493" class="html-sec">Section 3.2.2</a> for details) using the actual color maps employed in GeoAnnotateAssisted. Green ellipses are the labels an expert has placed as a result of the examination of these representations. Here, in panels, the following channels are shown: (<b>a</b>) <math display="inline"><semantics><mrow><mi>c</mi><mi>h</mi><msub><mn>9</mn><mi>n</mi></msub></mrow></semantics></math>, (<b>b</b>) <math display="inline"><semantics><mrow><mi>c</mi><mi>h</mi><msub><mn>5</mn><mi>n</mi></msub></mrow></semantics></math> and (<b>c</b>) <math display="inline"><semantics><mover accent="true"><mi>b</mi><mo>˜</mo></mover></semantics></math>. We do not show temperature color bars, since the values of the re-scaled features <math display="inline"><semantics><mrow><mi>c</mi><mi>h</mi><msub><mn>9</mn><mi>n</mi></msub></mrow></semantics></math>, <math display="inline"><semantics><mrow><mi>c</mi><mi>h</mi><msub><mn>5</mn><mi>n</mi></msub></mrow></semantics></math> and <math display="inline"><semantics><mover accent="true"><mi>b</mi><mo>˜</mo></mover></semantics></math> are unitless.</p>
Full article ">Figure 2
<p>High-level architecture of the annotation tool GeoAnnotateAssisted. The pictogram in the top-left corner stands for a user performing the annotation of MCS.</p>
Full article ">Figure 3
<p>Structure of the database used in GeoAnnotateAssisted to store the labels and tracking information. It matches completely with the structure of DaMesCoS-ETR. We provide a complete description of this schema in the GitHub repository of DaMesCoS (<a href="http://github.com/mkrinitskiy/damescos" target="_blank">http://github.com/mkrinitskiy/damescos</a>, accessed on 20 April 2023). Here, label “1” close to certain ID fields indicate that the corresponding identifiers are unique within the scope of its table (e.g., the field “id” of table “labels”) Also, label “*” close to “track_id” field of the table “track_labels” indicates that this identifier is not unique within the scope of the table.</p>
Full article ">Figure 4
<p>Lifecycle characteristics of MCS in DaMesCoS-ETR based on expert labeling: (<b>a</b>) the distribution of path lengths <span class="html-italic">L</span>, km, (<b>b</b>) the distribution of lifetime <math display="inline"><semantics><msub><mi>T</mi><mrow><mi>l</mi><mi>i</mi><mi>f</mi><mi>e</mi></mrow></msub></semantics></math>, h, and (<b>c</b>) the distribution of propagation velocity <span class="html-italic">V</span>, km/h. The distributions are presented in logarithmic scale of the variables. The lines are the kernel density estimation of the distributions.</p>
Full article ">Figure 5
<p>Focal loss <math display="inline"><semantics><msub><mi mathvariant="script">L</mi><mrow><mi>F</mi><mi>L</mi></mrow></msub></semantics></math> graphs vs. probability of ground truth class <math display="inline"><semantics><msub><mi>p</mi><mi>t</mi></msub></semantics></math> depending on focusing parameter <math display="inline"><semantics><mi>γ</mi></semantics></math> (from the original paper [<a href="#B63-remotesensing-15-03493" class="html-bibr">63</a>]). Here, we modified the captions for the sake of figure clarity. Note, that in the case of <math display="inline"><semantics><mrow><mi>γ</mi><mo>=</mo><mn>0</mn></mrow></semantics></math>, focal loss (FL) is equivalent to the cross-entropy (CE).</p>
Full article ">Figure 6
<p>RetinaNet architecture, from the original paper [<a href="#B63-remotesensing-15-03493" class="html-bibr">63</a>]. Here, we modified the captions of the first subnet to “backbone subnet” according to the common terminology established by 2023.</p>
Full article ">Figure 7
<p>Difference between knowledge extraction in traditional machine learning (<b>left</b>) and with transfer learning employed (<b>right</b>) (from Pan and Yang 2010 [<a href="#B110-remotesensing-15-03493" class="html-bibr">110</a>]).</p>
Full article ">Figure 8
<p>Stages of <math display="inline"><semantics><mrow><mi>B</mi><mi>T</mi><mi>D</mi></mrow></semantics></math> data transformations within domain adaptation: (<b>a</b>) the original <math display="inline"><semantics><mrow><mi>B</mi><mi>T</mi><mi>D</mi></mrow></semantics></math> distribution; (<b>b</b>) its normalized distribution acquired as a result of transformation in Equation (<a href="#FD4-remotesensing-15-03493" class="html-disp-formula">4</a>); (<b>c</b>) plot of the transfer function presented in Equation (<a href="#FD5-remotesensing-15-03493" class="html-disp-formula">5</a>), and (<b>d</b>) the empirical distribution of non-linearly scaled <math display="inline"><semantics><mrow><mi>B</mi><mi>T</mi><mi>D</mi></mrow></semantics></math> employed as one of the spatially distributed features in our study.</p>
Full article ">Figure 9
<p>Distributions of the features of remote sensing data transformed in accordance with Equations (<a href="#FD3-remotesensing-15-03493" class="html-disp-formula">3</a>) and (<a href="#FD5-remotesensing-15-03493" class="html-disp-formula">5</a>).</p>
Full article ">Figure 10
<p>Sampling strategy applied in our study for splitting the dataset into training and testing subsets. Here, we display individual Meteosat snapshots with corresponding MCS labels as “measurements” (thin vertical rectangles); days of observations are shown by thick rectangles; days of labeled observations contributing to a train subset are shown in light blue; and days of labeled observations contributing to a test subset are shown in light green.</p>
Full article ">Figure 11
<p>Learning rate schedule employed in our study.</p>
Full article ">Figure 12
<p>Plots of metrics <math display="inline"><semantics><mrow><mi>I</mi><mi>o</mi><mi>U</mi></mrow></semantics></math> and AP depending on <math display="inline"><semantics><msub><mi>t</mi><mi>p</mi></msub></semantics></math> as the supporting diagram for choosing the <math display="inline"><semantics><msub><mi>t</mi><mi>p</mi></msub></semantics></math>.</p>
Full article ">Figure 13
<p>Examples of MesCoSNet application in MCS identification problem using Meteosat remote sensing data. Here, source satellite imagery bands are presented using the color maps delivering the same visual experience as color maps exploited in GeoAnnotateAssisted. MCS labels detected by MesCoSNet are shown as pink rectangles; ground truth MCS labels from DaMesCoS-ETR are shown as yellow rectangles. We also display the certainty rate of MesCoSNet (class MCS probability) for the detected labels (see pink text in the top-left corner of pick rectangles.). Here, we present the cases of (<b>a</b>) our neural network MesCoSNet identifying an MCS close to ground truth and (<b>b</b>) our neural network MesCoSNet identifying one of three ground truth MCSs, thus, two MCSs were missed, whereas the identified one is located close to ground truth.</p>
Full article ">Figure 14
<p>Frequency of MCS occurrence over ETR in summer (May to September) as a result of MesCoSNet application with Meteosat SEVIRI imagery from 2014 till 2017: (<b>left</b>) frequency map; (<b>right</b>) the averaged over pink rectangle diurnal variation of MCS occurrence probability.</p>
Full article ">
25 pages, 11581 KiB  
Article
Machine Learning for Fog-and-Low-Stratus Nowcasting from Meteosat SEVIRI Satellite Images
by Driss Bari, Nabila Lasri, Rania Souri and Redouane Lguensat
Atmosphere 2023, 14(6), 953; https://doi.org/10.3390/atmos14060953 - 30 May 2023
Cited by 5 | Viewed by 2406
Abstract
Fog and low stratus (FLS) are meteorological phenomena that have a significant impact on all ways of transportation and public safety. Due to their similarity, they are often grouped together as a single category when viewed from a satellite perspective. The early detection [...] Read more.
Fog and low stratus (FLS) are meteorological phenomena that have a significant impact on all ways of transportation and public safety. Due to their similarity, they are often grouped together as a single category when viewed from a satellite perspective. The early detection of these phenomena is crucial to reduce the negative effects that they can cause. This paper presents an image-based approach for the short-term nighttime forecasting of FLS during the next 5 h over Morocco, based on geostationary satellite observations (Meteosat SEVIRI). To achieve this, a dataset of hourly night microphysics RGB product was generated from native files covering the nighttime cold season (October to April) of the 5-year period (2016–2020). Two optical flow techniques (sparse and dense) and three deep learning techniques (CNN, Unet and ConvLSTM) were used, and the performance of the developed models was assessed using mean squared error (MSE) and structural similarity index measure (SSIM) metrics. Hourly observations from Meteorological Aviation Routine Weather Reports (METAR) over Morocco were used to qualitatively compare the FLS existence in METAR, where it is also shown by the RGB product. Results analysis show that deep learning techniques outperform the traditional optical flow method with SSIM and MSE of about 0.6 and 0.3, respectively. Deep learning techniques show promising results during the first three hours. However, their performance is highly dependent on the number of filters and the computing resources, while sparse optical flow is found to be very sensitive to mask definition on the target phenomenon. Full article
(This article belongs to the Special Issue Decision Support System for Fog)
Show Figures

Figure 1

Figure 1
<p>Map showing the geographical position of Morocco in the northwestern corner of Africa and also the position of the used synoptic stations over the study domain. The topography is also shown on this map.</p>
Full article ">Figure 2
<p>Heatmap showing the frequency distribution of the onset time of FLS events as a function of their duration over the 17 synoptic stations that are working 24 h/24 h from a total of 43 stations. The period covers the cold season (October to April) of 5 years (2016–2020).</p>
Full article ">Figure 3
<p>Diagram of applied methodology for deep learning model development.</p>
Full article ">Figure 4
<p>MSE (<b>left</b>) and SSIM (<b>right</b>) distribution of sparse optical flow with and without mask definition as a function of the forecast hour for the case study of 21 November 2020. Persistence scores are also drawn as the benchmark.</p>
Full article ">Figure 5
<p>SSIM distribution of CNN, Unet and ConvLSTM with 64, 128 and 192 filters as a function of forecast hour for the case study of 18–19 February 2020. Persistence scores are also drawn as benchmarks.</p>
Full article ">Figure 6
<p>Experiment setup: schematic description of the four configurations used to predict the spatio-temporal evolution of the FLS event that occurred during the night 18–19 February 2020.</p>
Full article ">Figure 7
<p>Nowcasting output from the sparse optical flow model for the four configurations vs. observed satellite images. From top to bottom per two lines (forecast vs. observation): configuration 1, 2, 3 and 4. Case study of 18–19 February 2020. Red circles refer to the tracked fog areas.</p>
Full article ">Figure 8
<p>Nowcasting output from the dense optical flow model for the four configurations vs. observed satellite images. From top to bottom per two lines (forecast vs. observation): configuration 1, 2, 3 and 4. Case study of 18–19 February 2020. Red circles refer to the tracked fog areas.</p>
Full article ">Figure 9
<p>MSE (<b>left</b>) and SSIM (<b>right</b>) distribution of CNN, ConvLSTM, Unet, SparseOF and DenseOF as a function of forecast hour for the case study of 18–19 February 2020. Persistence scores are also drawn as benchmark.</p>
Full article ">Figure 10
<p>From top to down: (row 1) the input images, (row 2) the ground truth, (row 3) the CNN model predictions, (row 4) the ConvLSTM model predictions, and (row 5) the forecasts of Unet model. Use case of 18–19 February 2020.</p>
Full article ">Figure 11
<p>From top to down: (row 1) the input images, (row 2) the ground truth, (row 3) the CNN model predictions, (row 4) the ConvLSTM model predictions, and (row 5) the forecasts of Unet model. Use case of 18–19 February 2020.</p>
Full article ">Figure 12
<p>From top to down: (row 1) the input images, (row 2) the ground truth, (row 3) the CNN model predictions, (row 4) the ConvLSTM model predictions, and (row 5) the forecasts of Unet model. Use case of 18–19 February 2020.</p>
Full article ">Figure 13
<p>From top to down: (row 1) the input images, (row 2) the ground truth, (row 3) the CNN model predictions, (row 4) the ConvLSTM model predictions, and (row 5) the forecasts of Unet model. Use case of 18–19 February 2020.</p>
Full article ">Figure 14
<p>MSE (<b>left</b>) and SSIM (<b>right</b>) distribution as function of forecast hour for CNN, ConvLSTM, Unet, SparseOF and DenseOF over the testing set. Persistence scores are also drawn as benchmarks.</p>
Full article ">
Back to TopTop