Quantitative Reconstruction of Precipitation and Lake Areas During Early to Middle Holocene in Mu Us Desert, North China

Introduction

There were variations in monsoon precipitation in the geological data recorded in China during the Holocene (An et al., 2000). The data showed that the East Asian summer monsoon (EASM) strengthened, the monsoon rain belt moved northwestward, and precipitation increased in the early to middle Holocene (Yang et al., 2015a). Holocene precipitation variations in north China have been quantitatively reconstructed with multiple methods. In the Chinese Loess Plateau, Sun et al. (1995) established a susceptibility-rainfall conversion function by regression analysis of magnetic susceptibility of surface soils and contemporary precipitation data. The reconstructed precipitation was 100–150 mm higher than the present in the loess-desert transitional zone during the Holocene optimum. Precipitation reconstruction from the Gonghai Lake in the north loess plateau revealed gradually increasing rainfall since the last deglaciation, reaching a maximum of 30% higher than the present during the middle Holocene (Chen et al., 2015). In eastern Inner Mongolia Plateau, pollen records of Chagan Nuur showed that the East Asian monsoon precipitation began to increase in the early Holocene and reached a maximum in the middle Holocene, which was 30–50% higher than modern rainfall (Li et al., 2020). The lake hydrologic model and Holocene shoreline records were used to reconstruct the precipitation in Dali Lake, located near the northwestern limit of the EASM domain. The results showed that rainfall in the Dali Lake area increased by 11–5.5 ka, and the EASM expanded northwest (Goldsmith et al., 2017). The Bayanchagan Lake in southern Inner Mongolia also recorded a wet period in the early to middle Holocene, with annual precipitation approximately 32% higher than the present (Jiang et al., 2006). In the Alxa Plateau, the reconstructed average precipitation of the Jilantai playa during 8.5–3.5 ka was 246 mm, about 75% higher than the present (Wu et al., 2018). In the Badain Jaran Desert, Yang and Williams (2003) estimated precipitation of ∼200 mm during the early to middle Holocene, at least twice that of the modern period. The reconstructed precipitation in the Shiyang River drainage area during the middle Holocene was approximately 70% higher than the present (Liu and Li, 2017). The reconstruction results showed that the increase percentage in precipitation during the early to middle Holocene ranged from ∼30% to 100% compare with the present in north China, with significant differences near the northwestern limit of the EASM domain.

Lakes are treated as archives of earth history (Cohen, 2003). In structurally stable regions, lake level fluctuation is a sensitive response to catchment climatic change (Fritz, 2008). According to the quantitative reconstruction results, the climate in the early to middle Holocene was wetter than the present. Furthermore, the lacustrine landform and strata show the lakes’ highstand during this period. For instance, the lacustrine strata revealed by cores in the Ulan Buh Desert show that a unified paleolake developed in the northern part of the desert during the early to middle Holocene (Chen et al., 2014). The lake highstand in Badain Jaran Desert occurred at ∼10–6.6 cal ka BP (Wang et al., 2016). Evidence from ancient shorelines and lacustrine stratigraphy in Badain Jaran Desert indicate that the area of the lakes during this period was about five times that of modern lakes (Yang and Williams, 2003). The terrace of Baijian lake in Tengger Desert shows the lake level was ∼15 m higher than present at 8.5–5.1 ka (Zhang et al., 2004). Furthermore, lacustrine landform and strata recorded the existence of lake highstands during the early and middle Holocene in the south and east Inner Mongolia Plateau (Yang et al., 2011; Yang et al., 2015b; Zhang et al., 2016).

The Mu Us Desert (MUD) is one of the twelve sandy areas in China. It is located northwest edge area of the modern EASM (Figure 1A). Precipitation reconstruction is of great significance to the study of climate change in the EASM edge area and the evolution of the MUD. In addition, it is beneficial to understand the desert’s evolution and landform pattern by reconstructing paleolake area during the early to middle Holocene. Previous research has mainly focused on sediments in the southeast of the desert. Moreover, paleoclimatic records have been revealed from the aeolian sand-paleosol sequence, loess-paleosol sequence, and alluvial-lacustrine-peat deposition (Gao et al., 1993; Zhou et al., 1999; Jia et al., 2015; Liu et al., 2018a; Jia et al., 2018). Proxies such as magnetic susceptibility and pollen assemblage have been used to quantitatively reconstruct precipitation levels (Shi et al., 1988; Chen et al., 1993; Sun and Feng, 2013; Liu et al., 2017). The paleoclimatic research in the northwest of the desert has received relatively less attention. Moreover, current studies rarely involve the reconstruction of the paleolake area in the desert. Here, this study presents and discusses reconstructed precipitation and lake areas during the Holocene highstand based on the lakes in the closed drainage area of the west and north MUD.

FIGURE 1. Overview of the study area. (A) Map of northern China and location of study area. The modern isohyets are values averaged over 56 years (1961–2016); monsoon edge area is modified after Zhang et al. (1993); the red rectangle is the study area. (B) Topography of study area. The DEM data is from ASTER GDEM version2. (C) Diagrammatic section of groundwater flow system of SRBDY–DLTR Nuur drainage areas (modified after Hou et al., 2008); location of the section is shown by the yellow line in Figure 1B).

Regional Setting

The MUD, with an area of about 4 × 10⁴ km², is situated in the southeastern Ordos Plateau. Elevation ranges from 1600 m a.s.l. in the northwest to 1,200 m in the southeast. Bedrock hills from the arid denudation plateau in the northwest MUD extend to the southeast. Lakes and alluvial-lacustrine plains are distributed in valleys between bedrock hills. The aeolian landform presents a characteristically staggered distribution pattern of different forms of mobile, semi-fixed, and fixed dunes/land (Department of Geography of Peking University, 1983). The underlying bedrock in the MUD is Jurassic and Cretaceous sandstone, of which the weathered detritus provide one of the major provenances of the desert(Zhu et al., 1980). The mean annual precipitation decreased by about 450 mm to the northwest. Dominated by the East Asian monsoon, 60–70% of the annual rainfall occurs between July and September. The mean annual temperature ranges from 6.5 to 9°C.

There are 154 natural lakes in the MUD covering an area ≥0.1 km²(Wu et al., 2017). Most of the lakes are distributed in bedrock valleys or aeolian depressions. Jurassic and Cretaceous sandstone constitute the lake floor. These lakes are salt lakes with a mineral concentration ranging from 75 to 400 g/L, except for some freshwater or brackish water lakes such as Hongjian Nuur and Daotuhaizi in the southeast MUD(Zheng et al., 1992). The average depth of most lakes is within 2 m. Hongjian Nuur is the deepest lake with an average depth of 8.2 m at 1230 m a.s.l (Wang and Dou, 1998). Groundwater resources in MUD are abundant with a salinity less than 1 g/L. The groundwater table in the inter-dune depressions is within 1 m (Department of Geography of Peking University, 1983). The lakes are convergent centers of surface and groundwater runoff in local watershed (Hou et al., 2005).

The west and north MUD belong to the Ordos Plateau closed drainage area, showing neither surface nor groundwater hydraulic connection with the Yellow River and its tributaries (Hou et al., 2008). There are 108 natural lakes distributed in closed drainage area. The total area of these lakes is 301.1 km². These lakes can be divided into twenty-four drainage basins according to surface watershed (Figure 1B). The surface watershed is approximately consistent with the groundwater watershed (Hou et al., 2005). Each basin contains a single lake or catenulated lakes connected by seasonal river channels. Groundwater renewal ability decreases with depth in the closed drainage area. The annual renewal rate of medium and deep groundwater flow is less than 0.5% and 0.02%, respectively. The vertical hydraulic connection between shallow groundwater and medium to deep groundwater is weak (Wan et al., 2010). Therefore, the lakes in each drainage basin are independent discharge centers of surface runoff and shallow groundwater flow (Figure 1C).

Materials and Methods

Lacustrine Landform and Sample

Lacustrine plains, lake terraces, and lake shorelines are the landform relics of lake evolution. The lake terrace relics and corresponding lacustrine deposits, which represent the stable stage of lakes during their expansion period, are direct evidence of the evolutionary history of lakes and for reconstructing high lake levels. The lakes on the Ordos Plateau have been retreating since the late Quaternary. Hongjian Nuur and Bojianghaizi developed three lake terraces (Figure 2). T2 and T3 developed before the late Pleistocene, and T1 developed in the early and middle Holocene (Yuan et al., 1987; Shi, 1991). Holocene lake terrace and lacustrine strata in the MUD are generally no more than 10 m higher than the modern lake surface. The lacustrine sediments are gray-white or gray-green silty sand and fine sand with horizontal bedding. The sedimentary texture is loose with poor cementation. The lake terraces are stepped terrain with a flat top and a steep front edge. Lacustrine strata remain partially on the terrace. Based on remote sensing image interpretation, lake basins of Hongjian Nuur (HJN), Dalengturu Nuur (DLTR), Baozhai Nuur (BZN), Sharibuduyin Nuur (SRBDY), and Bojianghaizi (BJHZ) were selected for field investigation and sampling (Figure 1B, No.1–5).

FIGURE 2. Geomorphology and lacustrine strata of lakes in MUD. (A) Image of lake shorelines and T1in south HJN. (B) Sketch map of lake terraces of HJN (modified after Shi, 1991). (C) Image of lake terrace in the south DLTR. (D) Lacustrine strata in the terrace of the south DLTR. (E) Holocene lacustrine strata in the southeast beach of BZN. (F) Lake terrace and lacustrine strata in the east SRBDY. (G) Sketch map of lake terraces in BJHZ (modified after Long et al., 2007; Shi, 1991).The red star, circle, and square show OSL, TL, and ¹⁴C dates, respectively. The ¹⁴C dates were converted to calendar ages using Calib 7.02 based on the INTCAL 13 calibration (Reimer et al., 2013).

HJN is the largest lake in MUD. Shorelines and lake terraces are clearly visible in south HJN. The shorelines formed in recent decades are roughly parallel to the modern shoreline (Figures 2A,B). The top of the lacustrine strata in T1 is about 9 m above the modern lake. The sediment is gray-white very fine sand. The terrace of DLTR is in the south of the lake. The terrace is about 5.5 m higher than the modern lake, slightly inclined toward the lake (Figure 2C). The lacustrine strata is gray-green fine sand with clay, about 17 cm in thickness (Figure 2D). The Holocene terrace of BZN is located on the southeast beach of the lake about 9 m above the modern lake. The lacustrine layer is gray-white silt sand with well-developed horizontal bedding (Figure 2E). The lake terrace of SRBDY, located in the east of the lake, is about 8 m higher than the modern lake. The lacustrine strata at the top is about 40 cm, with gray-white fine sand and visible and well-developed horizontal bedding (Figure 2F). The T1 in the north of Bojianghaizi is about 7 m higher than the modern lake with gray-green and gray-black lacustrine clay (Figure 2G).

Two quartz Optical Stimulated Luminescence (OSL) dating samples were collected in the DLTR section. Pure quartz grains 38–63 μm size were extracted from the middle part of the sample tube using the procedure employed by the OSL Laboratory of the State Key Laboratory of Earth Surface Processes and Resource Ecology, Beijing Normal University. The SAR-SGC method was used to determine the equivalent dose (De) (Lai and Ou, 2013).The measurement equipment was a Risø DA-20 TL/OSL reader equipped with blue diodes (k = 470 ± 20 nm) and a 90Sr/90Y radioactive beta source. The concentrations of Uranium (U), Thorium (Th), and Potassium (K) were measured by inductively coupled plasma mass spectrometry (ICP-MS). The water content was measured by weighing the samples before and after drying. Finally, we obtained the boundary ages of the lacustrine deposit in the DLTR section.

Lake Hydrologic Model

The lake area is governed by lake inputs and outputs. The lake level remains stable when the water inputs equal outputs. Lake inputs are lake surface precipitation and drainage runoffs, whereas the outputs are water surface evaporation and groundwater discharge. Based on the reconstructed lake area, the lake hydrologic model can be used to calculate the precipitation of the corresponding period. The model used in this study followed Goldsmith et al. (2017):

where $\mathrm{P}$ is mean annual precipitation(m); $A_C$ is the catchment area (m²) where the area of the lake is subtracted from the area of the catchment; $F_R$ represents the fraction of the watershed precipitation that reaches the lake (including surface and groundwater runoff); $A_L$ is the area of the lake (m²); $E$ represents the mean annual evaporation of the lake surface (m); and $G$ represents groundwater discharge from the lake per year (m³).

The catchment area ($A_C$) was calculated from topographic maps with contour lines of 5 and 10 m. Considering the inter-annual fluctuations of the lake level (Xu et al., 2019), the modern area of the lake ($A_L$) was averaged from data of two periods: a topographic map drawn in the 1970s and a geomorphic map of the MUD (Wu et al., 2017). The precipitation was interpolated from the mean annual precipitation of meteorological stations in the MUD for 56 years (1961–2016). The highstand area of five lakes was calculated from topographic maps according to the elevation of the lacustrine strata. The watershed area of the lake highstand was the same as the modern catchment area because the elevation of each watershed was much higher than the modern lake and Holocene highstand. All the catchments are closed basins with no hydraulic connection between them. Thus, groundwater discharge from the lake ($G$) can be neglected for the hydrologic model. The meteorological data of each basin was obtained by interpolation.

The evaporation of the lake surface was estimated using a modified Penman equation (Penman, 1948; Sene et al., 1991):

where $E_d$ is evaporation per day (m/day); ∆ is the slope of the saturation vapor pressure curve at air temperature t (hpa/°C); H is net radiation, J/(m²∙day); $L$ represents latent heat of vaporization (J/kg); γ is psychrometric constant (hpa/°C); and $E_a$ is atmospheric conductance (m/day).

The climatological data (radiation, temperature, air pressure, vapor pressure, wind velocity, and sunshine duration) were obtained from the Chinese Meteorological Agency for the period 1961–2016 for meteorological stations located in the MUD. The mean annual evaporation of each meteorological station was calculated by summing evaporation per day ($E_d$). The evaporation of each lake surface ($E$) was mathematically interpolated from these stations. The runoff fraction ($F_R$) can be calculated using the closed-basin hydrological model. In addition, it can be estimated by the Budyko relationship (Budyko, 1974; Koster et al., 2006):

where ${\mathrm{E}}_c$ is catchment evaporation calculated as D = $R_{net}$/P $L$, where D is the dryness index, $R_{net}$ is mean annual net radiation, $P$ is annual precipitation, $L$ is latent heat of vaporization, and $R_{net}$/ $L$ is potential evaporation. Potential evaporation can be measured by the radiative balance or by pan evaporation (Koster et al., 2006; Goldsmith et al., 2017). It also can be approximately estimated using the Penman equation (Yang et al., 2012; Li et al., 2020).In our calculation, the Penman equation was used.

In the modern lake hydrologic model, the mean annual evaporation of meteorological stations in the MUD was calculated using the Penman equation. Then, the evaporation ($E$) of the five lake basins was calculated by mathematical interpolation. The modern runoff fraction ($F_R$) of each basin was calculated using the lake hydrologic model and Budyko relationship, respectively. The results show that the two runoff fractions ($F_{R1}$ and $F_{R2}$) are similar, indicating that both the modern lake hydrological model and Budyko relationship produced suitable estimates.

In this paper, first, the model was applied to HJN, DLTR, BZN, SRBDY, and BJHZ to reconstruct precipitation during the lake highstand. Second, the paleolake areas of other drainages (No. 6–24 in Figure 1B) were estimated using this model under the reconstructed precipitation pattern. In the paleolake hydrologic model, the evaporation was also calculated using the Penman equation. Considering the influence of wind speed and temperature on evaporation, we estimated that the evaporation during the lake’s highstand was at most 15% higher than the modern value. For the paleolakes, precipitation (or paleolake area) and runoff fraction could not be constrained independently. Therefore, we used both the closed-basin lake hydrologic model and the Budyko relationship to assess these parameters.

Extraction of Published Precipitation Data

A series of precipitation data was extracted from published articles to reconstruct the paleoprecipitation pattern in MUD combining with reconstructed precipitation by the model in section 3.2. The published precipitation data were extracted according to the following methods:1) the paleoprecipitation was extracted immediately from articles with exact data; 2) the paleoprecipitation was calculated using modern precipitation (1961–2016) multiplied by change percentage given in the articles; 3) with regard to data in precipitation series, data from the chosen period were digitalized and averaged. With the method above, forty precipitation data was extracted in the MUD and its surrounding region in north China.

Results

Lake Highstand Reconstructed by Lacustrine Landform and Strata

Based on field investigations, the highstand of five lakes was reconstructed using a topographic map and a digital elevation model (Figure 3). The area of the paleolakes is 1.65–3.56 times that of modern lakes, with the smallest ratio in HJN. The total area of paleolakes in the five catchments is 158.2 km², which is 2.05 times that of modern lakes. There was only one lake during the highstand period in the DLTR and SRBDY drainage. The drainage basins in HJN, BJHZ, and BZN contain more than one lake. In each basin, the elevation of the watershed between lakes is the upper limit of the upstream paleolake level. We also inferred that paleolakes existed in some depressions during the highstand period. For example, the Shejiahaizi (haizi means lake in the dialect of study area), located in depression northeast of HJN, was a lake that regressed decades ago according to the topographic map and remote sensing image. During the Holocene highstand period, Shejiahaizi was probably an independent lake that was ∼3 m higher than the HJN based on lacustrine landform. In the BZN drainage, lake Baozhai Nuur and the adjacent Baga Nuur were connected to a unified lake due to rising lake levels. In the west of the drainage, there are two depressions where paleolakes may develop. In the BJHZ drainage, two adjacent lakes were connected during the high lake level period.

FIGURE 3. Reconstructed highstand of five lake basins in the MUD (the DEM data is from ASTER GDEM version 2).

The quartz optically stimulated luminescence (OSL) ages of the lacustrine strata in DLTR are listed in Table 1. The quartz OSL ages confirmed the lake highstand occurred during the early to middle Holocene.

TABLE 1. Quartz OSL ages of lacustrine strata in DLTR.

Quantitative Precipitation Reconstruction Based on the Lake Hydrologic Model

The modern precipitation ($P$), modern and highstand evaporation ($E$), catchment area ($A_C$), and runoff fraction ($F_R$) required in the hydrologic model are listed in Table 2. The interpolated modern mean annual evaporation ranged from 995 to 1,065 mm with a rough trend of increasing from east to west. The evaporation of the paleolakes was calculated using modern values multiplied by 1.15. The runoff fraction ($F_R$) of the paleolake drainages during lake highstand was estimated by Budyko relationship. They ranged from 5% to 8%, higher than at present in each catchment. The runoff fraction ($F_R$) in both periods shows a decreasing trend from southeast to northwest.

TABLE 2. The hydrological parameters of each lake basins in the MUD during the modern and highstand period.

The precipitation of five lakes during the highstand was estimated using the lake hydrologic model (Table 2). The reconstructed precipitation of HJN in the eastern MUD was 500 mm. In the middle MUD, the precipitation of BZN and BJHZ was 471 and 472 mm, respectively. The precipitation of DLTR and SRBDY in the west was 443 and 444 mm, respectively. The reconstructed precipitation was 120–165 mm higher than the modern value. The increase in precipitation was 32–60%, and the increase in the west was significantly higher than that in the east.

Reconstruction of Precipitation Patterns

Precipitation patterns in the MUD during the early to middle Holocene were reconstructed by synthesizing precipitation data in section 4.2 and published data. A total of forty-five precipitation data from the early to middle Holocene in the MUD and surrounding areas were used (Supplementary Table S1). The isohyet of the MUD was reconstructed by interpolation (Figure 4). The results show that the precipitation during the early to middle Holocene in the MUD decreased from about 550 mm in the southeast to about 350 mm in the northwest, which is consistent with present levels. The 400 mm isohyet moved 130–170 km to the northwest, roughly consistent with the modern 250 mm isohyet. The moving distance in the southwest is slightly smaller than that in the northeast. The amplification of precipitation gradually increased from southeast to northwest. Due to the difference in amplification, the precipitation gradient in much of the MUD is lower than at present. The large gradient is distributed near the 400 mm isohyet in the west MUD.

FIGURE 4. Reconstructed precipitation and lake pattern in the MUD (the legends are same as Figure 3).

Reconstructed Paleolakes

The lake highstand of five drainage areas was reconstructed based on lacustrine landform and strata in section 4.1. Based on the reconstructed precipitation pattern in the MUD, the paleolake area of each drainage (No. 6–24) during the early to middle Holocene was estimated. First, the runoff fraction (F_R2) of each drainage basin was estimated using the Budyko relationship. Second, the lake area of each basin was calculated using the lake hydrologic model (Table 2). The total area of reconstructed paleolakes during the early to middle Holocene was 896.1 km², which is 2.96 times that of modern lakes. The lake area of each drainage is two to five times that of modern lakes with the largest increase in the southwest MUD. The area of the paleolake in the Wuhudong Nuur, Gouchi, and Beidachi basins (No. 22–24 in Table 2) are 4.57, 4.84, and 4.97 times that of modern lakes, respectively.

Based on the reconstructed lake area in each drainage (No. 6–24), we also attempted to reconstruct the paleolake pattern. For the catchment comprising one single lake, the paleolake expanded by topographic contour around the modern lake based on the calculated area. In the basin comprising a series of catenulated lakes, the elevation of the watershed between lakes was regarded as the altitude of the paleolake upstream. The upstream paleolake area was calculated using closed contour and the corresponding altitude. The downstream lake area was calculated as the total area minus the upstream lake area. Due to the rise in lake levels, adjacent lakes may have become unified lakes. Therefore, the modern 108 lakes decreased to eighty-four during the early to middle Holocene. In addition, there are small, closed depressions in certain catchments. These depressions were inferred as paleolakes during the Holocene highstand. A total of forty-three inferred lakes with an area more than 0.1 km² were identified. According to the method above, there were 127 paleolakes in the closed flow area of MUD during the early to middle Holocene. The paleolake pattern is shown in Figure 4.

Discussion

Estimation of Evaporation Used in the Lake Hydrologic Model

As one of the required parameters in the lake hydrologic model, open water evaporation is typically estimated by observing different pans. However, the observation data of a 20 m² evaporation pan, which is approximately equal to open water evaporation, is rare in the MUD. The observation data of D20 and 601B pans are significantly higher than open water evaporation in arid and semi-arid areas because of the pan wall. In certain regions, open water evaporation is estimated using observation data of pans multiplied by a conversion coefficient (Linacre, 1994). However, the conversion coefficient in the MUD was difficult to constrain. Therefore, a modified Penman equation was used to estimate the evaporation of the lake surface (Penman, 1948; Sene et al., 1991). The Penman equation has been used to calculate evaporation in arid and semi-arid China (Li et al., 2001; Yang et al., 2010). The difference between calculated values and measured data in 20 m² evaporation pans was no more than 10% according to experiments conducted in the arid regions of Xinjiang, in western China (Zhang and Zhou, 1992).

The observed evaporation of two pans and calculated evaporation using the Penman equation are listed in Table 3. The ratios $E_{D20}$/ $E$ and $E_{601B}$/ $E$ are 0.4–0.52 and 0.81–0.97, respectively. Both ratios decreased from the semi-arid southeast to the arid northeast. The average conversion coefficient between D20 pan and open water evaporation in arid northwest China and Inner Mongolia Plateau was estimated at 0.44 and 0.52, respectively (Xie and Wang, 2007). The $E_{D20}$/ $E$ ratio in the MUD was roughly between these two coefficients. The conversion coefficient between 601B pan and 20 m² pan in North China ranged from 0.83 to 0.93, which is close to the $E_{601B}$/ $E$ ratio in the MUD. The decreasing trend from semi-humid Henan to semi-arid Inner Mongolia was also consistent with the MUD (Qian et al., 1998). This indicates the evaporation calculated using the Penman equation was approximately equal to the open water evaporation.

TABLE 3. Evaporation of meteorological station in MUD.

In the paleolake hydrologic model, the evaporation of lakes was also estimated using the Penman equation. The temperature during the early to middle Holocene was about 1°C higher than modern times (Fang et al., 2011). Based on the sensitivity of evaporation to temperature change, a 1°C rise would produce evaporation rate changes of about 5–6% in the MUD. We also examined the sensitivity of evaporation to wind speed change. A 20% change in wind speed can lead to about a 10% change in evaporation, which is nearly consistent with the examination in Dali Lake in eastern Inner Mongolia (Goldsmith et al., 2017). Thus, considering the possible change in temperature and wind speed, the evaporation was assumed at 15% higher than the modern value at most.

Reconstructed Precipitation in the MUD

The lake hydrologic model is one of the common methods for reconstructing precipitation. The precipitation reconstructed in HJN was 500 mm. The result is approximate to the value 510 ± 25 mm estimated from pollen data from the Hejialiang section about 15 km east of the lake (Chen et al., 1993). The reconstructed precipitation in BZN was 471 mm. The rainfall calculated based on pollen assemblage in Taobao and Wushen profiles was 470 ± 25 and 461.8 mm, respectively (Shi, 1991; Chen et al., 1993). These two profiles are about 7 and 15 km to the southeast of BZN, respectively. This indicated that the lake hydrologic model and pollen assemblage produced suitable estimates. The mean annual precipitation of the HJN lake was ∼380 mm. The precipitation during the early to middle Holocene was ∼32% higher than present. The result is consistent with reconstructed results near the 400 mm isohyet in the MUD and southeast Inner Mongolia Plateau (Jiang et al., 2006; Xu et al., 2010). However, the result was much smaller than the estimated precipitation of Qinghai Lake (Shi et al., 1992). This may reflect that the migration distance of the monsoon rain belt varied in different regions of the monsoon edge area.

In the EASM area of China, the relative increase in precipitation during the early to middle Holocene varied from ∼10 to 15% in the humid southeastern coastal area to more than 70% in the arid northwest area (Yang and Williams, 2003; Stebich et al., 2015; Liu and Li, 2017; Li et al., 2018; Wu et al., 2018). The reconstructed relative increase in the MUD rose from ∼32% in the southeast to ∼60% in the northwest. The rising trend is consistent with the whole EASM area in China. However, the extent of relative increase ranges significantly from southeast to the northwest, indicating that the monsoon edge area was extremely sensitive to climatic change. During the early to middle Holocene, the precipitation gradient decreased significantly in the MUD due to the difference in precipitation increase. This probably reflects the monsoon rain belt and the monsoon edge area migration to the northwest. In consequence, the MUD was stably dominated by the EASM. The spatio-temporal pattern of climate showed a relatively stable, warm, and humid climate for millennia. It can also be confirmed from the peat and paleosol development process in the MUD. The development of peat and lacustrine in the MUD since the last deglaciation indicates the regional climate fluctuated considerably until 8.5 ka. Thereafter there were humid intervals until 6ka, during which lacustrine sediments were deposited (Liu et al., 2018b). The paleosol was generally regarded as the indicator of warm and humid climate. In the aeolian sand/paleosol sequence of the MUD, the duration of paleosol is much longer than the aeolian sand during the early to middle Holocene (Gao et al., 1993; Ding et al., 2021).

Paleolake Pattern in the Closed Drainage Area

The lakes in northern China generally reached their highstand due to the increase in precipitation in the early and middle Holocene. The lake landforms and lacustrine strata of HJN and BJHZ are direct evidence of the development of high lakes in the study area. The decrease in the grain size of lacustrine deposits in the lake cores of the Baahar Nuur and Qigai Nuur in north MUD indicates an increase in lake level during the early and middle Holocene (Guo et al., 2007; Sun and Feng, 2013). The clay deposits in Chagan Nuur in the western MUD indicate a stable deep-water sedimentary environment of ∼9–5 ka (Yang, 1997). Clay sediments in Hetongchahagan Nuur in the western MUD indicate a stable deep-water sedimentary environment at∼9—6 ka (Bernasconi et al., 1997). Grain size analysis of the DSGL profile in the southern MUD showed that the Salawusu Lake began to form in the early Holocene and reached its maximum in the middle Holocene (Liu et al., 2018b). These lake records above indicate that the lake highstands during the early to middle Holocene in the MUD were almost simultaneous and stable.

The total area of lakes in the closed drainage area is 2.96 times that of modern lakes. The relative height between the highstand and the modern lake surface is about 5–9 m. The increase of height and area was smaller than the lakes in Badain Jaran Desert and Tengger Desert (Yang and Williams, 2003; Zhang et al., 2004). Most of the lakes in the study area locate in bedrock valleys or aeolian bedrock depressions. The topography of the lake basins is mainly controlled by bedrock gullies. The modern aeolian activity has little impact on topography. The stable terrain is conducive to the reconstruction of paleolake pattern. According to the reconstruction results, the number of paleolakes increased along with the rise of the lake level and the expansion of the lake area. The adjacent lakes in certain basins connected and formed unified lakes due to the rise in water level. In addition, closed depressions were extracted from topographic maps and remote sensing images. Geomorphic relics, such as lake shorelines, can be distinguished around the depressions. Salt crystals on the surface of depressions indicated paleolakes may have existed (Figure 5). We inferred these closed depressions formed lakes due to the rise of precipitation during the early to middle Holocene. In each catchment, the altitude of the lake highstands was lower than the watershed between drainages. The pattern of closed drainage area is consistent with the present. There is no evidence from the lake landform, lacustrine strata, and climatic conditions supporting the formation of a unified megalake in the closed drainage area during the early to middle Holocene.

FIGURE 5. Inferred paleolake in the MUD during the early to middle Holocene (The white curves represent contour lines extracted from topographic maps, the number means elevation/m. The blue line with arrow indicates the possible flow direction. (A–D) correspond to locations in Figure 4).

Conclusion

The precipitation and lake areas of the MUD during the early and middle Holocene were quantitatively reconstructed based on the lacustrine landform, and by using the lake hydrologic model. A total of 127 paleolakes existed in the closed flow area during the early to middle Holocene. The paleolake area was 896.1 km², which is 2.96 times that of modern lakes. The relative height between the highstand and the modern lake surface is about 5–9 m. The precipitation during early to middle Holocene decreased from about 550 mm in the southeast to about 350 mm in the northwest. The 400 mm isohyet moved 130–170 km to the northwest, roughly consistent with the modern 250 mm isohyet. The relative increase in precipitation was ∼32–60%, and the increase in the west was significantly higher than the east. The precipitation gradient in much of the MUD was lower than the present. The results above indicate that the monsoon edge area and monsoon rain belt have migrated to the northwest. Furthermore, the MUD was stably dominated by the EASM. The spatio-temporal pattern of the climate during the early to middle Holocene was relatively warm and humid with a decreased precipitation gradient for millennia.

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 author.

Author Contributions

DL, YWu, and LT contributed to conception and design of the study. DL, LT, YWe, and TF conducted investigation and methodology. DL, YWe, and TF performed the statistical analysis. YWu provided the funding. All authors contributed to manuscript revision, read, and approved the submitted version.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2022.850633/full#supplementary-material

References