Refine
Year of publication
- 2018 (314) (remove)
Document Type
- Article (237)
- Doctoral Thesis (25)
- Postprint (20)
- Other (17)
- Review (10)
- Habilitation Thesis (2)
- Conference Proceeding (1)
- Master's Thesis (1)
- Report (1)
Keywords
- climate change (8)
- Himalaya (7)
- permafrost (7)
- Germany (5)
- landslides (4)
- Climate change (3)
- Kalahari (3)
- Naturgefahren (3)
- Palynology (3)
- Permafrost (3)
Institute
- Institut für Geowissenschaften (314) (remove)
TDS-1 GNSS Reflectometry
(2018)
This study presents the development and a systematic evaluation study of GNSS reflectometry wind speeds. After establishing a wind speed retrieval algorithm, UK TechDemoSat-1 (TDS-1) derived winds, from May 2015 to July 2017, are compared to the Advanced Scatterometer (ASCAT). ERA-Interim wind fields of the European Centre for Medium-range Weather Forecasts (ECMWF) and in situ observation from Tropical Atmosphere Ocean buoy array in the Pacific are taken as reference. One-year averaged TDS-1 global winds demonstrate small differences with ECMWF in a majority of areas as well as discuss under- and overestimations. The pioneering TDS-1 winds demonstrate a root-mean-squared error (RMSE) and bias of 2.77 and -0.33 m/s, which are comparable to the RMSE and bias derived by ASCAT winds, as large as 2.31 and 0.25 m/s, respectively. Using buoys measurements as reference, RMSE and bias of 2.23 and -0.03 m/s for TDS-1 as well as 1.40 and -0.68 m/s for ASCAT are obtained. Utilizing rain microwave-infrared estimates of the Tropical Rainfall Measuring Mission, rain-affected observation of both ASCAT and TDS-1 are collected and evaluated. Although ASCAT winds show a significant performance degradation resulting in an RMSE and bias of 3.16 and 1.03 m/s, respectively, during rain condition, TDS-1 shows a more reliable performance with an RMSE and bias of 2.94 and -0.21 m/s, respectively, which indicates the promising capability of GNSS forward scattering for wind retrievals during rain. A decrease in TDS-1-derived bistatic radar cross sections during rain events, at weak winds, is also demonstrated.
Dumortierite was synthesized in piston-cylinder experiments at 2.5-4.0 GPa, 650-700 degrees C in the Al2O3 -B2O3-SiO2-H2O (ABSH) system. Electron-microprobe (EMP) analyses reveal significant boron-excess (up to 0.26 B-[4] per formula unit, pfu) and silicon-deficiency relative to the ideal anhydrous dumortierite stoichiometry Al7BSi3O18 . The EMP data in conjunction with results from single-crystal Raman spectroscopy and powder X-ray diffraction provide evidence that silicon at the tetrahedral site is replaced by excess boron via the substitution Si-[4] <--> B-[4] + H. The Raman spectrum of synthetic dumortierite in the frequency region 2000 4000 cm(-1) comprises eight bands, of which six are located at frequencies below 3400 cm(-1). This points to strong hydrogen bonding, most likely O2-H center dot center dot center dot O7 and O7-H center dot center dot center dot O2, arising from a high number of octahedral vacancies at the All site and substitution of trivalent Al3+ and B3+ for Si4+ at Si1 and Si2 sites, causing decreasing acceptor-donor distances and lower incident valence at the acceptor oxygen. Contrary to the synthetic high-pressure ABSH-dumortierite, magnesiodumortierite from the Dora-Maira Massif, which is assumed to have formed at similar conditions (2.5-3.0 GPa, 700 degrees C), does not show any B-excess. Tourmaline shows an analogous behaviour in that magnesium-rich (e.g., dravitic) tourmaline formed at high pressure shows no or only minor amounts of tetrahedral boron, whereas natural aluminum-rich tourmaline and synthetic olenitic tourmaline formed at high pressures can incorporate significant amounts of tetrahedral boron. Two mechanisms might account for this discrepancy: (i) Structural avoidance of Mg-[6]-(OR3+)-R-[4] configurations in magnesiodumortierite due to charge deficieny at the oxygens O2 and O7 and strong local distortion of M1 due to decreased O2-O7 bond length, and/or (ii) decreasing fluid mobility of boron in Al-rich systems at high pressures.
Unravelling the spatiotemporal evolution of the Cenozoic Andean (Altiplano-Puna) plateau has been one of the most intriguing problems of South American geology. Despite a number of investigations, the early deformation and uplift history of this area remained largely enigmatic. This paper analyses the Paleogene tectono-sedimentary history of the Casa Grande Basin, in the present-day transition zone between the northern sector of the Puna Plateau and the northern part of the Argentine Eastern Cordillera. Our detailed mapping of synsedimentary structures records the onset of regional contractional deformation during the middle Eocene, revealing reactivation of Cretaceous extensional structures and the development of doubly vergent thrusts. This is in agreement with records from other southern parts of the Puna Plateau and the Eastern Cordillera. These observations indicate the existence of an Eocene broken foreland setting within the region, characterized by low-lying compressional basins and ranges with spatially disparate sectors of deformation, which was subsequently subjected to regional uplift resulting in the attainment of present-day elevations during the Neogene.
Coastal ecosystems in the Arctic are affected by climate change. As summer rainfall frequency and intensity are projected to increase in the future, more organic matter, nutrients and sediment could bemobilized and transported into the coastal nearshore zones. However, knowledge of current processes and future changes is limited. We investigated streamflow dynamics and the impacts of summer rainfall on lateral fluxes in a small coastal catchment on Herschel Island in the western Canadian Arctic. For the summer monitoring periods of 2014-2016, mean dissolved organic matter flux over 17 days amounted to 82.7 +/- 30.7 kg km(-2) and mean total dissolved solids flux to 5252 +/- 1224 kg km(-2). Flux of suspended sediment was 7245 kg km(-2) in 2015, and 369 kg km(-2) in 2016. We found that 2.0% of suspended sediment was composed of particulate organic carbon. Data and hysteresis analysis suggest a limited supply of sediments; their interannual variability is most likely caused by short-lived localized disturbances. In contrast, our results imply that dissolved organic carbon is widely available throughout the catchment and exhibits positive linear relationship with runoff. We hypothesize that increased projected rainfall in the future will result in a similar increase of dissolved organic carbon fluxes.
Nearly 13,000 years ago, the warming trend into the Holocene was sharply interrupted by a reversal to near glacial conditions. Climatic causes and ecological consequences of the Younger Dryas (YD) have been extensively studied, however proxy archives from the Mediterranean basin capturing this period are scarce and do not provide annual resolution. Here, we report a hydroclimatic reconstruction from stable isotopes (delta O-18, delta C-13) in subfossil pines from southern France. Growing before and during the transition period into the YD (12 900-12 600 cal BP), the trees provide an annually resolved, continuous sequence of atmospheric change. Isotopic signature of tree sourcewater (delta O-18(sw)) and estimates of relative air humidity were reconstructed as a proxy for variations in air mass origin and precipitation regime. We find a distinct increase in inter-annual variability of sourcewater isotopes (delta O-18(sw)), with three major downturn phases of increasing magnitude beginning at 12 740 cal BP. The observed variation most likely results from an amplified intensity of North Atlantic (low delta O-18(sw)) versus Mediterranean (high delta O-18(sw)) precipitation. This marked pattern of climate variability is not seen in records from higher latitudes and is likely a consequence of atmospheric circulation oscillations at the margin of the southward moving polar front.
Nearly 13,000 years ago, the warming trend into the Holocene was sharply interrupted by a reversal to near glacial conditions. Climatic causes and ecological consequences of the Younger Dryas (YD) have been extensively studied, however proxy archives from the Mediterranean basin capturing this period are scarce and do not provide annual resolution. Here, we report a hydroclimatic reconstruction from stable isotopes (delta O-18, delta C-13) in subfossil pines from southern France. Growing before and during the transition period into the YD (12 900-12 600 cal BP), the trees provide an annually resolved, continuous sequence of atmospheric change. Isotopic signature of tree sourcewater (delta O-18(sw)) and estimates of relative air humidity were reconstructed as a proxy for variations in air mass origin and precipitation regime. We find a distinct increase in inter-annual variability of sourcewater isotopes (delta O-18(sw)), with three major downturn phases of increasing magnitude beginning at 12 740 cal BP. The observed variation most likely results from an amplified intensity of North Atlantic (low delta O-18(sw)) versus Mediterranean (high delta O-18(sw)) precipitation. This marked pattern of climate variability is not seen in records from higher latitudes and is likely a consequence of atmospheric circulation oscillations at the margin of the southward moving polar front.
Pre-existing weakness zones in the lithosphere such as transform faults/fracture zones and extinct mid-oceanic ridges have been suggested to facilitate subduction initiation in an intra-oceanic environment. Here, we propose that the additional forcing coming from the mantle suction flow is required to trigger the conversion of a fracture zone/transform fault into a converging plate boundary. This suction flow can be induced either from the slab remnants of former converging plate boundaries or/and from slabs of neighbouring active subduction zones. Using 2-D coupled thermo-mechanical models, we show that a sufficiently strong mantle flow is able to convert a fracture zone/transform fault into a subduction zone. However, this process is feasible only if the fracture zone/transform fault is very close to the mid-oceanic ridge. Our numerical model results indicate that time of subduction initiation depends on the velocity, domain size and location of mantle suction flow and age of the oceanic plate.
Sub-seasonal thaw slump mass wasting is not consistently energy limited at the landscape scale
(2018)
Predicting future thaw slump activity requires a sound understanding of the atmospheric drivers and geomorphic controls on mass wasting across a range of timescales. On sub-seasonal timescales, sparse measurements indicate that mass wasting at active slumps is often limited by the energy available for melting ground ice, but other factors such as rainfall or the formation of an insulating veneer may also be relevant. To study the sub-seasonal drivers, we derive topographic changes from single-pass radar interferometric data acquired by the TanDEM-X satellites. The estimated elevation changes at 12m resolution complement the commonly observed planimetric retreat rates by providing information on volume losses. Their high vertical precision (around 30 cm), frequent observations (11 days) and large coverage (5000 km(2)) allow us to track mass wasting as drivers such as the available energy change during the summer of 2015 in two study regions. We find that thaw slumps in the Tuktoyaktuk coastlands, Canada, are not energy limited in June, as they undergo limited mass wasting (height loss of around 0 cm day 1) despite the ample available energy, suggesting the widespread presence of early season insulating snow or debris veneer. Later in summer, height losses generally increase (around 3 cm day 1), but they do so in distinct ways. For many slumps, mass wasting tracks the available energy, a temporal pattern that is also observed at coastal yedoma cliffs on the Bykovsky Peninsula, Russia. However, the other two common temporal trajectories are asynchronous with the available energy, as they track strong precipitation events or show a sudden speed-up in late August respectively. The observed temporal patterns are poorly related to slump characteristics like the headwall height. The contrasting temporal behaviour of nearby thaw slumps highlights the importance of complex local and temporally varying controls on mass wasting.
The lack of process-based classification procedures may lead to unrealistic hyetograph design due to complex oscillation of rainfall depths when assimilated at high temporal resolutions. Four consecutive years of sub-hourly rainfall data were assimilated in three study areas (Guaraira, GEB, Sao Joao do Cariri, CEB, and Aiuaba, AEB) under distinct climates (very hot semi-arid and tropical wet). This study aimed to define rainfall events (for Minimum Inter-event Time, MIT, and Minimum Rainfall Depth, MRD, equal to 30 min and 1.016 mm, respectively), classify their hyetograph types (rectangular, R, unimodal with left-skewed, UL, right-skewed, UR, and centred peaks, UC, bimodal, B, and shapeless, SL), and compare their key rainfall properties (frequency, duration, depth, rate and peak). A rain pulse aggregation process allowed for reshaping SL-events for six different time spans varying from 2 to 30 min. The results revealed that the coastal area held predominantly R-events (64% events and 49% rainfall depth), in western semi-arid prevailed UL-events (57% events and 63% rainfall depth), whereas in eastern semi-arid mostly were R-events (61% events and 30% rainfall depth) similar to coastal area. It is concluded that each cloud formation type had important effects on hyetograph properties, differentiating them even within the same climate.
Active and passive source data from two seismic experiments within the interdisciplinary project TIPTEQ (from The Incoming Plate to mega Thrust EarthQuake processes) were used to image and identify the structural and petrophysical properties (such as P- and S-velocities, Poisson's ratios, pore pressure, density and amount of fluids) within the Chilean seismogenic coupling zone at 38.25°S, where in 1960 the largest earthquake ever recorded (Mw 9.5) occurred. Two S-wave velocity models calculated using traveltime and noise tomography techniques were merged with an existing velocity model to obtain a 2D S-wave velocity model, which gathered the advantages of each individual model. In a following step, P- and S-reflectivity images of the subduction zone were obtained using different pre stack and post-stack depth migration techniques. Among them, the recent prestack line-drawing depth migration scheme yielded revealing results. Next, synthetic seismograms modelled using the reflectivity method allowed, through their input 1D synthetic P- and S-velocities, to infer the composition and rocks within the subduction zone. Finally, an image of the subduction zone is given, jointly interpreting the results from this work with results from other studies. The Chilean seismogenic coupling zone at 38.25°S shows a continental crust with highly reflective horizontal, as well as (steep) dipping events. Among them, the Lanalhue Fault Zone (LFZ), which is interpreted to be east-dipping, is imaged to very shallow depths. Some steep reflectors are observed for the first time, for example one near the coast, related to high seismicity and another one near the LFZ. Steep shallow reflectivity towards the volcanic arc could be related to a steep west-dipping reflector interpreted as fluids and/or melts, migrating upwards due to material recycling in the continental mantle wedge. The high resolution of the S-velocity model in the first kilometres allowed to identify several sedimentary basins, characterized by very low P- and S-velocities, high Poisson's ratios and possible steep reflectivity. Such high Poisson's ratios are also observed within the oceanic crust, which reaches the seismogenic zone hydrated due to bending-related faulting. It is interpreted to release water until reaching the coast and under the continental mantle wedge. In terms of seismic velocities, the inferred composition and rocks in the continental crust is in agreement with field geology observations at the surface along the proflle. Furthermore, there is no requirement to call on the existence of measurable amounts of present-day fluids above the plate interface in the continental crust of the Coastal Cordillera and the Central Valley in this part of the Chilean convergent margin. A large-scale anisotropy in the continental crust and upper mantle, previously proposed from magnetotelluric studies, is proposed from seismic velocities. However, quantitative studies on this topic in the continental crust of the Chilean seismogenic zone at 38.25°S do not exist to date.
Carbonate-rich silicate and carbonate melts play a crucial role in deep Earth magmatic processes and their melt structure is a key parameter, as it controls physical and transport properties. Carbon-rich melts can be strongly enriched in trace elements, but the structural incorporation mechanisms of these elements are difficult to study because such melts generally cannot be quenched to glasses. In this contribution we investigate the influence of CO2 on the local environments of trace elements contained in silicate glasses with variable CO2 concentrations and in silicate and carbonate melts. The melts were studied in-situ at high pressure and temperature conditions using the Paris-Edinburgh press (2.2 to 2.6 GPa and 1200 to 1500 degrees C). The compositions studied include sodium-rich peralkaline silicate melts and glasses and carbonate melts similar to those occurring naturally at Oldoinyo Lengai volcano. The local environments of yttrium (Y), lanthanum (La) and strontium (Sr) were investigated using extended X-ray absorption fine structure (EXAFS) spectroscopy. Main findings of the study suggest: (1) In peralkaline silicate glasses the local structure of Y is unaffected by the CO2 content. Contrary, a slight increase of oxygen bond lengths of Sr and La is inferred with increasing CO2 content in peralkaline glasses, while they remain constant in glasses of even higher peralkalinity independent of the CO2 content. (2) In silicate melts of different CO2 contents Y-O bond lengths are constant, while a slight increase within carbonate melt compositions is deduced. On the other hand, a steady bond lengths increase over the whole compositional range is inferred for La-O and Sr-O. This may well be explained by distinct preferences of these elements for specific local environments. Based on these new data, we suggest potential mechanisms for the structural incorporation of these elements, a key step towards understanding their partitioning behavior in natural magmatic systems.
The impressive number of stream gauges in Chile, combined with a suite of past and recent large earthquakes, makes Chile a unique natural laboratory to study several streams that recorded responses to multiple seismic events. We document changes in discharge in eight streams in Chile following two or more large earthquakes. In all cases, discharge increases. Changes in discharge occur for peak ground velocities greater than about 7-11cm/s. Above that threshold, the magnitude of both the increase in discharge and the total excess water do not increase with increasing peak ground velocities. While these observations are consistent with previous work in California, they conflict with lab experiments that show that the magnitude of permeability changes increases with increasing amplitude of ground motion. Instead, our study suggests that streamflow responses are binary. Plain Language Summary Earthquakes deform and shake the surface and the ground below. These changes may affect groundwater flows by increasing the permeability along newly formed cracks and/or clearing clogged pores. As a result, groundwater flow may substantially increase after earthquakes and remain elevated for several months. Here we document streamflow anomalies following multiple high magnitude earthquakes in multiple streams in one of the most earthquake prone regions worldwide, Chile. We take advantage of the dense monitoring network in Chile that recorded streamflow since the 1940s. We show that once a critical ground motion is exceeded, streamflow responses to earthquakes can be expected.
Germination marks a critical transition in plant life that is prone to high mortality. Strong selection pressure is therefore expected to finely tune it to environmental conditions. Our study on the common Mediterranean grass Brachypodium hybridum assessed whether germination behavior changes systematically along a steep natural rainfall gradient ranging from harsh desert to rather mild mesic-Mediterranean conditions, We specifically tested hypotheses that germination behavior confers greater risk-spreading in populations from drier, unpredictable environments, and that seeds from wetter populations are better competitors. In 14 populations (spanning 114-954 mm annual rainfall) we assessed three alternative key parameters of germination in a greenhouse experiment: between-year dormancy, days to emergence within a season, and temporal spread. Addition of neighbor seeds accounted for competition as another crucial environmental factor. In six of the 14 populations, we also compared seeds originating from corresponding north (more mesic) and south (more arid) exposed hill slopes to test whether germination patterns along the large-scale rainfall gradient are paralleled at this smaller scale. B. hybridum exhibited generally high germination fractions and rapid emergence with very little temporal spread, indicating overall little risk-spreading germination. Surprisingly, none of the three parameters changed systematically with increasing aridity, neither at large scale along the rainfall gradient nor at small scale between north and south exposures. Neighbor seeds, however, mildly suppressed germination. Germination of neighbor seeds, in turn, was more strongly suppressed by B. hybridum seeds from drier populations, and this effect was stronger for forb than for grass neighbor species. Our results provide strong evidence that increased risk-spreading germination is not a universal, essential strategy to persist in increasingly dry, unpredictable environments. They also highlight that competition with neighbors occurs even at the earliest plant life stage. Since neighbor effects were species-specific, competition among seeds can affect community composition at later plant stages.
To evaluate the spatiotemporal variations of ground motions in northern Chile, we built a high-quality rock seismic acceleration database and an interface earthquakes catalog. Two ground-motion prediction equation (GMPE) models for subduction zones have been tested and validated for the area. They were then used as backbone models to describe the time-space variations of earthquake frequency content (Fourier and response spectra). Consistent with previous studies of large subduction earthquakes, moderate interface earthquakes in northern Chile show an increase of the high-frequency energy released with depth. A regional variability of earthquake frequency content is also observed, which may be related to a lateral segmentation of the mechanical properties of the subduction interface. Finally, interface earthquakes show a temporal evolution of their frequency content in the earthquake sequence associated with the 2014 Iquique M-w 8.1 megathrust earthquake. Surprisingly, the change does not occur with the mainshock but is associated with an 8 month slow slip preceding the megathrust. Electronic Supplement: Strong-motion database.
The potential link between erosion rates at the Earth’s surface and changes in global climate has intrigued geoscientists for decades1,2 because such a coupling has implications for the influence of silicate weathering3,4 and organic-carbon burial5 on climate and for the role of Quaternary glaciations in landscape evolution1,6. A global increase in late-Cenozoic erosion rates in response to a cooling, more variable climate has been proposed on the basis of worldwide sedimentation rates7. Other studies have indicated, however, that global erosion rates may have remained steady, suggesting that the reported increases in sediment-accumulation rates are due to preservation biases, depositional hiatuses and varying measurement intervals8,9,10. More recently, a global compilation of thermochronology data has been used to infer a nearly twofold increase in the erosion rate in mountainous landscapes over late-Cenozoic times6. It has been contended that this result is free of the biases that affect sedimentary records11, although others have argued that it contains biases related to how thermochronological data are averaged12 and to erosion hiatuses in glaciated landscapes13. Here we investigate the 30 locations with reported accelerated erosion during the late Cenozoic6. Our analysis shows that in 23 of these locations, the reported increases are a result of a spatial correlation bias—that is, combining data with disparate exhumation histories, thereby converting spatial erosion-rate variations into temporal increases. In four locations, the increases can be explained by changes in tectonic boundary conditions. In three cases, climatically induced accelerations are recorded, driven by localized glacial valley incision. Our findings suggest that thermochronology data currently have insufficient resolution to assess whether late-Cenozoic climate change affected erosion rates on a global scale. We suggest that a synthesis of local findings that include location-specific information may help to further investigate drivers of global erosion rates.
Despite its societal relevance, the question whether fluctuations in flood occurrence or magnitude are coherent in space has hardly been addressed in quantitative terms. We investigate this question for Germany by analysing fluctuations in annual maximum series (AMS) values at 68 discharge gauges for the common time period 1932-2005. We find remarkable spatial coherence across Germany given its different flood regimes. For example, there is a tendency that flood-rich/-poor years in sub-catchments of the Rhine basin, which are dominated by winter floods, coincide with flood-rich/-poor years in the southern sub-catchments of the Danube basin, which have their dominant flood season in summer. Our findings indicate that coherence is caused rather by persistence in catchment wetness than by persistent periods of higher/lower event precipitation. Further, we propose to differentiate between event-type and non-event-type coherence. There are quite a number of hydrological years with considerable nonevent-type coherence, i.e. AMS values of the 68 gauges are spread out through the year but in the same magnitude range. Years with extreme flooding tend to be of event-type and non-coherent, i.e. there is at least one precipitation event that affects many catchments to various degree. Although spatial coherence is a remarkable phenomenon, and large-scale flooding across Germany can lead to severe situations, extreme magnitudes across the whole country within one event or within one year were not observed in the investigated period. (C) 2018 Elsevier B.V. All rights reserved.
Some like it hot
(2018)
Accumulating evidence has demonstrated considerable impact of climate change on biodiversity, with terrestrial ectotherms being particularly vulnerable. While climate-induced range shifts are often addressed in the literature, little is known about the underlying ecological responses at individual and population levels. Using a 30-yr monitoring study of the long-living nocturnal gecko Gehyra variegata in arid Australia, we determined the relative contribution of climatic factors acting locally (temperature, rainfall) or distantly (La Nina induced flooding) on ecological processes ranging from traits at the individual level (body condition, body growth) to the demography at population level (survival, sexual maturity, population sizes). We also investigated whether thermoregulatory activity during both active (night) and resting (daytime) periods of the day can explain these responses. Gehyra variegata responded to local and distant climatic effects. Both high temperatures and high water availability enhanced individual and demographic parameters. Moreover, the impact of water availability was scale independent as local rainfall and La Nina induced flooding compensated each other. When water availability was low, however, extremely high temperatures delayed body growth and sexual maturity while survival of individuals and population sizes remained stable. This suggests a trade-off with traits at the individual level that may potentially buffer the consequences of adverse climatic conditions at the population level. Moreover, hot temperatures did not impact nocturnal nor diurnal behavior. Instead, only cool temperatures induced diurnal thermoregulatory behavior with individuals moving to exposed hollow branches and even outside tree hollows for sun-basking during the day. Since diurnal behavioral thermoregulation likely induced costs on fitness, this could decrease performance at both individual and population level under cool temperatures. Our findings show that water availability rather than high temperature is the limiting factor in our focal population of G.variegata. In contrast to previous studies, we stress that drier rather than warmer conditions are expected to be detrimental for nocturnal desert reptiles. Identifying the actual limiting climatic factors at different scales and their functional interactions at different ecological levels is critical to be able to predict reliably future population dynamics and support conservation planning in arid ecosystems.
Signals stored in sediment
(2018)
Tectonic and climatic boundary conditions determine the amount and the characteristics (size distribution and composition) of sediment that is generated and exported from mountain regions. On millennial timescales, rivers adjust their morphology such that the incoming sediment (Qs,in) can be transported downstream by the available water discharge (Qw). Changes in climatic and tectonic boundary conditions thus trigger an adjustment of the downstream river morphology. Understanding the sensitivity of river morphology to perturbations in boundary conditions is therefore of major importance, for example, for flood assessments, infrastructure and habitats. Although we have a general understanding of how rivers evolve over longer timescales, the prediction of channel response to changes in boundary conditions on a more local scale and over shorter timescales remains a major challenge. To better predict morphological channel evolution, we need to test (i) how channels respond to perturbations in boundary conditions and (ii) how signals reflecting the persisting conditions are preserved in sediment characteristics. This information can then be applied to reconstruct how local river systems have evolved over time.
In this thesis, I address those questions by combining targeted field data collection in the Quebrada del Toro (Southern Central Andes of NW Argentina) with cosmogenic nuclide analysis and remote sensing data. In particular, I (1) investigate how information on hillslope processes is preserved in the 10Be concentration (geochemical composition) of fluvial sediments and how those signals are altered during downstream transport. I complement the field-based approach with physical experiments in the laboratory, in which I (2) explore how changes in sediment supply (Qs,in) or water discharge (Qw) generate distinct signals in the amount of sediment discharge at the basin outlet (Qs,out). With the same set of experiments, I (3) study the adjustments of alluvial channel morphology to changes in Qw and Qs,in, with a particular focus in fill-terrace formation. I transfer the findings from the experiments to the field to (4) reconstruct the evolution of a several-hundred meter thick fluvial fill-terrace sequence in the Quebrada del Toro. I create a detailed terrace chronology and perform reconstructions of paleo-Qs and Qw from the terrace deposits. In the following paragraphs, I summarize my findings on each of these four topics.
First, I sampled detrital sediment at the outlet of tributaries and along the main stem in the Quebrada del Toro, analyzed their 10Be concentration ([10Be]) and compared the data to a detailed hillslope-process inventory. The often observed non-linear increase in catchment-mean denudation rate (inferred from [10Be] in fluvial sediment) with catchment-median slope, which has commonly been explained by an adjustment in landslide-frequency, coincided with a shift in the main type of hillslope processes. In addition, the [10Be] in fluvial sediments varied with grain-size. I defined the normalized sand-gravel-index (NSGI) as the 10Be-concentration difference between sand and gravel fractions divided by their summed concentrations. The NSGI increased with median catchment slope and coincided with a shift in the prevailing hillslope processes active in the catchments, thus making the NSGI a potential proxy for the evolution of hillslope processes over time from sedimentary deposits. However, the NSGI recorded hillslope-processes less well in regions of reduced hillslope-channel connectivity and, in addition, has the potential to be altered during downstream transport due to lateral sediment input, size-selective sediment transport and abrasion.
Second, my physical experiments revealed that sediment discharge at the basin outlet (Qs,out) varied in response to changes in Qs,in or Qw. While changes in Qw caused a distinct signal in Qs,out during the transient adjustment phase of the channel to new boundary conditions, signals related to changes in Qs,in were buffered during the transient phase and likely only become apparent once the channel is adjusted to the new conditions. The temporal buffering is related to the negative feedback between Qs,in and channel-slope adjustments. In addition, I inferred from this result that signals extracted from the geochemical composition of sediments (e.g., [10Be]) are more likely to represent modern-day conditions during times of aggradation, whereas the signal will be temporally buffered due to mixing with older, remobilized sediment during times of channel incision.
Third, the same set of experiments revealed that river incision, channel-width narrowing and terrace cutting were initiated by either an increase in Qw, a decrease in Qs,in or a drop in base level. The lag-time between the external perturbation and the terrace cutting determined (1) how well terrace surfaces preserved the channel profile prior to perturbation and (2) the degree of reworking of terrace-surface material. Short lag-times and well preserved profiles occurred in cases with a rapid onset of incision. Also, lag-times were synchronous along the entire channel after upstream perturbations (Qw, Qs,in), whereas base-level fall triggered an upstream migrating knickzone, such that lag-times increased with distance upstream. Terraces formed after upstream perturbations (Qw, Qs,in) were always steeper when compared to the active channel in new equilibrium conditions. In the base-level fall experiment, the slope of the terrace-surfaces and the modern channel were similar. Hence, slope comparisons between the terrace surface and the modern channel can give insights into the mechanism of terrace formation.
Fourth, my detailed terrace-formation chronology indicated that cut-and-fill episodes in the Quebrada del Toro followed a ~100-kyr cyclicity, with the oldest terraces ~ 500 kyr old. The terraces were formed due to variability in upstream Qw and Qs. Reconstructions of paleo-Qs over the last 500 kyr, which were restricted to times of sediment deposition, indicated only minor (up to four-fold) variations in paleo-denudation rates. Reconstructions of paleo-Qw were limited to the times around the onset of river incision and revealed enhanced discharge from 10 to 85% compared to today. Such increases in Qw are in agreement with other quantitative paleo-hydrological reconstructions from the Eastern Andes, but have the advantage of dating further back in time.
Near-surface geophysical imaging of alluvial fan settings is a challenging task but crucial for understating geological processes in such settings. The alluvial fan of Ghor Al-Haditha at the southeast shore of the Dead Sea is strongly affected by localized subsidence and destructive sinkhole collapses, with a significantly increasing sinkhole formation rate since ca. 1983. A similar increase is observed also on the western shore of the Dead Sea, in correlation with an ongoing decline in the Dead Sea level. Since different structural models of the upper 50 m of the alluvial fan and varying hypothetical sinkhole processes have been suggested for the Ghor Al-Haditha area in the past, this study aimed to clarify the subsurface characteristics responsible for sinkhole development.
For this purpose, high-frequency shear wave reflection vibratory seismic surveys were carried out in the Ghor Al-Haditha area along several crossing and parallel profiles with a total length of 1.8 and 2.1 km in 2013 and 2014, respectively. The sedimentary architecture of the alluvial fan at Ghor Al-Haditha is resolved down to a depth of nearly 200 m at a high resolution and is calibrated with the stratigraphic profiles of two boreholes located inside the survey area.
The most surprising result of the survey is the absence of evidence of a thick (> 2–10 m) compacted salt layer formerly suggested to lie at ca. 35–40 m depth. Instead, seismic reflection amplitudes and velocities image with good continuity a complex interlocking of alluvial fan deposits and lacustrine sediments of the Dead Sea between 0 and 200 m depth. Furthermore, the underground section of areas affected by sinkholes is characterized by highly scattering wave fields and reduced seismic interval velocities. We propose that the Dead Sea mud layers, which comprise distributed inclusions or lenses of evaporitic chloride, sulfate, and carbonate minerals as well as clay silicates, become increasingly exposed to unsaturated water as the sea level declines and are consequently destabilized and mobilized by both dissolution and physical erosion in the subsurface. This new interpretation of the underlying cause of sinkhole development is supported by surface observations in nearby channel systems. Overall, this study shows that shear wave seismic reflection technique is a promising method for enhanced near-surface imaging in such challenging alluvial fan settings.
Near-surface geophysical imaging of alluvial fan settings is a challenging task but crucial for understating geological processes in such settings. The alluvial fan of Ghor Al-Haditha at the southeast shore of the Dead Sea is strongly affected by localized subsidence and destructive sinkhole collapses, with a significantly increasing sinkhole formation rate since ca. 1983. A similar increase is observed also on the western shore of the Dead Sea, in correlation with an ongoing decline in the Dead Sea level. Since different structural models of the upper 50 m of the alluvial fan and varying hypothetical sinkhole processes have been suggested for the Ghor Al-Haditha area in the past, this study aimed to clarify the subsurface characteristics responsible for sinkhole development. For this purpose, high-frequency shear wave reflection vibratory seismic surveys were carried out in the Ghor Al-Haditha area along several crossing and parallel profiles with a total length of 1.8 and 2.1 km in 2013 and 2014, respectively. The sedimentary architecture of the alluvial fan at Ghor Al-Haditha is resolved down to a depth of nearly 200 m at a high resolution and is calibrated with the stratigraphic profiles of two boreholes located inside the survey area. The most surprising result of the survey is the absence of evidence of a thick (>2-10 m) compacted salt layer formerly suggested to lie at ca. 35-40 m depth. Instead, seismic reflection amplitudes and velocities image with good continuity a complex interlocking of alluvial fan deposits and lacustrine sediments of the Dead Sea between 0 and 200 m depth. Furthermore, the underground section of areas affected by sinkholes is characterized by highly scattering wave fields and reduced seismic interval velocities. We propose that the Dead Sea mud layers, which comprise distributed inclusions or lenses of evaporitic chloride, sulfate, and carbonate minerals as well as clay silicates, become increasingly exposed to unsaturated water as the sea level declines and are consequently destabilized and mobilized by both dissolution and physical erosion in the subsurface. This new interpretation of the underlying cause of sinkhole development is supported by surface observations in nearby channel systems. Overall, this study shows that shear wave seismic reflection technique is a promising method for enhanced near-surface imaging in such challenging alluvial fan settings.
Despite remarkable tectonostratigraphic similarities along the Himalayan arc, pronounced topographic and exhumational variability exists in different morphotectonic segments. The processes responsible for this segmentation are debated. Of particular interest is a 30- to 40-km-wide orogen-parallel belt of rapid exhumation that extends from central Nepal to the western Himalaya and its possible linkage to a midcrustal ramp in the basal decollement, and the related growth of Lesser Himalayan duplex structures. Here we present 26 new apatite fission track cooling ages from the Beas-Lahul region, at the transition from the Central to the Western Himalaya (77 degrees-78 degrees E) to investigate segmentation in the Himalayan arc from a thermochronologic perspective. Together with previously published data from this part of the orogen, we document significant lateral changes in exhumation between the Dhauladar Range to the west, the Beas-Lahul region, and the Sutlej area to the east of the study area. In contrast to the Himalayan front farther east, exhumation in the far western sectors is focused at the frontal parts of the mountain range and associated with the hanging wall of the Main Boundary Thrust fault ramp. Our results allow us to spatially correlate the termination of the rapid exhumation belt with a midcrustal ramp to the west. We suggest that a plunging anticline at the northwestern edge of the Larji-Kullu-Rampur window represents the termination of the Central Himalayan segment, which is related to the evolution of the Lesser Himalayan duplex. Key Points
Annual greenhouse gas emissions have increased more than threefold between 1950 and 2014, posing a major threat to the integrity of the entire earth system and subsequently to humankind. Consequently, roadmaps towards low-carbon pathways are urgently needed. Our study contributes to a more detailed understanding of the dynamics of country based emission patterns and uses them to discuss prospective low-carbon pathways for countries. As availability of databases on sectoral emissions substantially increased, we employ machine learning techniques to classify emission features and pathways. By doing so, 18 representative emission patterns are derived. Overall emissions from seven sectors and for 167 countries covering the time span from 1950 to 2014 have been used in the analyses. The following significant trends can be observed: a) increasing per capita emissions due to growing fossil fuel use in many parts of the world, b) a decline in per capita emissions in some countries, and c) a shift in the emission shares, i.e., a reduction of agricultural and land use contributions in certain regions. Using the emission patterns, their dynamics, and best performing countries as role models, we show the possibility for gaining a decent human development without significantly increasing per capita emissions.
Runoff predictions in ungauged arctic basins using conceptual models forced by reanalysis data
(2018)
Due to global warming, the problem of assessing water resources and their vulnerability to climate drivers in the Arctic region has become a focus in the recent years. This study is aimed at investigating three lumped hydrological models to predict daily runoff of large-scale Arctic basins in the case of substantial data scarcity. All models were driven only by meteorological forcing reanalysis dataset without any additional information about landscape, soil, or vegetation cover properties of the studied basins. Model parameter regionalization based on transferring the whole parameter set showed good efficiency for predictions in ungauged basins. We run a blind test of the proposed methodology for ensemble runoff predictions on five sub-basins, for which only monthly observations were available, and obtained promising results for current water resources assessment for a broad domain of ungauged basins in the Russian Arctic.
Agricultural soil landscapes of hummocky ground moraines are characterized by 3D spatial patterns of soil types that result from profile modifications due to the combined effect of water and tillage erosion. We hypothesize that crops reflect such soil landscape patterns by increased or reduced plant and root growth. Root development may depend on the thickness and vertical sequence of soil horizons as well as on the structural development state of these horizons at different landscape positions. The hypotheses were tested using field data of the root density (RD) and the root lengths (RL) of winter wheat using the minirhizotron technique. We compared data from plots at the CarboZALF-D site (NE Germany) that are representing a non-eroded reference soil profile (Albic Luvisol) at a plateau position, a strongly eroded profile at steep slope (Calcaric Regosol), and a depositional profile at the footslope (Anocolluvic Regosol). At each of these plots, three Plexiglas access tubes were installed down to approx. 1.5 m soil depth. Root measurements were carried out during the growing season of winter wheat (September 2014-August 2015) on six dates. The root length density (RLD) and the root biomass density were derived from RD values assuming a mean specific root length of 100 m g(-1). Values of RD and RLD were highest for the Anocolluvic Regosol and lowest for the Calcaric Regosol. The maximum root penetration depth was lower in the Anocolluvic Regosol because of a relatively high and fluctuating water table at this landscape position. Results revealed positive relations between below-ground (root) and above-ground crop parameters (i.e., leaf area index, plant height, biomass, and yield) for the three soil types. Observed root densities and root lengths in soils at the three landscape positions corroborated the hypothesis that the root system was reflecting erosion-induced soil profile modifications. Soil landscape position dependent root growth should be considered when attempting to quantify landscape scale water and element balances as well as agricultural productivity.
Rock glaciers in semiarid mountains contain large amounts of ice and might be important water stores aside from glaciers, lakes, and rivers. Yet whether and how rock glaciers interact with river channels in mountain valleys remains largely unresolved. We examine the potential for rock glaciers to block or disrupt river channels, using a new inventory of more than 2000 intact rock glaciers that we mapped from remotely sensed imagery in the Karakoram (KR), Tien Shan (TS), and Altai (ALT) mountains. We find that between 5% and 14% of the rock glaciers partly buried, blocked, diverted or constricted at least 95 km of mountain rivers in the entire study area. We use a Bayesian robust logistic regression with multiple topographic and climatic inputs to discern those rock glaciers disrupting mountain rivers from those with no obvious impacts. We identify elevation and potential incoming solar radiation (PISR), together with the size of feeder basins, as dominant predictors, so that lower-lying and larger rock glaciers from larger basins are more likely to disrupt river channels. Given that elevation and PISR are key inputs for modelling the regional distribution of mountain permafrost from the positions of rock-glacier toes, we infer that river-blocking rock glaciers may be diagnostic of non-equilibrated permafrost. Principal component analysis adds temperature evenness and wet-season precipitation to the controls that characterise rock glaciers impacting on rivers. Depending on the choice of predictors, the accuracy of our classification is moderate to good with median posterior area-under-the-curve values of 0.71-0.89. Clarifying whether rapidly advancing rock glaciers can physically impound rivers, or fortify existing dams instead, deserves future field investigation. We suspect that rock-glacier dams are conspicuous features that have a polygenetic history and encourage more research on the geomorphic coupling between permafrost lobes, river channels, and the sediment cascades of semiarid mountain belts. (c) 2018 John Wiley & Sons, Ltd.
The formation and breaching of natural dammed lakes have formed the landscapes, especially in seismically active high-mountain regions. Dammed lakes pose both, potential water resources, and hazard in case of dam breaching. Central Asia has mostly arid and semi-arid climates. Rock glaciers already store more water than ice-glaciers in some semi-arid regions of the world, but their distribution and advance mechanisms are still under debate in recent research. Their impact on the water availability in Central Asia will likely increase as temperatures rise and glaciers diminish.
This thesis provides insight to the relative age distribution of selected Kyrgyz and Kazakh rock glaciers and their single lobes derived from lichenometric dating. The size of roughly 8000 different lichen specimens was used to approximate an exposure age of the underlying debris surface. We showed that rock-glacier movement differs signifcantly on small scales. This has several implications for climatic inferences from rock glaciers. First, reactivation of their lobes does not necessarily point to climatic changes, or at least at out-of-equilibrium conditions. Second, the elevations of rock-glacier toes can no longer be considered as general indicators of the limit of sporadic mountain permafrost as they have been used traditionally.
In the mountainous and seismically active region of Central Asia, natural dams, besides rock glaciers, also play a key role in controlling water and sediment infux into river valleys. However, rock glaciers advancing into valleys seem to be capable of infuencing the stream network, to dam rivers, or to impound lakes. This influence has not previously been addressed. We quantitatively explored these controls using a new inventory of 1300 Central Asian rock glaciers. Elevation, potential incoming solar radiation, and the size of rock glaciers and their feeder basins played key roles in predicting dam appearance. Bayesian techniques were used to credibly distinguish between lichen sizes on rock glaciers and their lobes, and to find those parameters of a rock-glacier system that are most credibly expressing the potential to build natural dams.
To place these studies in the region's history of natural dams, a combination of dating of former lake levels and outburst flood modelling addresses the history and possible outburst flood hypotheses of the second largest mountain lake of the world, Issyk Kul in Kyrgyzstan. Megafoods from breached earthen or glacial dams were found to be a likely explanation for some of the lake's highly fluctuating water levels. However, our detailed analysis of candidate lake sediments and outburst-flood deposits also showed that more localised dam breaks to the west of Issyk Kul could have left similar geomorphic and sedimentary evidence in this Central Asian mountain landscape. We thus caution against readily invoking megafloods as the main cause of lake-level drops of Issyk Kul. In summary, this thesis addresses some new pathways for studying rock glaciers and natural dams with several practical implications for studies on mountain permafrost and natural hazards.
The extent of continental rifts and subduction zones through deep geological time provides insights into the mechanisms behind supercontinent cycles and the long term evolution of the mantle. However, previous compilations have stopped short of mapping the locations of rifts and subduction zones continuously since the Neoproterozoic and within a self-consistent plate kinematic framework. Using recently published plate models with continuously closing boundaries for the Neoproterozoic and Phanerozoic, we estimate how rift and peri-continental subduction length vary from 1 Ga to present and test hypotheses pertaining to the supercontinent cycle and supercontinent breakup. We extract measures of continental perimeter-to-area ratio as a proxy for the existence of a supercontinent, where during times of supercontinent existence the perimeter-to-area ratio should be low, and during assembly and dispersal it should be high. The amalgamation of Gondwana is clearly represented by changes in the length of peri-continental subduction and the breakup of Rodinia and Pangea by changes in rift lengths. The assembly of Pangea is not clearly defined using plate boundary lengths, likely because its formation resulted from the collision of only two large continents. Instead the assembly of Gondwana (ca. 520 Ma) marks the most prominent change in arc length and perimeter-to-area ratio during the last billion years suggesting that Gondwana during the Early Palaeozoic could explicitly be considered part of a Phanerozoic supercontinent. Consequently, the traditional understanding of the supercontinent cycle, in terms of supercontinent existence for short periods of time before dispersal and re-accretion, may be inadequate to fully describe the cycle. Instead, either a two-stage supercontinent cycle could be a more appropriate concept, or alternatively the time period of 1 to 0 Ga has to be considered as being dominated by supercontinent existence, with brief periods of dispersal and amalgamation.
In public perception, abnormal animal behavior is widely assumed to be a potential earthquake precursor, in strong contrast to the viewpoint in natural sciences. Proponents of earthquake prediction via animals claim that animals feel and react abnormally to small changes in environmental and physico-chemical parameters related to the earthquake preparation process. In seismology, however, observational evidence for changes of physical parameters before earthquakes is very weak. In this study, we reviewed 180 publications regarding abnormal animal behavior before earthquakes and analyze and discuss them with respect to (1) magnitude-distance relations, (2) foreshock activity, and (3) the quality and length of the published observations. More than 700 records of claimed animal precursors related to 160 earthquakes are reviewed with unusual behavior of more than 130 species. The precursor time ranges from months to seconds prior to the earthquakes, and the distances from a few to hundreds of kilometers. However, only 14 time series were published, whereas all other records are single observations. The time series are often short (the longest is 1 yr), or only small excerpts of the full data set are shown. The probability density of foreshocks and the occurrence of animal precursors are strikingly similar, suggesting that at least parts of the reported animal precursors are in fact related to foreshocks. Another major difficulty for a systematic and statistical analysis is the high diversity of data, which are often only anecdotal and retrospective. The study clearly demonstrates strong weaknesses or even deficits in many of the published reports on possible abnormal animal behavior. To improve the research on precursors, we suggest a scheme of yes and no questions to be assessed to ensure the quality of such claims.
Natural extreme events are an integral part of nature on planet earth. Usually these events are only considered hazardous to humans, in case they are exposed. In this case, however, natural hazards can have devastating impacts on human societies. Especially hydro-meteorological hazards have a high damage potential in form of e.g. riverine and pluvial floods, winter storms, hurricanes and tornadoes, which can occur all over the globe. Along with an increasingly warm climate also an increase in extreme weather which potentially triggers natural hazards can be expected. Yet, not only changing natural systems, but also changing societal systems contribute to an increasing risk associated with these hazards. These can comprise increasing exposure and possibly also increasing vulnerability to the impacts of natural events. Thus, appropriate risk management is required to adapt all parts of society to existing and upcoming risks at various spatial scales. One essential part of risk management is the risk assessment including the estimation of the economic impacts. However, reliable methods for the estimation of economic impacts due to hydro-meteorological hazards are still missing. Therefore, this thesis deals with the question of how the reliability of hazard damage estimates can be improved, represented and propagated across all spatial scales. This question is investigated using the specific example of economic impacts to companies as a result of riverine floods in Germany.
Flood damage models aim to describe the damage processes during a given flood event. In other words they describe the vulnerability of a specific object to a flood. The models can be based on empirical data sets collected after flood events. In this thesis tree-based models trained with survey data are used for the estimation of direct economic flood impacts on the objects. It is found that these machine learning models, in conjunction with increasing sizes of data sets used to derive the models, outperform state-of-the-art damage models. However, despite the performance improvements induced by using multiple variables and more data points, large prediction errors remain at the object level. The occurrence of the high errors was explained by a further investigation using distributions derived from tree-based models. The investigation showed that direct economic impacts to individual objects cannot be modeled by a normal distribution. Yet, most state-of-the-art approaches assume a normal distribution and take mean values as point estimators. Subsequently, the predictions are unlikely values within the distributions resulting in high errors. At larger spatial scales more objects are considered for the damage estimation. This leads to a better fit of the damage estimates to a normal distribution. Consequently, also the performance of the point estimators get better, although large errors can still occur due to the variance of the normal distribution. It is recommended to use distributions instead of point estimates in order to represent the reliability of damage estimates.
In addition current approaches also mostly ignore the uncertainty associated with the characteristics of the hazard and the exposed objects. For a given flood event e.g. the estimation of the water level at a certain building is prone to uncertainties. Current approaches define exposed objects mostly by the use of land use data sets. These data sets often show inconsistencies, which introduce additional uncertainties. Furthermore, state-of-the-art approaches also imply problems of missing consistency when predicting the damage at different spatial scales. This is due to the use of different types of exposure data sets for model derivation and application. In order to face these issues a novel object-based method was developed in this thesis. The method enables a seamless estimation of hydro-meteorological hazard damage across spatial scales including uncertainty quantification. The application and validation of the method resulted in plausible estimations at all spatial scales without overestimating the uncertainty.
Mainly newly available data sets containing individual buildings make the application of the method possible as they allow for the identification of flood affected objects by overlaying the data sets with water masks. However, the identification of affected objects with two different water masks revealed huge differences in the number of identified objects. Thus, more effort is needed for their identification, since the number of objects affected determines the order of magnitude of the economic flood impacts to a large extent.
In general the method represents the uncertainties associated with the three components of risk namely hazard, exposure and vulnerability, in form of probability distributions. The object-based approach enables a consistent propagation of these uncertainties in space. Aside from the propagation of damage estimates and their uncertainties across spatial scales, a propagation between models estimating direct and indirect economic impacts was demonstrated. This enables the inclusion of uncertainties associated with the direct economic impacts within the estimation of the indirect economic impacts. Consequently, the modeling procedure facilitates the representation of the reliability of estimated total economic impacts. The representation of the estimates' reliability prevents reasoning based on a false certainty, which might be attributed to point estimates. Therefore, the developed approach facilitates a meaningful flood risk management and adaptation planning.
The successful post-event application and the representation of the uncertainties qualifies the method also for the use for future risk assessments. Thus, the developed method enables the representation of the assumptions made for the future risk assessments, which is crucial information for future risk management. This is an important step forward, since the representation of reliability associated with all components of risk is currently lacking in all state-of-the-art methods assessing future risk.
In conclusion, the use of object-based methods giving results in the form of distributions instead of point estimations is recommended. The improvement of the model performance by the means of multi-variable models and additional data points is possible, but small. Uncertainties associated with all components of damage estimation should be included and represented within the results. Furthermore, the findings of the thesis suggest that, at larger scales, the influence of the uncertainty associated with the vulnerability is smaller than those associated with the hazard and exposure. This leads to the conclusion that for an increased reliability of flood damage estimations and risk assessments, the improvement and active inclusion of hazard and exposure, including their uncertainties, is needed in addition to the improvements of the models describing the vulnerability of the objects.
Recent studies indicate that the tropical belt has been expanding during recent decades, which can significantly influence precipitation in subtropical climates. Often the location of the tropical border is identified using the Hadley cell edge (HCE) or the subtropical jet stream (STJ), but most studies concentrated on the zonal-mean state, thereby missing regional impacts. Here we detect longitudinal-resolved trends in STJ cores and HCEs over 1979-2016 in both hemispheres at a higher spatial and temporal resolution than previous studies. Besides pronounced regional trend differences in both sign and magnitude, we show that winter HCE and STJ variability increased in the Mediterranean region and decreased over the American and Asian continents. Rainfall variability in these regions changed likewise, and most of those changes can be explained by changes in HCE/STJ variability. This highlights the importance of understanding future tropical belt changes both regionally and in terms of variability. Plain Language Summary We applied a new network-based method to detect motion of the tropical climate border with longitudinal resolution. Depending on the longitudinal position, there are differences in both direction and magnitude of the border motion. In addition, we demonstrate that the rainfall variability is increasing in the Mediterranean region and decreasing over the American and Asian continents, which can be explained by the variability of the tropical belt location.
Both enstatite spectra are very bright in the VIS and NIR and show almost neutral to slightly bluish spectral slopes with a steep absorption in the UV. Very low iron in the enstatites (below similar to 0.04 wt% FeO) already results in weak albeit noticeable absorptions in the VNIR between 0.4 and 0.9 mu m. Orthoenstatite and clinoenstatite are not distinguishable based only on their spectra in the VIS and NIR. At the Reststrahlen bands in the MIR a systematic difference in the number and exact position of local minima at similar to 10 mu m between clinoenstatite and orthoenstatite is evident. This can be used to discern between the polymorphs in this wavelength range. Additionally, we can distinguish between Fe-free low- and high-Ca pyroxenes in the MIR.
Arctic tundra ecosystems have experienced unprecedented change associated with climate warming over recent decades. Across the Pan-Arctic, vegetation productivity and surface greenness have trended positively over the period of satellite observation. However, since 2011 these trends have slowed considerably, showing signs of browning in many regions. It is unclear what factors are driving this change and which regions/landforms will be most sensitive to future browning. Here we provide evidence linking decadal patterns in arctic greening and browning with regional climate change and local permafrost-driven landscape heterogeneity. We analyzed the spatial variability of decadal-scale trends in surface greenness across the Arctic Coastal Plain of northern Alaska (similar to 60,000 km(2)) using the Landsat archive (1999-2014), in combination with novel 30 m classifications of polygonal tundra and regional watersheds, finding landscape heterogeneity and regional climate change to be the most important factors controlling historical greenness trends. Browning was linked to increased temperature and precipitation, with the exception of young landforms (developed following lake drainage), which will likely continue to green. Spatiotemporal model forecasting suggests carbon uptake potential to be reduced in response to warmer and/or wetter climatic conditions, potentially increasing the net loss of carbon to the atmosphere, at a greater degree than previously expected.
The appropriate selection of recurrence thresholds is a key problem in applications of recurrence quantification analysis and related methods across disciplines. Here, we discuss the distribution of pairwise distances between state vectors in the studied system’s state space reconstructed by means of time-delay embedding as the key characteristic that should guide the corresponding choice for obtaining an adequate resolution of a recurrence plot. Specifically, we present an empirical description of the distance distribution, focusing on characteristic changes of its shape with increasing embedding dimension. Our results suggest that selecting the recurrence threshold according to a fixed percentile of this distribution reduces the dependence of recurrence characteristics on the embedding dimension in comparison with other commonly used threshold selection methods. Numerical investigations on some paradigmatic model systems with time-dependent parameters support these empirical findings.
Recurrence plots (RPs) provide an intuitive tool for visualizing the (potentially multi-dimensional) trajectory of a dynamical system in state space. In case only univariate observations of the system’s overall state are available, time-delay embedding has become a standard procedure for qualitatively reconstructing the dynamics in state space. The selection of a threshold distance 𝜀
, which distinguishes close from distant pairs of (reconstructed) state vectors, is known to have a substantial impact on the recurrence plot and its quantitative characteristics, but its corresponding interplay with the embedding dimension has not yet been explicitly addressed. Here, we point out that the results of recurrence quantification analysis (RQA) and related methods are qualitatively robust under changes of the (sufficiently high) embedding dimension only if the full distribution of pairwise distances between state vectors is considered for selecting 𝜀, which is achieved by consideration of a fixed recurrence rate.
The use of chironomids (Diptera: Chironomidae) as indicators of changes in natural and climatic settings allows one to reconstruct paleoclimatic and paleoenvironmental conditions. Based on the results of lithological and chironomid analyses, the loss on ignition (LOI) values, as well as the results of radiocarbon dating (C-14 AMS), new data on the paleoclimate in the Karelian Isthmus, in particular the quantitative reconstruction of the mean July temperature (T (July), A degrees C) in the Late Pleistocene and Holocene and the evolution of the ecosystem of Medvedevskoe Lake have been obtained.
We present the first high-resolution trace element (Mg/Ca, Sr/Ca, Ba/Ca) record from a stalagmite in southwestern Romania covering the last 3.6 ka, which provides the potential for quantitative climate reconstruction. Precise age control is based on three independent dating methods, in particular for the last 250 yr, where chemical lamina counting is combined with the identification of the 20th century radiocarbon bomb peak and Th-230/U dating. Long-term cave monitoring and model simulations of drip water and speleothem elemental variability indicate that precipitation-related processes are the main drivers of speleothem Mg/Ca ratios. Calibration against instrumental climate data shows a significant anti-correlation of speleothem Mg/Ca ratios with autumn/winter (October to March) precipitation (r = -0.61, p < 0.01), which is statistically robust when considering age uncertainties and auto-correlation. This relationship is used to develop a quantitative reconstruction of autumn/winter precipitation. During the late Holocene, our data suggest a heterogeneous pattern of past regional winter hydroclimate in the Carpathian/Balkan realm, along with intermittent weakening of the dominant influence of North Atlantic forcing. In agreement with other regional paleo-hydrological reconstructions, the observed variability reveals periodically occurring strong NW-SE hydro-climate gradients. We hypothesize, that this pattern is caused by shifts of the eastern edge of the area of influence of the NAO across central eastern Europe due to the confluence of North Atlantic forcing, and other climatic features such as the East Atlantic/Western Russia (EAWR) pattern. (C) 2018 Elsevier B.V. All rights reserved.
After a comprehensive geophysical prospecting the Quaternary MA 1/2 tina Maar, located on a line between the two Quaternary scoria cones Komorni could be revealed by a scientific drilling at the German-Czech border in 2007. Further geophysical field investigations led to the discovery of another geological structure about 2.5 km ESE of the small town Neualbenreuth (NE-Bavaria, Germany), inferred to be also a maar structure, being the fourth volcanic feature aligned along the NW-SE trending Tachov fault zone. It is only faintly indicated as a partial circular rim in the digital elevation model. Though not expressed by a clear magnetic anomaly, geoelectric and refraction seismic tomography strongly indicates a bowl-shaped depression filled with low-resistivity and low-velocity material, correlating well with the well-defined negative gravity anomaly of - 2.5 mGal. Below ca. 15 m-thick debris layer, successions of mostly laminated sediments were recovered in a 100 m-long sediment core in 2015. Sections of finely laminated layers, likely varves, rich in organic matter and tree pollen, were recognized in the upper (22-30 m) and lower (70-86 m) part of the core, respectively, interpreted as interglacials, whereas mostly minerogenic laminated deposits, poor in organic matter, and (almost) barren of tree pollen are interpreted as clastic glacial deposits. According to a preliminary age model based on magnetostratigraphy, palynology, radiocarbon dating, and cyclostratigraphy, the recovered sediments span the time window from about 85 ka back to about 270 ka, covering marine isotope stages 5-8. Sedimentation rates are in the range of 10 cm ka(-1) in interglacials and up to 100 cm ka(-1) in glacial phases. The stratigraphic record resembles the one from MA 1/2 tina Maar, with its eruption date being derived from a nearby tephra deposit at 288 +/- 17 ka, thus supporting the age model of the inferred Neualbenreuth Maar.
The Sesia zone in the Italian Western Alps is a piece of continental crust that has been subducted to eclogite-facies conditions and records a complex metamorphic history. The exact timing of events and the significance of geochronological information are debated due to the interplay of tectonic, metamorphic, and metasomatic processes. Here we present new geochronological data using Rb-Sr internal mineral isochrons and in situ Ar-40/Ar-39 laser ablation data to provide constraints on the relative importance of fluid-mediated mineral replacement reactions and diffusion for the interpretation of radiogenic isotope signatures, and on the use of these isotopic systems for dating metamorphic and variably deformed rocks. Our study focuses on the shear zone at the contact between two major lithological units of the Sesia zone, the eclogitic micaschists and the gneiss minuti. Metasedimentary rocks of the eclogitic micaschists unit contain phengite with step-like zoning in major element chemistry as evidence for petrologic disequilibrium. Distinct Ar-40/Ar-39 spot ages of relict phengite cores and over-printed rims demonstrate the preservation of individual age domains in the crystals. The eclogitic micaschists also show systematic Sr isotope disequilibria among different phengite populations, so that minimum ages of relict assemblage crystallization can be differentiated from the timing of late increments of deformation. The preservation of these disequilibrium features shows the lack of diffusive re-equilibration and underpins that fluid-assisted dissolution and recrystallization reactions are the main factors controlling the isotope record in these subduction-related metamorphic rocks. Blueschist-facies mylonites record deformation along the major shear zone that separates the eclogitic micaschists from the gneiss minuti. Two Rb-Sr isochrones that comprise several white mica fractions and glaucophane constrain the timing of this deformation and accompanying near-complete blueschist-facies re-equilibration of the Rb-Sr system to 60.1 +/- 0.9 Ma and 60.9 +/- 2.1 Ma, respectively. Overlapping ages in eclogitic micaschists of 60.1 +/- 1.1 (Rb-Sr isochron of sheared matrix assemblage), 58.6 +/- 0.8, and 60.9 +/- 0.4 Ma (white mica Ar-40/Ar-39 inverse isochron ages) support the significance of this age and show that fluid-rock interaction and partial re-equilibration occurred as much as several kilometers away from the shear zone. An earlier equilibration during high-pressure conditions in the eclogitic mica schists is recorded in minimum Rb-Sr ages for relict assemblages (77.2 +/- 0.8 and 72.4 +/- 1.1 Ma) and an Ar-40/Ar-39 inverse isochron age of 75.4 +/- 0.8 Ma for white mica cores, again demonstrating that the two isotope systems provide mutually supporting geochronological information. Local reactivation and recrystallization along the shear zone lasted >15 m.y., as late increments of deformation are recorded in a greenschist-facies mylonite by a Rb-Sr isochron age of 46.5 +/- 0.7 Ma.
The variation of Rayleigh ellipticity versus frequency is gaining popularity in site characterization. It becomes a necessary observable to complement dispersion curves when inverting shear wave velocity profiles. Various methods have been proposed so far to extract polarization from ambient vibrations recorded on a single three-component station or with an array of three-component sensors. If only absolute values were recovered 10 yr ago, new array-based techniques were recently proposed with enhanced efficiencies providing also the ellipticity sign. With array processing, higher-order modes are often detected even in the ellipticity domain. We suggest to explore the properties of a high-resolution beamforming where radial and vertical components are explicitly included. If N is the number of three-component sensors, 2N x 2N cross-spectral density matrices are calculated for all presumed directions of propagation. They are built with N radial and N vertical channels. As a first approach, steering vectors are designed to fit with Rayleigh wave properties: the phase shift between radial and vertical components is either -Pi/2 or Pi/2. We show that neglecting the ellipticity tilt due to attenuation has only minor effects on the results. Additionally, we prove analytically that it is possible to retrieve the ellipticity value from the usual maximization of the high-resolution beam power. The method is tested on synthetic data sets and on experimental data. Both are reference sites already analysed by several authors. A detailed comparison with previous results on these cases is provided.
Flash floods and debris flows are iconic hazards inmountainous regions with steep relief, high rainfall intensities, rapid snowmelt events, and abundant sediments. The cuesta landscapes of southern Germany hardly come to mind when dealing with such hazards. A series of heavy rainstorms dumping up to 140mm in 2 h caused destructive flash floods and debris flows in May 2016. The most severe damage occurred in the Braunsbach municipality, which was partly buried by 42,000 m(3) of boulders, gravel, mud, and anthropogenic debris from the small catchment of Orlacher Bach (similar to 6 km(2)). We analysed this event by combining rainfall patterns, geological conditions, and geomorphic impacts to estimate an average sediment yield of 14,000 t/km(2) that mostly (similar to 95%) came from some 50 riparian landslides and channel-bed incision of similar to 2 m. This specific sediment yield ranks among the top 20% globally, while the intensity-duration curve of the rainstormis similarly in the upper percentile range of storms that had triggered landslides. Compared to similar-sized catchments in the greater region hit by the rainstorms, we find that the Orlacher Bach is above the 95th percentile in terms of steepness, storm-rainfall intensity, and topographic curvatures. The flash flood transported a sediment volume equal to as much as 20-40% of the Pleistocene sediment volume stored in the Orlacher Bach fan, andmay have had several predecessors in the Holocene. River control structures from 1903 and records of a debris flow in the 1920s in a nearby catchment indicate that the local inhabitants may have been aware of the debris-flow hazards earlier. Such recurring and destructive events elude flood-hazard appraisals in humid landscapes of gentle relief, and broaden mechanistic views of how landslides and debris flows contribute to shaping small and deeply cut tributaries in the southern Germany cuesta landscape.
The equatorial electrojet is an enhanced eastward current in the dayside E region ionosphere flowing along the magnetic equator. The equatorial electrojet is highly variable as it is subject to various forcing mechanisms including atmospheric waves from the lower layers of the atmosphere. There are occasionally times when the intensity of the equatorial electrojet at a fixed longitude shows an oscillatory variation with a period of approximately 6days. We present case studies of such events based on the equatorial electrojet measurements from the CHAMP and Swarm satellites. The spatial and temporal variability of the equatorial electrojet intensity during these events reveals characteristics of a westward propagating wave with zonal wavenumber 1, consistent with the effect of the quasi-6-day planetary wave. Analyses of the geopotential height data from the Aura satellite confirm the presence of the quasi-6-day planetary wave in the lower thermosphere during the events. The amplitude of the quasi-6-day variation in the equatorial electrojet intensity depends on longitude, but no systematic longitudinal dependence is found for different events. During the event of August 2010, quasi-6-day variations are also observed by ground-based magnetometers and a radar in the Peruvian sector. The effect of the quasi-6-day wave accounts for up to +/- 5.9m/s in the equatorial vertical plasma velocity at noon, which is much larger than previously predicted by a numerical model. These results suggest that the quasi-6-day planetary wave is an important source of short-term variability in the equatorial ionosphere.
Quantifying the roles of single stations within homogeneous regions using complex network analysis
(2018)
Regionalization and pooling stations to form homogeneous regions or communities are essential for reliable parameter transfer, prediction in ungauged basins, and estimation of missing information. Over the years, several clustering methods have been proposed for regional analysis. Most of these methods are able to quantify the study region in terms of homogeneity but fail to provide microscopic information about the interaction between communities, as well as about each station within the communities. We propose a complex network-based approach to extract this valuable information and demonstrate the potential of our approach using a rainfall network constructed from the Indian gridded daily precipitation data. The communities were identified using the network-theoretical community detection algorithm for maximizing the modularity. Further, the grid points (nodes) were classified into universal roles according to their pattern of within- and between-community connections. The method thus yields zoomed-in details of individual rainfall grids within each community.
Subsurface lateral flow in hillslope soils depends on lower permeability or texture-contrasting soil horizons. In the arable hummocky soil landscape, erosion processes caused glacial till appearance closer to the soil surface at upslope positions. The objective of this work was to quantify the potential for subsurface lateral flow depending on the erosion-affected spatial hydropedological complexity. The eroded Haplic Luvisol profile was studied due to the presence of a relatively dense C horizon that varied in depth, thickness, and sloping angle. A two-dimensional numerical modeling and sensitivity analysis for the saturated hydraulic conductivity (K-s) of the C horizon and the depth to C horizon (i.e., soil solum thickness) was performed for rainstorms in 2011 and 2012 using HYDRUS-2D. A K-s value of <2.5 cm d(-1) for the C horizon was required for lateral flow initiation. Lateral flow was (i) increasing with decreasing solum thickness, indicating an erosion-induced feedback on subsurface lateral flow, and (ii) dependent on the soil moisture regime prior to rainstorms. The effect of lateral flow on the movement of a conservative tracer was simulated in the form of a "virtual experiment". Simulation scenarios revealed only a relatively small lateral shift of the tracer plume caused by a local decoupling of water (already lateral) from subsequent tracer movement (still vertical). Longer term simulations suggested that, for the present conditions, lateral flow was limited mostly to occasional summer storm events. Even without considering preferential flow contribution to lateral flow, highly complex hydropedologic interactions are present in erosion-affected heterogeneous soil profiles.
How information about sediment transport processes is transmitted to the sedimentary record remains a complex problem for the interpretation of fluvial stratigraphy. Alluvial fan deposits represent the condensed archive of sediment transport, which is at least partly controlled by tectonics and climate. For three coupled catchment-fan systems in northern Death Valley, California, we measure grain size across 12 well-preserved Holocene and late-Pleistocene surfaces, mapped in detail from field observations and remote sensing. Our results show that fan surfaces correlated to the late Pleistocene are, on average, 30-50% coarser than active or Holocene fan surfaces. We adopt a self-similar form of grain size distribution based on the observed stability of the ratio between mean grain size and standard deviation downstream. Using statistical analysis, we show that fan surface grain size distributions are self-similar. We derive a relative mobility function using our self-similar grain size distributions, which describes the relative probability of a given grain size being transported. We show that the largest mobile grain sizes are between 20 and 35mm, a value that varies over time and is clearly lower in the Holocene than in the Pleistocene; a change we suggest is due to a drier climate in the Holocene. These results support recent findings that alluvial fan sedimentology can record past environmental change and that these landscapes are potentially sensitive to climatic change over a glacial-interglacial cycle. We demonstrate that the self-similarity methodology offers a means to explore changes in relative mobility of grain sizes from preserved fluvial deposits. Plain Language Summary A key challenge in Earth Science is understanding how landscapes respond to climate. It may be possible to observe measurable differences in certain landscapes settings such as alluvial fans in desert regions. Alluvial fans are believed to be effective recorders of climate, representing a cumulative store of material transported downstream by rainfall-sensitive river systems. In northern Death Valley, California, we measure at high resolution grain size on three alluvial fans with surfaces that date from the Holocene and the arid climate of today to the 20-40% wetter late-Pleistocene epoch. We find that older late-Pleistocene surfaces are coarser on average than surfaces deposited during the modern and Holocene dry period, suggesting a changing sediment transport regime potentially in response to precipitation. We also show that measured grain size distributions within and between surfaces can be successfully normalized based on the decay in mean grain size and variance downstream, exhibiting a self-similar pattern. Finally, we employ a grain size relative mobility model using our field data to establish which grain sizes are likely to be in transport or locked in the substrate. This model predicts that during the wetter late-Pleistocene mobile grain sizes are up to 40% larger than during the Holocene.
Quantifying rock weakening due to decreasing calcite mineral content by numerical simulations
(2018)
The quantification of changes in geomechanical properties due to chemical reactions is of paramount importance for geological subsurface utilisation, since mineral dissolution generally reduces rock stiffness. In the present study, the effective elastic moduli of two digital rock samples, the Fontainebleau and Bentheim sandstones, are numerically determined based on micro-CT images. Reduction in rock stiffness due to the dissolution of 10% calcite cement by volume out of the pore network is quantified for three synthetic spatial calcite distributions (coating, partial filling and random) using representative sub-cubes derived from the digital rock samples. Due to the reduced calcite content, bulk and shear moduli decrease by 34% and 38% in maximum, respectively. Total porosity is clearly the dominant parameter, while spatial calcite distribution has a minor impact, except for a randomly chosen cement distribution within the pore network. Moreover, applying an initial stiffness reduced by 47% for the calcite cement results only in a slightly weaker mechanical behaviour. Using the quantitative approach introduced here substantially improves the accuracy of predictions in elastic rock properties compared to general analytical methods, and further enables quantification of uncertainties related to spatial variations in porosity and mineral distribution.
Quantifying rock weakening due to decreasing calcite mineral content by numerical simulations
(2018)
The quantification of changes in geomechanical properties due to chemical reactions is of paramount importance for geological subsurface utilisation, since mineral dissolution generally reduces rock stiffness. In the present study, the effective elastic moduli of two digital rock samples, the Fontainebleau and Bentheim sandstones, are numerically determined based on micro-CT images. Reduction in rock stiffness due to the dissolution of 10% calcite cement by volume out of the pore network is quantified for three synthetic spatial calcite distributions (coating, partial filling and random) using representative sub-cubes derived from the digital rock samples. Due to the reduced calcite content, bulk and shear moduli decrease by 34% and 38% in maximum, respectively. Total porosity is clearly the dominant parameter, while spatial calcite distribution has a minor impact, except for a randomly chosen cement distribution within the pore network. Moreover, applying an initial stiffness reduced by 47% for the calcite cement results only in a slightly weaker mechanical behaviour. Using the quantitative approach introduced here substantially improves the accuracy of predictions in elastic rock properties compared to general analytical methods, and further enables quantification of uncertainties related to spatial variations in porosity and mineral distribution.
Together with the gradual change of mean values, ongoing climate change is projected to increase frequency and amplitude of temperature and precipitation extremes in many regions of Europe. The impacts of such in most cases short term extraordinary climate situations on terrestrial ecosystems are a matter of central interest of recent climate change research, because it can not per se be assumed that known dependencies between climate variables and ecosystems are linearly scalable. So far, yet, there is a high demand for a method to quantify such impacts in terms of simultaneities of event time series.
In the course of this manuscript the new statistical approach of Event Coincidence Analysis (ECA) as well as it's R implementation is introduced, a methodology that allows assessing whether or not two types of event time series exhibit similar sequences of occurrences. Applications of the method are presented, analyzing climate impacts on different temporal and spacial scales: the impact of extraordinary expressions of various climatic variables on tree stem variations (subdaily and local scale), the impact of extreme temperature and precipitation events on the owering time of European shrub species (weekly and country scale), the impact of extreme temperature events on ecosystem health in terms of NDVI (weekly and continental scale) and the impact of El Niño and La Niña events on precipitation anomalies (seasonal and global scale).
The applications presented in this thesis refine already known relationships based on classical methods and also deliver substantial new findings to the scientific community: the widely known positive correlation between flowering time and temperature for example is confirmed to be valid for the tails of the distributions while the widely assumed positive dependency between stem diameter variation and temperature is shown to be not valid for very warm and very cold days. The larger scale investigations underline the sensitivity of anthrogenically shaped landscapes towards temperature extremes in Europe and provide a comprehensive global ENSO impact map for strong precipitation events.
Finally, by publishing the R implementation of the method, this thesis shall enable other researcher to further investigate on similar research questions by using Event Coincidence Analysis.
Mountain rivers respond to strong earthquakes by rapidly aggrading to accommodate excess sediment delivered by co-seismic landslides. Detailed sediment budgets indicate that rivers need several years to decades to recover from seismic disturbances, depending on how recovery is defined. We examine three principal proxies of river recovery after earthquake-induced sediment pulses around Pokhara, Nepal's second largest city. Freshly exhumed cohorts of floodplain trees in growth position indicate rapid and pulsed sedimentation that formed a fan covering 150 km2 in a Lesser Himalayan basin with tens of metres of debris between the 11th and 15th centuries AD. Radiocarbon dates of buried trees are consistent with those of nearby valley deposits linked to major medieval earthquakes, such that we can estimate average rates of re-incision since. We combine high-resolution digital elevation data, geodetic field surveys, aerial photos, and dated tree trunks to reconstruct geomorphic marker surfaces. The volumes of sediment relative to these surfaces require average net sediment yields of up to 4200 t km–2 yr–1 for the 650 years since the last inferred earthquake-triggered sediment pulse. The lithological composition of channel bedload differs from that of local bedrock, confirming that rivers are still mostly evacuating medieval valley fills, locally incising at rates of up to 0.2 m yr–1. Pronounced knickpoints and epigenetic gorges at tributary junctions further illustrate the protracted fluvial response; only the distal portions of the earthquake-derived sediment wedges have been cut to near their base. Our results challenge the notion that mountain rivers recover speedily from earthquakes within years to decades. The valley fills around Pokhara show that even highly erosive Himalayan rivers may need more than several centuries to adjust to catastrophic perturbations. Our results motivate some rethinking of post-seismic hazard appraisals and infrastructural planning in active mountain regions.
Evidence for Quaternary climate change in East Africa has been derived from outcrops on land and lake cores and from marine dust, leaf wax, and pollen records. These data have previously been used to evaluate the impact of climate change on hominin evolution, but correlations have proved to be difficult, given poor data continuity and the great distances between marine cores and terrestrial basins where fossil evidence is located. Here, we present continental coring evidence for progressive aridification since about 575 thousand years before present (ka), based on Lake Magadi (Kenya) sediments. This long-term drying trend was interrupted by many wet-dry cycles, with the greatest variability developing during times of high eccentricity-modulated precession. Intense aridification apparent in the Magadi record took place between 525 and 400 ka, with relatively persistent arid conditions after 350 ka and through to the present. Arid conditions in the Magadi Basin coincide with the Mid-Brunhes Event and overlap with mammalian extinctions in the South Kenya Rift between 500 and 400 ka. The 525 to 400 ka arid phase developed in the South Kenya Rift between the period when the last Acheulean tools are reported (at about 500 ka) and before the appearance of Middle Stone Age artifacts (by about 320 ka). Our data suggest that increasing Middle- to Late-Pleistocene aridification and environmental variability may have been drivers in the physical and cultural evolution of Homo sapiens in East Africa.
Preface
(2018)
Here we present high-resolution delta O-18 records obtained from speleothems collected in the eastern Bolivian Andes. The stable isotope records are related to the regional- to large-scale atmospheric circulation over South America and allow interpreting changes in delta O-18 during the last 1400 yr as a function of changes in precipitation regimes over the southern tropical Andes. Two distinct phases with more negative delta O-18 values, interpreted as periods of increased convective activity over the eastern Andean Cordillera in Bolivia are observed concomitantly with periods of global climate anomalies during the last millennium, such as the Medieval Climate Anomaly (MCA) and the Little Ice Age (LIA) respectively. Changes in the Bolivian delta O-18 record during the LIA are apparently related to a southward displacement of the Intertropical Convergence Zone (ITCZ), which acts as a main moisture driver to intensify convection over the tropical continent. During the MCA, however, the increased convective activity observed in the Bolivian record is likely the result of a different mechanism, which implies moisture sourced mainly from the southern tropical Atlantic. This interpretation is consistent with paleoclimate records further to the north in the tropical Andes that show progressively drier conditions during this time period, indicating a more northerly position of the ITCZ. The transition period between the MCA and the LIA shows a slight tendency toward increased delta O-18 values, indicating weakened convective activity. Our results also reveal a non-stationary anti-phased behavior between the delta O-18 reconstructions from Bolivia and northeastern Brazil that confirms a continental-scale east-west teleconnection across South America during the LIA.