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

Search Parameters:
Keywords = intercalibration

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
29 pages, 12604 KiB  
Article
The Characterization of the Railroad Valley Playa Test Site Using the DESIS Imaging Spectrometer from the Space Station Orbit
by Mohammad H. Tahersima, Kurtis Thome, Brian N. Wenny, Derrick Lampkin, Norvik Voskanian, Sarah Eftekharzadeh Kay and Mehran Yarahmadi
Remote Sens. 2025, 17(3), 396; https://doi.org/10.3390/rs17030396 - 24 Jan 2025
Viewed by 481
Abstract
The reflectance-based vicarious calibration approach uses measurements at well-understood test sites to provide top-of-atmosphere reference reflectance values suitable for inter-calibration approaches and does not require coincident views. The challenge is that results from such data may suffer from high variability from day to [...] Read more.
The reflectance-based vicarious calibration approach uses measurements at well-understood test sites to provide top-of-atmosphere reference reflectance values suitable for inter-calibration approaches and does not require coincident views. The challenge is that results from such data may suffer from high variability from day to day. Data from high-quality sensors, such as the imaging spectrometers on the International Space Station (ISS) platform, provide an opportunity to use improved fine spectral information about the test sites with various sun/sensor geometries and site surface and atmospheric conditions to improve the test sites’ characterization. The results here are based on data from the DLR Earth Sensing Imaging Spectrometer (DESIS) instrument installed on the ISS since 2018 combined with output from the Radiometric Calibration Network (RadCalNet) site at Railroad Valley Playa (RRV) to decouple the effects of sun/sensor geometry from the RadCalNet predictions. The approach here uses the precessing orbit of the ISS to allow similar sensor view zenith angles at varying sun angles over short periods that limit the impact of any sensor changes and highlight the bi-directional effects of the surface reflectance and atmospheric conditions. DESIS data collected at (i) similar solar angles but varying view angles, (ii) similar sensor angles and varying solar angles, and (iii) similar scatter angles are compared. The DESIS results indicate that the top-of-atmosphere reflectance spectra for RRV at similar solar zenith angles but with varying sensor viewing angles provide more consistent data than those with varying solar zenith but with similar sensor viewing angles. In addition, comparisons of reflectance spectra of the site performed in terms of the sensor view scatter angle show good agreement, indicating that a directional reflectance correction would be straightforward and could offer a significant improvement in the use of RadCalNet data. The work shows that observations from imaging spectroscopy data from DESIS, and eventually Earth Surface Mineral Dust Source Investigation (EMIT), Surface Biology and Geology (SBG), and the climate-quality sensor CLARREO Pathfinder (CPF), provide the opportunity for the development of a model-based, SI-traceable prediction of at-sensor radiance over the RRV site that would serve as the basis for similar site characterizations with error budgets valid for arbitrary view and illumination angles. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>A</b>) A sample Band 60 (552.91 center wavelength) DESIS image of the RRV captured on 22 August 2018 and (<b>B</b>) a zoomed-in version of the DESIS image showing the 34-pixel-by-34-pixel area of the test site. The white box in (<b>A</b>) has the same footprint as the image in (<b>B</b>). The white box of (<b>B</b>) shows the 34-pixel-by-34-pixel area of the test site.</p>
Full article ">Figure 2
<p>(<b>A</b>) Pixel-wise TOA spectral radiance (1156 blue spectra shown) for the DESIS image captured on 4 February 2022. (<b>B</b>) Site-averaged TOA spectral radiance (blue solid line) and site-averaged reflectance (red solid line) from the same DESIS image, and corresponding TOA reflectance reported by in situ-based RadCalNet data from the same date and within a few minutes of the DESIS image capture time.</p>
Full article ">Figure 3
<p>This figure shows the intercomparison of TOA reflectance spectra (the ratio of the red curves shown in <a href="#remotesensing-17-00396-f002" class="html-fig">Figure 2</a>) of the RRV site on 4 February 2022. Dashed lines represent the absolute uncertainty (k = 2) of the intercomparison.</p>
Full article ">Figure 4
<p>ISS is on non-SSO. Thus, the local overpass time of ISS-mounted instruments changes from scene to scene. A three-day orbit cycle of the ISS and swath for DESIS (<b>A</b>) is plotted using the Committee on Earth Observation Satellites (CEOS) Visualization Environment (COVE) [<a href="#B37-remotesensing-17-00396" class="html-bibr">37</a>]. DESIS has 30 m spatial resolutions with a 30 km swath width at a 400 km altitude. The data acquisition geometry of the data is shown in (<b>B</b>). Blue symbols represent the sensor geometry of DESIS, and red symbols represent sun geometry relative to the RRV site at the capturing times. The azimuth angles are given as 0°N, 90°E, 180°S, and 270°W. The zenith angles are 0° at the center and increase outwards.</p>
Full article ">Figure 5
<p>Geometry of incident and reflected beams to/from the surface of the radiometric calibration site. The in situ radiometers installed at the site always have a nadir view.</p>
Full article ">Figure 6
<p>(<b>A</b>) Seasonal and (<b>B</b>) temporal sun/sensor geometries for 41 DESIS images capture the RRV’s full view from December 2018 to December 2023. Sun zenith angles are marked in red, and sensor zenith angles are marked in blue. Clusters of scenes with similar sun or sensor angles captured close in time are identified from (<b>B</b>). The symbols facing up/down represent scenes captured when the ISS moved in ascending/descending orbits. The darker symbols are for DESIS scenes with associated RadCalNet data. Conversely, the lighter symbols represent DESIS scenes without reported RadCalNet data.</p>
Full article ">Figure 7
<p>(<b>A</b>) Time-agnostic sun/sensor geometries and (<b>B</b>) scatter angles of the 41 DESIS images with a clear and complete view of the RRV from December 2018 to December 2023. Time-agnostic clusters of scenes with similar sun or sensor angles are identified from (<b>A</b>). The symbols facing up/down represent scenes captured when the ISS moved in ascending/descending orbits. The darker symbols are for DESIS scenes that have associated RadCalNet data. Conversely, the lighter symbols represent DESIS scenes without reported RadCalNet data and that cannot be used for absolute or relative double-ratio comparisons.</p>
Full article ">Figure 8
<p>This figure shows (<b>A</b>) absolute ratios of DESIS’s TOA reflectance spectra to those of RadCalNet, (<b>B</b>) relative ratios of DESIS’s TOA reflectance spectra to a reference TOA reflectance spectrum of DESIS, and (<b>C</b>) relative double ratios of the absolute ratio spectra to a reference absolute ratio spectrum.</p>
Full article ">Figure 9
<p>Percent standard deviation of mean disagreement between at-sensor reflectance prediction and DESIS-derived reflectance.</p>
Full article ">Figure 10
<p>Ratios of TOA reflectance spectra of DESIS observations of the RRV site captured during June 2019 relative to a reference TOA reflectance spectrum of a DESIS image captured on 28 June 2019.</p>
Full article ">Figure 11
<p>The ratio of TOA reflectance spectra of DESIS observations of the RRV site captured during August 2019 relative to a reference TOA reflectance spectrum of a DESIS image captured on 19 August 2019.</p>
Full article ">Figure 12
<p>The ratio of TOA reflectance spectra of DESIS observations of the RRV site captured during October 2019 relative to a reference TOA reflectance spectrum of a DESIS image captured on 29 October 2019.</p>
Full article ">Figure 13
<p>The ratio of TOA reflectance spectra of DESIS observations of the RRV site captured in April–May 2022 relative to a reference TOA reflectance spectrum of a DESIS image captured on 11 April 2022.</p>
Full article ">Figure 14
<p>The ratio of site-averaged TOA reflectance spectra of DESIS observations of the RRV site captured with similar sun/sensor geometries relative to a reference site-averaged TOA reflectance spectrum captured on 28 June 2019.</p>
Full article ">Figure 15
<p>Absolute ratio (DESIS<sub>i</sub>/RCN<sub>i</sub>) spectra clustered based on view zenith angles in five degree increments from 0 to 5 degrees (<b>A</b>), 10 to 15 degrees (<b>B</b>), 5 to 10 degrees (<b>C</b>), and 15 to 20 degrees (<b>D</b>).</p>
Full article ">Figure 16
<p>Spectral solar irradiance terms for solar zenith angle sweep from 0 degrees (blue) to 80 degrees (red). Calculated for 30% surface reflectance using MODTRAN mid-latitude summer model. Direct (<b>A</b>) and Diffuse (<b>B</b>) components of the spectral solar irradiance are shown for a sweep of solar zenith angles from 0 (shown in blue) to 80 degrees (shown in red) in 10 degrees increments.</p>
Full article ">Figure 17
<p>Magnitudes of correction factors to account for MODTRAN view angle adjustments to match the exact sun/sensor geometries of each DESIS scene. There are 15 spectral lines corresponding to the 15 matchups (explained in <a href="#sec4-remotesensing-17-00396" class="html-sec">Section 4</a> and listed in <a href="#remotesensing-17-00396-t001" class="html-table">Table 1</a>). Results show a small effect of less than 0.5% between 600 nm and 1000 nm. The MODTRAN angle-adjustment correction factor can gradually increase up to ~2% toward the shorter wavelengths.</p>
Full article ">Figure 18
<p>Sun angle (<b>A</b>) and view angle (<b>B</b>) polar plots of the directional correction factors, which are calculated by taking a mean of the ratio of DESIS-observed TOA reflectance spectra to the predicted TOA reflectance spectra calculated by the radiative transfer codes that consider sun/sensor geometries of DESIS observations. The correction factor for each scene is marked with a distinctive marker on both the sun angle and view angle polar plots. The markers are color-coded with magnitudes of the correction factor according to the color bar shown.</p>
Full article ">Figure 19
<p>Bi-directional reflectance factors (BRFs), at 440 nm and 550 nm, 660 nm, and 860 nm, for sun angle <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>θ</mi> </mrow> <mrow> <mi>s</mi> </mrow> </msub> <mo>=</mo> <msup> <mrow> <mn>20</mn> </mrow> <mrow> <mo>∘</mo> </mrow> </msup> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>ϕ</mi> </mrow> <mrow> <mi>s</mi> </mrow> </msub> <mo>=</mo> <msup> <mrow> <mn>180</mn> </mrow> <mrow> <mo>∘</mo> </mrow> </msup> </mrow> </semantics></math>, computed based on the correction factors of the 15 DESIS scenes shown in <a href="#remotesensing-17-00396-f018" class="html-fig">Figure 18</a>. The principal planes are the <math display="inline"><semantics> <mrow> <msup> <mrow> <mn>0</mn> </mrow> <mrow> <mo>∘</mo> </mrow> </msup> <mo>−</mo> <msup> <mrow> <mn>180</mn> </mrow> <mrow> <mo>∘</mo> </mrow> </msup> </mrow> </semantics></math> directions with forward reflection at the top (less reflective compared to the nadir view) and the back reflection at the bottom (more reflective compared to the nadir view).</p>
Full article ">Figure 20
<p>BRFs of <a href="#remotesensing-17-00396-f019" class="html-fig">Figure 19</a> relative to respective BRFs at the nadir view. These relative BRF polar plots map view angle corrections needed for sun angles at <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>θ</mi> </mrow> <mrow> <mi>s</mi> </mrow> </msub> <mo>=</mo> <msup> <mrow> <mn>20</mn> </mrow> <mrow> <mo>∘</mo> </mrow> </msup> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>ϕ</mi> </mrow> <mrow> <mi>s</mi> </mrow> </msub> <mo>=</mo> <msup> <mrow> <mn>180</mn> </mrow> <mrow> <mo>∘</mo> </mrow> </msup> </mrow> </semantics></math>.</p>
Full article ">
17 pages, 7398 KiB  
Article
Inter-Calibration and Limb Correction of FY-3D/E MWTS for Long-Term Gridded Microwave Brightness Temperature Dataset
by Xinlu Xia, Mingjian Zeng, Xiaochun Luo, Xiao Shi, Rongsheng Jiang, Xinyi Yuan, Xiaozhuo Sang, Fei Tang and Xu Xu
Remote Sens. 2025, 17(1), 158; https://doi.org/10.3390/rs17010158 - 5 Jan 2025
Viewed by 488
Abstract
The Microwave Temperature Sounder-3 (MWTS-3) onboard the Chinese FengYun-3E (FY-3E) satellite is the third generation of Chinese microwave temperature sounder. Based on MWTS-2, the number of MWTS-3 channels has been increased from 13 to 17, which can observe the atmospheric temperature and water [...] Read more.
The Microwave Temperature Sounder-3 (MWTS-3) onboard the Chinese FengYun-3E (FY-3E) satellite is the third generation of Chinese microwave temperature sounder. Based on MWTS-2, the number of MWTS-3 channels has been increased from 13 to 17, which can observe the atmospheric temperature and water vapor profiles from the surface to the lower stratosphere. In this study, two generations of MWTSs onboard FY-3D/3E were inter-calibrated by the Double Difference (DD) method to eliminate bias. The results showed that the biases of tropospheric channels were stable (within 1 K) and the biases of stratospheric channels were relatively large (over 2 K). In addition, the weighting functions of all MWTS channels varied with fields of view (FOVs) due to different optical paths, causing the brightness temperature (TB) observations to display strong scan-dependent features, i.e., the limb effect. This work used a limb correction method to remove scan-dependent patterns so that the underlying weather signals could be uncovered. After inter-calibration and limb correction, this work converted the TB observations from MWTS-2/3 onto a global gridded dataset at 0.5° × 0.5° latitudinal and longitudinal resolutions using a method of nested interpolation. Based on this research, more long-term FengYun series satellite climate datasets can be established in the future. Full article
Show Figures

Figure 1

Figure 1
<p>Weighting functions (WFs) of FY-3E MWTS-3 channels 1–17 (including 13 channels same as MWTS-2 shown in black curves), calculated using Radiative Transfer for TOVS (RTTOV) with the U.S. standard atmosphere profile at nadir.</p>
Full article ">Figure 2
<p>Monthly mean local equator-crossing times of FY-3D and FY-3E calculated from nadir observations on ascending (solid curves) and descending (dashed curves) nodes from January 2019 to June 2023.</p>
Full article ">Figure 3
<p>Weighting functions of MWTS-3 channel 7 at different scan angles of FOVs 1 (black), 15 (blue), 30 (green), and 45 (red) calculated by RTTOV using the U.S. standard atmosphere, 1976 (COESA 1976) profile.</p>
Full article ">Figure 4
<p>MWTS-3 channel 7 TBs observed on the ascending nodes on 1 January 2023.</p>
Full article ">Figure 5
<p>(<b>a</b>) Ascending nodes from FY-3D (blue curve) and FY-3E (red curve) swaths on 1 January 2023 shown under a north-polar stereographic projection. (<b>b</b>) All SNO points between FY-3D and FY-3E from March 2022 to June 2023.</p>
Full article ">Figure 6
<p>DD bias for all channels between MWTS-2 and MWTS-3 (mean value and standard deviation are indicated by black curves). Channel numbers in the figure are all from MWTS-3. (<b>a</b>) channels 3, 4, 5; (<b>b</b>) channels 7, 9, 10; (<b>c</b>) channels 11, 12, 13; (<b>d</b>) channels 14, 15, 16, 17.</p>
Full article ">Figure 7
<p>Inter-sensor DD bias (blue solid curve) with standard deviation (red dashed curve) of MWTS channels between FY-3D and FY-3E.</p>
Full article ">Figure 8
<p>Flowchart of inter-calibration between MWTS-2 and MWTS-3.</p>
Full article ">Figure 9
<p>MWTS-3 channel 7 mean TBs averaged at every scan position and in every 2° latitudinal band over global ocean. Data are from January 2023.</p>
Full article ">Figure 10
<p>Limb-corrected values of MWTS-3 channel 7 <math display="inline"><semantics> <mrow> <mo>[</mo> <msubsup> <mrow> <mi mathvariant="normal">T</mi> </mrow> <mrow> <mi mathvariant="normal">b</mi> </mrow> <mrow> <mrow> <mi>Ch</mi> <mn>7</mn> </mrow> </mrow> </msubsup> <mfenced separators="|"> <mrow> <mrow> <mi mathvariant="normal">i</mi> <mo>,</mo> <mo> </mo> <mi mathvariant="normal">j</mi> </mrow> </mrow> </mfenced> <mo> </mo> <mo>-</mo> <msubsup> <mrow> <mi mathvariant="normal">T</mi> </mrow> <mrow> <mrow> <mi mathvariant="normal">b</mi> <mo>,</mo> <mo> </mo> <mi>LC</mi> </mrow> </mrow> <mrow> <mrow> <mi>Ch</mi> <mn>7</mn> </mrow> </mrow> </msubsup> <mfenced separators="|"> <mrow> <mrow> <mi mathvariant="normal">i</mi> <mo>,</mo> <mo> </mo> <mi mathvariant="normal">j</mi> </mrow> </mrow> </mfenced> <mo>]</mo> </mrow> </semantics></math> TBs (black contour lines at 2 K intervals) as a function of scan position and latitude at 2° intervals, as well as scan variations in the global mean TBs of channel 3, 7, 10, and 12 before (solid-colored curves) and after (dashed colored curves) the limb correction. MWTS-3 observations are collected from January 2023.</p>
Full article ">Figure 11
<p>MWTS-3 channel 7 TBs observed on the ascending nodes on 15 January 2023 (<b>a</b>) before and (<b>b</b>) after the limb correction.</p>
Full article ">Figure 12
<p>Global distributions of gridded MWTS-3 channel-7 TB after nested interpolation on 15 January 2023.</p>
Full article ">Figure 13
<p>Flowchart for establishing multi-satellite MWTS gridded dataset.</p>
Full article ">
30 pages, 60239 KiB  
Article
Retrieval and Evaluation of Global Surface Albedo Based on AVHRR GAC Data of the Last 40 Years
by Shaopeng Li, Xiongxin Xiao, Christoph Neuhaus and Stefan Wunderle
Remote Sens. 2025, 17(1), 117; https://doi.org/10.3390/rs17010117 - 1 Jan 2025
Viewed by 835
Abstract
In this study, the global land surface albedo namely GAC43 was retrieved for the years 1979 to 2020 using Advanced Very High Resolution Radiometer (AVHRR) global area coverage (GAC) data onboard National Oceanic and Atmospheric Administration (NOAA) and Meteorological Operational (MetOp) satellites. We [...] Read more.
In this study, the global land surface albedo namely GAC43 was retrieved for the years 1979 to 2020 using Advanced Very High Resolution Radiometer (AVHRR) global area coverage (GAC) data onboard National Oceanic and Atmospheric Administration (NOAA) and Meteorological Operational (MetOp) satellites. We provide a comprehensive retrieval process of the GAC43 albedo, followed by a comprehensive assessment against in situ measurements and three widely used satellite-based albedo products, the third edition of the CM SAF cLoud, Albedo and surface RAdiation (CLARA-A3), the Copernicus Climate Change Service (C3S) albedo product, and MODIS BRDF/albedo product (MCD43). Our quantitative evaluations indicate that GAC43 demonstrates the best stability, with a linear trend of ±0.002 per decade at nearly all pseudo invariant calibration sites (PICS) from 1982 to 2020. In contrast, CLARA-A3 exhibits significant noise before the 2000s due to the limited availability of observations, while C3S shows substantial biases during the same period due to imperfect sensors intercalibrations. Extensive validation at globally distributed homogeneous sites shows that GAC43 has comparable accuracy to C3S, with an overall RMSE of approximately 0.03, but a smaller positive bias of 0.012. Comparatively, MCD43C3 shows the lowest RMSE (~0.023) and minimal bias, while CLARA-A3 displays the highest RMSE (~0.042) and bias (0.02). Furthermore, GAC43, CLARA-A3, and C3S exhibit overestimation in forests, with positive biases exceeding 0.023 and RMSEs of at least 0.028. In contrast, MCD43C3 shows negligible bias and a smaller RMSE of 0.015. For grasslands and shrublands, GAC43 and MCD43C3 demonstrate comparable estimation uncertainties of approximately 0.023, with close positive biases near 0.09, whereas C3S and CLARA-A3 exhibit higher RMSEs and biases exceeding 0.032 and 0.022, respectively. All four albedo products show significant RMSEs around 0.035 over croplands but achieve the highest estimation accuracy better than 0.020 over deserts. It is worth noting that significant biases are typically attributed to insufficient spatial representativeness of the measurement sites. Globally, GAC43 and C3S exhibit similar spatial distribution patterns across most land surface conditions, including an overestimation compared to MCD43C3 and an underestimation compared to CLARA-A3 in forested areas. In addition, GAC43, C3S, and CLARA-A3 estimate higher albedo values than MCD43C3 in low-vegetation regions, such as croplands, grasslands, savannas, and woody savannas. Besides the fact that the new GAC43 product shows the best stability covering the last 40 years, one has to consider the higher proportion of backup inversions before 2000. Overall, GAC43 offers a promising long-term and consistent albedo with good accuracy for future studies such as global climate change, energy balance, and land management policy. Full article
Show Figures

Figure 1

Figure 1
<p>Local solar times and solar zenith angles of equator observations for all AVHRR-carrying NOAA and MetOp satellites used to generate GAC43 albedo products as shown in (<b>a</b>,<b>b</b>), respectively. SZA &gt; 90° indicates night conditions.</p>
Full article ">Figure 2
<p>Globally distributed sites with homogeneous characteristics and corresponding land cover types defined by the IGBP from the MCD12C1 product. Purple squares located in the desert are used to evaluate temporal stability, while other sites are utilized for direct validations.</p>
Full article ">Figure 3
<p>Flowchart for this study.</p>
Full article ">Figure 4
<p>The performance of full inversion and full and backup inversion at various IGBP land cover types.</p>
Full article ">Figure 5
<p>The performance of the GAC43 albedo with full inversions at various land cover types, where panels (<b>a</b>–<b>h</b>) represent the land cover types of BSV, CRO, DBF, EBF, ENF, GRA, OSH and WSA, respectively. In the plots, the red solid line represents the 1:1 line, and the green dotted line and purple solid lines represent the limits of deviation ±0.02 and ±0.04, respectively.</p>
Full article ">Figure 6
<p>Google Earth <sup>TM</sup> images were used to visually illustrate the heterogeneity surrounding selected homogeneous sites representing various land cover types: (<b>a</b>) EBF, (<b>b</b>) BSV, (<b>c</b>) CRO and (<b>d</b>) GRA, as defined by the MCD12C1 IGBP classification. The red circle in each image denotes a radius of 2.5 km.</p>
Full article ">Figure 7
<p>Inter-comparison performance among four satellite-based albedo products. The top four subfigures (<b>a</b>–<b>d</b>) show the accuracy of all available matching samples between in situ measurements and estimated albedo values derived from satellite products, while the bottom four subfigures (<b>e</b>–<b>h</b>) give the performance of that using same samples.</p>
Full article ">Figure 8
<p>The performance of four satellite-based albedo products using same samples across various land surface types, evaluated in terms of (<b>a</b>) RMSE and (<b>b</b>) bias, respectively. The <span class="html-italic">x</span>-axis represents the land cover type classified as forest, grassland or shrublands, cropland, and desert, and corresponding available samples.</p>
Full article ">Figure 9
<p>The temporal performance of four satellite-based albedo products related to in situ measurements, and each subplot represents one case of different land cover surface, including (<b>a</b>) EBF, (<b>b</b>) ENF, (<b>c</b>) DBF, (<b>d</b>) GRA, and (<b>e</b>) CRO, respectively. The grey shaded areas depict situations with snow cover.</p>
Full article ">Figure 10
<p>Spatial distributions of GAC43 BSA in July 2013 are shown in subgraph (<b>a</b>), with corresponding differences from (<b>b</b>) CLARA-A3, (<b>c</b>) C3S, and (<b>d</b>) MCD43C3 in the same month, respectively.</p>
Full article ">Figure 11
<p>Percentage difference in BSA values between (<b>a</b>) GAC43 and CLARA-A3, (<b>b</b>) GAC43 and C3S, and (<b>c</b>) GAC43 and MCD43C3 in July 2013.</p>
Full article ">Figure 12
<p>The scattering plots between GAC43 BSA and (<b>a</b>) CLARA-A3 BSA, (<b>b</b>) C3S BSA, and (<b>c</b>) MCD43C3 BSA using all snow-free monthly pixels in July 2013, where the red lines indicate 1:1.</p>
Full article ">Figure 13
<p>The monthly BSA for the four satellite-based products across various land cover types in July 2013, where panels (<b>a</b>–<b>i</b>) represent the land cover types of CRO, DBF, DNF, EBF, ENF, GRA, MF, SAV and WSA, respectively. In the plots, the bottom values of each albedo product are the median of all corresponding land cover estimates. The top values match available samples.</p>
Full article ">Figure 14
<p>Monthly BSA from GAC43, MCD43C3, C3S, and CALRA-A3 at three randomly selected PICS sites: (<b>a</b>) Arabia 2, 20.19°N, 51.63°E; (<b>b</b>) Libya 3, 23.22°N, 23.23°E; and (<b>c</b>) Sudan 1, 22.11°N, 28.11°E, all characterized by BSV land surfaces as defined by IGBP.</p>
Full article ">Figure 15
<p>Box plots of the slope per decade for GAC43, CLARA-A3, C3S, and MCD43C3 at all PICS sites, where (<b>a</b>–<b>d</b>) represent the corresponding statistics during 1982–1990, 1991–2000, 2001–2010 and 2011–2020, respectively, and three dashed grey lines represent the 75%, 50%, and 25% quantiles. Red dotted lines indicate the horizontal line where slope is 0.</p>
Full article ">Figure 16
<p>Percentage of full inversions for the years 2004, 2008, 2012, and 2016 based on GAC43 (<b>top</b>) and MCD43A3 (<b>bottom</b>).</p>
Full article ">Figure 17
<p>Percentage of full inversions of GAC43 at various continents from 1979 to 2020.</p>
Full article ">Figure A1
<p>Spatial distributions of GAC43 BSA in July 2004 are shown in subgraph (<b>a</b>), with corresponding differences from (<b>b</b>) CLARA-A3, (<b>c</b>) C3S, and (<b>d</b>) MCD43C3 in the same month, respectively.</p>
Full article ">Figure A2
<p>Spatial distributions of GAC43 BSA in July 2008 are shown in subgraph (<b>a</b>), with corresponding differences from (<b>b</b>) CLARA-A3, (<b>c</b>) C3S, and (<b>d</b>) MCD43C3 in the same month, respectively.</p>
Full article ">Figure A3
<p>Spatial distributions of GAC43 BSA in July 2012 are shown in subgraph (<b>a</b>), with corresponding differences from (<b>b</b>) CLARA-A3, (<b>c</b>) C3S, and (<b>d</b>) MCD43C3 in the same month, respectively.</p>
Full article ">Figure A4
<p>Spatial distributions of GAC43 BSA in July 2016 are shown in subgraph (<b>a</b>), with corresponding differences from (<b>b</b>) CLARA-A3, (<b>c</b>) C3S, and (<b>d</b>) MCD43C3 in the same month, respectively.</p>
Full article ">Figure A5
<p>Percentages of full inversions for the years between 1979 and 2020 based on GAC43 data record.</p>
Full article ">
22 pages, 3002 KiB  
Review
Overview of Operational Global and Regional Ocean Colour Essential Ocean Variables Within the Copernicus Marine Service
by Vittorio E. Brando, Rosalia Santoleri, Simone Colella, Gianluca Volpe, Annalisa Di Cicco, Michela Sammartino, Luis González Vilas, Chiara Lapucci, Emanuele Böhm, Maria Laura Zoffoli, Claudia Cesarini, Vega Forneris, Flavio La Padula, Antoine Mangin, Quentin Jutard, Marine Bretagnon, Philippe Bryère, Julien Demaria, Ben Calton, Jane Netting, Shubha Sathyendranath, Davide D’Alimonte, Tamito Kajiyama, Dimitry Van der Zande, Quinten Vanhellemont, Kerstin Stelzer, Martin Böttcher and Carole Lebretonadd Show full author list remove Hide full author list
Remote Sens. 2024, 16(23), 4588; https://doi.org/10.3390/rs16234588 - 6 Dec 2024
Viewed by 862
Abstract
The Ocean Colour Thematic Assembly Centre (OCTAC) of the Copernicus Marine Service delivers state-of-the-art Ocean Colour core products for both global oceans and European seas, derived from multiple satellite missions. Since 2015, the OCTAC has provided global and regional high-level merged products that [...] Read more.
The Ocean Colour Thematic Assembly Centre (OCTAC) of the Copernicus Marine Service delivers state-of-the-art Ocean Colour core products for both global oceans and European seas, derived from multiple satellite missions. Since 2015, the OCTAC has provided global and regional high-level merged products that offer value-added information not directly available from space agencies. This is achieved by integrating observations from various missions, resulting in homogenized, inter-calibrated datasets with broader spatial coverage than single-sensor data streams. OCTAC enhanced continuously the basin-level accuracy of essential ocean variables (EOVs) across the global ocean and European regional seas, including the Atlantic, Arctic, Baltic, Mediterranean, and Black seas. From 2019 onwards, new EOVs have been introduced, focusing on phytoplankton functional groups, community structure, and primary production. This paper provides an overview of the evolution of the OCTAC catalogue from 2015 to date, evaluates the accuracy of global and regional products, and outlines plans for future product development. Full article
(This article belongs to the Special Issue Oceans from Space V)
Show Figures

Figure 1

Figure 1
<p>Overview of the OCTAC catalogue evolutions of the single-sensor and multisensor global and regional OC products from 2015 to 2024. The blue lines mark the timelines of each product type; covered basins are marked in green and listed under each line; satellite sensors are marked in black; spatial resolution of products/datasets is marked in blue. The red dots mark the dates of the MY reprocessing.</p>
Full article ">Figure 2
<p>Spatial coverage of the Sentinel-3 OLCI 300 m and Sentinel-2 MSI 100 m datasets. (<b>A</b>) All European regional seas and a 200 km strip from the coastline in the global product for Sentinel-3 OLCI. (<b>B</b>) A 20 km strip from the coastline for the European coastal waters covered in 5 days with Sentinel-2 MSI.</p>
Full article ">Figure 3
<p>OC sensors and high-resolution imagers adopted upstream in OCTAC processing chains. Timelines of legacy, and current and forthcoming (approved and planned) sensors are displayed (source CEOS): red identifies science OC missions, blue identifies operational OC missions, and brown identifies high-resolution/land imagers.</p>
Full article ">Figure 4
<p>Mediterranean Sea satellite CHL trend over the period 1997-2023, based on the CMEMS product OCEANCOLOUR_MED_BGC_L4_MY_009_144. (<b>A</b>) Time series and linear trend of monthly regional average satellite CHL: the monthly regional average (weighted by pixel area) time series is shown in gray, with the de-seasonalized time series in green and the linear trend in blue. (<b>B</b>) Map of satellite CHL trend, expressed in % per year, with positive trends in red and negative trends in blue.</p>
Full article ">Figure 5
<p>Time series (1998–2023) of SDG 14.1.1a Level 2 sub-indicator for European countries. The potential eutrophication values for European waters are based on CMEMS OC regional products aggregated over the EEZ for each country. AL: Albania, BE: Belgium, BG: Bulgaria, CY: Cyprus, DE: Germany, DK: Denmark, EE: Estonia, EL: Greece, ES: Spain, FI: Finland, FO: Faroe Islands, FR: France, GE: Georgia, GL: Greenland, HR: Croatia, IE: Ireland, IS: Iceland, IT: Italy, LT: Lithuania, LV: Lat- via, MC: Monaco, ME: Montenegro, MT: Malta, NL: Netherlands, NO: Norway, PL: Poland, PT: Portugal, RO: Romania, SE: Sweden, SI: Slovenia, UK: United Kingdom.</p>
Full article ">
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 635
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 2955
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
Cited by 1 | Viewed by 1188
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 1468
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 3 | Viewed by 2125
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 1940
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 7 | Viewed by 1784
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
Cited by 2 | Viewed by 1713
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 1805
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 2742
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 2332
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 Special Issue Advances in Air Pollutant Monitoring through Low-Cost Sensors)
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 ">
Back to TopTop