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 (119)

Search Parameters:
Keywords = intercalibration

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
28 pages, 4077 KiB  
Article
Inter-Sensor Level 1 Radiometric Comparisons Using Deep Convective Clouds
by Louis Rivoire, Sébastien Clerc, Bahjat Alhammoud, Frédéric Romand and Nicolas Lamquin
Remote Sens. 2024, 16(23), 4445; https://doi.org/10.3390/rs16234445 - 27 Nov 2024
Viewed by 238
Abstract
To evaluate the radiometric performance of top-of-atmosphere reflectance images, Deep Convective Clouds (DCCs) can be used as temporally, spatially and spectrally stable targets. The DCCs method has been developed more than 20 years ago and applied recently to Sentinel-2 and Sentinel-3 sensors. In [...] Read more.
To evaluate the radiometric performance of top-of-atmosphere reflectance images, Deep Convective Clouds (DCCs) can be used as temporally, spatially and spectrally stable targets. The DCCs method has been developed more than 20 years ago and applied recently to Sentinel-2 and Sentinel-3 sensors. In this paper, among other developments, we built a new methodology upon those existing by using the bootstrap method and spectral band adjustment factors computed with the Hyper-Spectral Imager (HSI) from the Environmental Mapping and Analysis Program (EnMAP). This methodology is applied to the two Multi-Spectral Imager (MSI) instruments onboard Sentinel-2A and 2B, but also the two Operational Land Imager (OLI) instruments onboard Landsat 8 and 9, from visible wavelength at 442 nm to shortwave-infrared at 2200 nm, using images with a ground resolution spanning from 10 m to 60 m. The results demonstrate the good inter-calibration of MSI units A and B, which are within one percent of relative difference on average between January 2022 and June 2024 for all visible, near-infrared and shortwave-infrared bands, except for the band at 1375 nm for which saturation prevents the use of the method. Similarly, OLI and OLI-2 are found to have a relative difference on the same period lower than one percent for all 30 m resolution bands. Evaluation of the relative difference between the MSI sensors and the OLI sensors with the DCCs method gives values lower than three percent. Finally, these validation results are compared to those obtained with Pseudo-Invariant Calibration Sites (PICSs) over Libya-4: an agreement better than two percent is found between the DCCs and PICSs methods. Full article
Show Figures

Figure 1

Figure 1
<p>DCC detection steps for band B02 at 492 nm in product S2A_MSIL1C_20240505T144731_ N0510_R139_T20NKG_20240505T181400. (<b>a</b>) Raw Sentinel-2A TOA reflectance at 10 m resolution. (<b>b</b>) Subsampled TOA reflectance at 60 m resolution. (<b>c</b>) Subsampled TOA reflectance with detection thresholds applied using bands B8A and B10. (<b>d</b>) TOA reflectance with small DCC clusters removed. (<b>e</b>) TOA reflectance with morphological dilation applied to the DCC mask. (<b>f</b>) Top-Of-DCC reflectance after atmospheric correction has been applied.</p>
Full article ">Figure 2
<p>MSI-A histogram of a DCC in band 9 at 945 nm for all detectors in product S2A_ MSIL1C_20240505T144731_N0510_R139_T20NKG_20240505T181400.</p>
Full article ">Figure 3
<p>Sum of all MSI-A DCC histograms obtained in May 2024 along with the skewed Gaussian distribution fit and the associated first point of inflexion, mode and second point of inflexion.</p>
Full article ">Figure 4
<p>Distributions obtained with the bootstrap method (N equals 20 for readability) applied to MSI-A DCC histograms available in May 2024, with the second point of inflexion of each distribution.</p>
Full article ">Figure 5
<p>Example of DCC observed by EnMAP at 468 nm and used in this study. Product identifier is ENMAP-HSI-L1CDT0000060165_12-2024-02-01T08:30:08.938.</p>
Full article ">Figure 6
<p>Spectral response functions of Sentinel-2 MSI-A, MSI-B, Landsat 8 OLI and Lansat 9 OLI-2 band 3 at 560 nm at their highest resolution available and sampled at the resolution of EnMAP bands. (<b>a</b>) Sentinel-2 MSI-A and MSI-B. (<b>b</b>) Landsat 8 OLI and Lansat 9 OLI-2.</p>
Full article ">Figure 7
<p>Mean DCC spectrum measured with 15 EnMAP products, with Sentinel-2A bands overimposed. Error bars are a measure of the uncertainty in terms of the standard deviation of the mean.</p>
Full article ">Figure 8
<p>Location of products containing DCC pixels for Sentinel-2 units A and B, Landsat 8 and Landsat 9, in August 2023.</p>
Full article ">Figure 9
<p>Temporal evolution of the number of products containing DCC pixels for Sentinel-2 units A and B, Landsat 8 and Landsat 9, between January 2022 and June 2024.</p>
Full article ">Figure 10
<p>Temporal evolution of the reflectance indicators for Sentinel-2 units A and B, Landsat 8 and Landsat 9, between January 2022 and June 2024, for bands centered at 492 nm. (<b>a</b>) represents the temporal evolution of the absolute reflectance indicators while (<b>b</b>) represents the temporal evolution of the relative differences of a reflectance indicator of a sensor with that of Sentinel-2A. SBAFs are applied in (<b>a</b>,<b>b</b>); in (<b>a</b>) SBAFs are applied taking Sentinel-2A as reference. Error bars are a measure of the uncertainty in terms of the standard deviation of the mean. Dashed lines represent one and three percent relative difference with the reference.</p>
Full article ">Figure 11
<p>Same as <a href="#remotesensing-16-04445-f010" class="html-fig">Figure 10</a> for bands centred at 864 nm.</p>
Full article ">Figure 12
<p>Same as <a href="#remotesensing-16-04445-f010" class="html-fig">Figure 10</a> for bands centred at 1610 nm.</p>
Full article ">Figure 13
<p>Sentinel-2A image and associated histogram of a DCC over Barbados in band 10 on the 01/07/2024 (product S2A_MSIL1C_20240701T143751_N0510_R096_T20PRV_20240701T175847). Saturation over the DCC can be seen on the right-hand part of (<b>a</b>). (<b>a</b>) Sentinel-2A image. (<b>b</b>) Histogram of reflectances.</p>
Full article ">Figure 14
<p>Sentinel-2B over Sentinel-2A mean relative differences for the thirteen Sentinel-2 bands, averaged between January 2022 and June 2024; SBAFs are applied. Relative difference at 1375 nm is impacted by saturation. Error bars are a measure of the statistical dispersion of the relative difference over the averaging period. Dashed lines represent one and three percent relative difference with the reference.</p>
Full article ">Figure 15
<p>Landsat 9 over Landsat 8 mean relative differences for eight common bands, averaged between January 2022 and June 2024; SBAFs are applied. Error bars are a measure of the statistical dispersion of the relative difference over the averaging period. Dashed lines represent one and three percent relative difference with the reference.</p>
Full article ">Figure 16
<p>Sentinel-2B, Landsat 8 and Landsat 9 over Sentinel-2A mean relative differences for eight common bands, averaged between January 2022 and June 2024; SBAFs are applied. Relative differences at 1375 nm are impacted by MSI saturation. Error bars are a measure of the statistical dispersion of the relative difference over the averaging period. Dashed lines represent one and three percent relative difference with the reference.</p>
Full article ">Figure 17
<p>Mean Sentinel-2B, Landsat 8 and Landsat 9 over Sentinel-2A relative differences for bands at 442 nm, 492 nm, 560 nm, 665 nm and 864 nm, averaged between January 2022 and December 2023, for the PICSs method and the DCCs method. Error bars are a measure of the statistical dispersion of the relative difference over the averaging period. Dashed lines represent one and three percent relative difference with the reference.</p>
Full article ">Figure 18
<p>Relative differences of Sentinel-2B reflectance indicator over that of Sentinel-2A, for four geographical zones corresponding to strips of latitude and four months spanning over 2023. SBAFs are applied. Error bars are a measure of the uncertainty in terms of the standard deviation of the mean. Dashed lines represent one and three percent relative difference with the reference.</p>
Full article ">Figure A1
<p>Same as <a href="#remotesensing-16-04445-f010" class="html-fig">Figure 10</a> for bands centred at 442 nm.</p>
Full article ">Figure A2
<p>Same as <a href="#remotesensing-16-04445-f010" class="html-fig">Figure 10</a> for bands centred at 560 nm.</p>
Full article ">Figure A3
<p>Same as <a href="#remotesensing-16-04445-f010" class="html-fig">Figure 10</a> for bands centred at 665 nm.</p>
Full article ">Figure A4
<p>Same as <a href="#remotesensing-16-04445-f010" class="html-fig">Figure 10</a> for bands centred at 704 nm, with only Sentinel-2A and Sentinel-2B.</p>
Full article ">Figure A5
<p>Same as <a href="#remotesensing-16-04445-f010" class="html-fig">Figure 10</a> for bands centred at 740 nm, with only Sentinel-2A and Sentinel-2B.</p>
Full article ">Figure A6
<p>Same as <a href="#remotesensing-16-04445-f010" class="html-fig">Figure 10</a> for bands centred at 780 nm, with only Sentinel-2A and Sentinel-2B.</p>
Full article ">Figure A7
<p>Same as <a href="#remotesensing-16-04445-f010" class="html-fig">Figure 10</a> for bands centred at 833 nm, with only Sentinel-2A and Sentinel-2B.</p>
Full article ">Figure A8
<p>Same as <a href="#remotesensing-16-04445-f010" class="html-fig">Figure 10</a> for bands centred at 945 nm, with only Sentinel-2A and Sentinel-2B.</p>
Full article ">Figure A9
<p>Same as <a href="#remotesensing-16-04445-f010" class="html-fig">Figure 10</a> for bands centred at 1375 nm; y-axis scale increased to view error bars.</p>
Full article ">Figure A10
<p>Same as <a href="#remotesensing-16-04445-f010" class="html-fig">Figure 10</a> for bands centred at 2200 nm.</p>
Full article ">
35 pages, 16179 KiB  
Article
Vegetative Index Intercalibration Between PlanetScope and Sentinel-2 Through a SkySat Classification in the Context of “Riserva San Massimo” Rice Farm in Northern Italy
by Christian Massimiliano Baldin and Vittorio Marco Casella
Remote Sens. 2024, 16(21), 3921; https://doi.org/10.3390/rs16213921 - 22 Oct 2024
Viewed by 1992
Abstract
Rice farming in Italy accounts for about 50% of the EU’s rice area and production. Precision agriculture has entered the scene to enhance sustainability, cut pollution, and ensure food security. Various studies have used remote sensing tools like satellites and drones for multispectral [...] Read more.
Rice farming in Italy accounts for about 50% of the EU’s rice area and production. Precision agriculture has entered the scene to enhance sustainability, cut pollution, and ensure food security. Various studies have used remote sensing tools like satellites and drones for multispectral imaging. While Sentinel-2 is highly regarded for precision agriculture, it falls short for specific applications, like at the “Riserva San Massimo” (Gropello Cairoli, Lombardia, Northern Italy) rice farm, where irregularly shaped crops need higher resolution and frequent revisits to deal with cloud cover. A prior study that compared Sentinel-2 and the higher-resolution PlanetScope constellation for vegetative indices found a seasonal miscalibration in the Normalized Difference Vegetation Index (NDVI) and in the Normalized Difference Red Edge Index (NDRE). Dr. Agr. G.N. Rognoni, a seasoned agronomist working with this farm, stresses the importance of studying the radiometric intercalibration between the PlanetScope and Sentinel-2 vegetative indices to leverage the knowledge gained from Sentinel-2 for him to apply variable rate application (VRA). A high-resolution SkySat image, taken almost simultaneously with a pair of Sentinel-2 and PlanetScope images, offered a chance to examine if the irregular distribution of vegetation and barren land within rice fields might be a factor in the observed miscalibration. Using an unsupervised pixel-based image classification technique on SkySat imagery, it is feasible to split rice into two subclasses and intercalibrate them separately. The results indicated that combining histograms and agronomists’ expertise could confirm SkySat classification. Moreover, the uneven spatial distribution of rice does not affect the seasonal miscalibration object of past studies, which can be adjusted using the methods described here, even with images taken four days apart: the first method emphasizes accuracy using linear regression, histogram shifting, and histogram matching; whereas the second method is faster and utilizes only histogram matching. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Riserva San Massimo rice farm crops over 3 June 2021. SkySat image—EPSG: 4326.</p>
Full article ">Figure 2
<p>Process flowchart.</p>
Full article ">Figure 3
<p>SkySat unsupervised classification with vegetation in green and barren land in brown.</p>
Full article ">Figure 4
<p>SkySat image compared with the masks produced in MATLAB.</p>
Full article ">Figure 5
<p>Method 1 for NDVI.</p>
Full article ">Figure 6
<p>Method 1 for NDRE.</p>
Full article ">Figure 7
<p>Method 2 for NDVI—LR after HM is not useful for distribution and statistics.</p>
Full article ">Figure 8
<p>Method 2 for NDRE—LR after HM is not useful for distribution and statistics.</p>
Full article ">Figure A1
<p>Method 1 for NDVI—linear regression for full crops.</p>
Full article ">Figure A2
<p>Method 1 for NDVI—linear regression for vegetation subclass.</p>
Full article ">Figure A3
<p>Method 1 for NDVI—linear regression for barren land subclass.</p>
Full article ">Figure A4
<p>Method 1 for NDRE—linear regression for full crops.</p>
Full article ">Figure A5
<p>Method 1 for NDRE—linear regression for vegetation subclass.</p>
Full article ">Figure A6
<p>Method 1 for NDRE—linear regression for barren land subclass.</p>
Full article ">Figure A7
<p>Method 2 for NDVI—linear regression for full crops.</p>
Full article ">Figure A8
<p>Method 2 for NDVI—linear regression for vegetation subclass.</p>
Full article ">Figure A9
<p>Method 2 for NDVI—linear regression for barren land subclass.</p>
Full article ">Figure A10
<p>Method 2 for NDRE—linear regression for full crops.</p>
Full article ">Figure A11
<p>Method 2 for NDRE—linear regression for vegetation subclass.</p>
Full article ">Figure A12
<p>Method 2 for NDRE—linear regression for barren land subclass.</p>
Full article ">
14 pages, 7034 KiB  
Article
Macrophytes as Key Element to Determine Ecological Quality Changes in Transitional Water Systems: The Venice Lagoon as Study Case
by Adriano Sfriso, Alessandro Buosi, Yari Tomio, Giulia Silan, Marion Adelheid Wolf, Katia Sciuto and Andrea Augusto Sfriso
Environments 2024, 11(9), 209; https://doi.org/10.3390/environments11090209 - 22 Sep 2024
Viewed by 856
Abstract
According to European Union guidelines, the assessment of the ecological status of Transitional Water Systems (TWSs) should be based on the monitoring of biological communities rather than physico-chemical parameters and pollutants. Macrophytes, including aquatic angiosperms and macroalgae, are organisms that respond more quickly [...] Read more.
According to European Union guidelines, the assessment of the ecological status of Transitional Water Systems (TWSs) should be based on the monitoring of biological communities rather than physico-chemical parameters and pollutants. Macrophytes, including aquatic angiosperms and macroalgae, are organisms that respond more quickly to environmental changes by varying the structure and biomass of their assemblages. There are several ecological indices based on macrophytes, among them the Macrophyte Quality Index (MaQI), which has been intercalibrated with water and sediment parameters, nutrient concentrations, and pollutants and is used to determine the ecological status of Italian TWSs. In the Venice Lagoon, it was applied to 87 stations, showing a significant score increase over the last ten years of monitoring (2011–2021) due to progressive lagoon environmental recovery. The dominant taxa assemblages, previously dominated by Ulvaceae, were replaced by species of higher ecological value, with an increase in the number and distribution of sensitive species, as well as the spread and cover of aquatic angiosperms. The rise in the Ecological Quality Ratio (EQR) determined by the MaQI confirms the key role of macrophyte monitoring in detecting environmental changes in TWSs. In fact, a simple check of the presence or absence of aquatic angiosperms and sensitive species is sufficient for an initial rapid assessment of the ecological status of these environments. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Map of the Venice Lagoon.</p>
Full article ">Figure 2
<p>Macroalgal Quality Index (MaQI) scheme to assess the ecological status of transitional water systems (TWSs). As indicated by the Water Framework Directive (WFD 2000/60/EC) the Ecological Quality Ratio (EQR) values are highlighted in red (Bad conditions), ochre (Poor conditions), yellow (Moderate conditions), green (Good conditions) and light blue (High conditions).</p>
Full article ">Figure 3
<p>Ecological Quality Ratio (EQR) changes determined by the Macrophyte Quality Index (MaQI) in the water bodies (WBs) of the Venice Lagoon from 2011 to 2021. Legend: DV = Dogà Valley; CV = Cavallino Valley; ZV = Zappa Valley; ER = Euryhaline Restricted; ENR = Euryhaline Non-Restricted; PR = Polyhaline Restricted; PNR = Polyhaline Non-Restricted.</p>
Full article ">Figure 4
<p>Changes in (<b>A</b>) the mean number of total macroalgal taxa per station; (<b>B</b>) the mean number of sensitive macroalgal taxa per station; (<b>C</b>) the total number of stations colonized by aquatic angiosperms; (<b>D</b>) the mean percentage of aquatic angiosperm cover in the total stations during the 4 sampling years.</p>
Full article ">Figure 5
<p>(<b>A</b>) Number of colonized stations by each angiosperm species; (<b>B</b>) total mean cover of each angiosperm in the 87 stations between 2011 and 2021.</p>
Full article ">Figure 6
<p>(<b>A</b>) Mean values of Reactive Phosphorus (RP); (<b>B</b>) Dissolved Inorganic Nitrogen (DIN = sum of ammonium, nitrite, nitrate); (<b>C</b>) Total Chlorophyll-a (Chl-a); and (<b>D</b>) Total Suspended Solids (TSSs) in the 87 stations between 2011 and 2021.</p>
Full article ">
27 pages, 9855 KiB  
Article
Inter-Calibration of Passive Microwave Satellite Brightness Temperature Observations between FY-3D/MWRI and GCOM-W1/AMSR2
by Zuomin Xu, Ruijing Sun, Shuang Wu, Jiali Shao and Jie Chen
Remote Sens. 2024, 16(2), 424; https://doi.org/10.3390/rs16020424 - 22 Jan 2024
Viewed by 1206
Abstract
Microwave sensors possess the capacity to effectively penetrate through clouds and fog and are widely used in obtaining soil moisture, atmospheric water vapor, and surface temperature measurements. Long time-series datasets play a pivotal role in climate change studies. Unfortunately, the lifespan of operational [...] Read more.
Microwave sensors possess the capacity to effectively penetrate through clouds and fog and are widely used in obtaining soil moisture, atmospheric water vapor, and surface temperature measurements. Long time-series datasets play a pivotal role in climate change studies. Unfortunately, the lifespan of operational satellites often falls short of the needs of these extensive datasets. Hence, comparing and cross-calibrating sensors with similar configurations is paramount. The Microwave Radiation Imager (MWRI) onboard Fengyun-3D (FY-3D) is the latest generation of satellite-based microwave remote sensing instruments in China, and its data quality and application prospects have attracted widespread attention. To comprehensively assess the data quality of MWRI, a comparison of the orbital brightness temperature (TB) data between FY-3D/MWRI and Global Change Observation Mission 1st-Water (GCOM-W1)/Advanced Microwave Scanning Radiometer 2 (AMSR2) is conducted, and then a calibration model is established. The results indicate a strong correlation between the two sensors, with a correlation coefficient exceeding 0.9 across all channels. The mean bias ranges from −1.5 K to 0.15 K. Notably, the bias of vertical polarization is more pronounced than that of horizontal polarization. The TB distribution patterns and temporal evolutions are highly consistent for both sensors, particularly under snow and ice. The small intercepts and close-to-1 slopes obtained during calibration further demonstrate the minor data differences between the two sensors. However, the calibration process effectively reduces the existing errors, and the calibrated FY-3D/MWRI TB data are closer to GCOM-W1/AMSR2, with a mean bias approximately equal to 0 K and a correlation coefficient exceeding 0.99. The excellent consistency of the TB data between the two sensors provides a vital data basis for retrieving surface parameters and establishing long time-series datasets. Full article
(This article belongs to the Special Issue Remote Sensing Satellites Calibration and Validation)
Show Figures

Figure 1

Figure 1
<p>Spatial collocation location for FY-3D/MWRI and GCOM-W1/AMSR2 based on spatial, temporal, and geometry collocation criteria, as well as filtering with globally homogeneous pixels for the period of January–December 2022, on the 15th day and at the end of the month, for a total of 24 days (different days are shown in different colors).</p>
Full article ">Figure 2
<p>Differences in annual mean brightness temperature between FY−3D/MWRI and GCOM−W1/AMSR2 over the global land area: (<b>a</b>) 10.65 GHz (H; horizontal), (<b>b</b>) 10.65 GHz (V; vertical), (<b>c</b>) 18.7 GHz (H), (<b>d</b>) 18.7 GHz (V), (<b>e</b>) 23.8 GHz (H), (<b>f</b>) 23.8 GHz (V), (<b>g</b>) 36.5 GHz (H), (<b>h</b>) 36.5 GHz (V), (<b>i</b>) 89 GHz (H), and (<b>j</b>) 89 GHz (V).</p>
Full article ">Figure 2 Cont.
<p>Differences in annual mean brightness temperature between FY−3D/MWRI and GCOM−W1/AMSR2 over the global land area: (<b>a</b>) 10.65 GHz (H; horizontal), (<b>b</b>) 10.65 GHz (V; vertical), (<b>c</b>) 18.7 GHz (H), (<b>d</b>) 18.7 GHz (V), (<b>e</b>) 23.8 GHz (H), (<b>f</b>) 23.8 GHz (V), (<b>g</b>) 36.5 GHz (H), (<b>h</b>) 36.5 GHz (V), (<b>i</b>) 89 GHz (H), and (<b>j</b>) 89 GHz (V).</p>
Full article ">Figure 2 Cont.
<p>Differences in annual mean brightness temperature between FY−3D/MWRI and GCOM−W1/AMSR2 over the global land area: (<b>a</b>) 10.65 GHz (H; horizontal), (<b>b</b>) 10.65 GHz (V; vertical), (<b>c</b>) 18.7 GHz (H), (<b>d</b>) 18.7 GHz (V), (<b>e</b>) 23.8 GHz (H), (<b>f</b>) 23.8 GHz (V), (<b>g</b>) 36.5 GHz (H), (<b>h</b>) 36.5 GHz (V), (<b>i</b>) 89 GHz (H), and (<b>j</b>) 89 GHz (V).</p>
Full article ">Figure 2 Cont.
<p>Differences in annual mean brightness temperature between FY−3D/MWRI and GCOM−W1/AMSR2 over the global land area: (<b>a</b>) 10.65 GHz (H; horizontal), (<b>b</b>) 10.65 GHz (V; vertical), (<b>c</b>) 18.7 GHz (H), (<b>d</b>) 18.7 GHz (V), (<b>e</b>) 23.8 GHz (H), (<b>f</b>) 23.8 GHz (V), (<b>g</b>) 36.5 GHz (H), (<b>h</b>) 36.5 GHz (V), (<b>i</b>) 89 GHz (H), and (<b>j</b>) 89 GHz (V).</p>
Full article ">Figure 2 Cont.
<p>Differences in annual mean brightness temperature between FY−3D/MWRI and GCOM−W1/AMSR2 over the global land area: (<b>a</b>) 10.65 GHz (H; horizontal), (<b>b</b>) 10.65 GHz (V; vertical), (<b>c</b>) 18.7 GHz (H), (<b>d</b>) 18.7 GHz (V), (<b>e</b>) 23.8 GHz (H), (<b>f</b>) 23.8 GHz (V), (<b>g</b>) 36.5 GHz (H), (<b>h</b>) 36.5 GHz (V), (<b>i</b>) 89 GHz (H), and (<b>j</b>) 89 GHz (V).</p>
Full article ">Figure 3
<p>Comparison of TB in each channel of AMSR2 and MWRI. The different colors indicate the aggregation of data points. The subfigures represent the following frequency and polarization combinations: (<b>a</b>) 10.65 GHz (H; horizontal), (<b>b</b>) 10.65 GHz (V; vertical), (<b>c</b>) 18.7 GHz (H), (<b>d</b>) 18.7 GHz (V), (<b>e</b>) 23.8 GHz (H), (<b>f</b>) 23.8 GHz (V), (<b>g</b>) 36.5 GHz (H), (<b>h</b>) 36.5 GHz (V), (<b>i</b>) 89 GHz (H), and (<b>j</b>) 89 GHz (V).</p>
Full article ">Figure 3 Cont.
<p>Comparison of TB in each channel of AMSR2 and MWRI. The different colors indicate the aggregation of data points. The subfigures represent the following frequency and polarization combinations: (<b>a</b>) 10.65 GHz (H; horizontal), (<b>b</b>) 10.65 GHz (V; vertical), (<b>c</b>) 18.7 GHz (H), (<b>d</b>) 18.7 GHz (V), (<b>e</b>) 23.8 GHz (H), (<b>f</b>) 23.8 GHz (V), (<b>g</b>) 36.5 GHz (H), (<b>h</b>) 36.5 GHz (V), (<b>i</b>) 89 GHz (H), and (<b>j</b>) 89 GHz (V).</p>
Full article ">Figure 4
<p>The TB frequency statistics for each channel of AMSR2 (blue line) and MWRI (red line): (<b>a</b>–<b>j</b>) present the frequency statistical results for the ten channels of FY−3D/MWRI and GCOM−W1/AMSR2 from 10.65–89 GHz. The Min (maximum), Min (minimum), Mean (mean), and Std (standard deviation) of TB for two sensors are shown on the left-hand side.</p>
Full article ">Figure 4 Cont.
<p>The TB frequency statistics for each channel of AMSR2 (blue line) and MWRI (red line): (<b>a</b>–<b>j</b>) present the frequency statistical results for the ten channels of FY−3D/MWRI and GCOM−W1/AMSR2 from 10.65–89 GHz. The Min (maximum), Min (minimum), Mean (mean), and Std (standard deviation) of TB for two sensors are shown on the left-hand side.</p>
Full article ">Figure 5
<p>The total number of matching pixels, correlation coefficients (R<sup>2</sup>), the root mean square error (RMSE), and bias values for each channel in different TB segments. All the TB values are divided into 12 segments with a TB interval of 20 K.</p>
Full article ">Figure 5 Cont.
<p>The total number of matching pixels, correlation coefficients (R<sup>2</sup>), the root mean square error (RMSE), and bias values for each channel in different TB segments. All the TB values are divided into 12 segments with a TB interval of 20 K.</p>
Full article ">Figure 5 Cont.
<p>The total number of matching pixels, correlation coefficients (R<sup>2</sup>), the root mean square error (RMSE), and bias values for each channel in different TB segments. All the TB values are divided into 12 segments with a TB interval of 20 K.</p>
Full article ">Figure 6
<p>Land cover map with homogeneous pixels on a 25 km EASE-GRID projection grid.</p>
Full article ">Figure 7
<p>Distribution of TB differences (TB<sub>MWRI</sub> − TB<sub>AMSR2</sub>) and RMSEs between AMSR2 and MWRI data for horizontal and vertical polarizations across each channel for different land cover types. The vertical bars in the left box plot represent the upper to lower quartiles of the data, with the narrowing in the middle indicating the median value. The ends of the vertical dotted lines represent the maximum and minimum values. Points outside the vertical dotted line represent outliers.</p>
Full article ">Figure 7 Cont.
<p>Distribution of TB differences (TB<sub>MWRI</sub> − TB<sub>AMSR2</sub>) and RMSEs between AMSR2 and MWRI data for horizontal and vertical polarizations across each channel for different land cover types. The vertical bars in the left box plot represent the upper to lower quartiles of the data, with the narrowing in the middle indicating the median value. The ends of the vertical dotted lines represent the maximum and minimum values. Points outside the vertical dotted line represent outliers.</p>
Full article ">Figure 7 Cont.
<p>Distribution of TB differences (TB<sub>MWRI</sub> − TB<sub>AMSR2</sub>) and RMSEs between AMSR2 and MWRI data for horizontal and vertical polarizations across each channel for different land cover types. The vertical bars in the left box plot represent the upper to lower quartiles of the data, with the narrowing in the middle indicating the median value. The ends of the vertical dotted lines represent the maximum and minimum values. Points outside the vertical dotted line represent outliers.</p>
Full article ">Figure 8
<p>RMSE and bias distributions of pre−calibration MWRI data and after−calibration MWRI data with respect to AMSR2 data for each channel at different TB intervals. All the TB values are also divided into 12 segments with a TB interval of 20 K.</p>
Full article ">Figure 8 Cont.
<p>RMSE and bias distributions of pre−calibration MWRI data and after−calibration MWRI data with respect to AMSR2 data for each channel at different TB intervals. All the TB values are also divided into 12 segments with a TB interval of 20 K.</p>
Full article ">
12 pages, 5526 KiB  
Article
Adaptation of the European Fish Index (EFI+) to Include the Alien Fish Pressure
by Enric Aparicio, Carles Alcaraz, Rafel Rocaspana, Quim Pou-Rovira and Emili García-Berthou
Fishes 2024, 9(1), 13; https://doi.org/10.3390/fishes9010013 - 29 Dec 2023
Cited by 1 | Viewed by 1893
Abstract
The European Fish Index EFI+ is the only fish-based multimetric index for the assessment of the ecological status of running waters that is validated and thus applicable across most countries of the European Union. Metrics of the index rely on several attributes of [...] Read more.
The European Fish Index EFI+ is the only fish-based multimetric index for the assessment of the ecological status of running waters that is validated and thus applicable across most countries of the European Union. Metrics of the index rely on several attributes of the species present in the fish assemblage, irrespective of their native/alien status. The abundance of alien fish, together with other anthropogenic impacts, is one of the most important threats to the conservation of native fish and ecosystem health and is also an indicator of degraded stream conditions. Therefore, to improve the performance of the EFI+ in regions with high incidence of alien species, the EFI+ was adapted to include alien fish pressure as a new metric that reflects the number of alien species as well as the proportional abundance of alien individuals. The application of the adapted index (A-EFI+) is illustrated with data from several Iberian Mediterranean basins and showed similar or stronger correlations than the original EFI+ with anthropogenic pressure (land-use variables and alterations in hydrology and river morphology) and with other regional fish indices. EFI+ has been invaluable to intercalibrate fish indices across Europe, and A-EFI+ is similar but explicitly includes alien pressure, thus helping to provide a more comprehensive assessment of ecosystem health and to communicate it to society. Full article
(This article belongs to the Special Issue Biomonitoring and Conservation of Freshwater & Marine Fishes)
Show Figures

Figure 1

Figure 1
<p>Relationship between fish alien metrics, land-use variables and hydrological and morphological alterations. The panels above the diagonal show the Spearman rank correlation coefficients with significance level (*** <span class="html-italic">p</span> &lt; 0.001) and the panels below show the pairwise scatterplot with a smoothing curve (LOESS, red line). In the scatterplots, the <span class="html-italic">Y</span>-axis corresponds to the variable in the row diagonal and the <span class="html-italic">X</span>-axis to the column diagonal (e.g., the scatterplot on the bottom left has % alien individuals in the <span class="html-italic">Y</span>-axis and % forest cover in the <span class="html-italic">X</span>-axis).</p>
Full article ">Figure 2
<p>Relationship between the EFI+ and A-EFI+ indices (blue crosses and green circles, respectively) and the alien fish pressure (AFP) (average of % alien species and % alien individuals) in the study area. The simple regression lines are also shown (EFI+ = 0.803 − 0.435 AFP, <span class="html-italic">R</span><sup>2</sup><sub>adj</sub> = 0.465, <span class="html-italic">p</span> &lt; 0.001; A-EFI+ = 0.795 − 0.575 AFP, <span class="html-italic">R</span><sup>2</sup><sub>adj</sub> = 0.683, <span class="html-italic">p</span> &lt; 0.001).</p>
Full article ">Figure 3
<p>Relationship between the A-EFI+ and EFI+ indices and the effects of the alien fish pressure. The A-EFI+ index can be easily estimated from EFI+ with the following linear regression functions: A-EFI+ = −0.089 + 1.058 EFI+, <span class="html-italic">R</span><sup>2</sup><sub>adj</sub> = 0.940, <span class="html-italic">p</span> &lt; 0.001; A-EFI+ = 0.1308 + 0.827 EFI+ − 0.215 AFP, <span class="html-italic">R</span><sup>2</sup><sub>adj</sub> = 0.991, <span class="html-italic">p</span> &lt; 0.001.</p>
Full article ">Figure 4
<p>Relationship between fish indices and land-use variables. The panels above the diagonal show the Spearman rank correlation coefficients with significance level (*** <span class="html-italic">p</span> &lt; 0.001) and the panels below the pairwise scatterplot with a smoothing curve (LOESS, red line). In the scatterplots, the <span class="html-italic">Y</span>-axis corresponds to the variable in the row diagonal and the <span class="html-italic">X</span>-axis to the column diagonal (e.g., the scatterplot on the bottom left has A-EFI+ in the <span class="html-italic">Y</span>-axis and % forest cover in the <span class="html-italic">X</span>-axis).</p>
Full article ">Figure 5
<p>Relationship between biotic indices and land-use variables. The panels above the diagonal show the Spearman rank correlation coefficients with significance level (*** <span class="html-italic">p</span> &lt; 0.001) and the panels below the pairwise scatterplot with a smoothing curve (LOESS, red line). In the scatterplots, the <span class="html-italic">Y</span>-axis corresponds to the variable in the row diagonal and the <span class="html-italic">X</span>-axis to the column diagonal (e.g., the scatterplot on the bottom left has A-EFI+ in the <span class="html-italic">Y</span>-axis and % forest cover in the <span class="html-italic">X</span>-axis).</p>
Full article ">Figure 6
<p>Relationship between land-use variables and hydrological and morphological alterations with EFI+ and A-EFI+ indices. The panels above the diagonal show the Spearman rank correlation coefficients with significance level (*** <span class="html-italic">p</span> &lt; 0.001) and the panels below the pairwise scatterplot with a smoothing curve (LOESS, red line). In the scatterplots, the <span class="html-italic">Y</span>-axis corresponds to the variable in the row diagonal and the <span class="html-italic">X</span>-axis to the column diagonal (e.g., the scatterplot on the bottom left has A-EFI+ in the <span class="html-italic">Y</span>-axis and % forest cover in the <span class="html-italic">X</span>-axis).</p>
Full article ">
15 pages, 2602 KiB  
Article
Measurement of the Growth of the Main Commercial Rays (Raja clavata, Raja brachyura, Torpedo marmorata, Dipturus oxyrinchus) in European Waters Using Intercalibration Methods
by Andrea Bellodi, Pierluigi Carbonara, Kirsteen M. MacKenzie, Blondine Agus, Karen Bekaert, Eleanor S. I. Greenway, Maria C. Follesa, Manfredi Madia, Andrea Massaro, Michele Palmisano, Chiara Romano, Mauro Sinopoli, Francesca Ferragut-Perello and Kélig Mahé
Biology 2024, 13(1), 20; https://doi.org/10.3390/biology13010020 - 29 Dec 2023
Cited by 1 | Viewed by 1744
Abstract
The intercalibration of age readings represents a crucial step in the ageing procedure; the use of different sampling methods, structures, preparation techniques, and ageing criteria can significantly affect age and growth data. This study evaluated the precision and accuracy of ageing for the [...] Read more.
The intercalibration of age readings represents a crucial step in the ageing procedure; the use of different sampling methods, structures, preparation techniques, and ageing criteria can significantly affect age and growth data. This study evaluated the precision and accuracy of ageing for the most important North Atlantic (NA) and Mediterranean (M) ray species, Raja clavata, Raja brachyura, Torpedo marmorata, and Dipturus oxyrinchus, through exchange exercises carried out by readers from different laboratories. In addition, growth parameters were estimated from the obtained data. A total of 663 individual batoids were analysed. R. clavata and R. brachyura samples were obtained from both the NA and the M, while vertebral centra of T. marmorata and D. oxyrinchus were only available for the M. High reading variability was observed for all four evaluated species in terms of CV, APE, and PA. D. oxyrinchus and T. marmorata showed relatively slow growth and the von Bertalanffy model with fixed t0 and Gompertz’s model were, respectively, the most precise models for each of these species. In R. brachyura, females had a faster growth rate compared to combined sexes. The vbt0p proved the most precise model for describing growth in this species, and no statistical differences were found between the NO and the M. For R. clavata, the best-fitting model was the vbt0p for females and males in the NO and for females from the M, while the best-fitting model for males from the M and sexes combined for both areas was log.p. Distinct growth patterns were observed between the two study areas. Full article
(This article belongs to the Section Conservation Biology and Biodiversity)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Summary of the preparation methods used for the analysed species in the Atlantic Ocean and Mediterranean Sea. Dots and lines represent the structure interpretation by different readers.</p>
Full article ">Figure 2
<p>Coefficient of variation (CV), percentage of agreement (PA), and standard deviation values for each modal age of <span class="html-italic">R. clavata</span> (<b>A</b>), <span class="html-italic">R. brachyura</span> (<b>B</b>), <span class="html-italic">T. marmorata</span> (<b>C</b>), and <span class="html-italic">D. oxyrinchus</span> (<b>D</b>).</p>
Full article ">Figure 3
<p>Comparison of the estimated growth curves and the observed age-at-length data obtained for <span class="html-italic">R. clavata</span> (<b>left</b>) and <span class="html-italic">R. brachyura</span> (<b>right</b>) in the Atlantic Ocean and Mediterranean Sea.</p>
Full article ">Figure 4
<p>Age-at-length data for each species with the estimated growth curves from each best-fitting model: von Bertalanffy with constrained <span class="html-italic">t</span><sub>0</sub> (vbt0p) curves for <span class="html-italic">R. brachyura</span> and <span class="html-italic">D. oxyrinchus</span>, logistic curve (log.p) for <span class="html-italic">R. clavata</span>, and Gompertz’s curve (vbT1p) for <span class="html-italic">T. marmorata</span>.</p>
Full article ">
17 pages, 2997 KiB  
Article
Combining RadCalNet Sites for Radiometric Cross Calibration of Landsat 9 and Landsat 8 Operational Land Imagers (OLIs)
by Norvik Voskanian, Kurtis Thome, Brian N. Wenny, Mohammad H. Tahersima and Mehran Yarahmadi
Remote Sens. 2023, 15(24), 5752; https://doi.org/10.3390/rs15245752 - 15 Dec 2023
Cited by 5 | Viewed by 1512
Abstract
Combining images from multiple Earth Observing (EO) satellites increases the temporal resolution of the data, overcoming the limitations imposed by low revisit time and cloud coverage. However, this requires an intercalibration process to ensure that there is no radiometric difference in top-of-atmosphere (TOA) [...] Read more.
Combining images from multiple Earth Observing (EO) satellites increases the temporal resolution of the data, overcoming the limitations imposed by low revisit time and cloud coverage. However, this requires an intercalibration process to ensure that there is no radiometric difference in top-of-atmosphere (TOA) observations or to quantify any offset in the respective instruments. In addition, combining vicarious calibration processes to the intercalibration of instruments can provide a useful mechanism to validate and compare data from multiple sensors. The Radiometric Calibration Network (RadCalNet) provides automated surface and top-of-atmosphere reflectance data from multiple participating ground sites that can be used for instrument vicarious calibration. We present a comparative analysis of the Landsat 8 and Landsat 9 Operational Land Imagers (OLI) sensors and validate the data by comparing them to measurements from RadCalNet sites as a quantitative intercalibration approach. RadCalNet serves as a common reference for instrument radiometric calibration, providing SI-traceable TOA reflectance with its associated absolute uncertainties. This paper discusses the method of combining data from multiple sites and calculating the weighted average by comparing the TOA reflectance of the instruments and their associated uncertainties. The presented process provides a SI-traceable intercalibration methodology and quantifies the offset and uncertainty in the Landsat 8 and 9 OLI instruments, demonstrating that the two instruments are in good agreement with each other and the data can be reliably cross-correlated and used by the scientific community. Full article
Show Figures

Figure 1

Figure 1
<p>Relative spectral response of OLI sensor for Landsat 8 and Landsat 9 instruments.</p>
Full article ">Figure 2
<p>TOA reflectance (blue) and uncertainty in the TOA reflectance (orange) of a typical RVUS data set.</p>
Full article ">Figure 3
<p>Scatter plot of TOA reflectance and RadCalNet predicted TOA reflectance for Landsat 9 based on 2022 data set.</p>
Full article ">Figure 4
<p>Scatter plot of TOA reflectance and RadCalNet predicted TOA reflectance for Landsat 8 from 2013 to 2023.</p>
Full article ">Figure 5
<p>TOA reflectance ratios for all sites for Landsat 9 OLI band 4 for all dates in 2022.</p>
Full article ">Figure 6
<p>TOA reflectance ratios for all sites for Landsat 8 OLI band 4 for all dates in 2022.</p>
Full article ">Figure 7
<p>TOA reflectance ratios for all sites for Landsat 8 OLI band 4 for all dates since 2013.</p>
Full article ">Figure 8
<p>TOA reflectance ratio as a function of wavelength for bands 1–5 of Landsat 8 and Landsat 9 instruments based on weighted mean of the four RadCalNet sites. The dashed line represents the expected RadCalNet intercalibration uncertainty for an ideal instrument where <math display="inline"><semantics> <mrow> <mfrac> <mrow> <msub> <mrow> <mi>ρ</mi> </mrow> <mrow> <mi>T</mi> <mi>O</mi> <mi>A</mi> <mo>_</mo> <mi>O</mi> <mi>L</mi> <mi>I</mi> </mrow> </msub> </mrow> <mrow> <msub> <mrow> <mi>ρ</mi> </mrow> <mrow> <mi>T</mi> <mi>O</mi> <mi>A</mi> <mo>_</mo> <mi>R</mi> <mi>C</mi> <mi>N</mi> </mrow> </msub> </mrow> </mfrac> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 9
<p>Absolute uncertainties of TOA reflectance ratios as a function of TOA reflectance ratios for Landsat 9.</p>
Full article ">Figure 10
<p>Absolute uncertainties of TOA reflectance ratios as a function of TOA reflectance ratios for Landsat 8.</p>
Full article ">
17 pages, 5075 KiB  
Article
Validation of EMI-2 Radiometric Performance with TROPOMI over Dome C Site in Antarctica
by Jingming Su, Fuqi Si, Minjie Zhao, Haijin Zhou and Yan Hong
Remote Sens. 2023, 15(8), 2012; https://doi.org/10.3390/rs15082012 - 11 Apr 2023
Viewed by 1602
Abstract
(1) The Environmental Trace Gases Monitoring Instrument-2(EMI-2) is a high-quality spaceborne imaging spectrometer that launched in September 2021. To evaluate its radiometric calibration performance in-flight, the UV2 and VIS1 bands of EMI-2 were cross-calibrated by the corresponding bands (band3 and band4) of TROPOMI [...] Read more.
(1) The Environmental Trace Gases Monitoring Instrument-2(EMI-2) is a high-quality spaceborne imaging spectrometer that launched in September 2021. To evaluate its radiometric calibration performance in-flight, the UV2 and VIS1 bands of EMI-2 were cross-calibrated by the corresponding bands (band3 and band4) of TROPOMI over the pseudo-invariant calibration site Dome C. (2) After angle limitation and cloud filtering of the Earth radiance data measured by EMI-2 and TROPOMI over Dome C, the top of atmosphere (TOA) reflectance time series were calculated. The spectral adjustment factors (SAF) were derived from the solar spectrum measured by the sensor to minimize the uncertainties caused by the different spectral response functions (SRF) of sensors. In addition, a correction method based on the radiative transfer model (RTM) SCIATRAN was used to suppress unaccounted angular dependence of atmospheric scattering. The radiation performance of EMI-2 is evaluated using the TOA reflectance ratio of EMI-2 and TROPOMI, combining the SAF correction and RTM-based correction methods. (3) It was shown that the time series trending of the TOA reflectance ratio between EMI-2 measurements and TROPOMI demonstrate flat characteristics and strong correlation. The mean reflectance ratios range from 0.998 to 1.09. The standard deviation of the reflection ratio is less than 3%. For 328 nm, 335 nm, 340 nm, 460 nm, and 490 nm, the mean values are close to one, and the relative radiometric bias estimated through EMI-2 and TROPOMI intercalibration is less than 3%, and for other wavelengths, the biases are less than 6%, except for 416 nm, which behaves higher than 7%. The cross-calibration results show that the radiometric calibration of EMI-2 is within the relative accuracy requirement. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>EMI-2 optical layout.</p>
Full article ">Figure 2
<p>Geographical location of the calibration site.</p>
Full article ">Figure 3
<p>Flowchart for cross-calibration of EMI-2 to TROPOMI based on RTM method.</p>
Full article ">Figure 4
<p>Reflectance histograms of EMI-2 UV2 (340 nm) with and without anomaly removal.</p>
Full article ">Figure 5
<p>Polynomial relationship between the reflectance factor and SZA for EMI-2/TROPOMI data.</p>
Full article ">Figure 6
<p>TOA reflectance measured by EMI-2 on 8 February 2022 over Dome C.</p>
Full article ">Figure 7
<p>The irradiance spectra in ultraviolet (<b>a</b>) and visible (<b>b</b>) bands measured by EMI-2 (black) and TROPOMI (red) on 23 December 2021. The blue line is from SAO2010 solar irradiance reference spectrum and convolved with the ISRF of EMI-2.</p>
Full article ">Figure 8
<p>Comparison of spectral adjustment factors.</p>
Full article ">Figure 9
<p>Time series of atmospheric transmission correction factor between EMI-2 and TROPOMI for twelve bands.</p>
Full article ">Figure 9 Cont.
<p>Time series of atmospheric transmission correction factor between EMI-2 and TROPOMI for twelve bands.</p>
Full article ">Figure 10
<p>Mean and standard deviation series of atmospheric transmission correction factor of EMI-2 and TROPOMI in different months.</p>
Full article ">Figure 11
<p>Time series of reflectance ratios of EMI-2 to TROPOMI after spectral adjustment factor and RTM-based correction.</p>
Full article ">Figure 12
<p>The mean and standard deviation of reflectance ratio between EMI-2 and TROPOMI after spectral adjustment factor and RTM correction.</p>
Full article ">
24 pages, 11194 KiB  
Article
Characterization of the East—West Spatial Uniformity for GOES-16/17 ABI Bands Using the Moon
by Fangfang Yu, Xiangqian Wu, Xi Shao and Haifeng Qian
Remote Sens. 2023, 15(7), 1881; https://doi.org/10.3390/rs15071881 - 31 Mar 2023
Viewed by 1634
Abstract
The Advanced Baseline Imager (ABI) is the primary instrument onboard the NOAA Geostationary Operational Environmental Satellite-R Series (GOES-R) satellites, providing continuous weather imagery over the vast area in the Western Hemisphere. It is imperative to ensure consistent calibration accuracy within the instrument’s field [...] Read more.
The Advanced Baseline Imager (ABI) is the primary instrument onboard the NOAA Geostationary Operational Environmental Satellite-R Series (GOES-R) satellites, providing continuous weather imagery over the vast area in the Western Hemisphere. It is imperative to ensure consistent calibration accuracy within the instrument’s field of regard (FOR). This paper characterized the spatial uniformity in the east–west (EW) direction for the six ABI visible and near-infrared (VNIR) bands of the first two GOES-R satellites, GOES-16 (G16) and GOES-17 (G17), using a special collection of lunar chasing images during their post-launch testing and post-launch product testing (PLT/PLPT) periods. The EW response versus scan-angle (RVS) is examined with the normalized lunar irradiance ratios at varying scan angles combined from multiple lunar-chasing events. The impacts of straylight from the Earth were found in some of the B01–B03 lunar images. The straylight, including those scattered into the spacelook scenes near the polar regions and those leaked into space near the Moon, can cause RVS variation up to 1% for B01 and to a lesser magnitude for the other two bands. Straylight correction algorithms are applied for the accurate ABI lunar image irradiance calculation. After the corrections, the RVS variation is reduced to less than 0.3% for all the VNIR bands of both G16/17 in full-disk (FD) images. Results of this study also confirm that the Global Space-based Inter-Calibration System (GSICS) Implementation of the ROLO (GIRO) model has high relative accuracy for the ABI VNIR bands when the lunar images are collected within a relatively short time. The method described in this paper can be applied to validate the EW spatial uniformity for imagers on other geostationary satellites, including the recently launched GOES-18 and the future GOES-U satellites. Full article
Show Figures

Figure 1

Figure 1
<p>Illustration of ABI FOR (the circle in thin dotted line), with the Earth in the center and the Moon traversing along the red dashed line through part of ABI’s FOR.</p>
Full article ">Figure 2
<p>ABI field of regard (the outer circle in dashed gray line), Earth (the inner circle in solid black line), and the center locations of the MESO frames (dots with different colors) of the G16 (<b>left</b>) and G17 (<b>right</b>) lunar chasing events within the ABI FOR. G16 lunar chasing events: 11 February 2017 (in red), 12 February 2017 (in brown), 14 February 2017 (in pink), 4 December 2017 (in blue), and 12 July 2017 (in green); G17 lunar chasing events: 30 July 2018 (in red), 28 August 2018 (in brown), 29 August 2018 (in pink), and 20 October 2018 (in blue). Some of the G16 lunar images were not collected on 4 December 2017.</p>
Full article ">Figure 3
<p>(<b>a</b>) Locations of the G17 spacelook events (in red dots) for the timeline between 20:02 and 20:17 UTC on 30 July 2018. The spacelook positions are located east to the Earth for this timeline. (<b>b</b>) Time series of the NS scan angles for the spacelook scenes of this timeline.</p>
Full article ">Figure 4
<p>Subset GOES-16 ABI VNIR lunar images collected on 12 April 2017. The subsets are of comparable sizes in angle, with the number of samples in the B02 image about 2 × 2 times of that for B01, B03 and B05 and 4 × 4 times of that for B04 and B06.</p>
Full article ">Figure 5
<p>Same as <a href="#remotesensing-15-01881-f004" class="html-fig">Figure 4</a>, with lunar samples (albedo larger than 0.33%) marked as −0.1%.</p>
Full article ">Figure 6
<p>Lunar irradiance ratio versus the EW scan angles for the G16 VNIR bands.</p>
Full article ">Figure 7
<p>The EW scan-angle-dependent ABI scan-mirror reflectivity for the six G16 VNIR bands.</p>
Full article ">Figure 8
<p>The normalized lunar irradiance ratio versus scan angles for G16 VNIR bands. The dashed gray lines are ±0.5% from the mean (in solid gray).</p>
Full article ">Figure 9
<p>Same as <a href="#remotesensing-15-01881-f008" class="html-fig">Figure 8</a>, but for G17.</p>
Full article ">Figure 10
<p>Time series of the mean spacelook (SPLK) values for the G16 VNIR band during the lunar chasing event on 11 February 2017. The vertical dash line is the starting time for each 15 min timeline.</p>
Full article ">Figure 11
<p>As <a href="#remotesensing-15-01881-f010" class="html-fig">Figure 10</a>, but for G17 on 30 July 2018.</p>
Full article ">Figure 12
<p>Time series of the mean SPLK values for the G17 VNIR band during the lunar chasing event on 20 October 2018. The vertical dash line is the starting time for each 15 min timeline. The gray line is refers to the second y-axis for the FPM temperature.</p>
Full article ">Figure 13
<p>Linear relationship between the SPLK count and the VNIR FPM temperature for G17 on 20 October 2018 using detector #200 as an example. The open circles detected at B01 and B02 are outliers, and the solid circles are used for the final linear regression calculation.</p>
Full article ">Figure 14
<p>The normalized lunar irradiance ratio versus the scan angle for G16 B01–B03, after the spacelook correction with Equation (8).</p>
Full article ">Figure 15
<p>The normalized lunar irradiance ratio versus the scan angle for G17 B01–B03. The spacelook count on 30 July 2018 is corrected with Equation (8). The data on 28 August 2018, 29 August 2018, and 20 October 2018 are corrected with Equation (9).</p>
Full article ">Figure 16
<p>Albedo (%) of the sample from FD Swath #11 collected on 26 July 2018 at 19:15 UTC. The image is subsampled for every four samples in the EW direction.</p>
Full article ">Figure 17
<p>Change in the delta count with the scan angle for the samples at the northern part of FD Swath #11. Delta count is the count difference between a swath sample and the corresponding spacelook calculated with Equation (2). Each dot is the mean value for the samples from the northern part of Swath #11 (lines #350–636 in <a href="#remotesensing-15-01881-f012" class="html-fig">Figure 12</a>) at the same EW scan angle (vertical line). Data were collected at 19:15:38 UTC on 26 July 2018.</p>
Full article ">Figure 18
<p>Dependence of lunar irradiance ratio on EW scan angle for G16. The B01–B03 images are calibrated with the corrected spacelook and the L1alpha image straylight correction.</p>
Full article ">Figure 19
<p>Same as <a href="#remotesensing-15-01881-f018" class="html-fig">Figure 18</a>, but for G17.</p>
Full article ">
18 pages, 2574 KiB  
Article
Estimation of Leaf Area Index in a Typical Northern Tropical Secondary Monsoon Rainforest by Different Indirect Methods
by Xiansheng Xie, Yuanzheng Yang, Wuzheng Li, Nanyan Liao, Weihu Pan and Hongxin Su
Remote Sens. 2023, 15(6), 1621; https://doi.org/10.3390/rs15061621 - 17 Mar 2023
Cited by 6 | Viewed by 2437
Abstract
The leaf area index (LAI) is a crucial indicator for quantifying forest productivity and community ecological processes. Satellite remote sensing can achieve large-scale LAI monitoring, but it needs to be calibrated and validated according to the in situ measurements on the ground. In [...] Read more.
The leaf area index (LAI) is a crucial indicator for quantifying forest productivity and community ecological processes. Satellite remote sensing can achieve large-scale LAI monitoring, but it needs to be calibrated and validated according to the in situ measurements on the ground. In this study, we attempted to use different indirect methods to measure LAI in a tropical secondary forest. These methods included the LAI-2200 plant canopy analyzer (LAI-2200), Digital Hemispherical Photography (DHP), Tracing Radiation and Architecture of Canopies (TRAC), and Terrestrial Laser Scanning (TLS) (using single-station and multi-station measurements, respectively). Additionally, we tried to correct the measured LAI by obtaining indicators of woody components and clumping effects. The results showed that the LAI of this forest was large, with estimated values of 5.27 ± 1.16, 3.69 ± 0.74, 5.86 ± 1.09, 4.93 ± 1.33, and 3.87 ± 0.89 for LAI-2200, DHP, TRAC, TLS multi-station, and TLS single-station, respectively. There was a significant correlation between the different methods. LAI-2200 was significantly correlated with all other methods (p < 0.01), with the strongest correlation with DHP (r = 0.684). TRAC was significantly correlated with TLS single-station (p < 0.01, r = 0.283). TLS multi-station was significantly correlated with TLS single-station (p < 0.05, r = 0.266). With the multi-station measurement method, TLS could maximize the compensation for measurement bias due to the shadowing effects. In general, the clumping index of this forest was 0.94 ± 0.05, the woody-to-total area ratio was 3.23 ± 2.22%, and the total correction coefficient was 1.03 ± 0.07. After correction, the LAI estimates for all methods were slightly higher than before, but there was no significant difference among them. Based on the performance assessment of existing ground-based methods, we hope to enhance the inter-calibration between methods to improve their estimation accuracy under complex forest conditions and advance the validation of remote sensing inversion of the LAI. Moreover, this study also provided a practical reference to promote the application of LiDAR technology in tropical forests. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location of the study area.</p>
Full article ">Figure 2
<p>Measurement schemes for obtaining LAI by different indirect methods. (<b>a</b>) Measurement route of LAI-2200 B-value; (<b>b</b>) Camera shooting point for DHP; (<b>c</b>) Measurement route for TRAC; (<b>d</b>) Field measurement scenes for TLS; (<b>e</b>) Instrument and target ball site locations for TLS.</p>
Full article ">Figure 3
<p>The distribution of effective LAI data estimated by different indirect methods. (<b>a</b>) LAI-2200; (<b>b</b>) DHP; (<b>c</b>) TRAC; (<b>d</b>) TLS1 (TLS multi-station); (<b>e</b>) TLS2 (TLS single-station).</p>
Full article ">Figure 4
<p>The difference of effective LAI estimated by different indirect methods. (<b>a</b>) Box plots, the central line in the box indicates the median, and the bottom and top edges indicate the 25th and 75th percentiles. Outliers are plotted using the “○” symbol; (<b>b</b>) Estimating performance (Mean ± SD). Different lowercase letters (a, b, or c) indicate significant differences between the effective LAI estimated by different methods (<span class="html-italic">p</span> &lt; 0.05). TLS1 indicates TLS multi-station, and TLS2 indicates TLS single-station.</p>
Full article ">Figure 5
<p>Correction coefficients for different sample plots. (<b>a</b>) Clumping Index; (<b>b</b>) The woody-to-total area ratio; (<b>c</b>) The total correction coefficient. Different lowercase letters (a, b) indicate significant differences between the correction coefficients of different sample plots (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 6
<p>The estimated LAI before and after correction by different indirect methods (Mean ± SD). (<b>a</b>) LAI-2200; (<b>b</b>) DHP; (<b>c</b>) TRAC; (<b>d</b>) TLS1; (<b>e</b>) TLS2; (<b>f</b>) Total. TLS1 indicates TLS multi-station, and TLS2 indicates TLS single-station.</p>
Full article ">Figure 7
<p>Treatment of woody components based on PS method. (<b>a</b>) original image; (<b>b</b>) processed image.</p>
Full article ">
16 pages, 18130 KiB  
Article
Design and Implementation of a Low-Cost Air Quality Network for the Aburra Valley Surrounding Mountains
by Andrés Yarce Botero, Santiago Lopez Restrepo, Juan Sebastian Rodriguez, Diego Valle, Julian Galvez-Serna, Elena Montilla, Francisco Botero, Bas Henzing, Arjo Segers, Arnold Heemink, Olga Lucia Quintero and Nicolás Pinel
Pollutants 2023, 3(1), 150-165; https://doi.org/10.3390/pollutants3010012 - 1 Mar 2023
Cited by 3 | Viewed by 2161
Abstract
The densest network for measuring air pollutant concentrations in Colombia is in Medellin, where most sensors are located in the heavily polluted lower parts of the valley. Measuring stations in the higher elevations on the mountains surrounding the valley are not available, which [...] Read more.
The densest network for measuring air pollutant concentrations in Colombia is in Medellin, where most sensors are located in the heavily polluted lower parts of the valley. Measuring stations in the higher elevations on the mountains surrounding the valley are not available, which limits our understanding of the valley’s pollutant dynamics and hinders the effectiveness of data assimilation studies using chemical transport models such as LOTOS-EUROS. To address this gap in measurements, we have designed a new network of low-cost sensors to be installed at altitudes above 2000 m.a.s.l. The network consists of custom-built, solar-powered, and remotely connected sensors. Locations were strategically selected using the LOTOS-EUROS model driven by diverse meteorology-simulated fields to explore the effects of the valley wind representation on the transport of pollutants. The sensors transmit collected data to internet gateways for posterior analysis. Various tests to verify the critical characteristics of the equipment, such as long-range transmission modeling and experiments with an R score of 0.96 for the best propagation model, energy power system autonomy, and sensor calibration procedures, besides case exposure to dust and water experiments, to ensure IP certifications. An inter-calibration procedure was performed to characterize the sensors against reference sensors and describe the observation error to provide acceptable ranges for the data assimilation algorithm (<10% nominal). The design, installation, testing, and implementation of this air quality network, oriented towards data assimilation over the Aburrá Valley, constitute an initial experience for the simulation capabilities toward the system’s operative capabilities. Our solution approach adds value by removing the disadvantages of low-cost devices and offers a viable solution from a developing country’s perspective, employing hardware explicitly designed for the situation. Full article
(This article belongs to the Section Air Pollution)
Show Figures

Figure 1

Figure 1
<p>The figure shows two patterns for the hypothetical simulation for regional atmospheric pollutant transport. The panels in the image represent simulation results using the LOTOS-EUROS CTM for Medellín (<bold>a</bold>) and Rionegro (<bold>b</bold>) with a point source of nitrogen oxides equivalent to the emitted by the city in one day (1000 kg of nitrogen oxides/hour). These releases were over the background concentration (The background simulation was removed, which makes that the NO<sub>2</sub> was liberated from a specific point) during the hours of 06:00–20:00 on the third day of the simulation. The deposition was monitored from day 3 through day 9 (days 1 and 2 were used as model ramp-up time).</p>
Full article ">Figure 2
<p>Transversal cut comparisons, three-dimensional snapshot of the model output over the valley and assimilated outputs of a low-cost sensor network inside the valley.</p>
Full article ">Figure 3
<p>(<bold>a</bold>) Locations selected for the air quality network disposition around the surrounding mountains. The red cross is located at the point of view in the right direction from the arrow of the (<bold>c</bold>) photo. (<bold>b</bold>) table with the code and station coordinates. (<bold>c</bold>) The landscape of Medellín with the image of the east mountains of the valley and the disposition of two stations over the urban areas of this deep-seated valley.</p>
Full article ">Figure 4
<p>Data assimilation system with new measurements in the high part of the valley complementing the SIATA network.</p>
Full article ">Figure 5
<p>Flowchart with the various subsystems that comprise the Simple-4 partial centralized architecture. The Energy Power Supply (EPS) subsystem with the MPPT module controls the charge and discharge cycles of the battery and the load consumption. A dedicated microcontroller collects the payload information and, once the data is preprocessed, is delivered to the OB&amp;DH to command the communication subsystem. The other secondary system is the thermal monitoring system, a transversal support system conformed by the different thermal sensors in each PCB layer.</p>
Full article ">Figure 6
<p>(<bold>a</bold>) SimpleVital (Command On Board &amp; Data Handling (CB&amp;DH)), (<bold>b</bold>) SimplePayload (Gases/Pollution), (<bold>c</bold>) SimplePower (Energy Power Supply (EPS)), (<bold>d</bold>) SimpleCOMM (Communications Subsystems (COMM).</p>
Full article ">Figure 7
<p>(<bold>a</bold>) Vertical cut and Fresnel distance. (<bold>b</bold>) Spatial intensity over the Aburra Valley and (<bold>c</bold>) model generated from the different signal intensities and experimental locations.</p>
Full article ">Figure 8
<p>(<bold>a</bold>) Current, Voltage, and Power (Watts) before the battery and before the load measure for four days, (<bold>b</bold>) Energy supply simulation for the percentage of full, empty, not captured, and missing energy at the different point locations in the mountains around Medellín.</p>
Full article ">Figure 9
<p>(<bold>a</bold>) The case developed for the measurement device, (<bold>b</bold>) Transversal cut to see the inside disposition of the elements and the different chambers specially designed to orientate the air mass flux through the payload bay where the gas sensors are located, and (<bold>c</bold>) Render visualization of the case to appreciate the translucid top dome and the lateral airflow channels.</p>
Full article ">Figure 10
<p>(<bold>a</bold>) Infrared thermal images under operations for different operation routines of the module. IP5X IP6X test for the device. (<bold>b</bold>) Temperature and Humidity stationary chamber test. (<bold>c</bold>) dust and wind tests.</p>
Full article ">Figure 11
<p>Qualitative comparison of the device with the Vaisala AQT400 ®calibrated equipment.</p>
Full article ">Figure 12
<p>(<bold>a</bold>) Locations of a Simple 4 compared with two sensors of the Ciudadanos network at two altitudes. (<bold>b</bold>) PM10 concentrations. Interesting to see how the mountains can “contain” the contamination; it is clear these days that the concentration is greater within the valley and, on the other side, also greater than that recorded by the device. At the altitude of the module in the Columbus, that more wooded corridor of the slopes of the eastern mountains of the Aburrá Valley can also be seen.</p>
Full article ">Figure A1
<p>Different meteorology field input (ECMWF and WRF) the LOTOS-EUROS CTM simulating particulate matter 2.5 for the contingency period of the 3, 4, and 5 of March 2019. The red dots are the station network’s location. There is a noticeable difference in the particulate matter concentration fields that, in one situation, respect the other and is trapped by the valley structure.</p>
Full article ">
16 pages, 15503 KiB  
Technical Note
Soil Moisture Retrieval from Multi-GNSS Reflectometry on FY-3E GNOS-II by Land Cover Classification
by Cong Yin, Feixiong Huang, Junming Xia, Weihua Bai, Yueqiang Sun, Guanglin Yang, Xiaochun Zhai, Na Xu, Xiuqing Hu, Peng Zhang, Jinsong Wang, Qifei Du, Xianyi Wang and Yuerong Cai
Remote Sens. 2023, 15(4), 1097; https://doi.org/10.3390/rs15041097 - 17 Feb 2023
Cited by 8 | Viewed by 2460
Abstract
The reflected GNSS signals at the L-band is significantly advantageous in soil moisture monitoring as they are sensitive to the dielectric properties determined by the volumetric water content of topsoil, and they can penetrate vegetation, except in very dense forests. The Global Navigation [...] Read more.
The reflected GNSS signals at the L-band is significantly advantageous in soil moisture monitoring as they are sensitive to the dielectric properties determined by the volumetric water content of topsoil, and they can penetrate vegetation, except in very dense forests. The Global Navigation satellite system Occultation Sounder (GNOS-II) with a reflectometry technique onboard the Fengyun-3E (FY-3E) satellite, launched on 5 July 2021, is the first mission that can receive reflected Global Navigation Satellite System (GNSS) signals from GPS, BeiDou and Galileo systems. This paper presents the soil moisture retrieval results from the FY-3E GNOS-II mission using 16 months of data. In this study, the reflectivity observations from different GNSS systems were firstly intercalibrated with some differences analyzed. Observations were also corrected by considering vegetation attenuation for 16 different land cover classifications. Finally, an empirical model was constructed for volumetric soil moisture (VSM) estimation, where the reflectivity of GNOS-II was linearly related to the SMAP reference soil moisture for each 36 km ease grid. The overall root-mean-square error of the retrieved soil moisture is 0.049 compared with the SMAP product, and 0.054 compared with the in situ data. The results of the three GNSS systems show similar levels of accuracy. This paper, for the first time, demonstrates the feasibility of global soil moisture retrieval using multiple GNSS signals. Full article
Show Figures

Figure 1

Figure 1
<p>DDMs over the land surface from (<bold>a</bold>) BDS C29, (<bold>b</bold>) GPS PRN20 and (<bold>c</bold>) GAL E11.</p>
Full article ">Figure 2
<p>The overlap between GNOS-II tracks and SMAP tracks (<bold>a</bold>) The overlap between GNOS-II tracks and SMAP observations in one day (<bold>b</bold>) The overlap between GNOS-II tracks and SMAP observations in three days (the green tracks indicate the overlap of two measurements, whereas the red tracks represent the GNOS-R observations that are not coincident with the SMAP measurements).</p>
Full article ">Figure 2 Cont.
<p>The overlap between GNOS-II tracks and SMAP tracks (<bold>a</bold>) The overlap between GNOS-II tracks and SMAP observations in one day (<bold>b</bold>) The overlap between GNOS-II tracks and SMAP observations in three days (the green tracks indicate the overlap of two measurements, whereas the red tracks represent the GNOS-R observations that are not coincident with the SMAP measurements).</p>
Full article ">Figure 3
<p>Distribution of in situ soil moisture stations.</p>
Full article ">Figure 4
<p>Density scatter plots for the comparison between reflectivity over sea ice from different GNSS systems (7 July–31 October 2022). Bias is also displayed on each sub-figure.</p>
Full article ">Figure 5
<p>Density scatter plots for the comparison between reflectivity from different GNSS systems (7 July–31 October 2022): (<bold>a</bold>) before calibration; (<bold>b</bold>) after calibration.</p>
Full article ">Figure 6
<p>The GNOS-II reflectivity of land surface (10 July–9 September 2021).</p>
Full article ">Figure 7
<p>The estimated topography attenuation at an ease grid of 36 km.</p>
Full article ">Figure 8
<p>The histogram of the incidence angles for the observations of FY-3E GNOS-II.</p>
Full article ">Figure 9
<p>GNOS-II global soil moisture estimation (7 July–6 August 2022).</p>
Full article ">Figure 10
<p>GNOS-II soil moisture vs. SMAP soil moisture. (<bold>a</bold>) GPS; (<bold>b</bold>) BDS; (<bold>c</bold>) GAL.</p>
Full article ">Figure 11
<p>GNOS/SMAP soil moisture vs. ISMN in situ observations. (<bold>a</bold>) ARM Medford; (<bold>b</bold>) ARM Ashton; (<bold>c</bold>) ARM Waukomis; (<bold>d</bold>) SCAN Phillipsburg; (<bold>e</bold>) TxSON CR 200-18.</p>
Full article ">Figure 12
<p>The histogram of ubRMSE (<bold>a</bold>) GNOS vs. in situ (<bold>b</bold>) SMAP vs. in situ.</p>
Full article ">
15 pages, 2804 KiB  
Article
Imputation of Missing Parts in UAV Orthomosaics Using PlanetScope and Sentinel-2 Data: A Case Study in a Grass-Dominated Area
by Francisco R. da S. Pereira, Aliny A. Dos Reis, Rodrigo G. Freitas, Stanley R. de M. Oliveira, Lucas R. do Amaral, Gleyce K. D. A. Figueiredo, João F. G. Antunes, Rubens A. C. Lamparelli, Edemar Moro and Paulo S. G. Magalhães
ISPRS Int. J. Geo-Inf. 2023, 12(2), 41; https://doi.org/10.3390/ijgi12020041 - 28 Jan 2023
Cited by 1 | Viewed by 2347
Abstract
The recent advances in unmanned aerial vehicle (UAV)-based remote sensing systems have broadened the remote sensing applications for agriculture. Despite the great possibilities of using UAVs to monitor agricultural fields, specific problems related to missing parts in UAV orthomosaics due to drone flight [...] Read more.
The recent advances in unmanned aerial vehicle (UAV)-based remote sensing systems have broadened the remote sensing applications for agriculture. Despite the great possibilities of using UAVs to monitor agricultural fields, specific problems related to missing parts in UAV orthomosaics due to drone flight restrictions are common in agricultural monitoring, especially in large areas. In this study, we propose a methodological framework to impute missing parts of UAV orthomosaics using PlanetScope (PS) and Sentinel-2 (S2) data and the random forest (RF) algorithm of an integrated crop–livestock system (ICLS) covered by grass at the time. We validated the proposed framework by simulating and imputing artificial missing parts in a UAV orthomosaic and then comparing the original data with the model predictions. Spectral bands and the normalized difference vegetation index (NDVI) derived from PS, as well as S2 images (separately and combined), were used as predictor variables of the UAV spectral bands and NDVI in developing the RF-based imputation models. The proposed framework produces highly accurate results (RMSE = 6.77–17.33%) with a computationally efficient and robust machine-learning algorithm that leverages the wealth of empirical information present in optical satellite imagery (PS and S2) to impute up to 50% of missing parts in a UAV orthomosaic. Full article
(This article belongs to the Special Issue Geomatics in Forestry and Agriculture: New Advances and Perspectives)
Show Figures

Figure 1

Figure 1
<p>Location of the study area in the western region of the state of São Paulo, Brazil, and the UAV orthomosaic (true colour composite red-green-blue (RGB):321) of the study area in August 2019.</p>
Full article ">Figure 2
<p>Spectral (bandwidth and central wavelength) characteristics of UAV orthomosaic, PlanetScope, and Sentinel-2A imagery.</p>
Full article ">Figure 3
<p>Flowchart for UAV and satellite imagery processing and filling in the missing parts in UAV orthomosaics based on the RF algorithm and ancillary data derived from satellite image.</p>
Full article ">Figure 4
<p>Simulation of the missing parts (white strips) in the original UAV orthomosaic: (<b>a</b>) missing part of 10%—MP-10%, (<b>b</b>) missing part of 30%—MP-30%, (<b>c</b>) missing part of 50%—MP-50%.</p>
Full article ">Figure 5
<p>The relative importance of the predictor variables was measured by the variable importance metric in the eighteen RF-based imputation models using both PlanetScope and Sentinel-2 data.</p>
Full article ">Figure 6
<p>Comparison between the true colour composite (RGB:321) of the original UAV orthomosaic (<b>a</b>) and the imputed UAV orthomosaic—10% (<b>b</b>), 30% (<b>c</b>), and 50% (<b>d</b>) of image prediction. Squares #1, #2, and #3 show, in more detail, the original UAV orthomosaic, the introduced missing parts, and the imputed UAV orthomosaic in the highest level of missing parts assessed in this study (MP–50%).</p>
Full article ">Figure 7
<p>Original UAV NDVI image (<b>a</b>); predicted NDVI images with missing parts of 10% (<b>b</b>), 30% (<b>d</b>), and 50% (<b>f</b>); and the calculated NDVI images (based on the imputed R and NIR spectral bands) with missing parts of 10% (<b>c</b>), 30% (<b>e</b>), and 50% (<b>g</b>).</p>
Full article ">
17 pages, 4476 KiB  
Article
Using a Tandem Flight Configuration between Sentinel-6 and Jason-3 to Compare SAR and Conventional Altimeters in Sea Surface Signatures of Internal Solitary Waves
by Jorge M. Magalhaes, Ian G. Lapa, Adriana M. Santos-Ferreira, José C. B. da Silva, Fanny Piras, Thomas Moreau, Samira Amraoui, Marcello Passaro, Christian Schwatke, Michael Hart-Davis, Claire Maraldi and Craig Donlon
Remote Sens. 2023, 15(2), 392; https://doi.org/10.3390/rs15020392 - 8 Jan 2023
Cited by 8 | Viewed by 3771
Abstract
Satellite altimetry has been providing a continuous record of ocean measurements with numerous applications across the entire range of ocean sciences. A reference orbit has been used since 1992 with TOPEX/Poseidon, which was repeated in the Jason missions, and in the newly launched [...] Read more.
Satellite altimetry has been providing a continuous record of ocean measurements with numerous applications across the entire range of ocean sciences. A reference orbit has been used since 1992 with TOPEX/Poseidon, which was repeated in the Jason missions, and in the newly launched Sentinel-6 Michael Freilich (in November 2020) to continually monitor the trends of sea level rise and other properties of the sea surface. These multidecadal missions have evolved alongside major technological advances, whose measurements are unified into a single data record owing to continuous intercalibration and validation efforts. However, the new Sentinel-6 provides synthetic aperture radar (SAR) processing, which improves the along-track resolution of conventional altimeters from a few kilometres (e.g., for Jason-3) to about 300 m. This means a major leap in sampling towards higher frequencies of the ocean spectrum, which inevitably means reconciling the assumption of a uniform Brown surface between the footprints of the larger kilometre-scale conventional altimetry and those of the finer-scale SAR altimetry. To explore this issue, this study uses the vantage point of the Sentinel-6/Jason-3 tandem phase to compare simultaneous sea surface signatures of large-scale Internal Solitary Waves (ISWs) between SAR and conventional altimetry. These waves can modulate the sea surface into arrayed sections of increased and decreased roughness with horizontal scales up to 10 km, which inflict sharp transitions between increased and decreased backscatter in the radar altimeters. It is found that Sentinel-6 can provide more detailed structures of ISWs in standard level-2 products, when compared with those from the conventional Jason-3 (similarly to previous results reported from the SAR altimeter from Sentinel-3). However, a new and striking feature is found when comparing the radar backscatter between Sentinel-6 and Jason-3, which are in opposite phases in the ISWs. These intriguing results are discussed in light of the intrinsically different acquisition geometries of SAR and conventional altimeters as well as possible implications thereof. Full article
Show Figures

Figure 1

Figure 1
<p>Schematic representation of an ISW propagating with phase velocity <math display="inline"><semantics> <mover accent="true"> <mi>c</mi> <mo>→</mo> </mover> </semantics></math> and its sea surface roughness patterns with leading rough and trailing slick-like sections. Footprints are shown over the slick-like section for Jason-3 (round-shaped) and Sentinel-6 (rectangular)—assumed at nadir but shown at an angle for simplicity. Note that spatial scales are chosen to highlight that, when properly aligned, the smaller footprint in a SAR altimeter can sample consecutive echoes in each section of the ISW, while the larger footprint in a conventional altimeter can obtain mixed contributions from both sections simultaneously. However, both satellites can obtain mixed contributions as illustrated in the alternative view on the back left-hand side.</p>
Full article ">Figure 2
<p>Sentinel-1 image over the Banda Sea (see inset in top-left corner) acquired 7 March (2022) at 10h01m UTC. Three ISW packets are seen propagating to the north/northwest separated by typical semidiurnal periods and wavelengths (i.e., from −1M2 to +1M2, with packets separated approximately by λ<sub>packet</sub> ≈ 140 km). Note that, the width of the leading ISWs in each packet (λ<sub>ISW</sub>) is around 2 to 3 km. For reference, the ground-track of Sentinel-6/Jason-3 is shown in a blue line for pass 253 and the blue circle marks the location of the leading ISW-like signals highlighted in <a href="#remotesensing-15-00392-f003" class="html-fig">Figure 3</a>.</p>
Full article ">Figure 3
<p>Panels (<b>a</b>–<b>c</b>) show <span class="html-italic">σ</span><sub>0</sub>, SSHAs and SWHs obtained from Sentinel-6 (S6) and Jason-3 (J3) MLE4 level-2 products (at 20 Hz, except SSHAs for MLE4) in the Banda Sea (8 March 2022, see <a href="#app1-remotesensing-15-00392" class="html-app">Supplemental Materials S1</a>). The radar backscatter from Jason-3 (J3) is shown with an offset to highlight its opposite modulations in the vicinities of the ISWs (leading ISW marked with a blue rectangle). Panels (<b>d</b>–<b>f</b>), same as previous panels for Jason-3 processed with ALES and the Adaptive retrackers (at 20 Hz). Note that panel (<b>d</b>) also shows the radar backscatter from Sentinel-6 with an offset to highlight its correlation with the ALES and the Adaptive retrackers.</p>
Full article ">Figure 4
<p>Similar to <a href="#remotesensing-15-00392-f003" class="html-fig">Figure 3</a> in the Celebes Sea (5 January 2022, see <a href="#app1-remotesensing-15-00392" class="html-app">Supplemental Materials S1</a>). Note that in this case, the broken lines in panel (<b>a</b>,<b>b</b>) indicates missing data in MLE4 level-2 products.</p>
Full article ">Figure 5
<p>Similar to <a href="#remotesensing-15-00392-f003" class="html-fig">Figure 3</a> in the South China Sea (15 May 2021, see <a href="#app1-remotesensing-15-00392" class="html-app">Supplemental Materials S1</a>). Note that in this case, the broken lines indicate missing data in MLE4 level-2 products in panels (<b>a</b>–<b>c</b>) and in the ALES retracker in panels (<b>e</b>,<b>f</b>).</p>
Full article ">Figure 6
<p>Correlations coefficients (<span class="html-italic">R</span>) for all cases listed in <a href="#app1-remotesensing-15-00392" class="html-app">Supplementary Materials S1</a> and computed as described in <a href="#sec3-remotesensing-15-00392" class="html-sec">Section 3</a>. Green and blue circles represent correlations between Sentinel-6 and Jason-3 (ALES/Adaptive, respectively), and red circle represent correlations between Sentinel-6 and Jason-3 (MLE4). Two sets of correlation coefficients are shown, listed as follows: on the left for the leading ISWs (<span class="html-italic">R<sub>ISWs</sub></span>), and on the right for the waves’ background conditions (<span class="html-italic">R<sub>Background</sub></span>). A representative case is show for the Celebes Sea using data smoothed with a running mean of about 10 km (see also <a href="#remotesensing-15-00392-f004" class="html-fig">Figure 4</a>). Note that correlation coefficients for the South China Sea are only shown for the waves’ background (owing to missing data in the ISWs, see <a href="#remotesensing-15-00392-f005" class="html-fig">Figure 5</a>).</p>
Full article ">Figure 7
<p>Overview of selected cases listed in <a href="#app1-remotesensing-15-00392" class="html-app">Table S1 (in Supplementary Materials S1)</a> showing modulations in radar backscatter (<math display="inline"><semantics> <mrow> <mo>Δ</mo> <msub> <mi>σ</mi> <mn>0</mn> </msub> </mrow> </semantics></math>), SSHAs (<math display="inline"><semantics> <mo>Δ</mo> </semantics></math>SSHA) and SWHs (<math display="inline"><semantics> <mo>Δ</mo> </semantics></math>SWH) in the leading ISWs. Note that for consistency, all values assume the waves are travelling rightwards (with phase velocity <span class="html-italic">C</span>). In the case of <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msub> <mi>σ</mi> <mn>0</mn> </msub> </mrow> </semantics></math>, each pair represents the backscatter modulation in the rear/leading sections in comparison with an unperturbed background (taken ahead of the ISW). Yellow circles mark either missing values in level-2 products (NA) or estimates that are not in agreement with two-layer solitary wave theory.</p>
Full article ">Figure 8
<p>Similar to <a href="#remotesensing-15-00392-f003" class="html-fig">Figure 3</a>a, but comparing <math display="inline"><semantics> <mrow> <msub> <mi>σ</mi> <mn>0</mn> </msub> </mrow> </semantics></math> for Sentinel-6 (S6) level-2 products (at 20 Hz) in LRM (i.e., conventional altimeter) modes, and data from MLE4 in Jason-3 (J3).</p>
Full article ">
14 pages, 20752 KiB  
Article
Combining Artificial Neural Network and Driver–Pressure–State–Impact–Response Approach for Evaluating a Mediterranean Lake
by Christos Tsitsis, Dimitrios E. Alexakis, Konstantinos Moustris and Dimitra E. Gamvroula
Water 2023, 15(2), 266; https://doi.org/10.3390/w15020266 - 8 Jan 2023
Cited by 6 | Viewed by 2125
Abstract
The main objective of this research was to evaluate the surface water system of Lake Vegoritida (Region of Central Macedonia, Greece). The Driver–Pressure–State–Impact–Response (DPSIR) methodological approach was used. The analysis includes data from three (3) stations monitoring point source pollution and recording the [...] Read more.
The main objective of this research was to evaluate the surface water system of Lake Vegoritida (Region of Central Macedonia, Greece). The Driver–Pressure–State–Impact–Response (DPSIR) methodological approach was used. The analysis includes data from three (3) stations monitoring point source pollution and recording the most critical water quality measurement parameters in a time series data analysis from 1983 to 1997. The data will contribute to the analysis and was used to investigate, identify, and evaluate possible sources of chemical and ecological changes recorded in the lake. The artificial neural network (ANN) is a valuable tool for making predictions based on the water quality data set. The findings highlighted the increased concentration of nutrients that contribute to the presence of eutrophic conditions, while their seasonal variability is mainly due to factors, such as water level fluctuations and biological processes in the lake. The above, combined with the critical biotic indicators and factors alongside the reduction in biodiversity, indicated that only the most resistant species survive, confirming the previous finding. In Greece, systematic monitoring and reporting programs have recently been implemented, such as the ECOFRAME scheme and the guidelines proposed by the “Intercalibration Group for Mediterranean Lakes”. The water quality status could be classified as “High”, “High to Good”, and “High to Poor”, respectively, while the overall ecological assessment tends to change to poor conditions. The actions required at an early stage concern the planning of programs and actions that contribute to the sustainable management of land uses and the reduction in point sources of pollution, as well as the reduction of the applied quantities of agrochemicals on the cultivated land in the study area. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Map of Greece presenting the location of Lake Vegoritida; (<b>b</b>) Map of study area showing the location of monitoring stations (modified from Google Earth [<a href="#B32-water-15-00266" class="html-bibr">32</a>]).</p>
Full article ">Figure 2
<p>Typical architecture of a feed-forward multi-layer perceptron artificial neural (MLP-ANN) (modified from [<a href="#B47-water-15-00266" class="html-bibr">47</a>]).</p>
Full article ">Figure 3
<p>Determination coefficient between observed and predicted TP values by ANN.</p>
Full article ">Figure 4
<p>Time series of selected water quality parameters of Lake Vegoritida during the monitoring period 1983–1997: (<b>a</b>) water temperature; (<b>b</b>) pH; (<b>c</b>) CND; (<b>d</b>) DO; (<b>e</b>) TN; and (<b>f</b>) TP (data were provided by HMRDF 2019).</p>
Full article ">Figure 5
<p>Average values of water quality parameters of the studied surface water body and chemical conditions typology as stated by ECOFRAME and MED-GIG system (data were provided by HMRDF 2019, <sup>a</sup> Criteria given by ECOFRAME [<a href="#B15-water-15-00266" class="html-bibr">15</a>]; <sup>b</sup> Criteria given by Poikane [<a href="#B16-water-15-00266" class="html-bibr">16</a>]).</p>
Full article ">
Back to TopTop