Magnetic Properties of a Holocene Sediment Core from the Yeongsan Estuary, Southwest Korea: Implications for Diagenetic Effects and Availability as Paleoenvironmental Proxies

The sensitivity of magnetic properties, which characterize the mineralogy, concentration, and grain size distribution of magnetic minerals, to environmental processes may provide useful information on paleoenvironmental changes in estuarine environments. Magnetic property studies of estuaries are less common than other environments and, due to the west coast of South Korea having an abundance of estuaries, it provides a good place to study these processes. In this study, we analyzed a variety of magnetic properties based on magnetic susceptibility, hysteresis parameters, progressive acquisition of isothermal remanent magnetization and first-order reversal curve data from a Holocene muddy sediment core recovered from the Yeongsan Estuary on the west coast of South Korea. We examined diagenetic effects on magnetic properties and tested their availability as proxies of paleoenvironmental change. The presence of generally low magnetic susceptibility, ubiquitous greigite-like authigenic magnetic component, and very fine magnetic particle occurrence suggested that the analyzed sediments had undergone considerable early diagenetic alteration. Electron microscopic observations of magnetic minerals support this suggestion. Our results confirm that the use of initial bulk susceptibility as a stand-alone environmental change proxy is not recommended unless it is supported by additional magnetic analyses. We recognized the existence of ferromagnetic-based variabilities related to something besides the adverse diagenetic effects, and have examined possible relationships with sea-level and major climate changes during the Holocene. The most remarkable finding of this study is the two distinct intervals with high values in magnetic coercivity (Bc), coercivity of remanence (Bcr), and ratio of remanent saturation moment to saturation moment (Mrs/Ms) that were well coincident with the respective abrupt decelerations in the rate of sea-level rise occurred at around 8.2 and 7 thousand years ago. It is then inferred that such condition with abrupt drop in sea-level rise rate would be favorable for the abrupt modification of grain size distribution toward more single-domain-like content. We modestly propose consideration of the Bc, Bcr, and Mrs/Ms variability as a potential indicator for the initiation/occurrence of sea-level stillstand/slowstand or highstand during the Holocence, at least at estuarine environments in and around the studied area.


INTRODUCTION
Magnetic minerals are common constituents of a wide range of sediments and are sensitive to the physicochemical conditions of their surrounding environment (Thompson and Oldfield, 1986;Verosub and Roberts, 1995;Dekkers, 1997;Maher and Thompson, 1999;Evans and Heller, 2003;Torii, 2005;Liu et al., 2012). Magnetic measurements facilitate characterization of the mineralogy, concentration, and grain size distribution of magnetic minerals contained in sediments, providing information about sediment provenance, transportation, and deposition, as well as post-depositional diagenesis and input of urban and industrial sources into sediments (Karlin and Levi, 1983;Bloemendal et al., 1992;Robinson et al., 2000;Emiroglu et al., 2004;Rey et al., 2005;Kim et al., 2009;Liu et al., 2010;Szuszkiewicz et al., 2015;Tauxe et al., 2015;Pan et al., 2017).
It is generally thought that sedimentation during Holocene at estuaries and coastal lines are sensitive to not only their catchment (terrestrial) environmental changes and extreme hydrologic events such as heavy rainfalls and storms, which are closely associated with regional and global climate changes, but also regional sea-level change (e.g., Lim et al., 2017, Lim et al., 2019. Hence, such estuarine-coastal sediments are of potential use for reconstructing records of paleoenvironmental changes during Holocene. Along the west coast of South Korea, facing the eastern Yellow Sea, estuaries of large rivers such as the Han, Geum, and Yeongsan Rivers are widely distributed from north to south; these are tide-dominated depositional environments (Chough et al., 2004) forming a variety of deposits during the Late Pleistocene and Holocene (Park et al., 1998;Choi et al., 2003;Lim and Park, 2003;Nahm et al., 2008;Nahm and Hong, 2014;Moon et al., 2018). During the Holocene, these western coastal areas would have experienced dramatic environmental changes due to e.g., sea-level changes such as the Holocene transgression (Stanley and Warne, 1994;Kim and Kennett, 1998), the El Niño Southern Oscillation (ENSO) activity (e.g., Lim et al., 2017;Lu et al., 2018;Lim et al., 2019), the eastern Asian Monsoon (EAM; e.g., Chang, 2004;Dykoski et al., 2005;Selvaraj et al., 2007), and the Holocene climate optimum (HCO; e.g., An et al., 2000;Zhou et al., 2016;Park et al., 2019).
Bulk low-field magnetic susceptibility has been widely applied to interpret such sediment environments in South Korea (Park et al., 1998;Moon et al., 2018;Lim et al., 2004;Lim et al., 2014;Lim et al., 2015), yet comprehensive studies using other magnetic properties have been rare. Further studies are needed to comprehensively explore magnetic minerals in terms of their supply into the accommodation space, dissolution or alteration or generation by early diagenesis, abundance, and response to environmental changes. Additionally, characterization of such magnetic properties among fine materials (muds) in the estuaries of the western Korean Peninsula may provide fundamental information on the provenance and hydrodynamic transport of mud deposits in the Yellow Sea and East China Sea, which remain poorly constrained (Yang et al., 2003;Liu et al., 2010;Wang et al., 2010;Koo et al., 2018).
The objective of this study was to characterize variations in magnetic properties from a sediment core retrieved from an estuarine system in the Yeongsan Estuary on the west coast of South Korea and to identify diagenetic effects that can cause down-core variation of magnetic properties. We then tested the availability of such magnetic properties as proxies of paleoenvironmental changes during the Holocene.

STUDY AREA, MATERIALS, AND RADIOCARBON AGES
The Yeongsan Estuary is located within the city of Mokpo (population, ∼250,000) on the southwestern coast of South Korea; from land to sea, it is divided into Yeongsan Lake (fresh water lake), the inner and outer estuaries, and the coastal zone ( Figure 1; Williams et al., 2014). The estuary is under a macrotidal regime, with primarily semidiurnal tides of about 4.5 m tidal range (Byun et al., 2004). Across the estuary, the Yeongsan River (total length, 115 km; drainage basin area, 3,371 km 2 ) flows into the Yellow Sea through the relatively shallow (mostly <50 m water depth) ria coastline. The average depth and width of the Yeongsan River ranges from ∼10 to 19 m and from 0.6 to 1.3 km, respectively, in its lower reaches (Lee et al., 2009;Williams et al., 2014). The climate of the Korean Peninsula, including the study area, is dominated by seasonal monsoons: during winter, cold and dry north-northwesterly winds are accompanied by less precipitation; during summer, warmer and wetter south-southeasterly winds are accompanied by heavy precipitation and occasional typhoons (e.g., Chang, 2004). Due to the seasonality in precipitation, freshwater discharge of the Yeongsan River occurs mainly during summer (>80% of the annual mean; Ryu et al., 2004). The geology in and around the Yeongsan River basin is composed mainly of Precambrian gneiss, Paleozoic sedimentary rocks such as shale and mudstone, and Mesozoic granites, volcanic extrusive rocks, and tuffs (Choi et al., 2002).
The Yeongsan Estuary has experienced substantial coastal construction within the last 100 years, a dam construction (the Yeongsan Estuarine Dam, drawing boundary between the Inner Estuary and the Yeongsan Lake at present) by 1981, and land reclamation projects during the 1980s, which led to reduction of the estuarine area (e.g., Williams et al., 2014). Prior to the dam construction, tidally influenced environments spanned about 63 km upstream from the dam (Lee et al., 2009). During the land reclamation projects, seawalls/embankments were constructed within the area.
YAR-4 (34.82167°N, 126.55357°E; 0.5 m below mean sea level at the top) is a 20-m-long sediment core (inner diameter, ∼74 mm) that was recovered adjacent to the junction of the Yeongam River (also called the Yeongam tributary) and the main channel of the Yeongsan River ( Figure 1). Prior to the dam construction the core site would have been below mean sea level and the dam construction actually exposed the core location. The YAR-4 core site has had no tidal influence since the dam construction. The average accumulation rate of sediments around the site is presently estimated at ∼20 mm/yr (Williams et al., 2014). Weathered granite (basement rock) is apparent at the bottom 0.20 m of this core; from a core depth of 19.80 m to the top, various facies of sediments are visible ( Figure 2). The sediment descriptions, interpretation of the sedimentary environments, and accelerator-mass spectrometer (AMS) radiocarbon ( 14 C) age data for this core have been previously documented (Nakanishi et al., 2013), and are briefly described as follows. Nakanishi et al. (2013) divided the YAR-4 core sediment stratigraphy into six different sedimentary facies units (I, II, III, IV, V, and VI, in Figure 2) based on their analyses, including grain size, color, sedimentary features, and AMS 14 C age data. Units II and III are further subdivided into two subfacies units each, as IIa, IIb, IIIa, and IIIb. Figure 2 presents a summary of the core's sedimentary characteristics and the interpretation of associated sedimentary environments. Nakanishi et al. (2013) also estimated a number of AMS 14 C ages (i.e., only conventional age estimates, without 14 C calibration) from materials of terrestrial (plants) and marine (shells) origins to evaluate the marine reservoir effects on the sediments. Their results revealed that the shell 14 C ages were older by 0-500 years (i.e., the marine reservoir effects) than the plant 14 C ages.
In the present study, we selected part of the plant 14 C age data from the published age dataset to represent appropriate sediment ages (Table 1; Figure 3) to build up an age-depth relationship for this core and to estimate the sedimentation rates. These ages ranged from 8,130 (±70, 1σ) to 3,760 (±50, 1σ) yr BP in conventional age for units II to V, being converted into calibrated ages from 9,088 (±310, 2σ) to 4,114 (±179, 2σ) cal. yr BP by using the OxCal (version 4.4) calibration program (https://c14.arch.ox.ac.uk/oxcal.html) with the IntCal20 (Reimer et al., 2020). Unit I had no 14 C age determination due to the absence of samples suitable for 14 C dating. Based on these selected data points, we constructed an age-depth model ( Figure 3) by using the software "Undatable" (Lougheed and Obrochta, 2019), which ran 100,000 Monte Carlo iterations with setting parameters of xfactor 0.1 and bootpc 10 in treating uncertainty (for details about the xfactor, the bootpc, and the modeling, see Lougheed and Obrochta, 2019). This YAR-4 age-depth model indicated several drastic changes in the sedimentation rate, which ranged apparently between 0.4 and >25 mm/yr.

METHODS
All of the following subsample preparations and procedures for analyses of sediment particle size, total organic carbon (TOC), and a suite of magnetic properties were conducted at the laboratories of the Korea Institute of Geoscience and Mineral Resources (KIGAM, Republic of Korea). Only microscopic observations were conducted at the Central Research Facilities of Gyeongsang National University (Republic of Korea).
For the split half-core, subsamples for analyses were collected using cylindrical plastic tubes generally at ∼0.04-m intervals and a subsampling thickness of ∼0.02 m from the core depth interval of 3.20-20.00 m. The subsamples were oven-dried at 40°C for approximately 24 h prior to analyses.

Particle Size and Total Organic Carbon Measurements
The particle size distribution of bulk sediment was analyzed at ∼0.5-m intervals using a Mastersizer 2000 laser analyzer (Malvern Instrument Ltd., United Kingdom) using approximately 300 mg of each dried subsample, after treatment with 35% H 2 O 2 and 1 N HCl to dissolve organic matter and biogenic carbonates, respectively, and then with ultrasonic dispersion to facilitate complete disaggregation.
To measure the TOC content of bulk sediment, each subsample was treated with 1 N HCl at approximately 100°C for 1 h and then transferred to a tin combustion cup after rinsing with distilled water. The TOC content was analyzed at ∼0.2-m intervals using a CNS elemental analyzer (vario Micro cube; Elementar, Germany).

Measurements of Magnetic Properties
Magnetic properties were measured using three types of dried subsamples: bulk sediment sealed in a non-magnetic plastic box (volume, 7 cm 3 ), several tens of mg of bulk sediment sealed in a gelatin capsule and capped with glass wool, and several tens of mg of magnetic mineral extracts sealed in a gelatin capsule and capped with glass wool. Each extraction of magnetic minerals was obtained using a hand magnet. The numbers of each type of the subsample prepared were 153, 41, and 31, respectively. The "packed-in-plastic box" bulk subsamples were used in measurements with a MS2B magnetic susceptibility meter (Bartington Instruments Ltd., United Kingdom). Both of the bulk and extracts "packed-in-gelatin capsule" subsamples were used in measurements with a MicroMag Model 3,900 vibrating sample magnetometer (VSM; Princeton Measurements Corp., United States). By using the MS2B magnetic susceptibility meter, volume-specific magnetic susceptibilities with low-frequency (0.47 kHz) and highfrequency (4.7 kHz) (k LFQ and k HFQ , respectively) were obtained. By using the VSM, magnetic hysteresis loop measurement, progressive alternating field (AF) demagnetization of isothermal remanent magnetization (IRM), progressive IRM acquisition, and first-order reversal curves (FORCs) measurement were conducted. The MS2B measurements were applied on subsamples at about 0.04-0.10-m intervals, and the VSM measurements were applied on subsamples that were from a limited number of different horizons of the sediment core.
The k LFQ and k HFQ of subsample were established by the average of double measurements with the correction for diamagnetic contribution of the plastic box to magnetic susceptibility, which was made by subtracting the average value of measurements of five empty boxes. These allowed us to calculate frequency-dependent susceptibility (k FD in %) using the following formula: k FD (k LFQ − k HFQ )/k LFQ × 100.
Prior to the VSM measurements, the "packed-in-gelatin capsule" subsamples to be measured were weighed to permit calculation of mass-specific values. Each hysteresis loop was measured up to maximum fields of ±0.5 T, and then, after correction for paramagnetic contribution within the loop, saturation magnetization (M s ), saturation remanence (M rs ), and magnetic coercivity (B c ) were calculated. From each hysteresis loop, low-field magnetic susceptibility (χ lf ; total magnetic susceptibility) was calculated from the initial slope and high-field magnetic susceptibility (χ hf ), as a susceptibility estimate of paramagnetic (plus diamagnetic) minerals, was calculated from the high-field slope of the hysteresis loop. Note that for magnetic susceptibilities, k denotes volume-specific 1 | Accelerator mass spectrometry (AMS) 14 C age data for materials of terrestrial origin (plant fragments) of the YAR-4 core, selected from the dataset of Nakanishi et al. (2013). Age* denotes a conventional age estimate before 14 C calibration reported in Nakanishi et al. (2013). Age (±2σ) in cal. yr BP denotes the calibrated age estimate using the OxCal (version 4.4; https://c14.arch.ox.ac.uk/oxcal.html) and the IntCal20 (Reimer et al., 2020)  FIGURE 3 | Plot of the calibrated 14 C age data that are used for this study (also given in Table 1) with depth in core and elevation, and an age-depth model reconstructed by using the "Undatable" (Lougheed and Obrochta, 2019), with xfactor 0.1, bootpc 10, and 10 5 Monte Carlo iteration runs (see Supplementary Table S6). Each age data point is given by red-filled circle. The age-depth model median is denoted by green line. The 1σ and 2σ uncertainties of the age-depth model are denoted by dashed blues and black lines, respectively. The sedimentary unit divisions (I to VI) are supplementarily shown near the y-axis.
Frontiers in Earth Science | www.frontiersin.org April 2021 | Volume 9 | Article 593332 susceptibility in SI units (dimensionless), and χ denotes massspecific susceptibility (m 3 /kg in SI units). From each progressive AF demagnetization of IRM, the remanence coercivity (B cr ) was determined as the applied field at which the remanence becomes zero. Each progressive IRM acquisition was made until 1 T was reached. Among the IRM acquisition data, IRM obtained at 1 T (IRM 1T ; roughly considered as SIRM) was used as a magnetic property in this study.
To unmix the magnetic mineral components contributing to the total remanent magnetization, each progressive IRM acquisition dataset was processed using the Max Unmix software (Maxbauer et al., 2016), which obtains the best fit for the total IRM acquisition using the fewest possible components that are characterized by the SIRM, the median field (Bh) at which half of the SIRM is reached, and dispersion (Dp) of its corresponding cumulative lognormal distribution. For each FORC analysis, 232 or 266 FORCs were measured; these FORC data were processed to create a FORC diagram using the FORCinel software (Harrison and Feinberg, 2008) with a smoothing factor of 9 or 10. Magnetic susceptibility and IRM 1T generally reflect the abundance of the total of magnetic mineral types and ferromagnetic mineral types (refer to ferrimagnetic and antiferromagnetic, in this study), respectively. χ hf /χ lf reflects the contribution of paramagnetic components to the total susceptibility. k FD indicates the contribution to total susceptibility made by viscous superparamagnetic (SP; < ∼30 nm) grains. B c and B cr may be used as measures of the relative proportions between low-and high-coercivity ferromagnetic mineral components, with higher values corresponding to higher proportions of high-coercivity minerals. B cr /B c , M rs /M s and FORC diagram provide information related to grain size distribution, mainly for ferromagnetic minerals.

Electron Microscope Observations
Electron microscopic observations were conducted selectively on bulk sediment subsamples from four horizons of 7. 84, 8.70, 13.40, and 15.25 m depth in core, in order to identify existing magnetic minerals, especially low-and intermediate-coercivity minerals. Polished surfaces of the bulk sediment subsamples were observed using a JSM-7610 F field emission scanning electron microscope (FE-SEM, JEOL, Japan) equipped with an energy-dispersive X-ray spectroscopy (EDS). From the observations, back scattered electron images and EDS spectra of magnetic minerals were obtained.

RESULTS AND DISCUSSION
Particle Size and Total Organic Carbon for Bulk Sediments Table 2 and Figure 4 present the results of particle size and TOC analyses of bulk sediment subsamples, and statistics including minimum, maximum, mean and standard deviation values for sedimentary units, and their downcore variation. The mean FIGURE 4 | Downcore variations in mean grain size (MZ), total organic carbon (TOC) content and various magnetic properties for bulk sediment samples of the YAR-4 core. The magnetic properties shown include bulk magnetic susceptibility (k LFQ , in both linear and log scale), isothermal remanent magnetization (IRM) intensity imparted at 1 T (IRM 1T ), high-field to low-field magnetic susceptibility ratio (χ hf /χ lf ), frequency-dependent magnetic susceptibility (k FD ), coercive force (B c ), coercivity of remanence (B cr ), ratio of B cr to B c (B cr /B c ), and ratio of saturation remanence to saturation magnetization (M rs /M s ). In B c , B cr , B cr /B c , and M rs /M s plots, exceptionally, the respective values for magnetic extracts were also shown (assuming their values for the extracts having no significant differences with those for bulk sediments). In the k FD -depth diagram, the hatched area denotes lack of meaningful k FD values due to weak signals. Sedimentary unit divisions (I to VI) are shown by orange lines. On the right edge of the figure, the depths at which significant change in sedimentation rate appeared based on the age-depth model ( Figure 3) are shown: red triangle denotes increase in sedimentation rate at above the depth, and blue inverse triangle denotes its decrease at above the depth. Gray-shaded area denotes lack of analyzed data.
Frontiers in Earth Science | www.frontiersin.org April 2021 | Volume 9 | Article 593332 particle sizes from depths of 3.2-∼14 m was fairly constant at ∼10 μm (fine silt); below a depth of 14 m, the mean size tended to increase with depth, with fluctuations reaching 104 μm (very fine sand) at a depth of 18.70 m. The TOC content was mainly very low (<0.9%); its variation was minor to a depth of ∼14 m, except at 12.50 m, and appeared to be negatively correlated with the mean size variation throughout the entire sediment core. Table 2 and Figure 4 summarize the results of major magnetic property analyses. Table 2 provides summary statistics of the individual magnetic properties including k LFQ , IRM 1T , χ hf /χ lf , k FD , B c , B cr , B cr /B c and M rs /M s , for the sedimentary units. Examples of the measured and high-field slope corrected hysteresis loops are illustrated in Figure 5. VSM-derived magnetic properties for bulk subsamples were barely detectable due to the weak and quite noisy signals during measurements, such that most of these bulk samples did not produce meaningful values in IRM 1T , χ hf /χ lf , B c , B cr , B cr /B c and M rs /M s . Instead, B c , B cr , B cr /B c and M rs /M s were obtained mostly from magnetic extracts of subsamples. We assume that the extract subsamples are representative of the bulk materials for these B c , B cr , B cr /B c , and M rs /M s values. Although there was exactly no horizon where bulk and extract subsamples both were measured for these magnetic property values in this study, it may be supported in part by the fact that, in the interval of 7.90-8.00 m depth in core, their values on bulk and extract subsamples were apparently similar. Downcore variations in these magnetic properties are shown in Figure 4. k LFQ , which was obtained using the MS2B meter and was volume-corrected for its proper quantification, ranged between 18 and 857 × 10 −6 SI for all analyzed samples, and the means of the respective sedimentary units were from 94 to 358 × 10 −6 SI. Most of the individual values are much lower than 400 × 10 −6 SI.

Downcore Variation in Magnetic Properties
Considerably high values with >400 × 10 −6 SI are limited in the intervals corresponding to the gravel-rich (around 16.9 and 18.15-19.80 m depth) and the basement rock (19.80-20.00 m depth) lithology. In the downcore k LFQ variation, there are a number of significant susceptibility minima representing lowering by at most 65-80% of the overall mean (for example, at depths of approximately 8. 70, 9.70, and 12.20-12.50 m). There are some datasets of k values from Holocene sediment cores at FIGURE 5 | Examples of measured (uncorrected) and high-field slope (i.e., paramagnetic contribution) corrected hysteresis loops. The differences between uncorrected and corrected loops indicate substantial non-ferromagnetic (i.e., paramagnetic) contribution in the studied sediments. Note that, for the extract subsample of 17.65 m depth (in core), the corrected loop seems to be less well-constrained so that the hysteresis parameters (M rs , M s , and B c ) could not be determinable, but the determined χ hf /χ lf seems meaningful.
Frontiers in Earth Science | www.frontiersin.org April 2021 | Volume 9 | Article 593332 around the studied region, making available for comparison with our data: in the western coastal estuaries of South Korea, ∼0.1-0.2 × 10 −6 SI (Gyeonggi Bay, Moon et al., 2018), and ∼2-12 × 10 −6 SI (Namyang Bay, Lim et al., 2004); in the southern coast of South Korea, ∼300-650 × 10 −6 SI (Geoje Island, Lim et al., 2014), and ∼50-400 × 10 −6 SI (Yeoja Bay, Lim et al., 2015); in the inner shelf of the southwestern East China Sea, ∼200-250 × 10 −6 SI (Zheng et al., 2010), and ∼100-220 × 10 −6 SI (Zheng et al., 2011). The extremely low k values in the Gyeonggi Bay and the Namyang Bay may result from locally particular layers (e.g., organic-rich layer, siderite-rich layer). Besides this, the ranges of k values of the previous studies are comparable to that for the YAR-4 core. However, this comparison should be taken care because it is not straightforward to identify whether or not these previously reported k values were after applying volume corrections (including calibrations of the prepared sample volume and the sensor of instruments used).
Values of mass-specific susceptibility (χ) can be more proper in comparison with previously reported datasets. Unfortunately, our analyzed subsamples were not weighed, thereby not permitting to determine accurately χ values for the same subsamples. Nevertheless, the χ values could be approximated by converting the aforementioned volume-specific susceptibility (k) values, using a mean density of 1.17 g/cm 3 and a mean water content of 71.3 wt% for wet surface sediments in the Yeongsan Lake . Assuming the water density of 1.00 g/cm 3 and that the porosity of the studied sediments (silt-dominant, shallow depths) was little varied through the whole core (i.e., negligible porosity variation for shallow depths; e.g., Bahr et al., 2001), a mean value in dry density of the studied sediment samples could be approximated to 2.03 g/cm 3 (considered as a maximal value). Accordingly, the individual k LFQ values could be converted to the mass-specific susceptibility (χ LFQ ) values in the range of ∼0.9-∼42 × 10 −8 m 3 /kg, and the χ LFQ values were generally in the order of 5-10 × 10 −8 m 3 /kg. This range of mass magnetic susceptibility values is similar in the order of magnitude of or slightly lower than those of Holocene sediment cores of the Yangtze delta area (e.g., ∼10-90 × 10 −8 m 3 /kg, Chen et al., 2015; ∼35-130 × 10 −8 m 3 /kg, Pan et al., 2017), shallow shelf sediment cores with anoxic environments in the Korea Strait and off the west coast of the Korean Peninsula (∼18 × 10 −8 m 3 /kg, Liu et al., 2004; ∼14 × 10 −8 m 3 /kg, Liu et al., 2005).
χ hf /χ lf values ranged between ∼0.6 and 1.0, indicating significantly larger paramagnetic contributions than ferromagnetic contributions (see also Figure 5), presumably in most of the sampled sediments.
The k FD values, with exception of negative values, ranged between 0.0 and 37.2%, of which values >20% may not be reliable due to weak susceptibility signals. The mean values for respective sedimentary units were approximately between 3 and 10%, indicative of, in general, somewhat or considerable SP contributions in the sediments. Due to concern from the mixing with untrustworthy data, we did not used the k FD data in discussing on its detailed variability in association with the diagenetic, sea-level, and climate effects.
Individual B c , B cr , B cr /B c , and M rs /M s values ranged between 5 and 21 mT, between 24 and 45 mT, between 2.44 and 5.20, and between 0.06 and 0.30, respectively. Downcore variation in B c and M rs /M s created an apparent division into two zones at a depth of ∼8 m, with relatively constant low values in the lower zone and constant higher values in the upper zone. Interestingly, downcore variation in B c , B cr , and M rs /M s shared a distinct, drastic increase in the respective parameter at depths of 15.25 m and 7.84-7.96 m. At depths above 7.84 m, such high values were sustained after this abrupt increase, whereas low values appeared the same as the surrounding data points at depths immediately above and below 15.25 m (i.e., 15.05 and 15.45 m). Figure 6 presents plots of k LFQ vs. the other magnetic properties examined in this study, to assess possible relationships with k LFQ . Generally, we observed no or very week correlation with k LFQ , but B cr had a moderate positive correlation with k LFQ . This implies complexity of mineral magnetic factors controlling k LFQ (e.g., content of paramagnetic minerals, content, composition ratio, and grain size of ferromagnetic minerals).

Magnetic Mineralogy and Granulometry from Magnetic Measurements and Scanning Electron Microscope Observations
SEM observations provided back scattered electron images and EDS spectra of magnetic minerals in the sediments (Figures 7 and Supplementary Figure S1). From the observations, pyrite (FeS 2 ), greigite (Fe 3 S 4 ), and (titano)magnetite are identified. It is found that pyrite is ubiquitous and has the most dominance in magnetic minerals in all analyzed bulk sediment samples (i.e., 7.84, 8.70, 13.40, and 15.25 m depth in core; Figures 7A and Supplementary Figure S1). The pyrite exists as framboidal aggregates and separated euhedral crystals with variable sizes in the sediment matrix, in the vicinity of silicates, and in the voids of the matrix (Figures 7A,C, and Supplementary Figure S1). The pyrite grains forming the framboids are much smaller (roughly <0.3 μm in maximal length) than the separated euhedral ones (roughly 0.5-<10 μm in maximal length). These pyrite grains all are considered to be of post-depositional authigenic origin. The (titano)magnetite grains occur in the matrix, their shape is angular, and their size is between several microns and more than 10 microns ( Figure 7B). Such characteristics of the (titano) magnetite grains allows us to interpret it as detrital origin. The occurrence of (titano)magnetite grains is much less frequent than the pyrite grains. No sub-micron to nano-size (titano)magnetite grains are found. The greigite occurs in the shape of framboids or irregular aggregates within and around aggregates of pyrite crystals, and their size is much smaller (nano-scale; presumably in the SP to stable SD size range) than the pyrite grains ( Figures 7C,D). The gregite grains are interpreted as being of post-depositional origin.
As described in Downcore Variation in Magnetic Properties Section, paramagnetic minerals, rather than ferromagnetic minerals, appear to control the bulk magnetic properties (such as bulk magnetic susceptibility) of the studied sediments. Candidates for these major paramagnetic minerals are common iron-bearing clay minerals with detrital origin, pyrite, iron monosulfide (mackinawite, FeS), and siderite (FeCO 3 ) of diagenetic (authigenic) origin in estuarine sediments. The gray or olive-gray colors of most of the studied sediments may be associated with abundant pyrite, which is consistent with the SEM observations. Several studies of other Holocene tidal flat sediments on the west coast of South Korea have reported authigenic siderites that were formed by early diagenesis during the Holocene through interactions with freshwater (Khim et al., 2000;Choi et al., 2003;Lim et al., 2004). Lim et al. (2004) also documented the siderite-abundant sediment interval (the unit T1 in their study) that was characterized by reddish-or yellow-brown (10 YR 5/4) massive mud and 10-12 × 10 −6 SI in bulk magnetic susceptibility.
Meanwhile, ferromagnetic minerals in the studied sediments also contribute to magnetic property variation, and leave meaningful signals. The IRM component unmixing analyses revealed occurrence of magnetic mineral components with different coercivity spectra from the total IRM. ferrimagnetic), probably of detrital origin; hematite (α-Fe 2 O 3 ; antiferromagnetic); and goethite (α-FeOOH; antiferromagnetic), respectively (e.g., Maxbauer et al., 2016). Note that the components 3 and 4 having highcoercivity spectra were not recognized in the SEM observations of this study. The component 2 possibly imply greigite (Fe 3 S 4 ; ferrimagnetic) due to the similar remanence coercivity values for greigite-bearing sediments (mainly 60-95 mT, but 37-98 mT for a wider range; Roberts, 1995;Snowball, 1997;Peters and Dekkers, 2003). This can be supported by the occurrence of ultra-fine greigite grains in SEM observation. We find that the relative contributions of these components differ at different depths, but generally the contribution of the components 1 and 2 dominates throughout the sediments. Peters and Thompson (1998) discriminated between different magnetic minerals by constraining the typical ranges for the respective minerals in a biplot of IRM 1T /χ lf vs. B cr . As shown in Figure 9, part of our YAR-4 data (n 11) did not fall into the range of (titano-)magnetite or greigite, but instead were superimposed partly on the range of the A5 sediments of New Jersey Miocene clays reported in Nilsson et al. (2013) who interpreted them as mixtures of (titano-)magnetite and greigite. It might be also possible that the plotted range of our data was influenced partly by the presence of hematite (and occasionally goethite). For the respective biplots, a fit line and its corresponding correlation coefficient (R) is provided. Black (Cyan) circles denote the data where the magnetic property values were obtained from the bulk (extract) samples. For the k LFQ vs. k FD plot, the data that were in out of the range between −5 and 20 in k FD (gray circles), being regarded less reliable, were excluded when consdiering the correlation degree between k LFQ and k FD .
Frontiers in Earth Science | www.frontiersin.org April 2021 | Volume 9 | Article 593332 Figure 10 presents FORC diagrams for four selective subsamples at different depths to provide additional information in terms of magnetic mineralogy and its domain state (also associated to grain size distribution of the magnetic mineral assemblage). Although the analyzed FORC diagrams are somewhat noisy, probably due to low concentrations of magnetic minerals, they still offer meaningful information. The FORC diagram for the 5.16-m-depth subsample had a B c peak at ∼10 mT and vertically suppressed divergence in contour distribution that was partly closed but the others intersected  Figure 8). The 7.84-m-depth FORC diagram had strong B c (∼10 mT) and weaker B c (∼27 mT) peaks, a closed but vertically spreading contour distribution that also elongated to higher coercivities along the B c axis, and a negative region distributed along the lower half of the B u axis. This pattern might be indicative of the presence of PSD greigite (e.g., Figure 5 in Roberts et al., 2011) in addition to SD-to-PSD magnetite. The 13.40-m-depth FORC diagram contained two B c peaks at ∼10 mT and near the origin. Its contours were spreading vertically wider at lower B c values, and were not closed but intersected the B u axis, being vertically stretched along the lower half of the B u axis. At low B c values, the contours were also distended diagonally downward. Such features possibly indicate the presence of PSD detrital magnetite (Roberts et al., 2018a) and SP magnetite or another magnetic mineral (Roberts et al., 2000). In this case, the SP-related feature could be interpreted to originate from ultra-fine greigite grains as indicated by the SEM observations. The 15.25-m-depth FORC contours had a B c peak at ∼18 mT, and larger elongation to higher coercivities along the B c axis. The dense contour concentration up to at higher  Figure 4), and the depths above which the sedimentation rate changes significantly (the same as identified in Figure 3; the triangles are the same as in Figure 4) are supplementarily shown.
Frontiers in Earth Science | www.frontiersin.org April 2021 | Volume 9 | Article 593332 B c values relative to the other samples might be related to the higher component 2 content at the 15.25-m-depth indicated by the IRM unmixing result (Figure 8). The contours had a vertical spreading pattern that got wider gradually as B c decreased toward the origin, and were partly closed but the others intersected the B u axis. The contour pattern is interpreted to indicate the presence of SD and PSD grains (e.g., Roberts et al., 2000). Overall, the FORC contour patterns tell us that the 13.40-m-depth sample showed predominance of (relatively large) PSD and SP ferromagnetic grains, but the other three samples showed predominance of SD to (relatively small) PSD grains. The SP grains was interpreted as being of diagenetic greigite. Also, none of the FORC diagrams shown had such a contour pattern of central ridge with mainly suppressed vertical spreading, which is indicative of the presence of biogenic magnetite of authigenic origin . Thus, it is interpreted that biogenic magnetite cannot be a candidate for the middle-coercivity mineral component indicated by the IRM component unmixing (the component 2 in Figure 8). Figure 11 presents our hysteresis parameter ratio data (M rs / M s and B cr /B c ) on a Day diagram (Day et al., 1977). These data mostly fall in the PSD region of Dunlop (2002a), and almost all of them are distributed in the region between the SD + MD 3 and the SD + SP (10 nm) mixing lines of Dunlop (2002a) and Dunlop (2002b). As pointed out by Roberts et al. (2018b), such Day diagram result should be taken with care in diagnosing domain state (or grain size) particularly for mixed assemblage of multiple magnetic minerals as presented in this study. Nevertheless, it can be useful to recognize downcore relative variation in magnetic granulometric trend. The data at shallow depths down to 7.92 m fall within a relatively restricted region that is at the up left-hand side of the data distribution, whereas the others at deeper than 7.92 m, except the 15.25 m one, are distributed within the down right-hand region with low values in M rs /M s and B cr /B c . With a comprehensive view from the IRM unmixing and FORC results, such variation on Day diagram is interpreted to be caused mainly by difference in grain size distribution of ferromagnetic minerals, rather than by variation in mineralogical composition. Then it is interpreted that samples being plotted at more toward the up lefthand side on the Day diagram have more abundant ferromagnetic grains with close to stable SD to SD-PSD boundary in size but less SP grains, relative to the down righthand plotted samples.

Diagenetic Effects on Downcore Variation in Magnetic Properties
In a wide range of sediments, but especially marine sediments from shallow continental shelf to deep-sea environments, early post-depositional diagenesis generally causes selective dissolution of relatively coarse micron-size detrital iron oxides (e.g., magnetite, hematite), and authigenic formation and growth of nano-to sub-micron-size magnetic minerals (e.g., pyrite, greigite, siderite), leading to systematic changes in the concentration and grain size of magnetic minerals with depth (e.g., Karlin and Levi, 1983;Liu et al., 2004;Rowan et al., 2009;Hatfield, 2014;Roberts, 2015). Due to the degradation of primary magnetic signals associated with paleoenvironmental changes, identifying and understanding the effects of such early diagenesis on magnetic minerals is crucial for interpreting environmental magnetic records (Snowball and Thompson, 1990;Verosub and Roberts, 1995;Liu et al., 2004;Demory et al., 2005).
Using our results from the depth interval of 3.2-13.9 m, at which sediment mean particle sizes and TOC contents vary little, we attempted to assess relatively straightforwardly early diagenetic effects on downcore variations in magnetic properties in the YAR-4 core.
In a general model for sediment cores influenced by early diagenesis, downcore variation exhibits significantly decreasing magnetic susceptibility with increasing depth in the suboxic zone (below the top shallow sediments (oxic zone) down to the suboxic-sulfidic (anoxic) boundary) and relatively constant low-susceptibility values at depths below the suboxic-sulfidic boundary, under anoxic zones. The magnetic susceptibility range in the 3.2-13.9 m depth interval of the core YAR-4 was similar in the order of magnitude of those observed in anoxic environments of coastal and shallow shelf sediment cores around South Korea and east China (Liu et al., 2004;Liu et al., 2005;Chen et al., 2015;Pan et al., 2017; see also Downcore Variation in Magnetic Properties Section). This result implies that the target interval of the YAR-4 sediments experienced considerable early reductive diagenesis under anoxic conditions, causing considerable dissolution of detrital magnetic minerals. It is apparently in agreement with the current anoxic conditions at depths below 0.03 m, where methane concentrations are higher, at a site in Yeongsan Lake . This suggestion is also consistent with the widespread presence of the middle-coercivity component (component 2 in Figure 8), i.e., greigite (see Magnetic Mineralogy and Granulometry from Magnetic Measurements and FIGURE 9 | Biplot of IRM 1T /χ lf vs. B cr for the YAR-4 core sediments (red circles; n 11), with typical ranges (gray boxes) of the magnetic parameters for different minerals proposed by Peters and Thompson (1998) and data for the A5 sediments reported by Nilsson et al. (2013)  Scanning Electron Microscope Observations Section), which is a typical authigenic iron mineral, in addition to pyrite and iron monosulfide, which are generated in the sulfate-methane transition zone (SMTZ) (Roberts, 2015). Rowan et al. (2009) identified downcore variations producing counterclockwise loop trends in Day diagram from marine sediment cores around the world (Figure 9 of Rowan et al., 2009). They ascribed the systematic trends to the superposition of spatio-temporally progressive diagenetic changes of magnetic minerals. This changes include progressive reductive dissolution where finer pre-existing ferromagnetic grains (mostly detrital magnetite) preferentially dissolve to leave a coarser assemblage, and progressive sulfidization resulting in SP greigite nucleation and its progressive growth into the SD size range. The trends occur as the result of magnetic property combinations between survived coarse detrital grains and SP to SD authigenic greigite grains. The looping trend and the YAR-4 downcore variation in Day diagram (Figure 11) apparently share a similar pattern in part, but visible difference between them exists. The whole analyzed interval of the YAR-4 core should be correlated to the "zone 3" of the downcore magnetic property profile shown in Figures 2, 3 of Rowan et al. (2009), based on being the low magnetic susceptibility values throughout the analyzed interval ( Figure 4). However, the YAR-4 data in Day diagram lie within a wider B cr /B c range, being stretched toward the down right-hand region, than those of the "zone 3" (Rowan et al., 2009). Given this, the YAR-4 downcore variation is inferred to contain more complicated diagenetic effects than those indicated by Rowan et al. (2009) and/or possibly other effects caused by changes in depositional condition.
Magnetic mineral diagenesis is also generally thought to be controlled by depositional conditions, that is, mainly TOC and sedimentation rate (e.g., Hesse and Stolz, 1999). Changes in TOC and sedimentation rate often produce non-steady-state early diagenesis, where diagenetic reactions change with depth as the depositional conditions change (e.g., Thomson et al., 1984;Robinson et al., 2000;Emiroglu et al., 2004;Fu et al., 2008;Roberts, 2015). If magnetic property variations of the 3.2-13.9 m depth interval of the YAR-4 core are further ascribed to such nonsteady-state diagenesis, the YAR-4 downcore variations should be associated with changes in sedimentation rate because it is expected that the constant low TOC contents throughout the YAR-4 core would bring minimal or negligible effects. Riedinger et al. (2005) documented distinct minima in magnetic susceptibility within a distinct sediment interval, and interpreted them to result from a drastic change in sedimentation rate, which resulted in a fixation of the SMTZ at a specific depth causing substantial change of pre-existing abundant iron (oxyhydr)oxides to iron sulfides. Zheng et al. (2011) proposed that drastic changes in sedimentation rate resulted in vertical shifts of the SMTZ, of which each SMTZ eventually exhibited low magnetic susceptibility, high coercivity, and significant paramagnetic contribution as a result of dissolution of pre-existing detrital magnetite and hematite but strong resistance to the dissolution for hematite relative to magnetite, and their replacement mainly by paramagnetic pyrite. Abrajevitch and Kodama (2011) reported a relationship among mineral magnetic features, i.e., low magnetic susceptibility, coarsening of detrital magnetic grains (e.g., magnetite), increasing relative abundance of hematite, and decreasing relative abundance of goethite during periods with low sedimentation rate (i.e., sea-level highstand, in the case).
Indeed, the YAR-4 core appears to record at least four drastic changes in sedimentation rate and three distinct intervals with relatively low sedimentation rate, based on our age-depth model ( Figure 3). However, none of the relations between magnetic properties such as those identified in the above-stated previous studies are identified in the YAR-4 core (see e.g., Figures 4, 8). We thus infer that, in addition to the above-discussed possible diagenetic effects, additional effect(s) caused by another factor(s) should be contained in the downcore magnetic property variations for the YAR-4 core, and speculate the possible presence of certain paleoenvironmental change(s) or even lessknown diagenetic processes linked to that (in similar manner as reported in e.g., Larrasoaña et al., 2003, Blanchet et al., 2009, and  the review of Roberts, 2015) for the YAR-4 site as the additional controlling factor(s).

Searching for Magnetic Property Features Associated With Paleoenvironmental Changes
Sedimentation in the studied YAR-4 site during the Holocene would have been considerably influenced by sea-level change and hydrologic events that are regarded as closely associated with regional climate changes such as the ENSO activity, the EAM, and the HCO. Accordingly, types, concentration, and grain size distribution of the contained magnetic minerals would change through time, at least partly in response to the sea-level change and/or regional climate changes. Based on our age-depth model (Figure 3), the YAR-4 record of magnetic property variabilities presented in this study covers a period of approximately 4,800-9,000 cal. yr BP (Supplementary Figures S2, S3 and S4). In order to search for magnetic property proxies for paleoenvironmental changes, we compare temporal variations between some of the obtained magnetic properties (we choose k LFQ , B c , B cr , M rs /M s , and each of IRM-unmix components for this) and proxy indicating the change in sea-level (Lee and Chang, 2015;Song et al., 2018), the ENSO activity (Moy et al., 2002), the summer EAM (Dykoski et al., 2005), and the HCO with the 8.2 ka cooling event  around the Korean Peninsula, respectively ( Figure 12, Supplementary Figures S2, S3 and S4).
In advance of getting to the main discussion, we would like to mention that the use of only initial bulk magnetic susceptibility without detailed additional information about constituent magnetic minerals in interpreting a certain paleoenvironment change (as the same manner utilized in previous studies for Korea with similar purposes to this study) appears to be imprudent, given the k LFQ values that were substantially influenced by complicated magnetic mineral diagenesis as in the studied case.
Comparisons With El Niño Southern Oscillation Activity, Summer Eastern Asian Monsoon, and Holocene Climate Optimum Moy et al. (2002) presented a record of red color intensity variability through the Holocene as a proxy representing the ENSO activity variability from a sediment core of the lake Laguna Pallcacocha in the southern Ecuadorian Andes. There exist three extremely high peaks in ENSO activity during ∼4,800-9,000 cal. yr BP (∼4,750-5,000, ∼5,900-6,000, and ∼7,900-8,000 cal. yr BP; green shaded regions in Supplementary Figure S2), and such ENSO peaks are being considered to be related to high frequency and magnitude of flooding or freshwater input from inland at least in the southern part of the Korean Peninsula (Lim et al., 2017, Lim et al., 2019. However, we did not find any of our magnetic property variability data that could be well correlated to those ENSO peaks (Supplementary Figure S2). This implies no or, if any, meager regulation by extreme short-term hydrologic events (also associated with ENSO activity) on the magnetic property variabilities in the YAR-4 site.
FIGURE 11 | Day plot of our YAR-4 core sediment data (n 37) at different depths and sedimentary units. Circle and triangle symbols denotes the data measured from bulk and extract samples, respectively. Some data are labeled using the "depth in core" in meter, among which the data labeled in bold denote those where the associated FORC diagrams are provided in Figure 10. Single-domain (SD), pseudo-single-domain (PSD), and multi-domain (MD) boundaries, trend curves 1, 2, and 3 for SD + MD mixtures of magnetite and trend curves for SD + SP with 5-nm and SD + SP with 10-nm mixtures of magnetite were drawn by Dunlop (2002a), Dunlop (2002b). Colorshaded region represents a typical range for SD + SP mixtures of greigite indicated by Roberts et al. (2011). Division of the data into two groups by dashed enclosed lines was based on ranges of the B c , B cr , and M rs /M s values (see text for details).
FIGURE 12 | Comparison of our age-depth model for the YAR-4 core with two different suggested sea-level curves for the western coast of South Korea (Lee and Chang, 2015;Song et al., 2018). HBBM intervals (characterized by subtantially high B c , B cr , and M rs /M s values) and their corresponding ages are indicated by red lines and red shaded boxes. The YAR-4 age-depth model curve with the uncertainties are the same as in Figure 3. Sedimentary unit divisions are also supplementarily shown.
The oxygen isotope ratio variability from the Dongge Cave in the South China by Dykoski et al. (2005), adopted worldwide as representative of the summer EAM variability, seems to have similarity partly in centennial (short-term excursion; blue shaded regions shown in Supplementary Figure S3) and millennial (long-term trend; intervals between blue shaded individuals in Supplementary Figure S3) scales to our B c , B cr , and M rs /M s variabilities (Supplementary Figure S3). It shows apparently a relationship of low B c , B cr , and M rs /M s values with strong summer EAM. Strong summer EAM generally accompany high accumulation of rainfall precipitation that leads to high terrestrial material input, which can be possibly accompanied with increasing concentration and grain size of detrital magnetic minerals (e.g., magnetite). If it can be assumed that the low B c , B cr , and M rs /M s values were carried mainly by increasing grain size of detrital magnetite for the YAR-4 sediments, these magnetic properties can be dealt as summer EAM proxies. According to this inference, strong summer EAM conditions should be favorable for coarsening of detrital magnetite, via controlling the local hydraulic condition around the studied area. On the other hand, if so, lowering B c , B cr , and M rs /M s values should occur at the periods with extremely high ENSO peak (∼4,750-5,000, ∼5,900-6,000, and ∼7,900-8,000 cal. yr BP), through a similar process occurred in strong summer EAM periods with enhanced rainfall frequency and magnitude. However, such relationship appears not visible (Supplementary Figure S2). Therefore, we interpret that such potential hydraulic effect in relation to detrital magnetite grain size does not govern, but may control partly, the B c , B cr , and M rs /M s variabilities.
Additionally, the HCO period and the 8.2 ka (thousand years ago) abrupt cooling event are also the major phenomena of climate change during the Holocene. The HCO is a warm and humid period with relatively high precipitation, and occurred during ∼4,800-7,600 cal. yr BP for the southern Korean Peninsula cf., ∼5,100-8,200 cal. yr BP, a calibrated age range in this study with the OxCal and IntCal20 from ∼4,500 to 7,400 years BP in uncalibrated age originally reported by Yang et al., 2008, or ∼5,000-6,300 cal. yr BP by Nahm and Hong, 2014). The 8.2 ka event is a short-term climate excursion to cold and dry condition at ∼8,200 cal. yr BP, and has been identified also in South Korea (Park et al., 2018;Park et al., 2019). It might be good to look at the high (detrital) magnetite relevant abundance (i.e., the IRM-unmixed component 1 contribution) with less greigite relevant abundance (i.e., the IRM-unmixed component 2 contribution) during the HCO period, and vice versa around the 8.2 ka event (Supplementary Figure S4). Blanchet et al. (2009) reported that large flood deposits and glacial deposits are favorable for preservation of diagenetic greigite, and interpreted it to be related to dominance of reactive iron over organic matter and/ or HS − , i.e., enrichment in terrigenous sediments. According to this, it might be possible that the 8.2 ka cooling triggered relatively enhanced greigite preservation through the early diagenetic process at the studied YAR-4 site. This potential relationship will be worth being tested further in future investigations with more high-resolution magnetic analyses and the use of multiple cores at different sites.

Possible Magnetic Response to Abrupt Drop in the Rate of Sea-Level Rise
Many sea-level change curves have been suggested for the western coast of South Korea during the Holocene (e.g., Bloom and Park, 1985;Hwang et al., 1997;Chough et al., 2004;Lee et al., 2008;Lee and Chang, 2015;Song et al., 2018); however, no consensus has been reached (cf. the review by Choi, 2018). A major discrepancy among them would be the sea level relative to the present mean sea level (MSL) and its change pattern since ∼7,000 cal. yr BP: One group of the sea-level curves depicts that the sea level always was lower than the present MSL and continued a gradual rising toward the present MSL (e.g., Bloom and Park, 1985;Chough et al., 2004;Lee and Chang, 2015), whereas another group suggests a sea-level highstand (higher sea-level relative to the present MSL) during ∼4,000-7,000 cal. yr BP followed by a sealevel drop toward the present MSL (e.g., Hwang et al., 1997;Song et al., 2018). Unfortunately, the sedimentary characteristics and our age-depth model of the YAR-4 core (Figures 2, 12) independently do not allow to fully support one of the two different groups of suggested sea-level curves.
In any case, most suggested sea-level curves of these previous studies depict a drop in the rate of sea-level rise broadly at around 7,000 cal. yr BP. Besides this, the sea-level curve of Song et al. (2018) depicts an additional drastic drop in the sea-level rise rate at ∼8,200 cal. yr BP, then shortly followed by an increase of this rate again, inferring formation of a concomitant short-term slowstand/stillstand ( Figure 12). The existence of such "stepped" change in sea-level rise rate as indicated by Song et al. (2018) is being identified even in more recent studies (e.g., Tanabe, 2020). It also would be good to mention that attempts of the west coast sea-level reconstruction during the period back to prior to 8,000 cal. yr BP were rare in among the suggested sea-level curves (cf. Figure 1A in Choi, 2018). Lee and Chang (2015) generated the curve by fitting to the sea-level proxy data that were scarce but highly scattered in the sea-level value during prior to ∼7,800 cal. yr BP (Figure 2 of Lee and Chang, 2015). This may be the reason why the ∼8.2 ka potential drop in the rate of sea-level rise does not appear in the curve of Lee and Chang (2015) (and other previous studies depicting similar curves).
For the YAR-4 core, one remarkable, first-order magnetic feature is the abrupt increases in B c , B cr and M rs /M s (hereafter referred to as "HBBM feature") observed at depths of 15.25 and 7.84-7.96 m with their persistent high values at depths shallower than 7.84 m (Figure 4), corresponding to ∼8,100 cal. yr BP and since ∼6,700 cal. yr BP in inferred age (Figure 3; Supplementary Table S6). This observation seems not to be fully explained by solely the sulfidization-dominated diagenetic processes (even if sedimentation rate variability were considered) or correlated well with any of the ENSO activity variability (linked to the local hydrologic events), the summer EAM variability, and the HCO period, as interpreted above. Instead, the respective ages of the two abrupt B c , B cr , and M rs /M s increases are well coincident with those at which the rate of sea-level rise started to decrease suddenly according to the Song et al. (2018) sea-level curve, given the intrinsic uncertainty of the age-depth model construction (Figures 3,  12; Supplementary Table S6). These concordance allows us to infer that the HBBM feature is strongly related to change of the sea-level rise rate and possibly, sea-level stillstand/slowstand and highstand. The apparent absence of HBBM feature at the old stillstand/slowstand ( Figure 12) possibly might be attributed to its short-lived period.
It is worthwhile to discuss and constrain concerning possible mechanism to explain the relation between the change of the sea-level rise rate and the HBBM magnetic feature. A possible major factor enhancing B c and B cr could be among increasing relative abundance of hematite and goethite, increasing relative abundance of greigite, and decreasing relative abundance of detrital magnetite, or a combination of them. However, no similar characteristic concerning relative abundance of ferromagnetic minerals between the HBBM intervals ( Figure 8; see also Magnetic Mineralogy and Granulometry from Magnetic Measurements and Scanning Electron Microscope Observations Section) allows us to conclude its less probability of being major factor leading to enhancing B c and B cr . It, however, would like to be mentioned that, within the 3.21-7.04 m depth interval, there are decreasing magnetite with increasing more high-coercivity minerals in relative abundance, which may be attributable at least partly to the HBBM feature. On the other hand, the grain size distribution of ferromagnetic minerals between the HBBM-bearing intervals and the others could be clearly discriminated; the HBBM intervals showed more contents of stable SD to SD-PSD boundary grains but less SP contents (Figures 10, 11; see also Magnetic Mineralogy and Granulometry from Magnetic Measurements and Scanning Electron Microscope Observations Section). Moreover, the occurrence of the HBBM features during the HCO and the 8.2 ka event that had different climate conditions (warm and humid vs. cold and dry; Supplementary Figure S4) indicates its independence of both surface temperature and humidity. The HBBM occurrence also does not appear to be associated with the sedimentation rate ( Figure 12). Unfortunately, the exact mechanism driving the relation between the abrupt change in the grain size distribution and the abrupt drop in the sea-level rise rate is not identified by this study. However, we note that the varying sediment transport distance and pathway driven by the sea-level rise (e.g., Chen et al., 2017), leading to grain size change of detrital magnetic mineral sources, is not likely responsible for the observed abrupt change of B c , B cr , and M rs /M s values. We here modestly speculate one possibility that biochemical condition of the uppermost sediment involving the SMTZ and bottom water chemistry during the periods with low rate of sea-level rise might be favorable for rapid growth of greigite to the SD size range at the SMTZs (very close to the sediment top layer, in the studied site). Conversely, a certain condition during the periods with high rate of sea-level rise might permit forming SP greigite but inhibit its growth to the SD size range at the SMTZs that were shifted progressively upward. Consequently, the contrast between the grain size distributions driven by such sea-level change that controlled diagenetic modifications possibly might be the major cause of the occurrence of the HBBM features. This speculation will be verified in additional future works for firm establishment of a magnetic property proxy for the abrupt drop of the rate of sea-level rise or the initiation/occurrence of sea-level stillsatnd/slowstand.

CONCLUSION
In this study, we analyzed 3.2-to-19.8-m depth interval of a 20m-long Holocene muddy sediment core recovered from the Yeongsan Estuary, South Korea, using mineral magnetic measurements with sediment particle size and TOC content to characterize downcore variations in a variety of magnetic properties. We then evaluated the diagenetic effects on magnetic signals and tested their availability as proxies of paleoenvironmental change. The magnetic measurements included magnetic susceptibilities, hysteresis parameters, progressive IRM acquisition, and FORC analysis for each of the selected subsamples. The major findings of this study are as follows: (1) The analyzed sediments were generally characterized by relatively low bulk magnetic susceptibility values (∼19-245 × 10 −6 SI) and the predominance of paramagnetic (>60% of the total) rather than ferromagnetic contribution, but undoubtedly distinct downcore variations in ferromagnetic-related properties.
(2) The analyzed sediments would have undergone substantial early diagenetic alteration including dissolution and transformation of detrital magnetic minerals and authigenic magnetic mineral formation and growth, which eventually might led to the current complex variability in magnetic properties. For this reason, the raw bulk magnetic susceptibility values are not recommended as a stand-alone proxy of environmental change in (and around) the studied area.
(3) Abrupt increase in B c , B cr , and M rs /M s is coincide well with abrupt drop in the rate of sea-level rise, in other words, the initiation/occurrence of sea-level stillstand/slowstand or highstand, during the Holocene. It is worthwhile to clarify the exact mechanism causing their potential linkage and further consider the potential of these magnetic properties as a proxy, at least in and around the studied area (the western coast of South Korea), in future study. (4) We also have preliminarily explored possible relationships of the magnetic property variabilities with the ENSO activity involving local hydrologic events, the summer EAM, the HCO and the 8.2 ka cooling event. Among which, the possible relationship between the magnetite and greigite relative abundance and the HCO or the abrupt short-term cooling appears to be worth being tested further in future study.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.