1. Introduction
Coastal bays are broadly divided into estuary bays and semi-enclosed bays. An estuary bay connects an estuary with the open sea and is thus exposed to circulation of ocean currents and to freshwater inputs from river flows. A semi-enclosed bay receives little inflow of freshwater from land and the circulation of seawater to and from the open sea is limited by one or several restricted channels. Accordingly, seawater circulation, sedimentation processes, and ecosystem and sediment characteristics will differ depending on the type of bay [
1].
The west coast of the Korean Peninsula is an area of shallow water and a wide continental shelf. Its estuaries, bays, and tidal flats were formed by sea-level rise after the Pleistocene. This led to the formation of a ria coast that accounts for the various types of coastal environments, although meteorological events, such as storms and typhoons, have also determined the coastal features [
2,
3]. The large tidal flat consists of sediments that vary in their type and their characteristics [
4].
In the 1960s, a period that initiated the rapid growth of developing countries, including South Korea, large areas of the tidal flats were lost due to coastal development, including agriculture, industry, new town construction, and large-scale land reclamation projects. These environmental changes caused hydraulic changes in the tidal flats, adversely affecting their sedimentary environments and ecosystems, such as by the erosion/deposition of coastal sediments due to increased coastal pollution resulting from inflow of pollutants from land and changes in ocean currents and tides due to changes in the coastline [
1].
The geographical features of the study area include the estuaries of the Han River and the Geum River, and the interior of bays, which together with coastal development determine depositional mechanisms [
1]. The unique characteristics of the west coast of Korea have been examined in sedimentological studies [
5,
6,
7,
8], including those of surface sediments [
1,
9,
10], and long-term changes in the sediments and topography of the western and southern coasts [
8,
11,
12]. These studies have shown regional hydraulic changes that affect the transition environment, including the sedimentary environments of the tidal flats and therefore their biodiversity as well as their ability to supply nutrients to nearby waters and to remove pollutants [
13].
Cheonsu Bay is located on the central coast of the Yellow Sea. The Hwangdo tidal flat is an intertidal landform located in Cheonsu Bay, Taean-gun, Chungcheongnam-do, on the west coast of Korea. Seawater enters the bay and the waterway under the Anmyeon-do suspension (land) bridge. Prior to the 1980s, the bay was free of contaminants, as there was no inflow of large rivers, and biological productivity was high, making it suitable for aquaculture. However, the construction of seawalls near the study area from 1983 to 1985 resulted in the creation of Lake Ganwol and Lake Bunam, which caused significant changes in the sedimentary environment and ecosystem of Cheonsu Bay. After construction of the embankment, the absolute flow velocity in the northern part of the bay decreased greatly, and the exposure time of the northern tidal flat increased [
14].
The transition environment before and after construction of the embankment in Cheonsu Bay has been extensively investigated [
14,
15,
16,
17]. In 1979, Kim examined the sedimentation and sedimentary environment of Cheonsu Bay before construction of the embankment [
18]. Other researchers have focused on the post-construction characteristics and the distribution of the tidal flats, sediments, and waterways [
19,
20,
21,
22]. The surface sediments consisted of five sedimentary facies, as determined in June and October 2003, and sediment accumulation rates in the Hwangdo tidal flat over a period of 11 months indicated sediment deposition in the central part, whereas the eastern and western sides of the tidal flat had eroded [
1]. An analysis of high-spatial-resolution remotely sensed data (IKONOS and KOMPSAT-2, Korea multi-purpose satellite-2) showed that the percentage of sand grains larger than silt (0.0625 mm) at the Hwangdo tidal flat had increased considerably since the early 2000s. Mud-flat facies decreased by 5.81% in the late vs. the early 2000s, whereas mixed mud flats and sand flats increased by 4.46% and 2.14%, respectively [
23].
This study examined recent dynamic changes in the tidal flats in north Cheonsu Bay. Specifically, geochemical sediment characteristics (n = 107) measured in 2022 in the Hwangdo tidal flat were analyzed to determine the interaction of tidal-flat exposure with various environmental factors, and their correlation (e.g. grain-size distribution, mean size, gravel + sand + silt + clay compositions, sorting, skewness, kurtosis, dry bulk density, pore salinity, porosity, water content, TC, TOC, TN, TC, and TIC in surface sediments)
2. Materials and Methods
The Hwangdo tidal flat of Cheonsu Bay, Taean-gun, Chungcheongnam-do is located on the central coast of the Yellow Sea and is characterized by a diverse and complex environment and topography [
1]. Seawater enters the bay and the waterway under the Anmyeon-do suspension (land) bridge. Cheonsu Bay surrounds Anmyeondo Island in Taean County, Ganwoldo Island in Seosan City, and Yellow Sea line in Boryeong City, as shown in
Figure 1 (a) and (b). The total surface area was originally 380 km
2, but it was reduced to ~180 km
2 following the construction of an embankment as part of the Seosan A and B District reclamation project in 1984, the creation of Buam Lake and Ganwol Lake, and the construction of a second embankment as part of the Hongbo District reclamation project [
14]. The Cheonsu Bay has a semidiurnal tide, with a maximum high tide of 6.33 m, a minimum low tide of 2.86 m, and an average tide of 4.59 m. The 1.65-km wide × 5.15-km long Hwangdo tidal flat is adjacent to Anmyeon Island in Cheonsu Bay and it includes complex tidal channels. The sedimentary facies in the study area were composed of mud, mixed, and sand facies from the high tide to the low tide shoreline [
24]. The sampling stations for the surface sediments (n = 107) are shown in
Figure 1 (c) and (d).
The Hwangdo tidal flat is gentle and broad with an average slope of < 1 mm and a highly varied topographical structure that includes tidal channels ranging in width from 0.02 m to 2 m, chenier, and sand shoals. These topographic features affect the tidal flow and thereby the sediment-distribution characteristics and the topographical elevations. The latter affect the optical reflectivity by altering the exposure time depending on the state of the tide [
23]. In this study, the terrain altitude of the Hwangdo tidal flat was examined using Leica's Real Time Kinetic-Global Positioning System (RTK-GPS).
A drone survey was conducted using a full-frame sensor (Matrice 300 drone RTK-GPS + GLONASS + BeiDou + Galileo; DJI) as shown in
Figure 2 (a) and (b). Orthophotos were taken along the route and ground photos were taken from ground-control-point (GCP) stations (n = 17) using the Matrice 300 drone. In addition, tide gauges were installed at Hwangdo Bridge in
Figure 2 (c). The 3D error of each vertex was set to be within 0.02 m. RTK-GPS was operated using both a reference station and a mobile station. At the former, the reference point was set using the coordinates of the integrated reference point provided by the National Geographic Information Institute, South Korea. A moving RTK was obtained using a pole with a length of 2 m that supported a flat object in the center of the apex surface layer, with the end of the pole fixed to the flat object to measure the orthometric height (
Figure 2 (b)). The mean change in tidal height measured at Hwangdo Bridge during the sampling period was ~7.53 m (minimum: –3.86 m, maximum: +3.67 m).
Samples (~1 kg) were collected from the sediment surface (at depths of up to 10 mm) using a clean stain-less steel spoon and stored in the dark in an icebox (< 4°C) until delivery to the laboratory (spring range) for grain-size distribution of sediment particles, pore salinity, porosity, water content, TC, TOC, TN, TC, and TIC in surface sediments. The larger tidal range occurs during spring tides (spring range) at Hwangdo Bridge from October 25 to October 28, 2022 is shown in
Figure 3 (a), (c), and (d). Hwangdo Bridge was constructed in December 2011 by removing the embankment connecting Anmyeon-do and Hwang-do, Taean-gun in
Figure 3 (b). After the opening of the bridge, the seawater distribution resulting from the bridge construction improved the marine environment, including an increase in the benthos biotic productivity of the nearby tidal flats [
25].
The grain-size distribution of sediment particles was analyzed following a standard sieving procedure [
26] to obtain the sand fraction (particle size < 4 Φ). The mud fraction (> 4 Φ) was obtained by examining samples at 0.5-Φ intervals using a particle-size analyzer (SediGraph 5100; Micromeritics, Norcross, GA, USA) at the Korea Institute of Ocean Science and Technology (KIOST). An inclusive graphic method was used to determine sediment type, mean size (Φ), sorting (Φ), and skewness and kurtosis [
27].
For total organic carbon (TOC), 10 mL of 1 N hydrochloric acid (HCl) was added to 0.5 g of powdered sediment overnight to remove carbonate materials. The samples were then re-dried and re-weighed, with the difference recorded as the TOC for weight correction through acid (1 N HCl) treatment. For total carbon (TC) and total nitrogen (TN) measurements, unacidified samples were used. Acidified TOC, unacidified TC and TN were measured using an elemental analyzer (FlashEA 1112; Thermo Fisher Scientific, Waltham, MA, USA) at KIOST, by placing a few milligrams of treated sediments in an aluminum capsule. To ensure the reliability of the analytical data, two soil reference material (SRM) from the National Institute of Standards and Technology of the United States (nos. 338-40025 and 338-35210, Thermo Fisher Scientific) were included in the analysis. Carbon and nitrogen contents of the SRM analyzed were 2.30% and 0.21% with standard deviations of 0.02 and 0.01, respectively. And relative standard deviations (RSD) were 0.01 and 0.06, respectively, representing a high precision. The total inorganic carbon (TIC) content was calculated as TC (%) – TOC (%), and the total calcium carbonate content, based on the relative percent difference (%), as shown this is the example 1 and 2 of equation: [
28,
29].
To determine porewater salinity, a sediment sample of ~50 g was placed in a 50-mL Falcon tube and centrifuged at 2500 rpm for 10 min. The supernatant was then used to measure the electrical conductivity.
Since the sediments of tidal flats are always in contact with seawater, a sediment analysis must also consider the seawater in the pores between particles (interstitial water content, IWC) and the seawater on the surface of the tidal flat (surface water cover, SWC) [
30]. The nature of the seawater present in tidal flats will vary depending on the sediment particle size, the topography (slope) of the tidal flat, the exposure time,
etc., and will affect the optical reflectivity [
31]. In this study, samples for SWC analysis were not collected, while IWC was measured as the difference data between in weights change of the sediment before and after drying at 100°C. The moisture content was measured based on the weight change. Dry bulk density (g/cm
3) was defined as the ratio of the weight of the dry sample (W
d) and the total sample volume (V
t) and represents the mass of dry solids in a given bulk volume (5 ml) of sediment [
32,
33].
In flooded sediments, this is equivalent to the volumetric percentage of water in the sediment [
33]. Another important characteristic of sediments is their porosity (ρ), defined as the percentage of the total volume not occupied by solid particles, as shown in 3 of equation:
where V
w is the volume of water, V
t is the total volume of a sediment sample, and V
s is the sum of the volume of the solids.
The porosity (ρ) of the sediment can be also calculated based on the relationship between its dry bulk density and the particle density [
34], as shown in 4 of equation:
On October 26, 2022, a drone survey was conducted. The DEM was improved using the actual tidal data determined at Hwangdo Bridge, the drone data, and the GCP during the study period. The exposure time was calculated according to the tide-level change. The quality of the final processed data of the DEM was 1 m2 per DEM pixel. Data from 2022 were obtained by correcting the tidal-gauge data according to Korea Hydrographic and Oceanographic Agency (KHOA) in the Boryeong's sea level measurement tide area.
In the multivariate analysis (MVAs), environmental variables were square-root-transformed and a cluster analysis and a non-metric multidimensional scaling (MDS) analysis were performed using an Euclidean distance matrix [
35]. A MVAs is based on the multivariate statistics and is used to address situations in which multiple measurements of each experimental unit are made and the relationships among those measurements and their structures are important. MVAs include normal and general multivariate models and distribution theory, the study and measurement of relationships, probability computations of multidimensional regions, and exploration of data structures and patterns. In this study, a principal coordinate analysis (PCA) was applied to determine the major abiotic parameters influencing the environmental variables. MVAs were performed using Primer 7 + PERMANOVA statistical software (Plymouth Marine Laboratory, Plymouth, UK).
3. Results
In total, 13 surface facies were distributed throughout the study area. The sedimentary facies were classified following Folk and Ward [
27]. Based on 107 sample groups, four groups were identified. Group 1 (n = 10) consisted of slightly gravelly sand ((g)S), gravelly sand (gS), and sand (S); group 2 (n = 73) consisted of slightly gravelly muddy sand ((g)mS), gravelly muddy sand (gmS), muddy sand (mS), muddy sandy gravel (msG), and silty sand (zS); group 3 (n = 20) consisted of slightly gravelly sandy mud ((g)sM), sandy mud (sM), and sandy silt (sZ); and group 4 (n = 4) consisted of gravelly mud (gM) and mud (M). The sedimentary facies ranged from coarse sand (sand : mud ratio = 9:1) to sandy mud, in proportions of sandy silt (sZ) > muddy sand (mS) > slightly gravelly sandy mud ((g)sM). The respective mean values (range:) of gravel, sand, silt, and clay were 1.15% (0.00%–30.70%), 59.87% (6.98%–96.64%), 25.11% (0.15%–55.84%), and 13.86% (0.30%– 39.42%). The average sediment grain size (n=107) was +4.35 Φ (−0.31 to +7.33 Φ). The average sorting was +2.27 Φ (+0.61 to +4.52 Φ), indicative of moderately well sorted to extremely poorly sorted sediment grains. The average skewness was +0.55 (–0.26 to +0.76), with a positive value (tail to the right) when the proportion of coarse sediment was large. A range of kurtosis values of +1.90 (+0.51 to +5.59) is of the leptokurtic form if it is largely centrally concentrated. This shows medium to good and extreme characteristics in the mean sorting degree. The distribution of surface sediments was generally fine-grained sandy mud sediments in the west, muddy sand sediments in the central region, and coarse gravelly sand sediments in the east.
The average TOC, TC, TN, and TIC contents were 0.17% (range: 0.02–0.61%), 0.30% (0.13–1.16%), 0.02% (0.00–0.08%), and 1.06% (0.01–6.66%), respectively. The overall density was 1.34 g/mL (0.91–1.53 g/mL), the average porosity was 0.44% (0.18–0.63%), and the average water content was 25.87% (8.89–42.69%). The density tended to decrease as the average particle size decreased, whereas the moisture content and porosity showed the opposite trend. The TOC/TN ratio was consistently ≤ 10, which implies the influence of marine rather than terrestrial organisms. The average porewater salinity was 39.52 psu (range: 11.00–57.13 psu). The organic matter was characterized based on an analysis of particle size and the results confirmed the experimental values, the smaller the average particle size, the higher the concentration of organic matter.
The average elevation values based on data from the 107 stations was –0.66 m (–2.71 to +1.70) as shown in
Figure 4. The relationships between the sedimentary environmental factors and the topography of the intertidal zone were investigated using the Hwangdo tidal-flat intertidal DEM created by the GCP-drone orthogonal remote sensing method. Although this method is far less accurate than a leveling of sight lines, it is effective in assessing the overall topography of an intertidal zone [
36]. The topographical distribution was similar to that indicated by the orthometric height, with slopes rising slightly to the west (
Figure 4). The overall trend showed a topography with a high center, with peaks slightly skewed to the west and steep slopes in the west, while the eastern slope was relatively gentle. The two sandbars located to the east had very high elevation. In the western area, at around 0 m elevation, the features of the drone image were those of highly developed tidal channels in the western region. The overall trend showed a topography with a flat center, slightly tilted to the west, a gentle slope to the east and a sharp slope towards the western side in the morphology and a tidal channel. The two sandbars located in the east were of very high altitude. In the area ~2 m above sea level, the tidal channels were well-developed according to the DEM data, while the altitude of the rest of the area was relatively low. In terms of particle-size composition (mud ratio 36.7–95.2%) was predominant in the western region in the north-south direction, corresponding to the ridge of the upper intertidal zone of the Hwangdo tidal flat, and sand in the two sandbar sedimentary facies in the east (sand component 75.7–96.6%). This distribution reflects the development status of the main and minor channels, the presence or absence of surface residual water, and the topography [
37,
38,
39].
The above results show that the distribution of the tidal currents in the west of the Hwangdo tidal flat is closely related to topographical features. For reference in Ryu 2005, high-resolution images (IKONOS RGB (432 bands)) together with the DEM can be used to classify three types of sedimentary facies: mud, mixed, and sand sedimentary, considering the sedimentological characteristics of the intertidal zone, including the surface-layer size, surface residual water, tidal-channel distribution, and topography [
40].
Analysis of DEM profile lines showed that the changes in the elevation values of each of the eight latitudes decreased gradually as latitude decreased (H1 to H8 lines), and the overall topography distribution ranged from high (H1 line average -0.52 m) to low (H8 line average -1.55 m) in
Figure 5 (a) and
Figure 6 (a). The central region did not show a large change with latitude. However, as expected, the two sandbars located in the east had relatively high elevation. The 4.3-m elevation easting to 478 m on the H1 middle line was the topographical portion of the land. At the H4 line easting at 1,832 m, the elevation of 1.4 m represented a relatively high rock formation at the sandbar. The DEM-profile-line analysis showed that the elevation decreased rapidly with decreasing latitude. In addition, the elevation in the north increased irregularly from west to east, whereas the elevation in the south decreased, as shown in
Figure 5 (b) and
Figure 6 (b).
The area-change value according to the tidal-level change of the tidal flat was examined as well. The area at -4.0 m elevation was the maximum area range captured by the drone. At –4.0-m elevation, the exposed area was 4.841 km
2. With an increasing tidal level, only the high-elevation topography of the western part of the Hwangdo tidal flat remained at the 0.0 m mark, corresponding to an area of 0.800 km
2. The area with the largest increase in the rate of exposed area (Δ, %) at a tidal interval of 1 m or 0.5 m occurred at an elevation change from +2.0 m to 1.0 m, with an increase in the exposure area of 323%. Thereafter, the rate of increase gradually decreased as the altitude decreased. The tidal-flat exposure time is shown in
Figure 7 (a)-(l) and
Table 1 (25–28 Oct. 2022). The highest exposure-time ratio in the study area was within ~12 h, with an area of 2.47 km
2, and a rate of 50.8% on 26 Oct. in spring tides (spring range). Similarly, according to the 2022 KHOA, the highest exposure-time rate was within ~12 h, with an area of 2.48 km
2 and a rate of 51.1%. A calculation using the actual measurement data from the Boryeong’s sea level measurement tide data for the year 2022 showed that the highest exposure-time rate was within ~12 h, with an area of 1.43 km
2 and a rate of 29.5%. The drastic reduction during the period 25 to 28 Oct. 2022 compared to the control period can be explained by the tide level for the year 2022 according to Korea Hydrographic and Oceanographic Agency (KHOA), in which the greatest changes occurred twice a month on average throughout the year, when corrected by the Boryeong tide gauge for 2022 in
Figure 8 (a)-(c) and
Table 2.
4. Discussion and Conclusions
The spatial distribution of the Hwangdo tidal-flat stations was primarily determined by variations in sediment composition. Based on the cluster analysis and MDS the stations located at the tidal flat were divided into four groups (A, B, C, and D), representing upper, central, and lower stations, and regions of high organic-matter concentration, respectively (
Figure 9 and
Figure 10). In the dimension one (D1) in
Figure 10, the distance of two points present on the real line is the absolute value of the arithmetic difference of the coordinates (D1(p, q) = |p−q|). The major environmental variable in group A (upper stations) was the sand composition; in group B (central stations), mean particle size was silt; in group C (lower stations) it was the sand composition; and in group D it was the high organic carbon concentrations.
The degree to which the combinations of the various factors affected the formation of the tidal-flat groups in the ecliptic plane was determined. The closer the value to 1, the stronger the correlation (
Table 3). Among the analyzed factors, the single most influential was the average particle size. Combining the factors increased their degree of influence on group formation, but the differences were small and complicated explanations of the effect of each factor. The main factors influencing the formation of the Hwangdo tidal-flat group were sand, silt, and clay contents, which also collectively determined the average particle size.
The correlation between the exposure time and the spatial distribution of the sediment composition was low (
Figure 11). The tides near the tidal channel to the west of Hwangdo are of relatively low energy, which affects the tidal flats, with increasing sedimentation in the central area. At Cheonsu Bay, man-made structures such as seawalls have altered seawater circulation. The resulting changes in sedimentation patterns and flow have contributed to changes in the current and thus in the sedimentary phase type of the surface sediments of Hwangdo.
Figure 12 and
Table 4 shows the grouping by upper, middle, lower intertidal zones and the area of high organic-matter content.
An analysis of the oyster-reef area showed an increase in 2019 compared to 2015, based on the linear spectral unmixing classifications obtained with the KOMPSAT-2/3 data [
41]. By contrast at the H40, H27, and H28 stations, high tidal-channel-index values indicative of a dendritic pattern of tidal-channel distribution were identified in areas of high elevation and fine-grained sediments, as determined using remotely sensed data such as those from the KOMPSAT-2, Landsat Enhanced Thematic Mapper Plus (ETM+), and aerial photographs [
42]. A previous study [
43] confirmed our choice of a statistical approach to intertidal classification and topography to monitor the near-real-time spatiotemporal distribution changes in the tidal flats using continuous and stable Synthetic Aperture Radar (SAR) data obtained at local and regional scales. The dataset generated in this study can be used for field verification.
In October 2022, data on sedimentary-environment factors, such as geochemical data and topographical-elevation data, were obtained using drones, and data on the surface sedimentary facies vertex were obtained by analyzing the surface sediments. Thirteen types of sedimentary facies were identified. The distribution of surface sediments was generally fine-grained sandy mud sediments in the west, muddy sand sediments in the central region, and coarse gravelly sand sediments in the east. In the study by Kim et al. [
44], data on sedimentary environmental factors, such as surface sedimentary facies and elevation between 2004 and 2013, were obtained in field surveys, and the soil moisture content of each sediment was measured to analyze the correlation between seawater and optical reflectance in the tidal flat. The results showed 13 sedimentary facies distributed in the Hwangdo tidal flat, with a gradual increase in the sand content and elevation between 2004 and 2013 by Kim et al. [
44].
In a previous study [
40], the classification of sedimentary facies using IKONOS RGB (432 bands) images collected over the Hwangdo tidal flat in Cheonsu Bay was qualitatively investigated. The optical reflectance was compared with environmental and sedimentary factors such as grain size, tidal-channel patterns, and the presence of surface remnant water. The results showed that satellite images of high spatial resolution, such as those from IKONOS, can be used to classify surface sedimentary facies (mud, mixed and sand) if tidal sedimentary characteristics, including grain size, surface remnant water, and the tidal channel network, are considered as well.
In a study by Kim et al. [
45], red, green, and blue (RGB) orthoimagery from an unmanned aerial vehicle was used in combination with a field survey (232 samples) to produce a large-scale classification map for surface-sediment distribution, in accordance with sedimentology guide. The object-based method employed in that study was used to classify surface-sediment distribution based on correlations with spectral reflectance, grain size, and tidal channels. Six sediment types were distinguished according to their spectral reflectance and sediment properties, such as grain composition, and statistical parameters. In the present study, the classification was based on the four groups of Folk and Ward [
27], and the mapping results were similar to the four statistically processed groups.
A comparative evaluation of this study with previous studies showed that the seawater circulation in the northern part of Cheonsu Bay was changed by the embankment, thus affecting the sediment circulation. Changes in sedimentary patterns, rather than in the exposure time, accounted for the differences in the sediment compositions of the upper, middle, and lower stations, a response that is expected to continue. High organic carbon or benthic organisms will likewise be affected. Minimizing the impact on the ecosystem of the Hwangdo tidal flats in Cheonsu Bay caused by these changes, and preservation of the tidal flats and surrounding waters, will require implementation of a management plan for the sedimentary environment, one that takes into account the survey results of continuous monitoring. The results will contribute to improved management and conservation of tidal flats. It is judged that future research should focus on developing evaluation methods and procedures that will become key elements of marine spatial planning (MSP). In addition, it is necessary to develop through research such as evaluation target model through verification, evaluation elements, spatialization of each element, and verification of evaluation method. The MSP-based decision-making system should therefore be established. Several steps must be completed to achieve this, as well as to clarify issues associated with the Hwangdo tidal flat. With the ultimate aim of applying MSP throughout Korea, realizing a successful MSP system for the extensive and difficult to access Hwangdo tidal flat would be a useful starting point. The use of a spatial coordinate system, such as GPS, is a key element in MSP, and could be more effective than the sector-based management approach that predominated previously. Continuous efforts are needed to further advance and develop national MSP policies that take account of social-economic conditions, and the status of marine technologies.
In future studies, a dataset comprising geological environment characteristics, obtained based on the approach used in this study, will be used to assess changes in tidal-flat topography and the sedimentation environment.
Author Contributions
Conceptualization, J.-h.L. and H.J.W.; methodology, J.-h.L.; validation, J.-h.L, H.J.W. and H.-s.J.; formal analysis, J.-h.L, H.J.W. and H.-s.J.; investigation, J.-h.L, Y,J., J.B.J., J.S. and K.K.; data curation, J.-h.L.. Y. J., J.S. and K.K.; writing— original draft preparation, J.-h.L.; writing—review and editing, J.-h.L., H.J.W. and H.-s.J.; visualization, J.-h.L. and Y.J.; supervision, J.-h.L., H.J.W. and J.-h.R.; project administration, J.-h.L., K. K. and J.-h.R.; funding acquisition, J.-h.R. All authors have read and agreed to the published version of the manuscript.
Figure 1.
Map of the study area. (a) Location of the study area within South Korea. (b) Detailed map of the study area, showing Hwangdo tidal flats, Taean, Cheonsu Bay, west coast of Korea. (c) Red rhombi indicate the surface-sediment sampling stations (n = 107) and (d) Enlarged map of the sampling stations (n = 86) showing saltmarsh and tidal flat, tidal channels, small chenier, and sand shoals.
Figure 1.
Map of the study area. (a) Location of the study area within South Korea. (b) Detailed map of the study area, showing Hwangdo tidal flats, Taean, Cheonsu Bay, west coast of Korea. (c) Red rhombi indicate the surface-sediment sampling stations (n = 107) and (d) Enlarged map of the sampling stations (n = 86) showing saltmarsh and tidal flat, tidal channels, small chenier, and sand shoals.
Figure 2.
Field survey using drones for digital elevation model (DEM) processing. (a) Drone survey using a full-frame sensor (Matrice 300 drone Real-Time Kinematic-Global Positioning System (RTK-GPS) + GLONASS + BeiDou + Galileo; DJI). (b) RTK-GPS station for measuring tidal-flat elevations using a rover. (c) Orthophoto taken route and at a ground control point (GCP) station (n = 17) using a Matrice 300 drone. The location of Hwangdo Bridge, where the tide gauges were installed during the survey period, is also shown.
Figure 2.
Field survey using drones for digital elevation model (DEM) processing. (a) Drone survey using a full-frame sensor (Matrice 300 drone Real-Time Kinematic-Global Positioning System (RTK-GPS) + GLONASS + BeiDou + Galileo; DJI). (b) RTK-GPS station for measuring tidal-flat elevations using a rover. (c) Orthophoto taken route and at a ground control point (GCP) station (n = 17) using a Matrice 300 drone. The location of Hwangdo Bridge, where the tide gauges were installed during the survey period, is also shown.
Figure 3.
Hwangdo Bridge, the site of the installed tide gauges, and graphs of the tide-level changes during the survey period. (a) Actual tidal change at Hwangdo Bridge from October 25 to October 28, 2022. (b) The seawall constructed in 2012, before the construction of Hwangdo Bridge [
25]. (c) High tide on October 15, 2022 during spring tide. (d) Low tide on October 27, 2022 during spring tide (spring range).
Figure 3.
Hwangdo Bridge, the site of the installed tide gauges, and graphs of the tide-level changes during the survey period. (a) Actual tidal change at Hwangdo Bridge from October 25 to October 28, 2022. (b) The seawall constructed in 2012, before the construction of Hwangdo Bridge [
25]. (c) High tide on October 15, 2022 during spring tide. (d) Low tide on October 27, 2022 during spring tide (spring range).
Figure 4.
3D plot after DEM processing using the orthophoto taken route using a Matrice 300 drone.
Figure 4.
3D plot after DEM processing using the orthophoto taken route using a Matrice 300 drone.
Figure 5.
Grid line for elevation analysis. (a) Lines used in the analysis of elevation changes by latitude using DEM data. (b) Lines used in the analysis of elevation changes by longitude using DEM data in
Figure 4.
Figure 5.
Grid line for elevation analysis. (a) Lines used in the analysis of elevation changes by latitude using DEM data. (b) Lines used in the analysis of elevation changes by longitude using DEM data in
Figure 4.
Figure 6.
Elevation changes by longitude and latitude. (a) Altitude changes by latitude from H1 to H8. (b) Altitude changes by longitude from V1 to V5.
Figure 6.
Elevation changes by longitude and latitude. (a) Altitude changes by latitude from H1 to H8. (b) Altitude changes by longitude from V1 to V5.
Figure 7.
Changes in the exposed area of the tidal flats in the study area based on changes in the tide level (DEM –4.0 to 4.0 m).
Figure 7.
Changes in the exposed area of the tidal flats in the study area based on changes in the tide level (DEM –4.0 to 4.0 m).
Figure 8.
Comparison of exposure times in 2022 according to the research period processing method. Exposure-time images obtained (a) on October 26, 2022 in spring tides (spring range), (b) between October 25 to 28, 2022 in spring tides (spring range), and (c) for the year 2022 according to Korea Hydrographic and Oceanographic Agency (KHOA).
Figure 8.
Comparison of exposure times in 2022 according to the research period processing method. Exposure-time images obtained (a) on October 26, 2022 in spring tides (spring range), (b) between October 25 to 28, 2022 in spring tides (spring range), and (c) for the year 2022 according to Korea Hydrographic and Oceanographic Agency (KHOA).
Figure 9.
Cluster analysis of the square-root-transformed environmental variables. However, in the cluster analysis analysis, some samples did not belong to the correct group and were moved to other groups according to the results of multidimensional scaling (MDS) and principal components analysis (PCAs) (H25-1, wh-4, H29, OC3, H67, H16 -1, and A16).
Figure 9.
Cluster analysis of the square-root-transformed environmental variables. However, in the cluster analysis analysis, some samples did not belong to the correct group and were moved to other groups according to the results of multidimensional scaling (MDS) and principal components analysis (PCAs) (H25-1, wh-4, H29, OC3, H67, H16 -1, and A16).
Figure 10.
Multidimensional scaling (MDS) based on the D1 Euclidean distance using the square-root-transformed data.
Figure 10.
Multidimensional scaling (MDS) based on the D1 Euclidean distance using the square-root-transformed data.
Figure 11.
Principal components analysis (PCAs) of environmental variable factors.
Figure 11.
Principal components analysis (PCAs) of environmental variable factors.
Figure 12.
Grouping by upper, middle, and lower intertidal zone and the area of high organic matter. (a) Map of the station sites, (b) Orthogonal scale bar (50 × 50 cm) and landscape photos of the area in group A, (c) Shows a colour coded distribution of all the group samples on a map of the tidal flat, (d) group B, (e) group D, (f) group D, and (g) The seawater-logged road used to access the tidal flats.
Figure 12.
Grouping by upper, middle, and lower intertidal zone and the area of high organic matter. (a) Map of the station sites, (b) Orthogonal scale bar (50 × 50 cm) and landscape photos of the area in group A, (c) Shows a colour coded distribution of all the group samples on a map of the tidal flat, (d) group B, (e) group D, (f) group D, and (g) The seawater-logged road used to access the tidal flats.
Table 1.
Changes in the exposed areas (m2, km2, and rate of exposed area Δ%) of the tidal-flat in the study area based on changes in the tide level (DEM –4.0 to 4.0 m).
Table 1.
Changes in the exposed areas (m2, km2, and rate of exposed area Δ%) of the tidal-flat in the study area based on changes in the tide level (DEM –4.0 to 4.0 m).
Tidal level (m) |
Exposed area (m2) |
Exposed area (km2) |
Rate of exposed area (Δ, %) |
Tidal level (m) |
Exposed area (m2) |
Exposed area (km2) |
Rate of exposed area (Δ, %) |
+ 4.0 |
6,362 |
0.006 |
- |
-0.5 |
1,728,087 |
1.728 |
116 |
+ 3.0 |
8,274 |
0.008 |
30 |
-1.0 |
2,740,140 |
2.740 |
59 |
+ 2.0 |
13,494 |
0.013 |
63 |
-1.5 |
3,561,450 |
3.561 |
30 |
+ 1.0 |
57,126 |
0.057 |
323 |
-2,0 |
4,149,875 |
4.150 |
17 |
+ 0.5 |
215,435 |
0.215 |
277 |
-3.0 |
4,640,762 |
4.641 |
12 |
0.0 |
799,813 |
0.800 |
271 |
-4.0 |
4,841,078 |
4.841 |
4 |
Table 2.
Comparison of exposure times in 2022 according to the research period processing method (1 day, 3 days, and the year 2022 according to Korea Hydrographic and Oceanographic Agency (KHOA)).
Table 2.
Comparison of exposure times in 2022 according to the research period processing method (1 day, 3 days, and the year 2022 according to Korea Hydrographic and Oceanographic Agency (KHOA)).
Exposure time (h) (Guideline: 1 Jau. 2022-31 Dec. 2022, 1year) |
Exposed area (m2, km2) |
Exposed area (%) |
Exposure time (h) (Guideline: 25-28 Oct. 2022, 3 day) |
Exposed area (m2, km2) |
Exposed area (%) |
Exposure time (h) (Guideline: 26 Oct. 2022, 1 day) |
Exposed area (m2, km2) |
Exposed area (%) |
≤ 3 h |
557,753 0.56 |
11.5 |
≤ 3 h |
168,772 0.17 |
3.5 |
≤ 3 h |
151,850 0.15 |
3.1 |
~ 6 h |
1,251,492 1.25 |
25.7 |
~ 6 h |
289,550 0.29 |
6.0 |
~ 6 h |
347,895 0.35 |
7.2 |
~ 9 h |
1,048,798 1.05 |
21.6 |
~ 9 h |
1,153,271 1.15 |
23.7 |
~ 9 h |
1,226,007 1.23 |
25.2 |
~ 12 h |
1,433,997 1.43 |
29.5 |
~ 12 h |
2,467,386 2.47 |
50.8 |
~ 12 h |
2,484,270 2.48 |
51.1 |
~ 15 h |
500,274 0.50 |
10.3 |
~ 15 h |
755,080 0.76 |
15.5 |
~ 15 h |
628,871 0.63 |
12.9 |
~ 18 h |
52,503 0.05 |
1.1 |
~ 18 h |
16,386 0.02 |
0.3 |
~ 18 h |
12,224 0.01 |
0.3 |
~ 21 h |
6300 0.01 |
0.1 |
~ 21 h |
2632 0.00 |
0.1 |
~ 21 h |
2168 0.00 |
0.0 |
~ 24 hour |
4025 0.00 |
0.1 |
~ 24 hour |
1371 0.00 |
0.0 |
~ 24 hour |
1137 0.00 |
0.0 |
~ 24 hour < (at all times) |
6175 0.01 |
0.1 |
~ 24 hour < (at all times) |
6869 0.01 |
0.1 |
~ 24 hour < (at all times) |
,895 0.01 |
0.1 |
Sum |
4,861,317 |
100 |
Sum |
4,861,317 |
100 |
Sum |
4,861,317 |
100 |
Table 3.
Best result for each combination of variables among the textural parameters of the sediments based on Spearman rank (closer to 1 indicates higher correlation).
Table 3.
Best result for each combination of variables among the textural parameters of the sediments based on Spearman rank (closer to 1 indicates higher correlation).
Rank |
Factor and value |
1 |
Mean size (Mz) 0.780 |
2 |
Mean size (Mz), Sorting 0.855 |
3 |
Mean size (Mz), Sorting, TIC 0.888 |
4 |
Mean size (Mz), Sorting, TIC, Skewness 0.889 |
5 |
Mean size (Mz), Sorting, TIC, Skewness, Dry bulk density 0.891 |
Table 4.
Summary of sedimentary environments at re-presentive sampling sites.
Table 4.
Summary of sedimentary environments at re-presentive sampling sites.
Group |
Station No. |
Sedimentary facies type (Folk 1957) |
Gravel (%) |
Sand (%) |
Silt (%) |
Clay (%) |
Tidal channel pattern |
Elevation /Slope |
A |
A24 |
Mud flat gM (gravelly Mud) |
6.58 |
45.06 |
24.82 |
23.54 |
complex and meandering with oyster reefs |
Lower part /steep |
H27 |
Mud flat sM (sandy Mud) |
0.00 |
14.42 |
46.16 |
39.42 |
complex and meandering with oyster reefs |
Upper part /steep |
H28 |
Mud flat sZ (sandy Silt) |
0.00 |
26.15 |
49.76 |
24.09 |
complex and meandering with oyster reefs |
Upper part /steep |
H40 |
Mud flat M (Mud) |
0.00 |
15.9 |
49.33 |
34.77 |
complex and meandering |
Upper part /steep |
B |
A18 |
Mixed flat mS (muddy Sand) |
0.00 |
62.59 |
22.67 |
14.74 |
simple and straight |
Middle part /gentle |
H32 |
Mixed flat zS (silty Sand) |
0.00 |
54.87 |
31.36 |
13.77 |
simple and straight |
Middle part /gentle |
H57 |
Mixed flat mS (muddy Sand) |
0.00 |
78.59 |
11.54 |
9.87 |
simple and straight |
Middle part /gentle |
OC14 |
Sand flat gS (gravelly Sand) |
23.82 |
75.73 |
0.15 |
0.30 |
not present |
Upper part /steep |
C |
H25-1 |
Sand flat msG (muddy sandy gravel) |
30.70 |
61.39 |
2.89 |
5.02 |
simple and straight |
Upper part /steep |
H64 |
Sand flat (g)S (slightly gravelly Sand) |
0.34 |
89.82 |
4.58 |
5.26 |
not present |
Lower part /gentle |
D |
H23-1 |
Sand flat gmS (gravelly muddy Sand) |
14.31 |
57.24 |
14.49 |
13.96 |
simple and straight with shell framents |
Lower part /gentle |