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)
We present an algorithm that performs sequentially one-dimensional inversion of subsurface magnetic permeability and electrical conductivity by using multi-configuration electromagnetic induction sensor data. The presented method is based on the conversion of the in-phase and out-of-phase data into effective magnetic permeability and electrical conductivity of the equivalent homogeneous half-space. In the case of small-offset systems, such as portable electromagnetic induction sensors, for which in-phase and out-of-phase data are moderately coupled, the effective half-space magnetic permeability and electrical conductivity can be inverted sequentially within an iterative scheme. We test and evaluate the proposed inversion strategy using synthetic and field examples. First, we apply it to synthetic data for some highly magnetic environments. Then, the method is tested on real field data acquired in a basaltic environment to image a formation of archaeological interest. These examples demonstrate that a joint interpretation of in-phase and out-of-phase data leads to a better characterisation of the subsurface in magnetic environments such as volcanic areas.
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.
Here, we study the 3-D subduction initiation process induced by the interaction between a hot thermochemical mantle plume and oceanic lithosphere using thermo-mechanical viscoplastic finite difference marker-in-cell models. Our numerical modeling results show that self-sustaining subduction is induced by plume-lithosphere interaction when the plume is sufficiently buoyant, the oceanic lithosphere is sufficiently old and the plate is weak enough to allow the buoyant plume to. pass through it. Subduction initiation occurs following penetration of the lithosphere by the hot plume and the downward displacement of broken, nearly circular segments of lithosphere (proto-slabs) as a result of partially molten plume rocks overriding the proto-slabs. Our experiments show four different deformation regimes in response to plume-lithosphere interaction: a) self-sustaining subduction initiation, in which subduction becomes self-sustaining; b) frozen subduction initiation, in which subduction stops at shallow depths; c) slab break-off, in which the subducting circular slab breaks off soon after formation; and d) plume underplating, in which the plume does not pass through the lithosphere and instead spreads beneath it (i.e., failed subduction initiation). These regimes depend on several parameters, such as the size, composition, and temperature of the plume, the brittle/plastic strength and age of the oceanic lithosphere, and the presence/absence of lithospheric heterogeneities. The results show that subduction initiates and becomes self-sustaining when the lithosphere is older than 10 Myr and the non dimensional ratio of the plume buoyancy force and lithospheric strength above the plume is higher than approximately 2. The outcomes of our numerical experiments are applicable for subduction initiation in the modern and Precambrian Earth and for the origin of plume-related corona structures on Venus. (C) 2016 Elsevier B.V. All rights reserved.
Ice complex deposits are characteristic, ice-rich formations in northern East Siberia and represent an important part in the arctic carbon pool. Recently, these late Quaternary deposits are the objective of numerous investigations typically relying on outcrop and borehole data. Many of these studies can benefit from a 3D structural model of the subsurface for upscaling their observations or for constraining estimations of inventories, such as the local carbon stock. We have addressed this problem of structural imaging by 3D ground-penetrating radar (GPR), which, in permafrost studies, has been primarily used for 2D profiling. We have used a 3D kinematic GPR surveying strategy at a field site located in the New Siberian Archipelago on top of an ice complex. After applying a 3D GPR processing sequence, we were able to trace two horizons at depths below 20 m. Taking available borehole and outcrop data into account, we have interpreted these two features as interfaces of major lithologic units and derived a 3D cryostratigraphic model of the subsurface. Our data example demonstrated that a 3D surveying and processing strategy was crucial at our field site and showed the potential of 3D GPR to image geologic structures in complex ice-rich permafrost landscapes.
The importance of cultural ecosystem services in agricultural landscapes is increasingly recognized as agricultural scale enlargement and abandonment affect aesthetic and recreational values of agricultural landscapes. Landscape preference studies addressing these type of values often yield context-specific outcomes, limiting the applicability of their outcomes in landscape policy. Our approach measures the relative importance of landscape features across agricultural landscapes. This approach was applied in the agricultural landscapes of Winterswijk, The Netherlands (n=191) and the Markische Schweiz, Germany (n=113) among visitors in the agricultural landscape. We set up a parallel designed choice experiment, using regionally specific, photorealistic visualizations of four comparable landscape attributes. In the Dutch landscape visitors highly value hedgerows and tree lines, whereas groups of trees and crop diversity are highly valued in the German landscape. Furthermore, we find that differences in relative preference for landscape attributes are, to some extent, explained by socio-cultural background variables such as education level and affinity with agriculture of the visitors. This approach contributes to a better understanding of the cross-regional variation of aesthetic and recreational values and how these values relate to characteristics of the agricultural landscape, which could support the integration of cultural services in landscape policy. (C) 2015 Elsevier B.V. All rights reserved.
The molecular biomarker composition of two sediment cores from Sanabria Lake (NW Iberian Peninsula) and a survey of modern plants in the watershed provide a reconstruction of past vegetation and landscape dynamics since deglaciation. During a proglacial stage in Lake Sanabria (prior to 14.7 cal ka BP), very low biomarker concentration and carbon preference index (CPI) values similar to 1 suggest that the n-alkanes could have derived from eroded ancient sediment sources or older organic matter with high degree of maturity. During the Late glacial (14.7-11.7 cal ka BP) and the Holocene (last 11.7 cal ka BP) intervals with higher biomarker and triterpenoid concentrations (high %nC(29) , nC(31) alkanes), higher CPI and average carbon length (ACL), and lower P-aq (proportion of aquatic plants) are indicative of major contribution of vascular land plants from a more forested watershed (e.g. Mid Holocene period 7.0-4.0 cal ka BP). Lower biomarker concentrations (high %nC(27) alkanes), CPI and ACL values responded to short phases with decreased allochthonous contribution into the lake that correspond to centennial-scale periods of regional forest decline (e.g. 4-3 ka BP, Roman deforestation after 2.0 ka, and some phases of the LIA, seventeenth-nineteenth centuries). Human activities in the watershed were significant during early medieval times (1.3-1.0 cal ka BP) and since 1960 CE, in both cases associated with relatively higher productivity stages in the lake (lower biomarker and triterpenoid concentrations, high %nC(23) and %nC(31) respectively, lower ACL and CPI values and higher P-aq). The lipid composition of Sanabria Lake sediments indicates a major allochthonous (watershed-derived) contribution to the organic matter budget since deglaciation, and a dominant oligotrophic status during the lake history. The study constrains the climate and anthropogenic forcings and watershed versus lake sources in organic matter accumulation processes and helps to design conservation and management policies in mountain, oligotrophic lakes.
This manuscript proposes a method to assess hydrological drought in semi-arid environments under high impoundment rate and applies it to the semi-arid Jaguaribe River basin in Brazil. It analyzes droughts (1) in the largest reservoir systems; (2) in the Upper Basin, considering 4744 reservoirs, 800 wells and almost 18,000 cisterns; and (3) in reservoirs of different sizes during multiyear droughts. Results show that the water demand is constrained in the basin; hydrological and meteorological droughts are often out of phase; there is a negative correlation between storage level and drought severity; and the small systems cannot cope with long-term droughts.
A new view of Ecuador's complex geodynamics has been developed in the course of modeling seismic source zones for probabilistic seismic hazard analysis. This study focuses on two aspects of the plates' interaction at a continental scale: (a) age-related differences in rheology between Farallon and Nazca plates—marked by the Grijalva rifted margin and its inland projection—as they subduct underneath central Ecuador, and (b) the rapidly changing convergence obliquity resulting from the convex shape of the South American northwestern continental margin. Both conditions satisfactorily explain several characteristics of the observed seismicity and of the interseismic coupling. Intermediate-depth seismicity reveals a severe flexure in the Farallon slab as it dips and contorts at depth, originating the El Puyo seismic cluster. The two slabs position and geometry below continental Ecuador also correlate with surface expressions observable in the local and regional geology and tectonics. The interseismic coupling is weak and shallow south of the Grijalva rifted margin and increases northward, with a heterogeneous pattern locally associated to the Carnegie ridge subduction. High convergence obliquity is responsible for the North Andean Block northeastward movement along localized fault systems. The Cosanga and Pallatanga fault segments of the North Andean Block-South American boundary concentrate most of the seismic moment release in continental Ecuador. Other inner block faults located along the western border of the inter-Andean Depression also show a high rate of moderate-size earthquake production. Finally, a total of 19 seismic source zones were modeled in accordance with the proposed geodynamic and neotectonic scheme.
In the past, floods were basically managed by flood control mechanisms. The focus was set on the reduction of flood hazard. The potential consequences were of minor interest. Nowadays river flooding is increasingly seen from the risk perspective, including possible consequences. Moreover, the large-scale picture of flood risk became increasingly important for disaster management planning, national risk developments and the (re-) insurance industry. Therefore, it is widely accepted that risk-orientated flood management ap-proaches at the basin-scale are needed. However, large-scale flood risk assessment methods for areas of several 10,000 km² are still in early stages. Traditional flood risk assessments are performed reach wise, assuming constant probabilities for the entire reach or basin. This might be helpful on a local basis, but where large-scale patterns are important this approach is of limited use. Assuming a T-year flood (e.g. 100 years) for the entire river network is unrealistic and would lead to an overestimation of flood risk at the large scale. Due to the lack of damage data, additionally, the probability of peak discharge or rainfall is usually used as proxy for damage probability to derive flood risk. With a continuous and long term simulation of the entire flood risk chain, the spatial variability of probabilities could be consider and flood risk could be directly derived from damage data in a consistent way.
The objective of this study is the development and application of a full flood risk chain, appropriate for the large scale and based on long term and continuous simulation. The novel approach of ‘derived flood risk based on continuous simulations’ is introduced, where the synthetic discharge time series is used as input into flood impact models and flood risk is directly derived from the resulting synthetic damage time series.
The bottleneck at this scale is the hydrodynamic simu-lation. To find suitable hydrodynamic approaches for the large-scale a benchmark study with simplified 2D hydrodynamic models was performed. A raster-based approach with inertia formulation and a relatively high resolution of 100 m in combination with a fast 1D channel routing model was chosen.
To investigate the suitability of the continuous simulation of a full flood risk chain for the large scale, all model parts were integrated into a new framework, the Regional Flood Model (RFM). RFM consists of the hydrological model SWIM, a 1D hydrodynamic river network model, a 2D raster based inundation model and the flood loss model FELMOps+r. Subsequently, the model chain was applied to the Elbe catchment, one of the largest catchments in Germany. For the proof-of-concept, a continuous simulation was per-formed for the period of 1990-2003. Results were evaluated / validated as far as possible with available observed data in this period. Although each model part introduced its own uncertainties, results and runtime were generally found to be adequate for the purpose of continuous simulation at the large catchment scale.
Finally, RFM was applied to a meso-scale catchment in the east of Germany to firstly perform a flood risk assessment with the novel approach of ‘derived flood risk assessment based on continuous simulations’. Therefore, RFM was driven by long term synthetic meteorological input data generated by a weather generator. Thereby, a virtual time series of climate data of 100 x 100 years was generated and served as input to RFM providing subsequent 100 x 100 years of spatially consistent river discharge series, inundation patterns and damage values. On this basis, flood risk curves and expected annual damage could be derived directly from damage data, providing a large-scale picture of flood risk. In contrast to traditional flood risk analysis, where homogenous return periods are assumed for the entire basin, the presented approach provides a coherent large-scale picture of flood risk. The spatial variability of occurrence probability is respected. Additionally, data and methods are consistent. Catchment and floodplain processes are repre-sented in a holistic way. Antecedent catchment conditions are implicitly taken into account, as well as physical processes like storage effects, flood attenuation or channel–floodplain interactions and related damage influencing effects. Finally, the simulation of a virtual period of 100 x 100 years and consequently large data set on flood loss events enabled the calculation of flood risk directly from damage distributions. Problems associated with the transfer of probabilities in rainfall or peak runoff to probabilities in damage, as often used in traditional approaches, are bypassed.
RFM and the ‘derived flood risk approach based on continuous simulations’ has the potential to provide flood risk statements for national planning, re-insurance aspects or other questions where spatially consistent, large-scale assessments are required.
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.
Extra-tropical circulation systems impede poleward moisture advection by the Indian Summer Monsoon. In this context, the Himalayan range is believed to insulate the south Asian circulation from extra-tropical influences and to delineate the northern extent of the Indian Summer Monsoon in central Asia. Paleoclimatic evidence, however, suggests increased moisture availability in the Early Holocene north of the Himalayan range which is attributed to an intensification of the Indian Summer Monsoon. Nevertheless, mechanisms leading to a surpassing of the Himalayan range and the northern maximum extent of summer monsoonal influence remain unknown. Here we show that the Kunlun barrier on the northern Tibetan Plateau [similar to 36 degrees N] delimits Indian Summer Monsoon precipitation during the Holocene. The presence of the barrier relocates the insulation effect 1,000 km further north, allowing a continental low intensity branch of the Indian Summer Monsoon which is persistent throughout the Holocene. Precipitation intensities at its northern extent seem to be driven by differentiated solar heating of the Northern Hemisphere indicating dependency on energy-gradients rather than absolute radiation intensities. The identified spatial constraints of monsoonal precipitation will facilitate the prediction of future monsoonal precipitation patterns in Central Asia under varying climatic conditions.
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.
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.
We applied the geomorphic indices (hypsometry and stream length gradient) to evaluate the differential uplift of the central and southern Longmenshan, a mountain range characterized by rapid erosion, strong tectonic uplift, and devastating seismic hazards. The results of the geomorphic analysis indicate that the Beichuan-Yingxiu fault and the Shuangshi-Dachuan fault act as major tectonic boundaries separating areas experiencing rapid uplift from slow uplift. The results of the geomorphic analysis also suggest that the Beichuan-Yingxiu fault is the most active fault with the largest relative uplift rates compared to the rest of the faults in the Longmenshan fault system. We compared reflected relative uplift rates based on the hypsometry and stream length gradient indices with geological/geodetic absolute rates. Along-strike and across-strike variations in the hypsometry and stream length gradient correlate with the spatial patterns derived from the apatite fission track exhumation rates, the leveling-derived uplift rate, and coseismic vertical displacements during the 2008 Wenchuan earthquake. These data defined multiple fault relationships in a complex thrust zone and provided geomorphic evidence to evaluate the potential seismic hazards of the southern Longmenshan range.
Amides as thermo-sensitive tracers for investigating the thermal state of geothermal reservoirs
(2016)
The application of thermo-sensitive tracers is a promising technique for evaluating the thermal state of geothermal reservoirs. To extend the compound spectrum for hydrolyzable compounds to reservoir temperatures between 100 and 200 degrees C carboxamides were studied. The kinetic parameters of 17 self-synthesized amides were determined in hydrothermal batch and autoclave experiments. The influence of the molecular structure and the role of pH/pOH on hydrolysis kinetics were studied. Additionally, the thermal stabilities of the hydrolysis products were evaluated. The results demonstrate the high potential of tracers based on amide hydrolysis for use in medium enthalpy reservoirs. (C) 2016 Elsevier Ltd. All rights reserved.
Surveys for more than 9,500 households were conducted in the growing seasons 2002/2003 or 2003/2004 in eleven African countries: Burkina Faso, Cameroon, Ghana, Niger and Senegal in western Africa; Egypt in northern Africa; Ethiopia and Kenya in eastern Africa; South Africa, Zambia and Zimbabwe in southern Africa. Households were chosen randomly in districts that are representative for key agro-climatic zones and farming systems. The data set specifies farming systems characteristics that can help inform about the importance of each system for a country’s agricultural production and its ability to cope with short- and long-term climate changes or extreme weather events. Further it informs about the location of smallholders and vulnerable systems and permits benchmarking agricultural systems characteristics.
The resonance frequency of the transmission response in layered half-space model is important in the study of site effect because it is the frequency where the shake-ability of the ground is enhanced significantly. In practice, it is often determined by the H/V ratio technique in which the peak frequency of recorded H/V spectral ratio is interpreted as the resonance frequency. Despite of its importance, there has not been any formula of the resonance frequency of the layered half-space structure. In this paper, a simple approximate formula of the fundamental resonance frequency is presented after an exact formula in explicit form of the response function of vertically SH incident wave is obtained. The formula is in similar form with the one used in H/V ratio technique but it reflects several major effects of the model to the resonance frequency such as the arrangement of layers, the impedance contrast between layers and the half-space. Therefore, it could be considered as an improved formula used in H/V ratio technique. The formula also reflects the consistency between two approaches of the H/V ratio technique based on SH body waves or Rayleigh surface waves on the peak frequency under high impedance contrast condition. This formula is in explicit form and, therefore, may be used in the direct and inverse problem efficiently. A numerical illustration of the improved formula for an actual layered half-space model already investigated by H/V ratio technique is presented to demonstrate its new features and its improvement to the currently used formula.
An innovative approach to improve SRTM DEM using multispectral imagery and artificial neural network
(2016)
Although the Shuttle Radar Topography Mission [SRTM) data are a publicly accessible Digital Elevation Model [DEM) provided at no cost, its accuracy especially at forested area is known to be limited with root mean square error (RMSE) of approx. 14 m in Singapore's forested area. Such inaccuracy is attributed to the 5.6 cm wavelength used by SRTM that does not penetrate vegetation well. This paper considers forested areas of central catchment of Singapore as a proof of concept of an approach to improve the SRTM data set. The approach makes full use of (1) the introduction of multispectral imagery (Landsat 8), of 30 m resolution, into SRTM data; (2) the Artificial Neural Network (ANN) to flex its known strengths in pattern recognition and; (3) a reference DEM of high accuracy (1 m) derived through the integration of stereo imaging of worldview-1 and extensive ground survey points. The study shows a series of significant improvements of the SRTM when assessed with the reference DEM of 2 different areas, with RMSE reduction of ∼68% (from 13.9 m to 4.4 m) and ∼52% (from 14.2 m to 6.7 m). In addition, the assessment of the resulting DEM also includes comparisons with simple denoising methodology (Low Pass Filter) and commercially available product called NEXTMap® World 30™.
The analytical evaluation of diurnal temperature variation in riverbed sediments provides detailed information on exchange fluxes between rivers and groundwater. The underlying assumption of the stationary, one-dimensional vertical flow field is frequently violated in natural systems where subsurface water flow often has a significant horizontal component. In this paper, we present a new methodology for identifying the geometry of the subsurface flow field using vertical temperature profiles. The statistical analyses are based on model optimisation and selection and are used to evaluate the shape of vertical amplitude ratio profiles. The method was applied to multiple profiles measured around in-stream geomorphological structures in a losing reach of a gravel bed river. The predominant subsurface flow field was systematically categorised in purely vertical and horizontal (hyporheic, parafluvial) components. The results highlight that river groundwater exchange flux at the head, crest and tail of geomorphological structures significantly deviated from the one-dimensional vertical flow, due to a significant horizontal component. The geometry of the subsurface water flow depended on the position around the geomorphological structures and on the river level. The methodology presented in this paper features great potential for characterising the spatial patterns and temporal dynamics of complex subsurface flow geometries by using measured temperature time series in vertical profiles. (C) 2016 Elsevier B.V. All rights reserved.
Variational methods are employed in situations where exact Bayesian inference becomes intractable due to the difficulty in performing certain integrals. Typically, variational methods postulate a tractable posterior and formulate a lower bound on the desired integral to be approximated, e.g. marginal likelihood. The lower bound is then optimised with respect to its free parameters, the so-called variational parameters. However, this is not always possible as for certain integrals it is very challenging (or tedious) to come up with a suitable lower bound. Here, we propose a simple scheme that overcomes some of the awkward cases where the usual variational treatment becomes difficult. The scheme relies on a rewriting of the lower bound on the model log-likelihood. We demonstrate the proposed scheme on a number of synthetic and real examples, as well as on a real geophysical model for which the standard variational approaches are inapplicable.
This paper assesses the seasonality of the urban heat island (UHI) effect in the Greater London area (United Kingdom). Combining satellite-based observations and urban boundary layer climate modeling with the UrbClim model, the authors are able to address the seasonality of UHI intensity, on the basis of both land surface temperature (LST) and 2-m air temperature, for four individual times of the day (0130, 1030, 1330, and 2230 local time) and the daily means derived from them. An objective of this paper is to investigate whether the UHI intensities that are based on both quantities exhibit a similar hysteresis-like trajectory that is observed for LST when plotting the UHI intensity against the background temperature. The results show that the UrbClim model can satisfactorily reproduce both the observed urban rural LSTs and 2-m air temperatures as well as their differences and the hysteresis in the surface UHI. The hysteresis-like seasonality is largely absent in both the observed and modeled 2-m air temperatures, however. A sensitivity simulation of the UHI intensity to incoming solar radiation suggests that the hysteresis of the LST can mainly be attributed to the seasonal variation in incoming solar radiation.
Snowfall comprises a significant percentage of the annual water budget in High Mountain Asia (HMA), but snow water equivalent (SWE) is poorly constrained due to lack of in-situ measurements and complex terrain that limits the efficacy of modeling and observations. Over the past few decades, SWE has been estimated with passive microwave (PM) sensors with generally good results in wide, flat, terrain, and lower reliability in densely forested, complex, or high-elevation areas. In this study, we use raw swath data from five satellite - sensors the Special Sensor Microwave/Imager (SSMI) and Special Sensor Microwave Imager/Sounder (SSMIS) (1987-2015, F08, F11, F13, F17), Advanced Microwave Scanning Radiometer - Earth Observing System (AMSR-E, 2002-2011), AMSR2 (2012-2015), and the Global Precipitation Measurement (GPM, 2014-2015) - in order to understand the spatial and temporal structure of native sensor, topographic, and land cover biases in SWE estimates in HMA. We develop a thorough understanding of the uncertainties in our SWE estimates by examining the impacts of topographic parameters (aspect, relief, hillslope angle, and elevation), land cover, native sensor biases, and climate parameters (precipitation, temperature, and wind speed). HMA, with its high seasonality, large topographic gradients and low relief at high elevations provides an excellent context to examine a wide range of climatic, land-cover, and topographic settings to better constrain SWE uncertainties and potential sensor bias. Using a multi-parameter regression, we compare long-term SWE variability to forest fraction, maximal multiyear snow depth, topographic parameters, and long-term average wind speed across both individual sensor time series and a merged multi-sensor dataset. In regions where forest cover is extensive, it is the strongest control on SWE variability. In those regions where forest density is low (<5%), maximal snow depth dominates the uncertainty signal. In our regression across HMA, we find that forest fraction is the strongest control on SWE variability (75.8%), followed by maximal multi-year snow depth (7.82%), 90th percentile 10-m wind speed of a 10-year December-January-February (DJF) time series (5.64%), 25th percentile DJF 10-m wind speed (5.44%), and hillslope angle (5.24%). Elevation, relief, and terrain aspect show very low influence on SWE variability (<1%). We find that the GPM sensor provides the most robust regression results, and can be reliably used to estimate SWE in our study region. While forest cover and elevation have been integrated into many SWE algorithms, wind speed and long-term maximal snow depth have not. Our results show that wind redistribution of snow can have impacts on SWE, especially over large, flat, areas. Using our regression results, we have developed an understanding of sensor specific SWE uncertainties and their spatial patterns. The uncertainty maps developed in this study provide a first-order approximation of SWE-estimate reliability for much of HMA, and imply that high-fidelity SWE estimates can be produced for many high-elevation areas. (C) 2016 Elsevier Inc. All rights reserved.
Analysis of time-lapse ground-penetrating radar (GPR) data can provide information regarding subsurface hydrological processes, such as preferential flow. However, the analysis of time-lapse data is often limited by data quality; for example, for noisy input data, the interpretation of difference images is often difficult. Motivated by modern image-processing tools, we have developed two robust GPR attributes, which allow us to distinguish amplitude (contrast similarity) and time-shift (structural similarity) variations related to differences between individual time-lapse GPR data sets. We tested and evaluated our attributes using synthetic data of different complexity. Afterward, we applied them to a field data example, in which subsurface flow was induced by an artificial rainfall event. For all examples, we identified our structural similarity attribute to be a robust measure for highlighting time-lapse changes also in data with low signal-to-noise ratios. We determined that our new attribute-based workflow is a promising tool to analyze time-lapse GPR data, especially for imaging subsurface hydrological processes.
The water cycle of sites with shallow groundwater tables is characterized by complex interactions of hydrological and ecological processes. The water balance components, which are subject to diurnal fluctuations, are best measured with groundwater lysimeters. However, the lower boundary condition of such lysimeters affects most of the hydrological variables, particularly when considering short time scales, and has to be defined in such a way as to facilitate realistic simulations. In this paper, different means of controlling the lower boundary condition of groundwater lysimeters were compared with respect to their ability to simulate the behavior of the water balance components properly. Measurements of rain-free periods from a lysimeter station installed in the Spreewald wetland in north-east Germany were evaluated. The most common groundwater lysimeter type is controlled using a Mariotte bottle and sets the groundwater level in the soil monolith to a constant level, which here caused an alteration of the inflow to the lysimeter, with respect to both its value and diurnal behavior. Still, daily evapotranspiration values were realistic and this simple and robust approach may be used for time intervals not shorter than one day. High-resolution measurements can be gained from lysimeters that automatically adjust the groundwater level by a system of pumps and valves on an hourly basis. Still, reliable results were only obtained when the conditions in the lysimeter and the surrounding field, where the target groundwater level was measured, were in accordance. Otherwise (e.g., when the groundwater level differed) an unrealistic inflow behavior evolved. Reasonable results, even for slightly diverging conditions, were gained with a new approach that defined the lower boundary conditions by controlling the inflows and outflows of the lysimeter. This approach further enabled the groundwater level itself to be the study subject, thereby enlarging the field of possible applications of groundwater lysimeters. (C) 2015 Elsevier B.V. All rights reserved.
Rainfall-induced attenuation is a major source of underestimation for radar-based precipitation estimation at C-band. Unconstrained gate-by-gate correction procedures are known to be inherently unstable and thus not suited for unsupervised attenuation correction. In this study, we evaluate three different procedures to constrain gate-by-gate attenuation correction using reflectivity as the only input. These procedures are benchmarked against rainfall estimates from uncorrected radar data, using six years of radar observations from the single-polarized C-band radar in South-West Germany. The precipitation estimation error is obtained by comparing the radar-based estimates to rain gauge observations. All attenuation correction procedures benchmarked in this study lead to an effective improvement of precipitation estimation. The first method caps the corrections if the rain intensity increase exceeds a factor of two. The second method decreases the parameters of the attenuation correction iteratively for every radar beam calculation until attaining a stability criterion. The second method outperforms the first method and leads to a consistent distribution of path-integrated attenuation along the radar beam. As a third method, we propose a slight modification of Kraemer's approach which allows users to exert better control over attenuation correction by introducing an additional constraint that prevents unplausible corrections in cases of dramatic signal losses.
In March 2015, a new international blueprint for disaster risk reduction (DRR) was adopted in Sendai, Japan, at the end of the Third UN World Conference on Disaster Risk Reduction (WCDRR, 14-18 March 2015). We review and discuss the agreed commitments and targets, as well as the negotiation leading the Sendai Framework for DRR (SF-DRR) and discuss briefly its implication for the later UN-led negotiations on sustainable development goals and climate change.
Climate science today makes use of a variety of red globes to explore and communicate findings. These transform the iconography which informs this image: the idealised, even mythical vision of the blue, vulnerable and perfect marble is impaired by the application of the colours yellow and red. Since only predictions that employ a lot of red seem to exist, spectators are confronted with the message that the future Earth that might turn out as envisaged here is undesirable. Here intuitively powerful narrations of the end of the world may connect. By employing methods of art history and visual analysis, and building on examples from current Intergovernmental Panel on Climate Change reports and future scenario maps, this article explores how burning world images bear - intentionally or not - elements of horror and shock. My question explored here is as follows: should 'burning world' images be understood as a new and powerful cosmology?
The problem of estimating the maximum possible earthquake magnitude m(max) has attracted growing attention in recent years. Due to sparse data, the role of uncertainties becomes crucial. In this work, we determine the uncertainties related to the maximum magnitude in terms of confidence intervals. Using an earthquake catalog of Iran, m(max) is estimated for different predefined levels of confidence in six seismotectonic zones. Assuming the doubly truncated Gutenberg-Richter distribution as a statistical model for earthquake magnitudes, confidence intervals for the maximum possible magnitude of earthquakes are calculated in each zone. While the lower limit of the confidence interval is the magnitude of the maximum observed event, the upper limit is calculated from the catalog and the statistical model. For this aim, we use the original catalog which no declustering methods applied on as well as a declustered version of the catalog. Based on the study by Holschneider et al. (Bull Seismol Soc Am 101(4): 1649-1659, 2011), the confidence interval for m(max) is frequently unbounded, especially if high levels of confidence are required. In this case, no information is gained from the data. Therefore, we elaborate for which settings finite confidence levels are obtained. In this work, Iran is divided into six seismotectonic zones, namely Alborz, Azerbaijan, Zagros, Makran, Kopet Dagh, Central Iran. Although calculations of the confidence interval in Central Iran and Zagros seismotectonic zones are relatively acceptable for meaningful levels of confidence, results in Kopet Dagh, Alborz, Azerbaijan and Makran are not that much promising. The results indicate that estimating mmax from an earthquake catalog for reasonable levels of confidence alone is almost impossible.
To understand past flood changes in the Rhine catchment and in particular the role of anthropogenic climate change in extreme flows, an attribution study relying on a proper GCM (general circulation model) downscaling is needed. A downscaling based on conditioning a stochastic weather generator on weather patterns is a promising approach. This approach assumes a strong link between weather patterns and local climate, and sufficient GCM skill in reproducing weather pattern climatology. These presuppositions are unprecedentedly evaluated here using 111 years of daily climate data from 490 stations in the Rhine basin and comprehensively testing the number of classification parameters and GCM weather pattern characteristics. A classification based on a combination of mean sea level pressure, temperature, and humidity from the ERA20C reanalysis of atmospheric fields over central Europe with 40 weather types was found to be the most appropriate for stratifying six local climate variables. The corresponding skill is quite diverse though, ranging from good for radiation to poor for precipitation. Especially for the latter it was apparent that pressure fields alone cannot sufficiently stratify local variability. To test the skill of the latest generation of GCMs from the CMIP5 ensemble in reproducing the frequency, seasonality, and persistence of the derived weather patterns, output from 15 GCMs is evaluated. Most GCMs are able to capture these characteristics well, but some models showed consistent deviations in all three evaluation criteria and should be excluded from further attribution analysis.
Carbon and nutrient cycling in kettle hole sediments depending on hydrological dynamics: a review
(2016)
Kettle holes as a specific group of isolated, small lentic freshwater systems (LFS) often are (i) hot spots of biogeochemical cycling and (ii) exposed to frequent sediment desiccation and rewetting. Their ecological functioning is greatly determined by immanent carbon and nutrient transformations. The objective of this review is to elucidate effects of a changing hydrological regime (i.e., dry-wet cycles) on carbon and nutrient cycling in kettle hole sediments. Generally, dry-wet cycles have the potential to increase C and N losses as well as P availability. However, their duration and frequency are important controlling factors regarding direction and intensity of biogeochemical and microbiological responses. To evaluate drought impacts on sediment carbon and nutrient cycling in detail requires the context of the LFS hydrological history. For example, frequent drought events induce physiological adaptation of exposed microbial communities and thus flatten metabolic responses, whereas rare events provoke unbalanced, strong microbial responses. Different potential of microbial resilience to drought stress can irretrievably change microbial communities and functional guilds, gearing cascades of functional responses. Hence, dry-wet events can shift the biogeochemical cycling of organic matter and nutrients to a new equilibrium, thus affecting the dynamic balance between carbon burial and mineralization in kettle holes.
Carbonatites are peculiar magmatic rocks with mantle-related genesis, commonly interpreted as the products of melting of CO2-bearing peridotites, or resulting from the chemical evolution of mantle derived magmas, either through extreme "differentiation or secondary immiscibility. Here we report the first finding of anatectic carbonatites of crustal origin, preserved as calcite-rich polycrystalline inclusions in garnet from low-to-medium pressure migmatites of the Oberpfalz area, SW Bohemian Massif (Central Europe). These inclusions originally trapped a melt of calciocarbonatitic composition with a characteristic enrichment in Ba, Sr and LREE. This interpretation is supported by the results of a detailed microstructural and microchemical investigation, as well as re-melting experiments using a piston cylinder apparatus. Carbonatitic inclusions coexist in the same cluster with crystallized silicate melt inclusions (nanogranites) and COH fluid inclusions, suggesting conditions of primary immiscibility between two melts and a fluid during anatexis. The production of both carbonatitic and granitic melts during the same anatectic event requires a suitable heterogeneous protolith. This may be represented by a sedimentary sequence containing marble lenses of limited extension, similar to the one still visible in the adjacent central Moldanubian Zone. The presence of CO2-rich fluid inclusions suggests furthermore that high CO2 activity during anatexis may be required to stabilize a carbonate-rich melt in a silica-dominated, system. This natural occurrence displays a remarkable similarity with experiments on carbonate-silicate melt immiscibility, where CO2 saturation is a condition commonly imposed. In conclusion, this study shows how the investigation of partial melting through melt inclusion studies may unveil unexpected processes whose evidence, while preserved in stiff minerals such as garnet, is completely obliterated in the rest of the rock due to metamorphic re-equilibration. Our results thus provide invaluable new insights into the processes which shape the geochemical evolution of our planet, such as the redistribution of carbon and strategic metals during orogenesis. (C) 2016 Elsevier B.V. All rights reserved.
Carboniferous metagranites with U-Pb zircon crystallization ages of 331-315 Ma crop out in the Afyon zone in the northern margin of the Anatolide-Tauride Block, which is commonly regarded as part of Gondwana during the Late Palaeozoic. They are peraluminous, calc-alkaline and are characterized by increase in Rb and Ba, decrease in Nb-Ta, and enrichment in Sr and high LILE/HFSE ratios compatible with a continental arc setting. The metagranites intrude a metasedimentary sequence of phyllite, metaquartzite and marble; both the Carboniferous metagranites and metasedimentary rocks are overlain unconformably by Lower Triassic metaconglomerates, metavolcanics and Upper Triassic to Cretaceous recrystallized limestones. The low-grade metamorphism and deformation occurred at the Cretaceous-Tertiary boundary. There is no evidence for Carboniferous deformation and metamorphism in the region. Carboniferous arc-type granites and previously described Carboniferous subduction-accretion complexes on the northern margin of the Anatolide-Tauride Block suggest southward subduction of Paleotethys under Gondwana during the Carboniferous. Considering the Variscan-related arc granites in Pelagonian and Sakarya zones on the active southern margin of Laurasia, a dual subduction of Paleotethys can be envisaged between Early Carboniferous and Late Permian. However, the southward subduction was short-lived and by the Late Permian the Gondwana margin became passive. (C) 2016 Elsevier B.V. All rights reserved.
Strong waves in the mid-latitude circulation have been linked to extreme surface weather and thus changes in waviness could have serious consequences for society. Several theories have been proposed which could alter waviness, including tropical sea surface temperature anomalies or rapid climate change in the Arctic. However, so far it remains unclear whether any changes in waviness have actually occurred. Here we propose a novel meandering index which captures the maximum waviness in geopotential height contours at any given day, using all information of the full spatial position of each contour. Data are analysed on different time scale (from daily to 11 day running means) and both on hemispheric and regional scales. Using quantile regressions, we analyse how seasonal distributions of this index have changed over 1979-2015. The most robust changes are detected for autumn which has seen a pronounced increase in strongly meandering patterns at the hemispheric level as well as over the Eurasian sector. In summer for both the hemisphere and the Eurasian sector, significant downward trends in meandering are detected on daily timescales which is consistent with the recently reported decrease in summer storm track activity. The American sector shows the strongest increase in meandering in the warm season: in particular for 11 day running mean data, indicating enhanced amplitudes of quasi-stationary waves. Our findings have implications for both the occurrence of recent cold spells and persistent heat waves in the mid-latitudes.
Pollen influx (number of pollen grains cm−2 year−1) can objectively reflect the dispersal and deposition features of pollen within a certain time and space, and is often used as a basis for the quantitative reconstruction of palaeovegetation; however, little is known about the features and mechanisms of vertical dispersal of pollen. Here we present the results from a 5 year (2006–2010) monitoring program using pollen traps placed at different heights from ground level up to 60 m and surface soil samples in a mixed coniferous and deciduous broad-leaved woodland in the Changbai mountains, northeastern China. The pollen percentages and pollen influx from the traps have very similar characteristics to the highest values for Betula, Fraxinus, Quercus and Pinus, among the tree taxa and Artemisia, Chenopodiaceae and Asteraceae among the herb taxa. Pollen influx values vary significantly with height and show major differences between three distinct layers, above-canopy (≥32 m), within the trunk layer (8 ≤ 32 m) and on the ground (0 m). These differences in pollen influx are explained by differences in (i) the air flows in each of these layers and (ii) the fall speed of pollen of the various taxa. We found that the pollen recorded on the ground surface is a good representation of the major part of the pollen transported in the trunk space of the woodland. Comparison of the pollen influx values with the theoretical, calculated “characteristic pollen source area” (CPSA) of 12 selected taxa indicates that the pollen deposited on the ground surface of the woodland is a fair representation with 85–90 % of the total pollen deposited at a wind speed of 2.4 m s−1 coming from within ca. 1–5 km for Pinus and Quercus, ca. 5–10 km for Ulmus, Tilia, Oleaceae and Betula, ca. 20–40 km for Fraxinus, Poaceae, Chenopodiaceae, Populus and Salix, and ca. 30–60 km for Artemisia; it is also a good representation with 90–98 % of the total pollen deposited coming from within 60 km at a wind speed of 2.4 m s−1, or 100 km at a wind speed: 6 m s−1, for the 12 selected taxa used in the CPSA calculation. Furthermore, comparison with the vegetation map of the area around the sampling site shows that the pollen deposited on the ground represents all plant communities which grow in the study area within 70 km radius of the sampling site. In this study, the pollen percentages obtained from the soil surface samples are significantly biased towards pollen taxa with good preservation due to thick and robust pollen walls. Therefore, if mosses are available instead, soil samples should be avoided for pollen studies, in particular for the study of pollen-vegetation relationships, the estimation of pollen productivities and quantitative reconstruction of past vegetation. The results also indicate that the existing model of pollen dispersal and deposition, Prentice’s model, provides a fair description of the actual pollen dispersal and deposition in this kind of woodland, which suggests that the application of the landscape reconstruction algorithm would be relevant for reconstruction of this type of woodland in the past.
A link between chemical weathering and physical erosion exists at the catchment scale over a wide range of erosion rates(1,2). However, in mountain environments, where erosion rates are highest, weathering may be kinetically limited(3-5) and therefore decoupled from erosion. In active mountain belts, erosion is driven by bedrock landsliding(6) at rates that depend strongly on the occurrence of extreme rainfall or seismicity(7). Although landslides affect only a small proportion of the landscape, bedrock landsliding can promote the collection and slow percolation of surface runoff in highly fragmented rock debris and create favourable conditions for weathering. Here we show from analysis of surface water chemistry in the Southern Alps of New Zealand that weathering in bedrock landslides controls the variability in solute load of these mountain rivers. We find that systematic patterns in surface water chemistry are strongly associated with landslide occurrence at scales from a single hillslope to an entire mountain belt, and that landslides boost weathering rates and river solute loads over decades. We conclude that landslides couple erosion and weathering in fast-eroding uplands and, thus, mountain weathering is a stochastic process that is sensitive to climatic and tectonic controls on mass wasting processes.
Cities play a vital role in the global climate change mitigation agenda. City population density is one of the key factors that influence urban energy consumption and the subsequent GHG emissions. However, previous research on the relationship between population density and GHG emissions led to contradictory results due to urban/rural definition conundrum and the varying methodologies for estimating GHG emissions. This work addresses these ambiguities by employing the City Clustering Algorithm (CCA) and utilizing the gridded CO2 emissions data. Our results, derived from the analysis of all inhabited areas in the US, show a sub-linear relationship between population density and the total emissions (i.e. the sum of on-road and building emissions) on a per capita basis. Accordingly, we find that doubling the population density would entail a reduction in the total CO2 emissions in buildings and on-road sectors typically by at least 42%. Moreover, we find that population density exerts a higher influence on on-road emissions than buildings emissions. From an energy consumption point of view, our results suggest that on-going urban sprawl will lead to an increase in on-road energy consumption in cities and therefore stresses the importance of developing adequate local policy measures to limit urban sprawl. (C) 2016 Elsevier Ltd. All rights reserved.
The energy sector is both affected by climate change and a key sector for climate protection measures. Energy security is the backbone of our modern society and guarantees the functioning of most critical infrastructure. Thus, decision makers and energy suppliers of different countries should be familiar with the factors that increase or decrease the susceptibility of their electricity sector to climate change. Susceptibility means socioeconomic and structural characteristics of the electricity sector that affect the demand for and supply of electricity under climate change. Moreover, the relevant stakeholders are supposed to know whether the given national energy and climate targets are feasible and what needs to be done in order to meet these targets. In this regard, a focus should be on the residential building sector as it is one of the largest energy consumers and therefore emitters of anthropogenic CO 2 worldwide.
This dissertation addresses the first aspect, namely the susceptibility of the electricity sector, by developing a ranked index which allows for quantitative comparison of the electricity sector susceptibility of 21 European countries based on 14 influencing factors. Such a ranking has not been completed to date. We applied a sensitivity analysis to test the relative effect of each influencing factor on the susceptibility index ranking. We also discuss reasons for the ranking position and thus the susceptibility of selected countries. The second objective, namely the impact of climate change on the energy demand of buildings, is tackled by means of a new model with which the heating and cooling energy demand of residential buildings can be estimated. We exemplarily applied the model to Germany and the Netherlands. It considers projections of future changes in population, climate and the insulation standards of buildings, whereas most of the existing studies only take into account fewer than three different factors that influence the future energy demand of buildings. Furthermore, we developed a comprehensive retrofitting algorithm with which the total residential building stock can be modeled for the first time for each year in the past and future.
The study confirms that there is no correlation between the geographical location of a country and its position in the electricity sector susceptibility ranking. Moreover, we found no pronounced pattern of susceptibility influencing factors between countries that ranked higher or lower in the index. We illustrate that Luxembourg, Greece, Slovakia and Italy are the countries with the highest electricity sector susceptibility. The electricity sectors of Norway, the Czech Republic, Portugal and Denmark were found to be least susceptible to climate change. Knowledge about the most important factors for the poor and good ranking positions of these countries is crucial for finding adequate adaptation measures to reduce the susceptibility of the electricity sector. Therefore, these factors are described within this study.
We show that the heating energy demand of residential buildings will strongly decrease in both Germany and the Netherlands in the future. The analysis for the Netherlands focused on the regional level and a finer temporal resolution which revealed strong variations in the future heating energy demand changes by province and by month. In the German study, we additionally investigated the future cooling energy demand and could demonstrate that it will only slightly increase up to the middle of this century. Thus, increases in the cooling energy demand are not expected to offset reductions in heating energy demand. The main factor for substantial heating energy demand reductions is the retrofitting of buildings. We are the first to show that the given German and Dutch energy and climate targets in the building sector can only be met if the annual retrofitting rates are substantially increased. The current rate of only about 1 % of the total building stock per year is insufficient for reaching a nearly zero-energy demand of all residential buildings by the middle of this century. To reach this target, it would need to be at least tripled. To sum up, this thesis emphasizes that country-specific characteristics are decisive for the electricity sector susceptibility of European countries. It also shows for different scenarios how much energy is needed in the future to heat and cool residential buildings. With this information, existing climate mitigation and adaptation measures can be justified or new actions encouraged.
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.
Extreme hydro-meteorological events, such as severe droughts or heavy rainstorms, constitute primary manifestations of climate variability and exert a critical impact on the natural environment and human society. This is particularly true for high-mountain areas, such as the eastern flank of the southern Central Andes of NW Argentina, a region impacted by deep convection processes that form the basis of extreme events, often resulting in floods, a variety of mass movements, and hillslope processes. This region is characterized by pronounced E-W gradients in topography, precipitation, and vegetation cover, spanning low to medium-elevation, humid and densely vegetated areas to high-elevation, arid and sparsely vegetated environments. This strong E-W gradient is mirrored by differences in the efficiency of surface processes, which mobilize and transport large amounts of sediment through the fluvial system, from the steep hillslopes to the intermontane basins and further to the foreland. In a highly sensitive high-mountain environment like this, even small changes in the spatiotemporal distribution, magnitude and rates of extreme events may strongly impact environmental conditions, anthropogenic activity, and the well-being of mountain communities and beyond. However, although the NW Argentine Andes comprise the catchments for the La Plata river that traverses one of the most populated and economically relevant areas of South America, there are only few detailed investigations of climate variability and extreme hydro-meteorological events.
In this thesis, I focus on deciphering the spatiotemporal variability of rainfall and river discharge, with particular emphasis on extreme hydro-meteorological events in the subtropical southern Central Andes of NW Argentina during the past seven decades. I employ various methods to assess and quantify statistically significant trend patterns of rainfall and river discharge, integrating high-quality daily time series from gauging stations (40 rainfall and 8 river discharge stations) with gridded datasets (CPC-uni and TRMM 3B42 V7), for the period between 1940 and 2015. Evidence for a general intensification of the hydrological cycle at intermediate elevations (~ 0.5 – 3 km asl) at the eastern flank of the southern Central Andes is found both from rainfall and river-discharge time-series analysis during the period from 1940 to 2015. This intensification is associated with the increase of the annual total amount of rainfall and the mean annual discharge. However, most pronounced trends are found at high percentiles, i.e. extreme hydro-meteorological events, particularly during the wet season from December to February.An important outcome of my studies is the recognition of a rapid increase in the amount of river discharge during the period between 1971 and 1977, most likely linked to the 1976-77 global climate shift, which is associated with the North Pacific Ocean sea surface temperature variability. Interestingly, after this rapid increase, both rainfall and river discharge decreased at low and intermediate elevations along the eastern flank of the Andes. In contrast, during the same time interval, at high elevations, extensive areas on the arid Puna de Atacama plateau have recorded increasing annual rainfall totals. This has been associated with more intense extreme hydro-meteorological events from 1979 to 2014. This part of the study reveals that low-, intermediate, and high-elevation sectors in the Andes of NW Argentina respond differently to changing climate conditions.
Possible forcing mechanisms of the pronounced hydro-meteorological variability observed in the study area are also investigated. For the period between 1940 and 2015, I analyzed modes of oscillation of river discharge from small to medium drainage basins (102 to 104 km2), located on the eastern flank of the orogen. First, I decomposed the relevant monthly time series using the Hilbert-Huang Transform, which is particularly appropriate for non-stationary time series that result from non-linear natural processes. I observed that in the study region discharge variability can be described by five quasi-periodic oscillatory modes on timescales varying from 1 to ~20 years. Secondly, I tested the link between river-discharge variations and large-scale climate modes of variability, using different climate indices, such as the BEST ENSO (Bivariate El Niño-Southern Oscillation Time-series) index. This analysis reveals that, although most of the variance on the annual timescale is associated with the South American Monsoon System, a relatively large part of river-discharge variability is linked to Pacific Ocean variability (PDO phases) at multi-decadal timescales (~20 years). To a lesser degree, river discharge variability is also linked to the Tropical South Atlantic (TSA) sea surface temperature anomaly at multi-annual timescales (~2-5 years).
Taken together, these findings exemplify the high degree of sensitivity of high-mountain environments with respect to climatic variability and change. This is particularly true for the topographic transitions between the humid, low-moderate elevations and the semi-arid to arid highlands of the southern Central Andes. Even subtle changes in the hydro-meteorological regime of these areas of the mountain belt react with major impacts on erosional hillslope processes and generate mass movements that fundamentally impact the transport capacity of mountain streams. Despite more severe storms in these areas, the fluvial system is characterized by pronounced variability of the stream power on different timescales, leading to cycles of sediment aggradation, the loss of agriculturally used land and severe impacts on infrastructure.
Deciphering the response of sediment routing systems to climatic forcing is fundamental for understanding the impacts of climate change on landscape evolution. In the Kangra Basin (northwest Sub-Himalaya, India), upper Pleistocene to Holocene alluvial fills and fluvial terraces record periodic fluctuations of sediment supply and transport capacity on timescales of 10(3) to 10(5) yr. To evaluate the potential influence of climate change on these fluctuations, we compare the timing of aggradation and incision phases recorded within remnant alluvial fans and terraces with climate archives. New surface-exposure dating of six terrace levels with in-situ cosmogenic Be-10 indicates the onset of incision phases. Two terrace surfaces from the highest level (T1) sculpted into the oldest preserved alluvial fan (AF1) date back to 53.4 +/- 3.2 ka and 43.0 +/- 2.7 ka (1 sigma). T2 surfaces sculpted into the remnants of AF1 have exposure ages of 18.6 +/- 1.2 ka and 15.3 +/- 0.9 ka, while terraces sculpted into the upper Pleistocene-Holocene fan (AF2) provide ages of 9.3 +/- 0.4 ka (T3), 7.1 +/- 0.4 ka (T4), 5.2 +/- 0.4 ka (T5) and 3.6 +/- 0.2 ka (T6). Together with previously published OSL ages yielding the timing of aggradation, we find a correlation between variations in sediment transport with oxygen-isotope records from regions affected by the Indian Summer Monsoon. During periods of increased monsoon intensity and post-Last Glacial Maximum glacial retreat, aggradation occurred in the Kangra Basin, likely due to high sediment flux, whereas periods of weakened monsoon intensity or lower sediment supply coincide with incision. (C) 2016 Elsevier B.V. All rights reserved.
Lithic sandstones of braided-fluvial to supratidal facies in the Paleoarchean Moodies Group (similar to 3.22 Ga, Barberton Greenstone Belt, South Africa) include several regionally traceable units with common to abundant, in places rock-forming, nodular concretions of megaquartz pseudomorphs after gypsum, barite and calcite. Concretionary accumulations are stratiform and commonly associated with aqueously reworked, fine-grained, tuffaceous sediment of originally rhyodacitic composition and can grow to fist sized agglomerates in crusts tens of m in lateral extent. Weathering of tuffaceous material and feldspar delivered alkali cations such as Ca, Ba, and K, while carbonates were likely supplied by silicate weathering of mafic to ultramafic volcanic rocks during exposure to a CO2-rich atmosphere. Sulfate ions were partly delivered by oxidative pyrite dissolution which may have included microbial and abiotic disproportionation of volcanic S or SO2. Concretionary growth apparently took place under pedogenic to early diagenetic conditions within unconsolidated granular sediment in the vadose zone, dominated by seasonal fluctuations of the groundwater level under evaporitic conditions. The concretions likely represent the oldest terrestrial evaporites known to date and form part of the oldest known compound paleosols. Their formation and composition constrain the local occurrence of sulfate in the Archean atmo- and hydrosphere, their interaction with the emerging biosphere, Archean weathering regime, local climate, and vadose-zone hydrodynamics. (C) 2016 Elsevier B.V. All rights reserved.
Climatic controls on debris-flow activity and sediment aggradation: The Del Medio fan, NW Argentina
(2016)
In the Central Andes, several studies on alluvial terraces and valley fills have linked sediment aggradation to periods of enhanced sediment supply. However, debate continues over whether tectonic or climatic factors are most important in triggering the enhanced supply. The Del Medio catchment in the Humahuaca Basin (Eastern Cordillera, NW Argentina) is located within a transition zone between subhumid and arid climates and hosts the only active debris-flow fan within this intermontane valley. By combining Be-10 analyses of boulder and sediment samples within the Del Medio catchment, with regional morphometric measurements of nearby catchments, we identify the surface processes responsible for aggradation in the Del Medio fan and their likely triggers. We find that the fan surface has been shaped by debris flows and channel avulsions during the last 400 years. Among potential tectonic, climatic, and autogenic factors that might influence deposition, our analyses point to a combination of several favorable factors that drive aggradation. These are in particular the impact of occasional abundant rainfall on steep slopes in rock types prone to failure, located in a region characterized by relatively low rainfall amounts and limited transport capacity. These characteristics are primarily associated with the climatic transition zone between the humid foreland and the arid orogen interior, which creates an imbalance between sediment supply and sediment transfer. The conditions and processes that drive aggradation in the Del Medio catchment today may provide a modern analog for the conditions and processes that drove aggradation in other nearby tributaries in the past.
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.
Knowledge of the possible impacts of early diagenesis on lipid biomarkers in geologic settings is important for robust applications of lipid proxies for paleoclimate reconstructions. In this study, molecular distributions and carbon isotopic compositions of lipids were compared in two particle-size fractions (<0.3 mm and >0.3 mm) of twelve surface peat samples collected from Dajiuhu peatland, central China. The average chain length (ACL) values of long-chain n-alkanes, n-fatty alcohols, n-fatty acids and n-alkan-2-ones show no significant differences between the finer and coarser fractions. In contrast, the carbon preference index values of long-chain n-alkanes, n-fatty alcohols and n-alkan-2-ones have relatively smaller values in the finer fractions than in the coarser ones. Stanols were also more abundant in the finer fractions. In addition, the delta C-13 values of odd-numbered n-alkanes (C-23-C-33) were generally less negative in the finer fractions. Our results indicate that (1) the finer fractions probably experienced stronger degradation than the coarser fractions; (2) the less negative delta(CC)-C-13 values of odd-numbered n-alkanes (C-23-C-33) in the finer fractions were largely a result of greater heterotrophic reworking during degradation; (3) ACL values of long-chain n-alkyl lipids (n-alkanes, n-fatty alcohols and n-fatty acids, n-alkan-2-ones) appear to be reliable proxies to trace lipid sources and their associated paleoenvironmental signals in peat deposits.
The uppermost Cretaceous (upper Campanian-Maastrichtian) marine deposits of the central south Pyrenees host a rich larger benthic foraminiferal fauna and several rudist-rich levels. These marine deposits are directly overlain by the continental facies of the Aren and Tremp Formations, which are famous for their fossil dinosaur remains. Larger benthic foraminiferal distribution documents an important faunal turnover in all the carbonate platform environments within the photic zone, from open marine to littoral areas. Biostratigraphy indicates that this turnover occurred close to the Campanian-Maastrichtian boundary. This is also confirmed by strontium isotope stratigraphy which indicates an earliest Maastrichtian age for the appearance of the larger benthic foraminiferal assemblage constituted by Lepidorbitoides socialis, Clypeorbis mammillata, Wannierina cataluniensis, Orbitoides gruenbachensis, Siderolites aff. calcitrapoides, Fascispira colomi, Omphalocyclus macroporus and Laffiteina mengaudi. In particular, a numerical age of 71 Ma is obtained for the Hippurites radiosus level, just a few meters below the first continental deposits of the Aren sensu stricto Formation. The youngest marine sediments of the central south Pyrenees are early Maastrichtian in age. This is also an important constraint for the age of the end-Cretaceous dinosaur fossil localities of the Tremp basin. (C) 2015 Elsevier Ltd. All rights reserved.
The present-day stress state is a key parameter in numerous geoscientific research fields including geodynamics, seismic hazard assessment, and geomechanics of georeservoirs. The Taranaki Basin of New Zealand is located on the Australian Plate and forms the western boundary of tectonic deformation due to Pacific Plate subduction along the Hikurangi margin. This paper presents the first comprehensive wellbore-derived basin-scale in situ stress analysis in New Zealand. We analyze borehole image and oriented caliper data from 129 petroleum wells in the Taranaki Basin to interpret the shape of boreholes and determine the orientation of maximum horizontal stress (S-Hmax). We combine these data (151 S-Hmax data records) with 40 stress data records derived from individual earthquake focal mechanism solutions, 6 from stress inversions of focal mechanisms, and 1 data record using the average of several focal mechanism solutions. The resulting data set has 198 data records for the Taranaki Basin and suggests a regional S-Hmax orientation of N068 degrees E (22 degrees), which is in agreement with NW-SE extension suggested by geological data. Furthermore, this ENE-WSW average S-Hmax orientation is subparallel to the subduction trench and strike of the subducting slab (N50 degrees E) beneath the central western North Island. Hence, we suggest that the slab geometry and the associated forces due to slab rollback are the key control of crustal stress in the Taranaki Basin. In addition, we find stress perturbations with depth in the vicinity of faults in some of the studied wells, which highlight the impact of local stress sources on the present-day stress rotation.
Pluvial floods have caused severe damage to urban areas in recent years. With a projected increase in extreme precipitation as well as an ongoing urbanization, pluvial flood damage is expected to increase in the future. Therefore, further insights, especially on the adverse consequences of pluvial floods and their mitigation, are needed. To gain more knowledge, empirical damage data from three different pluvial flood events in Germany were collected through computer-aided telephone interviews. Pluvial flood awareness as well as flood experience were found to be low before the respective flood events. The level of private precaution increased considerably after all events, but is mainly focused on measures that are easy to implement. Lower inundation depths, smaller potential losses as compared with fluvial floods, as well as the fact that pluvial flooding may occur everywhere, are expected to cause a shift in damage mitigation from precaution to emergency response. However, an effective implementation of emergency measures was constrained by a low dissemination of early warnings in the study areas. Further improvements of early warning systems including dissemination as well as a rise in pluvial flood preparedness are important to reduce future pluvial flood damage.
Water scarcity, adaption on climate change, and risk assessment of droughts and floods are critical topics for science and society these days. Monitoring and modeling of the hydrological cycle are a prerequisite to understand and predict the consequences for weather and agriculture. As soil water storage plays a key role for partitioning of water fluxes between the atmosphere, biosphere, and lithosphere, measurement techniques are required to estimate soil moisture states from small to large scales.
The method of cosmic-ray neutron sensing (CRNS) promises to close the gap between point-scale and remote-sensing observations, as its footprint was reported to be 30 ha. However, the methodology is rather young and requires highly interdisciplinary research to understand and interpret the response of neutrons to soil moisture. In this work, the signal of nine detectors has been systematically compared, and correction approaches have been revised to account for meteorological and geomagnetic variations. Neutron transport simulations have been consulted to precisely characterize the sensitive footprint area, which turned out to be 6--18 ha, highly local, and temporally dynamic. These results have been experimentally confirmed by the significant influence of water bodies and dry roads. Furthermore, mobile measurements on agricultural fields and across different land use types were able to accurately capture the various soil moisture states. It has been further demonstrated that the corresponding spatial and temporal neutron data can be beneficial for mesoscale hydrological modeling. Finally, first tests with a gyrocopter have proven the concept of airborne neutron sensing, where increased footprints are able to overcome local effects.
This dissertation not only bridges the gap between scales of soil moisture measurements. It also establishes a close connection between the two worlds of observers and modelers, and further aims to combine the disciplines of particle physics, geophysics, and soil hydrology to thoroughly explore the potential and limits of the CRNS method.
Volunteered geographical information (VGI) and citizen science have become important sources data for much scientific research. In the domain of land cover, crowdsourcing can provide a high temporal resolution data to support different analyses of landscape processes. However, the scientists may have little control over what gets recorded by the crowd, providing a potential source of error and uncertainty. This study compared analyses of crowdsourced land cover data that were contributed by different groups, based on nationality (labelled Gondor and Non-Gondor) and on domain experience (labelled Expert and Non-Expert). The analyses used a geographically weighted model to generate maps of land cover and compared the maps generated by the different groups. The results highlight the differences between the maps how specific land cover classes were under-and over-estimated. As crowdsourced data and citizen science are increasingly used to replace data collected under the designed experiment, this paper highlights the importance of considering between group variations and their impacts on the results of analyses. Critically, differences in the way that landscape features are conceptualised by different groups of contributors need to be considered when using crowdsourced data in formal scientific analyses. The discussion considers the potential for variation in crowdsourced data, the relativist nature of land cover and suggests a number of areas for future research. The key finding is that the veracity of citizen science data is not the critical issue per se. Rather, it is important to consider the impacts of differences in the semantics, affordances and functions associated with landscape features held by different groups of crowdsourced data contributors.
Despite a long history of plate convergence at the western margin of the South American plate that has been ongoing since at least the Early Paleozoic, the southern Peruvian fore-arc displays little to no evidence of shortening. In the light of this observation, we assess the deformation history of the southern Peruvian fore-arc and its geodynamic implications. To accomplish this, we present a new structural and geo-thermochronological data set (zircon U-Pb, mica Ar-40/Ar-39, apatite and zircon fission-track and zircon (U-Th)/He analyses) for samples collected along a 400km long transect parallel to the trench. Our results show that the Mesoproterozoic gneissic basement was mainly at temperatures 350 degrees C since the Neoproterozoic and was later intruded by Jurassic volcanic arc plutons. Along the coast, a peculiar apatite fission-track age pattern, coupled with field observations and a synthesis of available geological maps, allows us to identify crustal-scale tilted blocks that span the coastal Peruvian fore-arc. These blocks, bounded by normal faults that are orthogonal to the trench, suggest post-60Ma trench-parallel extension that potentially accommodated oroclinal bending in this region. Block tilting is consistent with the observed and previously described switch in the location of sedimentary sources in the fore-arc basin. Our data set allows us to estimate the cumulative slip on these faults to be less than 2km and questions the large amount of trench-parallel extension suggested to have accommodated this bending.
Testate amoebae with self-secreted siliceous shell platelets ("idiosomes") play an important role in terrestrial silicon (Si) cycles. In this context, Si-dependent culture growth dynamics of idiosomic testate amoebae are of interest. Clonal cultures of idiosomic testate amoebae were analyzed under three different Si concentrations: low (50 mu mol L-1), moderate/site-specific (150 mu mol L-1) and high Si supply (500 mu mol L-1). Food (Saccharomyces cerevisiae) was provided in surplus. (i) Shell size of four different clones of idiosomic testate amoebae either decreased (Trinema galeata, Euglypha filifera cf.), increased (E. rotunda cf.), or did not change (E. rotunda) under the lowest Si concentration (50 mu mol Si L-1). (ii) Culture growth of idiosomic Euglypha rotunda was dependent on Si concentration. The more Si available in the culture medium, the earlier the entry into exponential growth phase. (iii) Culture growth of idiosomic Euglypha rotunda was dependent on origin of inoculum. Amoebae previously cultured under a moderate Si concentration revealed highest sustainability in consecutive cultures. Amoebae derived from cultures with high Si concentrations showed rapid culture growth which finished early in consecutive cultures. (iv) Si (diluted in the culture medium) was absorbed by amoebae and fixed in the amoeba shells resulting in decreased Si concentrations. (C) 2016 Elsevier GmbH. All rights reserved.
Ecosystem services have a significant impact on human wellbeing. While ecosystem services are frequently represented by monetary values, social values and underlying social benefits remain underexplored. The purpose of this study is to assess whether and how social benefits have been explicitly addressed within socio-economic and socio-cultural ecosystem services research, ultimately allowing a better understanding between ecosystem services and human well-being. In this paper, we reviewed 115 international primary valuation studies and tested four hypotheses associated to the identification of social benefits of ecosystem services using logistic regressions. Tested hypotheses were that (1) social benefits are mostly derived in studies that assess cultural ecosystem services as opposed to other ecosystem service types, (2) there is a pattern of social benefits and certain cultural ecosystem services assessed simultaneously, (3) monetary valuation techniques go beyond expressing monetary values and convey social benefits, and (4) directly addressing stakeholdeŕs views the consideration of social benefits in ecosystem service assessments. Our analysis revealed that (1) a variety of social benefits are valued in studies that assess either of the four ecosystem service types, (2) certain social benefits are likely to co-occur in combination with certain cultural ecosystem services, (3) of the studies that employed monetary valuation techniques, simulated market approaches overlapped most frequently with the assessment of social benefits and (4) studies that directly incorporate stakeholder's views were more likely to also assess social benefits.
Most climate change impacts manifest in the form of natural hazards. Damage assessment typically relies on damage functions that translate the magnitude of extreme events to a quantifiable damage. In practice, the availability of damage functions is limited due to a lack of data sources and a lack of understanding of damage processes. The study of the characteristics of damage functions for different hazards could strengthen the theoretical foundation of damage functions and support their development and validation. Accordingly, we investigate analogies of damage functions for coastal flooding and for wind storms and identify a unified approach. This approach has general applicability for granular portfolios and may also be applied, for example, to heat-related mortality. Moreover, the unification enables the transfer of methodology between hazards and a consistent treatment of uncertainty. This is demonstrated by a sensitivity analysis on the basis of two simple case studies (for coastal flood and storm damage). The analysis reveals the relevance of the various uncertainty sources at varying hazard magnitude and on both the microscale and the macroscale level. Main findings are the dominance of uncertainty from the hazard magnitude and the persistent behaviour of intrinsic uncertainties on both scale levels. Our results shed light on the general role of uncertainties and provide useful insight for the application of the unified approach.
The Greater and Lesser Caucasus mountains and their associated foreland basins contain similar rock types, experience a similar two-fold, along-strike variation in mean annual precipitation, and were affected by extreme base-level drops of the neighboring Caspian Sea. However, the two Caucasus ranges are characterized by decidedly different tectonic regimes and rates of deformation that are subject to moderate (less than an order of magnitude) gradients in climate, and thus allow for a unique opportunity to isolate the effects of climate and tectonics in the evolution of topography within active orogens. There is an apparent disconnect between modern climate, shortening rates, and topography of both the Greater Caucasus and Lesser Caucasus which exhibit remarkably similar topography along-strike despite the gradients in forcing. By combining multiple datasets, we examine plausible causes for this disconnect by presenting a detailed analysis of the topography of both ranges utilizing established relationships between catchment-mean erosion rates and topography (local relief, hillslope gradients, and channel steepness) and combining it with a synthesis of previously published low-temperature thermochronologic data. Modern climate of the Caucasus region is assessed through an analysis of remotely-sensed data (TRMM and MODIS) and historical streamflow data. Because along-strike variation in either erosional efficiency or thickness of accreted material fail to explain our observations, we suggest that the topography of both the western Lesser and Greater Caucasus are partially supported by different geodynamic forces. In the western Lesser Caucasus, high relief portions of the landscape likely reflect uplift related to ongoing mantle lithosphere delamination beneath the neighboring East Anatolian Plateau. In the Greater Caucasus, maintenance of high topography in the western portion of the range despite extremely low (<2-4 mm/y) modern convergence rates may be related to dynamic topography from detachment of the north-directed Greater Caucasus slab or to a recent slowing of convergence rates. Large-scale spatial gradients in climate are not reflected in the topography of the Caucasus and do not seem to exert any significant control on the tectonics or structure of either range. (C) 2016 Elsevier B.V. All rights reserved.
Deglaciation of Fennoscandia
(2016)
To provide a new reconstruction of the deglaciation of the Fennoscandian Ice Sheet, in the form of calendar-year time-slices, which are particularly useful for ice sheet modelling, we have compiled and synthesized published geomorphological data for eskers, ice-marginal formations, lineations, marginal meltwater channels, striae, ice-dammed lakes, and geochronological data from radiocarbon, varve, optically-stimulated luminescence, and cosmogenic nuclide dating. This is summarized as a deglaciation map of the Fennoscandian Ice Sheet with isochrons marking every 1000 years between 22 and 13 cal kyr BP and every hundred years between 11.6 and final ice decay after 9.7 cal kyr BP. Deglaciation patterns vary across the Fennoscandian Ice Sheet domain, reflecting differences in climatic and geomorphic settings as well as ice sheet basal thermal conditions and terrestrial versus marine margins. For example, the ice sheet margin in the high-precipitation coastal setting of the western sector responded sensitively to climatic variations leaving a detailed record of prominent moraines and other ice-marginal deposits in many fjords and coastal valleys. Retreat rates across the southern sector differed between slow retreat of the terrestrial margin in western and southern Sweden and rapid retreat of the calving ice margin in the Baltic Basin. Our reconstruction is consistent with much of the published research. However, the synthesis of a large amount of existing and new data support refined reconstructions in some areas. For example, the LGM extent of the ice sheet in northwestern Russia was located far east and it occurred at a later time than the rest of the ice sheet, at around 17-15 cal kyr BP. We also propose a slightly different chronology of moraine formation over southern Sweden based on improved correlations of moraine segments using new LiDAR data and tying the timing of moraine formation to Greenland ice core cold stages. Retreat rates vary by as much as an order of magnitude in different sectors of the ice sheet, with the lowest rates on the high-elevation and maritime Norwegian margin. Retreat rates compared to the climatic information provided by the Greenland ice core record show a general correspondence between retreat rate and climatic forcing, although a close match between retreat rate and climate is unlikely because of other controls, such as topography and marine versus terrestrial margins. Overall, the time slice reconstructions of Fennoscandian Ice Sheet deglaciation from 22 to 9.7 cal kyr BP provide an important dataset for understanding the contexts that underpin spatial and temporal patterns in retreat of the Fennoscandian Ice Sheet, and are an important resource for testing and refining ice sheet models. (C) 2015 The Authors. Published by Elsevier Ltd.
Das Wissen um die lokale Struktur von Seltenen Erden Elementen (SEE) in silikatischen und aluminosilikatischen Schmelzen ist von fundamentalem Interesse für die Geochemie der magmatischen Prozesse, speziell wenn es um ein umfassendes Verständnis der Verteilungsprozesse von SEE in magmatischen Systemen geht. Es ist allgemein akzeptiert, dass die SEE-Verteilungsprozesse von Temperatur, Druck, Sauerstofffugazität (im Fall von polyvalenten Kationen) und der Kristallchemie kontrolliert werden. Allerdings ist wenig über den Einfluss der Schmelzzusammensetzung selbst bekannt. Ziel dieser Arbeit ist, eine Beziehung zwischen der Variation der SEE-Verteilung mit der Schmelzzusammensetzung und der Koordinationschemie dieser SEE in der Schmelze zu schaffen.
Dazu wurden Schmelzzusammensetzungen von Prowatke und Klemme (2005), welche eine deutliche Änderung der Verteilungskoeffizienten zwischen Titanit und Schmelze ausschließlich als Funktion der Schmelzzusammensetzung zeigen, sowie haplogranitische bzw. haplobasaltische Schmelzzusammensetzungen als Vertreter magmatischer Systeme mit La, Gd, Yb und Y dotiert und als Glas synthetisiert. Die Schmelzen variierten systematisch im Aluminiumsättigungsindex (ASI), welcher bei den Prowatke und Klemme (2005) Zusammensetzungen einen Bereich von 0.115 bis 0.768, bei den haplogranitischen Zusammensetzungen einen Bereich von 0.935 bis 1.785 und bei den haplobasaltischen Zusammensetzungen einen Bereich von 0.368 bis 1.010 abdeckt. Zusätzlich wurden die haplogranitischen Zusammensetzungen mit 4 % H2O synthetisiert, um den Einfluss von Wasser auf die lokale Umgebung von SEE zu studieren. Um Informationen über die lokalen Struktur von Gd, Yb und Y zu erhalten wurde die Röntgenabsorptionsspektroskopie angewendet. Dabei liefert die Untersuchung der Feinstruktur mittels der EXAFS-Spektroskopie (engl. Extended X-Ray Absorption Fine Structure) quantitative Informationen über die lokale Umgebung, während RIXS (engl. resonant inelastic X-ray scattering), sowie die daraus extrahierte hoch aufgelöste Nahkantenstruktur, XANES (engl. X-ray absorption near edge structure) qualitative Informationen über mögliche Koordinationsänderungen von La, Gd und Yb in den Gläsern liefert. Um mögliche Unterschiede der lokalen Struktur oberhalb der Glastransformationstemperatur (TG) zur Raumtemperatur zu untersuchen, wurden exemplarisch Hochtemperatur Y-EXAFS Untersuchungen durchgeführt.
Für die Auswertung der EXAFS-Messungen wurde ein neu eingeführter Histogramm-Fit verwendet, der auch nicht-symmetrische bzw. nichtgaußförmige Paarverteilungsfunktionen beschreiben kann, wie sie bei einem hohen Grad der Polymerisierung bzw. bei hohen Temperaturen auftreten können. Die Y-EXAFS-Spektren für die Prowatke und Klemme (2005) Zusammensetzungen zeigen mit Zunahme des ASI, eine Zunahme der Asymmetrie und Breite der Y-O Paarverteilungsfunktion, welche sich in sich in der Änderung der Koordinationszahl von 6 nach 8 und einer Zunahme des Y-O Abstand um 0.13Å manifestiert. Ein ähnlicher Trend lässt sich auch für die Gd- und Yb-EXAFS-Spektren beobachten. Die hoch aufgelösten XANESSpektren für La, Gd und Yb zeigen, dass sich die strukturellen Unterschiede zumindest halb-quantitativ bestimmen lassen. Dies gilt insbesondere für Änderungen im mittleren Abstand zu den Sauerstoffatomen. Im Vergleich zur EXAFS-Spektroskopie liefert XANES jedoch keine Informationen über die Form und Breite von Paarverteilungsfunktionen. Die Hochtemperatur EXAFS-Untersuchungen von Y zeigen Änderungen der lokalen Struktur oberhalb der Glasübergangstemperatur an, welche sich vordergründig auf eine thermisch induzierte Erhöhung des mittleren Y-O Abstandes zurückführen lassen. Allerdings zeigt ein Vergleich der Y-O Abstände für Zusammensetzungen mit einem ASI von 0.115 bzw. 0.755, ermittelt bei Raumtemperatur und TG, dass der im Glas beobachtete strukturelle Unterschied entlang der Zusammensetzungsserie in der Schmelze noch stärker ausfallen kann, als bisher für die Gläser angenommen wurde.
Die direkte Korrelation der Verteilungsdaten von Prowatke und Klemme (2005) mit den strukturellen Änderungen der Schmelzen offenbart für Y eine lineare Korrelation, wohingegen Yb und Gd eine nicht lineare Beziehung zeigen. Aufgrund seines Ionenradius und seiner Ladung wird das 6-fach koordinierte SEE in den niedriger polymerisierten Schmelzen bevorzugt durch nicht-brückenbildende Sauerstoffatome koordiniert, um stabile Konfigurationen zu bilden. In den höher polymerisierten Schmelzen mit ASI-Werten in der Nähe von 1 ist 6-fache Koordination nicht möglich, da fast nur noch brückenbildende Sauerstoffatome zur Verfügung stehen. Die Überbindung von brückenbildenden Sauerstoffatomen um das SEE wird durch Erhöhung der Koordinationszahl und des mittleren SEE-O Abstandes ausgeglichen. Dies bedeutet eine energetisch günstigere Konfiguration in den stärker depolymerisierten Zusammensetzungen, aus welcher die beobachtete Variation des Verteilungskoeffizienten resultiert, welcher sich jedoch für jedes Element stark unterscheidet. Für die haplogranitischen und haplobasaltischen Zusammensetzungen wurde mit Zunahme der Polymerisierung auch eine Zunahme der Koordinationszahl und des durchschnittlichen Bindungsabstands, einhergehend mit der Zunahme der Schiefe und der Asymmetrie der Paarverteilungsfunktion, beobachtet. Dies impliziert, dass das jeweilige SEE mit Zunahme der Polymerisierung auch inkompatibler in diesen Zusammensetzungen wird. Weiterhin zeigt die Zugabe von Wasser, dass die Schmelzen depolymerisieren, was in einer symmetrischeren Paarverteilungsfunktion resultiert, wodurch die Kompatibilität wieder zunimmt.
Zusammenfassend zeigt sich, dass die Veränderungen der Schmelzzusammensetzungen in einer Änderung der Polymerisierung der Schmelzen resultieren, die dann einen signifikanten Einfluss auf die lokale Umgebung der SEE hat. Die strukturellen Änderungen lassen sich direkt mit Verteilungsdaten korrelieren, die Trends unterscheiden sich aber stark zwischen leichten, mittleren und schweren SEE. Allerdings konnte diese Studie zeigen, in welcher Größenordnung die Änderungen liegen müssen, um einen signifikanten Einfluss auf den Verteilungskoeffizenten zu haben. Weiterhin zeigt sich, dass der Einfluss der Schmelzzusammensetzung auf die Verteilung der Spurenelemente mit Zunahme der Polymerisierung steigt und daher nicht vernachlässigt werden darf.
In the last decades, an increasing number of studies analyzed spatial patterns in throughfall by means of variograms. The estimation of the variogram from sample data requires an appropriate sampling scheme: most importantly, a large sample and a layout of sampling locations that often has to serve both variogram estimation and geostatistical prediction. While some recommendations on these aspects exist, they focus on Gaussian data and high ratios of the variogram range to the extent of the study area. However, many hydrological data, and throughfall data in particular, do not follow a Gaussian distribution. In this study, we examined the effect of extent, sample size, sampling design, and calculation method on variogram estimation of throughfall data. For our investigation, we first generated non Gaussian random fields based on throughfall data with large outliers. Subsequently, we sampled the fields with three extents (plots with edge lengths of 25 m, 50 m, and 100 m), four common sampling designs (two grid-based layouts, transect and random sampling) and five sample sizes (50, 100, 150, 200, 400). We then estimated the variogram parameters by method-of-moments (non-robust and robust estimators) and residual maximum likelihood. Our key findings are threefold. First, the choice of the extent has a substantial influence on the estimation of the variogram. A comparatively small ratio of the extent to the correlation length is beneficial for variogram estimation. Second, a combination of a minimum sample size of 150, a design that ensures the sampling of small distances and variogram estimation by residual maximum likelihood offers a good compromise between accuracy and efficiency. Third, studies relying on method-of-moments based variogram estimation may have to employ at least 200 sampling points for reliable variogram estimates. These suggested sample sizes exceed the number recommended by studies dealing with Gaussian data by up to 100 %. Given that most previous through fall studies relied on method-of-moments variogram estimation and sample sizes <<200, currently available data are prone to large uncertainties. (C) 2016 Elsevier B.V. All rights reserved.
Im Graduiertenkolleg NatRiskChange der Universität Potsdam und anderen Forschungseinrichtungen werden beobachtete sowie zukünftig mögliche Veränderungen von Naturgefahren untersucht. Teil des strukturierten Doktorandenprogramms sind sogenannte Task-Force-Einsätze, bei denen die Promovierende zeitlich begrenzt ein aktuelles Ereignis auswerten. Im Zuge dieser Aktivität wurde die Sturzflut vom 29.05.2016 in Braunsbach (Baden-Württemberg) untersucht.
In diesem Bericht werden erste Auswertungen zur Einordnung der Niederschläge, zu den hydrologischen und geomorphologischen Prozessen im Einzugsgebiet des Orlacher Bachs sowie zu den verursachten Schäden beleuchtet.
Die Region war Zentrum extremer Regenfälle in der Größenordnung von 100 mm innerhalb von 2 Stunden. Das 6 km² kleine Einzugsgebiet hat eine sehr schnelle Reaktionszeit, zumal bei vorgesättigtem Boden. Im steilen Bachtal haben mehrere kleinere und größere Hangrutschungen über 8000 m³ Geröll, Schutt und Schwemmholz in das Gewässer eingetragen und möglicherweise kurzzeitige Aufstauungen und Durchbrüche verursacht. Neben den großen Wassermengen mit einer Abflussspitze in einer Größenordnung von 100 m³/s hat gerade die Geschiebefracht zu großen Schäden an den Gebäuden entlang des Bachlaufs in Braunsbach geführt.
Robust appraisals of climate impacts at different levels of global-mean temperature increase are vital to guide assessments of dangerous anthropogenic interference with the climate system. The 2015 Paris Agreement includes a two-headed temperature goal: "holding the increase in the global average temperature to well below 2 degrees C above pre-industrial levels and pursuing efforts to limit the temperature increase to 1.5 degrees C". Despite the prominence of these two temperature limits, a comprehensive overview of the differences in climate impacts at these levels is still missing. Here we provide an assessment of key impacts of climate change at warming levels of 1.5 degrees C and 2 degrees C, including extreme weather events, water availability, agricultural yields, sea-level rise and risk of coral reef loss. Our results reveal substantial differences in impacts between a 1.5 degrees C and 2 degrees C warming that are highly relevant for the assessment of dangerous anthropogenic interference with the climate system. For heat-related extremes, the additional 0.5 degrees C increase in global-mean temperature marks the difference between events at the upper limit of present-day natural variability and a new climate regime, particularly in tropical regions. Similarly, this warming difference is likely to be decisive for the future of tropical coral reefs. In a scenario with an end-of-century warming of 2 degrees C, virtually all tropical coral reefs are projected to be at risk of severe degradation due to temperature-induced bleaching from 2050 onwards. This fraction is reduced to about 90% in 2050 and projected to decline to 70% by 2100 for a 1.5 degrees C scenario. Analyses of precipitation-related impacts reveal distinct regional differences and hot-spots of change emerge. Regional reduction in median water availability for the Mediterranean is found to nearly double from 9% to 17% between 1.5 degrees C and 2 degrees C, and the projected lengthening of regional dry spells increases from 7 to 11%. Projections for agricultural yields differ between crop types as well as world regions. While some (in particular high-latitude) regions may benefit, tropical regions like West Africa, South-East Asia, as well as Central and northern South America are projected to face substantial local yield reductions, particularly for wheat and maize. Best estimate sea-level rise projections based on two illustrative scenarios indicate a 50cm rise by 2100 relative to year 2000-levels for a 2 degrees C scenario, and about 10 cm lower levels for a 1.5 degrees C scenario. In a 1.5 degrees C scenario, the rate of sea-level rise in 2100 would be reduced by about 30% compared to a 2 degrees C scenario. Our findings highlight the importance of regional differentiation to assess both future climate risks and different vulnerabilities to incremental increases in global-mean temperature. The article provides a consistent and comprehensive assessment of existing projections and a good basis for future work on refining our understanding of the difference between impacts at 1.5 degrees C and 2 degrees C warming.
Rivers draining the southern Himalaya provide most of the water supply for the densely populated Indo-Gangetic plains. Despite the importance of water resources in light of climate change, the relative contributions of rainfall, snow and glacier melt to discharge are not well understood, due to the scarcity of ground-based data in this complex terrain. Here, we quantify discharge sources in the Sutlej Valley, western Himalaya, from 2000 to 2012 with a distributed hydrological model that is based on daily, ground-calibrated remote-sensing observation. Based on the consistently good model performance, we analyzed the spatiotemporal distribution of hydrologic components and quantified their contribution to river discharge. Our results indicate that the Sutlej River's annual discharge at the mountain front is sourced to 55% by effective rainfall (rainfall reduced by evapotranspiration), 35% by snow melt and 10% by glacier melt. In the high-elevation orogenic interior glacial runoff contributes ∼30% to annual river discharge. These glacier melt contributions are especially important during years with substantially reduced rainfall and snowmelt runoff, as during 2004, to compensate for low river discharge and ensure sustained water supply and hydropower generation. In 2004, discharge of the Sutlej River totaled only half the maximum annual discharge; with 17.3% being sourced by glacier melt. Our findings underscore the importance of calibrating remote-sensing data with ground-based data to constrain hydrological models with reasonable accuracy. For instance, we found that TRMM (Tropical Rainfall Measuring Mission) product 3B42 V7 systematically overestimates rainfall in arid regions of our study area by a factor of up to 5. By quantifying the spatiotemporal distribution of water resources we provide an important assessment of the potential impact of global warming on river discharge in the western Himalaya. Given the near-global coverage of the utilized remote-sensing datasets this hydrological modeling approach can be readily transferred to other data-sparse regions.
Dynamics of mantle plumes
(2016)
Mantle plumes are a link between different scales in the Earth’s mantle: They are an important part of large-scale mantle convection, transporting material and heat from the core-mantle boundary to the surface, but also affect processes on a smaller scale, such as melt generation and transport and surface magmatism. When they reach the base of the lithosphere, they cause massive magmatism associated with the generation of large igneous provinces, and they can be related to mass extinction events (Wignall, 2001) and continental breakup (White and McKenzie, 1989).
Thus, mantle plumes have been the subject of many previous numerical modelling studies (e.g. Farnetani and Richards, 1995; d’Acremont et al., 2003; Lin and van Keken, 2005; Sobolev et al., 2011; Ballmer et al., 2013). However, complex mechanisms, such as the development and implications of chemical heterogeneities in plumes, their interaction with mid-ocean ridges and global mantle flow, and melt ascent from the source region to the surface are still not very well understood; and disagreements between observations and the predictions of classical plume models have led to a challenge of the plume concept in general (Czamanske et al., 1998; Anderson, 2000; Foulger, 2011). Hence, there is a need for more sophisticated models that can explain the underlying physics, assess which properties and processes are important, explain how they cause the observations visible at the Earth’s surface and provide a link between the different scales.
In this work, integrated plume models are developed that investigate the effect of dense recycled oceanic crust on the development of mantle plumes, plume–ridge interaction under the influence of global mantle flow and melting and melt migration in form of two-phase flow.
The presented analysis of these models leads to a new, updated picture of mantle plumes: Models considering a realistic depth-dependent density of recycled oceanic crust and peridotitic mantle material show that plumes with excess temperatures of up to 300 K can transport up to 15% of recycled oceanic crust through the whole mantle. However, due to the high density of recycled crust, plumes can only advance to the base of the lithosphere directly if they have high excess temperatures, high plume volumes and the lowermost mantle is subadiabatic, or plumes rise from the top or edges of thermo-chemical piles. They might only cause minor surface uplift, and instead of the classical head–tail structure, these low-buoyancy plumes are predicted to be broad features in the lower mantle with much less pronounced plume heads. They can form a variety of shapes and regimes, including primary plumes directly advancing to the base of the lithosphere, stagnating plumes, secondary plumes rising from the core–mantle boundary or a pool of eclogitic material in the upper mantle and failing plumes. In the upper mantle, plumes are tilted and deflected by global mantle flow, and the shape, size and stability of the melting region is influenced by the distance from nearby plate boundaries, the speed of the overlying plate and the movement of the plume tail arriving from the lower mantle. Furthermore, the structure of the lithosphere controls where hot material is accumulated and melt is generated. In addition to melting in the plume tail at the plume arrival position, hot plume material flows upwards towards opening rifts, towards mid-ocean ridges and towards other regions of thinner lithosphere, where it produces additional melt due to decompression. This leads to the generation of either broad ridges of thickened magmatic crust or the separation into multiple thinner lines of sea mount chains at the surface. Once melt is generated within the plume, it influences its dynamics, lowering the viscosity and density, and while it rises the melt volume is increased up to 20% due to decompression. Melt has the tendency to accumulate at the top of the plume head, forming diapirs and initiating small-scale convection when the plume reaches the base of the lithosphere. Together with the introduced unstable, high-density material produced by freezing of melt, this provides an efficient mechanism to thin the lithosphere above plume heads.
In summary, this thesis shows that mantle plumes are more complex than previously considered, and linking the scales and coupling the physics of different processes occurring in mantle plumes can provide insights into how mantle plumes are influenced by chemical heterogeneities, interact with the lithosphere and global mantle flow, and are affected by melting and melt migration. Including these complexities in geodynamic models shows that plumes can also have broad plume tails, might produce only negligible surface uplift, can generate one or several volcanic island chains in interaction with a mid–ocean ridge, and can magmatically thin the lithosphere.
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.
Ecosystems' exposure to climate change - Modeling as support for nature conservation management
(2016)
Effects of data and model simplification on the results of a wetland water resource management model
(2016)
This paper presents the development of a wetland water balance model for use in a large river basin with many different wetlands. The basic model was primarily developed for a single wetland with a complex water management system involving large amounts of specialized input data and water management details. The aim was to simplify the model structure and to use only commonly available data as input for the model, with the least possible loss of accuracy. Results from different variants of the model and data adaptation were tested against results from a detailed model. This shows that using commonly available data and unifying and simplifying the input data is tolerable up to a certain level. The simplification of the model has greater effects on the evaluated water balance components than the data adaptation. Because this simplification was necessary for large-scale use, we suggest that, for reasons of comparability, simpler models should always be applied with uniform data bases for large regions, though these should only be moderately simplified. Further, we recommend using these simplified models only for large-scale comparisons and using more specific, detailed models for investigations on smaller scales.
Effects of Data and Model Simplification on the Results of a Wetland Water Resource Management Model
(2016)
This paper presents the development of a wetland water balance model for use in a large river basin with many different wetlands. The basic model was primarily developed for a single wetland with a complex water management system involving large amounts of specialized input data and water management details. The aim was to simplify the model structure and to use only commonly available data as input for the model, with the least possible loss of accuracy. Results from different variants of the model and data adaptation were tested against results from a detailed model. This shows that using commonly available data and unifying and simplifying the input data is tolerable up to a certain level. The simplification of the model has greater effects on the evaluated water balance components than the data adaptation. Because this simplification was necessary for large-scale use, we suggest that, for reasons of comparability, simpler models should always be applied with uniform data bases for large regions, though these should only be moderately simplified. Further, we recommend using these simplified models only for large-scale comparisons and using more specific, detailed models for investigations on smaller scales.
Leaf wax n-alkane H-2 values carry important information about environmental and ecophysiological processes in plants. However, the physiological and biochemical drivers that shape leaf wax n-alkane H-2 values are not completely understood. It is particularly unclear why n-alkanes in grasses are typically H-2-depleted compared with plants from other taxonomic groups such as dicotyledonous plants and why C3 grasses are H-2-depleted compared with C4 grasses. To resolve these uncertainties, we quantified the effects of leaf water evaporative H-2-enrichment and biosynthetic hydrogen isotope fractionation on n-alkane H-2 values for a range of C3 and C4 grasses grown in climate-controlled chambers. We found that only a fraction of leaf water evaporative H-2-enrichment is imprinted on the leaf wax n-alkane H-2 values in grasses. This is interesting, as previous studies have shown in dicotyledonous plants a nearly complete transfer of this H-2-enrichment to the n-alkane H-2 values. We thus infer that the typically observed H-2-depletion of n-alkanes in grasses (as opposed to dicots) is because only a fraction of the leaf water evaporative H-2-enrichment is imprinted on the H-2 values. Our experiments also show that differences in n-alkane H-2 values between C3 and C4 grasses are largely the result of systematic differences in biosynthetic fractionation between these two plant groups, which was on average -198 and-159 parts per thousand for C3 and C4 grasses, respectively. We present novel and exciting data on how leaf wax n-alkane 2H values from grasses are affected by plant physiological (leaf water evaporative 2H-enrichment) or biochemical processes (biosynthetic hydrogen isotope fractionation). These results are very interesting because they shed new light on how naturally observed differences between 2H values from C3 and C4 grasses and dicots can be explained by systematic differences in the biosynthesis of n-alkanes between these plant group (i.e. largely driven by NADPH origins).
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.
Elevation-dependent changes in n-alkane delta D and soil GDGTs across the South Central Andes
(2016)
Surface uplift of large plateaus may significantly influence regional climate and more specifically precipitation patterns and temperature, sometimes complicating paleoaltimetry interpretations. Thus, understanding the topographic evolution of tectonically active mountain belts benefits from continued development of reliable proxies to reduce uncertainties in paleoaltimetry reconstructions. Lipid biomarker-based proxies provide a novel approach to stable isotope paleoaltimetry and complement authigenic or pedogenic mineral proxy materials, in particular outside semi-arid climate zones where soil carbonates are not abundant but (soil) organic matter has a high preservation potential. Here we present delta D values of soil-derived n-alkanes and mean annual air temperature (MAT) estimates based on branched glycerol dialkyl glycerol tetraether (brGDGT) distributions to assess their potential for paleoelevation reconstructions in the southern central Andes. We analyzed soil samples across two environmental and hydrological gradients that include a hillslope (26-28 degrees S) and a valley (22-24 degrees S) transect on the windward flanks of Central Andean Eastern Cordillera in NW Argentina. Our results show that present-day n-alkane delta D values and brGDGT-based MAT estimates are both linearly related with elevation and in good agreement with present-day climate conditions. Soil n-alkanes show a delta D lapse rate (A(delta D)) of -1.64 parts per thousand/100 m (R-2 = 0.91, p < 0.01) at the hillslope transect, within the range of delta D lapse rates from precipitation and surface waters in other tropical regions in the Andes like the Eastern Cordillera in Colombia and Bolivia and the Equatorial and Peruvian Andes. BrGDGT-derived soil temperatures are similar to monitored winter temperatures in the region and show a lapse rate of Delta T = -0.51 degrees C/100 m (R-2 = 0.91, p < 0.01), comparable with lapse rates from in situ soil temperature measurements, satellite derived land-surface temperatures at this transect, and weather stations from the Eastern Cordillera at similar latitude. As a result of an increasing leeward sampling position along the valley transect lapse rates are biased towards lower values and display higher scatter (Delta(delta D) = -0.9 parts per thousand/100 m, R-2 = 0.76, p < 0.01 and Delta T = -0.19 degrees C/100 m, R-2 = 0.48, p < 0.05). Despite this higher complexity, they are in line with lapse rates from stream-water samples and in situ soil temperature measurements along the same transect. Our results demonstrate that both soil n-alkane delta D values and MAT reconstructions based on brGDGTs distributions from the hillslope transect (Delta(delta D) = -1.64 parts per thousand/100 m, R-2 = 0.91, p < 0.01 and Delta T = -0.51 degrees C/100 m, R-2 = 0.91, p < 0.01) track the direct effects of orography on precipitation and temperature and hence the combined effects of local and regional hydrology as well as elevation. (C) 2016 Elsevier B.V. All rights reserved.
EnGeoMAP 2.0-Automated Hyperspectral Mineral Identification for the German EnMAP Space Mission
(2016)
Algorithms for a rapid analysis of hyperspectral data are becoming more and more important with planned next generation spaceborne hyperspectral missions such as the Environmental Mapping and Analysis Program (EnMAP) and the Japanese Hyperspectral Imager Suite (HISUI), together with an ever growing pool of hyperspectral airborne data. The here presented EnGeoMAP 2.0 algorithm is an automated system for material characterization from imaging spectroscopy data, which builds on the theoretical framework of the Tetracorder and MICA (Material Identification and Characterization Algorithm) of the United States Geological Survey and of EnGeoMAP 1.0 from 2013. EnGeoMAP 2.0 includes automated absorption feature extraction, spatio-spectral gradient calculation and mineral anomaly detection. The usage of EnGeoMAP 2.0 is demonstrated at the mineral deposit sites of Rodalquilar (SE-Spain) and Haib River (S-Namibia) using HyMAP and simulated EnMAP data. Results from Hyperion data are presented as supplementary information.
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.
Arctic coastal infrastructure and cultural and archeological sites are increasingly vulnerable to erosion and flooding due to amplified warming of the Arctic, sea level rise, lengthening of open water periods, and a predicted increase in frequency of major storms. Mitigating these hazards necessitates decision-making tools at an appropriate scale. The objectives of this paper are to provide such a tool by assessing potential erosion and flood hazards at Herschel Island, a UNESCO World Heritage candidate site. This study focused on Simpson Point and the adjacent coastal sections because of their archeological, historical, and cultural significance. Shoreline movement was analyzed using the Digital Shoreline Analysis System (DSAS) after digitizing shorelines from 1952, 1970, 2000, and 2011. For purposes of this analysis, the coast was divided in seven coastal reaches (CRs) reflecting different morphologies and/or exposures. Using linear regression rates obtained from these data, projections of shoreline position were made for 20 and 50 years into the future. Flood hazard was assessed using a least cost path analysis based on a high-resolution light detection and ranging (LiDAR) dataset and current Intergovernmental Panel on Climate Change sea level estimates. Widespread erosion characterizes the study area. The rate of shoreline movement in different periods of the study ranges from -5.5 to 2.7 mI double dagger a(-1) (mean -0.6 mI double dagger a(-1)). Mean coastal retreat decreased from -0.6 mI double dagger a(-1) to -0.5 mI double dagger a(-1), for 1952-1970 and 1970-2000, respectively, and increased to -1.3 mI double dagger a(-1) in the period 2000-2011. Ice-rich coastal sections most exposed to wave attack exhibited the highest rates of coastal retreat. The geohazard map combines shoreline projections and flood hazard analyses to show that most of the spit area has extreme or very high flood hazard potential, and some buildings are vulnerable to coastal erosion. This study demonstrates that transgressive forcing may provide ample sediment for the expansion of depositional landforms, while growing more susceptible to overwash and flooding.
The role of soil erosion in terrestrial carbon (C) sequestration and release remains one of the most important uncertainties in our attempts to determine the potential of soils to mediate climate change. Despite its widely recognized importance for terrestrial C sequestration, to date, no Earth System Model (ESM) implements soil erosion effects on carbon cycling in sufficient detail. So far, available studies have mostly investigated the magnitude of erosional C transport and in-situ measurements of vertical C fluxes on the catchment or regional scale. Recognizing the need to adequately represent C erosion processes and controls in ESMs, we provide a comprehensive cross-disciplinary review on lateral C redistribution in the landscape and discuss the implications for bio-geochemical cycling of carbon. We present current knowledge on the role of erosional C distribution in controlling the stabilization and release of C in soils, taking into consideration the important geomorphic, ecological, hydrologic, pedologic and micro-climatic processes and controls that affect soil organic carbon (SOC) stock, fluxes, and persistence in dynamic landscapes. Further, we provide an overview on latest experimental and modelling approaches that are being used to investigate the role of erosion in the carbon cycle. Finally, to advance our understanding of the role of soil redistribution in biogeochemical cycles of essential elements, we discuss the most promising topics for future research in this field. (C) 2015 Elsevier B.V. All rights reserved.
Flood damage has increased significantly and is expected to rise further in many parts of the world. For assessing potential changes in flood risk, this paper presents an integrated model chain quantifying flood hazards and losses while considering climate and land use changes. In the case study region, risk estimates for the present and the near future illustrate that changes in flood risk by 2030 are relatively low compared to historic periods. While the impact of climate change on the flood hazard and risk by 2030 is slight or negligible, strong urbanisation associated with economic growth contributes to a remarkable increase in flood risk. Therefore, it is recommended to frequently consider land use scenarios and economic developments when assessing future flood risks. Further, an adapted and sustainable risk management is necessary to encounter rising flood losses, in which non-structural measures are becoming more and more important. The case study demonstrates that adaptation by non-structural measures such as stricter land use regulations or enhancement of private precaution is capable of reducing flood risk by around 30 %. Ignoring flood risks, in contrast, always leads to further increasing losses-with our assumptions by 17 %. These findings underline that private precaution and land use regulation could be taken into account as low cost adaptation strategies to global climate change in many flood prone areas. Since such measures reduce flood risk regardless of climate or land use changes, they can also be recommended as no-regret measures.
Adequate irrigation inputs are essential for the application of hydrological models in irrigated catchments, but reliable data on both the amount and the frequency of irrigation applications are often missing at an appropriate spatial scale. In this paper, we demonstrate and test approaches to estimate irrigation inputs for distributed hydrological modelling. In this context, the Soil and Water Assessment Tool was applied to simulate water balances for an irrigated catchment in southeast Australia during the period 2008–2010. Two methods for estimating irrigation inputs were tested. One method was based on a fixed irrigation application rate, whereas the other one had variable irrigation rates depending on season and the irrigated crop. These two approaches were also compared with the ‘auto-irrigation’ method within the Soil and Water Assessment Tool model. The method with variable irrigation rates resulted in the most reasonable interpretation of the readily available irrigation data, consistent estimates of irrigation runoff coefficients throughout the year and the best fit to observed data on both drain flows at the catchment outlet and spatial evapotranspiration patterns. We also found that the different irrigation inputs significantly affected simulated water balances, in particular deep percolation under relatively dry climatic conditions. All these results suggest that it is possible to infer irrigation inputs from readily available data and local knowledge, adequate for hydrological modelling in irrigated catchments. Our study also demonstrates that, in order to predict reliable water balances in irrigated catchments, an accurate knowledge of irrigation scheduling and irrigation runoff is required. Copyright © 2015 John Wiley & Sons, Ltd.
estimating mean throughfall
(2016)
The selection of an appropriate spatial extent of a sampling plot is one among several important decisions involved in planning a throughfall sampling scheme. In fact, the choice of the extent may determine whether or not a study can adequately characterize the hydrological fluxes of the studied ecosystem. Previous attempts to optimize throughfall sampling schemes focused on the selection of an appropriate sample size, support, and sampling design, while comparatively little attention has been given to the role of the extent. In this contribution, we investigated the influence of the extent on the representativeness of mean throughfall estimates for three forest ecosystems of varying stand structure. Our study is based on virtual sampling of simulated throughfall fields. We derived these fields from throughfall data sampled in a simply structured forest (young tropical forest) and two heterogeneous forests (old tropical forest, unmanaged mixed European beech forest). We then sampled the simulated throughfall fields with three common extents and various sample sizes for a range of events and for accumulated data. Our findings suggest that the size of the study area should be carefully adapted to the complexity of the system under study and to the required temporal resolution of the throughfall data (i.e. event-based versus accumulated). Generally, event-based sampling in complex structured forests (conditions that favor comparatively long autocorrelations in throughfall) requires the largest extents. For event-based sampling, the choice of an appropriate extent can be as important as using an adequate sample size. (C) 2016 Elsevier B.V. All rights reserved.
This case study evaluates the suitability of radar-based quantitative precipitation estimates (QPEs) for the simulation of streamflow in the Marikina River Basin (MRB), the Philippines. Hourly radar-based QPEs were produced from reflectivity that had been observed by an S-band radar located about 90 km from the MRB. Radar data processing and precipitation estimation were carried out using the open source library wradlib. To assess the added value of the radar-based QPE, we used spatially interpolated rain gauge observations (gauge-only (GO) product) as a benchmark. Rain gauge observations were also used to quantify rainfall estimation errors at the point scale. At the point scale, the radar-based QPE outperformed the GO product in 2012, while for 2013, the performance was similar. For both periods, estimation errors substantially increased from daily to the hourly accumulation intervals. Despite this fact, both rainfall estimation methods allowed for a good representation of observed streamflow when used to force a hydrological simulation model of the MRB. Furthermore, the results of the hydrological simulation were consistent with rainfall verification at the point scale: the radar-based QPE performed better than the GO product in 2012, and equivalently in 2013. Altogether, we could demonstrate that, in terms of streamflow simulation, the radar-based QPE can perform as good as or even better than the GO product - even for a basin such as the MRB which has a comparatively dense rain gauge network. This suggests good prospects for using radar-based QPE to simulate and forecast streamflow in other parts of the Philippines where rain gauge networks are not as dense.
Event coincidence analysis for quantifying statistical interrelationships between event time series
(2016)
Studying event time series is a powerful approach for analyzing the dynamics of complex dynamical systems in many fields of science. In this paper, we describe the method of event coincidence analysis to provide a framework for quantifying the strength, directionality and time lag of statistical interrelationships between event series. Event coincidence analysis allows to formulate and test null hypotheses on the origin of the observed interrelationships including tests based on Poisson processes or, more generally, stochastic point processes with a prescribed inter-event time distribution and other higher-order properties. Applying the framework to country-level observational data yields evidence that flood events have acted as triggers of epidemic outbreaks globally since the 1950s. Facing projected future changes in the statistics of climatic extreme events, statistical techniques such as event coincidence analysis will be relevant for investigating the impacts of anthropogenic climate change on human societies and ecosystems worldwide.
There is increasing evidence for recent changes in the intensity and frequency of heavy precipitation and in the number of days with snow cover in many parts of Norway. The question arises as to whether these changes are also discernable with respect to their impacts on the magnitude and frequency of flooding and on the processes producing high flows. In this study, we tested up to 211 catchments for trends in peak flow discharge series by applying the Mann-Kendall test and Poisson regression for three different time periods (1962-2012, 1972-2012, 1982-2012). Field-significance was tested using a bootstrap approach. Over threshold discharge events were classified into rainfall vs. snowmelt dominated floods, based on a simple water balance approach utilizing a nationwide 1 x 1 km(2) gridded data set with daily observed rainfall and simulated snowmelt data. Results suggest that trends in flood frequency are more pronounced than trends in flood magnitude and are more spatially consistent with observed changes in the hydrometeorological drivers. Increasing flood frequencies in southern and western Norway are mainly due to positive trends in the frequency of rainfall dominated events, while decreasing flood frequencies in northern Norway are mainly the result of negative trends in the frequency of snowmelt dominated floods. Negative trends in flood magnitude are found more often than positive trends, and the regional patterns of significant trends reflect differences in the flood generating processes (FGPs). The results illustrate the benefit of distinguishing FGPs rather than simply applying seasonal analyses. The results further suggest that rainfall has generally gained an increasing importance for the generation of floods in Norway, while the role of snowmelt has been decreasing and the timing of snowmelt dominated floods has become earlier. (C) 2016 Elsevier B.V. All rights reserved.
The recognition and understanding of vegetated habitats in the fossil record are of crucial importance in order to investigate paleoecological responses and indirectly infer climate and sea-level changes. However, the low preservation potential of plants and macroalgae hampers a direct identification of these environments in the geological past. Here we present sedimentological and paleontological evidences as tool to identify the presence of different seagrass-vegetated environments in the shallow marine settings of the lower Eocene jafnayn platform of Oman and their responses to paleoenvironmental changes. The studied lower Eocene deposits consist of well bedded, nodular pacicstones dominated by encrusting acervulinid and alveolinid foraminifera passing upward to an alternance of packstones with echinoids and quartz grains and grainstones rich in Orbitolites, smaller miliolid foraminifera and quartz grains. The presence of seagrass is inferred by the occurrence of encrusting acervulinids and soritid Orbitolites, as well as by their test morphologies together with further sedimentological criteria. The clear shift observed in the faunal assemblages and sedimentary features may be related to a major reorganization of the carbonate system passing from a carbonate platform to a ramp-like platform with increased terrigenous sedimentation. Heterotroph tubular acervulinids and oligotroph alveolinids of the carbonate platform were replaced upward by more heterotroph organisms such as large, discoidal Orbitolites and smaller miliolids, most likely due to enhanced nutrient levels which would have led to a change of phytal substrate, from cylindrical-leaf dominated grasses into flat-leafed ones. (C) 2016 Elsevier B.V. All rights reserved.
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.
We present new experimental data of the low-temperature metastable region of liquid water derived from high-density synthetic fluid inclusions (996-916 kg m(-3)) in quartz. Microthermometric measurements include: (i) prograde (upon heating) and retrograde (upon cooling) liquid-vapour homogenisation. We used single ultrashort laser pulses to stimulate vapour bubble nucleation in initially monophase liquid inclusions. Water densities were calculated based on prograde homogenisation temperatures using the IAPWS-95 formulation. We found retrograde liquid-vapour homogenisation temperatures in excellent agreement with IAPWS-95. (ii) Retrograde ice nucleation. Raman spectroscopy was used to determine the nucleation of ice in the absence of the vapour bubble. Our ice nucleation data in the doubly metastable region are inconsistent with the low-temperature trend of the spinodal predicted by IAPWS-95, as liquid water with a density of 921 kg m(-3) remains in a homogeneous state during cooling down to a temperature of -30.5 degrees C, where it is transformed into ice whose density corresponds to zero pressure. (iii) Ice melting. Ice melting temperatures of up to 6.8 degrees C were measured in the absence of the vapour bubble, i.e. in the negative pressure region. (iv) Spontaneous retrograde and, for the first time, prograde vapour bubble nucleation. Prograde bubble nucleation occurred upon heating at temperatures above ice melting. The occurrence of prograde and retrograde vapour bubble nucleation in the same inclusions indicates a maximum of the bubble nucleation curve in the R-T plane at around 40 degrees C. The new experimental data represent valuable benchmarks to evaluate and further improve theoretical models describing the p-V-T properties of metastable water in the low-temperature region.
We present new experimental data of the low-temperature metastable region of liquid water derived from high-density synthetic fluid inclusions (996–916 kg m−3) in quartz. Microthermometric measurements include: (i) prograde (upon heating) and retrograde (upon cooling) liquid–vapour homogenisation. We used single ultrashort laser pulses to stimulate vapour bubble nucleation in initially monophase liquid inclusions. Water densities were calculated based on prograde homogenisation temperatures using the IAPWS-95 formulation. We found retrograde liquid–vapour homogenisation temperatures in excellent agreement with IAPWS-95. (ii) Retrograde ice nucleation. Raman spectroscopy was used to determine the nucleation of ice in the absence of the vapour bubble. Our ice nucleation data in the doubly metastable region are inconsistent with the low-temperature trend of the spinodal predicted by IAPWS-95, as liquid water with a density of 921 kg m−3 remains in a homogeneous state during cooling down to a temperature of −30.5 °C, where it is transformed into ice whose density corresponds to zero pressure. (iii) Ice melting. Ice melting temperatures of up to 6.8 °C were measured in the absence of the vapour bubble, i.e. in the negative pressure region. (iv) Spontaneous retrograde and, for the first time, prograde vapour bubble nucleation. Prograde bubble nucleation occurred upon heating at temperatures above ice melting. The occurrence of prograde and retrograde vapour bubble nucleation in the same inclusions indicates a maximum of the bubble nucleation curve in the ϱ–T plane at around 40 °C. The new experimental data represent valuable benchmarks to evaluate and further improve theoretical models describing the p–V–T properties of metastable water in the low-temperature region.
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.
Knowledge of the quality factor of near-surface materials is of fundamental interest in various applications. Attenuation can be very strong close to the surface and thus needs to be properly assessed. In recent years, several researchers have studied the retrieval of attenuation coefficients from the cross correlation of ambient seismic noise. Yet, the determination of exact amplitude information from noise-correlation functions is, in contrast to the extraction of traveltimes, not trivial. Most of the studies estimated attenuation coefficients on the regional scale and within the microseism band. In this paper, we investigate the possibility to derive attenuation coefficients from seismic noise at much shallower depths and higher frequencies (> 1 Hz). The Euroseistest area in northern Greece offers ideal conditions to study quality factor retrieval from ambient noise for different rock types. Correlations are computed between the stations of a small scale array experiment (station spacings < 2 km) that was carried out in the Euroseistest area in 2011. We employ the correlation of the coda of the correlation (C-3) method instead of simple cross correlations to mitigate the effect of uneven noise source distributions on the correlation amplitude. Transient removal and temporal flattening are applied instead of 1-bit normalization in order to retain relative amplitudes. The C-3 method leads to improved correlation results (higher signal-to-noise ratio and improved time symmetry) compared to simple cross correlations. The C-3 functions are rotated from the ZNE to the ZRT system and we focus on Love wave arrivals on the transverse component and on Love wave quality factors Q(L). The analysis is performed for selected stations being either situated on soft soil or on weathered rock. Phase slowness is extracted using a slant-stack method. Attenuation parameters are inferred by inspecting the relative amplitude decay of Love waves with increasing interstation distance. We observe that the attenuation coefficient gamma and Q(L) can be reliably extracted for stations situated on soft soil whereas the derivation of attenuation parameters is more problematic for stations that are located on weathered rock. The results are in acceptable conformance with theoretical Love wave attenuation curves that were computed using 1-D shear wave velocity and quality factor profiles from the Euroseistest area.
Injection of fluids into deep saline aquifers causes a pore pressure increase in the storage formation, and thus displacement of resident brine. Via hydraulically conductive faults, brine may migrate upwards into shallower aquifers and lead to unwanted salinisation of potable groundwater resources. In the present study, we investigated different scenarios for a potential storage site in the Northeast German Basin using a three-dimensional (3-D) regional-scale model that includes four major fault zones. The focus was on assessing the impact of fault length and the effect of a secondary reservoir above the storage formation, as well as model boundary conditions and initial salinity distribution on the potential salinisation of shallow groundwater resources. We employed numerical simulations of brine injection as a representative fluid. Our simulation results demonstrate that the lateral model boundary settings and the effective fault damage zone volume have the greatest influence on pressure build-up and development within the reservoir, and thus intensity and duration of fluid flow through the faults. Higher vertical pressure gradients for short fault segments or a small effective fault damage zone volume result in the highest salinisation potential due to a larger vertical fault height affected by fluid displacement. Consequently, it has a strong impact on the degree of shallow aquifer salinisation, whether a gradient in salinity exists or the saltwater-freshwater interface lies below the fluid displacement depth in the faults. A small effective fault damage zone volume or low fault permeability further extend the duration of fluid flow, which can persist for several tens to hundreds of years, if the reservoir is laterally confined. Laterally open reservoir boundaries, large effective fault damage zone volumes and intermediate reservoirs significantly reduce vertical brine migration and the potential of freshwater salinisation because the origin depth of displaced brine is located only a few decametres below the shallow aquifer in maximum. The present study demonstrates that the existence of hydraulically conductive faults is not necessarily an exclusion criterion for potential injection sites, because salinisation of shallower aquifers strongly depends on initial salinity distribution, location of hydraulically conductive faults and their effective damage zone volumes as well as geological boundary conditions.
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.
Freshwater ostracods (Crustacea, Ostracoda) are valuable biological indicators. In Arctic environments, their habitat conditions are barely known and the abundance and diversity of ostracods is documented only in scattered records with incomplete ecological characterization. To determine the taxonomic range of ostracod assemblages and their habitat conditions in polygon ponds in the Indigirka Lowland, north-east Siberia, we collected more than 100 living ostracod individuals per site with a plankton net (mesh size 65 mm) and an exhaustor system from 27 water bodies and studied them in the context of substrate and hydrochemical data. During the summer of 2011, a single pond site and its ostracod population was selected for special study. This first record of the ostracod fauna in the Indigirka Lowland comprises eight species and three additional taxa. Fabaeformiscandona krochini and F. groenlandica were documented for the first time in continental Siberia. Repeated sampling of a low-centre polygon pond yielded insights into the population dynamics of F. pedata. We identified air temperature and precipitation as the main external drivers of water temperatures, water levels, ion concentrations and water stable isotope composition on diurnal and seasonal scales.
Peatlands are scarce and threatened ecosystems in the semiarid region of Kyrgyzstan. Knowledge about the Kyrgyz peatlands is still poor and, especially, their genesis has hardly been investigated so far. Typically, the peatland substrates are characterised by the admixture of silt-sized particles in various quantities. In this work we report the abiotic properties and genesis of three peatlands within different altitudinal zones in southern Kyrgyzstan. We surveyed the stratification of the peatlands and their water chemistry. In addition, we investigated whether the silt found in the peatland substrates was deposited by wind, rivers or springs. The mineral constituents of the peatland substrates were analysed for particle size distribution and their elemental composition was compared with that of nearby loess, river and spring sediments using the immobile trace element titanium. One peatland shows a high abundance of different peatland substrates, indicating a frequent change of ecological conditions in the past. All three peatlands are fed by groundwater. Overgrazing and trampling by cattle has led to recent degradation of the upper peat layer. The resulting compaction of the peats prevents water from seeping into the substrates of the peatlands and subsequently changes their hydrology. Our results indicate that both wind and rivers have deposited silt in the peatlands, depending on their positions in the relief. Silts may also have been relocated by springs within the peatlands.
New age data have been obtained to time constrain the recent Quaternary volcanism of El Hierro (Canary Islands) and to estimate its recurrence rate. We have carried out Ar-40/Ar-39 geochronology on samples spanning the entire volcanostratigraphic sequence of the island and C-14 geochronology on the most recent eruption on the northeast rift of the island: 2280 +/- 30 yr BP. We combine the new absolute data with a revision of published ages onshore, some of which were identified through geomorphological criteria (relative data). We present a revised and updated chronology of volcanism for the last 33 ka that we use to estimate the maximum eruptive recurrence of the island. The number of events per year determined is 9.7 x 10(-4) for the emerged part of the island, which means that, as a minimum, one eruption has occurred approximately every 1000 years. This highlights the need of more geochronological data to better constrain the eruptive recurrence of El Hierro. (C) 2015 Elsevier Ltd. All rights reserved.
Lake Towuti is a tectonic basin, surrounded by ultramafic rocks. Lateritic soils form through weathering and deliver abundant iron (oxy)hydroxides but very little sulfate to the lake and its sediment. To characterize the sediment biogeochemistry, we collected cores at three sites with increasing water depth and decreasing bottom water oxygen concentrations. Microbial cell densities were highest at the shallow site a feature we attribute to the availability of labile organic matter (OM) and the higher abundance of electron acceptors due to oxic bottom water conditions. At the two other sites, OM degradation and reduction processes below the oxycline led to partial electron acceptor depletion. Genetic information preserved in the sediment as extracellular DNA (eDNA) provided information on aerobic and anaerobic heterotrophs related to Nitrospirae. Chloroflexi, and Therrnoplasmatales. These taxa apparently played a significant role in the degradation of sinking OM. However, eDNA concentrations rapidly decreased with core depth. Despite very low sulfate concentrations, sulfate-reducing bacteria were present and viable in sediments at all three sites, as confirmed by measurement of potential sulfate reduction rates. Microbial community fingerprinting supported the presence of taxa related to Deltaproteobacteria and Firmicutes with demonstrated capacity for iron and sulfate reduction. Concomitantly, sequences of Ruminococcaceae, Clostridiales, and Methanornicrobiales indicated potential for fermentative hydrogen and methane production. Such first insights into ferruginous sediments showed that microbial populations perform successive metabolisms related to sulfur, iron, and methane. In theory, iron reduction could reoxidize reduced sulfur compounds and desorb OM from iron minerals to allow remineralization to methane. Overall, we found that biogeochemical processes in the sediments can be linked to redox differences in the bottom waters of the three sites, like oxidant concentrations and the supply of labile OM. At the scale of the lacustrine record, our geomicrobiological study should provide a means to link the extant subsurface biosphere to past environments.