The mineralogic and isotopic fingerprint of equatorial carbonates: Kepulauan Seribu, Indonesia

Kepulauan Seribu is an isolated patch reef complex situated in the Java Sea (Indonesia) and is a typical example for a humid, equatorial carbonate system. We investigate the mineralogical and isotopic fingerprint of Panggang, one of the reef platforms of Kepulauan Seribu, to evaluate differences to other carbonate systems, using isotope in combination with XRD and SEM analysis. A characteristic property of shallow water (< 20 m) sediments from Kepulauan Seribu is their increased LMC content (~ 10%) derived from some genera of rotaliid foraminifers and bivalves. The relative abundance of these faunal elements in shallow waters might be related to at least temporary turbid conditions caused by sediment-laden river runoff. This influence is also evidenced by the presence of low amounts of siliciclastic minerals below the regional wave base. Kepulauan Seribu carbonates are characterized by very low δ13C and δ18O values. This is related to the isotopically depleted riverine input. The δ13CDIC in riverine water is reduced by the contribution of 12C from riverside mangroves. Deep atmospheric convection and intensive rains contribute 18O-depleted freshwater in the river catchments, finally reducing salinity in the Java Sea. The depleted δ13C signature in carbonates is further enhanced by the lack of green algae and inorganic carbonates and abundance of coral debris. Low δ18O values in carbonates are favored by the high water temperatures in the equatorial setting. Since equatorial carbonates in SE Asia, including the Java Sea, are typically influenced by high turbidity and/or river runoff, the observed distinctively low isotope values likely are characteristic for equatorial carbonate systems in the region.


Introduction
During the last decade the awareness has grown that carbonate systems in the humid equatorial zone have common characteristics that distinguishes them from carbonates in more arid regions (Tomascik et al. 1997;Wilson 2008). Features of these equatorial carbonates include the common occurrence of photozoan assemblages adapted to oligophotic (low light) conditions, common elements of the Heterozoan Association and mixed carbonate/clastic deposits. Humid equatorial carbonates are also characterized by a lack of coated grains, aggregates, chemically precipitated mud and the absence of associated evaporites (Wilson 2002). They typically occur in a zone between 15/20° N and 15/20° S in areas with warm (> 22 °C) surface waters and heavy precipitation (> 1500 mm per annum; Wilson, 2012). Due to the increased precipitation, carbonate deposition takes place in shelfal waters of normal to reduced salinity and high clastic influx. River runoff and upwelling cause increased nutrient concentrations in surface waters, which in turn lead to a shallow position of the photic zone.
A key region for the study of these humid equatorial carbonates is SE Asia, where the shallow shelfs are strongly affected by river runoff and/or upwelling (Tomascik et al. 1997;Wilson 2002). The diversity of environments and habitats renders SE Asia an ideal natural laboratory to investigate the variability of equatorial carbonate processes 1 3 and deposits (Tomascik et al. 1997;Wilson 2002Wilson , 2008. Sediment discharge from each of the four largest SE Asian islands varies from ca 300 to 1650 million tons/year (Milliman et al. 1999), leading to high turbidity and nutrient concentration in surface waters. Depth of abundant coral growth in SE Asian reefs is therefore limited to approximately 20-30 m and is considered among the shallowest in the world (Schlager 1992;Wilson 2008). This shallow photic zone corresponds with low light penetration in an oligo-mesotropic region characterized by nutrient input from river runoff, upwelling and plankton blooms (Wilson 2008). Due to these environmental conditions, SE Asian carbonates can be considered as typical equatorial carbonate systems.
The low salinity in the semi-enclosed seas of SE Asia also leads to a relatively low aragonite supersaturation compared to more arid regions, such as the tropical, western Atlantic (TWA). The development of Cenozoic carbonate systems in SE Asia is thought to be strongly influenced by regional changes in nutrient input and precipitation (Wilson 2008). Late Oligocene carbonates from SE Asia were often characterised by calcitic dominated assemblages. This was likely promoted by the decreased aragonite supersaturation of the low salinity waters in the semi-enclosed seas of SE Asia (Wilson 2012). An increase of river runoff and associated decrease in aragonite supersaturation is also held responsible for a shelf wide switch from inorganic aragonite to bioclastic calcite dominated sedimentation on the tropical NW Shelf of Australia during the early Holocene (Hallenberger et al. 2019;Ishiwa et al. 2019); a process which seems to have occurred repeatedly during the Pleistocene (Gallagher et al. 2018). Studies of modern carbonate systems from northern Australia and SE Asia suggest that faunal assemblages from the oligophotic zone are relatively calcite-rich (Wilson and Vecsei 2005;Wilson 2012). Generally, this change from coral dominated euphotic to foramol dominated oligophotic zones takes place in a water depth of around 20 m, but locally can be as shallow as 6 m in turbid coastal water (Wilson and Vecsei 2005).
Equatorial carbonates share some similarities with cool water systems including the abundance of calcitic components and heterozoan elements. In some cases, this can complicate the distinction between non-tropical and equatorial carbonates (Wilson and Vecsei 2005). Wilson (2012), therefore, recommended the use of additional geochemical and isotopic proxies to identify warm temperatures and low salinities associated with regional equatorial carbonates in the rock record. However, stable isotope studies from modern SE Asia that could be used as a template for comparison with fossil counterparts are still missing. It is also unclear if the predominance of calcitic bioclasts, that is apparent in the oligophotic zone of many carbonate systems in SE Asia, is visible in the mineralogy of the coral dominated shallow water area. In this study, we investigate how far the mineralogy and stable isotope signature of shallow water (< 20 m) equatorial carbonate from Kepulauan Seribu in the Java Sea is distinct from other carbonate systems. Stable isotope (δ 13 C and δ 18 O) analysis in combination with X-ray diffraction (XRD) and Scanning Electron Microscopy (SEM) is used to provide a fingerprint for equatorial carbonates from SE Asia. This data might be used as a template for comparison with carbonates in the rock record.

Study area
Situated in the Java Sea, Kepulauan Seribu (The Thousand Islands) consists of numerous coral reef islands ranging in size from a few meters across to length greater than one kilometer. From about 10 km off the northwest Java coastline, the island chain extends for about 50 km into the Java Sea. Kepulauan Seribu is situated on the NNE-SSW structural high of the Seribu platform, which is separated from the Sunda Basin towards the west and Arjuna basin to the east by the active N-S oriented Seribu Fault system (Scrutton 1976a). The north-south elongation of the platform is largely controlled by the regional fault pattern that was established during the early Cenozoic and reactivated at various times since (Park et al. 2010). The patch reef system of Kepulauan Seribu is located on the shallow Sunda Shelf where it rises from a water depth of 30-50 m to the sea surface. A deep and several kilometers wide, east-west oriented channel separates Pari Island to the south from most other islands to its north (Fig. 1). It appears that this channel forms an important hydrologic barrier, channeling the sediment-laden coastal waters off Java and Kalimantan away from the reef (Tomascik et al. 1997;Jordan 1998). The reefal buildups on the platform vary in size and shape but most of them are elongated in a roughly east-west orientation. This pronounced orientation is inherited from the Pleistocene drainage pattern which was abandoned and truncated by river capture prior to the Holocene transgression (Park et al. 2010). The islands of Kepulauan Seribu are considered to be located in a relatively sheltered position protected from severe storms and ocean swells due to the fact that it is surrounded by major landmasses, such as Sumatra, Java and Kalimantan (Fig. 1).
The controlling influence on modern reef growth and morphology in Kepulauan Seribu is believed to be the seasonal change in wind and current directions (Scrutton 1976a). The Java Sea is characterized by a monsoonal climate with a seasonal reversal between the West Monsoon (December-February), and the East Monsoon (April-October). The West Monsoon winds are known to be stronger than the East Monsoon winds, but less persistent (3 vs. 7 months). The East Monsoon winds therefore are the dominant factor for the net sediment transport (Poerbandono 1 3 2016). The climate in Kepulauan Seribu is tropical with relative humidity of 80-90% (Tomascik et al. 1997). The maximum sea surface temperature (SST) is observed in April (29.8 °C) and minimum in January (28.6 °C) (Cahyarini et al. 2016). Surface salinity in the Java Sea (30-34‰) is well below normal open ocean levels, reaching as low as 22‰ in the vicinity of some of the larger river mouths. The low salinity of the Java Sea is due to heavy precipitation and intensive river runoff from Java and Kalimantan (Sofian and Kozai 2007). The water column of the shallow western Java Sea is well mixed and shows little temperature or salinity stratification (Hamzah et al. 2020).
Water along Kepulauan Seribu is highly eutrophic (Baum et al. 2015). Seasonal variability of Chlorophyll-α distribution in the Java Sea is influenced by monsoon, through water exchange with adjacent seas and nutrient supply from river discharge (Koropitan and Ikeda 2016). In the eastern Java Sea, episodic upwelling leads to phytoplankton blooms during the East Monsoon (Koropitan and Ikeda 2016). The role of riverine nutrient input during the West Monsoon (rainy season) is limited to coastal regions (Koropitan and Ikeda 2016). The reefs in Kepulauan Seribu are directly influenced by intense anthropogenic pressure from local as well as regional sources (Baum et al. 2015). Polluted water masses from the Greater Jakarta Area are considered to affect reefs in nearshore areas close to river mouths around Jakarta Bay, while localized effects of anthropogenic stressors are more important for reef in the central and northern reefs of Kepulauan Seribu (Damar et al. 2012;Farhan and Lim 2012;Baum et al. 2015).
Each island and platform of Kepulauan Seribu is built up from the skeletal debris of coral reefs which rim the platform margins (Jordan 1998). The steep reef front extends seaward from the reef crest to a depth of about 15 m. Well-developed reef facies occur down to a depth of about 16-18 m (Jordan 1998), whereas reef development is restricted to an even shallower depth of 3-10 m around 40 km towards the south in Jakarta Bay (Baum et al. 2015). Reef framework builders at Kepulauan Seribu consist of a diverse fauna of Indo-Pacific corals and encrusting red algae. Below the base of continuous reef development, occurs a transition zone of small reef knolls or coralgal knobs that are surrounded by unconsolidated sediments. In water depths of about 24 m, the seafloor becomes entirely sediment and begins to slope less steeply. Below the reef and near reef zone, in a depth from 27 to 88 m occurs a foraminiferal molluscan packstone facies that is rich in terrigenous mud. The dominant bioclasts are benthic foraminifers, mollusc and coral fragments Fig. 1  1 3 (Jordan 1998). The sedimentation on the shelf further away from the direct influence of the reef system is dominated by siliciclastic mud.
Statistical analyses based on texture and composition reveal four sedimentary facies in the shallow waters (< 20 m) of the Kepulauan Seribu: coral grainstone, coral packstone/ grainstone, coral-mollusc packstone and mollusc wackestone (Utami et al. 2018). The occurrence of mollusc wackestone in the lagoon is controlled by water depth, while the sand apron and reef front do not show significant facies separation with water depth. Environmental facies maps were derived from remote sensing data (Utami et al. 2018): The sand apron can be subdivided into three facies with either no, sparse, or dense seagrass cover. One facies is further characterized by active coral growth. The deeper water zone can be separated into a shallow and a deep subtidal part of lagoons and platform margins. Distribution is therefore mainly controlled by water depth, density of seagrass cover and coral abundance (Utami et al. 2018). The satellite-derived environmental facies correlated only in the lagoon directly with sedimentary facies. No direct correlation of environmental facies to sedimentary facies exists in the sand apron due to the heterogeneity and complexity of the environment (Utami et al. 2018).
Carbonate sediments from the shallow water (< 20 m) zone of the Panggang reef platform, Kepulauan Seribu serve as the basis of this study (Fig. 1). Panggang is an atoll-like reef with an approximate area of 3.2 km 2 with a central lagoon surrounded by reef and sand apron facies. From the reef flat, the water depth increases from about 1 m to about 20 m in the deepest part of the lagoon. The sand apron on the eastern and western side of the lagoon is roughly two times wider compared to the sand apron on the northern and southern side. Only one small inlet traverses the sand aprons allowing very limited connection between the lagoon and the open sea. Panggang Island has the highest population density in Kepulauan Seribu, with more than 1000 people per km 2 .

Methods
A total of 45 surface sediment samples were collected from Panggang reef platform, Kepulauan Seribu. In shallow water, sediments were collected by hand and water depth was measured using a marked stick. Sampling in the lagoon and at the reef slope was conducted by scuba diving and water depth was determined using the diving gear. In the laboratory, samples were washed in fresh water and dried using an oven set to 70 °C for 7 h. Dried samples were first sieved through a 125 µm sieve and split of the coarse fraction was stored in a glass bottle for later point counting analysis. The remaining portion of the > 125 µm fraction was sieved through 250 µm, 500 µm, 1 mm, and 2 mm sieves. All individual grainsize fractions were weighed. Mean grainsize and sorting were calculated using the software Gradistat (Blott and Pye 2001). Point counting analysis was done on the split of the > 125 µm fraction by counting 300 grains under a binocular microscope. The results of these analyses were reported previously in Utami et al. (2018). Here, we present data on the mineralogy, texture and isotopic composition of the same sample set.
Mineralogy was analyzed using a Bruker D8 X-ray diffractometer equipped with a Cu-anode, a secondary graphite monochromator and a scintillation counter. Forty-five bulk surface sediment samples covering all sedimentary facies (Utami et al. 2018) were oven-dried, ground and mounted on sample holders. Additionally, the mineralogy of the < 125 µm fraction from five of these samples and five individual bivalve shells was measured separately. All selected bivalve shells were clean and devoid of incrustations or epibenthos. All measurements were conducted over an angle field of 48° (16°-64°) with a step size of 0.02° per five seconds. Identification and quantification of different mineral phases was achieved by standard Rietveld refinement using the software Profex (ver. 3.14.0). Standard analytical error is ± 2%. A total of 20 weight percent of corundum was added to each sample as an internal standard. This standard was systematically overestimated during Rietveld refinement (avg. 4%), which was attributed to the occurrence of amorphous silica in the form of sponge needles. The MgCO 3 content of calcite (in mol %) was determined using the shift in the position of the d(104) peak (Lumsden 1979). The standard error is ± 0.15% on the mean composition.
Measurements of δ 13 C and δ 18 O values on 15 bulk and five fine fraction samples were made with a Finnigan MAT 253 mass spectrometer connected to a Kiel IV carbonate preparation device at the Leibniz Laboratory, Kiel University. Additionally, nine individual benthic foraminifers (Calcarina, Amphistegina and Sorites) and five individual bivalve shells were measured separately. Bivalve shells were ground to powder before measurement. All Calcarina and Amphistegina, as well as one individual Sorites were split into two aliquots and measured as subsamples. Standard external error is better than ± 0.05‰ for δ 13 C and ± 0.08‰ for δ 18 O (standards NBS 19 and 20). δ 13 C and δ 18 O values are reported in ‰ versus the Vienna PeeDee belemnite (V-PDB). A core of 68 cm length was drilled in September 2005 from a massive Porites coral in the Kepulauan Seribu complex, inshore of Bidadari Island at site 6°01′55″S, 106°44′47″E ( Fig. 1) from 3 m water depth (Cahyarini et al. 2016). δ 13 C and δ 18 O values were measured at the Stable Isotope Laboratory, Vrije University, Amsterdam. For the analysis, the samples were reacted with 100% H 3 PO 4 at 75 °C in a Kiel Device connected to a Finnigan MAT 252 mass spectrometer. Standard external error is ± 0.05‰ for δ 13 C and ± 0.08‰ for δ 18 O (standards NBS 19 and GICS).
A total of 16 selected samples of the < 125 μm fraction were further analysed using scanning electron microscopy in the SEM laboratory of FMIPA (Institut Teknologi Bandung) and the SEM/EDS laboratory of UPP Chevron-Prodi Geologi (Institut Teknologi Bandung). Satellite images were acquired from Google Earth (http://www. googl e.com/earth ). The inverse distance weighting tool of QGIS (ver. 3.14.1-Pi) was used to interpolate data between the sampling points for the generation of maps. All statistical analyses were performed using PAST (Hammer et al. 2001). The Levenés test was used to check if homogeneity of variances was met. The non-parametric Mann-Whitney U test was used for pairwise comparisons of median values. The non-parametric Kruskal-Wallis test for independent samples was used to examine the medians of δ 13 C and δ 18 O values in different sedimentary facies. The significance level for all analysis was set at 95%.

Mineralogy
Sediments in Kepulauan Seribu are mostly composed of biogenic grains, most notably corals and molluscs, with minor components including benthic foraminifers, echinoderms, red algae and the green alga Halimeda. Non-skeletal grains are very rare and restricted to faecal pellets and aggregates cemented by red algae or serpulids (Table 1) (Scrutton 1976b;Utami et al. 2018). SEM investigation confirms that the fine fraction of sediments in Kepulauan Seribu is largely a product of skeletal breakdown (Fig. 2). X-ray diffraction results confirm that surface sediment on the Panggang reef platform in the central part of Kepulauan Seribu is composed mainly of carbonate minerals with small percentages of siliciclastic minerals. Aragonite is the dominant mineral and its contribution varies between 72 and 88%, with a mean of 82% (Table 1, Fig. 3). High magnesium calcite (HMC) and low magnesium calcite (LMC) contents are fairly equal throughout the island, for HMC, the value ranges from 5 to  Fig. 3). The HMC content is very similar in all facies, with an average value of 8.5%. The MgCO 3 content in HMC varies between 13 and 14.2% in mollusc wackestone and coral grainstone facies, respectively. The content of LMC is equal or higher compared to HMC in all facies. The highest aragonite (82%) and lowest LMC contents (9%) are found in the coral packstone/grainstone facies. The mollusc wackestone facies shows the lowest aragonite (77%) and highest LMC (11%) content. The highest siliciclastic content is present in the mollusc wackestone facies (quartz: 2%, smectite: 2%).

Stable isotopes of carbonates
Carbon isotope values (δ 13 C) range from − 1.3 to 0.9 ‰ and oxygen isotopes (δ 18 O) from -3.4 to -4.9 ‰ for bulk sediment (Fig. 6). The fine fraction (< 125 µm) samples did not show a systematic offset in their isotope values with respect to their bulk samples (Fig. 6). Coral abundance is negatively correlated with bulk sediment δ 13 C (r = − 0.52; p = 0.04) and δ 18 O (r = − 0.79; p < 0.01; Table 3), while mollusc abundance shows a positive correlation with δ 13 C (r = 0.58; p = 0.02) and δ 18 O (r = 0.79; p < 0.01). A negative correlation is also apparent for the relationship between the abundance of aragonite and δ 18 O (r = − 0.7; p ≤ 0.01). Aragonite and δ 13 C also show a weak negative correlation (r = − 0.43, p = 0.11). Coral isotope values from Bidadari ( Fig. 6) are centered around − 3 ‰ δ 13 C and − 6 ‰ δ 18 O and lower compared to the bulk samples. Isotope values from bivalves and benthic  Table 1. Mg: mole % MgCO 3 in Mg-Calcite foraminifers are generally higher compared to the bulk samples but partly overlap with samples from the mollusk wackestone facies (Fig. 6). Carbon and oxygen isotope values of bulk sediment from Kepulauan Seribu show more positive value in the inner lagoon and more negative value in the reef front (Fig. 7).

Spatial and compositional pattern in mineralogy
The dominant mineral phase in Panggang is aragonite (  Fig. 3), with slightly increased values close to areas with dense coral cover on the reef front (Fig. 4). The aragonitic red algae genera Galaxaura and Liagora also thrive on the Kepulauan Seribu reef slope (Atmadja 1977), likely also contributing to the high aragonite concentrations. HMC and LMC occur in fairly similar amounts in all facies (Table 1, Fig. 3), with slightly increased values of LMC in the mollusc wackestone facies of the inner lagoon and some isolated points on the western reef flat (Fig. 4). Corals as the main aragonite producers are less abundant in this area, while less light dependent carbonate producers, such as molluscs and echinoderms, are more common (Table 1). The highest LMC content is found in the mollusc wackestone facies (Table 1), which is known to dominate in the inner lagoon where coral grainstone and coral packstone/grainstone facies are absent (Utami et al. 2018). High LMC content in the inner lagoon is likely due to the contribution from molluscs, the most important component in the mollusc wackestone facies (Table 1). Bivalves with calcitic or mixed aragonitecalcite shells in Kepulauan Seribu belong to the families Ostreidae (LMC), Pectinidae (LMC + ARA) and Pinnidae (LMC + ARA) (van der Meij et al. 2009;Ueda 2020). Five individual bivalve shells were measured using XRD. All of them have a mixed aragonite-calcite shell containing on average 90% aragonite and 10% LMC. Locally, LMC is also abundant in the reef flat (Fig. 4), likely due to contributions from bivalves and larger benthic foraminifers (e.g. Amphistegina and Elphidium, Table 4) (Blackmon and Todd 1959) which are abundant in the coral grainstone   Table 2 facies (Table 1). SEM analysis shows that calcareous nannoplankton, as a potential LMC source, is near absent in the sediment. No clear correlation exists between HMC and individual biogenic grain types (Table 3). This is caused by the varied contribution of several different biota to the overall HMC content to the sediment. Red algae, such as Amphiroa, that are typically associated with sea grass in Kepulauan Seribu (Atmadja 1977), contribute HMC to the coral grainstone facies on the reef flat (Table 1). Echinoderms are more important in the deeper water coral-mollusc packstone facies and mollusc wackestone facies (Table 1) compared to other facies. Increased HMC contents in the inner lagoon ( Fig. 4) can therefore be interpreted as a higher contribution from echinoderms. Mg content is generally higher in coralline algae compared to other organisms with HMC skeleton (Chave 1954). This is consistent with the observation that the lowest Mg contents in HMC are found in the mollusc wackestone facies where coralline algae are rare (Table 1, Fig. 3).
The mineralogy of five fine fraction (< 125 µm) subsamples from the coral-mollusc packstone and mollusc wackestone facies was analyzed using XRD. The aragonite content in these fine fraction subsamples is 6% lower compared to their respective bulk samples (Table 2). In contrast, their HMC content is increased by 4% (Table 2). It has been proposed that fish can produce significant amounts of HMC mud (Salter et al. 2012). Most of the important (> 5% abundance) fish families in the area of Kepulauan Seribu (Madduppa 2013) produce Mg-rich amorphous calcium carbonate. This carbonate phase commonly dissolves within hours and therefore is unlikely to be preserved as sediment (Salter et al. 2017). Only one family (Apogonidae) produces mud-sized HMC ellipsoids with Mg content of around 30-45% (Salter et al. 2018). However, the typical crystal morphologies have not been observed during SEM investigations in our samples. It is, therefore, unlikely that fish-derived Mg calcite contributes significantly to the increased HMC values in the fine fraction. In contrast, SEM investigations confirm that the fine fraction is largely a product of skeletal breakdown (Fig. 2). Carbonate components that occur only in the mud to fine sand size (< 125 µm) fraction, such as tunicate spiculae and needles from Padina thalli, are exclusively aragonitic. Laboratory studies suggest  that abrasion rates of corals are much slower compared to dominant HMC components, such as coralline algae and echinoderms (Milliman et al. 1974). Increased HMC content in the fine fraction (Table 2) therefore could be related to preferential breakdown of HMC components. However, it remains unclear in how far these laboratory results can be applied to natural settings. Quantitative data on the composition of different grains size fractions from isolated carbonate platforms of Belize were presented by Gischler and Zingeler (2002). Our calculations based on this data indicate that coral fragments are actually increased in abundance relative to echinoderm and red algae fragments in the fine (< 125 µm) fraction as compared to the larger grain size fraction. This is opposite to what would be expected by the laboratory results on mechanical abrasion (Milliman et al. 1974), possibly indicating that other processes, such as bioerosion, are more important in reef settings. The preferential abrasion of HMC components therefore seems not to be the reason for the increased HMC content in the fine fraction. A known source of mud-sized HMC are red algae growing on sea grass (Land 1970;Perry et al. 2019). Patches of sea grasses are common on the Seribu islands and are associated with reduced grain sizes in the sediment, likely due to their baffling effect (Utami et al. 2018). They also might contribute to the increased HMC contents in the fine fraction. However, a comparison of the mineralogical composition between locations with dense and sparse sea grass cover on the sand apron of Panggang does not show a significant difference with respect to HMC content (Mann-Whitney U test, p = 0.91). Due to the abundance of sea grass, our dataset does not include samples from locations without sea grass.
Alternatively, micritization in northern Belize increases the HMC content of skeletal grains, which later are preferentially abraded and contribute to carbonate muds. Cementation in micro-bioerosion traces is an important mechanism for micritization (Reid and Macintyre 2000). High Mg calcite cement has been observed, besides other mineral phases, in endolithic microborings in corals (Schroeder 1972). Bioerosion and/or micritization are ubiquitous in reefs affected by terrestrial runoff or upwelling in SE Asia (Tomascik et al. 1997;Wilson 2012; Madden and Wilson 2012) including the Seribu Islands (Fig. 2) (Park et al. 1992(Park et al. , 2010. Both processes are most intense at the outer surface of carbonate grains. The higher surface-to-volume ratio of smaller grains could therefore lead to stronger alteration of its mineralogy compared to larger grains. Micritization and bioerosion could therefore have contributed to the increased HMC content in the fine fraction of Panggang (Table 2).
Small amounts of siliciclastic minerals are found in Kepulauan Seribu (Table 1). Quartz content increases with water depth (Table 3), in the inner lagoon and reef front. Quartz and smectite are strongly concentrated in the fine fraction (Table 2), indicating that quartz grains will settle from suspension in low-energy environments but are easily winnowed by waves and currents in high-energy environments. Smectite occurs only in > 6 m water depth within the inner lagoon (Fig. 4), where it likely settles below the wave base. The fact that siliciclastic minerals are found in almost all samples (Fig. 4) shows that Panggang is influenced by terrestrial runoff (Fig. 1). A possible source for the siliciclastics could be Sumatra, since Kepulauan Seribu is relatively isolated from a direct influence from Java and Kalimantan by the regional current system (Tomascik et al. 1997;Jordan 1998). Partially due to this siliciclastic influence, the depth of abundant coral growth in Kepulauan Seribu (0-18 m, Jordan 1998) is typical for equatorial carbonates in SE Asia, but very shallow compared to other regions (0-100 m) with clearer water and lower nutrient input (Wilson 2008).

Spatial patterns in stable isotope composition of sediment
Delineating facies based on their δ 13 C and δ 18 O values is not possible due to a lack of a significance difference between carbonate facies (Kruskal-Wallis test, p = 0.17 and p = 0.15, respectively). However, the average stable isotope values of the coral grainstone facies (δ 13 C = 0.0, δ 18 O = − 4.4) is lower compared to the mollusc wackestone facies (δ 13 C = 0.4, δ 18 O = − 4.0), suggesting some facies control (Fig. 6). Coral abundance shows a significantly negative correlation to δ 13 C and δ 18 O values (Table 3). This results from a strong kinetic control on isotope fractionation in corals (McConnaughey 2003). Coral skeletons therefore show a depletion of heavy isotopes and consequently a strong offset from equilibrium values towards more negative isotope ratios (Fig. 6). Corals exert a strong influence on the bulk isotope composition, since they are the predominant skeletal material and major aragonite producer (Fig. 6). Marine mollusc shells, as the second most important sedimentary component (Table 3), do not show kinetic effects in their isotopic composition (McConnaughey and Gillikin 2008). Their δ 18 O values are often close to equilibrium values (Grossman and Ku 1986) and therefore more positive compared to corals (Fig. 6) (Aharon 1991). Carbon isotope signatures in molluscs are complex, with near equilibrium values e.g. in giant clams (Aharon 1991) but often several ‰ lighter due to metabolic carbon contribution (McConnaughey and Gillikin 2008). In reefal carbonates, molluscs are typically less depleted in 13 C compared to cooccurring corals (Fig. 6) (González and Lohmann 1985). This results in a positive (negative) correlation (Table 3) between oxygen and carbon isotopes and molluscs (corals). Higher stable isotope values in the inner lagoon and sand apron of Panggang (Fig. 7) result from the increased abundance of molluscs in mollusc wackestone and coral packstone/grainstone facies (Table 1) and coincides with the overall lowest aragonite value (Fig. 4). The most negative δ 13 C and δ 18 O values occur at the reef front in coral grainstone and coral packstone/grainstone facies where coral cover is densest (Fig. 7). However, the compositional differences between individual facies types (Table 1) are not strong enough to result in separate isotope fields for each facies (Fig. 6).
The most important group of purely calcitic components are benthic foraminifers. Many shallow water benthic foraminifers from reefal environments precipitate their skeleton close to oxygen isotope equilibrium (± 1‰) (Saraswati et al. 2004;Maeda et al. 2017), although slightly stronger depletion of 18 O was observed as well (Wefer and Berger 1991). Carbon isotopes typically are lighter compared to equilibrium values with common 13 C depletion on the order of 1-3‰ for most miliolids and 2-5‰ for rotaliids (Wefer et al. 1981;Saraswati et al. 2004). This is consistent with lower δ 13 C values in rotaliids (Calcarina, Amphistegina) compared to miliolids (Sorites) in the study area (Fig. 6). Benthic foraminifers from reefal environments therefore are characterized by generally more positive stable isotope values compared to cooccurring corals (Fig. 6), but can show partial overlap in their δ 13 C ranges (González and Lohmann 1985;Gischler et al. 2012), especially if rotaliids are abundant. The observed negative correlation between aragonite content and stable isotope values (Table 3) likely results from this isotopic difference between corals and benthic foraminifers as the main aragonite and calcite producers, respectively. The fact that some molluscs contribute calcite to the sediment likely enhances this effect.
The positive correlation between stable isotopes and fines is a spurious relationship (Table 3), resulting from the strong positive (negative) correlation with molluscs (corals). This is supported by the observation that the fine fraction samples show no consistent isotopic difference to the bulk samples (Table 3). This similarity indicates that the fine fraction is genetically related to the coarse fraction and is formed by skeletal breakdown.

Comparison with other carbonate systems
To facilitate the comparison of mineralogical and stable isotope signatures between humid equatorial, arid (sub) tropical, and temperate carbonate systems we have grouped them geographically (Figs. 8 and 10). Carbonates from each geographical group share common characteristics and are deposited under similar environmental conditions (Fig. 9, Table 5). Carbonate systems from the tropical western Atlantic (TWA) and Persian Gulf represent mostly (sub)tropical carbonates (sensu Wilson 2002Wilson , 2012) characterized by at least seasonally arid conditions. These carbonates often contain facies with non-skeletal grains and are deposited in normal marine to hypersaline waters (Table 5). Exceptions are Florida Bay and Belize, which are characterized by riverine input leading to strong gradients in salinity ranging from brackish to hypersaline and absence or low abundance of non-skeletal grains. The carbonates from the Indo-Pacific Warm Pool (IPWP) (Figs. 8,9,10,Table 5) are deposited on isolated carbonate platforms under the highest average temperatures (> 28 °C) and high seasonal precipitation, but outside the direct influence of river runoff. Non-skeletal grains are uncommon. High latitude systems (Figs. 8,9,10,Table 5) are typical temperate water Heterozoan carbonates (Wilson 2012).
The δ 13 C signature of carbonates from Kepulauan Seribu is considerably more negative compared to most other platform and reef systems (Fig. 8). The carbon isotopic composition of marine carbonates is mainly influenced by the δ 13 C of the dissolved inorganic carbon (DIC) in water (Smith and Kroopnick 1981), the mineralogy (Rubinson and Clayton 1969;Romanek et al. 1992) and so called "vital effects" in carbonate secreting organisms (e.g. kinetic and metabolic effects) (McConnaughey 2003;McConnaughey and Gillikin 2008). Because of these vital effects skeletal carbonates are typically depleted in δ 13 C compared to inorganically precipitated carbonates from the same environment Gischler et al. 2012). Non-skeletal grains and inorganically precipitated muds are much more common in modern platforms and reefs from the tropical, western Atlantic (TWA) and Persian Gulf compared to carbonate systems from the Indo-Pacific Warm Pool or New Zealand (Table 5) (Milliman et al. 1974;Gischler 2011;Gallagher et al. 2018). Most carbonates from the TWA and Persian Gulf therefore show broader δ 13 C ranges encompassing more positive δ 13 C values (Fig. 8).
Due to mineralogical effects, inorganic aragonite is enriched in δ 13 C by about 1.8 ‰ compared to coprecipitated calcite (Rubinson and Clayton 1969;Romanek et al. 1992). This mineralogical effect should favor a positive correlation between aragonite content and δ 13 C. However, the actual relation depends on the source and therefore isotopic signature of the aragonite. Aragonitic non-skeletal grains and green algae (Fig. 8) are characterized by more positive δ 13 C values compared to most calcitic skeletal grains, resulting in a positive correlation between both parameters Gischler et al. 2012). Corals on the other hand are more depleted in δ 13 C compared to most calcite secreting organisms (Fig. 6), which explains the negative correlation between aragonite and δ 13 C in the coral dominated Kepulauan Seribu (Table 3) and the lower δ 13 C compared to the relatively calcite-rich cool water carbonates (Fig. 8). However, the δ 13 C signature of carbonate from Kepulauan Seribu is also more negative compared to most other coral dominated reef facies from the IPWP (Fig. 8).
Carbon isotope data from a Porites coral from Bidadari Island (Fig. 1) in the Java Sea, is on average more depleted compared to Porites corals from the Seychelles and Chagos both of which are located within the IPWP (Fig. 8) (Pfeiffer et al. 2004;Pfeiffer and Dullo 2006;Cahyarini et al. 2016). It is unlikely that this difference is solely due to kinetic related offsets, since all three coral colonies show similar growth rates (Pfeiffer et al. 2004;Pfeiffer and Dullo 2006;Cahyarini et al. 2016). More negative δ 13 C values of the Porites colony from Bidadari Island are more likely due to the influence of river run off. Corals have the potential to record land-ocean carbon transfer because they draw on dissolved inorganic carbon (DIC) for calcification. The δ 13 C of DIC in coastal water is influenced by the δ 13 C signature from adjacent rivers (Moyer and Grottoli 2011). Bidadari Island is located in the Java Sea (Fig. 1) and in close proximity to three major rivers in the Greater Jakarta Area (Farhan and Lim 2012). Estuaries in SE Asia show low δ 13 C DIC values due to the contribution of 12 C from riverside mangroves (Miyajima et al. 2009). The low δ 13 C of the coral skeleton in Kepulauan Seribu therefore likely reflects the enrichment in 12 C of ambient seawater due to the input of riverine DIC to the Java Sea.
To a somewhat lesser extent, light intensity (photosynthesis) also has been shown to influence the value of skeletal δ 13 C in corals (Moyer and Grottoli 2011). Corals are usually depleted in 13 C relative to ambient seawater as a result of kinetic and metabolic fractionation (Swart 1983;McConnaughey 1989;Allison et al. 1996;Grottoli and Wellington 1999). Environmental variables that influence coral metabolism should therefore also affect δ 13 C levels in coral skeletons (Grottoli 2002). Under low light conditions photosynthesis decreases and metabolic fractionation increases with the consequence of decreased skeletal δ 13 C (Grottoli 2002). Sediment input from rivers into the Java Sea potentially enhances turbidity in the semi-enclosed basin, leading to a reduction of light intensity and photosynthesis thus resulting in reduced skeletal δ 13 C values. This would be consistent with the relatively shallow depth of reef growth  (Swart et al. 1989;Andrews 1991), Kuwait ramp (Gischler and Lomando 2005), and an example from high latitudes; Wanganui shelf, New Zealand (Gillespie and Nelson 1997). For comparison the isotopes ranges of faunal elements are shown: Porites coral from Bidadari Island, Kepulauan Seribu (Cahyarini et al. 2016), Porites corals from Seychelles and Chagos (Pfeiffer et al. 2004;Pfeiffer and Dullo 2006), Halimeda bioherm, eastern Java Sea (Aharon 1991), bivalves and benthic foraminifers from Panggang reef platform. Location of Bidadari Island is shown in Fig. 1 (Jordan 1998) and the presence of siliciclastic minerals in sediment samples below the wave base (Fig. 4) in the study area. Turbidity can also be produced by algae blooms caused by increased nutrient input either through river water or untreated sewage from populated islands (Baum et al. 2015). Additionally, higher nutrient availability could influence the δ 13 C of coral skeletons by an increase in the heterotrophic feeding rate or an increase in the photosynthesis rate related to higher zooxanthellate concentrations (Grottoli 2002). However, these processes have opposing effect on the δ 13 C of coral skeletons. The net effect therefore will depend on the relative magnitudes of these various influences but is expected to be small (Heikoop et al. 2000).
The largest single anthropogenic influence on the isotope signature of recent corals likely is caused by the input of anthropogenic CO 2 to the global surface ocean. The accelerated admixture of fossil fuel derived CO 2 leads to a shift towards lower δ 13 C values in atmospheric CO 2 , the so called Suess effect (Keeling 1979). This Suess effect is often recorded in declining coral skeleton δ 13 C values over the twentieth century (Swart et al. 2010). Linsley et al. (2019) showed that the Suess effect can shift skeletal δ 13 C in corals by more than 1‰ since the 1950s. The carbon isotope composition of corals in Fig. 8, which are largely derived from the twentieth century, is therefore likely depleted relative to bulk sediments or sedimentary grains of mixed and/ or unknown age.
Florida Bay is the only location that shows a nearly complete overlap but broader total range in δ 13 C values compared to the sediments from Kepulauan Seribu (Fig. 8). In contrast to many other areas in the TWA, non-skeletal grains and inorganic mud are lacking in Florida Bay (Table 5). Average δ 13 C DIC of surface waters in Florida Bay is lowest in areas influenced by isotopically depleted river runoff from the Everglades, while the δ 13 C variability in other areas is controlled by remineralization of marine organic matter (Swart and Price 2002). Similar to Kepulauan Seribu, the lowest values in δ 13 C of carbonates deposited in Florida Bay  (Völpel et al. 2017). Location of carbonate sites: 1. Study area in Kepulauan Seribu, 2. Enewetak, 3. The Maldives, 4. Chagos, 5. Seychelles, 6. Kuwait ramp, 7. San Salvador, Bahamas, 8. Great Bahama Bank, 9. Florida Bay, 10. Belize, 11. Wanganui shelf, New Zealand seem to be a result of isotopically depleted river runoff and a lack of inorganic carbonate precipitation.
Since the pioneering work by Urey (1947), the δ 18 O values of carbonate minerals are known to be controlled by temperature of formation. Other controlling factors are the δ 18 O value (δ 18 O w ) of the precipitating fluid (Urey 1947;Epstein and Mayeda 1953), mineralogy (Tarutani et al. 1969), and "vital effects" in carbonate secreting organisms (McConnaughey 2003). The slope of the δ 18 O/temperature relationship in carbonates is commonly assumed to be in the range from − 0.2 to − 0.25 ‰/°C (Grossman and Ku 1986;Bemis et al. 1998), resulting in a decrease of 0.8-1 ‰ δ 18 O for a four degree increase in temperature.
Carbonates from Kepulauan Seribu are characterized by much lower δ 18 O values compared to other carbonate systems (Fig. 8). This can be explained by the combined effect of high sea surface temperatures and low δ 18 O w in the Java Sea (Fig. 9). Modelling results confirm that for carbonates precipitated in equilibrium with seawater, SE Asia is globally the region with the most negative δ 18 O signature (Fig. 9c). Within the IPWP the differences in temperature are small (Fig. 9a), but δ 18 O water is significantly reduced in SE Asia compared to e.g. the central equatorial Indian Ocean or Pacific (Fig. 9b).
Differences in hydrology are reflected in the lower δ 18 O values of a Porites coral from Kepulauan Seribu compared to those from the Seychelles and Chagos (Fig. 8). The difference of more than 1 ‰ in δ 18 O is more than can be explained by intercolony offsets (Sayani et al. 2019) or the temperature difference of ~ 1 °C between the Java Sea and  Gillespie and Nelson (1997) the other two sites. Most of the difference therefore can be attributed to the contrasting δ 18 O seawater values, which based on the difference in salinity should be around 0.9-1.1 ‰/δ 18 O assuming a slope of 0.22 ‰ for the salinity/ δ 18 O relationship (Fairbanks et al. 1997). The lower δ 18 O seawater values in the Java Sea result from heavy precipitation and increased runoff, and are also influenced by the δ 18 O w of the precipitation itself via the so called "amount effect" (Rozanski et al. 1993;Lau and Yang 2003). Facies also could influence the δ 18 O signature of bulk sediments due to differences in the "vital effects" of carbonate producing organisms (McConnaughey 2003;McConnaughey and Gillikin 2008) and mineralogy (Tarutani et al. 1969). However, the different δ 18 O values for individual grain types to a large extent seem to average out for bulk samples (Gischler et al. 2012). This results in relatively large overlaps between facies in most atolls, barrier reefs and platforms (Gischler et al. 2012), showing that the influence of facies in these settings is limited.
The carbonate mineralogy of shallow water (< 20 m) bulk sediments in the tropics and subtropics (IPWP, TWA) is predominantly composed of aragonite and HMC, with LMC usually present in low quantities < 10% (Fig. 10) (Ginsburg 1956;Purdy 1963;Weber and Schmalz 1968;Gischler 2006;Reijmer et al. 2009;Gischler et al. 2013;Schmitt and Gischler 2017;Husinec et al. 2019). The only exception is the Kuwait ramp system that shows LMC contents of up to ~ 40% (Gischler and Lomando 2005). In reef atoll and platform systems higher LMC content > 10% are usually spatially restricted and related to local dominance of rotaliid foraminifers, such as Amphistegina (Weber and Woodhead 1972), contribution from calcite precipitated in meteoric waters (Swart and Kramer 1997) or from microbial activity (Pusey 1975;Purdy and Gischler 2003). Otherwise, higher LMC (7-14%) contents generally occur only in deeper water between 20 and 70 m, mainly contributed from planktonic foraminifers, coccoliths and molluscs (Weber and Schmalz 1968;Gischler 2006). Aragonite also dominates the mineralogy of the study area in Kepulauan Seribu but HMC and LMC occur in nearly equal proportions. The main mineralogical difference to other reefal systems in the IPWP is the occurrence of clastic mineral phases (Table 1) and a higher LMC content (~ 10%) in the shallow water sediment (Fig. 10). Opposed to open ocean settings, the semi-enclosed Java Sea is influenced by river runoff. Terrestrial sediments reaching Kepulauan Seribu (Table 1) most likely are derived from Sumatran rivers, since oceanic currents isolate Kepulauan Seribu from sediment input from Java and Kalimantan (Tomascik et al. 1997;Jordan 1998). The suspended particulate matter promotes oligophotic conditions in relatively shallow water. Wilson and Vecsei (2005) argued that oligophotic conditions in water depth below 20 m are typical for humid equatorial carbonate systems in SE Asia and would favor the deposition of e.g. relatively calcite-rich foramol facies. Foramol facies in Kepulauan Seribu dominates in water depth below 27 m (Jordan 1998). In water depth < 20 m, rotaliid foraminifers contribute to LMC (Table 4) but their abundance is not high enough to explain the generally increased LMC content. Another important source of LMC in Kepulauan Seribu seems to be mixed aragonite-calcite bivalve shells. More mineralogical analysis from humid equatorial settings are needed to test if the increased LMC content observed in Kepulauan Seribu is a common feature of shallow water (< 20 m) reef systems in humid, equatorial settings.
Fine fraction sediment from the Panggang reef platform is showing a higher HMC compared to the bulk fraction, likely related to micritization and bioerosion. Higher HMC contents in the fine fraction are also found in Belize, the potential source is micritization of skeletal grains to Mg calcite (Weber and Schmalz 1968). Micritization of skeletal grains is also an important diagenetic process in other shallow marine carbonate platforms, e.g. The Bahamas where the grains are altered to aragonite (Purdy 1963;Bathurst 1966;Kendall and Skipwith 1969;Trumbull 1988). Preferential alteration of the carbonate grains into high Mg calcite is potentially controlled by lower aragonite supersaturation in water off Kepulauan Seribu. Fresh water influx from rivers promotes low salinities in the semi-enclosed seas of SE Asia, leading to relatively low aragonite supersaturation compared to more arid regions as the tropical, western Atlantic (Park et al. 2010). Seawater Mg:Ca ratios are another important control on carbonate mineralogy, with lower Mg:Ca favouring calcite precipitation at a given temperature (Morse et al. 1997). A recent compilation of Mg:Ca ratios in seawater shows that values are not spatially constant but influenced by river waters with low Mg:Ca values (Lebrato et al. 2020). No seawater data is available for the semi-enclosed seas of SE Asia, but the high riverine input could potentially also reduce the Mg:Ca ratios in the Java Sea. Depending on the magnitude of this effect, the reduced Mg:Ca ratios could facilitate the precipitation of calcite in the study area.
The mineralogical composition in the shallow water (< 20 m), photozoan dominated carbonates from Kepulauan Seribu overlaps with the most aragonite rich part of the cool water carbonate field (Fig. 10). However, based on the dominance of reefal material, the sediments can be clearly differentiated from heterozoan dominated cool water carbonates. This might be more difficult for the foramol facies that dominates at the toe of slope and below (Jordan 1998). Wilson (2012) pointed out that oligophotic, humid equatorial carbonates share many similarities with cool water carbonates. However, a marked difference is that HMC components in Panggang have much higher Mg-contents (11.7-14 mol% MgCO 3 ), compared to 4-12 mol% MgCO 3 in carbonate sediment of e.g. South Australia (James et al. 2005;O'Connell and James 2015). This difference can be best explained by the different water temperature between both regions. The MgCO 3 content in HMC of calcareous marine organism is strongly influenced by the water temperature, with an increase in Mg-content related to higher temperatures (Chave 1954). However, this difference in Mg-content will not be helpful in the distinction of ancient carbonates since HMC is losing Mg from the crystal lattice relatively early during diagenesis (Brachert and Dullo 2000).

Conclusion
The reef platforms of Kepulauan Seribu represent the shallow water, photozoan-dominated part of a humid equatorial carbonate system. The mineralogical and/or isotopic characteristic of their surface sediments are distinct from most other carbonate systems and likely are representative for many SE Asian reef systems affected by river runoff.
Biogenic debris from corals and molluscs are the most important grain types. Minor components include benthic foraminifers, echinoderms, red algae and the green alga Halimeda, while non-skeletal grains are very rare. Aragonite is the dominant mineral phase, mainly contributed from corals and to a lesser degree by molluscs. LMC and HMC occur in fairly similar amounts in all facies. Bivalves (calcitic and bimineralic) and some rotaliid foraminifers (e.g. Amphistegina and Elphidium) contribute LMC. Most red algae (e.g. Amphiroa), echinoderms, and benthic foraminifers (e.g. Calcarina and Operculina) are important sources for HMC. An increased HMC content in the fine fraction could be related to preferential breakdown of HMC components, possibly from red algae growing on sea grass. Alternatively, micritization and bioerosion could contribute to increased HMC. Small amounts of siliciclastic minerals (quartz and smectite) are preferentially deposited in deeper waters below the wave base. This shows that despite its relatively isolated position, Kepulauan Seribu is affected by sediment-laden river runoff.
The main mineralogical difference between carbonates from Kepulauan Seribu and most other reefal systems in the tropic western Atlantic and Indo-Pacific Warm Pool area, is the widespread presence of > 10% LMC content in the shallow water sediment. The LMC is contributed by organisms that are not or less light dependent (bivalves and some rotaliid foraminifers) and therefore can thrive under oligophotic conditions. The suspended particulate matter from terrestrial runoff likely limits the depth of the euphotic zone. This is reflected in the deposition of foramol facies below a depth of ~ 28 m (Jordan 1998), but seems to already affect the sediments in the upper 20 m. Compared to similar foramol facies from higher latitudes, the Mg content in HMC components from Kepulauan Seribu is much higher due to the different water temperature between both regions.
The δ 13 C signature of carbonates from Kepulauan Seribu is considerably more negative compared to nearly all other platform and reef systems. This is likely due to the combination of low δ 13 C DIC in surface waters and the lack of green algal grains and inorganically precipitated carbonates. The δ 13 C in surface waters likely is reduced by 12 C-enriched input of DIC to river water from riverside mangroves. Corals as the predominant skeletal material and major aragonite producer also exert an influence on the bulk isotope composition, resulting in depletion of heavy isotopes values due to kinetic "vital effects". Low light intensity due to turbid water in the semi-enclosed basin of the Java Sea potentially hampers photosynthesis, which further reduces skeletal δ 13 C values in corals due to metabolic effects. The Suess effect, which is caused by the combustion of fossil fuel, leads to a decline of skeletal δ 13 C over the twentieth century in corals. This signal might not be recorded to its full extent in the bulk sediment yet, which represents a grain mixture of different ages.
Carbonates from Kepulauan Seribu are also characterized by much lower δ 18 O values compared to other carbonate systems, due to high sea surface temperatures and low δ 18 O w in the Java Sea. Deep atmospheric convection leads to heavy precipitation and increased fluvial runoff of very 18 O-depleted water to the semi-enclosed Java Sea. The relatively low δ 18 O of benthic foraminifers and bivalves from Kepulauan Seribu shows that the oxygen isotope signature can be used to distinguish foramol facies from oligophotic equatorial carbonates and cool water carbonates in the rock record. commenting an earlier version of the manuscript. Dr. Nils Andersen (Kiel) is acknowledged for stable isotope analysis. Uwe Wollenberg is thanked for his help with XRD measurement. We thank two anonymous reviewers for insightful comments on the manuscript. This work was partially supported by National Geographic Explorer Grant No CP 087R 17. We acknowledge Google Earth for online imagery and QGIS (ver. 3.14.1-Pi) for GIS software.
Funding Open Access funding enabled and organized by Projekt DEAL.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.