Aspects of human physical and behavioural evolution during the last 1 million years

This paper reviews some of the main advances in our understanding of human evolution over the last 1 million years, presenting a holistic overview of a field defined by interdisciplinary approaches to studying the origins of our species. We begin by briefly summarizing the climatic context across the Old World for the last 1 million years before directly addressing the fossil and archaeological records. The main themes in this work explore (i) recent discoveries in the fossil record over the last 15 years, such as Homo naledi and Homo floresiensis; (ii) the implications of palaeogenetics for understanding the evolutionary history of, and relationships between, Neanderthals, Denisovans and Homo sapiens; (iii) the interplay between physiology and metabolic demand, landscape use, and behavioural adaptations in the evolution of morphological and behavioural innovation; and (iv) recent advances in archaeological understanding for the behavioural record, in particular that of the Neanderthals. This paper seeks to provide a broad‐scale, holistic perspective of our current understanding of human evolution for the last 1 Ma, providing a reference point for researchers that can be built upon as new discoveries continue to develop the landscapes of human evolution. © 2019 The Authors. Journal of Quaternary Science Published by Wiley Periodicals, Inc.     

Global chronostratigraphical correlation table for the last 2.7 million years

The table provides a correlation of chronostratigraphical subdivisions of late Cenozoic geological time, spanning the last 2. 7 million years. The formal division of the Quaternary is the responsibility of the International Commission on Stratigraphy＇s （ICS） Subcommission on Quaternary Stratigraphy （SQS）, in partnership with the International Union for Quaternary Research＇s （INQUA） Commission on Stratigraphy and Chronology （SACCOM）. This is the third published version of the chart. Earlier versions are Gibbard et al., 2004, 2005. See http：//www.quaternary.stratigraphy.org.uk/correlation/ for history.     

Modelling the long‐term fluvial erosion of the River Somme during the last million years

A process‐based model that simulates fluvial erosion in the River Somme Valley over the last million years is presented here. The model takes into account lithology and climatic influences and allows the simulating of undercapacity and overcapacity sediment transport behaviour. The model has been calibrated to a family of terraces within the River Somme Valley. When matched to this field data, simulation trials suggest that bedrock incision occurred principally from 120 to 60–40 kyr during the last climatic cycle and before the last glaciation. The impact of a progressive tectonic uplift (c. 60 m over c. 1 million years) on the River Somme has also been studied here. Extended over a longer period of time, the simulations suggest that 1 million years ago the profile of the River Somme had a lower slope gradient than today, with little relief throughout the Paris Basin.     

The last million years of deformation in part of the New Zealand plate-boundary zone

Tectonic deformation of Pleistocene sediments is occurring within the plate-boundary zone in the North Island of New Zealand. Folds, associated with an active fault, are forming mainly by flexural-slip, and contain comformable and unconformable relations which show that the folding started ca 1.0 Ma. Simple geometric models of fold growth indicate that the fold structure has shortened at an average rate of 0.76 ± 0.1 mm a⁻¹ across 4.5 km, though with an increase in the last 100 Ka of about 1.8 times the average rate, which is entirely compatible with the rate of shortening observed in the last 50 years from retriangulation studies. Areas for 15 km to the southeast of the fold structures show significantly lower horizontal shortening rates. The strain rates across the plate-boundary zone vary both with position and time, with a negative correlation between the uplift pattern and the inferred geological shortening rates. This can be explained by the presence of large semi-rigid blocks within the plate-boundary zone which are back-tilting and overriding weaker more rapidly deforming areas. The available evidence suggests that a large proportion of the deformation is occurring seismically.     

A numerical model of the evolution of ocean sulfate and sedimentary sulfur during the last 800 million years

A reservoir model describing the time evolution of the sedimentary cycle of sulfur over the past 800 my has been developed. As a first approximation, the ocean sulfate concentration is assumed to be time-independent. With this assumption, the model is integrated backward in time and a new initialization procedure is derived in order to calculate the present state of the system which must be compatible with both observational data and model equations. The effects of a variation of the present state of the cycle on its past evolution are investigated. It is found that, when the present gypsum reservoir content is too low or when the weathering rate constants are too high, no acceptable solution can be obtained for the evolution of the cycle, since one reservoir is forced to depletion. The sensitivity of the model to the mean isotopic composition of the sedimentary system and to the fractionation factor during pyrite formation is also studied.Moreover, a model with time-dependent ocean sulfate concentration was developed. The existence of an acceptable solution appears to be linked to the steady state hypothesis for ocean sulfate, since a model with no acceptable steady state solution may be integrated until t = −800 my without any problem of reservoir depletion when the time-dependent equations are used.A tentative evolution of the ocean sulfate concentration is calculated. It is shown that this concentration is negatively correlated to the δ³⁴S of seawater sulfate. The carbon cycle is modelled in order to compare the calculated δ¹³C of carbonate deposits to the observational data.     

Accretion of extraterrestrial matter during the last 80 million years and its effect on the marine osmium isotope record

The flux of particulate extraterrestrial (ET) matter to the deep-sea has been calculated using a four-component mixing model based on osmium concentrations and isotope ratios in slowly accumulating pelagic sediments from the Pacific Ocean. Nineteen published bulk-leach osmium isotope data pairs that cover the last 80 million years have been used for the calculation. The calculated annual particulate flux, averaged over several 100 ka by slow accumulation and bioturbation, ranges from ∼18,000 to ∼67,000 tons, with a mean value of ∼37,000 ± 13,000 tons. The data indicate no significant variability in the flux with time, except at the K-T boundary and are thus compatible with results based on Ir accumulation in deep-sea sediments. The inferred constancy of the ET matter flux is in contrast to recent results based on the ET ³He flux that integrates only a small size fraction of the entire ET matter flux. ³He data indicate variability by a factor of 6 over the past 70 Ma and a fivefold increase in the ET flux between 3 Ma and 1 Ma that is not seen in the bulk ET matter flux based on osmium isotope data.The apparent constancy of the cosmic matter flux over the past 80 Ma stands in marked contrast to the dynamic nature of the marine osmium isotope record and indicates that dissolution of cosmic matter does not drive changes in this record, except at the K-T boundary. A comparison between osmium isotope and Ir data from two pelagic sediment cores from the Pacific (DSDP Site 596 and LL44-GPC3) indicates that the seawater-soluble fraction of ET Os cannot exceed ∼36 kg/a and most likely is significantly smaller. The maximum value of 36 kg/a can account for approximately one half of the seawater-soluble ET matter flux necessary to balance the radiogenic continental runoff of Os without any additional unradiogenic source. Simple mass balance calculations indicate that an additional unradiogenic source of Os to the oceans, most likely alteration of oceanic crust, is required to balance the present-day seawater osmium-isotopic composition. This source is probably more important in balancing radiogenic continental runoff than is dissolution of cosmic matter in seawater.     

The last 1.35 million years at Tenaghi Philippon: revised chronostratigraphy and long-term vegetation trends

In addition to being of interest to ancient Greek and Roman historians, the site of Philippi, NE Greece, has long been noted in Quaternary circles for providing the longest continuous European pollen record, spanning the last one million years. Here the original age model is re-evaluated and a new marine-terrestrial correlation is proposed. An astronomical calibration procedure, based on a correspondence between changes in certain vegetation elements and March and June perihelion configurations, suggests that the base of the sequence extends back to 1.35 million years ago. The revised chronological framework for the Tenaghi Philippon sequence provides an opportunity to examine the long-term behaviour of individual taxa and vegetation trends within the context of global climate changes. Comparisons reveal a close correspondence between the terrestrial and marine records, in terms of orbital and suborbital variability. However, joint time-frequency analysis of the arboreal pollen record shows that the obliquity and eccentricity/precession signals persist into the ‘100-kyr’ and ‘41-kyr’ worlds, respectively, suggesting the operation of additional climate mechanisms that are independent of high-latitude glacial–interglacial effects. Unlike ice core and marine sequences, no change in the magnitude of interglacial tree population expansions is observed after the Mid-Brunhes Event. Instead, the Tenaghi Philippon record suggests a major shift in the vegetational composition of interglacials after MIS 16, with the establishment of forests of reduced diversity and a ‘modern’ appearance.     

Storminess variation during the last 6500 years as reconstructed from an ombrotrophic peat bog in Halland,southwest Sweden

Cores taken from an ombrotrophic peat bog in the coastal zone of Halland, southwest Sweden, were examined for wind transported mineral grains, pollen and humidity indicators. The core covers the period from 6500 cal. yr BP to present. Ombrotrophic conditions existed from ca. 4200 cal. yr BP onwards. Bog surface wetness fluctuated strongly until ca. 3700 cal. yr BP, with an apparent dominance of dry summer conditions from 4800–4500 cal. yr BP. Local wet shifts occurred around 4300, 2800, 2400 and 1500 cal. yr BP, whereas the most recent 600 years of the record show increasingly dry conditions. Mineral grain content, interpreted as aeolian sand influx (ASI), was used as a proxy for (winter) storm frequency and intensity until ca. 1500 cal. yr BP, after which increasing human impact, as reconstructed by pollen analysis, became a second important potential cause for increased sand drift. Strongly increased storminess occurred at 4800, 4200, 2800–2200, 1500, 1100 and 400–50 cal. yr BP, indicating a dominance of cold and stormy winters during these periods. Many of these storm periods apparently coincide with storm events in other sites in southwestern Scandinavia, suggesting that our ASI record reflects a regional scale climatic signal. Furthermore these stormy periods correlate to well‐known cold phases in the North Atlantic region, suggesting a link to large‐scale fluctuations in atmospheric circulation patterns. Copyright © 2006 John Wiley & Sons, Ltd.     

Unstable forearc stress in the eastern Nankai subduction zone for the last 2 million years

Transpressional tectonics characterizes the SW Japan arc. However, we will show in this article that offshore seismic profiles and onshore mesoscale faults indicate that the eastern part of the forearc was subject to transtensional tectonics since ca. 2.0 Ma. Offshore normal faults imaged on the profiles run parallel to the Nankai Trough, and started activity at 1.0 Ma, but transtensional tectonics commenced the onshore area earlier. In order to understand the stress history in the forearc region, we collected fault-slip data from onshore mesoscale faults in Plio-Pleistocene sedimentary rocks in the Kakegawa area at the northeastern extension of the offshore normal faults. Most of the mesoscale faults are oblique-normal, indicating that the area was subject to transtensional tectonics. The faults suggest that the compressional tectonic regime was followed by the transtensional one at 2.0 Ma, in agreement with regional tectonostratigraphic data, which indicate that folding ceased at that time. Present compressional stress followed the transtensional tectonic regime sometime in the late Pleistocene. Transtensional or extensional tectonic zone shifted from the Kakegawa area to the offshore region.These observations indicate that the state of stress just behind the accretionary prism of the eastern Nankai subduction zone has been unstable in the last 2 million years, suggesting that the forearc wedge has been at critical state in that gravitational force and basal shear traction on the wedge have been balanced, but the forearc tectonics has been susceptible to small perturbations. Possible factors compatible with the observed stress history include the change of subduction direction of the plate at 1.0 Ma, and the rapid uplift of Central Japan thereafter.     

Paleoceanography of the Santa Barbara Basin during the last 8000 years

The varved sediments of the Santa Barbara Basin off southern California, offer a unique opportunity to study the changes in oceanographic conditions of this nearshore area during the last 8000 yr. Quantitative analysis of Radiolaria found in recent surface sediment samples from the eastern North Pacific allows the identification of four “assemblages” which can be related to the physical oceanography of the California Current. Two assemblages are associated with the southerly flowing California Current, one with the main stream of the current (California Current Assemblage) and the other with the offshore flow along northern California (Central Assemblage). The two other assemblages are associated with the subtropical region of the eastern North Pacific (Subtropical Assemblage) and one found mostly off the coast of Baja California (Baja Assemblage). Analysis of the Radiolaria found in the varved sediments of a core from the Santa Barbara Basin give an 8000-yr continuous record of these four assemblages. The California and Baja Assemblages show only minor fluctuations in their importance in the sediments of the Santa Barbara Basin. The California Assemblage, however, shows a steady increase during this time period. Prior to 5400 yr B.P. the Radiolaria were predominately subtropical in character, whereas after 5400 yr B.P. the Central Assemblage becomes more important. Since 5400 yr B.P. most of the changes in the radiolarian fauna consist of fluctuations in the importance of these two assemblages. Past sea-surface temperatures for the month of February were calculated using the transfer function technique of J. Imbrie and N. G. Kipp (1971, In “The Late Cenozore Glacial Ages” L. K. Turekian, Ed.), (Chap. 5, Yale Univ. Press, New Haven, Conn.). The time series of paleotemperature estimates show major changes in the average February temperature of Santa Barbara Basin waters. The range of estimated temperatures (12°C) exceeds that of the historical observations of February temperatures in the Santa Barbara Basin but does not exceed the observed range for the California Current region. The intervals from 800 to 1800 yr B.P. 3600 to 3800 yr B.P. and 5400 to the end of the record appear to have been generally warmer than today. Comparison of the Holocene record of alpine glacial advances with the radiolarian assemblage and paleotemperature time series shows that the initiations of advances was coincident with a decrease in sea-surface temperatures and an increase in the importance of the Central Assemblage in the Santa Barbara Basin. The terminations of these advances were not marked by any consistent characteristic in the Santa Barbara Basin time series.     

The relation of the southern trade winds to upwelling processes during the last 75,000 years

In the equatorial Pacific, between the Galapagos Islands and the coast of South America, two kinds of upwelling of oceanic waters occur. One is related to coastal upwelling and the other to surfacing of the Equatorial Undercurrent. Both of those processes are associated with the development of the southeast trade winds blowing in this area. Coastal upwelling is increased when the trade winds are intensified, and the surfacing of the Equatorial Undercurrent occurs when the trades weaken. The development of coastal upwelling and the surfacing of the Equatorial Undercurrent are inferred from the radiolarian assemblages in the sediments. The abundance of quartz, opal, and radiolarian assemblages in the deep-sea sediments of this area, as well as the distance from the sample locations to land and to the quartz source, is correlated with the intensity of the trade winds (in February and August) through multiple regression analysis. The chronostratigraphy of core V1929 (3°35′S, 83°56′W), used in this study, is inferred on basis of its δ¹⁸0 record. During the last 75,000 years, the fluctuations in intensity of the trade winds have been concurrent with or preceded the fluctuations in the amount of ice stored on the continents. In general, the wind velocity of the winter trades has been intensified during cool climatic stages of the earth (δ¹⁸0 stages 4 and 2) and they have been relaxed during warm stages (δ¹⁸0 stages 3 and 1). Seasonal contrast of the trade winds has also fluctuated within time, having been relatively high during the upper part of δ¹⁸0 stage 3.     

The Dawson Cut Forest Bed in the Fairbanks area, Alaska, is about two million years old☆

The Dawson Cut Forest Bed lies in the lower part of thick, late Cenozoic loess deposits in the Fairbanks area. It is associated with several distal tephra beds that provide age control and offer the opportunity of its recognition elsewhere in central Alaska. EC tephra (named herein) occurs in the uppermost part of the Dawson Cut Forest Bed and its petrographic and chemical properties point to a co-magmatic relationship with PA tephra, which has not been found in direct association with the forest bed. Both tephra beds are pink and have unusually high Cl in their glass shards, which readily separates them from all other tephra beds in the Fairbanks area. They were produced by discrete eruptions, closely spaced in time. PA tephra has a glass-fission-track age of 2.02 ± 0.14 myr, indicating that the Dawson Cut Forest Bed must be about 2 million years old. The Palisades tephra (named herein) has very similar properties to these two tephra beds, suggesting that the buried forest bed just above it at the Palisades site on the Yukon River, about 250 km west of Fairbanks, correlates with the Dawson Cut Forest Bed.     

青藏高原40年来降水量时空变化趋势

利用青藏高原及其周边的97个气象台站1961-2000年的逐月降水量资料,以Mann-Kendall趋势检验方法结合主成分分析,分析了高原降水量40年来的时间变化趋势和空间演变特征,并探讨了其变化趋势的区域分异因素和可能的趋势突变时间。结果发现,高原站点40年年降水量大部分表现为增大趋势,只有青海东南部和南疆及西藏部分站点出现减小趋势,这一时间变化趋势的空间分布则大致表现为高原中东部和南北的反向变化,同时高原站点冷季降水量增大趋势明显;以冷、暖季降水量的第三载荷向量场分异进行的分区在10年年代际和逐年变化中都有明显体现,暖季区域分异因素的主要相关区域降水量变化趋势与相应载荷向量的时间变化相关显著,在1989年出现可信的突变时间点;冷季相关分异区域的降水量变化趋势在1987年出现明显增大趋势,同样与载荷向量的时间变化相关显著,但突变时间点检测不明显。     

Uranium-series disequilibrium in tuffs from Yucca Mountain,Nevada, as evidence of pore-fluid flow over the last million years

Samples of tuff from boreholes drilled into fault zones in the Exploratory Studies Facility (ESF) and relatively unfractured rock of the Cross Drift tunnels, at Yucca Mountain, Nevada, have been analysed by U-series methods. This work is part of a project to verify the finding of fast flow-paths through the tuff to ESF level, indicated by the presence of ‘bomb’ ³⁶Cl in pore fluids. Secular radioactive equilibrium in the U decay series, (i.e. when the radioactivity ratios ²³⁴U/²³⁸U, ²³⁰Th /²³⁴U and ²²⁶Ra/²³⁰Th all equal 1.00) might be expected if the tuff samples have not experienced radionuclide loss due to rock-water interaction occurring within the last million years. However, most fractured and unfractured samples were found to have a small deficiency of ²³⁴U (weighted mean ²³⁴U/²³⁸U=0.95±0.01) and a small excess of ²³⁰Th (weighted mean ²³⁰Th/²³⁴U 1.10±0.02). The ²²⁶Ra/²³⁰Th ratios are close to secular equilibrium (weighted mean=0.94±0.07). These data indicate that ²³⁴U has been removed from the rock samples in the last ∼350 ka, probably by pore fluids. Within the precision of the measurement, it would appear that ²²⁶Ra has not been mobilized and removed from the tuff, although there may be some localised ²²⁶Ra redistribution as suggested by a few ratio values that are significantly different from 1.0. Because both fractured and unfractured tuffs show approximately the same deficiency of ²³⁴U, this indicates that pore fluids are moving equally through fractured and unfractured rock. More importantly, fractured rock appears not to be a dominant pathway for groundwater flow (otherwise the ratio would be more strongly affected and the Th and Ra isotopic ratios would likely also show disequilibrium). Application of a simple mass-balance model suggests that surface infiltration rate is over an order of magnitude greater than the rate indicated by other infiltration models and that residence time of pore fluids at ESF level is about 400 a. Processes of U sorption, precipitation and re-solution are believed to be occurring and would account for these anomalous results but have not been included in the model. Despite the difficulties, the U-series data suggest that fractured rock, specifically the Sundance and Drill Hole Wash faults, are not preferred flow paths for groundwater flowing through the Topopah Spring tuff and, by implication, rapid-flow, within 50 a, from the surface to the level of the ESF is improbable.     

Constraints on ice-sheet thickness over Tibet during the last 40 000 years

The Late Quaternary glaciation of Tibet has received considerable attention in the last few decades due to its influence on the regional climate, especially the Asian summer monsoon. Recently, however, it has been argued that the Tibetan ice sheet also might have played an important role in initiating global-scale palaeoclimatic changes. Controversy, however, exists on the nature of Late Quaternary ice cover over Tibet due largely to the subjectivity in the interpretation of the sparse and complex geomorphological evidence. We have examined this problem in the light of δ ¹⁸O data (a temperature proxy) of ice cores-from the Dunde ice cap on the northern flank of Tibet. Considering only the gross features in the Dunde ice-core isotopic data, we have interpreted a temperature decrease of 4°–6°C and consequent lowering of equilibrium line altitude (ELA) in the range 700–850 m during the last glacial stage (LGS). This could have caused depression of the snow line below the mean altitutde of the Tibetan plateau, resulting in an areally extensive but marginally thick ice cover. However, if one also considers the possibility that precipitation on the Tibetan plateau during LGS may have been significantly lower than at present, the ELA depression would be much less than that estimated by considering the temperature effect alone.     

Model for paleoceanographic reconstructions of the California Current during the last 8000 years

Two independent data sets are used to develop a model for reconstructing sea-surface temperature and dynamic height anomaly distributions for the California Current during the last 8000 years. The first data set, all hydrographic data available for the California Current region, was used to determine the statistical relationships between the historical record of sea-surface conditions in the Santa Barbara Basin and all one-degree-square grid points of the California Current area. Given these relationships and the second data set, an 8000-year record of sea-surface temperatures and dynamic height anomalies from the Santa Barbara Basin, past sea-surface conditions throughout the California Current can be estimated for times before historical observations. The 8000-year record of sea-surface conditions was estimated by analysis of the radiolarian fauna found in a varved sediment core from the Santa Barbara Basin (Pisias, N. G., 1978, Quaternary Research 10, 366–384). The reconstructions of sea-surface temperature and dynamic height anomalies indicate that at times of cold sea-surface conditions in the Santa Barbara Basin, the flow of the California Current was much stronger than it is today or was during the times of the warmest sea-surface conditions during the last 8000 years. The atmospheric circulation during the winter of 1950, the period of the coldest recorded sea-surface temperatures in the Santa Barbara Basin based on the historical data set, contained a strong northerly component in the winds which is consistent with the inferred increase in the California Current at times of cold sea-surface temperatures. Times of warm sea-surface temperatures in the Santa Barbara Basin are characterized by decreased southward flow of the California Current and a marked increase in northward flow into the Santa Barbara Basin itself. In the historical record, times of warm sea-surface temperatures are often associated with high precipitation in southern California. The atmospheric circulation during the winter of 1968–1969 is characterized by strong eastward flow over southern California and a northward transport of warm humid air from the tropics into the region of southern California producing the high rainfall observed. The persistence of this atmospheric circulation could produce the more zonal flow predicted for the California Current during times of warmer average conditions in the Santa Barbara Basin.     

Fluctuation in the level of pluvial Lake Lahontan during the last 40,000 years

Samples of algal tufa, gastropods and calcite-cemented sand were collected from the Walker and Pyramid Lake areas of the Lahontan Basin, Nevada. X-ray diffraction petrographic and radiocarbon analyses show that massive forms of tufa such as the dendritic variety contain secondary carbon-bearing material and therefore yield unreliable radiocarbon dates. Dense coating of tufa (lithoid), however, gave radiocarbon ages in agreement with dates on coexisting aragonite gastropods. Radiocarbon data from the study were combined with previously dated noncarbonate materials [Born, S. M. (1972). “Lake Quaternary History, Deltaic Sedimentation, and Mudlump Formation at Pyramid Lake, Nevada”, Center for Water Resources, Desert Research Inst., Reno, Nevada] to give an internally consistent record of lake level fluctuations for the past 40,000 years. The main features of the Lahontan chronology are (1) extreme high stands (1330 m above sea level) 13,500 to 11,000 and 25,000 to 22,000 B.P., (2) a moderate high stand (1260 m above sea level) 20,000 to 15,000 B.P., (3) a low stand of unknown elevation 40,000 to 25,000 B.P., (4) an extremely low stand 9000 to 5000 B.P., and (5) an overall increase in the size of Walker and Pyramid Lakes during the past 5000 years, until the late 19th century. Pore fluid data indicate that Walker Lake desiccated sometime during the period 9050 to 6400 B.P. Salts deposited as a result of this dessication are still undergoing dissolution causing a flux of chloride, carbon, and other solute species from the sediments to the overlying lake water. Pore fluid data obtained from Pyramid Lake sediments do not indicate the presence of a concentrated brine at depth. This suggests that Pyramid Lake did not dry completely during this period although it may have been severely reduced in size. There has been considerable disagreement regarding the occurrence of extreme arid conditions (altithermal period) since 10,000 B.P. [Mehringer, P. J. (1977). “Models and Great Basin Prehistory”. Desert Research Inst. Pub, Reno, Nevada]. The data of this study suggest that such a climatic regime did occur in the western Great Basin during the period 9000 to 5000 B.P.     

Palaeoceanographic change in the northeastern South China Sea during the last 15000 years

The sedimentary succession of piston core RC26-16, dated by ¹⁴C accelerator mass spectrometry, provides a nearly continuous palaeoceanographic record of the northeastern South China Sea for the last 15000 yr. Planktic foraminiferal assemblages indicate that winter sea-surface temperatures (SSTs) rose from 18°C to about 24°C from the last glacial to the Holocene. A short-lived cooling of 1°C in winter temperature centred at about 11000 ¹⁴C yr ago may reflect the Younger Dryas cooling event in this area. Summer SSTs have remained between 27°C and 29°C throughout the record. The temperature difference between summer and winter was about ca. 9°C during the last glacial, much higher than the Holocene value of ca. 5°C. During the late Holocene a short-lived cooling event occurred at about 4000 ¹⁴C yr ago. Oxygen and carbon isotopic gradients between surface (0–50 m) and subsurface (50–100 m) waters were smaller during the last glacial than those in the Holocene. The fluctuation in the isotopic gradients are caused most likely by changes in upwelling intensity. Smaller gradients indicate stronger upwelling during the glacial winter monsoon. The fauna-derived estimates of nutrient content of the surface waters indicate that the upwelling induced higher fertility and biological productivity during the glacial. The winter monsoon became weaker during the Holocene. The carbonate compensation depth and foraminiferal lysocline were shallower during the Holocene, except for a short-lived deepening at about 5000 ¹⁴C yr ago. A preservation peak of planktic foraminifera and calcium carbonate occurred between 13400 and 12000 ¹⁴C yr ago, synchronous to the global preservation event of Termination I.     

Response of glacier area variation to climate change in the Kaidu-Kongque river basin,Southern Tianshan Mountains during the last 20 years

Glaciers are crucial water resources for arid inland rivers in Northwest China. In recent decades, glaciers are largely experiencing shrinkage under the climate-warming scenario, thereby exerting tremendous influences on regional water resources. The primary role of understudying watershed scale glacier changes under changing climatic conditions is to ensure sustainable utilization of regional water resources, to prevent and mitigate glacier-related disasters. This study maps the current (2020) distribution of glacier boundaries across the Kaidu-Kongque river basin, south slope of Tianshan Mountains, and monitors the spatial evolution of glaciers over five time periods from 2000–2020 through thresholded band ratios approach, using 25 Landsat images at 30 m resolution. In addition, this study attempts to understand the role of climate characteristics for variable response of glacier area. The results show that the total area of glaciers was 398.21 km² in 2020. The glaciers retreated by about 1.17 km²/a (0.26%/a) from 2000 to 2020. The glaciers were reducing at a significantly rapid rate between 2000 and 2005, a slow rate from 2005 to 2015, and an accelerated rate during 2015–2020. The meteorological data shows slight increasing trends of mean annual temperature (0.02°C/a) and annual precipitation (2.07 mm/a). The correlation analysis demonstrates that the role of temperature presents more significant correlation with glacier recession than precipitation. There is a temporal hysteresis in the response of glacier change to climate change. Increasing trend of temperature in summer proves to be the driving force behind the Kaidu-Kongque basin glacier recession during the recent 20 years.© 2021 China Geology Editorial Office.