Sedimentary Rock Magnetic Response to Holocene Environmental Instability in the Pearl River Delta

Located on the northern coast of the South China Sea, the densely populated Pearl River Delta has experienced the combined effects of sea-level change, monsoon-driven discharge, and especially human activity, since the late Holocene. However, how these factors have regulated the regional environmental and sedimentary evolution remains unclear. To better understand these processes, we conducted a high-resolution rock magnetic investigation of the Holocene sediments of core DS01, drilled in the vicinity of the West River channel in the head area of the Pearl River deltaic plain. The magnetic grain-size proxy of the ARM/κlf ratio (the ratio of anhysteretic remanent magnetization to low-field magnetic susceptibility) indicates a long-term fining trend of the magnetite grain size, which may be a response to an increase in the weathering intensity in the Asian monsoon region during the Holocene. An interval with an enhanced concentration of magnetic minerals (mainly magnetite and hematite) occurred during 7.7–4.8 kyr BP (calendar years before 1950), coinciding with a period of delta progradation. During the marine transgression in the early Holocene, two similar intervals of magnetic enrichment may reflect regional hydrodynamic shifts associated with cooling events at ∼9.5–9.3 kyr BP and 8.2 kyr BP. The subsequent 4.2 kyr BP cooling event possibly induced a cold and dry environment in the sediment source area. From ∼800 yr BP onward, there was a major increase in the sedimentary magnetic mineral content, likely in response to intensified agricultural and industrial activities.


INTRODUCTION
The formation of the sedimentary sequences of the Pearl River (Zhu Jiang) Delta in southern China likely began in the late Quaternary (e.g., Huang et al., 1982;Zhao, 1990). Developed on the northern coast of the South China Sea, in an area undergoing long-term tectonic subsidence, the regional sedimentary sequences of the deltaic plain have documented the complex history of marine transgression and regression (e.g., Zong et al., 2009;Wei and Wu, 2011;Zong et al., 2012;Zong et al., 2016;Wei X. et al., 2020). The modern Pearl River Delta is a densely populated and industrialized area, which is strongly influenced by the Asian summer monsoon (An, 2000;An et al., 2000). Determining the environmental and sedimentary evolution of the Pearl River Delta is a complex task, especially during the Holocene, when human activities became an increasingly important influence. Stalagmite records from China have provided detailed information on the paleoclimatic evolution of the Asian monsoon over the last 640 kyr (Cheng et al., 2016). For the latest 10 kyr, from the early Holocene onward, precisely dated proxy environmental records are available, such as the stalagmite records from southern China (Yuan et al., 2004;Dykoski et al., 2005;Wang et al., 2005). The δ 18 O record from Dongge Cave suggests a weakening trend of Asian monsoon intensity for most of the Holocene, which was punctuated by pronounced short-term monsoon events that can be correlated regionally or globally (Wang et al., 2005).
An integrated analysis of a transect of sediment cores recovered from along both sides of the current channel of the West River (Xi Jiang), a tributary extending west-east to the Pearl River catchment, has provided a Holocene record of the deltaic evolution of the head area of the Pearl River deltaic plain ( Figure 1) (Xiong et al., 2018;Fu et al., 2020). Within the area, lithostratigraphic units deposited above bedrock can be FIGURE 1 | Location map (A,B) and stratigraphic information (C) for sediment core DS01 (red circle) in the head area (dotted rectangle) of the Pearl River deltaic plain. Black circles mark the locations of sites discussed in the text: Dongge Cave (Wang et al., 2005), Dahu Swamp (Zhou et al., 2004), Huguangyan Maar Lake in south China (Wang et al., 2016), core HKUV11 (8.4 m water depth, Wu et al., 2017), core MD05-2905 (1,647 m water depth, Zhou et al., 2012) in the northern South China Sea, core NZ2 (4.6 m water depth) in Lingdingyang Bay (Pearl River estuary, Wu et al., 2021), and core PD and core SX97 from the banks of the Pearl River estuary Peng et al., 2014). The 10 radiocarbon dates for core DS01 shown in (C) are from Fu et al. (2020). Details of the methods used to produce the age data are given in Xiong et al. (2018). The maps were produced using Tencent Map App.
Frontiers in Earth Science | www.frontiersin.org June 2022 | Volume 10 | Article 882201 subdivided into three major groups: fluvial deposits, estuarinedeltaic deposits, and delta plain deposits, with decreasing depth. These sedimentary facies correspond to different sedimentary environments in the head area of the Pearl River Delta, and they provide an opportunity to obtain a relatively complete regional paleoenvironmental record throughout the Holocene.
In this study, we investigated a Holocene sedimentary sequence from the Pearl River Delta using rock magnetic measurements, which have been widely used for paleoenvironmental reconstruction (Thompson and Oldfield, 1986;Evans and Heller, 2003;Liu et al., 2012). Originating from the adjacent terrigenous areas of the Pearl River drainage, the detrital magnetic minerals in the sediments can potentially be used to decipher the paleoenvironmental history of the drainage area. Rock magnetic results have been reported for various drilling sites in the Pearl River estuary (e.g., Yang et al., 2008;Peng et al., 2014;Wu et al., 2021) (Figures 1A,B). These studies have provided abundant information about the interplay of sea-level changes and sedimentary evolution in the coastal area since at least the late Pleistocene. In the present study, we measured multiple rock magnetic parameters, at high stratigraphic resolution, from a sediment core from the delta, with the aim of reconstructing millennial-to centennial-scale records of the Holocene sedimentary evolution of the head area of the Pearl River Delta. We compare the records with previous reconstructions of past regional climatic and environmental changes and attempt to understand how changes in climate, sea-level, and human impacts influenced the Holocene evolution of this coastal area.

MATERIALS AND METHODS
We used sediment core DS01 ( Figure 1B; 23°13.47′N, 112°44.43′E; 25.61 m in length), obtained from the north side of the West River. The altitude of the coring site was~4.4 m above the sea level. The base of the core consists of compact, reddish weathered bedrock. Above, the sediments can be divided into three lithological units, as illustrated in Figure 1C. An age of 18,181-17,865 cal. yr Before Present (BP) was obtained for a sand layer (within Unit I) that unconformably covered the base rock. The three units are described as follows: Unit I: (25.61-23.6 m), sand and clayey silt with occasional plant residues. Unit II (23.6-5.6 m), mainly dark silt with occasional plant residues; soft. This unit comprises the main part of the core. Unit III: (5.6-0 m), mainly gray silt with occasional plant fragments; loose structure. Fu et al. (2020) identified four subunits that should comprise a complete Unit III within the vicinity of the head area of the Pearl River Delta (Figure 1). However, in core DS01, two of these subunits are missing (sub-units IIIa and IIIb). Thus, we infer the existence of another sedimentary discontinuity between Unit II and Unit III.
Although not continuously represented in core DS01, the sedimentary succession in the study area was found to record a progressive transition from an estuarine to a deltaic environment, associated with the upper part of Unit II ( Figure 1C). In this study, samples were collected every 2 cm by pressing plastic cubes (2 × 2 × 2 cm 3 ) into the center of the archive half of core DS01. The surface sediments (upper~110 cm) were not sampled to avoid sediments affected by recent agricultural and industrial activities. Another set of samples was collected at 10-20 cm intervals for analyses of grain size, organic carbon isotopes, and diatoms. These results were previously published (Xiong et al., 2018;Fu et al., 2020) and will be partially aiding the following discussions in our study.
Low-frequency (976 Hz) and high-frequency (15,616 Hz) volume magnetic susceptibilities (κ lf and κ hf , respectively) were measured on a total of 570 discrete cubic specimens (equivalent to a sampling resolution of 4.3 cm), using an AGICO Kappabridge MFK1-FA magnetic susceptibility meter. The volume frequencydependent susceptibility (κ fd ) was calculated as κ fd = κ lf -κ hf . κ lf generally reflects the concentration of paramagnetic and ferromagnetic components, and κ fd is sensitive to the presence of superparamagnetic grains (Thompson and Oldfield, 1986;Evans and Heller, 2003;Liu et al., 2012).
Anhysteretic remanent magnetization (ARM) was imparted in a peak alternating field (AF) of 100 mT with a direct current biasing field of 0.05 mT superimposed, using a 2G AF demagnetizer. ARM was measured using a 2G three-axis cryogenic superconducting rock magnetometer (755R) installed in a shielded room with residual fields of <300 nT, in the Paleomagnetism Laboratory of the South China Sea Institute of Oceanology, Chinese Academy of Sciences. ARM is sensitive to stable single-domain ferrimagnetic grains (e.g.,~40 nm for magnetite). We calculated the ARM/κ lf ratio as a proxy for detecting relative changes in the grain size of magnetite (Banerjee et al., 1981;King et al., 1982).
Saturation isothermal remanent magnetization (SIRM) was imparted in a field of 1T using an ASC Impulse Magnetizer (IM-10-30) and measured using an AGICO JR-6A rock magnetometer. Back fields of 100 mT (IRM -100 mT ) and 300 mT (IRM -300 mT ) were then imparted and measured sequentially. Back IRMs were imparted to every third sample (n = 190). The SIRM and backfield IRMs are used to calculate hard IRM (HIRM) and S-ratio (Thompson and Oldfield, 1986;Evans and Heller, 2003) as follows: SIRM mainly reflects the combined concentrations of ferrimagnetic (e.g., magnetite) and imperfect antiferromagnetic particles (e.g., hematite and goethite), predominantly the former, unless the magnetic mineral assemblages are dominated by imperfect antiferromagnetic material. HIRM and S-ratio are widely used as proxies for the absolute and relative concentrations of imperfect antiferromagnetic particles, respectively. The ARM/SIRM ratio was calculated as an additional magnetic grain-size proxy, since it is unaffected by paramagnetic and superparamagnetic minerals. In addition to the aforementioned routine magnetic measurements, selected samples were used for the following more detailed magnetic measurements. High-temperature magnetic susceptibility measurements (κ lf -T) were made on freeze-dried powder samples in an argon environment, using the Kappabridge MFK1-FA combined with a CS-4 high-temperature furnace. IRM acquisition curves were measured in magnetic fields ranging from 1 mT to 1.5 T (in 90 logarithmically distributed steps), followed by backfield demagnetization in a range of magnetic fields up to 100 mT. Hysteresis loops were obtained by cycling the magnetic field between ±1.0 T, with a step size of 5 mT. These measurements were made using a Lakeshore 8600 vibrating sample magnetometer (software version 1.3) at the Environmental Magnetism Laboratory of South China Normal University.

Magnetic Mineralogy
The thermomagnetic curves (κ lf -T) of the five representative samples are largely consistent, showing a slight decrease in κ lf from room temperature to~300°C during heating (see the lower panel in Figure 2). This may be caused by the gradual unblocking of paramagnetic (or superparamagnetic) minerals, such as siderite and pyrite (Dunlop and Özdemir, 2001;Roberts, 2015). Above~300°C, κ lf begins to increase and a "hump" is evident between 400 and 550°C. The cooling limb of the curve is not reversible because stoichiometric ferrimagnetic minerals have been produced during heating, possibly from the thermally induced alteration of paramagnetic iron-bearing silicates and/ or clays (e.g., Deng et al., 2001;Liu et al., 2020). All the curves show a sharp decrease near 580°C, the Curie temperature of magnetite, indicating that magnetite is likely the major magnetic carrier of the samples. The decreasing trend continues beyond 600°C, revealing the presence of a magnetic phase with a Curie temperature greater than that of magnetite. A weak convex pattern appears at~120°C during each of the cooling runs, likely arising from the presence of a minor amount of goethite. Figure 3A shows that the IRMs at 300 mT contributẽ 73-88% to the SIRM values, with S-ratios ranging from 11.5 to 95.3%, with a~61% median for the whole core (n = 192). The representative samples have a wide range of coercivities of remanence (~35-105 mT) ( Figure 3B). We conclude that the dominant remanence carriers are low-coercivity ferrimagnetic minerals, but with high-coercivity phases being abundant.
The hysteresis loops of the selected samples are slightly waspwaisted ( Figures 3C-G), which may result from the coexistence of two different magnetic components with contrasting coercivities (Roberts et al., 1995;Tauxe et al., 1996). Furthermore, the loops are almost closed when approaching 600 mT, which also demonstrates the coexistence of both lowand high-coercivity magnetic components in the sediments.

Down-Core Variations in Magnetic Properties
Down-core plots of multiple magnetic parameters are presented in Supplementary Figure S1. Using the pre-existing age data for the core ( Figure 1C), we built an age model based on linear FIGURE 2 | Temperature-dependent magnetic susceptibility (κ lf -T) curves for representative samples from core DS01. Arrows indicate heating (red) and cooling (blue) runs.
Frontiers in Earth Science | www.frontiersin.org June 2022 | Volume 10 | Article 882201 interpolation or extrapolation ( Figure 4). However, this was not possible for Unit I (25.6-23.6 m) because the evaluated ages for the upper part would be much younger than those of the base of Unit II, assuming a constant sedimentation rate (~16.5 cm/kyr) estimated from the pre-existing age data. We infer that Unit I was rapidly deposited in a high-energy freshwater environment. Two age data are available for Unit III, and they indicate a high sedimentation rate (~5 m/kyr) within the interval of 1308-420 yr BP, resulting in a clear age discrepancy between Unit III and Unit II. This sedimentary hiatus may have been caused by river channel migration or later human activities (Fu et al., 2020). Six age data are available for Unit II, and therefore the chronology for this lithologic unit is of a relatively high resolution. Within Unit II, four horizons are highlighted due to the relatively high values of most of the magnetic parameters (Supplementary Figure S1, Table 1 and Figure 4). Interval (a) (~22.0-23.6 m,~9.3-9.5 kyr BP) is a plateau-like zone of the magnetic parameters, at the base of Unit II, corresponding to the period when an estuarine environment began to dominate in the study area. In interval (b), around the middle of Unit II, a minor peak is evident in almost all the magnetic parameter FIGURE 3 | IRM acquisition curves and backfield IRM curves (A,B) and hysteresis loops (C-G) for representative samples from core DS01. Ms, saturation magnetization; Mrs, saturation remanence; Bc, coercivity; Bcr, coercivity of remanence. Hysteresis loop data of the samples shown in (C-F) were processed using the HystLab program (Paterson et al., 2018).
Frontiers in Earth Science | www.frontiersin.org June 2022 | Volume 10 | Article 882201 6 profiles. However, intervals (a) and (b) do not have comparably high HIRM values, which indicates the dominant contribution of ferrimagnetic minerals (magnetite) in these intervals. Except for intervals (a) and (b), Unit II shows increasing but slightly different values of the magnetic parameters, including κ lf , ARM, SIRM, HIRM, S-ratio, and ARM/κ lf (Supplementary Figures S1C,E-J). Therefore, these trends did not develop simultaneously. Within interval (c), κ lf , ARM, SIRM, HIRM, and S-ratio fluctuate within an interval of overall higher values. Interval (d), at the top of Unit II, is characterized by rapid changes in most of the magnetic parameters, with HIRM showing a significant peak, indicating a peak in the concentration of hard magnetic components (mainly hematite). In contrast, there are only minor peaks in κ lf , ARM, and SIRM.
There is a general upward-increasing trend in ARM/κ lf , as shown in Supplementary Figure S1I, suggesting the progressive fining of the magnetic grains (mainly magnetite). The ARM/ SIRM ratio shows relatively uniform values throughout most of the core, with no comparable fluctuations in the interval of 11.3-6.8 m [interval (c)], but with an increase thereafter (Supplementary Figure S1J).
Cross-plots showing the relationships between the magnetic parameters are shown in Supplementary Figure S2. For κ lf , κ fd , ARM, and SIRM (Supplementary Figures S2A-C,F, respectively) linear relationships are evident for most of the data, indicating a general consistency of the sedimentary magnetic assemblages; deviations in these linear relationships generally correspond to high-amplitude fluctuations of the respective parameters within the profiles. κ fd was not expected to be closely related with both ARM and SIRM (Supplementary Figures S2D,E) as they are indicative of distinctive magneticdomain states; however, the data points in both plots are generally clustered, as the magnetic minerals are seemingly concentrated within sedimentary intervals with a higher clay content (Figure 4).

Hard Isothermal Remanent Magnetization and S-Ratio
The rock magnetic results show that magnetite and hematite are the dominant magnetic minerals in core DS01 and that the parameters related to magnetic mineral concentration (κ lf , ARM, SIRM, HIRM, and S-ratio) show very similar stratigraphic patterns. The depositional context of the core location varied substantially since~18 kyr BP, and the material sources of the core sediments may have undergone substantial changes. This inference may undermine the use of HIRM and S-ratio, as they may be significantly affected by changes in magnetic mineral coercivity. Therefore, we used the L-ratio of Liu et al. (2007) to evaluate if these parameters were influenced by Al-substituted hematite or goethite, or both, since this would result in a wide range of coercivities which would invalidate their use. The L-ratio is calculated as follows: Figure 5A shows a nearly constant distribution for most of the data points. The outliers are related to the specific sediment layers in Figure 5B and Figure 4 and coincide with the transitional interval from Unit II to Unit III (Supplementary Figure S1), during which the depositional environment was very unstable. The horizontally shaded intervals in Figure 5B are characterized by magnetic minerals with substantially different coercivities, which may reflect different sources; hence, the use of HIRM and S-ratio may be inappropriate in these intervals. However, in the other intervals, both parameters can be used to quantify the absolute and relative concentrations of highcoercivity magnetic minerals Roberts et al., 2020).

DISCUSSION
The magnetic properties of sedimentary archives are often correlated with geological processes that control the provenance, transportation, deposition, and post-depositional modification (e.g., diagenesis) of sediments. The thermomagnetic analyses suggest a relatively uniform sedimentary magnetic mineral composition (Figure 2), likely suggesting a consistent sediment source. This inference is supported by the generally consistent coercivities of the magnetic minerals indicated by the L-ratio (Figure 5), as well as by the similar hysteresis loops of the representative samples ( Figures  3C-G). However, there are significant fluctuations in magnetic mineral concentration and in magnetic grain size, as shown by the stratigraphic variations of κ lf , ARM, SIRM, HIRM, and S-ratio and by the magnetic grain-size indicators, ARM/κ lf and ARM/ SIRM. We now focus on these fluctuations, especially on the intervals of enhanced magnetic concentration and their potential paleoenvironmental implications.

Episodes of Magnetic Mineral Enrichment
Intervals (a-d) in Unit II are dated to~9.5-9.3 kyr BP, 8.5-8.4 kyr BP, 7.7-4.8 kyr BP, and 4.4-4.2 kyr BP, respectively, while the uppermost part of the core may provide a record of the last~800 yr ( Figure 4, Table 1). Intervals (a), (b), and (d) in core DS01 correspond to three short Bond-like cooling events (Bond et al., 1997;Bond et al., 2001) that are recorded in the stalagmite record from Dongge Cave, in southern China, in which more enriched δ 18 O values represent intervals of a weakened monsoon ( Figure 4L). These rapid events may be dynamically linked to climate changes in the North Atlantic and Greenland (Dykoski et al., 2005;Wang et al., 2005). The linkage may be related to the strengthening of the Siberia High, which significantly affects the East Asian monsoon climate (Chen et al., 1991;Ding et al., 1995;Liu and Ding, 1998). Intervals (a), (b), and (c) in core DS01 correspond to enhanced concentrations of both hard and soft magnetic components ( Figure 4, Table 1, and Supplementary Figure S1). These intervals are related to quite different regional environments in terms of the degree of marine influence, representing the transition from an estuarine to a deltaic environment in the head area of the Pearl River Delta (Figure 1). In contrast, interval (d), at the top of Unit II, is characterized by a higher concentration of imperfect antiferromagnetic minerals. Combined with the evaluation of the HIRM and S-ratio records (below), we infer that interval (d) represents a cold climatic event, similar to interval (b); however, the two intervals may be the product of different fluvial processes and depositional environments. The uppermost part of the core may predominantly reflect anthropogenic influences, rather than natural processes. 4.1.1 Possible Brief Climatic Events at 9. 5-9.3, 8.5-8.4, and 4.4-4

.2 Thousand Years Before Present
During the early Holocene, the entire Pearl River Delta likely experienced a rapid sea-level rise ( Figure 4H). Xiong et al. (2018) found that the rates of sea-level rise around the southeast coast of China doubled from 10.5 to 9.5 kyr BP, increasing from 16 to 33 mm/yr. Global sea-level changes, reconstructed from fluctuations in global ice volume, show a near-uniform rise during~11.4-8.4 kyr BP ( Figure 4G, Lambeck et al., 2014). Core DS01 contains a prominent peak in the concentration of magnetic minerals after 9.5 kyr BP. At this time, the head area of the Pearl River Delta may have become an intertidal environment, corresponding to the lower part of Unit II ( Figures 1B,C). Therefore, the site location of core DS01 would have been strongly influenced by the marine transgression at~9.5 kyr BP. Subsequently, rapid sediment accumulation occurred at the core site, including fine-grained sediments. Interval (a) is well-defined in most of the magnetic profiles (Supplementary Figure S1) and has a comparable level of magnetic enhancement to interval (b). In both zones, there is a pronounced increase in the clay content and a corresponding low sand content (Figures 4D,E). Evidently, these fine-grained sediments had higher concentrations of magnetic minerals than the coarse fraction. Thus, the magnetic properties and bulk sediment grain size suggest that these two intervals represent an abrupt shift in hydrodynamic forces that may have been associated with changes in precipitation and runoff in the river catchment during a period of weakened summer monsoon intensity. Within the limits of the age model, intervals (b) (8.5-8.4 kyr BP) and (d) (4.4-4.2 kyr BP) are potentially correlative with the Holocene cooling events on the centennial scale, at~8.2 and 4.2 kyr BP, respectively. The magnetic parameters (Supplementary Figure S1), bulk sediment grain size, and carbon isotope ratios (Figure 4) collectively display coherent variations corresponding to intervals (a) and (b). By contrast, interval (d) is characterized by an increase in HIRM, hence showing an increase in the hematite content. This may reflect a cold, dry climate and oxidizing conditions in the source area. In a deltaic environment with a limited marine influence, the core site was likely subjected to weak hydrodynamic sorting during this interval ( Figures 4D,F). Adjacent sediment core HKUV11 likely recorded a rapid cooling event at~4.2 kyr BP ( Figure 4I; Wu et al., 2017), while core MD05-2905 from the northern South China Sea recorded an abrupt temperature decrease at~8.2 kyr BP ( Figure 4J; Zhou et al., 2012). However, there is no corresponding cooling signal in the global sea-level record ( Figure 4G, Lambeck et al., 2014). In southernmost mainland China, the total organic matter content of the sediments of Huguangyan Maar Lake recorded cooling events at 8.2 kyr BP and 9.2 kyr BP ( Figure 4K). Several other short-term cold/dry events are also evident in proxy environmental records from Holocene lake sediments from the region, including the 4.2 kyr BP cold event (Wang et al., 2016).

Environmental Changes During 7.7-4.8 Thousand Years Before Present
The diatom-derived salinity record (Supplementary Figure S1K) indicates that during the interval of 7.7-4.8 kyr BP, following the early-middle Holocene marine transgression, the study site experienced delta progradation and the accumulation of terrestrial sediments (or a marine regression) Zong et al., 2016;Xiong et al., 2018;Fu et al., 2020). Reconstructions of the Holocene sea-level history in the Pearl River Delta show that the sea level was approximately stable during the last 7 kyr Xiong et al., 2018). The diatom profile suggests that the head area of the Pearl River Delta experienced a reduced marine influence during this interval (Supplementary Figure S1K), which may be associated with the initiation of widespread sedimentation in the Pearl River deltaic basin (Figure 4). Sediments that were sourced from tropical and subtropical areas and had undergone intensive weathering were deposited in the basin area, leading to the overall sedimentary magnetic enhancement throughout this interval.

Environmental Changes Since~800 Years Before Present
Recent agricultural and industrial activities, such as cultivation, deforestation, and mining, may have accelerated the production of magnetic minerals in Earth surface sediments. Organic carbon isotope records (δ 13 C) from at least two drilling cores (including core DS01) from the head area of the Pearl River Delta reveal a general enrichment trend in the uppermost sedimentary unit (Fu et al., 2020). Numerous previous studies have provided evidence of agricultural activity in the Pearl River Delta since at least 2.5-2.2 kyr BP (e.g., Zong et al., 2010;Yang et al., 2012;Hu et al., 2013). Although the uppermost part of the core was not sampled for magnetic measurements, there is a pronounced increase in the contents of hematite and fine-grained magnetite since~800 yr BP ( Figures 4A,C), which is approximately coeval with the mid-Southern Song Dynasty. This interval corresponds to the peak period of Medieval warmth, represented in terrestrial climatic reconstructions (Mann et al., 2008). These changes in the magnetic profiles of core DS01 may be the result of intensified human activities in the Pearl River catchment, rather than to the phase of deltaic development over the area. Based on speleothem records, Zhang et al. (2008) hypothesized a close linkage between changes in monsoonal climate and Chinese cultural phases. A rapid increase in human disturbance from~900 yr BP is evident in the As record from a sediment core from the south of Hainan Island (Wan et al., 2015). The provision of high-resolution chronologies for the late Quaternary sediments of the study area remains challenging (e.g., Yim, 1999); moreover, there may be uncertainties in correlating records of monsoon intensity and rainfall, as argued by Zhang et al. (2010). Further research is needed to precisely determine the contribution of human activities to the properties of the recent sediments of the study area.

Holocene Climatic Optimum
The period before interval (c) (~7.7-4.8 kyr BP) corresponds to the early Holocene marine transgression and is regarded as the climatic optimum in East Asia. Early studies proposed quite different timings for the Holocene optimum (defined as the peak in monsoonal precipitation or effective moisture) across the domain of the East Asian summer monsoon: at~3 kyr BP in southern China and at~9 kyr BP in northern China (see An et al., 2000 and references therein). Zhou et al. (2022) proposed a gradual northward expansion of the Holocene optimum in the East Asian monsoon region (from 24.25°N to 48.74°N), based on pollen-based precipitation reconstructions. A depleted early Holocene stalagmite δ 18 O signal is evident at Dongge Cave ( Figure 4L), which likely indicates a wetter Asian summer monsoon climate, as the δ 18 O record of Dongge Cave is anticorrelated with the Asian summer monsoon intensity (Wang et al., 2005). Records of pollen (Zhou et al., 2004) and bulk organic carbon isotope ratios (Zhong et al., 2010) from Dahu Swamp, to the north of the East River in the Pearl River drainage (Figure 1), suggest a pronounced warm and wet period during the early Holocene (~10-6 kyr BP), which may be related to a strengthened East Asian summer monsoon. Wei Z. et al. (2020) conducted sedimentary magnetic analyses at Dahu Swamp and proposed a rainfall maximum during the mid-Holocene (~8-4 kyr BP). Although the magnetic record differs from other paleoenvironmental reconstructions from Dahu Swamp, the magnetic record is in phase with an equivalent record from Daping Swamp (26°10.18′N, 110°08′E), to the north-west of Dahu Swamp (Zhong et al., 2018). Both locations may be sensitive to local hydrological changes as they are in hydrologically closed lake basins (Zhong et al., 2018;Wei Z. et al., 2020). Huguangyan Maar Lake (Figure 1) may record the earliest timing of the Holocene optimum, as indicated by an elevated total organic matter content during 11.5-6.5 kyr BP (Wang et al., 2016). A pollen record from the lake sediments indicated that during the early Holocene (11.6-7.8 kyr BP), the area was dominantly occupied by trees and shrubs, with tropical trees recording maximum percentages during~9.5-8.0 kyr BP (Wang et al., 2007). Thus, there may be minor offsets in the timing of the climatic optimum between the reconstructions from Huguangyan Maar Lake and Dahu Swamp. However, it is possible to conclude that, at least in southeastern China, the Holocene climatic optimum was not significantly timetransgressive and that the Holocene optimum undoubtedly occurred during the early Holocene. In the vicinity of the site of core DS01, a rapid marine transgression occurred during this period, which was associated with relatively low and uniform magnetic mineral concentrations, except the interval (a) (Figure 4).

Long-Term Trend of Environmental Evolution
Both magnetic grain-size parameters in core DS01 (Supplementary Figures S1I,J) show no clear response to either marine transgression or regression, when compared, for example, with the salinity profile (Supplementary Figure S1K). This is quite different from the records from core SX97 (Peng et al., 2014) or core PD , which recorded frequent oscillations during the late Quaternary or Holocene. The ARM/κ lf ratio in core DS01 indicates the long-term fining of the magnetic grain size, except for the lower part of the core (below 22 m). Neither of the neighboring cores shows comparable trends in their magnetic records Peng et al., 2014;Wu et al., 2021). Seaward of the Pearl River estuary, the magnetic grains within the surface sediments become increasingly fine, which likely resulted from the dissolution of superparamagnetic grains (Ouyang et al., 2017). An increasing distance from the source areas could also result in a fining of the magnetic grain size, as finer grains are likely to be transported longer distances than coarser grains (Huang et al., 2021). However, this is not applicable here, as indicated by the L-ratio ( Figure 5). It is unlikely that sea-level changes were responsible for the long-term fining trend of magnetic grains evident in core DS01, although there is a close relationship between sea-level changes and sedimentary magnetic grain-size variations in the South China Sea (see Kissel et al., 2020 and references therein). Neither is there clear evidence for the effect of sea-level change on East Asian monsoon intensity (e.g., Wang et al., 2001). We tentatively attribute the magnetic grain-size trend to a slightly increasing weathering intensity in the material source areas. Within the East Asian monsoon domain, the monsoonal climate (as a source of heat and moisture) would be expected to be a major factor controlling erosion and weathering processes in the Pearl River catchment. In comparison, there was a long-term trend of decreasing monsoon precipitation in southern China since the early Holocene ( Figure 4L; Dykoski et al., 2005;Wang et al., 2005), which does not necessarily conflict with increased weathering during the same period because of the complex regional hydrothermal conditions. Additionally, sediment pore water may have affected the magnetic measurements. Since water is generally diamagnetic and the sedimentary water content generally decreases down-core, ARM/κ lf could be biased accordingly. However, this influence is likely to be minor, as the pore-water content was relatively low.

SUMMARY
We have obtained a high-resolution rock magnetic record for drilling core DS01 to reconstruct the Holocene environmental evolution of the Pearl River Delta. The results indicate a relatively uniform magnetic mineral assemblage. A sedimentary hiatus occurred during~3.9-1.3 kyr BP. Enhanced magnetic mineral concentrations have occurred at the top of the core since~800 yr BP, which we attribute to the intensified anthropogenic impact. The ARM/κ lf ratio indicates a long-term trend of decreasing magnetite grain size during the Holocene, which may be a response to intensified weathering in the East Asian monsoon region. However, we cannot exclude the influence of residual pore water on this trend. A prominent interval of increased magnetic mineral concentration occurred during~7.7-4.8 kyr BP, associated with delta progradation (or marine regression) in the study area. Three abrupt shifts in magnetic properties occurred at~9.5-9.3,~8.5-8.4, and 4.4-4.2 kyr BP, which are potentially correlated with Bondlike events. The two older events are coeval with the initiation of the Holocene optimum in southeastern China, while the youngest event occurred in the context of deltaic development and minimal direct marine influence. It should be noted that possible uncertainties in the age model of core DS01 hinder the definitive correlation of these events with records from elsewhere, and further investigations of the Holocene sedimentary sequences in the Pearl River Delta and adjacent regions of East Asian monsoon influence are needed.

DATA AVAILABILITY STATEMENT
New data from this study are available from the Mendeley Data (http://dx.doi.org/10.17632/wcdhtstwtc.1).

AUTHOR CONTRIBUTIONS
YW and SF conceived the study and collected sub-samples for measurements; HX conducted the core drilling; YZ initiated the study on the Pearl River Delta; YW, YH, and JC conducted environmental magnetic analyses; TO and SP performed the rock magnetic analyses; and YW wrote the manuscript with input Frontiers in Earth Science | www.frontiersin.org June 2022 | Volume 10 | Article 882201 from SF, HX, ZZ, and TO. All co-authors contributed to discussion and approved the submitted version.