Have a personal or library account? Click to login
Aboveground biomass growth mapping and CO2 uptake of multispecies seagrass meadow in Menjangan Besar Island, Karimunjawa Islands, using Sentinel-2 imagery Cover

Aboveground biomass growth mapping and CO2 uptake of multispecies seagrass meadow in Menjangan Besar Island, Karimunjawa Islands, using Sentinel-2 imagery

Open Access
|Dec 2025

Figures & Tables

Figure 1

Menjangan Besar Island, Karimunjawa Islands (True colour composite of Sentinel-2A imagery background acquired on 2 May 2024).
Menjangan Besar Island, Karimunjawa Islands (True colour composite of Sentinel-2A imagery background acquired on 2 May 2024).

Figure 2

Research flowchart. AGB, aboveground biomass; OCC, organic carbon content; PCA, principal component analysis; RMSE, root mean square error; SR, surface reflectance.
Research flowchart. AGB, aboveground biomass; OCC, organic carbon content; PCA, principal component analysis; RMSE, root mean square error; SR, surface reflectance.

Figure 3

Benthic habitat classification result of Menjangan Besar Island (SR image: OA: 67.9%, ntree: 300, mtry: log and impurity function: Gini index). SR, surface reflectance.
Benthic habitat classification result of Menjangan Besar Island (SR image: OA: 67.9%, ntree: 300, mtry: log and impurity function: Gini index). SR, surface reflectance.

Figure 4

Map of seagrass percent cover on Menjangan Besar Island (inversed PCA-derived SR image: ntree 200, mtry square root, nodesize 4, RMSE 12.5%, r 0.62, and R2 0.38). PCA, principal component analysis; RMSE, root mean square error; SR, surface reflectance.
Map of seagrass percent cover on Menjangan Besar Island (inversed PCA-derived SR image: ntree 200, mtry square root, nodesize 4, RMSE 12.5%, r 0.62, and R2 0.38). PCA, principal component analysis; RMSE, root mean square error; SR, surface reflectance.

Figure 5

Map of seagrass species composition of Menjangan Besar Island (PC SR image: OA 62.5%, ntree 100, mtry square root and impurity function entropy). PC, principal component; SR, surface reflectance.
Map of seagrass species composition of Menjangan Besar Island (PC SR image: OA 62.5%, ntree 100, mtry square root and impurity function entropy). PC, principal component; SR, surface reflectance.

Figure 6

The AGB growth map of multispecies seagrass meadows on Menjangan Besar Island. AGB, aboveground biomass.
The AGB growth map of multispecies seagrass meadows on Menjangan Besar Island. AGB, aboveground biomass.

Figure 7

The 1:1 plot between reference and estimated AGB growth values for classes: (A) Ea, (B) Th-Cr, (C) Ea-Th, (D) BS. AGB, aboveground biomass; BS, bare substrate.
The 1:1 plot between reference and estimated AGB growth values for classes: (A) Ea, (B) Th-Cr, (C) Ea-Th, (D) BS. AGB, aboveground biomass; BS, bare substrate.

Figure 8

The CO2 uptake estimation map of the multispecies seagrass meadow of Menjangan Besar Island.
The CO2 uptake estimation map of the multispecies seagrass meadow of Menjangan Besar Island.

Figure 9

The 1:1 plot of field CO2 uptake reference and estimated values for: (A) Ea, (B) Th-Cr, (C) Ea-Th and (D) BS. BS, bare substrate.
The 1:1 plot of field CO2 uptake reference and estimated values for: (A) Ea, (B) Th-Cr, (C) Ea-Th and (D) BS. BS, bare substrate.

Figure 10

Relationship between seagrass percent cover and AGB growth at: (A) community, (B) species Ea, and (C) species Th-Cr. AGB, aboveground biomass.
Relationship between seagrass percent cover and AGB growth at: (A) community, (B) species Ea, and (C) species Th-Cr. AGB, aboveground biomass.

Overall accuracy of benthic habitat classification experimental results using SR, PC SR, deglint and DII images

ntreemtryImpurityOverall accuracy (%)
SRPC SRDeglintDII
100Square rootGini67.6661.4762.3953.44
Entropy66.9760.5563.5354.36
LogGini67.2060.3261.2455.96
Entropy66.5162.1662.1655.96
200Square rootGini66.9761.0163.3054.36
Entropy66.2862.1662.8455.05
LogGini67.4361.7062.3953.67
Entropy66.9762.3960.7855.05
300Square rootGini66.7461.9362.8454.36
Entropy67.6661.7061.0155.28
LogGini67.8961.2461.7054.36
Entropy66.9761.0161.7055.50
400Square rootGini67.6661.4762.3954.59
Entropy67.2061.7061.7054.59
LogGini66.9761.7061.9355.05
Entropy66.7462.6161.4754.82

AGB growth calculation results for Ea and Th species

No.SpeciesPercent cover (%)Shoot density (shoot · m−2)Average growth rate (cm · leaf−1 · 5 days)DW (gDW · shoot−1 · 5 days)AGB growth (gDW · m−2 · year)−1
1.Ea90.91756.86n = 20Sum = 157.7 cm · 5 daysMin = 2.6Max = 12.80.117641.34
2.54.5552444.66
3.52.1773624.24
4.76.4784718.30
5.36.3643367.70
6.45.83105897.88
7.66.6755470.32
8.20.0034290.74
9.30.0079675.55
10.37.5075641.34
11.Th, Cr30.431675.63n = 14Sum = 78.8 cm · 5 daysMin = 3.9Max = 7.70.03543.09
12.50.0024563.21
13.86.36718185.25
14.56.5227871.72
15.31.8217043.86
16.57.1419049.02
17.59.0937496.49
18.54.55434111.97
19.39.1320953.92
20.45.4524462.95

The R2 value and resultant regression function from: seagrass at the community level, Ea, and Th-Cr AGCG from AGBG leaf-marking samples_

SpeciesR2Regression equation
Community0.88AGCGc = ((2E–05) × (AGBGc2)) + (0.073 × AGBGc) – 0.3845
Ea1AGCGEa = ((4E–19) × (AGBGEa2)) + (0.0621 × AGBGEa) + (3E–13)
Th-Cr1AGCGTh-cr = ((–2E–19) × (AGBGTh-Cr2)) + (0.1131 × AGBGTh-Cr) – (1E–13)

The R2 value and modelling regression equation of seagrass at the community level, Ea and Th-Cr AGB growth_

SpeciesR2Regression equation
Community0.42AGB growth = (–0.0291 × (pCv2)) + (13.793 × PCv) + 57.885
Ea0.58AGB growthEa = (–0.1449 × (pCv2)) + (19.152 × PCv) + 36.171
Th-Cr0.87AGB growthTh.Cr = (θ.191 × (pCv2)) + (3.3648 × PCv) + 60.809

Error matrix of benthic habitat classification using SR image, ntree 300, mtry log and impurity function Gini index_

Reference class
BSDCDMCBSMCMaMSBSC + MTotalUA (%)
Classified classBSD149121031020273.8
CD271135049574.7
MCBS1712300226347.6
MCMa00000220
MSBS2501 4507163.4
C + M010101333.39
Total19385656789436
PA (%)77.283.546.2057.711.1 OA = 67.9%

Experimental results of RF regression modelling seagrass AGB growth

ClassEaTh-CrEa-ThBS
ntreemtryNodesizeRMSE (gDW · m−2 · year−1rR2RMSE (gDW · m−2 · year−1RR2RMSE (gDW · m−2 · year−1rR2RMSE (gDW · m−2 · year−1rR2
100Square root262.580.930.86136.840.760.5897.410.760.58106.800.820.67
362.160.930.86129.220.800.6496.540.800.64105.400.820.67
462.480.930.86130.560.790.63104.520.820.66100.990.840.70
562.900.930.86127.830.800.65112.180.780.6294.930.880.78
Log262.480.930.87128.910.800.6498.820.770.59109.290.810.65
362.090.930.87134.420.770.59102.370.790.62106.940.810.66
463.020.930.86133.720.770.60104.970.820.6795.030.870.76
563.210.930.86125.850.820.6695.920.870.7699.890.850.73
200Square root261.040.930.87132.940.780.60102.630.740.55101.430.840.71
362.110.930.86132.320.780.61100.140.780.61104.070.830.68
461.770.930.87128.830.800.6494.220.830.7097.220.860.74
564.600.930.86125.930.810.66102.250.820.6797.970.870.75
Log261.340.930.87135.830.760.58100.430.760.58103.970.830.69
360.830.930.87130.380.790.62102.760.770.59101.200.840.71
461.930.930.87130.450.790.6297.180.870.7698.900.850.73
563.240.930.86125.550.820.6798.600.870.7597.310.880.77

Rules for developing benthic habitat classification scheme for mapping seagrass AGB growth mapping_

CategoryRules
DominantThe percent cover of a class is >80% or <80%, while other classes are <20%.
IrrelevantThe percent cover of a class is <20% so the class is not included in the classification.
MixThe difference in percent cover between classes is <20%.
Sub-classThe difference in percent cover between classes is >20% and the subclass name is added with a ‘+’ sign.

Classification scheme for seagrass species composition mapping (Wicaksono & Hafizt, 2013)

ClassSpecies compositionExplanation
EaEnhalus acoroidesGrows and extends vertically in the water column and sometimes reaches the surface. Other life-forms may exist in this class but are not significant.(Ea species ≥ 70%, other species ≤ 30%)
Th-CrThalassia hemprichii, Cymodocea rotundata, Halodule uninervis, Syringodium isoetifolium, Halophila ovalisGrows to cover the substrate but does not extend vertically in the water column. Other life-forms may exist in this class but are not significant.(Ea species ≤ 70%, other species ≥ 30%)
Ea-ThEnhalus acoroides, Thalassia hemprichii, Cymodocea rotundata, Halodule uninervis, Syringodium isoetifolium, Halophila ovalisA mix of Ea and Th, Cr species with a significant percent cover ratio.

Number of training area and accuracy assessment samples for each benthic habitat class_

Benthic habitat classTraining areaAccuracy assessment samples
‘Bare substrate dominated’219194
‘Coral dominated’109108
‘Mix of coral and bare substrate’6766
‘Mix of coral and macroalgae’66
‘Mix of seagrass and bare substrate’7978
‘Coral + macroalgae’1010
Total490462

The R2 value and resultant regression function from: seagrass at the community level, Ea, and Th-Cr AGCG from photo quadrat samples of AGBG_

ClassR2Regression equation
Ea1AGCGEa = ((2E–18) × (AGBGEa2)) + (0.0621 × AGBGEa) + (4E–13)
Th-Cr1AGCGTh-cr = ((–2E–18) × (AGBGTh-cr2)) + (0.1131 × AGBGTh-cr) + (2E–13)
Ea-Th1AGCGEa-Th = ((2E-05) × (AGBGEa-Th2)) + (0.073 × AGBGEa-Th) – 0.3845
BS1AGCGBS = ((2E–05) × (AGBGBS2)) + (0.073 × AGBGBS) – 0.3845
DOI: https://doi.org/10.26881/oahs-2025.1.27 | Journal eISSN: 1897-3191 | Journal ISSN: 1730-413X
Language: English
Page range: 303 - 327
Submitted on: May 17, 2025
|
Accepted on: Nov 21, 2025
|
Published on: Dec 31, 2025
Published by: University of Gdańsk
In partnership with: Paradigm Publishing Services
Publication frequency: 4 issues per year

© 2025 Aurellia Devina Putri, Pramaditya Wicaksono, published by University of Gdańsk
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.