[go: up one dir, main page]

 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

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
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
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Article Types

Countries / Regions

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
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (6,555)

Search Parameters:
Keywords = Sentinel-2A

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
16 pages, 8588 KiB  
Article
A Novel Approach for Farmland Size Estimation in Small-Scale Agriculture Using Edge Counting and Remote Sensing
by Jingnan Du, Sucheng Xu, Jinshan Li, Jiakun Duan and Wu Xiao
Remote Sens. 2024, 16(16), 2981; https://doi.org/10.3390/rs16162981 - 14 Aug 2024
Abstract
Accurate and timely information on farmland size is crucial for agricultural development, resource management, and other related fields. However, there is currently no mature method for estimating farmland size in smallholder farming areas. This is due to the small size of farmland plots [...] Read more.
Accurate and timely information on farmland size is crucial for agricultural development, resource management, and other related fields. However, there is currently no mature method for estimating farmland size in smallholder farming areas. This is due to the small size of farmland plots in these areas, which have unclear boundaries in medium and high-resolution satellite imagery, and irregular shapes that make it difficult to extract complete boundaries using morphological rules. Automatic farmland mapping algorithms using remote sensing data also perform poorly in small-scale farming areas. To address this issue, this study proposes a farmland size evaluation index based on edge frequency (ECR). The algorithm utilizes the high temporal resolution of Sentinel-2 satellite imagery to compensate for its spatial resolution limitations. First, all Sentinel-2 images from one year are used to calculate edge frequencies, which can divide farmland areas into low-value farmland interior regions, medium-value non-permanent edges, and high-value permanent edges (PE). Next, the Otsu’s thresholding algorithm is iteratively applied twice to the edge frequencies to first extract edges and then permanent edges. The ratio of PE to cropland (ECR) is then calculated. Using the North China Plain and Northeast China Plain as study areas, and comparing with existing farmland size datasets, the appropriate estimation radius for ECR was determined to be 1600 m. The study found that the peak ECR value for the Northeast China Plain was 0.085, and the peak value for the North China Plain was 0.105. The overall distribution was consistent with the reference dataset. Full article
(This article belongs to the Special Issue Advances in Remote Sensing for Crop Monitoring and Food Security)
Show Figures

Figure 1

Figure 1
<p>The red borders outline the six major grain-producing provinces in China, which serve as our study area. From top to bottom, left to right, they are Inner Mongolia and Liaoning Province in the Northeast China Plain, and Shandong Province, Henan Province, Hubei Province, and Anhui Province in the North China Plain. There are a total of 1792 sample points from the Geo-Wiki plot size dataset that fall within the study area. Both color and size are used to display the points based on plot size for better visualization.</p>
Full article ">Figure 2
<p>Technical flowchart of this study. Abbreviations: NDVI stands for normalized difference vegetation index; Otsu, Otsu binary segmentation algorithm; ECR, edge cropland ratio [<a href="#B6-remotesensing-16-02981" class="html-bibr">6</a>].</p>
Full article ">Figure 3
<p>The figure illustrates the Edge count generation process within a 3200 m radius of a sample point (Sample ID: 962700, Latitude: 45.417702, Longitude: 121.545998) located in the Inner Mongolia Autonomous Region. Edge count represents the number of times each pixel is marked as an edge. In the grayscale image, brighter areas indicate higher counts, while pure black areas represent non-farmland regions.</p>
Full article ">Figure 4
<p>(<b>a</b>) shows the buffer zones of different radii for the sample point, ranging from 200 m to 3200 m. (<b>b</b>) displays the edge count images obtained at different radii. The edge count result for the 3200 m radius is shown in (<b>c</b>), with non-farmland areas set as transparent. (<b>d</b>–<b>f</b>), respectively, represent the extracted edges, permanent edges, and edge frequency distribution of (<b>c</b>).</p>
Full article ">Figure 5
<p>The boxplot illustrates the distribution of estimated ECR values at different radii. The numerical statistics are shown in the right figure.</p>
Full article ">Figure 6
<p>The distribution of ECR for each parcel size group is shown in separate figures with different radii.</p>
Full article ">Figure 7
<p>The line chart on the left shows the probability density distribution of ECRs at different radii, while the bar chart on the right shows the number of Field size labels at each level, categorized by province.</p>
Full article ">Figure 8
<p>The figure presents three samples of different farmland sizes in separate columns, from left to right: XL, M, and XS. The three rows from top to bottom are: Google Satellite basemap and sample’s metadata, edge count image and grayscale histogram (upper right corner), and identified permanent edges (in red).</p>
Full article ">Figure 9
<p>The probability density curves of the 1600 m ECRs are displayed, grouped by parcel size. The three vertical dashed lines in the left figure represent the optimal thresholds for classifying farmland size based on ECRs, determined using Spearman’s rank correlation coefficient. The confusion matrix on the right shows the comparison between our ECR-predicted farmland sizes and the manually interpreted labeling results.</p>
Full article ">
26 pages, 8634 KiB  
Article
New Insights on the Information Content of the Normalized Difference Vegetation Index Sentinel-2 Time Series for Assessing Vegetation Dynamics
by César Sáenz, Víctor Cicuéndez, Gabriel García, Diego Madruga, Laura Recuero, Alfonso Bermejo-Saiz, Javier Litago, Ignacio de la Calle and Alicia Palacios-Orueta
Remote Sens. 2024, 16(16), 2980; https://doi.org/10.3390/rs16162980 - 14 Aug 2024
Abstract
The Sentinel-2 NDVI time series information content from 2017 to 2023 at a 10 m spatial resolution was evaluated based on the NDVI temporal dependency in five scenarios in central Spain. First, time series were interpolated and then filtered using the Savitzky–Golay, Fast [...] Read more.
The Sentinel-2 NDVI time series information content from 2017 to 2023 at a 10 m spatial resolution was evaluated based on the NDVI temporal dependency in five scenarios in central Spain. First, time series were interpolated and then filtered using the Savitzky–Golay, Fast Fourier Transform, Whittaker, and Maximum Value filters. Temporal dependency was assessed using the Q-Ljung-Box and Fisher’s Kappa tests, and similarity between raw and filtered time series was assessed using Correlation Coefficient and Root Mean Square Error. An Interpolating Efficiency Indicator (IEI) was proposed to summarize the number and temporal distribution of low-quality observations. Type of climate, atmospheric disturbances, land cover dynamics, and management were the main sources of variability in five scenarios: (1) rainfed wheat and barley presented high short-term variability due to clouds (lower IEI in winter and spring) during the growing cycle and high interannual variability due to precipitation; (2) maize showed stable summer cycles (high IEI) and low interannual variability due to irrigation; (3) irrigated alfalfa was cut five to six times during summer, resulting in specific intra-annual variability; (4) beech forest showed a strong and stable summer cycle, despite the short-term variability due to clouds (low IEI); and (5) evergreen pine forest had a highly variable growing cycle due to fast responses to temperature and precipitation through the year and medium IEI values. Interpolation after removing non-valid observations resulted in an increase in temporal dependency (Q-test), particularly a short term in areas with low IEI values. The information improvement made it possible to identify hidden periodicities and trends using the Fisher’s Kappa test. The SG filter showed high similarity values and weak influence on dynamics, while the MVF showed an overestimation of the NDVI values. Full article
(This article belongs to the Special Issue Crops and Vegetation Monitoring with Remote/Proximal Sensing II)
Show Figures

Figure 1

Figure 1
<p>Image in true color from © Google Earth (Landsat/Copernicus). (<b>a</b>) Map of Spain and (<b>b</b>) study area.</p>
Full article ">Figure 2
<p>Sentinel-2 image processing flowchart of time series filtering.</p>
Full article ">Figure 3
<p>Spatial distribution of the vegetation species in the study area.</p>
Full article ">Figure 4
<p>Valid observations according to the SCL band of Sentinel-2. (<b>a</b>) Spatial distribution in percentage and (<b>b</b>) number of pixels for each level of valid observations.</p>
Full article ">Figure 5
<p>(<b>a</b>) Length of the longest gap within time series; (<b>b</b>) season in which the longest gap occurs.</p>
Full article ">Figure 6
<p>Interpolating Efficiency Indicator for Tile 30TUN: (<b>a</b>) spatial distribution and (<b>b</b>) frequency distribution.</p>
Full article ">Figure 7
<p>Time series of pure pixels for a wheat (<b>a</b>), barley (<b>c</b>), maize (<b>e</b>), and alfalfa (<b>g</b>) crop. Time series of the interpolated NDVI (<b><span style="color:#00B0F0">−</span></b>) and its four filters: SG (<b><span style="color:red">−</span></b>), FFT (<b><span style="color:#92D050">−</span></b>), WHI (<b>−</b>), and MVF (<b><span style="color:yellow">−</span></b>), and (<b>b</b>,<b>d</b>,<b>f</b>,<b>h</b>) the average year, i.e., Buys-Ballot. The colors of each filter and interpolation are the same as those of the entire time series.</p>
Full article ">Figure 8
<p>Time series of pure pixels for a beech forest (<b>a</b>) and a pine forest (<b>c</b>). Time series of the interpolated NDVI (<b><span style="color:#00B0F0">−</span></b>) and its four filters: SG (<b><span style="color:red">−</span></b>), FFT (<b><span style="color:#92D050">−</span></b>), WHI (<b>−</b>), and MVF (<b><span style="color:yellow">−</span></b>), and (<b>b</b>,<b>d</b>) the average year, i.e., Buys-Ballot. The colors of each filter and interpolation are the same as those of the entire time series.</p>
Full article ">Figure 9
<p>Representation of the evolution at the image level of the interpolation and filtering processes of the Sentinel-2 NDVI time series using the Whittaker filter. (<b>a</b>) NDVI image (10/06/20) without interpolation, (<b>b</b>) NDVI image (15/06/20) without interpolation, (<b>c</b>) NDVI image (20/06/20) without interpolation, (<b>d</b>) NDVI image (15/06/20) interpolated, and (<b>e</b>) NDVI image (15/06/20) filtered using the Whittaker filter.</p>
Full article ">Figure 10
<p>Values of the Q-test for the short term (lags 1, 2, 3, 4, 5, 6, and 7).</p>
Full article ">Figure 11
<p>Values of the Q-test at lag 36 (6 months), lag 73 (one year), and lag 146 (two years).</p>
Full article ">Figure 12
<p>Average Fk value for period 73 (one year) for each vegetation type (left axis) and percentage of increase after interpolation and filtering (right axis).</p>
Full article ">Figure 13
<p>Time series of a pixel declared as alfalfa implemented during the first year.</p>
Full article ">
40 pages, 19379 KiB  
Article
Evaluation of Sentinel-5P TROPOMI Methane Observations at Northern High Latitudes
by Hannakaisa Lindqvist, Ella Kivimäki, Tuomas Häkkilä, Aki Tsuruta, Oliver Schneising, Michael Buchwitz, Alba Lorente, Mari Martinez Velarte, Tobias Borsdorff, Carlos Alberti, Leif Backman, Matthias Buschmann, Huilin Chen, Darko Dubravica, Frank Hase, Pauli Heikkinen, Tomi Karppinen, Rigel Kivi, Erin McGee, Justus Notholt, Kimmo Rautiainen, Sébastien Roche, William Simpson, Kimberly Strong, Qiansi Tu, Debra Wunch, Tuula Aalto and Johanna Tamminenadd Show full author list remove Hide full author list
Remote Sens. 2024, 16(16), 2979; https://doi.org/10.3390/rs16162979 - 14 Aug 2024
Abstract
The Arctic and boreal regions are experiencing a rapid increase in temperature, resulting in a changing cryosphere, increasing human activity, and potentially increasing high-latitude methane emissions. Satellite observations from Sentinel-5P TROPOMI provide an unprecedented coverage of a column-averaged dry-air mole fraction of methane [...] Read more.
The Arctic and boreal regions are experiencing a rapid increase in temperature, resulting in a changing cryosphere, increasing human activity, and potentially increasing high-latitude methane emissions. Satellite observations from Sentinel-5P TROPOMI provide an unprecedented coverage of a column-averaged dry-air mole fraction of methane (XCH4) in the Arctic, compared to previous missions or in situ measurements. The purpose of this study is to support and enhance the data used for high-latitude research through presenting a systematic evaluation of TROPOMI methane products derived from two different processing algorithms: the operational product (OPER) and the scientific product (WFMD), including the comparison of recent version changes of the products (OPER, OPER rpro, WFMD v1.2, and WFMD v1.8). One finding is that OPER rpro yields lower XCH4 than WFMD v1.8, the difference increasing towards the highest latitudes. TROPOMI product differences were evaluated with respect to ground-based high-latitude references, including four Fourier Transform Spectrometer in the Total Carbon Column Observing Network (TCCON) and five EM27/SUN instruments in the Collaborative Carbon Column Observing Network (COCCON). The mean TROPOMI–TCCON GGG2020 daily median XCH4 difference was site-dependent and varied for OPER rpro from −0.47 ppb to 22.4 ppb, and for WFMD v1.8 from 1.2 ppb to 19.4 ppb with standard deviations between 13.0 and 20.4 ppb and 12.5–15.0 ppb, respectively. The TROPOMI–COCCON daily median XCH4 difference varied from −26.5 ppb to 5.6 ppb for OPER rpro, with a standard deviation of 14.0–28.7 ppb, and from −5.0 ppb to 17.2 ppb for WFMD v1.8, with a standard deviation of 11.5–13.0 ppb. Although the accuracy and precision of both TROPOMI products are, on average, good compared to the TCCON and COCCON, a persistent seasonal bias in TROPOMI XCH4 (high values in spring; low values in autumn) is found for OPER rpro and is reflected in the higher standard deviation values. A systematic decrease of about 7 ppb was found between TCCON GGG2014 and GGG2020 product update highlighting the importance of also ensuring the reliability of ground-based retrievals. Comparisons to atmospheric profile measurements with AirCore carried out in Sodankylä, Northern Finland, resulted in XCH4 differences comparable to or smaller than those from ground-based remote sensing. Full article
(This article belongs to the Special Issue Remote Sensing of Greenhouse Gas Emissions II)
Show Figures

Figure 1

Figure 1
<p>Locations of the high-latitude TCCON, COCCON, and AirCore sites overlaid on a map of the regional permafrost extent re-gridded from ESA CCI Permafrost data. Regions with &gt;90% permafrost extent are considered continuous permafrost.</p>
Full article ">Figure 2
<p>Effect of the COCCON prior correction on TROPOMI OPER rpro (blue) and WFMD v1.8 (red) daily median XCH<sub>4</sub> at (<b>a</b>) Kiruna, Sweden, (<b>b</b>) SN039 and (<b>c</b>) SN122 in Sodankylä, Finland, (<b>d</b>) Fairbanks, Alaska, USA, and (<b>e</b>) St. Petersburg, Russia.</p>
Full article ">Figure 3
<p>Spatial variability of monthly-averaged total column methane for TROPOMI OPER rpro (<b>left</b>) and WFMD v1.8 (<b>right</b>) for April, August, and October in 2020. The grid size is 0.25° × 0.2°.</p>
Full article ">Figure 4
<p>Spatial variability of the difference in monthly-averaged total column methane for TROPOMI OPER rpro and WFMD v1.8 in 2020. The grid size is 0.25° × 0.2°.</p>
Full article ">Figure 5
<p>Temporal dependence of (<b>a</b>) TROPOMI OPER rpro–WFMD v1.8 XCH<sub>4</sub> difference, (<b>b</b>) TROPOMI OPER–OPER rpro XCH<sub>4</sub> difference, and (<b>c</b>) TROPOMI WFMD v1.2–WFMD v1.8 XCH<sub>4</sub> difference, colored based on 5-degree latitude bands north of 50°N. Gaps represent missing data.</p>
Full article ">Figure 6
<p>Averages of (<b>a</b>) TROPOMI OPER rpro–WFMD v1.8 XCH<sub>4</sub> difference, (<b>b</b>) TROPOMI OPER–OPER rpro XCH<sub>4</sub> difference, and (<b>c</b>) TROPOMI WFMD v1.2–WFMD v1.8 XCH<sub>4</sub> difference for different latitude bands. The error bars denote the standard deviation. Colored circles denote the mean difference of the satellite product to the ground-based TCCON instrument (GGG2020 retrieval) located in that latitude band.</p>
Full article ">Figure 7
<p>Seasonal dependence of the number of TROPOMI observations over different permafrost regions for (<b>a</b>) OPER, (<b>b</b>) OPER rpro, (<b>c</b>) WFMD v1.2, and (<b>d</b>) WFMD v1.8. The colors refer to different permafrost classes, determined using the ESA Permafrost CCI Level 4 product.</p>
Full article ">Figure 8
<p>Daily median XCH<sub>4</sub> from TCCON GGG2020 (blue) and co-located TROPOMI (red) (<b>a</b>) OPER rpro and (<b>c</b>) WFMD v1.8 retrievals, and TROPOMI–TCCON differences (<b>b</b>,<b>d</b>) at Eureka, Canada. The dashed line in (<b>b</b>,<b>d</b>) shows the mean difference at the site and shaded area the standard deviation of the mean.</p>
Full article ">Figure 9
<p>Daily median XCH<sub>4</sub> from TCCON GGG2020 (blue) and co-located TROPOMI (red) (<b>a</b>) OPER rpro and (<b>c</b>) WFMD v1.8 retrievals, and TROPOMI–TCCON differences (<b>b</b>,<b>d</b>) at Ny-Ålesund, Norway. The dashed line in (<b>b</b>,<b>d</b>) shows the mean difference at the site and shaded area the standard deviation of the mean.</p>
Full article ">Figure 10
<p>Daily median XCH<sub>4</sub> from TCCON GGG2020 (blue) and co-located TROPOMI (red) (<b>a</b>) OPER rpro and (<b>c</b>) WFMD v1.8 retrievals, and TROPOMI–TCCON differences (<b>b</b>,<b>d</b>) at Sodankylä, Finland. The dashed line in (<b>b</b>,<b>d</b>) shows the mean difference at the site and shaded area the standard deviation of the mean.</p>
Full article ">Figure 11
<p>Daily median XCH<sub>4</sub> from TCCON GGG2020 (blue) and co-located TROPOMI (red) (<b>a</b>) OPER rpro and (<b>c</b>) WFMD v1.8 retrievals, and TROPOMI–TCCON differences (<b>b</b>,<b>d</b>) at East Trout Lake, Canada. The dashed line in (<b>b</b>,<b>d</b>) shows the mean difference at the site and shaded area the standard deviation of the mean.</p>
Full article ">Figure 12
<p>Comparison of TROPOMI–TCCON mean XCH<sub>4</sub> differences for TCCON GGG2014 and TCCON GGG2020, considering TROPOMI (<b>a</b>) OPER rpro and (<b>b</b>) WFMD v1.8 products. The TCCON data have been selected to include only the days for which both processing versions yield good-quality retrievals, so their temporal sampling is identical.</p>
Full article ">Figure 13
<p>Daily median XCH<sub>4</sub> from COCCON (red) and co-located TROPOMI (blue) (<b>a</b>) OPER rpro and (<b>c</b>) WFMD v1.8 retrievals, and TROPOMI–COCCON differences (<b>b</b>,<b>d</b>) at Kiruna, Sweden. The dashed line in (<b>b</b>,<b>d</b>) shows the mean differences and the shaded area denotes one standard deviation from the mean.</p>
Full article ">Figure 14
<p>Daily median XCH<sub>4</sub> from COCCON (red) and co-located TROPOMI (blue) (<b>a</b>) OPER rpro and (<b>c</b>) WFMD v1.8 retrievals, and TROPOMI–COCCON differences (<b>b</b>,<b>d</b>) at Sodankylä, Finland (SN039). The dashed line in (<b>b</b>,<b>d</b>) shows the mean differences and the shaded area denotes one standard deviation from the mean.</p>
Full article ">Figure 15
<p>Daily median XCH<sub>4</sub> from COCCON (red) and co-located TROPOMI (blue) (<b>a</b>) OPER rpro and (<b>c</b>) WFMD v1.8 retrievals, and TROPOMI–COCCON differences (<b>b</b>,<b>d</b>) at Sodankylä, Finland (SN122). The dashed line in (<b>b</b>,<b>d</b>) shows the mean differences and the shaded area denotes one standard deviation from the mean.</p>
Full article ">Figure 16
<p>Daily median XCH<sub>4</sub> from COCCON (green) and co-located TROPOMI (orange) (<b>a</b>) OPER rpro and (<b>c</b>) WFMD v1.8 retrievals, and TROPOMI–COCCON differences (<b>b</b>,<b>d</b>) at Fairbanks, Alaska, USA. The dashed line in (<b>b</b>,<b>d</b>) shows the mean differences and the shaded area denotes one standard deviation from the mean.</p>
Full article ">Figure 17
<p>Daily median XCH<sub>4</sub> from COCCON (ref) and co-located TROPOMI (blue) (<b>a</b>) OPER rpro and (<b>c</b>) WFMD v1.8 retrievals, and TROPOMI–COCCON differences (<b>b</b>,<b>d</b>) at St. Petersburg, Russia. The dashed line in (<b>b</b>,<b>d</b>) shows the mean differences and the shaded area denotes one standard deviation from the mean.</p>
Full article ">Figure 18
<p>Comparison of measured AirCore profiles (black curves), TCCON GGG2020 prior profiles (red), and TROPOMI prior profiles: (<b>a</b>) TROPOMI OPER rpro (blue), and (<b>c</b>) TROPOMI WFMD v1.8 (orange) prior profiles. The TROPOMI prior profiles are also scaled (dashed lines) so that the profile corresponds to each retrieved XCH<sub>4</sub>. The difference of the scaled profile compared to the measured concentration is also shown: (<b>b</b>) OPER rpro–AirCore, and (<b>d</b>) WFMD v1.8–AirCore.</p>
Full article ">Figure 19
<p>Collection of the AirCore XCH<sub>4</sub> results with co-located TCCON GGG2020 and TROPOMI observations.</p>
Full article ">Figure 20
<p>Collection of the TROPOMI high-latitude evaluation results obtained in this paper. The bar plots show the average difference of TROPOMI (OPER rpro or WFMD v1.8 product) and ground-based references, and the error bars depict the standard deviation. The references from left to right are TCCON stations at East Trout Lake, Sodankylä, Ny-Ålesund, and Eureka, then COCCON stations at St. Petersburg, Fairbanks, Sodankylä (SN039 and SN122), and Kiruna, and AirCore measurements at Sodankylä. It should be noted that the temporal sampling differs between the references.</p>
Full article ">Figure A1
<p>TCCON GGG2020 (blue) and co-located TROPOMI daily median XCH<sub>4</sub> (red) with (<b>a</b>) OPER and (<b>c</b>) WFMD v1.2 retrievals, and TROPOMI−TCCON differences (<b>b</b>,<b>d</b>) at Eureka (Canada) TCCON site. The dashed line in (<b>b</b>,<b>d</b>) show the mean difference at the site and shaded area the standard deviation of the mean.</p>
Full article ">Figure A2
<p>TCCON GGG2020 (blue) and co-located TROPOMI daily median XCH<sub>4</sub> (red) with (<b>a</b>) OPER and (<b>c</b>) WFMD v1.2 retrievals, and TROPOMI−TCCON differences (<b>b</b>,<b>d</b>) at Ny-Ålesund (Norway) TCCON site. The dashed line in (<b>b</b>,<b>d</b>) show the mean difference at the site and shaded area the standard deviation of the mean.</p>
Full article ">Figure A3
<p>TCCON GGG2020 (blue) and co-located TROPOMI daily median XCH<sub>4</sub> (red) with (<b>a</b>) OPER and (<b>c</b>) WFMD v1.2 retrievals, and TROPOMI−TCCON differences (<b>b</b>,<b>d</b>) at Sodankylä (Finland) TCCON site. The dashed line in (<b>b</b>,<b>d</b>) show the mean difference at the site and shaded area the standard deviation of the mean.</p>
Full article ">Figure A4
<p>TCCON GGG2020 (blue) and co-located TROPOMI daily median XCH<sub>4</sub> (red) with (<b>a</b>) OPER and (<b>c</b>) WFMD v1.2 retrievals, and TROPOMI−TCCON differences (<b>b</b>,<b>d</b>) at East Trout Lake (Canada) TCCON site. The dashed line in (<b>b</b>,<b>d</b>) show the mean difference at the site and shaded area the standard deviation of the mean.</p>
Full article ">
21 pages, 19442 KiB  
Article
Pasture Quality Assessment through NDVI Obtained by Remote Sensing: A Validation Study in the Mediterranean Silvo-Pastoral Ecosystem
by João Serrano, Shakib Shahidian, Luís Paixão, José Marques da Silva and Luís Lorenzo Paniágua
Agriculture 2024, 14(8), 1350; https://doi.org/10.3390/agriculture14081350 - 13 Aug 2024
Viewed by 187
Abstract
Monitoring the evolution of pasture availability and quality throughout the growing season is the basis of grazing management in extensive Mediterranean livestock systems. Remote sensing (RS) is an innovative tool that, among many other applications, is being developed for detailed spatial and temporal [...] Read more.
Monitoring the evolution of pasture availability and quality throughout the growing season is the basis of grazing management in extensive Mediterranean livestock systems. Remote sensing (RS) is an innovative tool that, among many other applications, is being developed for detailed spatial and temporal pasture quality assessment. The aim of the present study is to evaluate the potential of satellite images (Sentinel-2) to assess indicators of pasture quality (pasture moisture content, PMC, crude protein, CP and neutral detergent fiber, NDF) using the normalized difference vegetation index (NDVI). Field measurements were conducted over three years at eight representative fields of the biodiversity and variability of dryland pastures in Portugal. A total of 656 georeferenced pasture samples were collected and processed in the laboratory. The results show a significant correlation between pasture quality parameters (PMC, CP and NDF) obtained in standard laboratory methods and NDVI satellite-derived data (R2 of 0.72, 0.75, and 0.50, respectively). The promising findings obtained in this large-scale validation study (three years and eight fields) encourage further research (i) to test and develop other vegetation indexes for monitoring pasture nutritive value; (ii) to extend this research to pastures of the other Mediterranean countries, building large and representative datasets and developing more robust and accurate monitoring models based on freely available Sentinel-2 images; (iii) to implement an extension program for agricultural managers to popularize the use of these technological tools as the basis of grazing and pasture management. Full article
(This article belongs to the Section Digital Agriculture)
Show Figures

Figure 1

Figure 1
<p>Location of the eight experimental fields (seven in Portugal and one in Spain).</p>
Full article ">Figure 2
<p>Location of the eight sampling areas in each experimental field and illustrative photography: (<b>a</b>) AZI; (<b>b</b>) CUB; (<b>c</b>) GRO; (<b>d</b>) MIT.</p>
Full article ">Figure 3
<p>Location of the eight sampling areas in each experimental field and illustrative photography: (<b>a</b>) MUR; (<b>b</b>) PAD; (<b>c</b>) QF; (<b>d</b>) TAP.</p>
Full article ">Figure 4
<p>Flowchart of the methodology used.</p>
Full article ">Figure 5
<p>Evolution of mean values of pasture crude protein (CP), neutral detergent fiber (NDF) and moisture content (PMC), in MIT experimental field (11 pasture collection campaigns in three years).</p>
Full article ">Figure 6
<p>Thermo-pluviometric diagram of the Meteorological Station of Mitra (Évora, Portugal) between July 2018 and June 2021.</p>
Full article ">Figure 7
<p>Regression analysis between NDVI obtained by remote sensing and pasture moisture content (PMC).</p>
Full article ">Figure 8
<p>Regression analysis between NDVI and pasture crude protein (CP).</p>
Full article ">Figure 9
<p>Regression analysis between NDVI obtained by remote sensing and pasture neutral detergent fiber (NDF).</p>
Full article ">
10 pages, 771 KiB  
Brief Report
Prevalence of Co-Infections in Primary Care Patients with Medically Attended Acute Respiratory Infection in the 2022/2023 Season
by Maja Sočan, Katarina Prosenc and Maja Mrzel
Viruses 2024, 16(8), 1289; https://doi.org/10.3390/v16081289 - 13 Aug 2024
Viewed by 213
Abstract
In the post-pandemic period, an endemic circulation of respiratory viruses has been re-established. Respiratory viruses are co-circulating with SARS-CoV-2. We performed a retrospective analysis of co-infections in primary care patients with medically attended acute respiratory infections (MAARI) who consulted from week 40/2022 to [...] Read more.
In the post-pandemic period, an endemic circulation of respiratory viruses has been re-established. Respiratory viruses are co-circulating with SARS-CoV-2. We performed a retrospective analysis of co-infections in primary care patients with medically attended acute respiratory infections (MAARI) who consulted from week 40/2022 to week 39/2023 and were tested for a panel of respiratory viruses. Out of 2099 samples tested, 1260 (60.0%) were positive for one virus. In 340 samples, co-infection was detected: two viruses in 281 (13.4%), three viruses in 51 (2.4%), and four viruses in eight (0.4%) samples. Respiratory viruses co-infected the patients with MAARI at very different rates. The lowest rates of co-infections were confirmed for influenza B (13.8%) and influenza A (22.9%) and the highest for human bocaviruses (84.0%) and human parechoviruses (82.1%). Co-infections were detected in 28.2% of SARS-CoV-2 positive samples. SARS-CoV-2 has never been co-infected with influenza B virus, enterovirus or adenovirus, although the latter was found as a co-infecting virus with all other respiratory viruses tested. The rate of co-infections decreased significantly with increasing age (p-value 0.000), and no difference was found regarding gender (p-value 0.672). It is important to understand the epidemiology of respiratory co-infections for prevention and management decisions in patients with MAARI. Full article
Show Figures

Figure 1

Figure 1
<p>The number of ARI cases infected with one respiratory virus or in combination with other respiratory viruses.</p>
Full article ">Figure 2
<p>The number of monthly detections of respiratory viruses in the 2022/2023 season in Slovenia.</p>
Full article ">
21 pages, 28441 KiB  
Article
MSSFNet: A Multiscale Spatial–Spectral Fusion Network for Extracting Offshore Floating Raft Aquaculture Areas in Multispectral Remote Sensing Images
by Haomiao Yu, Yingzi Hou, Fangxiong Wang, Junfu Wang, Jianfeng Zhu and Jianke Guo
Sensors 2024, 24(16), 5220; https://doi.org/10.3390/s24165220 (registering DOI) - 12 Aug 2024
Viewed by 250
Abstract
Accurately extracting large-scale offshore floating raft aquaculture (FRA) areas is crucial for supporting scientific planning and precise aquaculture management. While remote sensing technology offers advantages such as wide coverage, rapid imaging, and multispectral capabilities for FRA monitoring, the current methods face challenges in [...] Read more.
Accurately extracting large-scale offshore floating raft aquaculture (FRA) areas is crucial for supporting scientific planning and precise aquaculture management. While remote sensing technology offers advantages such as wide coverage, rapid imaging, and multispectral capabilities for FRA monitoring, the current methods face challenges in terms of establishing spatial–spectral correlations and extracting multiscale features, thereby limiting their accuracy. To address these issues, we propose an innovative multiscale spatial–spectral fusion network (MSSFNet) designed specifically for extracting offshore FRA areas from multispectral remote sensing imagery. MSSFNet effectively integrates spectral and spatial information through a spatial–spectral feature extraction block (SSFEB), significantly enhancing the accuracy of FRA area identification. Additionally, a multiscale spatial attention block (MSAB) captures contextual information across different scales, improving the ability to detect FRA areas of varying sizes and shapes while minimizing edge artifacts. We created the CHN-YE7-FRA dataset using Sentinel-2 multispectral remote sensing imagery and conducted extensive evaluations. The results showed that MSSFNet achieved impressive metrics: an F1 score of 90.76%, an intersection over union (IoU) of 83.08%, and a kappa coefficient of 89.75%, surpassing those of state-of-the-art methods. The ablation results confirmed that the SSFEB and MSAB modules effectively enhanced the FRA extraction accuracy. Furthermore, the successful practical applications of MSSFNet validated its generalizability and robustness across diverse marine environments. These findings highlight the performance of MSSFNet in both experimental and real-world scenarios, providing reliable, precise FRA area monitoring. This capability provides crucial data for scientific planning and environmental protection purposes in coastal aquaculture zones. Full article
(This article belongs to the Special Issue Application of Satellite Remote Sensing in Geospatial Monitoring)
Show Figures

Figure 1

Figure 1
<p>Geographical locations of and images derived from the study areas ((<b>A</b>). Changhai County in Liaoning Province, (<b>B</b>). Jinshitan Bay in Liaoning Province, (<b>C</b>). Rongcheng Bay in Shandong Province, (<b>D</b>). Haizhou Bay in Jiangsu Province, (<b>E</b>). Dayu Bay in Jiangsu Province, (<b>F</b>). Sansha Bay in Fujian Province, and (<b>G</b>). Zhaoan Bay in Fujian Province).</p>
Full article ">Figure 2
<p>Structure of MSSFNet.</p>
Full article ">Figure 3
<p>Structure of the SSFEB.</p>
Full article ">Figure 4
<p>Structure of the MSAB.</p>
Full article ">Figure 5
<p>Visual comparison between MSSFNet and other mainstream networks (Subfigures (<b>a</b>–<b>g</b>) represent various scenarios. The white area in the prediction result map indicates the FRA and the black area indicates the background area. The red rectangle indicates the false detection region, the yellow rectangle indicates the missed detection region, and the green rectangle indicates the overlapping region).</p>
Full article ">Figure 6
<p>Comparison among the attention visualizations produced for the modules integrated into MSSFNet (Subfigures (<b>a</b>–<b>f</b>) represent various scenarios. The results of the extraction are plotted from blue to red to indicate low to high levels of attention, respectively).</p>
Full article ">Figure 7
<p>The FRA distribution areas of Changhai County, Haizhou Bay, and Sansha Bay in 2024 (the white areas in the figure represent the FRA regions, with T1 denoting Changhai County, T2 signifying Haizhou Bay, and T3 denoting Sansha Bay).</p>
Full article ">Figure 8
<p>Geographical locations and FRA extraction results outside the sampling areas (T1 denotes Longwangtang Bay in Liaoning Province, T2 represents Zhangjia Bay in Shandong Province, T3 signifies Taizhou Bay in Zhejiang Province, and T4 denotes Qingchuan Bay in Fujian Province). The white areas in the figure represent the FRA extraction results.</p>
Full article ">
21 pages, 42176 KiB  
Article
Application of Getis-Ord Correlation Index (Gi) for Burned Area Detection Improvement in Mediterranean Ecosystems (Southern Italy and Sardinia) Using Sentinel-2 Data
by Antonio Lanorte, Gabriele Nolè and Giuseppe Cillis
Remote Sens. 2024, 16(16), 2943; https://doi.org/10.3390/rs16162943 - 12 Aug 2024
Viewed by 284
Abstract
This study collects the results obtained using the Getis-Ord local spatial autocorrelation index (Gi) with the aim of improving the classification of burned area detection maps generated from spectral indices (i.e., dNBR index) derived from Sentinel-2 satellite data. Therefore, the work proposes an [...] Read more.
This study collects the results obtained using the Getis-Ord local spatial autocorrelation index (Gi) with the aim of improving the classification of burned area detection maps generated from spectral indices (i.e., dNBR index) derived from Sentinel-2 satellite data. Therefore, the work proposes an adaptive thresholding approach that also includes the application of a similarity index (Sorensen–Dice Similarity Index) with the aim of adaptively correcting classification errors (false-positive burned pixels) related to the spectral response of burned/unburned areas. In this way, two new indices derived from the application of the Getis-Ord local autocorrelation analysis were created to test their effectiveness. Three wildfire events were considered, two of which occurred in Southern Italy in the summer of 2017 and one in Sardinia in the summer of 2019. The accuracy assessment analysis was carried out using the CEMS (Copernicus Emergency Management Service) on-demand maps. The results show the remarkable performance of the two new indices in terms of their ability to reduce the false positives generated by dNBR. In the three sites considered, the false-positive reduction percentage was around 95–96%. The proposed approach seems to be adaptable to different vegetation contexts, and above all, it could be a useful tool for mapping burned areas to support post-fire management activities. Full article
Show Figures

Figure 1

Figure 1
<p>Location and perimeter of the burned areas analysed as provided by CEMS.</p>
Full article ">Figure 2
<p>Workflow of the proposed approach.</p>
Full article ">Figure 3
<p>NBR pre-fire, NBR post-fire, and dNBR maps for Brienza fire.</p>
Full article ">Figure 4
<p>NBR pre-fire, NBR post-fire, and dNBR maps for San Fili-Rende fire.</p>
Full article ">Figure 5
<p>NBR pre-fire, NBR post-fire. and dNBR maps for Tanca-Altara fire.</p>
Full article ">Figure 6
<p>Area of Interest as reported by CEMS (inside the white line) and Region of interest used in the present study (red areas with black squares). On the right are the same views but at a higher zoom level.</p>
Full article ">Figure 7
<p>Comparison of burned areas (red) from CEMS (on the left) with dNBRGi, dGiNBR, and dNBR for the three study cases.</p>
Full article ">Figure 8
<p>Comparison of the dNBR indices dNBRGi and dGiNBR with the reference burned area as reported by CEMS. Highlighting of false positives and negatives.</p>
Full article ">Figure 9
<p>The red colour highlights the burned areas as reported by CEMS (San Fili) and as calculated in the study. The figure shows the improvement obtained by applying one of the indices compared to dNBR.</p>
Full article ">Figure 10
<p>An example of a persistent false-positive has been highlighted within the white circle corresponding to an area of dry vegetation (San Fili). This area was not classified as burned by CEMS (<b>left</b>) but was present in the dNBR (<b>centre</b>) and indices (<b>right</b>) developed in this study.</p>
Full article ">
23 pages, 3670 KiB  
Article
Modelling Soil Moisture Content with Hydrus 2D in a Continental Climate for Effective Maize Irrigation Planning
by Nxumalo Gift Siphiwe, Tamás Magyar, János Tamás and Attila Nagy
Agriculture 2024, 14(8), 1340; https://doi.org/10.3390/agriculture14081340 - 10 Aug 2024
Viewed by 472
Abstract
In light of climate change and limited water resources, optimizing water usage in agriculture is crucial. This study models water productivity to help regional planners address these challenges. We integrate CROPWAT-based reference evapotranspiration (ETo) with Sentinel 2 data to calculate daily [...] Read more.
In light of climate change and limited water resources, optimizing water usage in agriculture is crucial. This study models water productivity to help regional planners address these challenges. We integrate CROPWAT-based reference evapotranspiration (ETo) with Sentinel 2 data to calculate daily evapotranspiration and water needs for maize using soil and climate data from 2021 to 2023. The HYDRUS model predicted volumetric soil moisture content, validated against observed data. A 2D hydrodynamic model within HYDRUS simulated temporal and spatial variations in soil water distribution for maize at a non-irrigated site in Hungary. The model used soil physical properties and crop evapotranspiration rates as inputs, covering crop development stages from planting to harvest. The model showed good performance, with R² values of 0.65 (10 cm) and 0.81 (60 cm) in 2021, 0.51 (10 cm) and 0.50 (60 cm) in 2022, and 0.38 (10 cm) and 0.72 (60 cm) in 2023. RMSE and NRMSE values indicated reliability. The model revealed water deficits and proposed optimal irrigation schedules to maintain soil moisture between 32.2 and 17.51 V/V%. This integrated approach offers a reliable tool for monitoring soil moisture and developing efficient irrigation systems, aiding maize production’s adaptation to climate change. Full article
(This article belongs to the Section Agricultural Water Management)
Show Figures

Figure 1

Figure 1
<p>Study site in Kondoros, Hungary.</p>
Full article ">Figure 2
<p>Maize at the end of June in 2021, 2022, and 2023.</p>
Full article ">Figure 3
<p>Overview of the study Hydrus-2D model soil water content simulation process (adopted from [<a href="#B65-agriculture-14-01340" class="html-bibr">65</a>]).</p>
Full article ">Figure 4
<p>Meteorological conditions at the study site.</p>
Full article ">Figure 5
<p>Maize NDVI and <span class="html-italic">ET<sub>c</sub></span> as a function of elapsed time after sowing from 2021 to 2023 with (<b>a</b>), (<b>b</b>), and (<b>c</b>), respectively.</p>
Full article ">Figure 6
<p>Soil water content at 10 cm and 60 cm depths was both measured and simulated during the 2021–2023 growing season.</p>
Full article ">
26 pages, 14290 KiB  
Article
Exploratory Analysis Using Deep Learning for Water-Body Segmentation of Peru’s High-Mountain Remote Sensing Images
by William Isaac Perez-Torres, Diego Armando Uman-Flores, Andres Benjamin Quispe-Quispe, Facundo Palomino-Quispe, Emili Bezerra, Quefren Leher, Thuanne Paixão and Ana Beatriz Alvarez
Sensors 2024, 24(16), 5177; https://doi.org/10.3390/s24165177 - 10 Aug 2024
Viewed by 488
Abstract
High-mountain water bodies represent critical components of their ecosystems, serving as vital freshwater reservoirs, environmental regulators, and sentinels of climate change. To understand the environmental dynamics of these regions, comprehensive analyses of lakes across spatial and temporal scales are necessary. While remote sensing [...] Read more.
High-mountain water bodies represent critical components of their ecosystems, serving as vital freshwater reservoirs, environmental regulators, and sentinels of climate change. To understand the environmental dynamics of these regions, comprehensive analyses of lakes across spatial and temporal scales are necessary. While remote sensing offers a powerful tool for lake monitoring, applications in high-mountain terrain present unique challenges. The Ancash and Cuzco regions of the Peruvian Andes exemplify these challenges. These regions harbor numerous high-mountain lakes, which are crucial for fresh water supply and environmental regulation. This paper presents an exploratory examination of remote sensing techniques for lake monitoring in the Ancash and Cuzco regions of the Peruvian Andes. The study compares three deep learning models for lake segmentation: the well-established DeepWaterMapV2 and WatNet models and the adapted WaterSegDiff model, which is based on a combination of diffusion and transformation mechanisms specifically conditioned for lake segmentation. In addition, the Normalized Difference Water Index (NDWI) with Otsu thresholding is used for comparison purposes. To capture lakes across these regions, a new dataset was created with Landsat-8 multispectral imagery (bands 2–7) from 2013 to 2023. Quantitative and qualitative analyses were performed using metrics such as Mean Intersection over Union (MIoU), Pixel Accuracy (PA), and F1 Score. The results achieved indicate equivalent performance of DeepWaterMapV2 and WatNet encoder–decoder architectures, achieving adequate lake segmentation despite the challenging geographical and atmospheric conditions inherent in high-mountain environments. In the qualitative analysis, the behavior of the WaterSegDiff model was considered promising for the proposed application. Considering that WatNet is less computationally complex, with 3.4 million parameters, this architecture becomes the most pertinent to implement. Additionally, a detailed temporal analysis of Lake Singrenacocha in the Vilcanota Mountains was conducted, pointing out the more significant behavior of the WatNet model. Full article
(This article belongs to the Special Issue Application of Satellite Remote Sensing in Geospatial Monitoring)
Show Figures

Figure 1

Figure 1
<p>Location of the study area.</p>
Full article ">Figure 2
<p>Landsat-8 scenes selected for study.</p>
Full article ">Figure 3
<p>Combining process from B2 to B7 into a single 6-channel image.</p>
Full article ">Figure 4
<p>From left to right: <math display="inline"><semantics> <mrow> <mo>(</mo> <mi>θ</mi> <mo>,</mo> <mi>ρ</mi> <mo>)</mo> </mrow> </semantics></math> parameter space, deskwed image, cropped image, and division of the image into 256 × 256 pixel patches.</p>
Full article ">Figure 5
<p>Mask creation process.</p>
Full article ">Figure 6
<p>WatNet model architecture.</p>
Full article ">Figure 7
<p>DeepWaterMapV2 model architecture based on 3 primary blocks.</p>
Full article ">Figure 8
<p>General architecture of WaterSegDiff based on a conditioning model and a diffusion model that integrate their information through two conditioning mechanisms, <math display="inline"><semantics> <mi mathvariant="script">U</mi> </semantics></math>-SA and SS-Former.</p>
Full article ">Figure 9
<p>SS-Former internal architecture consisting of two symmetrical cross-attention modules.</p>
Full article ">Figure 10
<p>Qualitative analysis of 5 selected samples that represent large lakes with compact structures. Showing the RGB image, ground truth, NDWI, WatNet, DeepWaterMapV2, and WaterSegDiff results. (<b>a</b>) Large and irregular lake, (<b>b</b>) two lakes with compact structure, (<b>c</b>) scene with river crossing, (<b>d</b>) large lake in mountainous region, (<b>e</b>) lake surrounded by dense vegetation.</p>
Full article ">Figure 11
<p>Qualitative analysis of 5 selected samples that represent small and dispersed lakes. Showing the RGB image, ground truth, NDWI, WatNet, DeepWaterMapV2, and WaterSegDiff results. (<b>a</b>,<b>b</b>) Snowy scene with shadows with presence of clear and turbid lakes, (<b>c</b>) completely snowy scene, (<b>d</b>,<b>e</b>) partially snowy area with scattered lakes.</p>
Full article ">Figure 12
<p>The edges extracted from Lake Singrenacocha based on NDWI, WatNet, DeepWaterMapV2, and WaterSegDiff. Highlights in yellow, green, blue, and red for the years 2014, 2016, 2018, and 2020, respectively.</p>
Full article ">Figure 13
<p>Graphical representation of the segmentation performance of Lake Singrenacocha during the years 2014, 2016, 2018, and 2020.</p>
Full article ">
18 pages, 4684 KiB  
Article
Monitoring Water Quality Parameters Using Sentinel-2 Data: A Case Study in the Weihe River Basin (China)
by Tieming Liu, Zhao Guo, Xiaoping Li, Teng Xiao, Jiaxin Liu and Yuanzhi Zhang
Sustainability 2024, 16(16), 6881; https://doi.org/10.3390/su16166881 (registering DOI) - 10 Aug 2024
Viewed by 514
Abstract
Based on Sentinel-2 multispectral image data and existing research results, the comprehensive water quality index (CWQI), NH4+-N, and total phosphorus (TP) in the Weihe River and its tributaries were estimated. Furthermore, a verified model was obtained by fitting the regression [...] Read more.
Based on Sentinel-2 multispectral image data and existing research results, the comprehensive water quality index (CWQI), NH4+-N, and total phosphorus (TP) in the Weihe River and its tributaries were estimated. Furthermore, a verified model was obtained by fitting the regression using the measured and inverted data. The verified model results show that the average relative error of the CWQI is only 9.80%, the goodness of fit of NH4+-N and TP concentrations is 0.62 and 0.61, respectively, and the average relative errors are 19.40% and 24.70%, respectively. The accuracy of the verified model is relatively high, and it can approximately invert the distribution of the three parameters of the Weihe River and its tributaries. In December 2023, except for the Bahe River between Puhua Town and Sanli Town in Lantian County, most of the water bodies in the Weihe River and its tributaries had good water quality. The study can provide an example of how to monitor water quality information using Sentinel-2 data in similar river basins. Full article
(This article belongs to the Section Sustainable Water Management)
Show Figures

Figure 1

Figure 1
<p>The distribution map of Xi’an City, Xianyang City, and the Weihe River and its tributaries.</p>
Full article ">Figure 2
<p>Water body distribution of the Weihe River and its tributaries in December 2023.</p>
Full article ">Figure 3
<p>Locations of measured sections (S1–S12), towns, and sewage treatment plants (red points indicate measured sections, green points indicate small towns, and black points indicate sewage treatment plants).</p>
Full article ">Figure 4
<p>Line chart of measured and inverted CWQI values.</p>
Full article ">Figure 5
<p>Line chart of measured and inverted NH<sub>4</sub><sup>+</sup>-N concentrations.</p>
Full article ">Figure 6
<p>Line chart of measured and inverted TP concentrations.</p>
Full article ">Figure 7
<p>Line chart of measured and verified CWQI for sections S8–S12.</p>
Full article ">Figure 8
<p>NH<sub>4</sub><sup>+</sup>-N regression results.</p>
Full article ">Figure 9
<p>TP regression results.</p>
Full article ">Figure 10
<p>Distribution of CWQI for the Weihe River and its tributaries in December 2023 estimated by the inversion model.</p>
Full article ">Figure 11
<p>Distribution of CWQI for the Weihe River and its tributaries in December 2023 obtained from the verified model (“a” represents the overall distribution of the river, and “b” represents the river segment in the area of Lantian County).</p>
Full article ">Figure 12
<p>Distribution of NH<sub>4</sub><sup>+</sup>-N for the Weihe River and its tributaries in December 2023 obtained from the inversion model.</p>
Full article ">Figure 13
<p>Distribution of NH<sub>4</sub><sup>+</sup>-N for the Weihe River and its tributaries in December 2023 obtained from the verified model (“c” represents the overall distribution of the river, and “d” represents the river segment in the area of Lantian County).</p>
Full article ">Figure 14
<p>Distribution of TP for the Weihe River and its tributaries in December 2023 obtained from the inversion model.</p>
Full article ">Figure 15
<p>Distribution of TP for the Weihe River and its tributaries in December 2023 obtained from the verified model (“e” represents the overall distribution of the river, and “f” represents the river segment in the area of Lantian County).</p>
Full article ">
20 pages, 2575 KiB  
Article
Recent Developments in Satellite Remote Sensing for Air Pollution Surveillance in Support of Sustainable Development Goals
by Dimitris Stratoulias, Narissara Nuthammachot, Racha Dejchanchaiwong, Perapong Tekasakul and Gregory R. Carmichael
Remote Sens. 2024, 16(16), 2932; https://doi.org/10.3390/rs16162932 - 9 Aug 2024
Viewed by 246
Abstract
Air pollution is an integral part of climatic, environmental, and socioeconomic current affairs and a cross-cutting component of certain United Nations Sustainable Development Goals (SDGs). Hence, reliable information on air pollution and human exposure is a crucial element in policy recommendations and decisions. [...] Read more.
Air pollution is an integral part of climatic, environmental, and socioeconomic current affairs and a cross-cutting component of certain United Nations Sustainable Development Goals (SDGs). Hence, reliable information on air pollution and human exposure is a crucial element in policy recommendations and decisions. At the same time, Earth Observation is steadily gaining confidence as a data input in the calculation of various SDG indicators. The current paper focuses on the usability of modern satellite remote sensing in the context of SDGs relevant to air quality. We introduce the socioeconomic importance of air quality and discuss the current uptake of geospatial information. The latest developments in Earth Observation provide measurements of finer spatial, temporal, and radiometric resolution products with increased global coverage, long-term continuation, and coherence in measurements. Leveraging on the two latest operational satellite technologies available, namely the Sentinel-5P and the Geostationary Environment Monitoring Spectrometer (GEMS) missions, we demonstrate two potential operational applications for quantifying air pollution at city and regional scales. Based on the two examples and by discussing the near-future anticipated geospatial capabilities, we showcase and advocate that the potential of satellite remote sensing as a, complementary to ground station networks, source of air pollution information is gaining confidence. As such, it can be an invaluable tool for quantifying global air pollution and deriving robust population exposure estimates. Full article
Show Figures

Figure 1

Figure 1
<p>Settlements with available data on PM<sub>2.5</sub> concentrations between 2010 and 2019. Adapted from the World Health Organization [<a href="#B31-remotesensing-16-02932" class="html-bibr">31</a>] with permission.</p>
Full article ">Figure 2
<p>Global gridded map of the adjusted population count. Adapted from the Gridded Population of the World (GPWv4) dataset. Source: Center for International Earth Science Information Network—CIESIN—Columbia University [<a href="#B33-remotesensing-16-02932" class="html-bibr">33</a>].</p>
Full article ">Figure 3
<p>Global annual mean of geophysical PM<sub>2.5</sub> estimates for the year 2015 based on advances in satellite observations. Black dots represent ground stations. Adapted from Hammer et al. [<a href="#B23-remotesensing-16-02932" class="html-bibr">23</a>]. Source: <a href="https://pubs.acs.org/doi/10.1021/acs.est.0c01764" target="_blank">https://pubs.acs.org/doi/10.1021/acs.est.0c01764</a> (accessed on 30 July 2024). Further permissions related to the material excerpted should be directed to the ACS.</p>
Full article ">Figure 4
<p>Annual mean of the tropospheric vertical column of NO<sub>2</sub> for the year 2021 retrieved from Sentinel-5P satellite over Bangkok, Thailand. The blue dots represent the locations of the regulatory-grade ground stations available in this region. The layers are superimposed over a natural-color satellite image of the city of Bangkok.</p>
Full article ">Figure 5
<p>An operational product from the GEMS instrument: estimated surface PM<sub>2.5</sub> concentrations over Asia acquired on 25 February 2022 (retrieved from the NIER) (<b>left</b>) and monthly mean GEMS AOD (550 nm) image for March 2023 (<b>right</b>).</p>
Full article ">Figure 6
<p>Reprocessed monthly mean of NO<sub>2</sub> (<b>left</b>) and monthly maximum for SO<sub>2</sub> (<b>right</b>) from the operational data provided by GEMS for the month of November 2023.</p>
Full article ">
15 pages, 613 KiB  
Article
A Technological Perspective of Bringing Climate Change Adaptation, Disaster Risk Reduction, and Food Security Together in South Africa
by Annegrace Zembe, Livhuwani David Nemakonde, Paul Chipangura, Christo Coetzee and Fortune Mangara
Sustainability 2024, 16(16), 6844; https://doi.org/10.3390/su16166844 - 9 Aug 2024
Viewed by 784
Abstract
As disasters and climate change risks, particularly droughts and floods, continue to affect food security globally, most governments, including South Africa, have resorted to the use of technology to incorporate climate change adaptation and disaster risk reduction to address FS issues. This is [...] Read more.
As disasters and climate change risks, particularly droughts and floods, continue to affect food security globally, most governments, including South Africa, have resorted to the use of technology to incorporate climate change adaptation and disaster risk reduction to address FS issues. This is because most institutions and policies that address climate change adaptation, disaster risk reduction, and food security operate in parallel, which usually leads to the polarisation of interventions and conflicting objectives, thus leaving the issue of FS unresolved. The study aimed to investigate how food security projects are incorporating climate change adaptation and disaster risk reduction using technology. A qualitative research design was applied, whereby in-depth interviews were conducted with ten project participants from two projects, while 24 key informants were purposively selected from government and research institutions. The study’s main findings revealed that both projects incorporate climate change adaptation and disaster risk reduction measures in most of their food value chains. Although the projects are different, they still face similar challenges, such as a lack of expertise, resources, and funding, and an inadequate regulatory environment to improve their farming practices. The study brings in the practical side of addressing the coherence between food security, climate change adaptation, and disaster risk reduction through technology. Full article
(This article belongs to the Section Sustainable Agriculture)
Show Figures

Figure 1

Figure 1
<p>Rooftop farm situated on Chamber of Mines building in Johannesburg [<a href="#B64-sustainability-16-06844" class="html-bibr">64</a>].</p>
Full article ">
18 pages, 3947 KiB  
Article
Potential of the Bi-Static SAR Satellite Companion Mission Harmony for Land-Ice Observations
by Andreas Kääb, Jérémie Mouginot, Pau Prats-Iraola, Eric Rignot, Bernhard Rabus, Andreas Benedikter, Helmut Rott, Thomas Nagler, Björn Rommen and Paco Lopez-Dekker
Remote Sens. 2024, 16(16), 2918; https://doi.org/10.3390/rs16162918 - 9 Aug 2024
Viewed by 313
Abstract
The EarthExplorer 10 mission Harmony by the European Space Agency ESA, scheduled for launch around 2029–2030, consists of two passive C-band synthetic-aperture-radar companion satellites flying in a flexible constellation with one Sentinel-1 radar satellite as an illuminator. Sentinel-1 will serve as transmitter and [...] Read more.
The EarthExplorer 10 mission Harmony by the European Space Agency ESA, scheduled for launch around 2029–2030, consists of two passive C-band synthetic-aperture-radar companion satellites flying in a flexible constellation with one Sentinel-1 radar satellite as an illuminator. Sentinel-1 will serve as transmitter and receiver of radar waves, and the two Harmonys will serve as bistatic receivers without the ability to transmit. During the first and last year of the 5-year mission, the two Harmony satellites will fly in a cross-track interferometric constellation, such as that known from TanDEM-X, about 350 km ahead or behind the assigned Sentinel-1. This constellation will provide 12-day repeat DEMs, among other regions, over most land-ice and permafrost areas. These repeat DEMs will be complemented by synchronous lateral terrain displacements from the well-established offset tracking method. In between the cross-track interferometry phases, one of the Harmony satellites will be moved to the opposite side of the Sentinel-1 to form a symmetric bistatic “stereo” constellation with ±~350 km along-track baseline. In this phase, the mission will provide opportunity for radar interferometry along three lines of sight, or up to six when combining ascending and descending acquisitions, enabling the measurement of three-dimensional surface motion, for instance sub- and emergence components of ice flow, or three-dimensional deformation of permafrost surfaces or slow landslides. Such measurements would, for the first time, be available for large areas and are anticipated to provide a number of novel insights into the dynamics and mass balance of a range of mass movement processes. Full article
(This article belongs to the Special Issue Remote Sensing of the Cryosphere II)
Show Figures

Figure 1

Figure 1
<p>Graphical summary of Harmony land-ice measurement objectives.</p>
Full article ">Figure 2
<p>The Harmony mission consists of two satellites with one passive SAR instrument each, Harmony-A and Harmony-B. They will fly in two alternating configurations in convoy with Sentinel-1, which serves as a radar transmitter for the two receive-only Harmony satellites. (<b>a</b>) The stereo configuration is optimized to measure surface motion vectors on land and ocean and is foreseen for years 2–4 of the mission. (<b>b</b>) The cross-track interferometric (XTI) configuration is optimized to measure land surface topography every 12 days during years 1 and 5.</p>
Full article ">Figure 3
<p>Schematics of Harmony time series for the cases of (<b>a</b>) small elevation changes (few meters over several months to years) and of (<b>b</b>) large elevation changes (tens of meters) for a hypothetical glacier point. (<b>b</b>) shows combined elevation changes and horizontal speeds. During the full mission years 1 and 5, Harmony is foreseen to measure dozens of DEMs over glaciers and ice-sheet margins globally (black points with error bars). The blue and brown curves indicate hypothetical idealized glacier elevation variations over time. Note that the vertical axes of panel (<b>a</b>,<b>b</b>) have scales that are different by an order of magnitude. (<b>a</b>) The mission will be able to deliver glacier volume changes Δh from differencing the DEM stacks from years 1 and 5. For small glacier thickness changes, the penetration bias between real surface and a radar-interferometric DEM is substantial, relative to the expected elevation changes, and needs to be dealt with. (<b>b</b>) Harmony’s repeat DEMs, potentially combined with lateral surface displacements from radar offset tracking, can be used to study the interplay between changes in ice dynamics and ice thickness, such as for calving glaciers or glacier surges. Compared to the large elevation changes expected for such cases, radar penetration and DEM errors are relatively small and of less concern than in case (<b>a</b>).</p>
Full article ">Figure 4
<p>Visualization and comparison of some expected Harmony results using Tunabreen (78.46°N, 17.4°E), Svalbard, as an example. (<b>a</b>) Sentinel-1 amplitude image (gamma0, 12 March 2020). (<b>b</b>) Sentinel-1 12-day offsets over 19–31 January 2019. (<b>c</b>) Sentinel-2 image (infrared false color, for orientation only, 2 August 2019). (<b>d</b>) Elevation differences from Arctic-DEM [<a href="#B29-remotesensing-16-02918" class="html-bibr">29</a>] strips of 17 March 2015 and a mosaic of 11 and 15 March 2020 resampled to 100 m resolution with random noise of ±0.5 m/yr added. (<b>e</b>) Same as (<b>d</b>), but resampled to 50 m and with ±0.2 m/yr noise added. (<b>f</b>) Elevation trends from ASTER stacks 2015–2019 [<a href="#B24-remotesensing-16-02918" class="html-bibr">24</a>]. (<b>g</b>) TanDEM-X topographic change product, computed between TanDEM-X elevations compiled over 2010–2014 and elevations from 2017. Note, panels (<b>d</b>,<b>e</b>) are visualizations (not simulations) of Harmony XTI 5-year DEM differences for the threshold (panel (<b>d</b>)) and goal requirements (panel (<b>e</b>)). They do not include potential errors from SAR XTI processing (e.g., steep terrain, penetration), but rather errors from optical stereo processing (e.g., small clouds, lack of visual contrast). Black glacier outlines are from Randolph glacier inventory v7. Note also that the ASTER DEM differences are just shown for visual comparison, and that the potential reasons for the differences between the DEM differences from ASTER and Arctic-DEM are not the subject of the present contribution.</p>
Full article ">Figure 5
<p>Harmony’s interferometric line-of-sight diversity. (<b>a</b>) Sentinel-1 alone is able to provide radar interferometry along one line-of-sight (two if both ascending and descending orbits are combined; grey column to the right). (<b>b</b>) In XTI configuration, Harmony will provide two lines of sight (strictly speaking three, but the lines of sight of both Harmony satellites are very similar), and four if both ascending and descending orbits are combined. (<b>c</b>) In stereo configuration, Harmony will provide interferometry along three, or six, lines of sight, respectively. The three lines of sight per orbit in stereo configuration lie, though, approximately in one oblique plane. (<b>d</b>) Harmony interferograms from one orbit can also be combined with interferograms from an opposite orbit of the other Sentinel-1, providing four lines of sight.</p>
Full article ">Figure 6
<p>Interferometric measurements by the Harmony mission along two or more lines of sight will facilitate the measurement of three-dimensional glacier flow close to the surface. Such measurements will connect between the vertical component of ice flow, thickness changes over time, and local mass balance (i.e., directly measure components of the so-called kinematic boundary condition). Changes in the penetrated snow and firn pack could lead, though, to offsets between the true ice particle displacement and the displacement between the phase centers of the radar waves.</p>
Full article ">Figure 7
<p>2D schemes of different idealized cross-sections and kinematics of slow landslides. (<b>a</b>) surface slope steeper than bedding slope, (<b>b</b>) both similar, (<b>c</b>) bedding slope steeper than surface, (<b>d</b>) rotational landslide. Interferometric measurements from repeat Harmony data would enable estimating two- or three-dimensional surface velocities (red arrows). DEM stacks from Harmony’s XTI phases (or other DEMs) can give elevation changes (blue arrows). Surface at time 1, black line; surface at time 2, black dashed line; bedding plane, grey line.</p>
Full article ">Figure 8
<p>(<b>a</b>) Normalized differences between the two unit-less amplitude images of a TanDEM-X bi-static acquisition of 13 April 2015 over Aletsch Glacier (<b>b</b>) and Bernese Alps, Switzerland, with approx. 500 m along-track baseline and 2 km cross-track baseline, i.e., approx. 2060 m total baseline. The noise-filtered color-coded normalized differences are transparently laid over one of the amplitude images. The more blue or red, resp., the stronger the differences are. Image in raw radar geometry, flying direction from top to bottom (descending), look direction from right to left. Strongest variations in differences between the two amplitude images are found for higher elevations, likely related to variations in snow cover type. (For other examples, see Stefko et al. [<a href="#B81-remotesensing-16-02918" class="html-bibr">81</a>].)</p>
Full article ">
19 pages, 7874 KiB  
Article
Mapping the Fraction of Vegetation Coverage of Potamogeton crispus L. in a Shallow Lake of Northern China Based on UAV and Satellite Data
by Junjie Chen, Quanzhou Yu, Fenghua Zhao, Huaizhen Zhang, Tianquan Liang, Hao Li, Zhentan Yu, Hongli Zhang, Ruyun Liu, Anran Xu and Shaoqiang Wang
Remote Sens. 2024, 16(16), 2917; https://doi.org/10.3390/rs16162917 - 9 Aug 2024
Viewed by 437
Abstract
Under the background of global change, the lake water environment is facing a huge threat from eutrophication. The rapid increase in curly-leaf pondweed (Potamogeton crispus L.) in recent years has seriously threatened the ecological balance and the water diversion safety of the [...] Read more.
Under the background of global change, the lake water environment is facing a huge threat from eutrophication. The rapid increase in curly-leaf pondweed (Potamogeton crispus L.) in recent years has seriously threatened the ecological balance and the water diversion safety of the eastern route of China’s South-to-North Water Diversion Project. The monitoring and control of curly-leaf pondweed is imperative in shallow lakes of northern China. Unmanned Aerial Vehicles (UAVs) have great potential for monitoring aquatic vegetation. However, merely using satellite remote sensing to detect submerged vegetation is not sufficient, and the monitoring of UAVs on aquatic vegetation is rarely systematically evaluated. In this study, taking Nansi Lake as a case, we employed Red–Green–Blue (RGB) UAV and satellite datasets to evaluate the monitoring of RGB Vegetation Indices (VIs) in pondweed and mapped the dynamic patterns of the pondweed Fractional Vegetation Coverage (FVC) in Nansi Lake. The pondweed FVC values were extracted using the RGB VIs and the machine learning method. The extraction of the UAV RGB images was evaluated by correlations, accuracy assessments and separability. The correlation between VIs and FVC was used to invert the pondweed FVC in Nansi Lake. The RGB VIs were also calculated using Gaofen-2 (GF-2) and were compared with UAV and Sentinel-2 data. Our results showed the following: (1) The RGB UAV could effectively monitor the FVC of pondweed, especially when using Support Vector Machine that (SVM) has a high ability to recognize pondweed in UAV RGB images. Two RGB VIs, RCC and RGRI, appeared best suited for monitoring aquatic plants. The correlations between four RGB VIs based on GF-2, i.e., GCC, BRI, VDVI, and RGBVI and FVCSVM calculated by the UAV (p < 0.01) were better than those obtained with other RGB VIs. Thus, the RGB VIs of GF-2 were not as effective as those of the UAV in pondweed monitoring. (2) The binomial estimation model constructed by the Normalized Difference Water Index (NDWI) of Sentinel-2 showed a high accuracy (R2 = 0.7505, RMSE = 0.169) for pondweed FVC and can be used for mapping the FVC of pondweed in Nansi Lake. (3) Combined with the Sentinel-2 time-series data, we mapped the dynamic patterns of pondweed FVC in Nansi Lake. It was determined that the flooding of pondweed in Nansi Lake has been alleviated in recent years, but the rapid increase in pondweed in part of Nansi Lake remains a challenging management issue. This study provides practical tools and methodology for the innovative remote sensing monitoring of submerged vegetation. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the study area with the distribution of the sample sites. (<b>a</b>) Location of the study area; (<b>b</b>) detailed location of the study area; (<b>c</b>) sample points in the study area; (<b>d</b>) sample points in Dushan Lake; (<b>e</b>) sample points in Weishan Lake.</p>
Full article ">Figure 2
<p>A flowchart of this study.</p>
Full article ">Figure 3
<p>Illustration of RGB VI separability determination.</p>
Full article ">Figure 4
<p>FVC box plots extracted by different methods, including SVM, the Dimidiate Pixel Model and dynamic thresholding.</p>
Full article ">Figure 5
<p>The correlation coefficients between the different FVC values extracted by SVM, the Dimidiate Pixel Model and dynamic thresholding.</p>
Full article ">Figure 6
<p>(<b>a</b>) Correlation between remote sensing VIs, FVC<sub>NDVI</sub>, FVCsvm and the mean RGB VIs, (<b>b</b>) correlation between remote sensing VIs, FVC<sub>NDVI</sub>, FVCsvm and FVC values by UAV. “*” and “**” represent significant differences, with <span class="html-italic">p</span> &lt; 0.05 and <span class="html-italic">p</span> &lt; 0.01, respectively.</p>
Full article ">Figure 7
<p>Correlation analysis between RGB VIs by GF-2 and the means of RGB VIs by the UAV (<b>a</b>), FVC by the UAV (<b>b</b>) and remote sensing VIs by Sentinel-2 (<b>c</b>). “*” and “**” represent significant differences, with <span class="html-italic">p</span> &lt; 0.05 and <span class="html-italic">p</span> &lt; 0.01, respectively.</p>
Full article ">Figure 8
<p>Accuracy assessment results for RGB VIs. (<b>a</b>) Overall accuracy, (<b>b</b>) F1 score.</p>
Full article ">Figure 9
<p>Statistical results of separability in the acquired images for RGB VIs.</p>
Full article ">Figure 10
<p>Comparison between estimated and measured pondweed FVC.</p>
Full article ">Figure 11
<p>Mapping pondweed FVC in Nansi Lake based on the NDWI binomial estimation model (14 May 2023).</p>
Full article ">Figure 12
<p>Seasonal change in pondweed FVC in Nansi Lake, 2023.</p>
Full article ">Figure 13
<p>Inter-annual changes in pondweed FVC in Nansi Lake, 2018–2023.</p>
Full article ">Figure 14
<p>Different growth periods of pondweed imaged by the RGB UAV.</p>
Full article ">
23 pages, 11056 KiB  
Article
Co-Kriging-Guided Interpolation for Mapping Forest Aboveground Biomass by Integrating Global Ecosystem Dynamics Investigation and Sentinel-2 Data
by Yingchen Wang, Hongtao Wang, Cheng Wang, Shuting Zhang, Rongxi Wang, Shaohui Wang and Jingjing Duan
Remote Sens. 2024, 16(16), 2913; https://doi.org/10.3390/rs16162913 - 9 Aug 2024
Viewed by 414
Abstract
Mapping wall-to-wall forest aboveground biomass (AGB) at large scales is critical for understanding global climate change and the carbon cycle. In previous studies, a regression-based method was commonly used to map the spatially continuous distribution of forest AGB with the aid of optical [...] Read more.
Mapping wall-to-wall forest aboveground biomass (AGB) at large scales is critical for understanding global climate change and the carbon cycle. In previous studies, a regression-based method was commonly used to map the spatially continuous distribution of forest AGB with the aid of optical images, which may suffer from the saturation effect. The Global Ecosystem Dynamics Investigation (GEDI) can collect forest vertical structure information with high precision on a global scale. In this study, we proposed a collaborative kriging (co-kriging) interpolation-based method for mapping spatially continuous forest AGB by integrating GEDI and Sentinel-2 data. First, by fusing spectral features from Sentinel-2 images with vertical structure features from GEDI, the optimal estimation model for footprint-level AGB was determined by comparing different machine-learning algorithms. Second, footprint-level predicted AGB was used as the main variable, with rh95 and B12 as covariates, to build a co-kriging guided interpolation model. Finally, the interpolation model was employed to map wall-to-wall forest AGB. The results showed the following: (1) For footprint-level AGB, CatBoost achieved the highest accuracy by fusing features from GEDI and Sentinel-2 data (R2 = 0.87, RMSE = 49.56 Mg/ha, rRMSE = 27.06%). (2) The mapping results based on the interpolation method exhibited relatively high accuracy and mitigated the saturation effect in areas with higher forest AGB (R2 = 0.69, RMSE = 81.56 Mg/ha, rRMSE = 40.98%, bias = −3.236 Mg/ha). The mapping result demonstrates that the proposed method based on interpolation combined with multi-source data can be a promising solution for monitoring spatially continuous forest AGB. Full article
(This article belongs to the Special Issue Remote Sensing and Lidar Data for Forest Monitoring)
Show Figures

Figure 1

Figure 1
<p>Location of the study area. (<b>a</b>) The administrative boundary of the State of California and the location of the study area within the state (i.e., where the red box is); (<b>b</b>) ground elevation distribution of the study area; (<b>c</b>) local distribution of GEDI footprints.</p>
Full article ">Figure 2
<p>Technology road map of this study.</p>
Full article ">Figure 3
<p>Filtered distribution of GEDI footprints: (<b>a</b>) urban area; (<b>b</b>) vegetation area. The green area represents vegetation, the purple area represents urban areas, and the blue area represents the ocean.</p>
Full article ">Figure 4
<p>The correlation between the features and forest AGB, with all significance levels of the selected features at 0.01.</p>
Full article ">Figure 5
<p>Feature selection outcomes. (<b>a</b>) Features selected using the stepwise method. (<b>b</b>) Features selected using the random forest method.</p>
Full article ">Figure 6
<p>Scatter plots between reference AGB and predicted AGB derived from the combined optical and GEDI data. The red dotted lines represent the fitted trend-lines.</p>
Full article ">Figure 7
<p>Histograms of AGB: (<b>a</b>) ALS-derived and (<b>b</b>) interpolation-method-derived.</p>
Full article ">Figure 8
<p>Boxplot of AGB distribution. IQR = Q3 − Q1, where Q1 is the first quartile and Q3 is the third quartile.</p>
Full article ">Figure 9
<p>Forest AGB mapping results. (<b>a</b>) Satellite images of the study area. (<b>b</b>) Forest AGB map derived from the co-kriging interpolation model. The red line represents the boundary of the study area (i.e., Sonoma County).</p>
Full article ">Figure 10
<p>Comparison of interpolation accuracy for different covariate combinations.</p>
Full article ">Figure 11
<p>Accuracy assessment of the wall-to-wall forest AGB products. (<b>a</b>) The predicted AGB derived from co-kriging interpolation. (<b>b</b>) The predicted AGB derived from the regression method. The red dotted lines represent the fitted trend-lines.</p>
Full article ">
Back to TopTop