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

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

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (418)

Search Parameters:
Keywords = MODIS-Aqua

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
24 pages, 6494 KiB  
Article
Reconstruction of Fine-Spatial-Resolution FY-3D-Based Vegetation Indices to Achieve Farmland-Scale Winter Wheat Yield Estimation via Fusion with Sentinel-2 Data
by Xijia Zhou, Tao Wang, Wei Zheng, Mingwei Zhang and Yuanyuan Wang
Remote Sens. 2024, 16(22), 4143; https://doi.org/10.3390/rs16224143 - 6 Nov 2024
Viewed by 565
Abstract
The spatial resolution (250–1000 m) of the FY-3D MERSI is too coarse for agricultural monitoring at the farmland scale (20–30 m). To achieve the winter wheat yield (WWY) at the farmland scale, based on FY-3D, a method framework is developed in this work. [...] Read more.
The spatial resolution (250–1000 m) of the FY-3D MERSI is too coarse for agricultural monitoring at the farmland scale (20–30 m). To achieve the winter wheat yield (WWY) at the farmland scale, based on FY-3D, a method framework is developed in this work. The enhanced deep convolutional spatiotemporal fusion network (EDCSTFN) was used to perform a spatiotemporal fusion on the 10 day interval FY-3D and Sentinel-2 vegetation indices (VIs), which were compared with the enhanced spatial and temporal adaptive reflectance fusion model (ESTARFM). In addition, a BP neural network was built to calculate the farmland-scale WWY based on the fused VIs, and the Aqua MODIS gross primary productivity product was used as ancillary data for WWY estimation. The results reveal that both the EDCSTFN and ESTARFM achieve satisfactory precision in the fusion of the Sentinel-2 and FY-3D VIs; however, when the period of spatiotemporal data fusion is relatively long, the EDCSTFN can achieve greater precision than ESTARFM. Finally, the WWY estimation results based on the fused VIs show remarkable correlations with the WWY data at the county scale and provide abundant spatial distribution details about the WWY, displaying great potential for accurate farmland-scale WWY estimations based on reconstructed fine-spatial-temporal-resolution FY-3D data. Full article
Show Figures

Figure 1

Figure 1
<p>Overview of the study region: (<b>a</b>) location of the Weihe Plain; (<b>b</b>) FY-3D false colour composite image for 3 May 2020; and (<b>c</b>) locations of the county-scale WWY data points used in the WWY estimation.</p>
Full article ">Figure 2
<p>Flowchart of the 10 day interval VI imagery reconstruction and farmland-scale WWY estimation.</p>
Full article ">Figure 3
<p>Flowchart of farmland-scale WWY estimation: (<b>a</b>) Y estimation model based on the cumulative GPP; and (<b>b</b>) farmland-scale Y estimation model based on multiple parameters.</p>
Full article ">Figure 4
<p>Results of the consistency analysis of the Sentinel-2 and FY-3D VIs: (<b>a</b>) R<sup>2</sup> values between the aggregated Sentinel-2 VI imagery and the FY-3D VI imagery at an SR of 250 m; and (<b>b</b>) average deviations and RMSE values of the fitting results between the aggregated Sentinel-2 VI imagery and FY-3D VI imagery. The error line in (<b>b</b>) denotes the RMSE of the fitting results.</p>
Full article ">Figure 5
<p>EVI at each WW growing stage from 2020 to 2022.</p>
Full article ">Figure 6
<p>Y estimation model based on the cumulative GPP for the main WW growing period and the Y estimation precision evaluation results in 2020 and 2021. The dotted lines in the figures denote the fitted linear functions, which are close to the diagonal solid lines, indicating that the systematic deviation in the Y estimation results is small. (<b>a</b>) Linear regression model between the cumulative GPP data for the main WW growing period and the county-scale WWY from 2014 to 2018, (<b>b</b>) linear regression results between the WWY estimation results from 2020 based on the cumulative GPP and county-scale Y statistical data, and (<b>c</b>) linear regression results between the WWY estimation results in 2021 based on the cumulative GPP and county-scale Y statistical data.</p>
Full article ">Figure 7
<p>WWY estimation results for 2020 to 2022 based on the MODIS cumulative GPP data.</p>
Full article ">Figure 8
<p>Farmland-scale WWY estimation results for the Weihe Plain from 2020 to 2022 based on multiple parameters.</p>
Full article ">Figure 9
<p>Linear regression results between the farmland-scale WWY estimation results and the Y statistical data in 2020 and 2021. The dotted lines in the figures denote the fitted linear functions, which are close to the diagonal solid lines, indicating that the systematic deviation of the Y estimation results is small.</p>
Full article ">
28 pages, 14472 KiB  
Article
Characteristics of R2019 Processing of MODIS Sea Surface Temperature at High Latitudes
by Chong Jia, Peter J. Minnett and Malgorzata Szczodrak
Remote Sens. 2024, 16(21), 4102; https://doi.org/10.3390/rs16214102 - 2 Nov 2024
Viewed by 436
Abstract
Satellite remote sensing is the best way to derive sea surface skin temperature (SSTskin) in the Arctic. However, as surface temperature retrieval algorithms in the infrared (IR) part of the electromagnetic spectrum are designed to compensate for atmospheric effects mainly due [...] Read more.
Satellite remote sensing is the best way to derive sea surface skin temperature (SSTskin) in the Arctic. However, as surface temperature retrieval algorithms in the infrared (IR) part of the electromagnetic spectrum are designed to compensate for atmospheric effects mainly due to water vapor, MODIS SSTskin retrievals have larger uncertainties at high latitudes where the atmosphere is very dry and cold, which is an extreme in the distribution of global conditions. MODIS R2019 SSTskin fields are currently derived using latitudinally and monthly dependent algorithm coefficients, including an additional band above 60°N to better represent the effects of Arctic atmospheres. However, the R2019 processing of MODIS SSTskin still has some unrevealed error characteristics. This study uses 21 years (2002–2022) of collocated, simultaneous satellite brightness temperature (BT) data from Aqua MODIS and in situ buoy-measured subsurface temperature data from iQuam for validation. Unlike elsewhere over the oceans, the 11 μm and 12 μm BT differences are poorly related to the column water vapor at high latitudes, resulting in poor atmospheric water vapor correction. Anomalous BT difference signals are identified, caused by the temperature and humidity inversions in the lower troposphere, which are especially significant during the summer. Although the existence of negative BT differences is physically reasonable, this makes the retrieval algorithm lose its effectiveness. Moreover, the statistics of the MODIS SSTskin data when compared with the iQuam buoy temperature data show large differences (in terms of mean and standard deviation) for the matchups at the Northern Atlantic and Pacific sides of the Arctic due to the disparity of in situ measurements and distinct surface and vertical atmospheric conditions. Therefore, it is necessary to further improve the retrieval algorithms to obtain more accurate MODIS SSTskin data to study surface ocean processes and climate change in the Arctic. Full article
(This article belongs to the Section Ocean Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Distribution of global drifters, moored buoys, and some other platforms in June 2024 from <a href="https://www.ocean-ops.org/dbcp/network/maps.html" target="_blank">https://www.ocean-ops.org/dbcp/network/maps.html</a> (accessed on 2 September 2024).</p>
Full article ">Figure 2
<p>(<b>a</b>) Time series of the Aqua MODIS-iQuam SST difference (ΔSST) above 60°N from 2002 to 2022, including both QL = 0 and QL = 1 data. The black dashed horizontal line indicates the value of −0.17 K. (<b>b</b>) Histogram of the ΔSST with fitted normal distribution as blue dashed curve. (<b>c</b>) Scatter plot of Aqua MODIS SST<sub>skin</sub> and in situ buoy SST measurements from iQuam colored according to data density. (<b>d</b>) Map of Aqua MODIS ΔSST. (<b>e</b>) Map of Aqua MODIS ΔSST &gt; 2 K or &lt;−2 K.</p>
Full article ">Figure 3
<p>Aqua MODIS ΔSST as a function of time difference between the satellite and in situ buoy measurements for the (<b>a</b>) daytime and (<b>b</b>) nighttime data. The scatter plots are colored according to the data density. The dashed black horizontal lines are at −0.17 K, and the linear regressions are plotted as black solid lines with the functional expressions given in the plots.</p>
Full article ">Figure 4
<p>Similar to <a href="#remotesensing-16-04102-f003" class="html-fig">Figure 3</a>, but for the Aqua MODIS ΔSST as a function of spatial distance between the satellite and in situ buoy measurements.</p>
Full article ">Figure 5
<p>BT difference (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="normal">B</mi> <mi mathvariant="normal">T</mi> </mrow> <mrow> <mn>11</mn> </mrow> </msub> <mo>−</mo> <msub> <mrow> <mi mathvariant="normal">B</mi> <mi mathvariant="normal">T</mi> </mrow> <mrow> <mn>12</mn> </mrow> </msub> </mrow> </semantics></math>) as a function of total column water vapor in the MUDB during the (<b>a</b>) day and (<b>b</b>) at night. The scatter plots are colored according to data density, and the black dashed lines indicate a BT difference of zero. The correlation coefficients are also given on the upper right.</p>
Full article ">Figure 6
<p>(<b>a</b>–<b>l</b>) Monthly scatter plots of Aqua MODIS ΔSST as a function of BT difference colored according to data density and fitted by linear regression in black dashed lines. (<b>m</b>) Monthly regression slope (blue) and intercept (red) variation, with the error bars indicating the uncertainties.</p>
Full article ">Figure 7
<p>(<b>a</b>) Map distribution of the anomalous BT differences in the R2019 MUDB above 60°N. Histograms of (<b>b</b>) latitude and (<b>c</b>) longitude for Aqua MODIS measurements in the MUDB when the BT differences are &gt;0.036 K (light blue) and are ≤0.036 K (light red).</p>
Full article ">Figure 8
<p>MERRA-2 reanalysis atmospheric profiles matched with Aqua MODIS R2019 MUDB north of 60°N including (<b>a</b>) air temperature and (<b>b</b>) specific humidity, plotted as the mean (dashed line) ± 1 STD (envelope) beneath the 500 hPa level for the measurements of BT difference &gt; 0.036 K (blue) and BT difference ≤ 0.036 K (red).</p>
Full article ">Figure 9
<p>(<b>a</b>–<b>l</b>) Monthly distributions (mean ± 1 STD) of the MERRA-2 air temperature profiles beneath the 500 hPa level for the normal and abnormal BT difference cases. The number of data points in each dataset is given in the lower left section of each panel.</p>
Full article ">Figure 10
<p>Similar to <a href="#remotesensing-16-04102-f009" class="html-fig">Figure 9</a>a–1, but for the MERRA-2 specific humidity profiles beneath the 500 hPa level. The number of each dataset is given in the upper right section of each panel.</p>
Full article ">Figure 11
<p>Monthly occurrence frequency of temperature inversion (TI) only (blue), humidity inversion (HI) only (red), simultaneous TI and HI (yellow), and neither TI nor HI (green) for (<b>a</b>) BT difference &gt; 0.036 K and (<b>b</b>) BT difference ≤ 0.036 K. The highest frequency among the TI/HI situations in each month is shown in the corresponding part of the column.</p>
Full article ">Figure 12
<p>Monthly variations of (<b>a</b>) TI strength, (<b>b</b>) HI strength, and (<b>c</b>) BT difference between 11 μm and 12 μm for BT difference &gt; 0.036 K (blue) and BT difference ≤ 0.036 K (red) only when TI or HI exists. For the boxplots, the ends of the boxes, the ends of the whiskers, and the line in the box represent the 25th and 75th percentiles, the minimum and maximum values that are not outliers, and the median, respectively. Outliers are values beyond 75th percentile +1.5* (interquartile range) or less than 25th percentile −1.5* (interquartile range), which are not shown here. (<b>d</b>) Scatter plot of BT difference as a function of TI strength colored according to the HI strength. The black dotted line indicates the BT difference value of 0.036 K, the threshold of positive BT differences greater than the NEΔT effects. The red dots and error bars are the mean and STD of the BT difference, calculated at 1 K intervals of TI strength.</p>
Full article ">Figure 13
<p>Monthly boxplots of the Aqua MODIS ΔSST in the Atlantic Sector (blue) and the Pacific Sector (red) of the Arctic. Outliers are not shown.</p>
Full article ">Figure 14
<p>(<b>a</b>–<b>l</b>) Monthly distributions (mean ± 1 STD) of the MERRA-2 air temperature profiles beneath the 500 hPa level matched with the data in the MUDB at the Atlantic and Pacific Sectors of the Arctic. The number of each dataset is given in the lower left section of each panel.</p>
Full article ">Figure 15
<p>Similar to <a href="#remotesensing-16-04102-f014" class="html-fig">Figure 14</a>a–l, but for the MERRA-2 specific humidity profiles beneath the 500 hPa level. The number of each dataset is given in the upper right section of each panel.</p>
Full article ">Figure 16
<p>Monthly occurrence frequency of TI only (blue), HI only (red), simultaneous TI and HI (yellow), and neither TI nor HI (green) for the data in MUDB for the (<b>a</b>) Atlantic Sector and (<b>b</b>) Pacific Sector of the Arctic. The highest frequency among the inversion situations in each month is shown in the corresponding part of the column. Monthly variations of (<b>c</b>) TI strength and (<b>d</b>) HI strength for the Atlantic Sector (blue) and Pacific Sector (red) are plotted only when TI or HI exists. Monthly boxplots of (<b>e</b>) BT difference between 11 μm and 12 μm and (<b>f</b>) total column water vapor at the Atlantic side (blue) and Pacific side (red) are also given. The black dotted line in (<b>e</b>) is the horizon of 0.036 K. Outliers are not shown in the boxplots.</p>
Full article ">Figure 17
<p>Maps of monthly MERRA-2 total column water vapor north of 60°N in (<b>a</b>) January, (<b>b</b>) April, (<b>c</b>) July, and (<b>d</b>) October. The data are averaged from 2002 to 2022.</p>
Full article ">Figure 18
<p>Monthly boxplots of (<b>a</b>) latitude of MODIS pixels and (<b>b</b>) in situ subsurface SST in the MUDB for the Atlantic (blue) and Pacific (red) Sectors. Outliers are not shown here.</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 988
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 ">
20 pages, 16133 KiB  
Article
Changes in Vegetation Cover and the Relationship with Surface Temperature in the Cananéia–Iguape Coastal System, São Paulo, Brazil
by Jakeline Baratto, Paulo Miguel de Bodas Terassi and Emerson Galvani
Remote Sens. 2024, 16(18), 3460; https://doi.org/10.3390/rs16183460 - 18 Sep 2024
Viewed by 779
Abstract
The objective of this article is to investigate the possible correlations between vegetation indices and surface temperature in the Cananéia–Iguape Coastal System (CICS), in São Paulo (Brazil). Vegetation index data from MODIS orbital products were used to carry out this work. The Normalized [...] Read more.
The objective of this article is to investigate the possible correlations between vegetation indices and surface temperature in the Cananéia–Iguape Coastal System (CICS), in São Paulo (Brazil). Vegetation index data from MODIS orbital products were used to carry out this work. The Normalized Difference Vegetation Index (NDVI) and the Enhanced Vegetation Index (EVI) were acquired from the MODIS/Aqua sensor (MYD13Q1) and the leaf area index (LAI) from the MODIS/Terra (MOD15A2H). Surface temperature data were acquired from MODIS/Aqua (MYD11A2). The data were processed using Google Earth Engine and Google Colab. The data were collected, and spatial and temporal correlations were applied. Correlations were applied in the annual and seasonal period. The annual temporal correlation between vegetation indices and surface temperature was positive, but statistically significant for the LAI, with r = 0.43 (90% significance). In the seasonal period, positive correlations occurred in JFM for all indices (95% significance). Spatially, the results of this research indicate that the largest area showed a positive correlation between VI and LST. The hottest and rainiest periods (OND and JFM) had clearer and more significant correlations. In some regions, significant and clear correlations were observed, such as in some areas in the north, south and close to the city of Iguape. This highlights the complexity of the interactions between vegetation indices and climatic attributes, and highlights the importance of considering other environmental variables and processes when interpreting changes in vegetation. However, this research has significantly progressed the field, by establishing new correlations and demonstrating the importance of considering climate variability, for a more accurate understanding of the impacts on vegetation indices. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the study area (<b>A</b>,<b>B</b>) and land use mapping (<b>C</b>).</p>
Full article ">Figure 2
<p>Variation in surface temperature and monthly (<b>A</b>) and annual (<b>B</b>) rainfall for the Cananéia-Iguape Coastal System for the 20032022 period. Source: MODIS/Aqua and CHIRPS, 2024.</p>
Full article ">Figure 3
<p>Annual variation in vegetation indices for the 2003–2022 period in the Cananéia–Iguape Coastal System.</p>
Full article ">Figure 4
<p>Scatter plot of annual NDVI (<b>a</b>), EVI (<b>b</b>) and LAI (<b>c</b>) values and surface temperature from 2003 to 2022.</p>
Full article ">Figure 5
<p>Scatter plot of seasonal values of VI–NDVI (<b>a</b>,<b>d</b>), EVI (<b>b</b>,<b>e</b>) and LAI (<b>c</b>,<b>f</b>)—and surface temperature for the JFM (<b>a</b>–<b>c</b>) and AMJ (<b>d</b>–<b>f</b>) quarter.</p>
Full article ">Figure 6
<p>Scatter plot of seasonal values of VI–NDVI (<b>a</b>,<b>d</b>), EVI (<b>b</b>,<b>e</b>) and LAI (<b>c</b>,<b>f</b>)—and climate variables for the JAS (<b>a</b>–<b>c</b>) and OND (<b>d</b>–<b>f</b>) quarter.</p>
Full article ">Figure 7
<p>Annual linear correlation between surface temperature and NDVI (<b>A</b>), EVI (<b>B</b>) and LAI (<b>C</b>) between 2003 and 2022.</p>
Full article ">Figure 8
<p>Seasonal linear correlation between surface temperature and VI between 2004 and 2022 for the JFM (<b>A</b>–<b>C</b>) and AMJ (<b>D</b>–<b>F</b>) periods.</p>
Full article ">Figure 9
<p>Seasonal linear correlation between surface temperature and VI between 2004 and 2022 for the JAS (<b>A</b>–<b>C</b>) and OND (<b>D</b>–<b>F</b>) periods.</p>
Full article ">
17 pages, 3602 KiB  
Article
Understanding Two Decades of Turbidity Dynamics in a Coral Triangle Hotspot: The Berau Coastal Shelf
by Faruq Khadami, Ayi Tarya, Ivonne Milichristi Radjawane, Totok Suprijo, Karina Aprilia Sujatmiko, Iwan Pramesti Anwar, Muhamad Faqih Hidayatullah and Muhamad Fauzan Rizky Adisty Erlangga
Water 2024, 16(16), 2300; https://doi.org/10.3390/w16162300 - 15 Aug 2024
Viewed by 946
Abstract
Turbidity serves as a crucial indicator of coastal water health and productivity. Twenty years of remote sensing data (2003–2022) from the Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) satellite were used to analyze the spatial and temporal variations in turbidity, as measured by total [...] Read more.
Turbidity serves as a crucial indicator of coastal water health and productivity. Twenty years of remote sensing data (2003–2022) from the Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) satellite were used to analyze the spatial and temporal variations in turbidity, as measured by total suspended matter (TSM), in the Berau Coastal Shelf (BCS), East Kalimantan, Indonesia. The BCS encompasses the estuary of the Berau River and is an integral part of the Coral Triangle, renowned for its rich marine and coastal habitats, including coral reefs, mangroves, and seagrasses. The aim of this research is to comprehend the seasonal and interannual patterns of turbidity and their associations with met-ocean parameters, such as wind, rainfall, and climate variations like the El Niño–Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD). The research findings indicate that the seasonal spatial pattern of turbidity is strongly influenced by monsoon winds, while its temporal patterns are closely related to river discharge and rainfall. The ENSO and IOD climate cycles exert an influence on the interannual turbidity variations, with turbidity values decreasing during La Niña and negative IOD events and conversely increasing during El Niño and positive IOD events. Furthermore, the elevated turbidity during negative IOD and La Niña coincides with rising temperatures, potentially acting as a compound stressor on marine habitats. These findings significantly enhance our understanding of turbidity dynamics in the BCS, thereby supporting the management of marine and coastal ecosystems in the face of changing climatic and environmental conditions. Full article
(This article belongs to the Section Oceans and Coastal Zones)
Show Figures

Figure 1

Figure 1
<p>Berau Coastal Shelf. The red boxes mark regions selected for spatial averaging, which encompass the north, river mouth, and south regions. The triangle symbol shows the river discharge station. The color shading in the image represents the varying bathymetry depths.</p>
Full article ">Figure 2
<p>Mean climatological TSM (mg/L) in the BCS.</p>
Full article ">Figure 3
<p>Monthly climatological TSM and wind rose.</p>
Full article ">Figure 4
<p>TSM variability at three locations in the BCS, showing median, interquartile ranges, max/min, and outliers. The black dots indicating outliers.</p>
Full article ">Figure 5
<p>Monthly differences in average TSM (mg/L): average TSM in the (<b>a</b>) north, (<b>b</b>) river mouth, and (<b>c</b>) south areas of the BCS from 2003 to 2022. The black dots indicating outliers.</p>
Full article ">Figure 6
<p>Power spectral density of average TSM in (<b>a</b>) north, (<b>b</b>) river mouth, and (<b>c</b>) south areas of the BCS.</p>
Full article ">Figure 7
<p>Berau River discharge data from May 2006 to January 2008. The black line represents the observed discharge values. The red line shows the data smoothed by a 30-day low-pass filter.</p>
Full article ">Figure 8
<p>The variation in the 15-month moving average of (<b>a</b>) the TSM at river mouth, (<b>b</b>) the precipitation, the (<b>c</b>) Dipole Mode Index, (<b>d</b>) El Niño–Southern Oscillation Index (NINO 3.4), and (<b>e</b>) Sea Surface Temperature (SST) anomaly. The red (blue) color indicating positive (negative) anomaly.</p>
Full article ">
19 pages, 9008 KiB  
Article
The Carpathian Agriculture in Poland in Relation to Other EU Countries, Ukraine and the Environmental Goals of the EU CAP 2023–2027
by Marek Zieliński, Artur Łopatka, Piotr Koza and Barbara Gołębiewska
Agriculture 2024, 14(8), 1325; https://doi.org/10.3390/agriculture14081325 - 9 Aug 2024
Viewed by 923
Abstract
This study discusses the issue of determining the direction and strength of changes taking place in the structure of agricultural land in the mountain and foothill areas of the Carpathians in Poland in comparison with Slovakia, Romania and Ukraine. The most important financial [...] Read more.
This study discusses the issue of determining the direction and strength of changes taking place in the structure of agricultural land in the mountain and foothill areas of the Carpathians in Poland in comparison with Slovakia, Romania and Ukraine. The most important financial institutional measures dedicated to the protection of the natural environment in Polish agriculture in the Areas facing Natural and other specific Constraints (ANCs) mountain and foothill in the first year of the CAP 2023–2027 were also established. Satellite data from 2001 to 2022 were used. The analyses used the land use classification MCD12Q1 provided by NASA and were made on the basis of satellite imagery collections from the MODIS sensor placed on two satellites: TERRA and AQUA. In EU countries, a decreasing trend in agricultural areas has been observed in areas below 350 m above sea level. In areas above 350 m, this trend weakened or even turned into an upward trend. Only in Ukraine was a different trend observed. It was found that in Poland, the degree of involvement of farmers from mountain and foothill areas in implementing financial institutional measures dedicated to protecting the natural environment during the study period was not satisfactory. Full article
Show Figures

Figure 1

Figure 1
<p>Scheme of the analysis of agriculture within separate groups of communes due to the fact and nuisance of ANCs mountain and foothill in Poland. Source: own study.</p>
Full article ">Figure 2
<p>Distribution of communes with different shares of ANCs mountain and foothill in Poland. Source: own study ISSPC SRI; IAFE NRI.</p>
Full article ">Figure 3
<p>Land use in the Carpathians in 2001 and 2022. Source: own study based on MODIS.</p>
Full article ">Figure 4
<p>Trends in the percentage share [%] of the total agricultural area and cropland in the total area of land in the Carpathians in 2001–2022. Source: own study based on MODIS.</p>
Full article ">Figure 5
<p>Number of farms participating in practices under eco-schemes, in organic and agri–environment–climate measures in communes with different shares of ANCs mountain and foothill in Poland in 2023. Source: own study based on ARMA.</p>
Full article ">Figure 6
<p>Share of [%] farms with eco-schemes in total number of farms in communes with ANCs mountain and foothill in 2023. Source: own study based on ARMA.</p>
Full article ">Figure 7
<p>Share of [%] farms with organic and agri–environmental–climate measure in total number of farms in communes with ANCs mountain and foothill in 2023.</p>
Full article ">Figure 8
<p>Agricultural area covered by practices under eco-schemes, ecological and agri–environment–climate measures in communes with different shares of ANCs mountain and foothill in Poland in 2023. Source: own study based on ARMA.</p>
Full article ">Figure 9
<p>Share [%] of UAA in farms with eco-schemes in total UAA in communes with ANCs mountain and foothill in 2023.</p>
Full article ">Figure 10
<p>Share [%] of UAA covered by organic and agri–environmental–climate measures in total UAA in communes with ANCs mountain and foothill in 2023.</p>
Full article ">
25 pages, 6444 KiB  
Article
Long-Term Evaluation of Aerosol Optical Properties in the Levantine Region: A Comparative Analysis of AERONET and Aqua/MODIS
by Ayse Gokcen Isik, S. Yeşer Aslanoğlu and Gülen Güllü
Remote Sens. 2024, 16(14), 2651; https://doi.org/10.3390/rs16142651 - 20 Jul 2024
Viewed by 944
Abstract
The focus on aerosol analysis in the Levantine Region is driven by climate-change impacts, the region’s increasing urban development and industrial activities, and its geographical proximity to major dust-source areas. This study conducts a comparative analysis of aerosol optical depth data from Aqua/MODIS [...] Read more.
The focus on aerosol analysis in the Levantine Region is driven by climate-change impacts, the region’s increasing urban development and industrial activities, and its geographical proximity to major dust-source areas. This study conducts a comparative analysis of aerosol optical depth data from Aqua/MODIS and AERONET during different periods between 2003 and 2023 at four stations: IMS-METU-ERDEMLI (Mersin/Türkiye) (2004–2019), CUT-TEPAK (Limassol/Cyprus) (2010–2023), Cairo_EMA_2 (Cairo/Egypt) (2010–2023), and SEDE_BOKER (Sede Boker/Israel) (2003–2023). The objective is to evaluate the variability and reliability of AOD measurements between satellite and ground-based observations and to determine how well they represent regional climatology. The highest percentage of measurements within the expected error envelope was observed at the IMS-METU-ERDEMLI station, indicating the best agreement between MODIS and AERONET data at this location. The Seasonal-Trend Decomposition using Loess (STL) method revealed consistent spring and summer peaks influenced by dust transport from the Sahara and the Middle East, with lower values in winter. The study also considers the influence of cloud fraction on MODIS measurements and includes aerosol classification. A statistically significant slight positive trend in AOD values was identified at the IMS-METU-ERDEMLI station. Conversely, no significant trends were detected at the other stations. The results of this study agree with those of previous research on the impact of long-range dust transport on regional aerosol loadings, emphasizing the importance of integrating satellite and ground-based observations. Full article
Show Figures

Figure 1

Figure 1
<p>Levantine Region with the AERONET stations used for the study.</p>
Full article ">Figure 2
<p>MODIS AOD versus AERONET AOD at Erdemli Station from 2004 to 2019. The red line is the regression line, and the green dashed lines define the envelope of expected error (EE). Colors of the points represent the cloud-fraction range associated with MODIS measurements.</p>
Full article ">Figure 3
<p>MODIS AOD versus AERONET AOD at Cyprus station from 2010 to 2023. The red line is the regression line, and the green dashed lines define the envelope of expected error (EE). Colors of the points represent the cloud-fraction range associated with MODIS measurements.</p>
Full article ">Figure 4
<p>MODIS AOD versus AERONET AOD at Cairo station from 2010 to 2023. The red line is the regression line, and the green dashed lines define the envelope of expected error (EE). Colors of the points represent the cloud-fraction range associated with MODIS measurements.</p>
Full article ">Figure 5
<p>MODIS AOD versus AERONET AOD at Sede Boker station from 2003 to 2023. The red line is the regression line, and the green dashed lines define the envelope of expected error (EE). Colors of the points represent the cloud-fraction range associated with MODIS measurements.</p>
Full article ">Figure 6
<p>Comparison of AOD measurements from MODIS to those from AERONET at IMS-METU-ERDEMLI, CUT-TEPAK, Cairo_EMA_2, and SEDE_BOKER stations over a decade, from 2010 to 2019.</p>
Full article ">Figure 7
<p>Seasonal comparison between MODIS and AERONET AOD measurements across four locations. The seasons are abbreviated as (<b>a</b>) DJF (December, January, February), (<b>b</b>) MAM (March, April, May), (<b>c</b>) JJA (June, July, August), and (<b>d</b>) SON (September, October, November).</p>
Full article ">Figure 8
<p>Histograms of the MODIS and AERONET measurements at 550 nm Across Different AERONET Stations: (<b>a</b>) IMS-METU-ERDEMLI, (<b>b</b>) CUT-TEPAK, (<b>c</b>) Cairo_EMA_2, (<b>d</b>) SEDE_BOKER.</p>
Full article ">Figure 9
<p>Aerosol classification for the IMS-METU-ERDEMLI station in Erdemli/Mersin (Time period: 2004–2019).</p>
Full article ">Figure 10
<p>Aerosol classification for the CUT-TEPAK station in Limassol/Cyprus (Time period: 2010–2023).</p>
Full article ">Figure 11
<p>Aerosol classification for the Cairo_EMA_2 station in Cairo/Egypt (Time period: 2010–2023).</p>
Full article ">Figure 12
<p>Aerosol classification for the SEDE_BOKER station in Sede Boker/Israel (Time period: 2003–2023).</p>
Full article ">Figure 13
<p>Trends in MODIS and AERONET AOD at 550 nm for four different stations: (<b>a</b>) IMS-METU-ERDEMLI (2004–2019), (<b>b</b>) CUT-TEPAK (2010–2023), (<b>c</b>) Cairo_EMA_2 (2010–2023), and (<b>d</b>) SEDE_BOKER (2003–2023).</p>
Full article ">Figure A1
<p>Scatterplot of MODIS AOD versus AERONET AOD for different cloud-fraction intervals (CF) across four different sites: (<b>a</b>–<b>c</b>) IMS-METU-ERDEMLI, (<b>d</b>–<b>f</b>) CUT-TEPAK, (<b>g</b>–<b>i</b>) Cairo_EMA, and (<b>j</b>–<b>l</b>) SEDE_BOKER. Each station was analyzed for three cloud-fraction intervals: CF (0–0.3), CF (0.3–0.6), and CF (0.6–1.0).</p>
Full article ">Figure A2
<p>STL decomposition of MODIS AOD time-series data at 550 nm for four different sites: (<b>a</b>) IMS-METU-ERDEMLI, (<b>b</b>) CUT-TEPAK, (<b>c</b>) Cairo_EMA_2, and (<b>d</b>) SEDE_BOKER. Each subplot includes the original data, trend component, seasonal component, and residuals from 2003 to 2023.</p>
Full article ">
16 pages, 1280 KiB  
Article
Are Harmful Algal Blooms Increasing in the Great Lakes?
by Karl R. Bosse, Gary L. Fahnenstiel, Cal D. Buelo, Matthew B. Pawlowski, Anne E. Scofield, Elizabeth K. Hinchey and Michael J. Sayers
Water 2024, 16(14), 1944; https://doi.org/10.3390/w16141944 - 10 Jul 2024
Viewed by 977
Abstract
This study used satellite remote sensing to investigate trends in harmful algal blooms (HABs) over the last 21 years, focusing on four regions within the Laurentian Great Lakes: western Lake Erie, Green Bay, Saginaw Bay, and western Lake Superior. HABs in the water [...] Read more.
This study used satellite remote sensing to investigate trends in harmful algal blooms (HABs) over the last 21 years, focusing on four regions within the Laurentian Great Lakes: western Lake Erie, Green Bay, Saginaw Bay, and western Lake Superior. HABs in the water column were identified from remote sensing-derived chlorophyll concentrations, and surface HAB scums were classified based on the Normalized Difference Vegetation Index (NDVI) band ratio index. Using imagery from the Moderate Resolution Imaging Spectroradiometer sensor on the Aqua satellite (MODIS-Aqua) from 2002 to 2022, we generated daily estimates of the HAB and surface scum extents for each region, which were then averaged to generate mean annual extents. We observed a significant decline in the Saginaw Bay mean annual HAB extents over the 21-year study period. Otherwise, no significant changes were observed over this period in any region for either the HAB or surface scum mean annual extents, thus suggesting that HABs are not increasing in the Great Lakes. Despite the lack of increasing trends, the blooms are still recurring annually and causing a negative impact on the nearby communities; thus, we believe that it is crucial to continue studying Great Lakes HABs to monitor the impact of current and future abatement strategies. Full article
Show Figures

Figure 1

Figure 1
<p>The regions studied for historical HAB extents include western Lake Erie (WLE), Saginaw Bay, Green Bay, and western Lake Superior.</p>
Full article ">Figure 2
<p>Mean annual HAB extents for WLE (panel (<b>A</b>)), Saginaw Bay (panel (<b>B</b>)), and Green Bay (panel (<b>C</b>)). Error bars represent the standard error of the mean. Only Saginaw Bay had a significant trend over the data record (slope = −2.1 km<sup>2</sup>/yr; <math display="inline"><semantics> <mrow> <mi>p</mi> <mo>&lt;</mo> <mn>0.05</mn> </mrow> </semantics></math>).</p>
Full article ">Figure 3
<p>Mean annual surface HAB scum extents for WLE (panel (<b>A</b>)), Saginaw Bay (panel (<b>B</b>)), and Green Bay (panel (<b>C</b>)). Error bars represent the standard error of the mean.</p>
Full article ">Figure 4
<p>Mean annual HAB extent plotted against mean annual surface HAB scum extent for WLE (panel (<b>A</b>)), Saginaw Bay (panel (<b>B</b>)), and Green Bay (panel (<b>C</b>)). Only Green Bay showed a significant correlation between these indices (<span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">Figure 5
<p>Mean annual algal bloom extent for Lake Superior. Error bars represent the standard error of the mean.</p>
Full article ">
29 pages, 17646 KiB  
Article
Dust Events over the Urmia Lake Basin, NW Iran, in 2009–2022 and Their Potential Sources
by Abbas Ranjbar Saadat Abadi, Karim Abdukhakimovich Shukurov, Nasim Hossein Hamzeh, Dimitris G. Kaskaoutis, Christian Opp, Lyudmila Mihailovna Shukurova and Zahra Ghasabi
Remote Sens. 2024, 16(13), 2384; https://doi.org/10.3390/rs16132384 - 28 Jun 2024
Viewed by 781
Abstract
Nowadays, dried lake beds constitute the largest source of saline dust storms, with serious environmental and health issues in the surrounding areas. In this study, we examined the spatial–temporal distribution of monthly and annual dust events of varying intensity (dust in suspension, blowing [...] Read more.
Nowadays, dried lake beds constitute the largest source of saline dust storms, with serious environmental and health issues in the surrounding areas. In this study, we examined the spatial–temporal distribution of monthly and annual dust events of varying intensity (dust in suspension, blowing dust, dust storms) in the vicinity of the desiccated Urmia Lake in northwestern (NW) Iran, based on horizontal visibility data during 2009–2022. Dust in suspension, blowing dust and dust storm events exhibited different monthly patterns, with higher frequencies between March and October, especially in the southern and eastern parts of the Urmia Basin. Furthermore, the intra-annual variations in aerosol optical depth at 500 nm (AOD550) and Ångström exponent at 412/470 nm (AE) were investigated using Terra/Aqua MODIS (Moderate Resolution Imaging Spectroradiometer) data over the Urmia Lake Basin (36–39°N, 44–47°E). Monthly distributions of potential coarse aerosol (AE < 1) sources affecting the lower troposphere over the Urmia Basin were reconstructed, synergizing Terra/Aqua MODIS AOD550 for AE < 1 values and HYSPLIT_4 backward trajectories. The reconstructed monthly patterns of the potential sources were compared with the monthly spatial distribution of Terra MODIS AOD550 in the Middle East and Central Asia (20–70°E, 20–50°N). The results showed that deserts in the Middle East and the Aral–Caspian arid region (ACAR) mostly contribute to dust aerosol load over the Urmia Lake region, exhibiting higher frequency in spring and early summer. Local dust sources from dried lake beds further contribute to the dust AOD, especially in the western part of the Urmia Basin during March and April. The modeling (DREAM8-NMME-MACC) results revealed high concentrations of near-surface dust concentrations, which may have health effects on the local population, while distant sources from the Middle East are the main controlling factors to aerosol loading over the Urmia Basin. Full article
Show Figures

Figure 1

Figure 1
<p>The topography of the study area and locations of 26 weather stations in the Urmia Basin, NW Iran. Black circles represent the main synoptic stations (MS) with 8 observations per day, while open circles symbolize the intermediate synoptic stations (IS) with 5 observations per day.</p>
Full article ">Figure 2
<p>A snapshot of the Urmia Lake region (September 2003). The cells of a 1° × 1° size are shown and were used to extract data series of the Terra/Aqua MODIS AOD<sub>550</sub> and Ångström exponent from 2009 to 2022. Parts of the dried bottom covered with salt are seen at the east (in N-cell) and southeast (CNTR-cell) boundaries of the lake.</p>
Full article ">Figure 3
<p>Monthly number of dust events with horizontal visibilities less than or equal to 5 km for widespread suspended dust (code 06) during 2009–2022.</p>
Full article ">Figure 4
<p>Monthly number of dust events with horizontal visibilities less than or equal to 5 km for blowing dust and dust storms (BDSS, as defined in <a href="#remotesensing-16-02384-t001" class="html-table">Table 1</a>) during 2009–2022.</p>
Full article ">Figure 5
<p>Number of horizontal visibilities reduced due to severe (<b>a</b>,<b>d</b>), moderate (<b>b</b>,<b>e</b>) and weak (<b>c</b>,<b>f</b>) intensity of WSD (WW06) and BDSS (WW_Oth).</p>
Full article ">Figure 6
<p>The 80th, 90th and 95th percentiles, as well as the maximum DSC (<math display="inline"><semantics> <mrow> <mrow> <mrow> <mi mathvariant="sans-serif">μ</mi> <mi mathvariant="normal">g</mi> </mrow> <mo>/</mo> <mrow> <msup> <mrow> <mi mathvariant="normal">m</mi> </mrow> <mrow> <mn>3</mn> </mrow> </msup> </mrow> </mrow> </mrow> </semantics></math>) values of the surface concentration of suspended particles, predicted by the DREAM8-NMME-MACC model for urban and rural points in the study area, based on the range of changes in the AQI for PM<sub>10</sub> (<a href="#remotesensing-16-02384-t003" class="html-table">Table 3</a>).</p>
Full article ">Figure 7
<p>Intra-annual variations in the monthly average Ångström exponent (lines) and AOD<sub>550</sub> values (bars) from Terra MODIS over grid cells around the Urmia Basin during 2009–2022. (<b>a</b>) NW-cell. (<b>b</b>) W-cell. (<b>c</b>) SW-cell. (<b>d</b>) S-cell. (<b>e</b>) SE-cell. (<b>f</b>) E-cell. (<b>g</b>) NE-cell. (<b>h</b>) N-cell (evaluative graph). (<b>i</b>). CNTR-cell (evaluative graph). Red and light pink (for N and CNTR cells) lines indicate the annual mean AOD<sub>550</sub> in 2009–2022.</p>
Full article ">Figure 8
<p>Intra-annual variations in the monthly average AOD<sub>550</sub> values for the cases of AE &lt; 1 (gray bars) and monthly total number of AE &lt; 1 cases (blue bars) in 2009–2022 by Terra MODIS data. (<b>a</b>) NW-cell. (<b>b</b>) W-cell. (<b>c</b>) SW-cell. (<b>d</b>) S-cell. (<b>e</b>) SE-cell. (<b>f</b>) E-cell. (<b>g</b>) NE-cell. (<b>h</b>) N-cell (evaluative graph). (<b>i</b>) CNTR-cell (evaluative graph). Red lines represent the mean annual AOD<sub>550</sub> due to coarse aerosols only.</p>
Full article ">Figure 9
<p>Monthly mean regional ABL contribution to AOD<sub>550</sub> at AE &lt; 1.0 (CWT-AOD<sub>550</sub>) according to Aqua and Terra MODIS data for the Urmia Lake region (marked with a black rectangle) in 2009–2022: The spatial resolution is 0.5° × 0.5°. Here and after, the Aral Sea is shown in its boundaries of the 1960s.</p>
Full article ">Figure 10
<p>The spatial distributions of monthly mean Terra MODIS Deep Blue AOD<sub>550</sub> during 2009–2022. The Urmia Lake region is marked with a black rectangle. The spatial resolution is 1° × 1°.</p>
Full article ">Figure A1
<p>Intra-annual variations in the monthly average Ångström exponent (lines) and the monthly average AOD<sub>550</sub> (bars) over the following cells in 2009–2021 by Aqua MODIS data. (<b>a</b>) NW-cell. (<b>b</b>) W-cell. (<b>c</b>) SW-cell. (<b>d</b>) S-cell. (<b>e</b>) SE-cell. (<b>f</b>) E-cell. (<b>g</b>) NE-cell. (<b>h</b>) N-cell (evaluative graph). (<b>i</b>) CNTR-cell (evaluative graph). Blue and red (light pink for N and CNTR cells) lines indicate, respectively, the annual mean Ångström exponent and AOD<sub>550</sub> values in 2009–2022.</p>
Full article ">Figure A2
<p>Intra-annual variations in the monthly average AOD<sub>550</sub> (gray bars) and monthly total number of coarse aerosol (Ångström exponent &lt; 1) events N (blue bars) in 2009–2022 by Aqua MODIS data. (<b>a</b>) NW-cell. (<b>b</b>) W-cell. (<b>c</b>) SW-cell. (<b>d</b>) S-cell. (<b>e</b>) SE-cell. (<b>f</b>) E-cell. (<b>g</b>) NE-cell. (<b>h</b>) N-cell (evaluative graph). (<b>i</b>) CNTR-cell (evaluative graph). Red (and light pink for N and CNTR cells) lines represent the mean annual AOD<sub>550</sub> averaged by coarse aerosol events.</p>
Full article ">Figure A3
<p>The spatial distributions of the monthly mean Terra-MODIS AOD<sub>550</sub> (Combined Deep Blue/Dark Target) during 2009–2022. The Urmia Lake region is marked with a black rectangle. The spatial resolution is 1° × 1°.</p>
Full article ">
19 pages, 8018 KiB  
Article
Characteristics of Yellow Sea Fog under the Influence of Eastern China Aerosol Plumes
by Jiakun Liang and Jennifer D. Small Griswold
Remote Sens. 2024, 16(13), 2262; https://doi.org/10.3390/rs16132262 - 21 Jun 2024
Viewed by 669
Abstract
Sea fog is a societally relevant phenomenon that occurs under the influence of specific oceanic and atmospheric conditions including aerosol conditions. The Yellow Sea region in China regularly experiences sea fog events, of varying intensity, that impact coastal regions and maritime activities. The [...] Read more.
Sea fog is a societally relevant phenomenon that occurs under the influence of specific oceanic and atmospheric conditions including aerosol conditions. The Yellow Sea region in China regularly experiences sea fog events, of varying intensity, that impact coastal regions and maritime activities. The occurrence and structure of fog are impacted by the concentration of aerosols in the air where the fog forms. Along with industrial development, air pollution has become a serious environmental problem in Northeastern China. These higher pollution levels are confirmed by various satellite remote sensing instruments including the Moderate Resolution Imaging Spectroradiometer (MODIS) aboard the Aqua satellite that observes aerosol and cloud properties. These observations show a clear influence of aerosol loading over the Yellow Sea region, which can impact regional sea fog. In this study, high-resolution data sets from MODIS Aqua L2 are used to investigate the relationships between cloud properties and aerosol features. Using a bi-variate comparison method, we find that, for most cases, larger values of COT (cloud optical thickness) are related to both a smaller DER (droplet effective radius) and higher CTH (cloud top height). However, in the cases where fog is thinner with many zero values in CTH, the larger COT is related to both a smaller DER and CTH. For fog cases where the aerosol type is dominated by smoke (e.g., confirmed fire activities in the East China Plain), the semi-direct effect is indicated and may play a role in determining fog structure such that a smaller DER corresponds with thinner fog and smaller COT values. Full article
(This article belongs to the Section Atmospheric Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>MODIS Aqua L1B Granule Images highlighting different fog case scenarios. (<b>a</b>) Fog case on 2 May 2020, red box: “incomplete” fog area, the upper portion of the Yellow Sea is not included in the MODIS granule. (<b>b</b>) Fog case on 31 July 2020, cyan box: fog area covered by high cloud. (<b>c</b>) Fog case on 28 March 2012, yellow box: pollution (aerosol) band visible on and offshore.</p>
Full article ">Figure 2
<p>An example of CTH modification for the fog case was on 13 May 2018. (<b>a</b>) The original CTH of 5 km resolution from MODIS Aqua L2 cloud data product, (<b>b</b>) the modified CTH.</p>
Full article ">Figure 3
<p>MODIS Aqua L1B Granule Image of a fog case on 13 May 2018 (<b>a</b>), and CTH of (<b>b</b>) mean temperature inversion height 633 m, (<b>c</b>) 700 m, (<b>d</b>) 800 m. (<b>e</b>) DER at 1 km resolution from MODIS Aqua L2 cloud data product, (<b>f</b>) result of the CTH for the selected fog area after applying the DER mask and land-sea mask.</p>
Full article ">Figure 4
<p>Terrestrial aerosol types surrounding the Yellow Sea region from the MODIS Aqua L2 aerosol data product. (<b>a</b>) Fog case on 23 May 2006, main aerosol type: sulfate and dust. (<b>b</b>) Fog case on 8 June 2007, main aerosol type: heavy absorbing smoke and sulfate. (<b>c</b>) Fog case on 2 May 2008, main aerosol type: sulfate. (<b>d</b>) Fog case on 3 May 2009, main aerosol type: sulfate and dust. (<b>e</b>) Fog case on 4 May 2009, main aerosol type: sulfate and dust. (<b>f</b>) Fog case on 17 May 2011, main aerosol type: sulfate. (<b>g</b>) Fog case on 1 June 2011, main aerosol type: heavy absorbing smoke, dust, and sulfate. (<b>h</b>) Fog case on 28 March 2012, main aerosol type: sulfate. (<b>i</b>) Fog case on 8 April 2014, main aerosol type: sulfate. (<b>j</b>) Fog case on 9 April 2014, main aerosol type: sulfate and dust. (<b>k</b>) Fog case on 10 April 2016, main aerosol type: sulfate and dust. (<b>l</b>) Fog case on 13 April 2016, main aerosol type: sulfate and dust. (<b>m</b>) Fog case on 14 April 2016, main aerosol type: sulfate and dust. (<b>n</b>) Fog case on 13 May 2018, main aerosol type: sulfate and dust. (<b>o</b>) Fog case on 6 June 2018, main aerosol type: heavy absorbing smoke and sulfate.</p>
Full article ">Figure 5
<p>Fog cases with fire occurrences around the Shandong Peninsula on 8 June 2007 (the first column, (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>,<b>m</b>)), 1 June 2011 (the second column (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>,<b>n</b>)), and 6 June 2018 (the third column (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>,<b>o</b>)). (<b>a</b>–<b>c</b>) Satellite RGB visible image from MODIS L2B Granule Image, red box: pollution band. (<b>d</b>–<b>f</b>) Thermal indicators of fire from NASA World View. (<b>g</b>–<b>i</b>) Vertical structures of air temperature (blue line) and dew point temperature (red line) from Sounding files at Qingdao Station. (<b>j</b>–<b>l</b>) Temperature advection calculated from the NECP/NCAR reanalysis data. The black line indicates the geopotential height at 1000 mb, the black arrows indicate the wind direction at the speed of 10 m/s unit, and the red (blue) areas indicate the warm (cold) temperature advection. (<b>m</b>–<b>o</b>) AOD from MODIS Aqua L2 aerosol data product.</p>
Full article ">Figure 6
<p>Bi-variate comparison for 15 fog cases. Diagonal Pattern (<b>a</b>,<b>d</b>–<b>j</b>,<b>l</b>–<b>n</b>) refers to distributions with larger COT values corresponding to smaller DER values and larger CTH values. Left-Right Pattern (<b>c</b>,<b>k</b>) refers to distributions with larger COT values corresponding to larger DER values and smaller CTH values. Inverse-Diagonal Pattern (<b>b</b>,<b>o</b>) refers to distributions with larger COT values corresponding to both larger DER values and larger CTH values.</p>
Full article ">Figure 7
<p>Aerosol, wind conditions, and cloud properties for the sea fog case on 28 March 2012, from MODIS Aqua L2 cloud data. (<b>a</b>) AOD form MODIS Aqua L2 aerosol data product. (<b>b</b>) Surface wind from NCEP/NCAR reanalysis dataset. (<b>c</b>) DER. (<b>d</b>) CTH. (<b>e</b>) COT.</p>
Full article ">Figure 8
<p>CTH from the MODIS Aqua L2 cloud data product. (<b>a</b>) Fog case on 2 May 2008. (<b>b</b>) Fog case on 10 April 2016.</p>
Full article ">Figure 9
<p>Cloud properties and aerosol for the sea fog case on 8 June 2007, from the MODIS Aqua L2 cloud data. (<b>a</b>) DER. (<b>b</b>) CTH. (<b>c</b>) COT.</p>
Full article ">Figure 10
<p>The sum bi-variate comparison of the 15 fog cases.</p>
Full article ">
18 pages, 4918 KiB  
Article
Assessment of Accuracy of Moderate-Resolution Imaging Spectroradiometer Sea Surface Temperature at High Latitudes Using Saildrone Data
by Chong Jia, Peter J. Minnett and Malgorzata Szczodrak
Remote Sens. 2024, 16(11), 2008; https://doi.org/10.3390/rs16112008 - 3 Jun 2024
Viewed by 1171
Abstract
The infrared (IR) satellite remote sensing of sea surface skin temperature (SSTskin) is challenging in the northern high-latitude region, especially in the Arctic because of its extreme environmental conditions, and thus the accuracy of SSTskin retrievals is questionable. Several Saildrone [...] Read more.
The infrared (IR) satellite remote sensing of sea surface skin temperature (SSTskin) is challenging in the northern high-latitude region, especially in the Arctic because of its extreme environmental conditions, and thus the accuracy of SSTskin retrievals is questionable. Several Saildrone uncrewed surface vehicles were deployed at the Pacific side of the Arctic in 2019, and two of them, SD-1036 and SD-1037, were equipped with a pair of IR pyrometers on the deck, whose measurements have been shown to be useful in the derivation of SSTskin with sufficient accuracy for scientific applications, providing an opportunity to validate satellite SSTskin retrievals. This study aims to assess the accuracy of MODIS-retrieved SSTskin from both Aqua and Terra satellites by comparisons with collocated Saildrone-derived SSTskin data. The mean difference in SSTskin from the SD-1036 and SD-1037 measurements is ~0.4 K, largely resulting from differences in the atmospheric conditions experienced by the two Saildrones. The performance of MODIS on Aqua and Terra in retrieving SSTskin is comparable. Negative brightness temperature (BT) differences between 11 μm and 12 μm channels are identified as being physically based, but are removed from the analyses as they present anomalous conditions for which the atmospheric correction algorithm is not suited. Overall, the MODIS SSTskin retrievals show negative mean biases, −0.234 K for Aqua and −0.295 K for Terra. The variations in the retrieval inaccuracies show an association with diurnal warming events in the upper ocean from long periods of sunlight in the Arctic. Also contributing to inaccuracies in the retrieval is the surface emissivity effect in BT differences characterized by the Emissivity-introduced BT difference (EΔBT) index. This study demonstrates the characteristics of MODIS-retrieved SSTskin in the Arctic, at least at the Pacific side, and underscores that more in situ SSTskin data at high latitudes are needed for further error identification and algorithm development of IR SSTskin. Full article
Show Figures

Figure 1

Figure 1
<p>The cruise trajectories of two NASA-funded Saildrone vehicles, SD-1036 (white) and SD-1037 (magenta), deployed during the 2019 Arctic Cruise from 15 May to 11 October. The background SST map is taken from the Multiscale Ultrahigh Resolution (MUR) Level-4 SST analysis data [<a href="#B37-remotesensing-16-02008" class="html-bibr">37</a>] on 16 September 2019. The subplot is a picture of the Saildrones at the starting point, which is courtesy of Saildrone Inc.</p>
Full article ">Figure 2
<p>The Aqua MODIS–Saildrone SST<sub>skin</sub> difference as a function of the (<b>a</b>) distance and (<b>b</b>) time difference in the allowed spatial–temporal intervals in the matchup criteria. Data include both SD-1036 and SD-1037 measurements. The one-to-one matchups were determined based on the smallest separation between the Saildrone measurement and MODIS pixel. The black linear fitted lines are given with the expression in the top right corner. (<b>c</b>,<b>d</b>) are similar to (<b>a</b>,<b>b</b>), but for the one-to-one matchups determined by the closest timestamp. The regressions in (<b>a</b>,<b>c</b>) pass the significance test at a 95% confidence level, but slopes in (<b>b</b>,<b>d</b>) are not significantly different from zero.</p>
Full article ">Figure 3
<p>(<b>a</b>) Histogram (normal distribution fitted curve in blue) of the Aqua MODIS–Saildrone SST<sub>skin</sub> difference and (<b>b</b>) the scatter plot of Saildrone- and MODIS-derived SST<sub>skin</sub> colored by the data density. (<b>c</b>,<b>d</b>) are similar to (<b>a</b>,<b>b</b>), but for Terra MODIS–Saildrone matchups.</p>
Full article ">Figure 4
<p>Histograms of the (<b>a</b>) Aqua MODIS BT difference between 11 μm and 12 μm channels and (<b>b</b>) air–sea temperature difference (ASTD) for the matchup data during the SD-1036 (light blue) and SD-1037 (light red) cruises. (<b>c</b>) The data density scatter plot of the BT difference and ASTD in Aqua MODIS–Saildrone matchups with the fitted dashed line.</p>
Full article ">Figure 5
<p>Maps of ASTD for Aqua MODIS–Saildrone matchups for (<b>a</b>) SD-1036 and (<b>b</b>) SD-1037. (<b>c</b>) The bivariate histogram for the longitude and latitude of the Aqua MODIS pixels matched with the SD-1036 (blue) and SD-1037 (red) measurements. The marks on some red columns indicate the heights of corresponding blue bars overwhelmed by the red ones.</p>
Full article ">Figure 6
<p>Reanalysis data from MERRA-2 matched with the Aqua MODIS–Saildrone matchups during the SD-1036 and SD-1037 cruises showing the vertical profiles of (<b>a</b>) specific humidity and (<b>b</b>) air temperature plotted as the mean (line and dots) ±1 standard deviation (envelope), as well as (<b>c</b>) the histogram of the total column water vapor.</p>
Full article ">Figure 7
<p>Histogram of RTTOV-simulated BT difference between 11 μm and 12 μm for Aqua MODIS pixels matched with Saildrone measurements during SD-1036 and SD-1037 cruises.</p>
Full article ">Figure 8
<p>Scatter plots (colored by data density) of the Aqua MODIS–Saildrone SST<sub>skin</sub> difference as a function of (<b>a</b>) the amplitude of diurnal warming with a fitted black dashed line when diurnal warming exists and (<b>b</b>) the Emissivity-introduced BT difference (EΔBT) with red dots and error bars indicating the mean and STD of temperature differences, calculated at 0.16 K intervals. The histogram distributions of diurnal warming and EΔBT are also plotted as the background for the data during SD-1036 and SD-1037 cruises separately.</p>
Full article ">
17 pages, 3972 KiB  
Article
Quantitative Assessment of Volcanic Thermal Activity from Space Using an Isolation Forest Machine Learning Algorithm
by Claudia Corradino, Arianna Beatrice Malaguti, Micheal S. Ramsey and Ciro Del Negro
Remote Sens. 2024, 16(11), 2001; https://doi.org/10.3390/rs16112001 - 1 Jun 2024
Cited by 2 | Viewed by 1337
Abstract
Understanding the dynamics of volcanic activity is crucial for volcano observatories in their efforts to forecast volcanic hazards. Satellite imager data hold promise in offering crucial insights into the thermal behavior of active volcanoes worldwide, facilitating the assessment of volcanic activity levels and [...] Read more.
Understanding the dynamics of volcanic activity is crucial for volcano observatories in their efforts to forecast volcanic hazards. Satellite imager data hold promise in offering crucial insights into the thermal behavior of active volcanoes worldwide, facilitating the assessment of volcanic activity levels and identifying significant changes during periods of volcano unrest. The Moderate Resolution Imaging Spectroradiometer (MODIS) sensor, aboard NASA’s Terra and Aqua satellites, provides invaluable data with high temporal and spectral resolution, enabling comprehensive thermal monitoring of eruptive activity. The accuracy of volcanic activity characterization depends on the quality of models used to relate the relationship between volcanic phenomena and target variables such as temperature. Under these circumstances, machine learning (ML) techniques such as decision trees can be employed to develop reliable models without necessarily offering any particular or explicit insights. Here, we present a ML approach for quantifying volcanic thermal activity levels in near real time using thermal infrared satellite data. We develop an unsupervised Isolation Forest machine learning algorithm, fully implemented in Google Colab using Google Earth Engine (GEE) which utilizes MODIS Land Surface Temperature (LST) data to automatically retrieve information on the thermal state of volcanoes. We evaluate the algorithm on various volcanoes worldwide characterized by different levels of volcanic activity. Full article
Show Figures

Figure 1

Figure 1
<p>Overview and representative image of the volcanoes studied. All images were captured via Google Earth Pro [<a href="http://www.earth.google.com" target="_blank">http://www.earth.google.com</a>; accessed on 27 March 2024] and QGIS [<a href="https://qgis.org/it/site/" target="_blank">https://qgis.org/it/site/</a>; accessed on 27 March 2024].</p>
Full article ">Figure 2
<p>Workflow of the proposed three-step approach.</p>
Full article ">Figure 3
<p>NLS adaptation phase to the anomaly detection phase for Etna. The reconstructed signal over the real thermal signal (<b>a</b>), the reconstructed error (<b>b</b>), the LST temperature of the detected anomalies (<b>c</b>), and the real thermal signal identified as anomalous over the real thermal signal (<b>d</b>).</p>
Full article ">Figure 4
<p>Time series of the temperature above average of the detected volcanic anomalies for Etna, Klyuchevskoy, and Lascar. Blue arrows indicate the starting of an eruption, while the gray shaded area indicates time windows characterized by an increase in activity preceding the eruption.</p>
Full article ">Figure 5
<p>Time series of the temperature above average of the detected volcanic anomalies for Fuego, Popocatépetl, and Stromboli. Blue arrows indicate the starting of an eruption, while the gray shaded area indicates time windows characterized by an increase in activity preceding the eruption.</p>
Full article ">Figure 6
<p>Time series of the temperature above average of the detected volcanic anomalies and the activity levels low (green bars), moderate (orange bars) and high (red bars) for Etna, Klyuchevskoy, Lascar, Fuego, Popocatépetl, and Stromboli.</p>
Full article ">
25 pages, 8968 KiB  
Article
Consistency of Aerosol Optical Properties between MODIS Satellite Retrievals and AERONET over a 14-Year Period in Central–East Europe
by Lucia-Timea Deaconu, Alexandru Mereuță, Andrei Radovici, Horațiu Ioan Ștefănie, Camelia Botezan and Nicolae Ajtai
Remote Sens. 2024, 16(10), 1677; https://doi.org/10.3390/rs16101677 - 9 May 2024
Viewed by 1134
Abstract
Aerosols influence Earth’s climate by interacting with radiation and clouds. Remote sensing techniques aim to enhance our understanding of aerosol forcing using ground-based and satellite retrievals. Despite technological advancements, challenges persist in reducing uncertainties in satellite remote sensing. Our study examines retrieval biases [...] Read more.
Aerosols influence Earth’s climate by interacting with radiation and clouds. Remote sensing techniques aim to enhance our understanding of aerosol forcing using ground-based and satellite retrievals. Despite technological advancements, challenges persist in reducing uncertainties in satellite remote sensing. Our study examines retrieval biases in MODIS sensors on Terra and Aqua satellites compared to AERONET ground-based measurements. We assess their performance and the correlation with the AERONET aerosol optical depth (AOD) using 14 years of data (2010–2023) from 29 AERONET stations across 10 Central–East European countries. The results indicate discrepancies between MODIS Terra and Aqua retrievals: Terra overestimates the AOD at 16 AERONET stations, while Aqua underestimates the AOD at 21 stations. The examination of temporal biases in the AOD using the calculated estimated error (ER) between AERONET and MODIS retrievals reveals a notable seasonality in coincident retrievals. Both sensors show higher positive AOD biases against AERONET in spring and summer compared to fall and winter, with few ER values for Aqua indicating poor agreement with AERONET. Seasonal variations in correlation strength were noted, with significant improvements from winter to summer (from R2 of 0.58 in winter to R2 of 0.76 in summer for MODIS Terra and from R2 of 0.53 in winter to R2 of 0.74 in summer for MODIS Aqua). Over the fourteen-year period, monthly mean aerosol AOD trends indicate a decrease of −0.00027 from AERONET retrievals and negative monthly mean trends of the AOD from collocated MODIS Terra and Aqua retrievals of −0.00023 and −0.00025, respectively. An aerosol classification analysis showed that mixed aerosols comprised over 30% of the total aerosol composition, while polluted aerosols accounted for more than 22%, and continental aerosols contributed between 22% and 24%. The remaining 20% consists of biomass-burning, dust, and marine aerosols. Based on the aerosol classification method, we computed the bias between the AERONET AE and MODIS AE, which showed higher AE values for AERONET retrievals for a mixture of aerosols and biomass burning, while for marine aerosols, the MODIS AE was larger and for dust the results were inconclusive. Full article
(This article belongs to the Special Issue Remote Sensing of Aerosols, Planetary Boundary Layer, and Clouds)
Show Figures

Figure 1

Figure 1
<p>Spatial distribution of averaged aerosol optical depth (AOD) derived from MODIS Terra dataset over Central–East Europe for four seasons, (<b>a</b>) winter, (<b>b</b>) spring, (<b>c</b>) summer, and (<b>d</b>) autumn, between 2010 and 2023. The map also indicates the location of the 29 AERONET stations used for comparison. The size of the circles for the AERONET stations represents the number of data points available at each station, while the color bars indicate AOD values, cutoff at 0.25 for both MODIS and AERONET measurements.</p>
Full article ">Figure 2
<p>Spatial distribution of averaged AOD derived from MODIS Aqua dataset over the same region and time period as in <a href="#remotesensing-16-01677-f001" class="html-fig">Figure 1</a>, (<b>a</b>) winter, (<b>b</b>) spring, (<b>c</b>) summer, and (<b>d</b>) autumn. The figure layout, including the AERONET dataset and colour bar scale, remains consistent with <a href="#remotesensing-16-01677-f001" class="html-fig">Figure 1</a>.</p>
Full article ">Figure 3
<p>Comparison between the AERONET aerosol optical depth (AOD) at 550 nm and collocated MODIS Terra AOD at 550 nm for the four seasons (<b>a</b>) winter—DJF, (<b>b</b>) spring—MAM, (<b>c</b>) summer—JJA, and (<b>d</b>) autumn—SON, spanning from 2010 to 2023 over Central–East Europe. The scatter plots depict the statistical correlation between the two datasets for each season across all AERONET stations, with marginal histograms illustrating the distribution of AOD values for both datasets. The colour map indicates the corresponding AERONET Ångström exponent (AE) calculated with AOD at 440 and 870 nm. Linear equation values (slope and intercept), R values, R-square values, RMSE values, and number of points are shown in the table for each season.</p>
Full article ">Figure 4
<p>Same as <a href="#remotesensing-16-01677-f003" class="html-fig">Figure 3</a>, comparing the AERONET AOD at 550 nm with the collocated MODIS Aqua AOD at 550 nm, (<b>a</b>) winter—DJF, (<b>b</b>) spring—MAM, (<b>c</b>) summer—JJA, and (<b>d</b>) autumn—SON. The scatter plots, marginal histograms, colour map, and statistical metrics remain consistent with <a href="#remotesensing-16-01677-f003" class="html-fig">Figure 3</a>.</p>
Full article ">Figure 5
<p>Time series of monthly mean error ratios (ER in Equation (1)) (<b>a</b>) and the number of collocations (<b>b</b>) for the collocated dataset from 29 selected AERONET stations in CEE, time period 2010–2023. The Terra record is in red, and the Aqua is in blue. The number of collocations is season dependent (low number during winter months). Horizontal red lines are the ER value at −1 and 1.</p>
Full article ">Figure 6
<p>The distribution of AOD differences between AERONET and MODIS (<b>a</b>) Terra and (<b>b</b>) Aqua presented as probability density functions (PDF). Each of the 29 coloured plots represents the PDF of the ΔAOD (difference in daily AOD at 550 nm) between AERONET and MODIS retrievals for a specific location spanning from 2010 to 2023. The x-axis represents the ΔAOD values, ranging from negative to positive, while values on the colour bar represent the mean MODIS AOD bias at each individual location. The blue shades indicate a negative mean ΔAOD, while the red shades correspond to a positive mean ΔAOD. The vertical black line indicates the value 0 of ΔAOD.</p>
Full article ">Figure 7
<p>MODIS Terra (<b>a</b>) and Aqua (<b>b</b>) mean biases with confidence intervals (95%) specific for each AERONET location. Colour code same as <a href="#remotesensing-16-01677-f004" class="html-fig">Figure 4</a>.</p>
Full article ">Figure 8
<p>(<b>a</b>) The monthly mean aerosol optical depth (AOD) at 550 nm (blue line) with standard deviation and confidence interval at 95% for monthly means, spanning January 2010 to June 2023, alongside the monthly AOD trend (red line). (<b>b</b>) Similar for MODIS Terra and (<b>c</b>) MODIS Aqua. The multiannual mean AOD and standard deviation with confidence intervals are also shown.</p>
Full article ">Figure 9
<p>Relationship between Ångström exponent (AE) at wavelengths 440/870 nm and aerosol optical depth (AOD) at 440 nm. Six aerosol classes are distinguished based on thresholds of AE and AOD: biomass (green), continental (red), dust (orange), marine (blue), mixed (purple), and polluted (black).</p>
Full article ">Figure 10
<p>Pie charts illustrating aerosol type classification for total (<b>a</b>), urban (<b>b</b>), and rural (<b>c</b>) locations. Associated bar charts depict monthly aerosol type distribution spanning 2010–2023 across Central–East Europe for the three categories: total (<b>a</b>), urban (<b>b</b>), and rural (<b>c</b>). Aerosol types represented include marine (blue), dust (yellow), continental (red), mixed (purple), polluted (black), and biomass burning (green).</p>
Full article ">Figure 11
<p>Ångström exponent biases between AERONET and MODIS Terra retrievals for 6 aerosol types: biomass burning, continental, dust, marine, mixed, and polluted. The colour bar shows the points density. Descriptive statistics for these categories is also presented.</p>
Full article ">
27 pages, 2083 KiB  
Article
A Wide-Angle Hyperspectral Top-of-Atmosphere Reflectance Model for the Libyan Desert
by Fuxiang Guo, Xiaobing Zheng, Yanna Zhang, Wei Wei, Zejie Zhang, Quan Zhang and Xin Li
Remote Sens. 2024, 16(8), 1406; https://doi.org/10.3390/rs16081406 - 16 Apr 2024
Viewed by 874
Abstract
Reference targets with stability, uniformity, and known reflectance on the Earth’s surface, such as deserts, can be used for the absolute radiometric calibration of satellite sensors. A wide-angle hyperspectral reflectance model at the top of atmosphere (TOA) over such a reference target will [...] Read more.
Reference targets with stability, uniformity, and known reflectance on the Earth’s surface, such as deserts, can be used for the absolute radiometric calibration of satellite sensors. A wide-angle hyperspectral reflectance model at the top of atmosphere (TOA) over such a reference target will expand the applicability of on-orbit calibration to different spectral bands and angles. To achieve the long-term, continuous, and high-precision absolute radiometric calibration of remote sensors, a wide-angle hyperspectral TOA reflectance model of the Libyan Desert was constructed based on spectral reflectance data, satellite overpass parameters, and atmospheric parameters from the Terra/Aqua and Earth Observation-1 (EO-1) satellites between 2003 and 2012. By means of angle fitting, viewing angle grouping, and spectral extension, the model is applicable for absolute radiometric calibration of the visible to short-wave infrared (SWIR) bands for sensors within viewing zenith angles of 65 degrees. To validate the accuracy and precision of the model, a total of 3120 long-term validations of model accuracy and 949 cross-validations with the Landsat 8 Operational Land Imager (OLI) and Suomi National Polar-Orbiting Partnership (NPP) Visible Infrared Imaging Radiometer Suite (VIIRS) satellite sensors between 2013 and 2020 were conducted. The results show that the TOA reflectance calculated by the model had a standard deviation (SD) of relative differences below 1.9% and a root-mean-square error (RMSE) below 0.8% when compared with observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Landsat 8 OLI. The SD of the relative differences and the RMSE were within 2.7% when predicting VIIRS data. Full article
(This article belongs to the Special Issue Remote Sensing Satellites Calibration and Validation)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The Libyan Desert in a satellite image. The overview map at the upper-right corner indicates the site’s location.</p>
Full article ">Figure 2
<p>Long-term clear-sky TOA reflectance in the red band of Landsat satellites.</p>
Full article ">Figure 3
<p>Long-term TOA reflectance of Aqua MODIS bands 1–7.</p>
Full article ">Figure 4
<p>Band coverage of the sensors.</p>
Full article ">Figure 5
<p>Daily nonuniformity distribution of the TOA reflectance and the iterative shrinking method.</p>
Full article ">Figure 6
<p>Distributions and sensitivity analyses of atmospheric parameters: AOD, water vapor, and ozone. (<b>a</b>) Distribution of AOD. (<b>b</b>) Relative difference analysis for AOD. (<b>c</b>) Distribution of water vapor. (<b>d</b>) Sensitivity analysis for water vapor. (<b>e</b>) Distribution of ozone. (<b>f</b>) Relative difference analysis for ozone. A positive relative difference corresponds to a brighter TOA reflectance.</p>
Full article ">Figure 7
<p>TOA reflectance as a function of the SZA and VZA in the blue band.</p>
Full article ">Figure 8
<p>Viewing angle distributions of Terra and Aqua.</p>
Full article ">Figure 9
<p>Hyperion hyperspectral TOA reflectance profiles.</p>
Full article ">Figure 10
<p>Spectral extension of the model results to Landsat 8 OLI bands <math display="inline"><semantics> <msubsup> <mi>ρ</mi> <mi>j</mi> <mrow> <mo>″</mo> </mrow> </msubsup> </semantics></math>.</p>
Full article ">Figure 11
<p>Spectral extension process.</p>
Full article ">Figure 12
<p>Observed and predicted TOA reflectance data for Terra MODIS.</p>
Full article ">Figure 13
<p>Relative differences between predictions and observations for Terra MODIS.</p>
Full article ">Figure 14
<p>Relative differences between predictions and observations for Aqua MODIS.</p>
Full article ">Figure 15
<p>Relative differences between predictions and observations for Landsat 8 OLI.</p>
Full article ">Figure 16
<p>Relative differences between predictions and observations for Suomi NPP VIIRS.</p>
Full article ">Figure 17
<p>Statistics of the model results. The positions of the circles along the y-axis represent the mean values of the relative differences. The error bars have a width of <math display="inline"><semantics> <mrow> <mo>±</mo> <mn>1</mn> </mrow> </semantics></math> SD around the circles.</p>
Full article ">
15 pages, 4657 KiB  
Article
Estimating Cotton Yield in the Brazilian Cerrado Using Linear Regression Models from MODIS Vegetation Index Time Series
by Daniel A. B. de Siqueira, Carlos M. P. Vaz, Flávio S. da Silva, Ednaldo J. Ferreira, Eduardo A. Speranza, Júlio C. Franchini, Rafael Galbieri, Jean L. Belot, Márcio de Souza, Fabiano J. Perina and Sérgio das Chagas
AgriEngineering 2024, 6(2), 947-961; https://doi.org/10.3390/agriengineering6020054 - 9 Apr 2024
Cited by 1 | Viewed by 1193
Abstract
Satellite remote sensing data expedite crop yield estimation, offering valuable insights for farmers’ decision making. Recent forecasting methods, particularly those utilizing machine learning algorithms like Random Forest and Artificial Neural Networks, show promise. However, challenges such as validation performances, large volume of data, [...] Read more.
Satellite remote sensing data expedite crop yield estimation, offering valuable insights for farmers’ decision making. Recent forecasting methods, particularly those utilizing machine learning algorithms like Random Forest and Artificial Neural Networks, show promise. However, challenges such as validation performances, large volume of data, and the inherent complexity and inexplicability of these models hinder their widespread adoption. This paper presents a simpler approach, employing linear regression models fitted from vegetation indices (VIs) extracted from MODIS sensor data on the Terra and Aqua satellites. The aim is to forecast cotton yields in key areas of the Brazilian Cerrado. Using data from 281 commercial production plots, models were trained (167 plots) and tested (114 plots), relating seed cotton yield to nine commonly used VIs averaged over 15-day intervals. Among the evaluated VIs, Enhanced Vegetation Index (EVI) and Triangular Vegetation Index (TVI) exhibited the lowest root mean square errors (RMSE) and the highest determination coefficients (R2). Optimal periods for in-season yield prediction fell between 90 and 105 to 135 and 150 days after sowing (DAS), corresponding to key phenological phases such as boll development, open boll, and fiber maturation, with the lowest RMSE of about 750 kg ha−1 and R2 of 0.70. The best forecasts for early crop stages were provided by models at the peaks (maximum value of the VI time series) for EVI and TVI, which occurred around 80–90 DAS. The proposed approach makes the yield predictability more inferable along the crop time series just by providing sowing dates, contour maps, and their respective VIs. Full article
Show Figures

Figure 1

Figure 1
<p>Spatial location of the 281 plots cultivated with cotton in the states of Mato Grosso (MT), Goiás (GO), and Bahia (BA) (typical Brazilian Cerrado biome) used to train (yellow dots) and test the models (blue dots) and three cotton plots zoomed (Plots 1, 2, and 3).</p>
Full article ">Figure 2
<p>MODIS product MCD43A4 V6.1: (<b>a</b>) NIR images obtained for two commercial cotton production plots; (<b>b</b>) original extracted time series (daily data) for red, NIR, blue, green, and SWIR spectral bands; (<b>c</b>) time series in 15-day interval averages. Plot on left presented low cotton yield on right high cotton yield.</p>
Full article ">Figure 3
<p>Average time series of reflectance for seven MODIS spectral bands for five cotton yield classes (<b>a</b>) and the linear correlation coefficients (r) between cotton yield and reflectance for the different 15-day days after sowing (DAS) intervals in the training dataset (167 plots) (<b>b</b>).</p>
Full article ">Figure 4
<p>Time series of NDVI, EVI, SAVI, and TVI from MODIS sensor for five cotton yield classes (<b>a</b>) and determination coefficients (R<sup>2</sup>) of linear correlations between cotton yield and 9 vegetation indices (<a href="#agriengineering-06-00054-t002" class="html-table">Table 2</a>) (<b>b</b>) for the different 15-day days after sowing (DAS) classes for the training dataset (167 plots).</p>
Full article ">Figure 5
<p>Linear correlations between TVI and seed cotton for different 15-day DAS for the training dataset, from 75–90 to 180–195 DAS, 75–195 DAS (averaged period) and the time series peak value. Linear equations shown in <a href="#agriengineering-06-00054-t003" class="html-table">Table 3</a>.</p>
Full article ">Figure 6
<p>Comparisons between observed and predicted seed cotton yield for the testing dataset (114 plots) for different 15-day DAS (75–90 to 180–195 DAS), 75–195 DAS (averaged period), the time series peak value, and their respective root mean square errors (RMSE), obtained using TVI linear regression models (<a href="#agriengineering-06-00054-t003" class="html-table">Table 3</a>).</p>
Full article ">Figure 7
<p>Root mean square error (RMSE) (<b>a</b>) and linear determination coefficients (R<sup>2</sup>) (<b>b</b>) between predicted and observed yield for the 114 plots of the testing dataset for different 15-day DAS; comparison of RMSE (<b>c</b>) and predicted yield (<b>d</b>) for 105–120 DAS, 75–195 DAS, and peak equations for TVI, EVI, SAVI, and NDVI. Dotted horizontal line in (<b>d</b>) represents the average observed yield for the 114 plots.</p>
Full article ">Figure 8
<p>Relationships between the residual error (Yield<sub>estimated</sub> − Yield<sub>observed</sub>; estimated by TVI peak model) and observed cotton yield (<b>a</b>); observed yield and sawing day from 1 November (<b>b</b>); observed yield and cotton cycle (<b>c</b>); and correlation coefficients (r) between observed yield and monthly accumulated precipitation, monthly averaged minimum (Tmin) and maximum (Tmax) temperatures, and Tmax − Tmin (<b>d</b>) for month 1 to 6 after sowing.</p>
Full article ">
Back to TopTop