Refine
Year of publication
- 2016 (292) (remove)
Document Type
- Article (221)
- Doctoral Thesis (37)
- Other (14)
- Postprint (13)
- Review (5)
- Habilitation Thesis (1)
- Report (1)
Keywords
- climate change (5)
- erosion (5)
- water balance (5)
- Pollen (4)
- August 2002 flood (3)
- Body waves (3)
- Central Europe (3)
- Floods Directive (3)
- Himalaya (3)
- Himalayas (3)
Institute
- Institut für Geowissenschaften (292) (remove)
Mineral topsoils possess large organic carbon (OC) contents but there is only limited knowledge on the mechanisms controlling the preservation of organic matter (OM) against microbial decay. Samples were taken from the uppermost mineral topsoil horizon (0 to 5 cm) of seven sites under mature deciduous forest showing OC contents between 69 and 164 g kg(-1) and a wide range in mineral characteristics. At first, organic particles and the water-extractable OM were removed from the soil samples. Thereafter, Na-pyrophosphate extractable organic matter (OM(PY)), assumed to be indicative for OM bound via cation mediated interactions, and the OM remaining in the extraction residue (OM(ER)), supposed to be indicative for OM occluded in mechanically highly stable micro-aggregates, were sequentially separated and quantified. The composition of OM(PY) and OM(ER) was analyzed by FTIR and their stability by C-14 measurements. The OC remaining in the extraction residues accounted for 38 to 59% of the bulk soil OC (SOC) suggesting a much larger relevance of OM(ER) for the OM dynamic in the analyzed soils as compared with OM(PY) that accounted for 1.6 to 7.5% of the SOC. The FUR analyses revealed a lower relative proportion of C=O groups in OM(ER) compared to OM(PY) indicating differences in the degree of microbial processing between these fractions. Correlation analyses suggest an increase in the stability of OM(PY) with the soil pH and contents of Na-pyrophosphate soluble Fe, Al, and Mg and an increase in the stability of OM(ER) with the soil pH and the contents of clay and oxalate-soluble Fe and Al. Despite the detected influence of soil mineral characteristics on the turnover of OM(PY) and OM(ER), the Delta C-14 signatures indicated mean residence times less than 100 years. The presence of less stabilized OM in these fractions can be derived from methodological uncertainties and/or the fast cycling compartment of mineral-associated OM. (C) 2015 Elsevier B.V. All rights reserved.
The spring of Nadaouiyeh Ain Askar - Paleoecology of a Paleolithic oasis in arid central Syria
(2016)
The site Nadaouiyeh Aïn Askar, an ancient artesian spring near the village of El Kowm, Central Syria, is an example of long lasting human occupation in a desert environment throughout the Middle and Late Pleistocene. The excavations expose a succession of sedimentary units, containing an artifact assemblage assigned to the Acheulean techno-complex. Unit VI, attributed to the Marine Isotope Stage 13, is rich in ostracod valves and was chosen for the present environmental study. From these sediments Heterocypris salina, H. incongruens, Cyprideis torosa, Ilyocypris cf. bradyi, I. inermis, I. cf. gibba, Darwinula stevensoni, Plesiocypridopsis newtoni, Pseudocandona compressa, Candona cf. neglecta, Pseudocandona sp., Trajancypris sp., Physocypria sp. and Mixtacandona sp. are documented for the first time in the Middle Pleistocene of the arid environment of central Syria. Data from these microfossils as well as geochemical proxies implicate three phases, turning the wetland from a palustrine setting into a spring supplied pond with increasing salinity. The high mineralization of the spring waters enables a discussion about early hominin adaptability to brackish waters as drinking water resources, common within the steppe and desert environments along the “out-of-Africa”—corridor in the eastern Mediterranean.
The aim of the present study was to investigate nutritional and metabolic parameters during the dry and early postpartum periods of ovulatory and anovulatory cows, as well as their postpartum reproductive performance. Blood samples from 20 multiparous Holstein cows were collected once a week from the far-off dry period to 3 weeks postpartum. Early postpartum (0-3 weeks) ovulation was confirmed using plasma progesterone concentration profiles, and cows were considered ovulatory if they had resumed luteal activity by this point (n = 9), whereas cows that had not were considered anovulatory (n = 11). Data from the ovulatory and anovulatory cows were analyzed separately for the far-off dry period (7-4 weeks prepartum), the close-up dry period (3-1 weeks prepartum), and the early postpartum period (0-3 weeks). Serum gamma-glutamyl transpeptidase activity (far-off, P = 0.065; close-up, P = 0.051; and early postpartum, P = 0.030) and aspartate aminotransferase (close-up, P = 0.050 and early postpartum, P = 0.087) activities were higher in anovulatory than in ovulatory cows. The days open period was longer (P = 0.019) in anovulatory than in ovulatory cows, and the number of artificial inseminations per conception (P = 0.025) was greater. In conclusion, we found that continuously high gamma-glutamyl transpeptidase activities in serum, which may be induced by liver disorders, prevent subsequent ovulation and affect subsequent fertility, even if cows obtain sufficient ovulation-related energy and beta-carotene.
Meteorological extreme events have great potential for damaging railway infrastructure and posing risks to the safety of train passengers. In the future, climate change will presumably have serious implications on meteorological hazards in the Alpine region. Hence, attaining insights on future frequencies of meteorological extremes with relevance for the railway operation in Austria is required in the context of a comprehensive and sustainable natural hazard management plan of the railway operator. In this study, possible impacts of climate change on the frequencies of so-called critical meteorological conditions (CMCs) between the periods 1961-1990 and 2011-2040 are analyzed. Thresholds for such CMCs have been defined by the railway operator and used in its weather monitoring and early warning system. First, the seasonal climate change signals for air temperature and precipitation in Austria are described on the basis of an ensemble of high-resolution Regional Climate Model (RCM) simulations for Europe. Subsequently, the RCM-ensemble was used to investigate changes in the frequency of CMCs. Finally, the sensitivity of results is analyzed with varying threshold values for the CMCs. Results give robust indications for an all-season air temperature rise, but show no clear tendency in average precipitation. The frequency analyses reveal an increase in intense rainfall events and heat waves, whereas heavy snowfall and cold days are likely to decrease. Furthermore, results indicate that frequencies of CMCs are rather sensitive to changes of thresholds. It thus emphasizes the importance to carefully define, validate, andif neededto adapt the thresholds that are used in the weather monitoring and warning system of the railway operator. For this, continuous and standardized documentation of damaging events and near-misses is a pre-requisite.
Experience has shown that river floods can significantly hamper the reliability of railway networks and cause extensive structural damage and disruption. As a result, the national railway operator in Austria had to cope with financial losses of more than EUR 100 million due to flooding in recent years. Comprehensive information on potential flood risk hot spots as well as on expected flood damage in Austria is therefore needed for strategic flood risk management. In view of this, the flood damage model RAIL (RAilway Infrastructure Loss) was applied to estimate (1) the expected structural flood damage and (2) the resulting repair costs of railway infrastructure due to a 30-, 100- and 300-year flood in the Austrian Mur River catchment. The results were then used to calculate the expected annual damage of the railway subnetwork and subsequently analysed in terms of their sensitivity to key model assumptions. Additionally, the impact of risk aversion on the estimates was investigated, and the overall results were briefly discussed against the background of climate change and possibly resulting changes in flood risk. The findings indicate that the RAIL model is capable of supporting decision-making in risk management by providing comprehensive risk information on the catchment level. It is furthermore demonstrated that an increased risk aversion of the railway operator has a marked influence on flood damage estimates for the study area and, hence, should be considered with regard to the development of risk management strategies.
Tectonic blocks and slabs of mafic-ultramafic rocks are distributed discontinuously in the Yuli metamorphic belt of Taiwan. The blocks include rare omphacite metagabbros and gamet-epidote blueschists in the Wanjung and Juisui (Tamayen) areas, respectively. Such high-pressure (HP) mineral assemblages have been attributed to a mid-Miocene subduction event. However, the surrounding psammitic, pelitic and chloritic schists are the dominant greenschist-facies lithologies of the Yuli belt. In the Chinshuichi area, tectonic blocks are enclosed in garnet-bearing metapelites, suggesting elevated pressures. In this area, we recently discovered meta-plagiogranite containing the assemblage glaucophane + omphacite (X-Jd up to 039) + rutile + quartz, indicating P-T conditions near 13 kbar/550 degrees C. New equilibrium phase modeling of a garnet-paragonite mica schist and compositional isopleths for peak assemblage minerals garnet and phengite (Si = 3.33-3.37 pfu) indicate metamorphic conditions of 15.5-17 kbar/530-550 degrees C. These P-T estimates are higher than previously reported in the Yuli belt and suggest that both tectonic blocks and host metapelites underwent HP metamorphism. The juxtaposition of tectonic blocks and metapelites apparently occurred during the formation of a subduction-accretionary complex, followed by exhumation facilitated by a collisional event. These new findings imply that HP metamorphism was not limited to tectonic blocks, and instead played a significant role attending orogenesis in eastern Taiwan. (C) 2015 Elsevier B.V. All rights reserved.
Storm runoff from the Marikina River Basin frequently causes flood events in the Philippine capital region Metro Manila. This paper presents and evaluates a system to predict short-term runoff from the upper part of that basin (380km(2)). It was designed as a possible component of an operational warning system yet to be installed. For the purpose of forecast verification, hindcasts of streamflow were generated for a period of 15 months with a time-continuous, conceptual hydrological model. The latter was fed with real-time observations of rainfall. Both ground observations and weather radar data were tested as rainfall forcings. The radar-based precipitation estimates clearly outperformed the raingauge-based estimates in the hydrological verification. Nevertheless, the quality of the deterministic short-term runoff forecasts was found to be limited. For the radar-based predictions, the reduction of variance for lead times of 1, 2 and 3hours was 0.61, 0.62 and 0.54, respectively, with reference to a no-forecast scenario, i.e. persistence. The probability of detection for major increases in streamflow was typically less than 0.5. Given the significance of flood events in the Marikina Basin, more effort needs to be put into the reduction of forecast errors and the quantification of remaining uncertainties.
Quantitative thermodynamic and geochemical modeling is today applied in a variety of geological environments from the petrogenesis of igneous rocks to the oceanic realm. Thermodynamic calculations are used, for example, to get better insight into lithosphere dynamics, to constrain melting processes in crust and mantle as well as to study fluid-rock interaction. The development of thermodynamic databases and computer programs to calculate equilibrium phase diagrams have greatly advanced our ability to model geodynamic processes from subduction to orogenesis. However, a well-known problem is that despite its broad application the use and interpretation of thermodynamic models applied to natural rocks is far from straightforward. For example, chemical disequilibrium and/or unknown rock properties, such as fluid activities, complicate the application of equilibrium thermodynamics.
One major aspect of the publications presented in this Habilitationsschrift are new approaches to unravel dynamic and chemical histories of rocks that include applications to chemically open system behaviour. This approach is especially important in rocks that are affected by element fractionation due to fractional crystallisation and fluid loss during dehydration reactions. Furthermore, chemically open system behaviour has also to be considered for studying fluid-rock interaction processes and for extracting information from compositionally zoned metamorphic minerals. In this Habilitationsschrift several publications are presented where I incorporate such open system behaviour in the forward models by incrementing the calculations and considering changing reacting rock compositions during metamorphism. I apply thermodynamic forward modelling incorporating the effects of element fractionation in a variety of geodynamic and geochemical applications in order to better understand lithosphere dynamics and mass transfer in solid rocks.
In three of the presented publications I combine thermodynamic forward models with trace element calculations in order to enlarge the application of geochemical numerical forward modeling. In these publications a combination of thermodynamic and trace element forward modeling is used to study and quantify processes in metamorphic petrology at spatial scales from µm to km. In the thermodynamic forward models I utilize Gibbs energy minimization to quantify mineralogical changes along a reaction path of a chemically open fluid/rock system. These results are combined with mass balanced trace element calculations to determine the trace element distribution between rock and melt/fluid during the metamorphic evolution. Thus, effects of mineral reactions, fluid-rock interaction and element transport in metamorphic rocks on the trace element and isotopic composition of minerals, rocks and percolating fluids or melts can be predicted.
One of the included publications shows that trace element growth zonations in metamorphic garnet porphyroblasts can be used to get crucial information about the reaction path of the investigated sample. In order to interpret the major and trace element distribution and zoning patterns in terms of the reaction history of the samples, we combined thermodynamic forward models with mass-balance rare earth element calculations. Such combined thermodynamic and mass-balance calculations of the rare earth element distribution among the modelled stable phases yielded characteristic zonation patterns in garnet that closely resemble those in the natural samples. We can show in that paper that garnet growth and trace element incorporation occurred in near thermodynamic equilibrium with matrix phases during subduction and that the rare earth element patterns in garnet exhibit distinct enrichment zones that fingerprint the minerals involved in the garnet-forming reactions.
In two of the presented publications I illustrate the capacities of combined thermodynamic-geochemical modeling based on examples relevant to mass transfer in subduction zones. The first example focuses on fluid-rock interaction in and around a blueschist-facies shear zone in felsic gneisses, where fluid-induced mineral reactions and their effects on boron (B) concentrations and isotopic compositions in white mica are modeled. In the second example, fluid release from a subducted slab and associated transport of B and variations in B concentrations and isotopic compositions in liberated fluids and residual rocks are modeled. I show that, combined with experimental data on elemental partitioning and isotopic fractionation, thermodynamic forward modeling unfolds enormous capacities that are far from exhausted.
In my publications presented in this Habilitationsschrift I compare the modeled results to geochemical data of natural minerals and rocks and demonstrate that the combination of thermodynamic and geochemical models enables quantification of metamorphic processes and insights into element cycling that would have been unattainable so far.
Thus, the contributions to the science community presented in this Habilitatonsschrift concern the fields of petrology, geochemistry, geochronology but also ore geology that all use thermodynamic and geochemical models to solve various problems related to geo-materials.
The subduction of hydrated slab mantle is the most important and yet weakly constrained factor in the quantification of the Earth's deep geologic water cycle. The most critical unknowns are the initial hydration state and the dehydration behavior of the subducted oceanic mantle. Here we present a combined thermomechanical, thermodynamic, and geochemical model of the Kamchatka subduction zone that indicates significant dehydration of subducted slab mantle beneath Kamchatka. Evidence for the subduction of hydrated oceanic mantle comes from across-arc trends of boron concentrations and isotopic compositions in arc volcanic rocks. Our thermodynamic-geochemical models successfully predict the complex geochemical patterns and the spatial distribution of arc volcanoes in Kamchatka assuming the subduction of hydrated oceanic mantle. Our results show that water content and dehydration behavior of the slab mantle beneath Kamchatka can be directly linked to compositional features in arc volcanic rocks. Depending on hydration depth of the slab mantle, our models yield water recycling rates between 1.1 × 103 and 7.4 × 103 Tg/Ma/km corresponding to values between 0.75 × 106 and 5.2 × 106 Tg/Ma for the entire Kamchatkan subduction zone. These values are up to one order of magnitude lower than previous estimates for Kamchatka, but clearly show that subducted hydrated slab mantle significantly contributes to the water budget in the Kamchatkan subduction zone.
The Dead Sea region has faced substantial environmental challenges in recent decades, including water resource scarcity, similar to 1 m annual decreases in the water level, sinkhole development, ascending-brine freshwater pollution, and seismic disturbance risks. Natural processes are significantly affected by human interference as well as by climate change and tectonic developments over the long term. To get a deep understanding of processes and their interactions, innovative scientific approaches that integrate disciplinary research and education are required. The research project DESERVE (Helmholtz Virtual Institute Dead Sea Research Venue) addresses these challenges in an interdisciplinary approach that includes geophysics, hydrology, and meteorology. The project is implemented by a consortium of scientific institutions in neighboring countries of the Dead Sea (Israel, Jordan, Palestine Territories) and participating German Helmholtz Centres (KIT, GFZ, UFZ). A new monitoring network of meteorological, hydrological, and seismic/geodynamic stations has been established, and extensive field research and numerical simulations have been undertaken. For the first time, innovative measurement and modeling techniques have been applied to the extreme conditions of the Dead Sea and its surroundings. The preliminary results show the potential of these methods. First time ever performed eddy covariance measurements give insight into the governing factors of Dead Sea evaporation. High-resolution bathymetric investigations reveal a strong correlation between submarine springs and neo-tectonic patterns. Based on detailed studies of stratigraphy and borehole information, the extension of the subsurface drainage basin of the Dead Sea is now reliably estimated. Originality has been achieved in monitoring flash floods in an arid basin at its outlet and simultaneously in tributaries, supplemented by spatio-temporal rainfall data. Low-altitude, high resolution photogrammetry, allied to satellite image analysis and to geophysical surveys (e.g. shear-wave reflections) has enabled a more detailed characterization of sinkhole morphology and temporal development and the possible subsurface controls thereon. All the above listed efforts and scientific results take place with the interdisciplinary education of young scientists. They are invited to attend joint thematic workshops and winter schools as well as to participate in field experiments. (C) 2015 The Authors. Published by Elsevier B.V.
Flood loss modeling is an important component for risk analyses and decision support in flood risk management. Commonly, flood loss models describe complex damaging processes by simple, deterministic approaches like depth-damage functions and are associated with large uncertainty. To improve flood loss estimation and to provide quantitative information about the uncertainty associated with loss modeling, a probabilistic, multivariable Bagging decision Tree Flood Loss Estimation MOdel (BT-FLEMO) for residential buildings was developed. The application of BT-FLEMO provides a probability distribution of estimated losses to residential buildings per municipality. BT-FLEMO was applied and validated at the mesoscale in 19 municipalities that were affected during the 2002 flood by the River Mulde in Saxony, Germany. Validation was undertaken on the one hand via a comparison with six deterministic loss models, including both depth-damage functions and multivariable models. On the other hand, the results were compared with official loss data. BT-FLEMO outperforms deterministic, univariable, and multivariable models with regard to model accuracy, although the prediction uncertainty remains high. An important advantage of BT-FLEMO is the quantification of prediction uncertainty. The probability distribution of loss estimates by BT-FLEMO well represents the variation range of loss estimates of the other models in the case study.
Water research is introduced from the combined perspectives of natural and social science and cases of citizen and stakeholder coproduction of knowledge. Using the overarching notion of transdisciplinarity, we examine how interdisciplinary and participatory water research has taken place and could be developed further. It becomes apparent that water knowledge is produced widely within society, across certified disciplinary experts and noncertified expert stakeholders and citizens. However, understanding and management interventions may remain partial, or even conflicting, as much research across and between traditional disciplines has failed to integrate disciplinary paradigms due to philosophical, methodological, and communication barriers. We argue for more agonistic relationships that challenge both certified and noncertified knowledge productively. These should include examination of how water research itself embeds and is embedded in social context and performs political work. While case studies of the cultural and political economy of water knowledge exist, we need more empirical evidence on how exactly culture, politics, and economics have shaped this knowledge and how and at what junctures this could have turned out differently. We may thus channel the coproductionist critique productively to bring perspectives, alternative knowledges, and implications into water politics where they were not previously considered; in an attempt to counter potential lock-in to particular water policies and technologies that may be inequitable, unsustainable, or unacceptable. While engaging explicitly with politics, transdisciplinary water research should remain attentive to closing down moments in the research process, such as framings, path-dependencies, vested interests, researchers’ positionalities, power, and scale.
A partially non-ergodic ground-motion prediction equation is estimated for Europe and the Middle East. Therefore, a hierarchical model is presented that accounts for regional differences. For this purpose, the scaling of ground-motion intensity measures is assumed to be similar, but not identical in different regions. This is achieved by assuming a hierarchical model, where some coefficients are treated as random variables which are sampled from an underlying global distribution. The coefficients are estimated by Bayesian inference. This allows one to estimate the epistemic uncertainty in the coefficients, and consequently in model predictions, in a rigorous way. The model is estimated based on peak ground acceleration data from nine different European/Middle Eastern regions. There are large differences in the amount of earthquakes and records in the different regions. However, due to the hierarchical nature of the model, regions with only few data points borrow strength from other regions with more data. This makes it possible to estimate a separate set of coefficients for all regions. Different regionalized models are compared, for which different coefficients are assumed to be regionally dependent. Results show that regionalizing the coefficients for magnitude and distance scaling leads to better performance of the models. The models for all regions are physically sound, even if only very few earthquakes comprise one region.
The Ms ∼ 7.7 Sarez-Pamir earthquake of 1911 February 18 is the largest instrumentally recorded earthquake in the Pamir region. It triggered one of the largest landslides of the past century, building a giant natural dam and forming Lake Sarez. As for many strong earthquakes from that time, information about source parameters of the Sarez-Pamir earthquake is limited due to the sparse observations. Here, we present the analysis of analogue seismic records of the Sarez-Pamir earthquake. We have collected, scanned and digitized 26 seismic records from 13 stations worldwide to relocate the epicentre and determine the event's depth (∼26 km) and magnitude (mB7.3 and Ms7.7). The unusually good quality of the digitized waveforms allowed their modelling, revealing an NE-striking sinistral strike-slip focal mechanism in accordance with regional tectonics. The shallow depth and magnitude (Mw7.3) of the earthquake were confirmed. Additionally, we investigated the possible contribution of the landslide to the waveforms and present an alternative source model assuming the landslide and earthquake occurred in close sequence.
Precision horticulture encompasses site- or tree-specific management in fruit plantations. Of decisive importance is spatially resolved data (this means data from each tree) from the production site, since it may enable customized and, therefore, resource-efficient production measures.
The present thesis involves an examination of the apparent electrical conductivity of the soil (ECa), the plant water status spatially measured by means of the crop water stress index (CWSI), and the fruit quality (e.g. fruit size) for Prunus domestica L. (plums) and Citrus x aurantium, Syn. Citrus paradisi (grapefruit). The goals of the present work were i) characterization of the 3D distribution of the apparent electrical conductivity of the soil and variability of the plant’s water status; ii) investigation of the interaction between ECa, CWSI, and fruit quality; and iii) an approach for delineating management zones with respect to managing trees individually.
To that end, the main investigations took place in the plum orchard. This plantation got a slope of 3° grade on Pleistocene and post-Pleistocene substrates in a semi-humid climate (Potsdam, Germany) and encloses an area of 0.37 ha with 156 trees of the cultivar ˈTophit Plusˈ on a Wavit rootstock. The plantation was laid in 2009 with annual and biannual trees spaced 4 m distance along the irrigation system and 5 m between the rows. The trees were watered three times a week with a drip irrigation system positioned 50 cm above ground level providing 1.6 l per tree per event. With the help of geoelectric measurements, the apparent electrical conductivity of the upper soil (0.25 m) was measured for each tree with an electrode spacing of 0.5 m (4-point light hp). In this manner, the plantation was spatially charted with respect to the soil’s ECa. Additionally, tomography measurements were performed for 3D mapping of the soil ECa and spot checks of drilled cores with a profile of up to 1 m. The vegetative, generative, and fruit quality data were collected for each tree. The instantaneous plant water status was comprehensively determined in spot checks with the established Scholander method for water potential analysis (Scholander pressure bomb) as well as thermal imaging. An infrared camera was used for the thermal imaging (ThermaCam SC 500), mounted on a tractor 3.3 m above ground level. The thermal images (320 x 240 px) of the canopy surface were taken with an aperture of 45° and a geometric resolution of 8.54 x 6.41 mm. With the aid of the canopy temperature readings from the thermal images, cross-checked with manual temperature measurements of a dry and a wet reference leaf, the crop water stress index (CWSI) was calculated. Adjustments in CWSI for measurements in a semi-humid climate were developed, whereas the collection of reference temperatures was automatically collected from thermal images.
The bonitur data were transformed with the help of a variance stabilization process into a normal distribution. The statistical analyses as well as the automatic evaluation routine were performed with several scripts in MATLAB® (R2010b and R2016a) and a free program (spatialtoolbox). The hot spot analysis served to check whether an observed pattern is statistically significant. The method was evaluated with an established k-mean analysis. To test the hot-spot analysis by comparison, data from a grapefruit plantation (Adana, Turkey) was collected, including soil ECa, trunk circumference, and yield data. The plantation had 179 trees on a soil of type Xerofkuvent with clay and clay-loamy texture. The examination of the interaction between the critical values from the soil and plant water status information and the vegetative and generative plant growth variables was performed with the application from ANOVA.
The study indicates that the variability of the soil and plant information in fruit production is high, even considering small orchards. It was further indicated that the spatial patterns found in the soil ECa stayed constant through the years (r = 0.88 in 2011-2012 and r = 0.71 in 2012-2013). It was also demonstrated that CWSI determination may also be possible in semi-humid climate. A correlation (r = - 0.65, p < 0.0001) with the established method of leaf water potential analysis was found. The interaction between the ECa from various depths and the plant variables produced a highly significant connection with the topsoil in which the irrigation system was to be found. A correlation between yield and ECatopsoil of r = 0.52 was determined. By using the hot-spot analysis, extreme values in the spatial data could be determined. These extremes served to divide the zones (cold-spot, random, hot-spot). The random zone showed the highest correlation to the plant variables.
In summary it may be said that the cumulative water use efficiency (WUEc) was enhanced with high crop load. While the CWSI had no effect on fruit quality, the interaction of CWSI and WUEc even outweighed the impact of soil ECa on fruit quality in the production system with irrigation. In the plum orchard, irrigation was relevant for obtaining high quality produce even in the semi-humid climate.
The warm water geothermal reservoir below the village of Waiwera in New Zealand has been known by the native Maori for centuries. Development by the European immigrants began in 1863. Until the year 1969, the warm water flowing from all drilled wells was artesian. Due to overproduction, water up to 50 A degrees C now needs to be pumped to surface. Further, between 1975 and 1976, all warm water seeps on the beach of Waiwera ran dry. Within the context of sustainable water management, hydrogeological models must be developed as part of a management plan. Approaches of varying complexity have been set-up and applied since the 1980s. However, none of the models directly provide all results required for optimal water management. Answers are given simply to parts of the questions, nonetheless improving resource management of the geothermal reservoir.
The northern Tien Shan of Kyrgyzstan and Kazakhstan has been affected by a series of major earthquakes in the late 19th and early 20th centuries. To assess the significance of such a pulse of strain release in a continental interior, it is important to analyze and quantify strain release over multiple time scales. We have undertaken paleoseismological investigations at two geomorphically distinct sites (Panfilovkoe and Rot Front) near the Kyrgyz capital Bishkek. Although located near the historic epicenters, both sites were not affected by these earthquakes. Trenching was accompanied by dating stratigraphy and offset surfaces using luminescence, radiocarbon, and Be-10 terrestrial cosmogenic nuclide methods. At Rot Front, trenching of a small scarp did not reveal evidence for surface rupture during the last 5000 years. The scarp rather resembles an extensive debris-flow lobe. At Panfilovkoe, we estimate a Late Pleistocene minimum slip rate of 0.2 +/- 0.1 mm/a, averaged over at least two, probably three earthquake cycles. Dip-slip reverse motion along segmented, moderately steep faults resulted in hanging wall collapse scarps during different events. The most recent earthquake occurred around 3.6 +/- 1.3 kyr ago (1 sigma), with dip-slip offsets between 1.2 and 1.4 m. We calculate a probabilistic paleomagnitude to be between 6.7 and 7.2, which is in agreement with regional data from the Kyrgyz range. The morphotectonic signals in the northern Tien Shan are a prime example of deformation in a tectonically active intracontinental mountain belt and as such can help understand the longer-term coevolution of topography and seismogenic processes in similar structural settings worldwide.
Waterbodies such as lakes and ponds are abundant in vast Arctic landscapes and strongly affect the thermal state of the surrounding permafrost. In order to gain a better understanding of the impact of small-and medium-sized waterbodies on permafrost and the formation of thermokarst, a land surface model was developed that can represent the vertical and lateral thermal interactions between waterbodies and permafrost. The model was validated using temperature measurements from two typical waterbodies located within the Lena River delta in northern Siberia. Impact simulations were performed under current climate conditions as well as under a moderate and a strong climate-warming scenario. The performed simulations demonstrate that small waterbodies can rise the sediment surface temperature by more than 10 degrees C and accelerate permafrost thaw by a factor of between 4 and 5. Up to 70% of this additional heat flux into the ground was found to be dissipated into the surrounding permafrost by lateral ground heat flux in the case of small, shallow, and isolated waterbodies. Under moderate climate warming, the lateral heat flux was found to reduce permafrost degradation underneath waterbodies by a factor of 2. Under stronger climatic warming, however, the lateral heat flux was too small to prevent rapid permafrost degradation. The lateral heat flux was also found to strongly impede the formation of thermokarst. Despite this stabilizing effect, our simulations have demonstrated that underneath shallow waterbodies (<1 m), thermokarst initiation happens 30 to 40 years earlier than in simulations without preexisting waterbody.
Climate change increases riverine carbon outgassing, while export to the ocean remains uncertain
(2016)
Any regular interaction of land and river during flooding affects carbon pools within the terrestrial system, riverine carbon and carbon exported from the system. In the Amazon basin carbon fluxes are considerably influenced by annual flooding, during which terrigenous organic material is imported to the river. The Amazon basin therefore represents an excellent example of a tightly coupled terrestrial-riverine system. The processes of generation, conversion and transport of organic carbon in such a coupled terrigenous-riverine system strongly interact and are climate-sensitive, yet their functioning is rarely considered in Earth system models and their response to climate change is still largely unknown. To quantify regional and global carbon budgets and climate change effects on carbon pools and carbon fluxes, it is important to account for the coupling between the land, the river, the ocean and the atmosphere. We developed the RIVerine Carbon Model (RivCM), which is directly coupled to the well-established dynamic vegetation and hydrology model LPJmL, in order to account for this large-scale coupling. We evaluate RivCM with observational data and show that some of the values are reproduced quite well by the model, while we see large deviations for other variables. This is mainly caused by some simplifications we assumed. Our evaluation shows that it is possible to reproduce large-scale carbon transport across a river system but that this involves large uncertainties. Acknowledging these uncertainties, we estimate the potential changes in riverine carbon by applying RivCM for climate forcing from five climate models and three CO2 emission scenarios (Special Report on Emissions Scenarios, SRES). We find that climate change causes a doubling of riverine organic carbon in the southern and western basin while reducing it by 20% in the eastern and northern parts. In contrast, the amount of riverine inorganic carbon shows a 2- to 3-fold increase in the entire basin, independent of the SRES scenario. The export of carbon to the atmosphere increases as well, with an average of about 30 %. In contrast, changes in future export of organic carbon to the Atlantic Ocean depend on the SRES scenario and are projected to either decrease by about 8.9% (SRES A1B) or increase by about 9.1% (SRES A2). Such changes in the terrigenous-riverine system could have local and regional impacts on the carbon budget of the whole Amazon basin and parts of the Atlantic Ocean. Changes in riverine carbon could lead to a shift in the riverine nutrient supply and pH, while changes in the exported carbon to the ocean lead to changes in the supply of organic material that acts as a food source in the Atlantic. On larger scales the increased outgassing of CO2 could turn the Amazon basin from a sink of carbon to a considerable source. Therefore, we propose that the coupling of terrestrial and riverine carbon budgets should be included in subsequent analysis of the future regional carbon budget.
Hydrological and climatological controls on radiocarbon concentrations in a tropical stalagmite
(2016)
Precisely-dated stalagmites are increasingly important archives for the reconstruction of terrestrial paleoclimate at very high temporal resolution. In-depth understanding of local conditions at the cave site and of the processes driving stalagmite deposition is of paramount importance for interpreting proxy signals incorporated in stalagmite carbonate. Here we present a sub-decadally resolved dead carbon fraction (DCF) record for a stalagmite from Yok Balum Cave (southern Belize). The record is coupled to parallel stable carbon isotope (delta C-13) and U/Ca measurements, as well as radiocarbon (C-14) measurements from soils overlying the cave system. Using a karst carbon cycle model we disentangle the importance of soil and karst processes on stalagmite DCF incorporation, revealing a dominant host rock dissolution control on total DCF. Covariation between DCF, delta C-13, and U/Ca indicates that karst processes are a common driver of all three parameters, suggesting possible use of delta C-13 and trace element ratios to independently quantify DCF variability. A statistically significant multi-decadal lag of variable length exists between DCF and reconstructed solar activity, suggesting that solar activity influenced regional precipitation in Mesoamerica over the past 1500 years, but that the relationship was non-static. Although the precise nature of the observed lag is unclear, solar-induced changes in North Atlantic oceanic and atmospheric dynamics may play a role. (C) 2016 Elsevier Ltd. All rights reserved.
Widespread landscape changes are presently observed in the Arctic and are most likely to
accelerate in the future, in particular in permafrost regions which are sensitive to climate warming. To assess current and future developments, it is crucial to understand past
environmental dynamics in these landscapes. Causes and interactions of environmental variability can hardly be resolved by instrumental records covering modern time scales. However, long-term
environmental variability is recorded in paleoenvironmental archives. Lake sediments are important archives that allow reconstruction of local limnogeological processes as well as past environmental changes driven directly or indirectly by climate dynamics. This study aims at
reconstructing Late Quaternary permafrost and thermokarst dynamics in central-eastern Beringia,
the terrestrial land mass connecting Eurasia and North America during glacial sea-level low stands. In order to investigate development, processes and influence of thermokarst dynamics, several sediment cores from extant lakes and drained lake basins were analyzed to answer the
following research questions:
1. When did permafrost degradation and thermokarst lake development take place and what were enhancing and inhibiting environmental factors?
2. What are the dominant processes during thermokarst lake development and how are
they reflected in proxy records?
3. How did, and still do, thermokarst dynamics contribute to the inventory and properties of organic matter in sediments and the carbon cycle?
Methods applied in this study are based upon a multi-proxy approach combining
sedimentological, geochemical, geochronological, and micropaleontological analyses, as well as
analyses of stable isotopes and hydrochemistry of pore-water and ice. Modern field observations of water quality and basin morphometrics complete the environmental investigations.
The investigated sediment cores reveal permafrost degradation and thermokarst dynamics on different time scales. The analysis of a sediment core from GG basin on the northern Seward
Peninsula (Alaska) shows prevalent terrestrial accumulation of yedoma throughout the Early to
Mid Wisconsin with intermediate wet conditions at around 44.5 to 41.5 ka BP. This first wetland
development was terminated by the accumulation of a 1-meter-thick airfall tephra most likely originating from the South Killeak Maar eruption at 42 ka BP. A depositional hiatus between 22.5 and 0.23 ka BP may indicate thermokarst lake formation in the surrounding of the site which forms a yedoma upland till today. The thermokarst lake forming GG basin initiated 230 ± 30 cal a
BP and drained in Spring 2005 AD. Four years after drainage the lake talik was still unfrozen below 268 cm depth.
A permafrost core from Mama Rhonda basin on the northern Seward Peninsula preserved a
full lacustrine record including several lake phases. The first lake generation developed at 11.8 cal ka BP during the Lateglacial-Early Holocene transition; its old basin (Grandma Rhonda) is still partially preserved at the southern margin of the study basin. Around 9.0 cal ka BP a shallow and more dynamic thermokarst lake developed with actively eroding shorelines and potentially intermediate shallow water or wetland phases (Mama Rhonda). Mama Rhonda lake drainage at 1.1 cal ka BP was followed by gradual accumulation of terrestrial peat and top-down refreezing of the lake talik. A significant lower organic carbon content was measured in Grandma Rhonda deposits (mean TOC of 2.5 wt%) than in Mama Rhonda deposits (mean TOC of 7.9 wt%) highlighting the impact of thermokarst dynamics on biogeochemical cycling in different lake generations by thawing and mobilization of organic carbon into the lake system.
Proximal and distal sediment cores from Peatball Lake on the Arctic Coastal Plain of Alaska revealed young thermokarst dynamics since about 1,400 years along a depositional gradient based on reconstructions from shoreline expansion rates and absolute dating results. After its initiation as a remnant pond of a previous drained lake basin, a rapidly deepening lake with increasing oxygenation of the water column is evident from laminated sediments, and higher Fe/Ti and Fe/S ratios in the sediment. The sediment record archived characterizing shifts in depositional regimes and sediment sources from upland deposits and re-deposited sediments from drained thaw lake basins depending on the gradually changing shoreline configuration. These changes are evident from alternating organic inputs into the lake system which highlights the potential for thermokarst lakes to recycle old carbon from degrading permafrost deposits of its catchment.
The lake sediment record from Herschel Island in the Yukon (Canada) covers the full Holocene period. After its initiation as a thermokarst lake at 11.7 cal ka BP and intense thermokarst activity until 10.0 cal ka BP, the steady sedimentation was interrupted by a depositional hiatus at 1.6 cal ka BP which likely resulted from lake drainage or allochthonous slumping due to collapsing shore lines. The specific setting of the lake on a push moraine composed of marine deposits is reflected in the sedimentary record. Freshening of the maturing lake is indicated by decreasing electrical conductivity in pore-water. Alternation of marine to freshwater ostracods and foraminifera confirms decreasing salinity as well but also reflects episodical re-deposition of allochthonous marine sediments.
Based on permafrost and lacustrine sediment records, this thesis shows examples of the Late Quaternary evolution of typical Arctic permafrost landscapes in central-eastern Beringia and the complex interaction of local disturbance processes, regional environmental dynamics and global climate patterns. This study confirms that thermokarst lakes are important agents of organic matter recycling in complex and continuously changing landscapes.
Permafrost-related processes drive regional landscape dynamics in the Arctic terrestrial system. A better understanding of past periods indicative of permafrost degradation and aggradation is important for predicting the future response of Arctic landscapes to climate change. Here, we used a multi-proxy approach to analyse a4m long sediment core from a drained thermokarst lake basin on the northern Seward Peninsula in western Arctic Alaska (USA). Sedimentological, biogeochemical, geochronological, micropalaeontological (ostracoda, testate amoebae) and tephra analyses were used to determine the long-term environmental Early-Wisconsin to Holocene history preserved in our core for central Beringia. Yedoma accumulation dominated throughout the Early to Late-Wisconsin but was interrupted by wetland formation from 44.5 to 41.5ka BP. The latter was terminated by the deposition of 1m of volcanic tephra, most likely originating from the South Killeak Maar eruption at about 42ka BP. Yedoma deposition continued until 22.5ka BP and was followed by a depositional hiatus in the sediment core between 22.5 and 0.23ka BP. We interpret this hiatus as due to intense thermokarst activity in the areas surrounding the site, which served as a sediment source during the Late-Wisconsin to Holocene climate transition. The lake forming the modern basin on the upland initiated around 0.23ka BP and drained catastrophically in spring 2005. The present study emphasises that Arctic lake systems and periglacial landscapes are highly dynamic and that permafrost formation as well as degradation in central Beringia was controlled by regional to global climate patterns as well as by local disturbances. Copyright (c) 2015 John Wiley & Sons, Ltd.
Permafrost degradation influences the morphology, biogeochemical cycling and hydrology of Arctic landscapes over a range of time scales. To reconstruct temporal patterns of early to late Holocene permafrost and thermokarst dynamics, site-specific palaeo-records are needed. Here we present a multi-proxy study of a 350-cm-long permafrost core from a drained lake basin on the northern Seward Peninsula, Alaska, revealing Lateglacial toHolocene thermokarst lake dynamics in a central location of Beringia. Use of radiocarbon dating, micropalaeontology (ostracods and testaceans), sedimentology (grain-size analyses, magnetic susceptibility, tephra analyses), geochemistry (total nitrogen and carbon, total organic carbon, C-13(org)) and stable water isotopes (O-18, D, dexcess) of ground ice allowed the reconstruction of several distinct thermokarst lake phases. These include a pre-lacustrine environment at the base of the core characterized by the Devil Mountain Maar tephra (22800 +/- 280cal. a BP, Unit A), which has vertically subsided in places due to subsequent development of a deep thermokarst lake that initiated around 11800cal. a BP (Unit B). At about 9000cal. a BP this lake transitioned from a stable depositional environment to a very dynamic lake system (Unit C) characterized by fluctuating lake levels, potentially intermediate wetland development, and expansion and erosion of shore deposits. Complete drainage of this lake occurred at 1060cal. a BP, including post-drainage sediment freezing from the top down to 154cm and gradual accumulation of terrestrial peat (Unit D), as well as uniform upward talik refreezing. This core-based reconstruction of multiple thermokarst lake generations since 11800cal. a BP improves our understanding of the temporal scales of thermokarst lake development from initiation to drainage, demonstrates complex landscape evolution in the ice-rich permafrost regions of Central Beringia during the Lateglacial and Holocene, and enhances our understanding of biogeochemical cycles in thermokarst-affected regions of the Arctic.
Abrupt climate changes and fluctuations over short time scales are superimposed on long-term climate changes. Understanding rapid climate fluctuations at the decadal time scale over the past millennium will enhance our understanding of patterns of climate variability and aid in forecasting climate changes in the future. In this study, climate changes on the southeastern Tibetan Plateau over the past millennium were determined from a 4.82-m-long sediment core from Basomtso Lake. At the centennial time scale, the Medieval Climate Anomaly (MCA), Little Ice Age (LIA) and Current Warm Period (CWP) are distinct in the Basomtso region. Rapid climate fluctuations inferred from five episodes with higher sediment input and likely warmer conditions, as well as seven episodes with lower sediment input and likely colder conditions, were well preserved in our record. These episodes with higher and lower sediment input are characterized by abrupt climate changes and short time durations. Spectral analysis indicates that the climate variations at the centennial scale on the southeastern Tibetan Plateau are influenced by solar activity during the past millennium.
The onset of modern central Asian atmospheric circulation is traditionally linked to the interplay of surface uplift of the Mongolian and Tibetan-Himalayan orogens, retreat of the Paratethys sea from central Asia and Cenozoic global cooling. Although the role of these players has not yet been unravelled, the vast dust deposits of central China support the presence of arid conditions and modern atmospheric pathways for the last 25 million years (Myr). Here, we present provenance data from older (42-33 Myr) dust deposits, at a time when the Tibetan Plateau was less developed, the Paratethys sea still present in central Asia and atmospheric pCO(2) much higher. Our results show that dust sources and near-surface atmospheric circulation have changed little since at least 42 Myr. Our findings indicate that the locus of central Asian high pressures and concurrent aridity is a resilient feature only modulated by mountain building, global cooling and sea retreat.
The loess and paleosol sequences of the Chinese Loess Plateau are composed of Quaternary dust, the origin of which has been the subject of considerable debate. Some recent U-Pb geochronological studies of eolian zircons have proposed the existence of two major wind pathways: from the north and northwest, through the Badan Jaran, Tengger, and Mu Us Deserts during interglacials, and from the west, through the Qaidam Basin during glacials. Others have emphasized the importance of Yellow River sediment supply in the Chinese Loess Plateau sediment budget. However, tracking dust source regions through U-Pb dating in a statistically robust manner is particularly complex given the similar age peaks in the age probability distributions of potential source regions in Central Asia. This paper presents 2410 new U-Pb ages of detrital zircons from wind-eroded strata, Quaternary eolian deposits, and modern river sands in central China in order to increase the robustness and the spatial resolution of zircon age distributions in dust source regions. We then propose a new mixture modeling technique to statistically address the contribution of these different sources to the Chinese Loess Plateau sedimentary budget. Our contribution estimates indicate that eolian supply to the Chinese Loess Plateau is dominated (60%-70%) by reworking of Yellow River sediment. Moreover, evidence of Qaidam Basin-sourced zircons (15%-20%) in both loess (glacial) and paleosol (interglacial) layers corroborates the existence of an erosive wind pathway through the Qaidam Basin during glacials and implies that a substantial portion of the interglacial dust is recycled from older glacial loess. We propose that sediment reworking of Yellow River sediment and older loess deposits by wind on the Chinese Loess Plateau homogenized the eolian zircon populations toward a glacial provenance due to higher (2-20 times) dust accumulation rates during glacials. These findings suggest that the Chinese Loess Plateau has evolved as a more dynamic landform than previous thought, where wind deflation, fluvial input, lateral transport, and accumulation of sediment are equally important. These internal reworking effects would then significantly bias the paleoclimatic interpretations based on eolian dust properties of the Chinese Loess Plateau.
Monitoring the phase space of ecosystems: Concept and examples from the Quillow catchment, Uckermark
(2016)
Ecosystem research benefits enormously from the fact that comprehensive data sets of high quality, and covering long time periods are now increasingly more available. However, facing apparently complex interdependencies between numerous ecosystem components, there is urgent need rethinking our approaches in ecosystem research and applying new tools of data analysis. The concept presented in this paper is based on two pillars. Firstly, it postulates that ecosystems are multiple feedback systems and thus are highly constrained. Consequently, the effective dimensionality of multivariate ecosystem data sets is expected to be rather low compared to the number of observables. Secondly, it assumes that ecosystems are characterized by continuity in time and space as well as between entities which are often treated as distinct units. Implementing this concept in ecosystem research requires new tools for analysing large multivariate data sets. This study presents some of them, which were applied to a comprehensive water quality data set from a long-term monitoring program in Northeast Germany in the Uckermark region, one of the LTER-D (Long Term Ecological Research network, Germany) sites. Short-term variability of the kettle hole water samples differed substantially from that of the stream water samples, suggesting different processes generating the dynamics in these two types of water bodies. However, again, this seemed to be due to differing intensities of single processes rather than to completely different processes. We feel that research aiming at elucidating apparently complex interactions in ecosystems could make much more efficient use from now available large monitoring data sets by implementing the suggested concept and using corresponding innovative tools of system analysis. (C) 2015 Elsevier Ltd. All rights reserved.
Mediterranean climate is characterized by highly irregular rainfall patterns with marked differences between wet and dry seasons which lead to highly variable hydrological fluvial regimes. As a result, and in order to ensure water availability and reduce its temporal variability, a high number of large dams were built during the 20th century (more than 3500 located in Mediterranean rivers). Dams modify the flow regime but also interrupt the continuity of sediment transfer along the river network, thereby changing its functioning as an ecosystem. Within this context, the present paper aims to assess the suspended sediment loads and dynamics of two climatically contrasting Mediterranean regulated rivers (i.e. the Esera and Siurana) during a 2-yr period. Key findings indicate that floods were responsible for 92% of the total suspended sediment load in the River Siurana, while this percentage falls to 70% for the Esera, indicating the importance of baseflows on sediment transport in this river. This fact is related to the high sediment availability, with the Esera acting as a non-supply-limited catchment due to the high productivity of the sources (i.e. badlands). In contrast, the Siurana can be considered a supply-limited system due to its low geomorphic activity and reduced sediment availability, with suspended sediment concentration remaining low even for high magnitude flood events. Reservoirs in both rivers reduce sediment load up to 90%, although total runoff is only reduced in the case of the River Esera. A remarkable fact is the change of the hydrological character of the River Lem downstream for the dam, shifting from a humid mountainous river regime to a quasi-invariable pattern, whereas the Siurana experiences the opposite effect, changing from a flashy Mediterranean river to a more constant flow regime below the dam. (C) 2015 Elsevier B.V. All rights reserved.
The Earth’s shallow subsurface with sedimentary cover acts as a waveguide to any incoming wavefield. Within the framework of my thesis, I focused on the characterization of this shallow subsurface within tens to few hundreds of meters of sediment cover. I imaged the seismic 1D shear wave velocity (and possibly the 1D compressional wave velocity). This information is not only required for any seismic risk assessment, geotechnical engineering or microzonation activities, but also for exploration and global seismology where site effects are often neglected in seismic waveform modeling.
First, the conventional frequency-wavenumber (f - k) technique is used to derive the dispersion characteristic of the propagating surface waves recorded using distinct arrays of seismometers in 1D and 2D configurations. Further, the cross-correlation technique is applied to seismic array data to estimate the Green’s function between receivers pairs combination assuming one is the source and the other the receiver. With the consideration of a 1D media, the estimated cross-correlation Green’s functions are sorted with interstation distance in a virtual 1D active seismic experiment. The f - k technique is then used to estimate the dispersion curves. This integrated analysis is important for the interpretation of a large bandwidth of the phase velocity dispersion curves and therefore improving the resolution of the estimated 1D Vs profile.
Second, the new theoretical approach based on the Diffuse Field Assumption (DFA) is used for the interpretation of the observed microtremors H/V spectral ratio. The theory is further extended in this research work to include not only the interpretation of the H/V measured at the surface, but also the H/V measured at depths and in marine environments. A modeling and inversion of synthetic H/V spectral ratio curves on simple predefined geological structures shows an almost perfect recovery of the model parameters (mainly Vs and to a lesser extent Vp). These results are obtained after information from a receiver at depth has been considered in the inversion.
Finally, the Rayleigh wave phase velocity information, estimated from array data, and the H/V(z, f) spectral ratio, estimated from a single station data, are combined and inverted for the velocity profile information. Obtained results indicate an improved depth resolution in comparison to estimations using the phase velocity dispersion curves only. The overall estimated sediment thickness is comparable to estimations obtained by inverting the full micortremor H/V spectral ratio.
We present an integrated approach for deriving the 1D shear wave velocity (Vs) information at few tens to hundreds of meters down to the first strong impedance contrast in typical sedimentary environments. We use multiple small aperture seismic arrays in 1D and 2D configuration to record active and passive seismic surface wave data at two selected geotechnical sites in Germany (Horstwalde & Lobnitz). Standard methods for data processing include the Multichannel Analysis of Surface Waves (MASW) method that exploits the high frequency content in the active data and the sliding window frequency-wavenumber (f-k) as well as the spatial autocorrelation (SPAC) methods that exploit the low frequency content in passive seismic data. Applied individually, each of the passive methods might be influenced by any source directivity in the noise wavefield. The advantages of active shot data (known source location) and passive microtremor (low frequency content) recording may be combined using a correlation based approach applied to the passive data in the so called Interferometric Multichannel Analysis of Surface Waves (IMASW). In this study, we apply those methods to jointly determine and interpret the dispersion characteristics of surface waves recorded at Horstwalde and Lobnitz. The reliability of the dispersion curves is controlled by applying strict limits on the interpretable range of wavelengths in the analysis and further avoiding potentially biased phase velocity estimates from the passive f-k method by comparing to those derived from the SPatial AutoCorrelation method (SPAC). From our investigation at these two sites, the joint analysis as proposed allows mode extraction in a wide frequency range (similar to 0.6-35 Hz at Horstwalde and similar to 1.5-25 Hz at Lobnitz) and consequently improves the Vs profile inversion. To obtain the shear wave velocity profiles, we make use of a global inversion approach based on the neighborhood algorithm to invert the interpreted branches of the dispersion curves. Within the uncertainty given by the apparent spread of forward models we find that besides a well defined sediment velocity range also a reasonable minimum estimate of bedrock depth and bedrock velocity can be achieved. The Vs estimate for the best model in Horstwalde ranges from similar to 190 m/s at the surface up to similar to 390 m/s in the bottom of the soft sediment column. The bedrock starts earliest around 200 m depth and bedrock velocities are higher than 1000 m/s. In Lobnitz, we observe slightly lower velocities for the sediments (similar to 165-375 m/s for the best model) and a minimum thickness of 75 m. (C) 2016 Elsevier B.V. All rights reserved.
We compute seismic velocity profiles by a combined inversion of surface-wave phase-velocity dispersion curves together with the full spectrum of the microtremor horizontal-to-vertical (H/V) spectral ratio at two sediment-covered sites in Germany. The sediment deposits are approximately 100 m thick at the first test site and approximately 400 m thick at the second test site. We have used an extended physical model based on the diffuse wavefield assumption for the interpretation of the observed microtremor H/V spectral ratio. The extension includes the interpretation of the microtremor H/V spectral ratio observed at depth (in boreholes). This full-wavefield approach accounts for the energy contribution from the body and surface waves, and thus it allows for inverting the properties of the shallow subsurface. We have obtained the multimode phase velocity dispersion curves from an independent study, and a description of the extracted branches and their interpretation was developed. The inversion results indicate that the combined approach using seismic ambient noise and actively generated surface-wave data will improve the accuracy of the reconstructed near-surface velocity model, a key step in microzonation, geotechnical engineering, seismic statics corrections, and reservoir imaging.
The reservoir sediments are important sinks for organic carbon (OC), the OC burial being dependent on two opposite processes, deposition and mineralization. Hence factors such as severe water level fluctuations are expected to influence the rate of OC accumulation as they may affect both deposition and mineralization. The Barasona Reservoir has been historically threatened by siltation, whilst the use of water for irrigation involves a drastic decrease of the water level. In this context, we have studied the physical and chemical characteristics (grain size, major and minor elemental compositions, organic and inorganic carbon, and nitrogen) of the recent sediments of the Barasona Reservoir and the relationships among them in order to: a) elucidate the main processes governing OC accumulation, b) evaluate the rate of OC mineralization and c) approach the effect of drought on the sediment characteristics in this system. Our results indicated that Barasona sediments were dominated by fine silts (>60%) and clays (>20%), the mean particle size decreasing from tail to dam. Desiccation increased particle sorting and size distribution became bimodal, but no effect on average size was observed. Attending to the composition, Barasona sediments were very homogeneous with low concentrations of nitrogen (TN) and phosphorus (<12 g kg(-1) dw and <0.6 g kg(-1) dw, respectively) and high concentration of OC (approximate to 36 kg(-1) dw). TN was negatively related to dry weight Sediment mixing due to drastic changes in water level may have favoured the observed homogeneity of Barasona sediments affecting carbon, major ions and grain size. The high amount of OC deposited in Barasona sediment suggested that the adsorption of OC onto fine particles was more important than in boreal lakes. The rate of oxygen consumption by wet sediment ranged from 2.26 to 3.15 mg O-2 m(-2) day(-1), values close to those compiled for Mediterranean running waters. (C) 2015 Elsevier B.V. All rights reserved.
Subsurface water flows play a key role in the distribution of water and solutes and thereby in the water availability for plants. However, the characterization of different flow processes (i.e.matrix and preferential flow), the frequency and factors that cause them, is relatively rare. This characterization enables a better understanding of spatio-temporal variability of water resources and allows for the design of models to be improved. Using a method based on the time derivative of soil moisture variation known as maximum wetting slope, types of soil wetting processes were classified and quantified. For this, capacitance sensors, which registered the volumetric water content at high temporal resolution (30min) for more than two hydrological years, were installed at different depths and placed in soil moisture stations with different vegetation covers, lithology and topographic position. Results indicated that there is a general behaviour or pattern of soil moisture dynamics in the catchment with a dominant occurrence of slower soil wetting processes (>50%), caused by matrix flows, and a low occurrence of those faster processes (<30%), originated by preferential flows. Nevertheless, when the total volume of water is considered, preferential flow becomes the dominant process, so that the ecological role of both flow types becomes prominent in water-limited environments. Statistical multivariate analyses based on data-mining techniques proved that although both flow types depend on variables associated with precipitation and antecedent soil moisture conditions, faster soil wetting processes are mainly related to variables such as rainfall intensity and topography, while slower soil wetting processes are related to flow velocity, soils or vegetation. Copyright (c) 2015 John Wiley & Sons, Ltd.
Regulation alters the characteristics of riversty transforming parts of them into lakes, affecting their hydrology and also the physical, chemical, and biological characteristics and dynamics. Reservoirs have proven to be very effective retaining particulate materials, thereby avoiding the downstream transport of suspended sediment and the chemical substances associated with it (e.g. Carbon, C or Nitrogen, N). The study of fluvial transport of C and N is of great interest since river load represents a major link to the global C and N cycles. Moreover, reservoirs are the most important sinks for organic carbon among inland waters and have a potential significance as nitrogen sinks. In this respect, this paper investigates the effects of a Pyrenean reservoir on the runoff, suspended sediment, C and N derived from the highly active Esera and Isabena rivets. Key findings indicate that the reservoir causes a considerable impact on the Esera-Isabena river fluxes, reducing them dramatically as almost all the inputs are retained within the reservoir. Despite the very dry study year (2011-2012), it can be calculated that almost 300,000 t of suspended sediment were deposited into the Barasona Reservoir, from which more than 16,000 were C (i.e. 2200 t as organic C) and 222 t were N. These values may not be seen as remarkable in a wider global context but, assuming that around 30 hm(3) of sediment are currently stored in the reservoir, figures would increase up to ca. 2.6 x 10(6) t of C (i.e. 360,000 t of organic C) and 35,000 t of N. Nevertheless, these values are indicative and should be treated with caution as there is incomplete understanding of all the processes which affect C and N. Further investigation to establish a more complete picture of C and N yields and budgets by monitoring the different processes involved is essential. (C) 2015 Elsevier B.V. All rights reserved.
A broad array of new provenance and stable isotope data are presented from two magnetostratigraphically dated sections in the south-eastern Issyk Kul basin of the Central Kyrgyz Tien Shan. The results presented here are discussed and interpreted for two plausible magnetostratigraphic age models. A combination of zircon U-Pb provenance, paleocurrent and conglomerate clast count analyses is used to determine sediment provenance. This analysis reveals that the first coarse-grained, syntectonic sediments (Dzhety Oguz formation) were sourced from the nearby Terskey Range, supporting previous thermochronology-based estimates of a ca. 25-20 Ma onset of deformation in the range. Climate variations are inferred using carbonate stable isotope (delta O-18 and delta C-13) data from 53 samples collected in the two sections and are compared with the oxygen isotope compositions of modern water from 128 samples. Two key features are identified in the stable isotope data set derived from the sediments: (1) isotope values, in particular delta C-13, decrease between ca. 26.0 and 23.6 or 25.6 and 21.0 Ma, and (2) the scatter of delta O-18 values increased significantly after ca. 22.6 or 16.9 Ma. The first feature is interpreted to reflect progressively wetter conditions. Because this feature slightly post-dates the onset of deformation in the Terskey Range, we suggest that it has been caused by orographically enhanced precipitation, implying that surface uplift accompanied late Cenozoic deformation and rock uplift in the Terskey Range. The increased scatter could reflect variable moisture source or availability caused by global climate change following the onset of Miocene glaciations at ca. 22.6 Ma, or enhanced evaporation during the Mid-Miocene climatic optimum at ca. 17-15 Ma.
In this study, we provide a comprehensive analysis of trends in the extremes during the Indian summer monsoon (ISM) months (June to September) at different temporal and spatial scales. Our goal is to identify and quantify spatiotemporal patterns and trends that have emerged during the recent decades and may be associated with changing climatic conditions. Our analysis primarily relies on quantile regression that avoids making any subjective choices on spatial, temporal, or intensity pattern of extreme rainfall events. Our analysis divides the Indian monsoon region into climatic compartments that show different and partly opposing trends. These include strong trends toward intensified droughts in Northwest India, parts of Peninsular India, and Myanmar; in contrast, parts of Pakistan, Northwest Himalaya, and Central India show increased extreme daily rain intensity leading to higher flood vulnerability. Our analysis helps explain previously contradicting results of trends in average ISM rainfall.
Knowledge about the magnitude of localised flooding of riverine areas is crucial for appropriate land management and administration at regional and local levels. However, detection and delineation of localised flooding with remote sensing techniques are often hampered on floodplains by the presence of herbaceous vegetation. To address this problem, this study presents the application of full waveform airborne laser scanning (ALS) data for detection of floodwater extent. In general, water surfaces are characterised by low values of backscattered energy due to water absorption of the infrared laser shots, but the exact strength of the recorded laser pulse depends on the area covered by the targets located within a laser pulse footprint area. To account for this we analysed the physical quantity of radio metrically calibrated ALS data, the backscattering coefficient, in relation to water and vegetation coverage within a single laser footprint. The results showed that the backscatter was negatively correlated to water coverage, and that of the three distinguished classes of water coverage (low, medium, and high) only the class with the largest extent of water cover (>70%) had relatively distinct characteristics that can be used for classification of water surfaces. Following the laser footprint analysis, three classifiers, namely AdaBoost with Decision Tree, Naive Bayes and Random Forest, were utilised to classify laser points into flooded and non-flooded classes and to derive the map of flooding extent. The performance of the classifiers is highly dependent on the set of laser points features used. Best performance was achieved by combining radiometric and geometric laser point features. The accuracy of flooding maps based solely on radiometric features resulted in overall accuracies of up to 70% and was limited due to the overlap of the backscattering coefficient values between water and other land cover classes. Our point-based classification methods assure a high mapping accuracy (similar to 89%) and demonstrate the potential of using full-waveform ALS data to detect water surfaces on floodplain areas with limited water surface exposition through the vegetation canopy. (C) 2016 International Society for Photogrammetry and Remote Sensing, Inc. (ISPRS). Published by Elsevier B.V. All rights reserved.
Earthquakes deform Earth's surface, building long-lasting topographic features and contributing to landscape and mountain formation.
However, seismic waves produced by earthquakes may also destabilize hillslopes, leading to large amounts of soil and bedrock moving downslope. Moreover, static deformation and shaking are suspected to damage the surface bedrock and therefore alter its future properties, affecting hydrological and erosional dynamics. Thus, earthquakes participate both in mountain building and stimulate directly or indirectly their erosion. Moreover, the impact of earthquakes on hillslopes has important implications for the amount of sediment and organic matter delivered to rivers, and ultimately to oceans, during episodic catastrophic seismic crises, the magnitude of life and property losses associated with landsliding, the perturbation and recovery of landscape properties after shaking, and the long term topographic evolution of mountain belts. Several of these aspects have been addressed recently through individual case studies but additional data compilation as well as theoretical or numerical modelling are required to tackle these issues in a more systematic and rigorous manner.
This dissertation combines data compilation of earthquake characteristics, landslide mapping, and seismological data interpretation with physically-based modeling in order to address how earthquakes impact on erosional processes and landscape evolution. Over short time scales (10-100 s) and intermediate length scales (10 km), I have attempted to improve our understanding and ability to predict the amount of landslide debris triggered by seismic shaking in epicentral areas. Over long time scales (1-100 ky) and across a mountain belt (100 km) I have modeled the competition between erosional unloading and building of topography associated with earthquakes. Finally, over intermediate time scales (1-10 y) and at the hillslope scale (0.1-1 km) I have collected geomorphological and seismological data that highlight persistent effects of earthquakes on landscape properties and behaviour.
First, I compiled a database on earthquakes that produced significant landsliding, including an estimate of the total landslide volume and area, and earthquake characteristics such as seismic moment and source depth. A key issue is the accurate conversion of landslide maps into volume estimates. Therefore I also estimated how amalgamation - when mapping errors lead to the bundling of multiple landslide into a single polygon - affects volume estimates from various earthquake-induced landslide inventories and developed an algorithm to automatically detect this artifact. The database was used to test a physically-based prediction of the total landslide area and volume caused by earthquakes, based on seismological scaling relationships and a statistical description of the landscape properties. The model outperforms empirical fits in accuracy, with 25 out of 40 cases well predicted, and allows interpretation of many outliers in physical terms. Apart from seismological complexities neglected by the model I found that exceptional rock strength properties or antecedent conditions may explain most outliers.
Second, I assessed the geomorphic effects of large earthquakes on landscape dynamics by surveying the temporal evolution of precipitation-normalized landslide rate. I found strongly elevated landslide rates following earthquakes that progressively recover over 1 to 4 years, indicating that regolith strength drops and recovers. The relaxation is clearly non-linear for at least one case, and does not seem to correlate with coseismic landslide reactivation, water table level increase or tree root-system recovery. I suggested that shallow bedrock is damaged by the earthquake and then heals on annual timescales. Such variations in ground strength must be translated into shallow subsurface seismic velocities that are increasingly surveyed with ambient seismic noise correlations. With seismic noise autocorrelation I computed the seismic velocity in the epicentral areas of three earthquakes where I constrained a change in landslide rate. We found similar recovery dynamics and timescales, suggesting that seismic noise correlation techniques could be further developed to meaningfully assess ground strength variations for landscape dynamics. These two measurements are also in good agreement with the temporal dynamics of post-seismic surface displacement measured by GPS. This correlation suggests that the surface healing mechanism may be driven by tectonic deformation, and that the surface regolith and fractured bedrock may behave as a granular media that slowly compacts as it is sheared or vibrated.
Last, I compared our model of earthquake-induced landsliding with a standard formulation of surface deformation caused by earthquakes to understand which parameters govern the competition between the building and destruction of topography caused by earthquakes. In contrast with previous studies I found that very large (Mw>8) earthquakes always increase the average topography, whereas only intermediate (Mw ~ 7) earthquakes in steep landscapes may reduce topography. Moreover, I illustrated how the net effect of earthquakes varies with depth or landscape steepness implying a complex and ambivalent role through the life of a mountain belt. Further I showed that faults producing a Gutenberg-Richter distribution of earthquake sizes, will limit topography over a larger range of fault sizes than faults producing repeated earthquakes with a characteristic size.
Large, compressional earthquakes cause surface uplift aswell as widespread mass wasting. Knowledge of their trade-off is fragmentary. Combining a seismologically consistent model of earthquake-triggered landsliding and an analytical solution of coseismic surface displacement, we assess how the mass balance of single earthquakes and earthquake sequences depends on fault size and other geophysical parameters. We find that intermediate size earthquakes (M-w 6-7.3) may cause more erosion than uplift, controlled primarily by seismic source depth and landscape steepness, and less so by fault dip and rake. Such earthquakes can limit topographic growth, but our model indicates that both smaller and larger earthquakes (M-w < 6, M-w > 7.3) systematically cause mountain building. Earthquake sequences with a Gutenberg-Richter distribution have a greater tendency to lead to predominant erosion, than repeating earthquakes of the same magnitude, unless a fault can produce earthquakes with M-w > 8 or more.
We present a new, seismologically consistent expression for the total area and volume of populations of earthquake-triggered landslides. This model builds on a set of scaling relationships between key parameters, such as landslide spatial density, seismic ground acceleration, fault length, earthquake source depth, and seismic moment. To assess the model we have assembled and normalized a catalog of landslide inventories for 40 shallow, continental earthquakes. Low landscape steepness causes systematic overprediction of the total area and volume of landslides. When this effect is accounted for, the model predicts the total landslide volume of 63% of 40 cases to within a factor 2 of the volume estimated from observations (R-2 = 0.76). The prediction of total landslide area is also sensitive to the landscape steepness, but less so than the total volume, and it appears to be sensitive to controls on the landslide size-frequency distribution, and possibly the shaking duration. Some outliers are likely associated with exceptionally strong rock mass in the epicentral area, while others may be related to seismic source complexities ignored by the model. However, the close match between prediction and estimate for about two thirds of cases in our database suggests that rock mass strength is similar in many cases and that our simple seismic model is often adequate, despite the variety of lithologies and tectonic settings covered. This makes our expression suitable for integration into landscape evolution models and application to the anticipation or rapid assessment of secondary hazards associated with earthquakes.
Successively smaller glacial extents have been proposed for continental Eurasia during the stadials of the last glacial period leading up to the Last Glacial Maximum (LGM). At the same time the large mountainous region east of Lake Baikal, Transbaikalia, has remained unexplored in terms of glacial chronology despite clear geomorphological evidence of substantial past glaciations. We have applied cosmogenic Be-10 exposure dating and optically stimulated luminescence to establish the first quantitative glacial chronology for this region. Based on eighteen exposure ages from five moraine complexes, we propose that large mountain ice fields existed in the Kodar and Udokan mountains during Oxygen Isotope Stage 2, commensurate with the global LGM. These ice fields fed valley glaciers (>100 km in length) reaching down to the Chara Depression between the Kodar and Udokan mountains and to the valley of the Vitim River northwest of the Kodar Mountains. Two of the investigated moraines date to the Late Glacial, but indications of incomplete exposure among some of the sampled boulders obscure the specific details of the post-LGM glacial history. In addition to the LGM ice fields in the highest mountains of Transbaikalia, we report geomorphological evidence of a much more extensive, ice-cap type glaciation at a time that is yet to be firmly resolved. (C) 2015 Elsevier Ltd. All rights reserved.
We analyze the spatiotemporal distribution of fault geometries from seismicity induced by fluid injection at The Geysers geothermal field. The consistency of these faults with the local stress field is investigated using (1) the fault instability coefficient I comparing the orientation of a fault with the optimal orientation for failure in the assumed stress field and (2) the misfit angle beta between slip vectors observed from focal mechanisms and predicted from stress tensor. A statistical approach is applied to calculate the most likely fault instabilities considering the uncertainties from focal mechanisms and stress inversion. We find that faults activated by fluid injection may display a broad range in orientations. About 72% of the analyzed seismicity occurs on faults with favorable orientation for failure with respect to the stress field. However, a number of events are observed either to occur on severely misoriented faults or to slip in a different orientation than predicted from stress field. These events mostly occur during periods of high injection rates and are located in proximity to the injection wells. From the stress inversion, the friction coefficient providing the largest overall instability is mu = 0.5. About 91% of the events are activated with an estimated excess pore pressure <10 MPa, in agreement with previous models considering the combined effect of thermal and poroelastic stress changes from fluid injection. Furthermore, high seismic activity and largest magnitudes occur on favorably oriented faults with large instability coefficients and low slip misfit angles.
Reconstructions of the deposition rate of windblown mineral dust in ocean sediments offer an important means of tracking past climate changes and of assessing the radiative and biogeochemical impacts of dust in past climates. Dust flux estimates in ocean sediments have commonly been based on the operationally defined lithogenic fraction of sediment samples. More recently, dust fluxes have been estimated from measurements of helium and thorium, as rare isotopes of these elements (He-3 and Th-230) allow estimates of sediment flux, and the dominant isotopes (He-4 and Th-232) are uniquely associated with the lithogenic fraction of marine sediments. In order to improve the fidelity of dust flux reconstructions based on He and Th, we present a survey of He and Th concentrations in sediments from dust source areas in East Asia, Australia and South America. Our data show systematic relationships between He and Th concentrations and grain size, with He concentrations decreasing and Th concentrations increasing with decreasing grain size. We find consistent He and Th concentrations in the fine fraction (<5 mu m) of samples from East Asia, Australia and Central South America (Puna-Central West Argentina), with Th concentrations averaging 14 mu g/g and He concentrations averaging 2 mu cc STP/g. We recommend use of these values for estimating dust fluxes in sediments where dust is dominantly fine-grained, and suggest that previous studies may have systematically overestimated Th-based dust fluxes by 30%. Source areas in Patagonia appear to have lower He and Th contents than other regions, as fine fraction concentrations average 0.8 mu cc STP/g and 9 mu g/g for He-4 and Th-232, respectively. The impact of grain size on lithogenic He and Th concentrations should be taken into account in sediments proximal to dust sources where dust grain size may vary considerably. Our data also have important implications for the hosts of He in long-traveled dust and for the He-3/He-4 ratio used for terrigenous He in studies of extraterrestrial He in sediments and ice. We also investigate the use of He/Th ratios as a provenance tracer. Our results suggest differences in fine fraction He/Th ratios between East Asia, Australia, central South America and Patagonia, with ratios showing a positive relationship with the geological age of source rocks. He/Th ratios may thus provide useful provenance information, for example allowing separation of Patagonian sources from Puna-Central West Argentina or Australian dust sources. He/Th ratios in open-ocean marine sediments are similar to ratios in the fine fraction of upwind dust source areas. He/Th ratios in mid-latitude South Atlantic sediments suggest that dust in this region primarily derives from the Puna-Central West Argentina region (23-32 degrees S) rather than Patagonia (>38 degrees S). In the equatorial Pacific, He/Th ratios are much lower than in extratropical Pacific sediments or potential source areas measured as a part of this study (East Asia, South America, Australia) for reasons that are at present unclear, complicating their use as provenance tracers in this region. (C) 2015 Elsevier Ltd. All rights reserved.
The current view regarding the timing of regionally developed penetrative tectonic fabrics in sedimentary rocks is that their development postdates lithification of those rocks. In this case, fabric development is achieved by a number of deformation mechanisms, including grain rigid body rotation, crystal-plastic deformation, and pressure solution. The latter is believed to be the primary mechanism responsible for the domainal structure of cleavage in low-grade metamorphic rocks. In this study we combine field observations with strain studies to characterize considerable (>50%) Acadian crustal shortening in a Devonian clastic sedimentary sequence from southwest Ireland. Despite these high levels of shortening there is a marked absence of the domainal cleavage structure and intraclast deformation that are expected with this level of deformation. Fabrics in these rocks are predominantly a product of rigid body rotation and repacking of extraformational clasts during deformation of a clastic sedimentary sequence before lithification was complete.
Surface movements during the largest subduction zone earthquakes commonly drown coastlines. Yet, on geological timescales, coastlines above subduction zones uplift. Here I use a morphometric analysis combined with a numerical model of landscape evolution to estimate uplift rates along the central Andean rasa-a low-relief coastal surface bounded by a steep cliff formed by wave erosion. I find that the rasa has experienced steady uplift of 0.13 +/- 0.04 mm per year along a stretch of more than 2,000 km in length, during the Quaternary. These long-term uplift rates do not correlate with Global Positioning System (GPS) measurements of interseismic movements over the decadal scale, which implies that permanent uplift is not predominantly accumulated during the interseismic period. Instead, the rate of rasa uplift correlates with slip during earthquakes straddling the crust-mantle transition, the Moho. Such deeper earthquakes with magnitude 7 to 8 that occurred between 1995 and 2012 resulted in decimetres of coastal uplift. Slip during these earthquakes is located below the locked portion of the plate interface, and therefore may translate into permanent deformation of the overlying plate, where it causes uplift of the coastline. Thus, lower parts of the plate boundary are stably segmented over hundreds to millions of years. I suggest the coastline marks the surface expression of the transition between the shallow, locked seismogenic domain and the deeper, conditionally stable domain where modest earthquakes build up topography.
The repeated occurrence of exceptional floods within a few years, such as the Rhine floods in 1993 and 1995 and the Elbe and Danube floods in 2002 and 2013, suggests that floods in Central Europe may be organized in flood-rich and flood-poor periods. This hypothesis is studied by testing the significance of temporal clustering in flood occurrence (peak-over-threshold) time series for 68 catchments across Germany for the period 1932-2005. To assess the robustness of the results, different methods are used: Firstly, the index of dispersion, which quantifies the departure from a homogeneous Poisson process, is investigated. Further, the time-variation of the flood occurrence rate is derived by non-parametric kernel implementation and the significance of clustering is evaluated via parametric and non-parametric tests. Although the methods give consistent overall results, the specific results differ considerably. Hence, we recommend applying different methods when investigating flood clustering. For flood estimation and risk management, it is of relevance to understand whether clustering changes with flood severity and time scale. To this end, clustering is assessed for different thresholds and time scales. It is found that the majority of catchments show temporal clustering at the 5% significance level for low thresholds and time scales of one to a few years. However, clustering decreases substantially with increasing threshold and time scale. We hypothesize that flood clustering in Germany is mainly caused by catchment memory effects along with intra- to inter-annual climate variability, and that decadal climate variability plays a minor role. (C) 2016 Elsevier B.V. All rights reserved.