Late Tertiary tectonic evolution of northern Iran: A case for simple crustal folding

Here we present a crustal folding or buckling mechanism to explain the rootless 3–5 km high Alborz Mountains in northern Iran as well as  10 km of Late Miocene to recent subsidence in the south Caspian basin and  3–6 km of subsidence in the central Iranian basin in the context of the middle Miocene to recent Arabia–Eurasia collision. A key element of the mechanism is the presence of lateral and vertical lithospheric strength contrasts between the north Iranian continental and south Caspian oceanic crusts: when compression from the collision is applied across the region, the strong south Caspian oceanic crust, buried under > 10 km of premiddle Miocene sediment, interacts with the bottom of the mechanically strong continental upper crust of northern Iran, resulting in upward buckling of the continental crust and downward buckling of the oceanic crust. We test this mechanism using a finite-element numerical model with a Maxwell rheology and obtain results that are consistent with the geological and geophysical observations. The observations compiled here and the model results demonstrate the potential for using this region as a natural laboratory for studying the early stages of continent–oceanic collision, including processes like basin inversion, fault localization and, potentially, subduction initiation.     

A peat core based estimate of Late-glacial and Holocene methane emissions from northern peatlands

We estimate the intensity of Late-glacial and Holocene methane emissions from peatlands based on their paleo net primary production (PNPP). The PNPP is derived from the carbon accumulation rates of the studied bog profile (Etang de la Gruère, Switzerland), which are corrected for the degree of peat degradation. The obtained PNPP curve is taken as a proxy for methane emissions. It shows relatively high values (90 g C m^− 2 yr^− 1) early in the Bolling/Allerod and drops to low values (40 g C m^− 2 yr^− 1) during the Younger Dryas cold period. With the onset of the Holocene the PNPP increases strongly up to 150 g C m^− 2 yr^− 1 around ca. 10,000 Cal. yr bp. This is followed by a decline to minimum values (30 to 40 g C m^− 2 yr^− 1) between 6500 and 4000 Cal. yr bp. Thereafter, the PNPP starts to increase again to reach its highest value (175 g C m^− 2 yr^− 1) around 1000 Cal. yr bp.The PNPP curve correlates well with the evolution of the atmospheric methane concentrations as derived from Greenland ice-cores. For example, minima in atmospheric methane reported during the Younger Dryas and around 5200 Cal. yr bp are coinciding with the lowest values of PNPP and the negative atmospheric methane peak at 8200 Cal. yr bp corresponds to a marked decrease in PNPP.Our PNPP curve suggests that the methane emissions from northern peatlands evolved similar to those of low latitude wetlands and together they largely determined the evolution of atmospheric methane throughout the Late-glacial and the Holocene. The abruptness of the rise of atmospheric methane at the end of the Younger Dryas probably points to an additional source (e.g. marine gas hydrates), but very early in the Holocene the peatlands have likely become the dominant source of atmospheric methane.     

Cryospheric change in China

This paper provides an overview of the current status of the cryosphere in China and its changes. Up-to-date statistics of the cryosphere in China are summarized based on the latest available data. There are 46,377 glaciers in China, covering an area of 59,425 km². The glacier ice reserve is estimated to be about 5600 km³ and the annual glacier runoff is about 61.6 × 10⁹ m³. The continuous snow cover extent (> 60 days) in China is about 3.4 × 10⁶ km² and the maximum water equivalent is 95.9 × 10⁹ m³ yr^− 1. The permafrost area in China is about 1.72 × 10⁶ km². The total ground ice reserve on the Qinghai–Tibetan Plateau is estimated to be about 10,923 km³. Recent investigations indicated that glacier areas in China have shrunk about 2–10% over the past 45 yr. Total glacier area has receded by about 5.5%. Snow mass has increased slightly. Permafrost is clearly degrading, as indicated by shrinking areas of permafrost, increasing depth of the active layer, rising of lower limit of permafrost, and thinning of the seasonal frost depth. Some models predict that glacier area shrinkage could be as high as 26.7% in 2050, with glacier runoff increasing until its maximum in about 2030. Although snow mass shows an increasing trend in western China, in eastern China the trend is toward decreasing snow mass, with increasing interannual fluctuations. Permafrost degradation is likely to continue, with one-third to one-half of the permafrost on the Qinghai–Tibetan Plateau anticipated to degrade by 2100. Most of the high-temperature permafrost will disappear by then. The permafrost in northeastern China will retreat further northward.     

Quaternary uplift of northern England

Upland flats, attributable to erosion, have long been recognised in the landscape of the Lake District region of NW England, at altitudes of up to ~ 800 m O.D. Extrapolation using uplift rates derived from dated Pleistocene sites (karstic caves and other features) in the adjacent Pennine uplands suggests that if this succession of flats formed close to sea-level they date from the Middle Pliocene onwards, indicating a subsequent time-averaged uplift rate of almost 0.3 mm a^− 1. Numerical modelling indicates that erosion of surrounding areas at a typical rate of 0.2 mm a^− 1 since 3.1 Ma could have caused this uplift, as well as constraining the local effective viscosity of the lower crust as ~ 4 × 10¹⁸ Pa s and the typical local Moho temperature as ~ 650 °C. It is thus feasible that most of the topography of northern England has developed since the Middle Pliocene, as a consequence of coupling between erosion and the resulting induced flow in the lower continental crust. The much faster vertical crustal motions indicated in this part of northern England, compared with SE England, are thus mainly a consequence of much greater mobility of the lower crust in the north, due to its younger thermal age and the heating effect of radioactive Palaeozoic granites. Uplift of this magnitude, which has previously gone unrecognised, may have affected post-Pliocene global climate.     

Dynamic topography of the East European craton: Shedding light upon lithospheric structure, composition and mantle dynamics

Most of the East European Craton lacks surface relief; however, the amplitude of topography at the top of the basement exceeds 20 km, the amplitude of topography undulations at the crustal base reaches almost 30 km with an amazing amplitude of ca. 50 km in variation in the thickness of the crystalline crust, and the amplitude of topography variations at the lithosphere–asthenosphere boundary exceeds 200 km. This paper examines the relative contributions of the crust, the subcrustal lithosphere, and the dynamic support of the sublithospheric mantle to maintain surface topography, using regional seismic data on the structure of the crystalline crust and the sedimentary cover, and thermal and large-scale P- and S-wave seismic tomography data on the structure of the lithospheric mantle. For the Precambrian lithosphere, an analysis of Vp/Vs ratio at 100, 150, 200, and 250 km depths does not show any age-dependence, suggesting that while Vp/Vs ratio can be effectively used to outline the cratonic margins, it is not sensitive to compositional variations within the cratonic lithosphere.Statistical analysis of age-dependence of velocity, density, and thermal structure of the continental crust and subcrustal lithosphere in the study area (0–62E, 45–72N) allows to link lithospheric structure with the tectonic evolution of the region since the Archean. Crustal thickness decreases systematically with age from 42–44 km in regions older than 1.6 Ga to 37–40 km in the Paleozoic–Mesoproterozoic structures, and to ca. 31 km in the Meso-Cenozoic regions. However, the isostatic contribution of the crust to the surface topography of the East European Craton is almost independent of age (ca. 4.5 km) due to an interplay of age-dependent crustal and sedimentary thicknesses and lithospheric temperatures.On the contrary, the contribution of the subcrustal lithosphere to the surface topography strongly depends on the age, being slightly positive (+ 0.3 + 0.7 km) for the regions older than 1.6 Ga and negative (− 0.5–1 km) for younger structures. This leads to age-dependent variations in the residual topography, i.e. the topography which cannot be explained by the assumed thermal and density structure of the lithosphere, and which can (at least partly) originate from the dynamic component caused by the mantle flow. Positive dynamic topography at the cratonic margins, which exceeds 2 km in the Norwegian Caledonides and in the Urals, clearly links their on-going uplift with deep mantle processes. Negative residual topography beneath the Archean-Paleoproterozoic cratons (− 1–2 km) indicates either a smaller density deficit (ca. 0.9%) in their subcrustal lithosphere than predicted by global petrologic data on mantle-derived xenoliths or the presence of a strong convective downwelling in the mantle. Such mantle downflows can effectively divert heat from the lithospheric base, leading to a long-term survival of the Archean-Paleoproterozoic lithosphere.     

Climate changes and recent glacier behaviour in the Chilean Lake District

Atmospheric temperatures measured at the Chilean Lake District (38°–42°S) showed contrasting trends during the second half of the 20th century. The surface cooling detected at several meteorological stations ranged from − 0.014 to − 0.021 °C a^− 1, whilst upper troposphere (850–300 gpm) records at radiosonde of Puerto Montt (41°26′S/73°07′W) revealed warming between 0.019 and 0.031 °C a^− 1. Regional rainfall data collected from 1961 to 2000 showed the overall decrease with a maximum rate of − 15 mm a^− 2 at Valdivia st. (39°38′S/73°05′W). These ongoing climatic changes, especially the precipitation reduction, seem to be related to El Niño–Southern Oscillation (ENSO) phenomena which has been more frequent after 1976. Glaciers within the Chilean Lake District have significantly retreated during recent decades, in an apparent out-of-phase response to the regional surface cooling. Moreover, very little is known about upper troposphere changes and how they can enhance the glacier responses. In order to analyse their behaviour in the context of the observed climate changes, Casa Pangue glacier (41°08′S/71°52′W) has been selected and studied by comparing Digital Elevation Models (DEMs) computed at three different dates throughout the last four decades. This approach allowed the determination of ice elevation changes between 1961 and 1998, yielding a mean thinning rate of − 2.3 ± 0.6 m a^− 1. Strikingly, when ice thinning is computed for the period between 1981 and 1998, the resulting rate is 50% higher (− 3.6 ± 0.6 m a^− 1). This enhanced trend and the related area loss and frontal retreat suggests that Casa Pangue might currently be suffering negative mass balances in response to the upper troposphere warming and decreased precipitation of the last 25–30 yr, as well as debris cover would not prevent the glacier from a fast reaction to climate forcing. Most of recent glaciological studies regarding Andean glaciers have concentrated on low altitude changes, namely frontal variations, however, in order to better understand the regional glacier changes, new data are necessary, especially from the accumulation areas.     

Active crustal deformation beyond the SE margin of the Tibetan Plateau: Constraints from the evolution of fluvial systems

An integrated explanation is proposed for the Late Cenozoic crustal deformation in Yunnan, SW China, using sedimentary and geomorphological evidence from the Yangtze and Red River systems. The observed fluvial incision indicates up to ~ 15 km of crustal thickening, associated with ~ 3 km of uplift, apparently triggered at ~ 8 Ma by monsoon-induced erosion drawing mobile lower crust from beneath Tibet to the northwest. The mobile lower-crustal layer beneath Yunnan was initially very thin, but a positive feedback loop developed, whereby each incremental influx of lower-crust widened and heated this layer, facilitating the next increment. At ~ 5 Ma, the shear tractions exerted on the brittle upper-crust by this flowing lower crust became sufficient to reactivate pre-existing lines of weakness, dragging blocks of the brittle layer southward and creating the region′s modern active fault systems. This region thus provides a dramatic example of crustal deformation induced by Late Cenozoic climate change, notwithstanding its location adjoining the India–Eurasia plate boundary.     

Are sea-level-rise trends along the coasts of the north Indian Ocean consistent with global estimates?

Mean-sea-level data from coastal tide gauges in the north Indian Ocean were used to show that low-frequency variability is consistent among the stations in the basin. Statistically significant trends obtained from records longer than 40 years yielded sea-level-rise estimates between 1.06–1.75 mm yr^− 1, with a regional average of 1.29 mm yr^− 1, when corrected for global isostatic adjustment (GIA) using model data. These estimates are consistent with the 1–2 mm yr^− 1 global sea-level-rise estimates reported by the Intergovernmental Panel on Climate Change.     

Mid Pleistocene climate instability in the subtropical northwestern Atlantic

We present Globigerinoides ruber, G. sacculifer and Neogloboquadrina dutertrei oxygen isotope records from northwestern subtropical Atlantic Site 1058 spanning the mid Pleistocene ( 600 to 400 ka). The high temporal resolution of these records ( 800 yr) allows us to compare millennial-scale climate signals during one of the most extreme glacial periods of the Pleistocene (Marine Isotope Stage (MIS) 12) to an earlier, less extreme glacial (MIS 14), as well as to two full interglacial intervals (MIS 13 and MIS 15). We observe excellent agreement in the timing and amplitude of variations between the surface-most dwelling species G. ruber and Northern Hemisphere insolation during the two interglacial periods. There is some expression of Northern Hemisphere insolation during glacial MIS 14; however, during the more extreme glacial MIS 12 Northern Hemisphere insolation patterns are not apparent in any of the planktonic foraminiferal δ¹⁸O records. Insolation remains relatively high, but δ¹⁸O values increase toward the characteristic δ¹⁸O maximum of MIS 12 in all three of the records. On the millennial-scale, all three species display their highest amplitude δ¹⁸O variations (with a period between 4–6 kyr) during glacial MIS 12. Suborbital-scale variability is also statistically significant during glacial MIS 14, but the amplitude is smaller. These results support hypotheses linking millennial-scale climate fluctuations to the extent of continental glaciation. We propose that the relatively high degree of sea surface instability during one of the most extreme glacial periods of the Pleistocene arises from the competing effects of strong atmospheric winds related to the presence of a large ice sheet to the north and persistently high incident solar radiation during this interval of time.     

Recent evolution and mass balance of Cordón Martial glaciers, Cordillera Fueguina Oriental

Past and present glacier changes have been studied at Cordón Martial, Cordillera Fueguina Oriental, Tierra del Fuego, providing novel data for the Holocene deglaciation history of southern South America and extrapolating as well its future behavior based on predicted climatic changes. Regional geomorphologic and stratigraphic correlations indicate that the last glacier advance deposited the ice-proximal (“internal”) moraines of Cordón Martial, around 330 ¹⁴C yr BP, during the Late Little Ice Age (LLIA). Since then glaciers have receded slowly, until 60 years ago, when major glacier retreat started. There is a good correspondence for the past 100 years between the surface area variation of four small cirque glaciers at Cordón Martial and the annual temperature and precipitation data of Ushuaia. Between 1984 and 1998, Martial Este Glacier lost 0.64 ± 0.02 × 10⁶ m³ of ice mass (0.59 ± 0.02 × 10⁶ m³ w.e.), corresponding to an average ice thinning of 7.0 ± 0.2 m (6.4 ± 0.2 m w.e), according to repeated topographic mapping. More detailed climatic data have been obtained since 1998 at the Martial Este Glacier, including air temperature, humidity and solar radiation. These records, together with the monthly mass balance measured since March 2000, document the annual response of the Martial Este Glacier to the climate variation. Mass balances during hydrological years were positive in 2000, negative in 2001 and near equilibrium in 2002. Finally, using these data and the regional temperature trend projections, modeled for different future scenarios by the Atmosphere-Ocean Model (GISS-NASA/GSFC), potential climatic-change effects on this mountain glacier were extrapolated. The analysis shows that only the Martial Este Glacier may survive this century.     

UBVRI photometry of the eclipsing symbiotic star CI Cyg during 1991–1996

UBVRI observations of the eclipsing symbiotic star CI Cyg made during 1991–1995 are analyzed, the results of which indicate that the system is in the same quiescent state as during 1985–1990. Variations in the "extra-eclipse" color index U-B with an amplitude of about 0^m.3–0^m.4 and a characteristic time of about 9–10 yr have been detected for the first time. A more thorough analysis of all available data in the UBVRI bands of the spectrum will be required to investigate the existence of possible long-period variations in the other color indices and brightness of the system.Translated fromAstrofizika, Vol. 39, No. 2, pp. 211–216, April–June, 1996.     

Reconstruction of Mediterranean sea level fields for the period 1945–2000

The distribution of sea level in the Mediterranean Sea is recovered for the period 1945–2000 by using a reduced space optimal interpolation analysis. The method involves estimating empirical orthogonal functions from satellite altimeter data spanning the period 1993–2005 that are then combined with tide gauge data to recover sea level fields over the period 1945–2000. The reconstruction technique is discussed and its robustness is checked through different tests. For the altimetric period (1993–2000) the prediction skill is quantified over the whole domain by comparing the reconstructed fields with satellite altimeter observations. For past times the skill can only be tested locally, by validating the reconstruction against independent tide gauge records. The reconstructed distribution of sea level trends for the period 1945–2000 shows a positive peak in the Ionian Sea (up to 1.5 mm yr^− 1) and a negative peak of − 0.5 mm yr^− 1 in a small area to the south-east of Crete. Positive trends are found nearly everywhere, being larger in the western Mediterranean (between 0.5 and 1 mm yr^− 1) than in the eastern Mediterranean (between 0 and 0.5 mm yr^− 1). The estimated rate of mean sea level rise for the period 1945–2000 is 0.7 ± 0.2 mm yr^− 1, i.e. about a half of the rate estimated for global mean sea level. These overall results do not appear to be very sensitive to the distribution of tide gauges. The poorest results are obtained in open-sea regions with intense mesoscale variability not correlated with any tide gauge station, such as the Algerian Basin.     

Long-chain alkenones and U37′ variability along a south–north transect in the Western Pacific Ocean

In this paper, we present 50 surface water samples collected during the IMAGES III cruise (June–July 1997) along a transect from New Zealand to the China Sea (42°S–178°E, 21°N–120°E) covering a temperature range from 13.3 to 30.4 °C. A very worthwhile aspect of this study is a coupling of both biomarker (alkenone) and coccolithophorid counting. We show that the U₃₇^k′–temperature relationship is very similar to the Prahl et al. [Geochim. Cosmochim. Acta 52 (1988) 2203] culture calibration and to the global core top calibration of Müller et al. [Geochim. Cosmochim. Acta 62 (1998) 1757]. However, in the warmest surface waters of the Western Pacific ocean (>26.4 °C) where Gephyrocapsa oceanica is likely the most widespread species, the associated U₃₇^k′ has a constant value of 1.0. The consequence is that above this temperature threshold, U₃₇^k′ cannot be used as an accurate paleothermometer.     

Temperature changes on the Tibetan Plateau during the past 600 years inferred from ice cores and tree rings

The climatological signal of δ¹⁸O variations preserved in ice cores recovered from Puruogangri ice field in the central Tibetan Plateau (TP) was calibrated with regional meteorological data for the past 50 years. For the period AD 1860–2000, 5-yearly averaged ice core δ¹⁸O and a summer temperature reconstruction derived from pollen data from the same ice core were compared. The statistical results provide compelling evidence that Puruogangri ice core δ¹⁸O variations represent summer temperature changes for the central TP, and hence regional temperature history during the past 600 years was revealed. A comparison of Puruogangri ice core δ¹⁸О with several other temperature reconstructions shows that broad-scale climate anomalies since the Little Ice Age occurred synchronously across the eastern and southern TP, and the Himalayas. Common cold periods were identified in the 15th century, 1625–1645 AD, 1660–1700 AD, 1725–1775 AD, 1795–1830 AD, 1850–1870 AD, 1890–1920 AD, 1940–1950 AD, and 1975–1985 AD. The period 1725–1775 AD was one of the most prolonged cool periods during the past 400 years and corresponded to maximum Little Ice Age glacier advance of monsoonal temperate glaciers of the TP.     

Sea level forcing in the Mediterranean Sea between 1960 and 2000

Sea level trends and inter-annual variability in the Mediterranean Sea for the period 1960–2000 is explored by comparing observations from tide gauges with sea level hindcasts from a barotropic 2D circulation model, and two full primitive equation 3D ocean circulation models, a regional one and the Mediterranean component of a global one,. In the 2D model, 50% of the sea level variance was found to result from the wind and atmospheric pressure forcing. In the 3D models, 20% of the sea level variance was explained by the steric effects. The sea level residuals at the tide gauges locations, calculated by subtraction of the 2D model output from the sea level observations are significantly correlated (r = 0.4) with the steric signals from the 3D models. After the removal of the atmospheric and the steric contributions the tide-gauge sea level records indicate a period where sea level was stable (1960–1975) and a period where sea level was rising (1975–2000) with rates in the range 1.1–1.8 mm/yr. A part of the residual trend can be explained by the contribution of local land movements (0.3 mm/yr) while its major part indicates a global signal, probably mass addition, appearing after 1975.     

Morphological and ice-dynamical changes on the Tasman Glacier, New Zealand, 1990–2007

This paper presents data concerning recent (1990–2007) surface morphological and ice-dynamical changes on the Tasman Glacier, New Zealand. We use remote-sensing data to derive rates of lake growth, glacier velocities and rates of glacier surface lowering. Between 1990 and 2007, the glacier terminus receded ~ 3.5 km and a large ice-contact proglacial lake developed behind the outwash head. By 2007 the lake area was ~ 6 km² and had replaced the majority of the lowermost 4 km of the glacier tongue. There is evidence that lake growth is proceeding at increasing rates — the lake area doubled between 2000 and 2007 alone. Measured horizontal glacier velocities decline from 150 m a^− 1 in the upper glacier catchment to almost zero at the glacier terminus and there is a consequent down-glacier increase in surface debris cover. Surface debris mapping shows that a large catastrophic rockfall onto the glacier surface in 1991 is still evident as a series of arcuate debris ridges below the Hochstetter icefall. Calculated glacier surface lowering is most clearly pronounced around the terminal area of the glacier tongue, with down-wasting rates of 4.2 ± 1.4 m a^− 1 in areas adjacent to the lateral moraine ridges outside of the current lake extent. Surface lowering rates of approximately 1.9 ± 1.4 m a^− 1 are common in the upper areas of the glacier. Calculations of future lake expansion are dependent on accurate bathymetric and bed topography surveys, but published data indicate that a further 8–10 km of the glacier is susceptible to calving and further lake development in the future.     

Geodetic GPS-based analysis of recent crustal motions in Victoria Land (Antarctica)

GPS measurements were collected within the framework of the VLNDEF (Victoria Land Network for DEFormation control) project, which was started in 1999 with the aim of detecting crustal deformation in Northern Victoria Land (Antarctica). The network was established in 1999 and is composed of one permanent station (TNB1), which has been observing since 1998, and 28 periodically surveyed control points. Three complete campaigns and some partial surveying of the network have been carried out to date.Data processing and analysis have been performed using an undifferenced approach for the network position within the ITRF. A double-differences-based strategy has been applied for movement detection. The data processing and analysis of results have been carried out for all available data, both periodically acquired and long time series.GPS measurements collected between December 1999 and February 2006 indicate a mean “absolute” motion of the region of v_e = 11.3 mm/yr and v_n = − 11.1 mm/yr and rock uplift rates of v_u = 2.8 mm/yr. These values are consistent with Antarctic plate motion and the general postglacial rebound models of the region. The relative motions within VLNDEF are small and only few points show velocities greater than the confidence levels.     

Estimation of net accumulation rate at a Patagonian glacier by ice core analyses using snow algae

Snow algae in a 45.97-m-long ice core from the Tyndall Glacier (50°59′05″S, 73°31′12″W, 1756 m a.s.l.) in the Southern Patagonian Icefield were examined for potential use in ice core dating and estimation of the net accumulation rate. The core was subjected to visual stratigraphic observation and bulk density measurements in the field, and later to analyses of snow algal biomass, water isotopes (¹⁸O, D), and major dissolved ions. The ice core contained many algal cells that belonged to two species of snow algae growing in the snow near the surface: Chloromonas sp. and an unknown green algal species. Algal biomass and major dissolved ions (Na⁺, K⁺, Mg²⁺, Ca²⁺, Cl⁻, SO₄²⁻) exhibited rapid decreases in the upper 3 m, probably owing to melt water elution and/or decomposition of algal cells. However, seasonal cycles were still found for the snow algal biomass, ¹⁸O, D-excess, and major ions, although the amplitudes of the cycles decreased with depth. Supposing that the layers with almost no snow algae were the winter layers without the melt water essential to algal growth, we estimated that the net accumulation rate at this location was 12.9 m a^− 1 from winter 1998 to winter 1999, and 5.1 m from the beginning of winter to December 1999. These estimates are similar to the values estimated from the peaks of ¹⁸O (17.8 m a^− 1 from summer 1998 to summer 1999 and 11.0 m from summer to December 1999) and those of D-excess (14.7 m a^− 1 from fall 1998 to fall 1999 and 8.6 m a^− 1 from fall to December 1999). These values are much higher than those obtained by past ice core studies in Patagonia, but are of the same order of magnitude as those predicted from various observations at ablation areas of Patagonian glaciers.     

Late Pleistocene and Holocene environments in the Nile basin

Owing to the very gently sloping nature of the flood plain in the lower White Nile valley, which is underlain by a former lake-bed, the depositional record in that area is unusually well preserved. In Egypt and along the Blue Nile phases of erosion have destroyed segments of the sedimentary record, but the White Nile sequence is a good proxy for both the main Nile and the Blue Nile. During the last 15 ka, at least, times of high flow in the Blue Nile and main Nile were synchronous with those in the White Nile.Not all the White Nile flood deposits have been preserved but calibrated radiocarbon dates obtained on fossil freshwater and amphibious Pila shells and fish bones indicate that White Nile levels were high around 14.7–13.1 ka, 9.7–9.0 ka, 7.9–7.6 ka, 6.3 ka and 3.2–2.8 ka. The Blue Nile record is more fragmentary and that of the main Nile even more so except for the Holocene Nile delta. Calibrated radiocarbon ages for high Blue Nile flows indicate very high flood levels towards 13.9–13.2 ka, 8.6 ka, 7.7 ka and 6.3 ka.Incision by the Blue Nile and main Nile has caused progressive incision in the White Nile amounting to at least 4 m since the terminal Pleistocene  15 ka ago and at least 2 m over the past 9 ka. The Blue Nile seems to have cut down at least 10 m since  15 ka and at least 4 m since 9 ka. The time-transgressive and relatively late inception of plant domestication in the Nile valley may partly reflect this history of incision. Nile incision would propagate upstream into the White Nile valley, draining previously swampy areas along the valley floor, which would then become accessible to cultivation.     

Modeled impacts of changes in tundra snow thickness on ground thermal regime and heat flow to the atmosphere in Northernmost Alaska

Seasonal snow covers the tundra surface for up to nine months of each year on the Alaskan North Slope. Variations in the snow thickness could strongly influence the thermal regime of the underlying soil and permafrost, and the surface energy balance. The impacts of increases and decreases in the tundra snow thickness on the thermal regime of snow surface, active layer, and permafrost, and on the conductive heat flow to the atmosphere were investigated numerically, by using an improved surface energy balance approach based one-dimensional heat transfer model. The baseline inputs for the numerical model are mean daily meteorological data and surface albedos collected at Barrow, Alaska from 1995 through 1999. Based on a study for the long-term mean daily maximum and minimum snow thickness distributions at Barrow in the snow season of 1948 through 1997, a snow thickness factor was defined and five simulation cases were run for the snow season of 1997–1998 by changing the snow thickness factor. The modeled results indicate that changes in snow thickness have significant impacts on ground thermal regimes and conductive heat flow to the atmosphere. Decreasing the snow thickness by 50% led to the maximum ground temperature decrease of 1.48 °C at 0.29 m depth, and 0.72 °C at 3.0 m depth; the magnitude of the mean conductive heat flow to the atmosphere for December increase of 4.3 Wm^− 2. Increasing the snow thickness by 50% resulted in the maximum ground temperature increase of 1.44 °C at 0.29 m depth, and 0.66 °C at 3.0 m depth; the magnitude of the mean conductive heat flow to the atmosphere for December decrease of 1.57 W m^− 2. On an annual basis, variation in the snow thickness by 50%, the ground temperature variations of more than 0.25 °C occurred as deep as 8.0 m below the ground surface. The modeled results also show that changes in snow thickness have a relatively small influence on the snow surface temperature.