Impact of interacting bark structure and rainfall conditions on stemflow variability in a temperate beech-oak forest,central Germany

ABSTRACT

Trees concentrate rainfall to near-stem soils via stemflow. When canopy structures are organized appropriately, stemflow can even induce preferential flow through soils, transporting nutrients to biogeochemically active areas. Bark structure significantly affects stemflow, yet bark-stemflow studies are primarily qualitative. We used a LaserBark to compute bark microrelief (MR), ridge-to-furrow amplitude (R) and slope (S) metrics per American Society of Mechanical Engineering standards (ASME-B46.1–2009) for two morphologically contrasting species (Fagus sylvatica L. (European beech), Quercus robur L. (pendunculate oak)) under storm conditions with strong bark water storage capacity (BWSC) influence in central Germany. Smaller R and S for F. sylvatica significantly lowered BWSC, which strongly and inversely correlated to maximum funnelling ratios and permitted stemflow generation at lower rain magnitudes. Larger R and S values in Q. robur reduced funnelling, diminishing stemflow drainage for larger storms. Quercus robur funnelling and stemflow was more reliant on intermediate rain intensities and intermittency to maintain bark channel-dependent drainage pathways. Shelter provided by Q. robur’s ridged bark also appears to protect entrained water, lengthening mean intrastorm dry periods necessary to affect stemflow. Storm conditions where BWSC plays a major role in stemflow accounted for much of 2013’s rainfall at the nearest meteorological station (Wulferstedt).

Editor M.C. Acreman; Associate editor not assigned     

A method to map riparian exotic vegetation (Salix spp.) area to inform water resource management

Hydrological processes within riparian environments worldwide are impacted when introduced species invade. Monitoring and management at substantial expense, are subsequently required to combat deleterious effects on the environment and stream hydrology. Willow species (Salicaceae: Salix spp.) introduced into Australia have spread throughout many riparian systems causing adverse environmental impacts, with high rates of water extraction when located within stream beds (in‐stream willows) thus altering hydrology. Strategies exist to manage willows; however, simpler, more cost‐effective methods are required to map and monitor spatial and temporal distributions. A method is presented to discriminate willow stands from surrounding native riparian vegetation using a single, very high 2 m resolution multispectral WorldView‐2 satellite image. A combination of spectral bands ‘coastal blue’ (400–450 nm), ‘red’ (630–690 nm), ‘red edge’ (705–745 nm) and ‘near infrared2’ (860–1040 nm), minimum noise fraction transformation, median filtering and maximum likelihood supervised classification provided the highest discriminatory power within a 25 km² study area. Of the spectral bands, coastal blue, red edge and near infrared2 are new bands that are not available in other multispectral sensors. These bands proved critical to the success of discriminating willows from other land cover categories (overall accuracy of 97%). Stream channels were defined by incorporating a LiDAR‐derived digital elevation model to discriminate between willows on stream banks and within stream beds. Canopy area estimates of in‐stream willows, coupled with water savings estimates from willow removal, suggest removal of 10.4 ha of Salix fragilis canopy from within river channels in the study area will potentially return 41 ML year^?1 to the environment. The method presented improves our understanding of willow impacts on hydrology and aids decisions regarding willow removal for water resource management. Copyright © 2013 John Wiley & Sons, Ltd.     

Peat depth as a control on Sphagnum moisture stress during seasonal drought

Peatlands are globally important long-term sinks of carbon, however there is concern that enhanced peat decomposition and moss moisture stress due to climate change mediated drought will reduce moss productivity making these ecosystems vulnerable to carbon loss and associated long-term degradation. Peatlands are resilient to summer drought moss stress because of negative ecohydrological feedbacks that generally maintain a wet peat surface, but where feedbacks may be contingent on peat depth. We tested this ‘survival of the deepest’ hypothesis by examining water table (WT) position, near-surface moisture content, and soil water tension in peatlands that differ in size, peat depth, and catchment area during a summer drought. All shallow sites (<40 cm depth) lost their WT (i.e., the groundwater well was dry) for considerable time during the drought period. Near-surface soil water tension increased dramatically at shallow sites following WT loss, increasing ~5–7.5× greater at shallow sites compared to deep sites (≥40 cm depth). During a mid-summer drought intensive field survey, we found that 60–67% of plots at shallow sites exceeded a 100 mb tension threshold used to infer moss water stress. Unlike the shallow sites, tension typically did not exceed this 100 mb threshold at the deep sites. Using species dependent water content – chlorophyll fluorescence thresholds and relations between volumetric water content and WT depth, Monte Carlo simulations suggest that moss had nearly twice the likelihood of being stressed at shallow sites (0.38 ± 0.24) compared to deep sites (0.22 ± 0.18). This study provides evidence that mosses in shallow peatland may be particularly vulnerable to warmer and drier climates in the future, but where species composition may play an important role. We argue that a critical ‘threshold’ peat depth specific for different hydrogeological and hydroclimatic regions can be used to assess what peatlands are especially vulnerable to climate change mediated drought.     

Changes in bark properties and hydrology following prescribed fire in Pinus taeda and Quercus montana

In the eastern United States, the use of prescribed fire as a silvicultural technique to manage for desirable upland tree species is increasing in popularity. Bark physical properties such as thickness, density, and porosity have known associations with fire tolerance among species. These physical properties simultaneously influence rainfall interception and canopy storage and thus are of interest across a range of disciplines. Furthermore, while these characteristics are innate to a species, it is unknown whether repeated exposure to fire facilitates physical change in bark structure and whether these changes are consistent among species. To answer these questions, bark samples were collected from mature pine (Pinus taeda L.) and oak (Quercus montana Willd.) trees from sites across the Bankhead National Forest in Alabama, USA under three different burn regimes: 3-year cycle, 9-year cycle, and no fire. Samples were analysed in the laboratory for bulk density, porosity, water storage capacity, and hygroscopicity (the amount of atmospheric water vapour absorbed by bark during non-rainfall conditions). Drying rates of saturated samples under simulated wetting conditions were also assessed. Oak bark had higher bulk density, lower porosity, and dried slower than pine bark. Interestingly, bark from both species had lower bulk density, higher porosity, greater water storage capacity, and dried faster in stands that were burned every 3 years compared to other fire regimes (p < 0.001). In summary, this study demonstrates that prescribed fire regimes in an eastern US forest alter bark structure and thus influence individual tree control on hydrological processes. The increase in bark water storage capacity, coupled with faster bark evaporation times may lead to less water inputs to the forest floor and drier overall conditions. Further investigation of this fire-bark-water feedback loop is necessary to understand the extent of these mechanisms controlling landscape-scale conditions.     

Estimation of freshwater and material fluxes from the Yangtze River into the East China Sea by using TOPEX/Poseidon altimeter data

In order to determine material fluxes in rivers by non‐contact methods, it is essential to estimate river discharge first. Although developed and optimized for open oceans, satellite radar altimetry has the potential to monitor variations in the levels of inland waters such as lakes and rivers. Making use of the concept of an ‘assumed reference point’, we converted TOPEX/Poseidon satellite altimetry data on water level variations in the Yangtze River (Changjiang) to ‘water level’ data. We also used ‘water level’ time‐series data and in situ river discharge to establish a rating curve. By use of the rating curve, we converted data on ‘water level’ derived from 7 years (1993–99) of TOPEX/Poseidon data to actual river discharge. On the basis of statistical correlation between discharge and nutrient concentration data collected in 1987–88 and in 1998–99, we estimated the total amounts of freshwater and material fluxes transferred by the Yangtze River during the 1990s. The result reveals that an overall, but very slight, increase in freshwater and material fluxes occurred during the 1990s. Copyright © 2005 John Wiley & Sons, Ltd.     

Tritium in ground water of the georgia piedmont: Implications for recharge and flow paths

Environmental tritium was measured in 33 natural water samples representative of precipitation, stream runoff, and groundwater (derived principally from production wells) within the Georgia Piedmont Province. Major ion analyses were used to assist in the interpretation of the tritium results. Tritium concentrations were significantly greater within shallow groundwater derived from the regolith (28–34 TU) and stream runoff (25–30 TU) than within recent rainfall (4–17 TU). Based upon the decay-corrected tritium input function, this probably indicates that at least some of the shallow water is stored within the regolith for a period of approximately 25 years. A ‘post-bomb’ component of recharge was present in all groundwater derived from production wells in the study area. Groundwater sampled from the bedrock aquifers was commonly less tritiated than either stream runoff or shallow water stored in the regolith. the lower tritium concentrations May, have resulted from the mixing of ‘pre-bomb’ water stored within the fractures or the transitional zone directly above the bedrock and modern water stored in the shallow regolith. the preponderance of modern water provides evidence that groundwater flow paths are areally restricted within this setting, probably confined to local surface water drainage basins. the residence time of groundwater in the Piedmont is limited by the lack of deep, gravity-driven regional flow and the localized vertical flow induced by pumping. the results of this study indicate that relatively small tritium concentration variations (10-20 TU) May, have regional hydrological significance in the southeastern Piedmont Province and similar settings.     

Implications of measurement metrics on soil freezing curves: A simulation of freeze–thaw hysteresis

Soil freeze–thaw events have important implications for water resources, flood risk, land productivity, and climate change. A property of these phenomena is the relationship between unfrozen water content and sub-freezing temperature, known as the soil freezing characteristic curve (SFC). It is documented that this relationship exhibits hysteretic behaviour when frozen soil thaws, leading to the definition of the soil thawing characteristic curve (STC). Although explanations have been given for SFC/STC hysteresis, the effect that ‘scale’ – particularly ‘measurement scale’ – may have on these curves has received little attention. The most commonly used measurement scale metric is the ‘support’, which is the spatial (or temporal) unit within which the measured variable is integrated or soil volume sampled. We show (a) measurement support can influence the range and shape of the SFC and (b) hysteresis can be attributed, in part, to the support and location of the measurements comprising the SFC/STC. We simulated lab measured temperature, volumetric water content (VWC), and permittivity from soil samples undergoing freeze–thaw transitions using Hydrus-1D and a modified Dobson permittivity model. To assess the effect of measurement support and location on SFC/STC, we masked the simulated temperature and VWC/permittivity extent to match the instrument's support and location. By creating a detailed simulation of the intra- and inter-support variability associated with the penetration of a freezing front, we demonstrate how measurement support and location can influence the temperature range over which water freezing events are captured. We show it is possible to simulate hysteresis in homogenous media with purely geometric considerations, suggesting that SFC/STC hysteresis may be more of an apparent phenomenon than mechanistically real. Lastly, we develop an understanding of how the location and support of soil temperature and VWC/permittivity measurements influence the temperature range over which water freezing events are captured.     

Sensitivity of nitrate concentration-discharge patterns to soil nitrate distribution and drainage properties in the vertical dimension

We propose a conceptual model that examines the ‘variable source area’ (VSA) and ‘nitrate flushing’ hypothesis in the vertical direction, and use this approach to explain nitrate concentration–discharge relationships in a semi-arid watershed. We use an eco-hydrology simulation model (RHESSys) to show that small changes in the vertical distribution of nitrate mass and their interaction with soil hydraulic conductivity can result in abrupt changes in the nitrate concentration–discharge relationship. We show that the estimated concentration–discharge relationship is sensitive to the parameters governing soil vertical nitrate distribution and soil hydraulic conductivity, at both patch scale and watershed scale, where lateral redistribution of water and nitrate is also accounted for. Given heterogeneity in nitrogen inputs, uptake processes, soil drainage and storage processes, substantial variation in parameters that describe rate of changes in vertical distribution of soil nitrate and hydraulic properties is likely both within and between watersheds. Thus, we argue that vertical ‘variable source area’ processes may be as important as lateral VSA in determining concentration discharge relationships.     

MID-WINTER STEMFLOW DRAINAGE FROM BIGTOOTH ASPEN (POPULUS GRANDIDENTATA MICHX.) IN CENTRAL MASSACHUSETTS

Under winter conditions, stemflow drainage in forested ecosystems is often assumed to be a negligible component of the hydrological cycle. This paper reports on mid-winter stemflow drainage from the broadleaved deciduous tree species Populus grandidentata. Stemflow volumes from this species at air temperatures of < 0°C were found to be comparable to rainfall-generated stemflow during summer. Over the three-month period January–March 1993, stemflow ranged from 5.4 to 9.9% of the incident gross precipitation. Expressed as depth equivalents per unit trunk basal area, these stemflow inputs ranged from 1.8 to 4.9 m. These concentrated mid-winter inputs of liquid water to the bases of canopy trees were attributable to: (1) snow interception by the leafless woody frame of each tree; (2) snow retention by glazed ice precipitation associated with the snowfall event; (3) increased temperature at the bark/snow interface caused by the low albedo of the bark tissue; and (4) convergence of snowmelt drainage from steeply inclined upthrust primary branches. The hydrological and ecological significance of liquid water inputs to the forest floor under sub-zero conditions are discussed. © 1997 by John Wiley & Sons Ltd.     

The geomorphological unit hydrograph from a historical‐critical perspective

In this paper we present a brief overview of geomorphological instantaneous unit hydrograph (GIUH) theories and analyze their successful path without hiding their limitations. The history of the GIUH is subdivided into three major sections. The first is based on the milestone works of Rodríguez‐Iturbe and Valdés (Water Resources Research 1979; 15 (6): 1409–1420) and Gupta et al. (Water Resources Research 1980; 16 (5): 855–862), which recognized that a treatment of water discharges with ‘travel times’ could provide a rich interpretation of the theory of the instantaneous unit hydrograph (IUH). We show how this was possible, what assumptions were made, which of these assumptions can be relaxed, and which have become obsolete and been discarded. The second section focuses on the width‐function‐based IUH (WFIUH) approach and its achievements in assessing the interplay of the topology and geometry of the network with water dynamics. The limitations of the WFIUH approach are described, and a way to work around them is suggested. Finally, a new formal approach to estimating the water budget by ‘travel times’, which derives from a suitable use of the water budget equation and some hypotheses, has been introduced and disentangled. Copyright © 2015 John Wiley & Sons, Ltd.     

Evaluation of vertical seismic response for a large‐scale beam‐supported aqueduct

By the theories of potential flow and structural vibration, the formulae for evaluating the ‘wet’ (with water) frequencies and mode shapes of the beam‐supported aqueduct are derived through a simplified fluid‐structure interaction analysis. The time‐history formulae of structural responses to the vertical seismic excitation are obtained. Applying the response‐spectrum principle, the equivalent vertical earthquake load exerted on the beam and the corresponding effects are also derived. Several illustrative examples are conducted. The analytical results show that: (i) The ‘wet’ frequencies of the structure are lower than the corresponding ‘dry’ (without water) frequencies due to the participating water mass, but the ‘wet’ mode shapes are identical to the corresponding ‘dry’ ones. (ii) The water mass plays an important role in the vertical seismic response, which varies with the different geological sites. For the different seismic inputs, the deeper the water is, the greater are the structural responses. (iii) The vertical seismic effects on the beam are generally not too small to be neglected and should be considered in the structural designs of a beam‐supported aqueduct. Copyright © 2002 John Wiley & Sons, Ltd.     

The identification of runoff-production mechanisms using environmental isotopes in a tussock grassland catchment,eastern otago,New Zealand

A previous hydrometric study of runoff production in tussock grassland drainage basins in Otago (45°50′S, 169°45′E), New Zealand, revealed a marked change of slope in storm hydrograph recessions. An environmental isotope study was initiated to investigate the runoff mechanisms operating and to test specific hypotheses to explain this break in the hydrograph recession. The results indicated that for quickflow volumes in excess of 10mm, the first part of the storm hydrograph can be attributed to two separate sources, namely, ‘old’ water from a shallow, unconfined groundwater reservoir and ‘new’ water from saturation overland flow on the lower wetlands of concave slopes. Despite the extensive area of wetlands, ‘old’ water runoff from the unconfined groundwater reservoir is delivered more rapidly to the stream than ‘new’ water from saturation overland flow. Substantial surface storage in the wetlands has first to be exceeded before rain becomes a significant part of stream discharge. For quickflow volumes less than 10mm, only ‘old’ water from groundwater contributes to the first part of the hydrograph recession. This means that only the largest 7 per cent of storms (in terms of quickflow volume) generate quickflow containing significant amounts of ‘new water’. The second part of the recession of the storm hydrograph consists of ‘old’ water derived from a remarkably well-mixed shallow unconfined groundwater body.     

An isotopic and geochemical study of snowmelt runoff in a small arctic watershed

‘Old’ water contributions to snowmelt runoff in a stream can be defined as water which was stored in the catchment prior to the start of the runoff event in question. We used mass balance techniques for natural oxygen-18 and several chemical parameters (electrical conductivity, Ca and Mg) to estimate the magnitude and timing of ‘old’ water contributions to snowmelt runoff in the Apex River watershed during the 1983 field season. The Apex River catchment is located in the southern part of Baffin Island, Northwest Territories, Canada. The watershed has an area of 60 km², it is in the zone of continuous permafrost, and the geology is dominantly Precambrian gneiss with sporadic, thin, glacial overburden. The isotopic data indicate that for the snowmelt season of 1983 (third highest peak discharge of 11 years of record), approximately 50 per cent of the peak stream discharge consisted of ‘old’ water. Our data also suggest that about 60 per cent of the entire 1983 hydrograph was ‘old’ water. The chemical parameters give old water contributions which are at least 10 per cent less than the isotopically-derived estimate, but they are consistent with the isotopic estimate during peak flow.     

Modelling stage–discharge relationships in anastomosed bedrock-influenced sections of the Sabie River system

Flow dynamics in a bedrock-influenced river system, the Sabie River, South Africa, have been found to be significantly different from those in temperate alluvial systems. The lack of lateral water connectivity leads to multiple bedrock distributaries with varying water surface elevations across a cross-section. Distributary activation is dependent on upstream breaching of bedrock barriers between distributaries by rising discharge. Where measurement of individual stage–discharge relationships in each distributary was not possible, a ‘Multiple Stage’ model was developed to predict hydraulic conditions in each distributary, using a single measured rating curve and knowledge of individual distributary water surface elevations at a low flow. Use of the ‘Multiple Stage’ model has enabled realistic prediction of channel geometry and hydraulic variables, that accounts for the different stages found in bedrock-influenced sections, yet is not prohibitively data intensive. Predicted ‘Multiple Stage’ results for maximum depth and velocity demonstrate the vast improvement on modelling flow dynamics, when compared to the conventional assumption of a single stage representing the whole cross-section. © 1998 John Wiley & Sons, Ltd.     

Submersible surveys of benthos near a turbidity cloud

Seabed impacted by settling solids from a turbidity plume can be surveyed by submersible using a ‘bounce’ technique to determine depth of the turbidity front, and hence depth above which benthos observations are practical. The protocol adopted was that after frontal depth was determined by a preliminary dive, the site for the benthic survey was selected. On arrival on the seabed a reconnaisance was made for the observers to agree on the identification of visible benthos. Then 100 × 1 m² transects were surveyed using a 1 m² quadrat attached to the submersible in view of the port-side observer. At the site investigated, two shallow water stations (depths 13–15 m, and 16–18 m) showing light deposits of mine tailings were seen to support associations of large epifauna and infauna. The first station had an almost single-species stand of the burrowing anenome Pachycerianthus fimbriatus at approximately 100 · 100 m^?2, and the second a diverse species association with several starfish and other species at 1–2 · 100 m^?2. Burrow holes (up to 3–5 cm diameter) of various sizes indicated that large infaunal species were present. At a deeper station (36–27 m) with heavy tailings and almost no burrow holes, a single Dungeness crab, Cancer magister, and an unidentified shrimp were seen in the 100 m² transect. Each dive site was inshore (shallower) from a routinely monitored benthos station shown repeatedly to support a population of small infauna in the tailings. Burrow holes have the potential for determining the identity and abundance of large infaunal species present if an identification system can be developed. The technique of epifaunal and burrow hole surveys (by submersible and scuba diver), combined with infaunal and sediment core surveys and contemporary theory on infaunal succession, provides the potential for a procedure to map the distribution in shallow water of tailings impact and benthic recovery (succession).     

Integrated cave drip monitoring for epikarst recharge estimation in a dry Mediterranean area,Sif Cave,Israel

Understanding recharge mechanisms and controls in karst regions is extremely important for managing water resources because of the dynamic nature of the system. The objective of this study was to evaluate water percolation through epikarst by monitoring water flow into a cave and conducting artificial irrigation and tracer experiments, at Sif Cave in Wadi Sussi, Israel from 2005 through 2007. The research is based on continuous high‐resolution direct measurements of both rainfall and water percolation in the cave chamber collected by three large PVC sheets which integrate drips from three different areas (17, 46, and 52 m²). Barrels equipped with pressure transducers record drip rate and volume for each of the three areas. The combined measured rainfall and cave data enables estimation of recharge into the epikarst and to better understand the relationship of rainfall‐recharge. Three distinct types of flow regimes were identified: (1) ‘Quick flow’ through preferential flow paths (large fractures and conduits); (2) ‘Intermediate flow’ through a secondary crack system; and (3) ‘Slow flow’ through the matrix. A threshold of ～100 mm of rain at the beginning of the rainy season is required to increase soil water content allowing later rainfall events to percolate deeper through the soil and to initiate dripping in the cave. During winter, as the soil water content rises, the lag time between a rain event and cave drip response decreases. Annual recharge (140–160 mm in different areas in the cave) measured represents 30–35% of annual rainfall (460 mm). Copyright © 2011 John Wiley & Sons, Ltd.     

Development of macrofaunal communities on dredged material used for mudflat enhancement: a comparison of three beneficial use schemes after one year

In recent years, dredged material has become regarded as a potential resource and used to create and/or improve intertidal habitats ('beneficial use' schemes). This paper presents the results of a sampling programme to investigate the short-term macrofaunal recovery of three beneficial use schemes in south-east England in terms of species and functional diversity. Environmental parameters (sediment redox potential, and water, organic carbon and silt/clay contents) and univariate community attributes (total individuals and species, diversity, evenness and biomass) at the recharge sites had attained reference levels at two schemes while assemblages differed significantly in terms of species composition at all three schemes. While trophic group proportionality had re-established at one scheme, an increased grazer dominance was apparent at another while the proportion of sub-surface deposit feeders decreased at the third. Total individuals and species number of the developing communities were negatively correlated with sediment redox potential at 4 cm and % silt/clay, respectively. The implications of these results for monitoring the recovery of future fine-grained beneficial use schemes are discussed.     

Evidence of field-scale shifts in transpiration dynamics following bark beetle infestation: Stomatal conductance responses

Amplified eruptive outbreaks of bark beetles as a consequence of climate change can cause tree mortality that significantly affects terrestrial water and carbon fluxes. However, the lack of field-scale observations of underlying physiological mechanisms currently hampers the expression of such ecosystem disturbances in predictive modelling. Based on a unique flux tower dataset from a subalpine forest located in the Rocky Mountains, mechanisms of stomatal response to an extensive bark beetle outbreak were investigated using various models and parametrizations. The datasets cover a decade, including the periods of pre-infestation, infestation, and post-infestation. Field measurements showed considerable decreases in evapotranspiration (ET), transpiration (T), and leaf area index (LAI) during the two-year infestation period compared to the pre-infestation period. Model interpretations of observed water and carbon fluxes indicated that the overall reductions in T were not solely due to decreased LAI, but also to changes in physiological behaviours. The summer season's canopy-scale stomatal conductance was significantly reduced during the infestation period, from 0.0018 to 0.0011 m s⁻¹. One primary reason for the observed variations is likely that the bark beetle infestation hampers the water transport in the xylem. The damage of xylem has important implications for water use efficiency (WUE), which also significantly influences the parameterization of stomatal conductance. When using stomatal conductance models to forecast ecosystem dynamics, it is crucial to recalibrate the model's parameters to ensure the accurate depiction of stomatal dynamics during various infestation periods. The neglect of the temporal variability of canopy-scale stomatal conductance under ecosystem disturbances (e.g., bark beetle infestations) in current earth system models, therefore, requires specific attention in assessments of large-scale water and carbon balances.     

Seasonality intensification and long-term winter cooling as a part of the Late Pliocene climate development

A mutual climatic range method is applied to the Mediterranean marine pollen record of Semaforo (Vrica section, Calabria, Italy) covering the period from ∼2.46 Ma to ∼2.11 Ma. The method yields detailed information on summer, annual and winter temperatures and on precipitation during the nine obliquity and precession-controlled ‘glacial’ periods (marine isotope stages 96 to 80) and eight ‘interglacial’ periods (marine isotope stages 95 to 81) characterising this time interval. The reconstruction reveals higher temperatures of at least 2.8 °C in mean annual and 2.2 °C in winter temperatures, and 500 mm in precipitation during the ‘interglacials’ as compared to the present-day climate in the study area. During the ‘glacials’, temperatures are generally lower as compared to the present-day climate in the region, but precipitation is equivalent. Along the consecutive ‘interglacials’, a trend toward a reduction in annual and winter temperatures by more than 2.3 °C, and toward a higher seasonality is observed. Along the consecutive ‘glacials’, a trend toward a strong reduction in all temperature parameters of at least 1.6 °C is reconstructed. Climatic amplitudes of ‘interglacial–glacial’ transitions increase from the older to the younger cycles for summer and annual temperatures. The cross-spectral analyses suggest obliquity related warm/humid–cold/dry ‘interglacial–glacial’ cycles which are superimposed by precession related warm/dry– cold/humid cycles. A time displacement in the development of temperatures and precipitation is indicated for the obliquity band by temperatures generally leading precipitation change at ∼4 kyr, and on the precession band of ∼9.6 kyr in maximum.