Abstract
High-resolution mineral magnetic and pollen records from overlapping piston cores from Lake Xiarinur (42°37′N, 115°28′E) document detailed changes in environment and vegetation since the last deglaciation in the Hunshandake Sandy Land, Inner Mongolia. The formation of Lake Xiarinur commenced during the Bølling–Allerød warming, as reflected by an abrupt shift in sedimentary facies from eolian sands to lacustrine sediments at a core depth of 3.5 m (~14.1 cal. kyr BP). The pollen records demonstrate that desert vegetation with sparse herbs occurred before 14.1 cal. kyr BP, which was succeeded by meadow grassland vegetation from ~14.1 –13.3 cal. kyr BP (the Bølling -Allerød). A dramatic decrease in pollen concentration occurred between ~13.3 and 11.7 cal. kyr BP corresponding to the Younger Dryas. The remarkable increases in pollen and Pediastrum concentrations at ~11.7 cal. kyr BP suggest that the environment began to ameliorate, and the most humid conditions prevailed until ~8.5 cal. kyr BP. The maxima of magnetic concentration–dependent and magnetic-grain-size-dependent parameters, together with high pollen concentrations, between ~10 and 9 cal. kyr BP are interpreted as a brief interval of high annual precipitation. Our results indicate a stepwise decrease in effective moisture in the mid-Holocene and an accelerated drying trend in the late Holocene, characteristics typical of an East Asian summer monsoon–influenced region. The broad similarities between our data and other lake records from central-eastern Inner Mongolia, well-dated speleothem δ18O records from southern and central China, and summer insolation at 45°N, support the idea that vegetation and climate changes in the Hunshandake Sandy Land were controlled mainly by fluctuations in the response of the East Asian monsoon to Northern Hemisphere summer insolation.
Keywords
Introduction
Knowledge of climatic variability during the Holocene and the last deglacial transition is important for a better understanding of the natural background variability of paleoclimate and paleoenvironmental change (e.g. An et al., 2000; Stott et al., 2004; Wang et al., 2005; Yuan et al., 2004). Climatic transition zones are particularly sensitive to climate change, and thus have been extensively used for paleoclimatic reconstructions (e.g. Clarke and Rendell, 1998; Peck et al., 2002). The semi-arid and arid areas of the Mongolian Plateau and surrounding areas are of interest for paleoenvironmental and paleoclimatic studies because the environments in these regions are sensitive to climatic change (e.g. Feng et al., 2006a, 2006b; Sun and Ding, 1998; Sun et al., 1999; Yang et al., 2013; Zhang et al., 2011). In recent decades, a number of geological records from these regions, including eolian deposits, lake sediments, and soils, have been examined to reconstruct the environmental history since the last deglaciation, particularly during the Holocene (e.g. An et al., 2000; Chen et al., 2003; Feng et al., 2006a, 2006b; Jiang et al., 2010; Shi et al., 1993; Xiao et al., 2004, 2006, 2009; Yang et al., 2011, 2013). While some studies proposed that conditions were arid during the early Holocene in the western part of arid central Asia dominated by the Westerlies (Chen et al., 2008; Feng et al., 2005; Yang et al., 2011; Yang and Scuderi, 2010), the central-eastern part of arid central Asia, in the Westerlies–monsoon transition zone, was humid during the early Holocene (Jiang et al., 2010; Lee et al., 2013; Yang et al., 2013). Controversies over both the timing and duration of the Holocene climatic optimum in this region, in addition to spatial differences in climatic variations on centennial to decadal time-scales, have hampered detailed regional comparisons and caused difficulty for the prediction of the future effects of climate change on Mongolian ecosystems. These discrepancies may be partly because of the uncertainties in proxy interpretation and dating and by spatial differences in the Asian monsoon system. Thus, more independent multi-proxy data sets are needed to improve understanding of climatic change in the monsoon–Westerlies transition zone.
The Hunshandake Sandy Land in the semi-arid zone of eastern Inner Mongolia is located at the northern margin of the East Asian summer monsoon (EASM) and at the eastern edge of the Northern Hemisphere Westerlies (NHW; Figure 1). The dune landscape of the Sandy Land is therefore jointly influenced by the East Asian monsoon (EAM) and the NHW (Jiang et al., 2010; Yang et al., 2013). Earlier studies, primarily based on the ages and stratigraphic relationships between paleosols and eolian sand layers in the eolian sequences and on geochemical and pollen-derived moisture proxies in lake sediments, reported that significant environmental change has occurred in this region (e.g. Jiang et al., 2010; Xiao et al., 2008; Yang et al., 2013). It is therefore important to perform multi-proxy paleoenvironmental studies at different sites in order to provide a comprehensive understanding of climatic patterns in central-eastern Inner Mongolia.

Topography and atmospheric circulation systems of China. Arrows indicate the wind directions of the East Asian summer monsoon (red), the East Asian winter monsoon (blue), the Indian monsoon (orange), and the mid-latitude Westerlies (lilac, dashed). Crosses mark the locations of Lake Xiarinur (red) and other sites referred to in the text (yellow): Lake Bayancagan (41°39′N, 115°13′E), Lake Angulinur (41°18′N, 114°20′E), Lake Dalinur (43°13–23′N, 116°29–45′E), Lake Daihai (40°35′N, 112°41′E), Lake Hulun (49°8′N, 117°30′E), Lake Ulaan (44°30′N, 103°39′E), Sanbao Cave (31°40′N, 110°26′E), and Dongge Cave (25°17′N, 108°5′E). The inset shows the bathymetry of Lake Xiarinur and the location of the coring sites (solid red circle).
In previous studies, a magnetic approach was seldom used to reconstruct short-term environmental fluctuations from different archives developed in the Hunshandake Sandy Land. By measuring variations in the composition, concentration, and grain size of magnetic minerals, mineral magnetism has been widely used to reconstruct paleoclimatic and paleoenvironmental processes recorded in marine sediments (Robinson, 1986; Rousse et al., 2006), lake sediments (Chen et al., 1999a, 1999b; Thouveny et al., 1994; Wang et al., 2008), and eolian sediments (Fang et al., 1999; Heller and Liu, 1984; Kukla et al., 1988). Of these natural archives, lake sediments have been widely used for studying environmental processes on a wide range of spatial and temporal scales (Li et al., 2006; Thompson and Oldfield, 1986; Verosub and Roberts, 1995). This is partly because lake sediments can provide data with higher temporal resolution than marine sediments, and they also offer abundant information on diverse aspects of environmental change, such as erosion history and soil-forming processes, in addition to the climatic history of continental regions and regional responses to large-scale climate change (Dearing et al., 2001; Oldfield et al., 1979; Rosenbaum and Reynolds, 2004).
Lake Xiarinur, located in the heart of the Hunshandake Sandy Land, is likely to be sensitive to the intensity and fine-scale shifts of the climatic drivers in eastern Inner Mongolia. Here, we present combined records of high-resolution magnetic, pollen, and particle size analyses of sediments from Lake Xiarinur to improve understanding of paleomonsoon variability since the last deglaciation of this region. The main aims of the study are as follows: (1) to assess the potential of mineral magnetic properties for reconstructing paleoenvironmental changes in the semi-arid Hunshandake Sandy Land, (2) to understand the processes of formation and evolution of the lake sediments using combined pollen and magneto-climatological analyses, and (3) to evaluate the variable contributions of the Holocene EAM and NHW in this climatically sensitive region.
Study site
Lake Xiarinur (42°37′N, 115°28′E) is a small closed saline lake, situated in the central area of the Hunshandake Sandy Land in central-eastern Inner Mongolia, ~260 km north of Beijing (Figure 1). The lake has a surface area of 3.1 km2 and a maximum water depth of 1–2 m. The catchment is surrounded by semi-fixed dunes of ~2–20 m in height, on which the vegetation is dominated by meadow-steppe species with small forest patches surrounding the catchment area (Sun et al., 2007; Zheng et al., 2006). The mean annual temperature and precipitation in the area are ~1.5°C and ~365 mm, respectively. About 80–90% of the precipitation occurs between July and September, a distribution that is typical of regions influenced by the EAM. From the middle of November to early April, the ground is covered with ice. During winter, the Mongolian High generates the prevailing cold, dry, northwesterly airflow which brings strong winds and cold air to the region. During summer, warm and moist air masses are driven northwards by the pressure gradient between the Subtropical High and the Continental Low; they interact with cold air from the northwest, producing the majority of the annual precipitation.
Sampling and methods
Sampling
A suite of overlapping piston cores has been retrieved from Lake Xiarinur since 2007 (Figure 2). Core XN-2007 did not reach the lake basement, while core XRD-2012 reached eolian sands at the bottom of the lake sediment sequence, at depths of <4 m. In this study, we focused on the upper 5-m-long composite section of XRD-2012 that represents the complete sedimentary succession. The cores were split into half longitudinally; one half was sampled at 1-cm intervals for determination of water content and dry density and for pollen, and grain-size analyses, while the other half was used for mineral magnetic analysis. A total of 240 magnetic samples were extracted from the split face of the core sections using a continuous sequence of cubic plastic boxes (external dimensions 2 cm × 2 cm × 2 cm with an internal volume of 7 cm3). All samples were stored at 4°C. Radiocarbon ages were determined by accelerator mass spectrometry in the Poznań Radiocarbon Laboratory, Poland, and Peking University, China. All AMS 14C ages were calibrated using the program CALIB 4.3 (Stuiver et al., 1998).

Lithological column and age–depth model for Lake Xiarinur. The black crosses and blue dots represent calibrated AMS 14C ages of cores XN-2007 and XRD-2012, respectively. The red dots denote 137Cs-derived ages of core XRD-2012.
The upper 12 cm was sampled at a 0.5-cm interval for determination of 137Cs and 210Pb activity (Figure 3). Activity measurements of 137Cs, 210Pb, and 226Ra were carried out by gamma spectrometry using a low-background well-type germanium detector (EGPC 100P-15R). Each sample was packed in a 15-mm polyethylene tube for 3 weeks of storage in sealed containers to allow radioactive equilibration. Total 210Pb activity was determined by gamma spectrometry via its energy at 46.5 keV. The short-lived 226Ra daughter nuclides 214Pb (241.9, 295.2, and 351.9 keV) and 214Bi (609.3 keV) were measured to determine the supported 210Pb content in order to calculate the unsupported 210Pb content. 137Cs was measured by its emissions at 662 keV.

137Cs, 210Pb, and 226Ra activities versus depth for core XRD-2012.
Mineral magnetic analyses
The magnetic susceptibility (χ, normalized by mass) was determined for all samples using an AGICO MFK1-FA Kappabridge magnetic susceptibility meter at frequencies of 976, 3904, and 15,616 Hz (with sensitivities of 2 × 10−8, 6 × 10−8, and 12 × 10−8 SI, respectively) in a peak magnetic field of 200 A/m. Here, χfd is defined as χ976 Hz – χ15,616 Hz. Natural Remanent Magnetization (NRM), Anhysteretic Remanent Magnetization (ARM), and Saturation Isothermal Remanent Magnetization (SIRM) of all samples were measured using a 2G Enterprises Model 755-4K superconducting rock magnetometer housed in a magnetically shielded room. ARM was imparted in a 2G alternating field demagnetizer in a DC field of 0.05 mT and a maximum AC field of 100 mT; it is expressed as an anhysteretic susceptibility (χARM) by dividing by the strength of the DC field. SIRM was imparted in a DC field of 1 T using an ASC Model IM-10-30 pulse magnetizer. Back-field IRM was imparted in a DC field of 0.3 T to obtain the ratio S−0.3 T (−IRM−0.3 T/SIRM).
A total of 19 representative samples were selected for hysteresis and first-order reversal curve (FORC) measurements using a MicroMag 3900 AGM. The hysteresis loops were corrected for high field magnetic susceptibility, enabling the calculation of the ratio of saturation remanent magnetization to saturation magnetization (Mrs/Ms) and coercivity (Bc). DC demagnetization fields were used to determine the coercivity of remanence (Bcr). A total of 120 reversal curves were measured with an averaging time of 0.2 s per data point, and a smoothing factor of 5 was applied to suppress noise in the contours of the FORC diagram.
Pollen and particle size analyses
Pollen analysis was carried out on 129 samples taken from the core at intervals of 2–4 cm. Lycopodium spore tablets were added to samples to calculate pollen concentrations (Moore et al., 1991). Sediment samples (1–3 g) were treated with 10% HCl, 10% KOH, and 36% HF to remove carbonate, humid acid, and silicate, respectively, before mounting the fossil pollen samples on microscopic slides using water-free glycerol. The treated samples were examined using an Olympus BX51 transmitted light microscope. The remains of Pediastrum were also counted. Particle size analysis was performed on 165 samples taken from the core at 2-cm intervals using a Malvern Mastersizer 2000 laser particle size analyzer after the carbonate-free samples had been ultrasonically treated in a 10% (NaPO3)6 solution.
Results
Lithology and chronology
The cores XN-2007 and XRD-2012 are well correlated, except for obvious compression in the upper 90 cm of core XRD-2012 because of the freezing of the uppermost sediment (Figure 2). Freezing of the uppermost part of the sediment column is the result of shallowing of the lake (the water depth was only 0.2 m at the sampling site) caused by the extraction of ground water for agriculture in recent years. The lithological sequence can be divided into two parts: a lower interval (>3.5-m depth) consisting of homogeneous coarse-grained eolian sands, and an upper interval of lake sediments (0–3.5 m) consisting mainly of gray, silty marls intercalated with dark, clayey silts and/or brownish-red clay (Figure 2). The change from eolian sand to lake sediments at a depth of 3.5 m, in addition to a sudden increase in water content and a corresponding decrease in dry density, marks the formation of the lake at that level.
The 137Cs, 210Pb, and 226Ra activities versus depth in core XRD-2012 are shown in Figure 3. An abrupt increase in 137Cs concentration at 5-cm depth is associated with the abrupt increase in nuclear testing in 1951/1952 (Crusius and Anderson, 1995), and provides satisfactory age control for the uppermost sediments. The absence of a prominent 137Cs peak in the upper 5 cm may be attributed to disturbance of the surface sediments either during coring or by human activity. 210Pb dating cannot be used because of the low levels of unsupported 210Pb in the sediments, coarse-grained sediments, and large variations in accumulation rates, as suggested by Lamoureux and Gilbert (2004). The final age–depth model was constructed using linear interpolation of five AMS 14C dates from terrestrial plant remains (Table 1) together with the 137Cs-derived ages for the uppermost sediments (Figure 2).
Radiocarbon ages from Lake Xiarinur.
All AMS 14C ages were calibrated using the atmospheric data set from the calibration program CALIB 4.3 (Stuiver et al., 1998). The AMS 14C age at the depth of 20 cm is too old and was not used in the age model.
Pollen assemblage zones
A total of 50 terrestrial fossil pollen types were identified, including 4 coniferous trees, 16 broadleaf trees, and 29 shrubs and herbs. In addition Pediastrum, one of the genera of green algae that are most abundant in very shallow lakes, was identified. Using the concentrations of the major pollen taxa and the affinity scores of meadow-steppe, broadleaf, and coniferous trees, the composite core can be divided into six pollen assemblage zones, as follows (Figure 4).

Pollen concentration diagram for the Lake Xiarinur sediment core, including sediment dates. Six pollen assemblage zones are indicated to the right.
Zone 1 (3.9–3.5 m; ~15.3–14.1 cal. kyr BP)
The pollen concentration is relatively low in this zone (138–7,315 grains/g). The pollen assemblages are dominated by herbs, with Artemisia (55.2–72.6%), Ranunculaceae (1.28–13.04%), and Chenopodiaceae (2.2–14%) being the dominant taxa (Figure 4). The consistently low percentages and concentrations of pollen from coniferous trees, and from plants from aquatic and wetland environments, suggest relatively dry and cold conditions (Figure 4). Pediastrum was not detected, indicating that the lake did not exist at this time.
Zone 2 (3.5–3.2 m; ~14.1–13.3 cal. kyr BP)
The pollen concentrations are significantly higher in this zone (23,460–164,214 grains/g). Shrubs and herbs (Artemisia, Chenopodiaceae, Ranunculaceae, and Gramineae) are the dominant taxa, comprising 77.9–90% of total pollen, with Artemisia as the major contributor. Broadleaf trees comprise the majority of the remaining total pollen (8–20.3%), including Betula (3.4–14.1%), Carpinus (0.4–2.2%), Quercus (<1.7%), and Ulmus (<1%). Conifer pollen occurs at relatively low levels. Percentages of wetland herbs (Cyperaceae <4.7%) and aquatic Typha are significantly higher than in Zone 1. Pteridophyte spores are absent. Meadow grassland is the inferred vegetation type. Pediastrum is present in this zone at ~14.1 cal. kyr BP, indicating the formation of the lake.
Zone 3 (3.2–2.8 m; ~13.3–11.7 cal. kyr BP)
Zone 3 is characterized by a dramatic decrease in pollen concentrations (12,264–94,697 grains/g). Shrubs and herbs comprise 80–88.5% of the total pollen, with Artemisia as the dominant taxon (56.1–66.9%). Chenopodiaceae (4.4–9.3%), Ranunculaceae (4.7–9%), and Gramineae (2.0–8.2%) are significantly represented. Broadleaf trees comprise 8.7–18.5% of total pollen, with Betula (5.1–10.1%) the most significant component. Conifer pollen (e.g. Pinus, Abies, and Picea) is present at a consistently low level (0.7%) and probably represent long-distance transport. Percentages of aquatic and wetland pollen types are much lower (0.5–2.8%) than in Zone 2; Cyperaceae (0.5–2.5%) is the major component. Pediastrum occurs intermittently. Meadow grassland continues to dominate in this zone; however, the decreased pollen concentration indicates a trend toward sparser vegetation coverage.
Zone 4 (2.8–1.6 m; ~11.7–8.5 cal. kyr BP)
This zone is marked by a significant rise in total pollen concentration (76,846–486,550 grains/g). Although shrubs and herbs are dominant (76.7–90.6%), Artemisia percentages (31.2–60.3%) are much lower than in Zone 3. The representation of Chenopodiaceae (4.9–16.8%) and Gramineae (2.2–9.6%) increases steadily. Broadleaf trees (6.1–17.7%), coniferous trees (0.2–3.3%), and aquatic and wetland taxa (2.1–8.5%) are represented. Pediastrum reaches a peak at 2.3-m depth (~10.3 cal. kyr BP), indicating the highest lake level recorded during the interval of study. The high pollen concentrations and assemblages indicate a continuation of herb- and shrub-rich meadow grassland but with much denser vegetation coverage under an inferred relatively warm and moist climate.
Zone 5 (1.6–1.0 m; ~8.5–4.2 cal. kyr BP)
The pollen concentration in Zone 5 (72,630–513,375 grains/g) is generally lower than in Zone 4. Artemisia comprises 45.8–57.5% of total pollen, while Gramineae and Chenopodiaceae contribute 1.3–7.5% and 2.6–16%, respectively. The percentage of Betula increases slightly, while the percentages of both Ulmus (0.8–4.6%) and Carpinus (2.1–4.5%) decrease relative to Zone 4. The maximum total pollen and Pediastrum concentration occur at ~5.7 cal. kyr BP, reflecting a short interval of optimum climatic conditions.
Zone 6 (1–0 m; ~4.2 cal. kyr BP–present)
The most significant feature of Zone 6 is the sharp decrease in total pollen concentration (1,216–72,983 grains/g) from ~4.2 cal. kyr BP. Shrubs and herbs comprise a relatively low proportion of the total pollen (32.4–49.9%), while the percentages of conifers (10.2–70.4%) increased, and those of broadleaf trees (48–6,243 grains/g) declined. Pollen concentrations of wetland taxa, including Cyperaceae and Typha, are significantly lower than in Zone 5, although they still have relatively high percentages. Pediastrum exhibits a decreasing trend in this zone, indicating the continued shrinkage of the lake.
Mineral magnetic properties
The combination of the Day plot and the FORC diagram provides a powerful tool to identify both the coercivity and grain-size distribution of magnetic materials (Day et al., 1977; Roberts et al., 2000). The hysteresis loops are all closed below 500 mT, suggesting the predominance of low-coercivity ferrimagnetic minerals (Figure 5a and b). The average values (n = 19) of Bc and Bcr are 13.2 mT and 30.6 mT, respectively, both of which are close to the theoretical values of single-domain (SD) magnetite at 10 mT and 33 mT, respectively. In all FORC diagrams, the contours along the Bu axis represent the axis of the interaction field among magnetic grains between 5 and 10 mT (Figure 5c and d), indicating the dominance of SD magnetite in the samples (Muxworthy and Dunlop, 2002; Roberts et al., 2000). Compared with larger interacting grains, weakly or noninteracting uniaxial single-domain (UNISD) magnetite particles are characterized by relatively closed contours along the Bc axis and a narrow spread along the Bu axis (Egli et al., 2010; Muxworthy and Roberts, 2007). The FORC diagrams exhibit a narrow (central) ridge along the Bc axis, which is characteristic of UNISD particles (Figure 5c and d). Progressive acquisition of IRM in magnetic fields of <1 T shows that for both high- and low-χ samples, saturation is achieved at ~300 mT, consistent with the presence of a dominant ferrimagnetic mineral phase (Figure 5e). The hysteresis parameters (Bcr/Bc and Mrs/Ms), calculated after removal of the paramagnetic component of the magnetization, are plotted on a Day plot (Day et al., 1977) in Figure 5f. All samples fall in the pseudo single-domain (PSD) range regardless of the differences in the two ratios. Both S−0.3 T and IRM@100 mT/SIRM values reflect the relative proportion of anti-ferromagnetic minerals in the magnetic assemblage (Bloemendal et al., 1988; Hesse, 1997; Robinson, 1986). Values of S−0.3 T are close to 1 when virtually no high-coercivity minerals are present, while lower values indicate a higher content of high-coercivity minerals (Bloemendal et al., 1988). The S−0.3 T values are generally >0.85 throughout the core, implying the significant presence of low-coercivity minerals.

Mineral magnetic parameters for the Lake Xiarinur sediment core. (a and b) Hysteresis loops, corrected for paramagnetism, and (c and d) FORC diagrams for two representative samples. (e) IRM acquisition curves, and (f) hysteresis ratios plotted on a Day diagram (Day et al., 1977).
Magnetic susceptibility (χ) commonly reflects the concentration of ferrimagnetic minerals in sediments (Thompson and Oldfield, 1986; Verosub and Roberts, 1995). ARM is notably sensitive to small SD particles, while SIRM indicates relatively coarse magnetic particles (Evans and Heller, 2003). The largely concordant variations among these magnetic concentration–dependent parameters throughout the core suggest that they are controlled mainly by variations in the concentration of magnetic minerals (Figure 6). The mass-specific χfd is a reliable proxy for the concentration of viscous SP particles (~20–25 nm for magnetite and maghemite; Liu et al., 2005; Zhou et al., 1990). χARM/χ and χARM/SIRM are widely employed to detect variations in magnetic mineral grain size (e.g. Evans and Heller, 2003; Maher and Thompson, 1991; Zhou et al., 1990). Figure 6 shows that higher χARM/χ and χARM/SIRM values generally correspond to higher χ, χfd, and ARM values, implying that the magnetic mineral grain sizes are generally finer in the intervals with higher concentrations of magnetic minerals. The ARM/SIRM ratio is also diagnostic of bacterial magnetosomes, and ARM/SIRM ratios between 0.15 and 0.25 are characteristic of intact magnetosomes (Moskovitz et al., 1993). ARM/SIRM values throughout the core are between 0.001 and 0.064, consistently lower than 0.15, suggesting a minimal contribution from bacterial magnetite.

Plots of mineral magnetic parameters, water content, dry density, and grain-size distribution against age. (a) Mass-specified magnetic susceptibility (χ), (b) frequency-dependent susceptibility (χfd), (c) saturation isothermal remanent magnetization (SIRM), (d) anhysteretic remanent magnetization (ARM), (e) χARM/χ, (f) χARM/SIRM, (g) IRM@100 mT/SIRM, (h) S−0.3 T (−IRM−0.3 T/SIRM), (i) water content, (j) dry density, and (k) grain-size fraction > 35 µm. The blue bar denotes the Younger Dryas (YD) event, and the yellow bar indicates the most humid interval in the early Holocene. Here, χfd is defined as (χ976 Hz - χ15,616 Hz)/χ976 Hz. IRM@100 mT is the remanence after alternating field (AF) demagnetization of an IRM imparted in a 1-T field with a peak AF of 100 mT.
Discussion
Formation and evolution of Lake Xiarinur
Accurate AMS 14C dating of the formation and evolution of Lake Xiarinur is required for paleoclimatic reconstruction. In sandy land areas, it is difficult to obtain suitable material for macrofossil AMS 14C dating. In this study, macrofossils (e.g. the stems from shrubs) found in the sediment provide a reliable time control for the sedimentary sequence. Moreover, the core ages project to the sediment surface at an age that is very close to zero (Figure 2), also implying that the carbon reservoir effect is not significant in Lake Xiarinur. As shown in Figure 2, the sedimentary facies changes from eolian sands to lacustrine sediments at the depth of 3.5 m (~14.1 cal. kyr BP). A sharp increase in water content and pollen concentration at 3.5-m depth indicates that the lake formed at this time. Lake formation may have occurred in response to the Bølling–Allerød warming (Denton et al., 2010; Kienast et al., 2003), a period characterized by a rapid rise in temperature following the Last Glacial Maximum. Subsequently, Lake Xiarinur recorded the following additional climatic stages as inferred from the pollen assemblages, variations in lithology, sediment dry density, and mineral magnetic characteristics: Stage 1 (~13.3–11.7 cal. kyr BP), Stage 2 (~11.7–8.5 cal. kyr BP), Stage 3 (~8.5–4.2 cal. kyr BP), and Stage 4 (~4.2 cal. kyr BP–present), corresponding to Pollen Zones 3–6, respectively. These stages may correspond to the Younger Dryas (YD) and the early, middle, and late Holocene, respectively (Figures 4 and 5). During the YD, the vegetation was dominated by herbs and shrubs, including mainly Artemisia and Ranunculaceae, indicating relatively cold and dry conditions (Figure 4). This time interval is also characterized by a significant decrease in both pollen concentration and water content and by an increase in bulk and magnetic mineral grain sizes, as reflected by the grain-size fraction >35 µm (%) and the χARM/SIRM ratio (Figure 6). We propose that this cold stage may represent a regional response to the YD event (Fairbanks, 1989; Lowe et al., 2008), during which Greenland Summit was 15°C colder, as inferred from thermally fractionated nitrogen and argon isotope ratios from the Greenland Ice Sheet Project Two (GISP2; Alley et al., 1993; Dansgaard et al., 1989). The period between ~11.7 and 8.5 cal. kyr BP was characterized by dramatic increases in both pollen concentration and water content, suggesting a dense vegetation cover (Figure 4). There was a remarkable increase in Pediastrum assemblages during the early Holocene, between 11.7 and 10.2 cal. kyr BP. Since Pediastrum, a genus of green algae, favors freshwater and has been regarded as a proxy for lake level change (Alhonen and Ristiluoma, 1973; Jiang et al., 2010; Komárek and Jankovská, 2001), the significant increase in its concentration clearly indicates a rapid increase in lake level at the onset of the Holocene. Total pollen concentrations began to decrease gradually between ~8.5 and 4.2 cal. kyr BP but still remained at relatively high levels. A significant increase in the concentrations of Pediastrum and various pollen taxa between 5.8 and 4.9 cal. kyr BP probably suggests a rise in lake level during the mid-Holocene (Figure 4). After ~4.2 cal. kyr BP, Pinus and other coniferous trees, which can survive in an arid environment, began to replace other broadleaf forest taxa. The significant reduction in the total pollen concentration during this period clearly indicates a decline in vegetation cover and general environmental degradation, although the concentrations of Pinus and other coniferous taxa were still relatively high (Figure 4). The ‘4.2 ka event’, which is characterized by cold and dry conditions, has been identified globally as a period of severe and abrupt climate change (Bond et al., 2001; Cullen et al., 2000; Stott et al., 2004; Thompson et al., 2002). At ~1 cal. kyr BP, a slight increase in the concentration of the pollen of Chenopodiaceae and aquatic plants implies a slight amelioration of the climate. In summary, the biostratigraphic data from Lake Xiarinur suggest humid conditions in the early and mid-Holocene, and a dry climate in the late Holocene. Our results provide a clear indication of the regional effects of global climate change on vegetation variations in the Hunshandake Sandy Land since the last deglaciation.
Precipitation-controlled magnetic records at Lake Xiarinur
Both wind-blown dust and a combination of sediment derived from surface runoff, rain splash, and remobilization of marginal lake sediments are potential sources of detrital material entering Lake Xiarinur. Fluctuations of the mineral magnetic properties of the lake sediments may, to some extent, reflect varying contributions from these two sources of sediment from within and beyond the catchment area (Thompson and Oldfield, 1986). Considering the location of Lake Xiarinur in relation to dust-transporting wind systems (Figure 1), both the East Asian winter monsoon and the Westerlies may have promoted inputs of eolian dust from the arid Gobi and other deserts in northwestern China when northwesterly/westerly surface winds predominated. These sandy materials may have also been transported to the lake either by strong winds during cold–dry periods, or by intensified rain splash/wave erosion under warm–humid conditions, particularly during the Holocene climatic optimum.
The vegetation changes at Lake Xiarinur since the last deglaciation indicate that maximum moisture levels occurred during the early Holocene, as indicated by the highest pollen concentrations, including abundant broadleaf trees, and of Pediastrum (Figure 4). During this period, the marked increase in annual precipitation resulted in enhanced vegetation coverage, which probably stabilized the dunes and enabled the formation of pedogenically produced soil-like material (Yang et al., 2013; Zhu et al., 1980). The occurrence of brownish-red clay in the depth interval 2.25–2.1 m (corresponding to ~10.2–9.9 cal. kyr BP) confirms that abundant and regular precipitation induced strong pedogenesis.
All magnetic concentration–dependent parameters (e.g. χ, ARM, and SIRM) show stepwise increases at the onset of the early Holocene, reaching peak values between ~10 and 9 cal. kyr BP (Figure 6). Higher χfd and χARM/χ values during this period (Figure 6) indicate the presence of abundant fine-grained ferrimagnetic particles (Evans and Heller, 2003; Maher and Thompson, 1991; Zhou et al., 1990). This is consistent with the prevailing view that pedogenic processes result in the formation of fine-grained ferrimagnetic minerals and, therefore, significant magnetic enhancement (Banerjee et al., 1993; Maher and Thompson, 1991; Zhou et al., 1990). Previous studies have shown that ferrimagnetic minerals are usually relatively abundant in soils, particularly in the topsoil, compared with wind-blown dust and other materials originating from dunes (Oldfield et al., 1985; Peck et al., 2004). We propose that rain splash or wave erosion associated with increased rainfall and increasing lake levels eroded material surrounding the lake, causing increased detrital inputs. This eroded material probably contained higher concentrations of partly pedogenized loess-derived magnetic minerals with higher χ values, consistent with the magnetic properties of the sediment cores.
The period between ~8.5 and 4.2 cal. kyr BP is characterized by low χ, ARM, and SIRM values, implying low concentrations of ferrimagnetic minerals. Both χARM/χ and χARM/SIRM ratios show stepwise decreasing trends, indicative of a coarsening of the ferrimagnetic grains. A significant increase in IRM@100 mT/SIRM and a corresponding decrease in S−0.3 T suggest decreasing contributions of low-coercivity ferrimagnetic grains. We attribute this to decreasing rates of rain splash/wave erosion and consequently declining surface-erosion-related input of ferrimagnetic minerals to the lake.
The period from ~4.2 cal. kyr BP to present is characterized by gradually decreasing trends of χ, ARM, and SIRM values, implying further climatic deterioration. Lower χARM/χ and χfd values, compared with the early and middle Holocene, reflect further coarsening of the ferrimagnetic particles (Figure 6). Lower S−0.3 T (Figure 6h), but higher IRM@100 mT/SIRM values (Figure 6g), suggest a further decrease in the percentages of ferrimagnetic minerals. χARM/χ, χARM/SIRM, and IRM@100 mT/SIRM, together with pollen concentrations, show comparable high-amplitude variations between ~2 and 0.7 cal. kyr BP, reflecting a short time interval inferred to have been wet and warm (Figure 6).
In summary, pollen concentrations, the magnetic-grain-size-dependent parameter χARM/χ, and the magnetic-composition-dependent parameter IRM@100 mT/SIRM exhibit comparable trends of variation (Figure 7). We suggest that the enhanced magnetic properties of the Lake Xiarinur sediment cores are probably controlled by precipitation delivered by the Asian summer monsoon that resulted in detrital inputs dominated by pedogenic ferrimagnetic grains, and to a lesser extent by a varying influx of eolian dust carried by the northwesterly/westerly winds of the winter monsoon. Both the pollen and mineral magnetic data suggest high precipitation and surface erosion during the early Holocene, decreased effective moisture in the mid-Holocene, and a clear trend to arid conditions in the late Holocene.

Comparison of sedimentary records from Lake Xiarinur with other lake records from Inner Mongolia, well-dated speleothem δ18O records from southern and central China, and summer insolation at 45°N. (a) Pollen concentrations from the Lake Xiarinur sediment core, (b) IRM@100 mT/SIRM ratios of the Lake Xiarinur sediments, (c) χARM/χ ratios of the Lake Xiarinur sediments, (d) annual precipitation (Pann) at Lake Bayancagan (mm; Jiang et al., 2010), (e) total organic carbon (TOC; %) of Lake Ulaan sediments (Lee et al., 2013), (f) speleothem δ18O values (‰) from Sanbao Cave (red; Dong et al., 2010), (g) speleothem δ18O values (‰) from Dongge Cave (green; Dykoski et al., 2005; Wang et al., 2005), and (h) summer insolation at 45°N (W/m2; Berger, 1978). The dark gray, pink, blue, red, yellow, and light gray bars correspond to the Mystery Interval (Denton et al., 2006), the Bølling–Allerød warming, the Younger Dryas (YD) event, the early Holocene, the mid-Holocene, and the late Holocene, respectively.
Influence of the EASM on the climate of central-eastern Inner Mongolia
Effective moisture determines the areal extent of deserts in arid regions, the water level of lakes, and the density of vegetation cover, and therefore has a strong influence on the ecological environment (Lee et al., 2013). The effective moisture in arid Central Asia is dominated by the Westerlies and was lowest in the early Holocene and highest in the mid-Holocene (Feng et al., 2005; Xiao et al., 2004; Yang and Scuderi, 2010; Yang et al., 2011; Zhang et al., 2012; Zhao and Yu, 2012). However, in EASM-dominated regions, effective moisture was strongest in the early and mid-Holocene, and weakened after the mid-Holocene (Chen et al., 2008; Jiang et al., 2010; Lee et al., 2013; Yang et al., 2013; Zhang et al., 2011). Feng et al. (2006a) comprehensively summarized the evidence indicating that the Inner Mongolian Plateau was warm and dry during the early Holocene (from ~10.5–9.5 to ~8–7.5 cal. kyr BP), and was warm and wet from ~7.5 to ~3.5 cal. kyr BP. Using multi-proxy analyses of sediments from Lake Angulinur, Wang et al. (2010) proposed that the environment became more benign at ~10.9 cal. kyr BP, with relatively wet and warm conditions during ~10.9–8.9 cal. kyr BP, followed by environmental deterioration from ~7.4 cal. kyr BP; the driest and coolest conditions occurred during ~2.2–0.48 cal. kyr BP. Based on high-resolution, well-dated paleo-records (pollen, Pediastrum, and δ18O of authigenic carbonate) from Lake Bayancagan, Jiang et al. (2006, 2010) concluded that the most humid conditions occurred between ~10.5 and ~6.5 cal. kyr BP, and that the mid-Holocene cooling began at ~8 cal. kyr BP, with a drying phase commencing ~1.5 kyr later. At Lake Ulaan, to the east of the southeastern Gobi Altai Mountains, records of total organic carbon (TOC), C/N ratios, and weathering intensity also indicated that the effective moisture was highest during the early Holocene, and was followed by a decreasing trend, reaching the lowest values during the late Holocene (Lee et al., 2013). Lake Dalinur, located 120 km northeast of Lake Xiarinur, experienced the highest lake level during the early Holocene (~11.5–7.6 cal. kyr BP); dramatic fluctuations in lake level during the mid-Holocene (~7.6–3.4 cal. kyr BP) with three intervals of lower lake stands at ~6.6–5.8, ~5.1–4.8, and ~4.45–3.75 cal. kyr BP, and a general shrinking trend of the lake during the late Holocene (~3.45 cal. kyr BP to present; Xiao et al., 2009). The most recent study of eolian sequences in the Hunshandake Sandy Land clearly indicates the formation of paleosols at ~9.6–3 cal. kyr BP, confirming that the dune landscape was stabilized by vegetation to a large degree during the early and middle Holocene (Yang et al., 2013). Syntheses of moisture changes in monsoon-influenced regions also demonstrate that a humid climate dominated during the early and middle Holocene, and that a drier climate prevailed during the late Holocene, with an abrupt shift at ~4.5 cal. kyr BP (Herzschuh, 2006; Zhao et al., 2009).
The persistent deterioration of the climate at Lake Xiarinur since ~4.2 cal. kyr BP is also consistent with several climate records from central-eastern Inner Mongolia and adjacent regions. At Lake Dalinur, low organic carbon and high inorganic carbon concentrations during ~4.45–3.75 cal. kyr BP suggest a major interval of decreased monsoonal precipitation (Xiao et al., 2009). At Lake Hulun, a dramatic fall in tree and shrub pollen percentages indicates an extremely dry event at ~4.4–4.2 cal. kyr BP (Wen et al., 2010a, 2010b). This event is also recorded at Lake Daihai, as indicated by the decline in woody plants at ~4.45–2.9 cal. kyr BP (Xiao et al., 2004). The upper boundary of a dark paleosol layer at a well-exposed section at Saihanxili, in the western part of the Hunshandake Sandy Land, was dated to ~4.2 cal. kyr BP, indicating the onset of drier conditions in the late Holocene (Yang et al., 2013).
The location of Lake Xiarinur at the margin of the EASM and the mid-latitude Westerlies implies that both of these atmospheric circulation systems probably played key roles in determining the course of environmental change in this area. Figure 7 shows that the pollen and magnetic records from Lake Xiarinur, annual precipitation reconstructions at Lake Bayancagan, and TOC (%) records of Lake Ulaan, all exhibit very similar trends of variation during the Holocene, typical of a summer-monsoon-influenced domain and in broad agreement with the synthesized EAM strength (Wang and Feng, 2013). A striking resemblance between pollen concentrations from Lake Xiarinur and two absolute-dated stalagmite δ18O records from Sanbao Cave (Dong et al., 2010) and Dongge Cave (Dykoski et al., 2005; Hu et al., 2008; Wang et al., 2005), both of which are characterized by rapid decreases in δ18O values during the early Holocene that gradually increase toward the present, further confirms that the lake records from central-eastern Inner Mongolia provide a consistent pattern of Asian summer monsoon precipitation-controlled Holocene climate change (Figure 7). This conclusion can be reached regardless of the controversies over the contributions of EASM and Indian summer monsoon to speleothem oxygen isotope records from southern and central China (Dong et al., 2010; Dykoski et al., 2005; Wang et al., 2005). Comparison of the lake proxy records with summer insolation at 45°N (Figure 7) supports the idea that vegetation patterns and regional climate changes in central-eastern Inner Mongolia during the Holocene were controlled by fluctuations in the response of the EASM to Northern Hemisphere summer insolation, which acts as a major driver of global temperature change (Kutzbach and Gallimore, 1988).
Conclusion
A high-resolution, radiocarbon and 137Cs dated, pollen and mineral magnetic record from Lake Xiarinur in the Hunshandake Sandy Land, central-eastern Inner Mongolia, reveals a detailed sequence of environmental changes during the last deglaciation and the Holocene. A change from eolian sand deposits to a lacustrine sequence at ~14.1 cal. kyr BP records the formation of Lake Xiarinur during the Bølling–Allerød warming, in response to rapid Northern Hemisphere deglaciation. The pollen data demonstrate that conditions began to ameliorate at ~11.7 cal. kyr BP, marking the end of the YD stadial, and maximum moisture levels were reached during the early Holocene. A gradual decrease in effective moisture occurred during the mid-Holocene (~8.5–4.2 cal. kyr BP), and an environmental deterioration toward increased aridity commenced at ~4.2 cal. kyr BP. Increases in the concentration of Pediastrum at ~11.7–10.2 and ~5.8–4.9 cal. kyr BP are inferred to represent two wet intervals with substantial rises in lake level. Mineral magnetic analysis reveals that the highest concentrations of magnetic minerals occurred at ~10–9 cal. kyr BP, which is interpreted as reflecting high detrital input of loess-derived magnetic minerals with a significant component of ferrimagnetic minerals of pedogenic origin by either rain splash or wave erosion associated with increased rainfall and/or a rise in lake level. The subsequent decline in the concentration of magnetic minerals is thus attributed to declining surface-erosion-related input of ferrimagnetic minerals into the lake. The stepwise decline in effective humidity through the Holocene, as inferred from parallel vegetation and magneto-climatological variations at Lake Xiarinur, is broadly consistent with other lake records from central-eastern Inner Mongolia, well-dated speleothem δ18O records from southern and central China, and summer insolation at 45°N. We therefore propose that regional climatic changes and associated variations in the density and composition of vegetation cover in the Hunshandake Sandy Land during the last deglaciation and throughout the Holocene were controlled mainly by fluctuations in the response of the EASM to Northern Hemisphere summer insolation.
Footnotes
Acknowledgements
We sincerely thank Professor Zhaoyan Gu, Dr Yongfu Chen, Bing Xu, and Zihua Tang for help in field sampling and in the laboratory. This manuscript benefits considerably from thorough and thoughtful reviews by J Bloemendal and JW Zhang.
Funding
This study was supported by the China Geological Survey (grant number 1212011120045) and the National Natural Science Foundation of China (grants 41274074 and 41371219).
