The precise determination of the stable C and O isotope fractionation between water and calcite (CC) and water and aragonite (AR) is of special interest for climate reconstructions, e.g. paleotemperatures. Previous studies reported results from both laboratory and field experiments, but their results are only partly consistent. Here we present C and O isotope data of a stalagmite from the Swiss Alps, which shows CC-AR transitions along individual growth layers. Using detailed analyses both laterally and perpendicular to such layers we examined the difference in the C and O isotope fractionation factor of the HCO3- - CC and the HCO3- - AR system. For O this difference is similar to the water-CC and water-AR offset provided in experimental studies. The O isotope fractionation difference in the water-CC and water-AR system is comparable to those determined in laboratory studies but shows a statistically significant correlation with the CaCO3 precipitation rate. For C we found a fractionation difference, which is independent of CaCO3 precipitation rate and with slightly smaller values for the fractionation offset between HCO3- - CC and HCO3- - AR compared to literature values. However, we also found an unexpected decrease in delta C-13 along growth layers, which contradicts the widely used concept of Rayleigh fractionation during CO2 degassing and CaCO3 precipitation. The results of this study can be used e.g., to correct stable isotope time series of stalagmites showing CC-AR transitions along their growth axes. (C) 2018 Elsevier Ltd. All rights reserved.

Time series of groundwater and stream water quality often exhibit substantial temporal and spatial variability, whereas typical existing monitoring data sets, e.g. from environmental agencies, are usually characterized by relatively low sampling frequency and irregular sampling in space and/or time. This complicates the differentiation between anthropogenic influence and natural variability as well as the detection of changes in water quality which indicate changes in single drivers. We suggest the new term "dominant changes" for changes in multivariate water quality data which concern (1) multiple variables, (2) multiple sites and (3) long-term patterns and present an exploratory framework for the detection of such dominant changes in data sets with irregular sampling in space and time. Firstly, a non-linear dimension-reduction technique was used to summarize the dominant spatiotemporal dynamics in the multivariate water quality data set in a few components. Those were used to derive hypotheses on the dominant drivers influencing water quality. Secondly, different sampling sites were compared with respect to median component values. Thirdly, time series of the components at single sites were analysed for long-term patterns. We tested the approach with a joint stream water and groundwater data set quality consisting of 1572 samples, each comprising sixteen variables, sampled with a spatially and temporally irregular sampling scheme at 29 sites in northeast Germany from 1998 to 2009. The first four components were interpreted as (1) an agriculturally induced enhancement of the natural background level of solute concentration, (2) a redox sequence from reducing conditions in deep groundwater to post-oxic conditions in shallow groundwater and oxic conditions in stream water, (3) a mixing ratio of deep and shallow groundwater to the streamflow and (4) sporadic events of slurry application in the agricultural practice. Dominant changes were observed for the first two components. The changing intensity of the first component was interpreted as response to the temporal variability of the thickness of the unsaturated zone. A steady increase in the second component at most stream water sites pointed towards progressing depletion of the denitrification capacity of the deep aquifer.

The AlpArray programme is a multinational, European consortium to advance our understanding of orogenesis and its relationship to mantle dynamics, plate reorganizations, surface processes and seismic hazard in the Alps-Apennines-Carpathians-Dinarides orogenic system. The AlpArray Seismic Network has been deployed with contributions from 36 institutions from 11 countries to map physical properties of the lithosphere and asthenosphere in 3D and thus to obtain new, high-resolution geophysical images of structures from the surface down to the base of the mantle transition zone. With over 600 broadband stations operated for 2 years, this seismic experiment is one of the largest simultaneously operated seismological networks in the academic domain, employing hexagonal coverage with station spacing at less than 52 km. This dense and regularly spaced experiment is made possible by the coordinated coeval deployment of temporary stations from numerous national pools, including ocean-bottom seismometers, which were funded by different national agencies. They combine with permanent networks, which also required the cooperation of many different operators. Together these stations ultimately fill coverage gaps. Following a short overview of previous large-scale seismological experiments in the Alpine region, we here present the goals, construction, deployment, characteristics and data management of the AlpArray Seismic Network, which will provide data that is expected to be unprecedented in quality to image the complex Alpine mountains at depth.

As flood impacts are increasing in large parts of the world, understanding the primary drivers of changes in risk is essential for effective adaptation. To gain more knowledge on the basis of empirical case studies, we analyze eight paired floods, that is, consecutive flood events that occurred in the same region, with the second flood causing significantly lower damage. These success stories of risk reduction were selected across different socioeconomic and hydro-climatic contexts. The potential of societies to adapt is uncovered by describing triggered societal changes, as well as formal measures and spontaneous processes that reduced flood risk. This novel approach has the potential to build the basis for an international data collection and analysis effort to better understand and attribute changes in risk due to hydrological extremes in the framework of the IAHSs Panta Rhei initiative. Across all case studies, we find that lower damage caused by the second event was mainly due to significant reductions in vulnerability, for example, via raised risk awareness, preparedness, and improvements of organizational emergency management. Thus, vulnerability reduction plays an essential role for successful adaptation. Our work shows that there is a high potential to adapt, but there remains the challenge to stimulate measures that reduce vulnerability and risk in periods in which extreme events do not occur.

The shape of trait distributions may inform about the selective forces that structure ecological communities. Here, we present a new moment-based approach to classify the shape of observed biomass-weighted trait distributions into normal, peaked, skewed, or bimodal that facilitates spatio-temporal and cross-system comparisons. Our observed phytoplankton trait distributions exhibited substantial variance and were mostly skewed or bimodal rather than normal. Additionally, mean, variance, skewness und kurtosis were strongly correlated. This is in conflict with trait-based aggregate models that often assume normally distributed trait values and small variances. Given these discrepancies between our data and general model assumptions we used the observed trait distributions to test how well different aggregate models with first- or second-order approximations and different types of moment closure predict the biomass, mean trait, and trait variance dynamics using weakly or moderately nonlinear fitness functions. For weakly non-linear fitness functions aggregate models with a second-order approximation and a data-based moment closure that relied on the observed correlations between skewness and mean, and kurtosis and variance predicted biomass and often also mean trait changes fairly well and better than models with first-order approximations or a normal-based moment closure. In contrast, none of the models reflected the changes of the trait variances reliably. Aggregate model performance was often also poor for moderately nonlinear fitness functions. This questions a general applicability of the normal-based approach, in particular for predicting variance dynamics determining the speed of trait changes and maintenance of biodiversity. We evaluate in detail how and why better approximations can be obtained.

Carbonatite magmatism is a highly efficient transport mechanism from Earth’s mantle to the crust, thus providing insights into the chemistry and dynamics of the Earth’s mantle. One evolving and promising tool for tracing magma interaction are stable iron isotopes, particularly because iron isotope fractionation is controlled by oxidation state and bonding environment. Meanwhile, a large data set on iron isotope fractionation in igneous rocks exists comprising bulk rock compositions and fractionation between mineral groups. Iron isotope data from natural carbonatite rocks are extremely light and of remarkably high variability. This resembles iron isotope data from mantle xenoliths, which are characterized by a variability in δ56Fe spanning three times the range found in basalts, and by the extremely light values of some whole rock samples, reaching δ56Fe as low as -0.69 ‰ in a spinel lherzolite. Cause to this large range of variations may be metasomatic processes, involving metasomatic agents like volatile bearing high-alkaline silicate melts or carbonate melts. The expected effects of metasomatism on iron isotope fractionation vary with parameters like melt/rock-ratio, reaction time, and the nature of metasomatic agents and mineral reactions involved. An alternative or additional way to enrich light isotopes in the mantle could be multiple phases of melt extraction. To interpret the existing data sets more knowledge on iron isotope fractionation factors is needed. To investigate the behavior of iron isotopes in the carbonatite systems, kinetic and equilibration experiments in natro-carbonatite systems between immiscible silicate and carbonate melts were performed in an internally heated gas pressure vessel at intrinsic redox conditions at temperatures between 900 and 1200 °C and pressures of 0.5 and 0.7 GPa. The iron isotope compositions of coexisting silicate melt and carbonate melt were analyzed by solution MC-ICP-MS. The kinetic experiments employing a Fe-58 spiked starting material show that isotopic equilibrium is obtained after 48 hours. The experimental studies of equilibrium iron isotope fractionation between immiscible silicate and carbonate melts have shown that light isotopes are enriched in the carbonatite melt. The highest Δ56Fesil.m.-carb.melt (mean) of 0.13 ‰ was determined in a system with a strongly peralkaline silicate melt composition (ASI ≥ 0.21, Na/Al ≤ 2.7). In three systems with extremely peralkaline silicate melt compositions (ASI between 0.11 and 0.14) iron isotope fractionation could analytically not be resolved. The lowest Δ56Fesil.m.-carb.melt (mean) of 0.02 ‰ was determined in a system with an extremely peralkaline silicate melt composition (ASI ≤ 0.11 , Na/Al ≥ 6.1). The observed iron isotope fractionation is most likely governed by the redox conditions of the system. Yet, in the systems, where no fractionation occurred, structural changes induced by compositional changes possibly overrule the influence of redox conditions. This interpretation implicates, that the iron isotope system holds the potential to be useful not only for exploring redox conditions in magmatic systems, but also for discovering structural changes in a melt. In situ iron isotope analyses by femtosecond laser ablation coupled to MC-ICP-MS on magnetite and olivine grains were performed to reveal variations in iron isotope composition on the micro scale. The investigated sample is a melilitite bomb from the Salt Lake Crater group at Honolulu (Oahu, Hawaii), showing strong evidence for interaction with a carbonatite melt. While magnetite grains are rather homogeneous in their iron isotope compositions, olivine grains span a far larger range in iron isotope ratios. The variability of δ56Fe in magnetite is limited from - 0.17 ‰ (± 0.11 ‰, 2SE) to +0.08 ‰ (± 0.09 ‰, 2SE). δ56Fe in olivine range from -0.66‰ (± 0.11 ‰, 2SE) to +0.10 ‰ (± 0.13 ‰, 2SE). Olivine and magnetite grains hold different informations regarding kinetic and equilibrium fractionation due to their different Fe diffusion coefficients. The observations made in the experiments and in the in situ iron isotope analyses suggest that the extremely light iron isotope signatures found in carbonatites are generated by several steps of isotope fractionation during carbonatite genesis. These may involve equilibrium and kinetic fractionation. Since iron isotopic signatures in natural systems are generated by a combination of multiple factors (pressure, temperature, redox conditions, phase composition and structure, time scale), multi tracer approaches are needed to explain signatures found in natural rocks.

Great megathrust earthquakes arise from the sudden release of energy accumulated during centuries of interseismic plate convergence. The moment deficit (energy available for future earthquakes) is commonly inferred by integrating the rate of interseismic plate locking over the time since the previous great earthquake. But accurate integration requires knowledge of how interseismic plate locking changes decades after earthquakes, measurements not available for most great earthquakes. Here we reconstruct the post-earthquake history of plate locking at Guafo Island, above the seismogenic zone of the giant 1960 (M-w = 9.5) Chile earthquake, through forward modeling of land-level changes inferred from aerial imagery (since 1974) and measured by GPS (since 1994). We find that interseismic locking increased to similar to 70% in the decade following the 1960 earthquake and then gradually to 100% by 2005. Our findings illustrate the transient evolution of plate locking in Chile, and suggest a similarly complex evolution elsewhere, with implications for the time- and magnitude-dependent probability of future events.

The Southern Central Andes (33°-36°S) are an excellent natural laboratory to study orogenic deformation processes, where boundary conditions, such as the geometry of the subducted plate, impose an important control on the evolution of the orogen. On the other hand, the South American plate presents a series of heterogeneities that additionally impart control on the mode of deformation. This thesis aims to test the control of this last factor over the construction of the Cenozoic Andean orogenic system. From the integration of surface and subsurface information in the southern area (34-36°S), the evolution of Andean deformation over the steeply dipping subduction segment was studied. A structural model was developed evaluating the stress state from the Miocene to the present-day and its influence in the migration of magmatic fluids and hydrocarbons. Based on these data, together with the data generated by other researchers in the northern zone of the study area (33-34°S), geodynamic numerical modeling was performed to test the hypothesis of the decisive role of upper-plate heterogeneities in the Andean evolution. Geodynamic codes (LAPEX-2D and ASPECT) which simulate the behavior of materials with elasto-visco-plastic rheologies under deformation, were used. The model results suggest that upper-plate contractional deformation is significantly controlled by the strength of the lithosphere, which is defined by the composition of the upper and lower crust, and by the proportion of lithospheric mantle, which in turn is determined by previous tectonic events. In addition, the previous regional tectono-magmatic events also defined the composition of the crust and its geometry, which is another factor that controls the localization of deformation. Accordingly, with more felsic lower crustal composition, the deformation follows a pure-shear mode, while more mafic compositions induce a simple-shear deformation mode. On the other hand, it was observed that initial lithospheric thickness may fundamentally control the location of deformation, with zones characterized by thin lithosphere are prone to concentrate it. Finally, it was found that an asymmetric lithosphere-astenosphere boundary resulting from corner flow in the mantle wedge of the eastward-directed subduction zone tends to generate east-vergent detachments.

We used a change point analysis on a late Pleistocene-Holocene lake-sediment record from the Chew Bahir basin in the southern Ethiopian Rift to determine the amplitude and duration of past climate transitions. The most dramatic changes occurred over 240 yr (from similar to 15,700 to 15,460 yr) during the onset of the African Humid Period (AHP), and over 990 yr (from similar to 4875 to 3885 yr) during its protracted termination. The AHP was interrupted by a distinct dry period coinciding with the high-latitude Younger Dryas stadial, which had an abrupt onset (less than similar to 100 yr) at similar to 13,260 yr and lasted until similar to 11,730 yr. Wet-dry-wet transitions prior to the AHP may reflect the high-latitude Dansgaard-Oeschger cycles, as indicated by cross-correlation of the potassium record with the NorthGRIP ice core record between similar to 45-20 ka. These findings may contribute to the debates regarding the amplitude, and duration and mechanisms of past climate transitions, and their possible influence on the development of early modern human cultures.

Cosmic-ray neutron sensing (CRNS) is a promising proximal soil sensing technique to estimate soil moisture at intermediate scale and high temporal resolution. However, the signal shows complex and non-unique response to all hydrogen pools near the land surface, providing some challenges for soil moisture estimation in practical applications. Aims of the study were 1) to assess the uncertainty of CRNS as a stand-alone approach to estimate volumetric soil moisture in cropped field 2) to identify the causes of this uncertainty 3) and possible improvements. Two experimental sites in Germany were equipped with a CRNS probe and point-scale soil moisture network. Additional monitoring activities were conducted during the crop growing season to characterize the soil-plant systems. This data is used to identify and quantify the different sources of uncertainty (factors). An uncertainty analysis, based on Monte Carlo approach, is applied to propagate these uncertainties to CRNS soil moisture estimations. In addition, a sensitivity analysis based on the Sobol’ method is performed to identify the most important factors explaining this uncertainty. Results show that CRNS soil moisture compares well to the soil moisture network when these point-scale values are weighted to account for the spatial sensitivity of the signal and other sources of hydrogen (lattice water and organic carbon) are added to the water content. However, the performance decreases when CRNS is considered as a stand-alone method to retrieve the actual (non-weighted) volumetric soil moisture. The support volume (penetration depth and radius) shows also a considerable uncertainty, especially in relatively dry soil moisture conditions. Four of the seven factors analyzed (the vertical soil moisture profile, bulk density, incoming neutron correction and the calibrated parameter N0) were found to play an important role. Among the possible improvements identified, a simple correction factor based on vertical point-scale soil moisture profiles shows to be a promising approach to account for the sensitivity of the CRNS signal to the upper soil layers.

The Siwalik sedimentary rocks of the Himalayan foreland basin preserve a record of Himalayan orogenesis, paleo-drainage evolution, and erosion. This study focuses on the still poorly studied easternmost Himalaya Siwalik record located directly downstream of the Namche Barwa syntaxis. We use luminescence, palaeomagnetism, magnetostratigraphy, and apatite fission-track dating to constrain the depositional ages of three Siwalik sequences: the Sibo outcrop (Upper Siwalik sediments at ca. 200-800 ka), the Remi section (Middle and Upper Siwalik rocks at >0.8-<8.8 +/- 2.4 Ma), and the Siang section (Middle Siwalik rocks at <9.3 +/- 1.5 to <13.5 +/- 1.5 Ma). Cretaceous-Paleogene detrital zircon and apatite U-Pb ages, characteristic of the Transhimalayan Gangdese Batholiths that crop out northwest of the syntaxis, are present throughout the Sibo, Remi, and Siang successions, confirming the existence of a Yarlung-Brahmaputra connection since at least the Late Miocene. A ca. 500 Ma zircon population increases up section, most strikingly sometime between 3.6 to 6.6 Ma, at the expense of Transhimalayan grains. We consider the ca. 500 Ma population to be derived from the Tethyan or Greater Himalaya, and we interpret the up-section increase to reflect progressive exhumation of the Namche Barwa syntaxis. Early Cretaceous zircon and apatite U-Pb ages are rare in the Sibo, Remi, and Siang successions, but abundant in modern Siang River sediments. Zircons of this age range are characteristic of the Transhimalayan Bomi-Chayu batholiths, which crop out east of the syntaxis and are eroded by the Parlung River, a modern tributary of the Siang River. We interpret the difference in relative abundance of Early Cretaceous zircons between the modern and ancient sediments to reflect capture of the Parlung by the Siang after 800 ka.

Current mitigation efforts and existing future commitments are inadequate to accomplish the Paris Agreement temperature goals. In light of this, research and debate are intensifying on the possibilities of additionally employing proposed climate geoengineering technologies, either through atmospheric carbon dioxide removal or farther-reaching interventions altering the Earth’s radiative energy budget. Although research indicates that several techniques may eventually have the physical potential to contribute to limiting climate change, all are in early stages of development, involve substantial uncertainties and risks, and raise ethical and governance dilemmas. Based on present knowledge, climate geoengineering techniques cannot be relied on to significantly contribute to meeting the Paris Agreement temperature goals.

The Red Indian Line (RIL) in central Newfoundland is the suture, where the main tract of the Iapetus Ocean was closed at similar to 452 Ma during accretion of the peri-Gondwanan Victoria arc with the composite active Laurentian margin. The protracted deformation history of this soft collision started at similar to 471 Ma with accretion of oceanic terranes to the active composite Laurentian margin. After Iapetus closure both colliding active margins were progressively deformed and metamorphosed during Silurian and Devonian (Salinic, Acadian and Neoacadian orogenic cycles). Peak conditions of the very low- to medium-grade, heterogeneously distributed metamorphism were determined by pseudosection techniques within the range of 2-7 kbar, 230-450 degrees C during increase of the metamorphic field gradient from similar to 12 degrees C/km to similar to 32 degrees C/km over time. Multiple metamorphic crystallisation stages were dated by white mica Ar-40/Ar-39 spot and plateau ages, additional Rb-Sr mineral isochrons involving white mica and one U/Pb age of titanite. All resulting ages between 439 +/- 4 Ma and 356 +/- 16 Ma postdate the closure of Iapetus. Results differ along two transects: The oldest ages of 443-421 Ma (Salinic orogenic cycle) were observed along the northern transect through the RIL zone with minimal younger overprint. Hence low temperature, intermediate to high pressure conditions (4.0-7.0 kbar, 230-340 degrees C) achieved during Taconic-Salinic underthrusting are well preserved. During Acadian dextral transpression the Taconic-Salinic structural wedge was tilted subvertically. In contrast, rocks along the southern transect through the RIL zone mainly show Acadian ages of 408-390 Ma with local preservation of older ages. Acadian deformation occurred under low temperature/low pressure conditions (similar to 250-450 degrees C, 2.5-4.6 kbar). Also Silurian terrestrial cover rocks were buried under these conditions. Acadian-Neoacadian deformation (393-340 Ma) becomes younger towards the northwest and progressively localized in transcurrent fault zones. This final foreland deformation at shallow crustal level established the Acadian/Neoacadian orogenic front in central Newfoundland slightly northwest of the RIL.

Both enstatite spectra are very bright in the VIS and NIR and show almost neutral to slightly bluish spectral slopes with a steep absorption in the UV. Very low iron in the enstatites (below similar to 0.04 wt% FeO) already results in weak albeit noticeable absorptions in the VNIR between 0.4 and 0.9 mu m. Orthoenstatite and clinoenstatite are not distinguishable based only on their spectra in the VIS and NIR. At the Reststrahlen bands in the MIR a systematic difference in the number and exact position of local minima at similar to 10 mu m between clinoenstatite and orthoenstatite is evident. This can be used to discern between the polymorphs in this wavelength range. Additionally, we can distinguish between Fe-free low- and high-Ca pyroxenes in the MIR.

Plain Language Summary The 2015 Gorkha earthquake in Nepal caused severe losses in the hydropower sector. The country temporarily lost similar to 20% of its hydropower capacity, and >30 hydropower projects were damaged. The projects hit hardest were those that were affected by earthquake-triggered landslides. We show that these projects are located along very steep rivers with towering sidewalls that are prone to become unstable during strong seismic ground shaking. A statistical classification based on a topographic metric that expresses river steepness and earthquake ground acceleration is able to approximately predict hydropower damage during future earthquakes, based on successful testing of past cases. Thus, our model enables us to estimate earthquake damages to hydropower projects in other parts of the Himalayas. We find that >10% of the Himalayan drainage network may be unsuitable for hydropower infrastructure given high probabilities of high earthquake damages.

More than 41% of the Earth’s land area is covered by permanent or seasonally arid dryland ecosystems. Global development and human activity have led to an increase in aridity, resulting in ecosystem degradation and desertification around the world. The objective of the present work was to investigate and compare the microbial community structure and geochemical characteristics of two geographically distinct saline pan sediments in the Kalahari Desert of southern Africa. Our data suggest that these microbial communities have been shaped by geochemical drivers, including water content, salinity, and the supply of organic matter. Using Illumina 16S rRNA gene sequencing, this study provides new insights into the diversity of bacteria and archaea in semi-arid, saline, and low-carbon environments. Many of the observed taxa are halophilic and adapted to water-limiting conditions. The analysis reveals a high relative abundance of halophilic archaea (primarily Halobacteria), and the bacterial diversity is marked by an abundance of Gemmatimonadetes and spore-forming Firmicutes. In the deeper, anoxic layers, candidate division MSBL1, and acetogenic bacteria (Acetothermia) are abundant. Together, the taxonomic information and geochemical data suggest that acetogenesis could be a prevalent form of metabolism in the deep layers of a saline pan.

Tomato (Solanum lycopersicum L.) being an important vegetable is cultivated and used throughout the world. It not only contributes in fulfilling the basic nutritional requirements of the human body but also has many health benefits due to its rich biochemical composition. However, its production at large scale is hampered by many limiting factors such as biotic and abiotic stresses. Among the different abiotic stresses, drought poses drastic impact on tomato yield. Drought stress is genetically regulated by many transcription factors that not only regulate the stress responsive mechanism but also facilitate the growth and development of tomato plants. NAC is an important stress related transcription factor genes family, and the ATAF1 gene, a member of this family, is involved in ABA signaling and stress response. In this study, tomato variety Rio Drande was transformed with drought resistant ATAF1 gene via Agrobacterium mediated gene transformation method. The ATAF1 gene was first cloned in the pK7WFG2 vector having kanamycin selectable marker and then it was introduced in the Agrobacterium tumefaciens strain GV3101 through heat shock method. The tomato cotyledon and hypocotyl ex-plants of variety "Rio Ggrande" were cultured on callus induction medium (MS + 2.5 mg/L IAA + 2 mg/L BAP). The calli were then infected with Agrobacterium tumefaciens strain GV3101 containing ATAF1 gene and selection was carried out on the kanamycin selectable medium (MS + 100 mg/L Kan), and were regenerated on MS medium with 1 mg/L IAA + 1 mg/L BAP. Out of 216 putative transformed calli, 13 calli were able to regenerate on the selection medium. Of the 13 calli, three transgenic tomato plantlets were recovered, and these were confirmed through PCR analysis for the presence of 432 bp fragment of ATAF1 gene. The transformation protocol reported here can be used to generate drought resistant tomato plants in future.

Risk-based insurance is a commonly proposed and discussed flood risk adaptation mechanism in policy debates across the world such as in the United Kingdom and the United States of America. However, both risk-based premiums and growing risk pose increasing difficulties for insurance to remain affordable. An empirical concept of affordability is required as the affordability of adaption strategies is an important concern for policymakers, yet such a concept is not often examined. Therefore, a robust metric with a commonly acceptable affordability threshold is required. A robust metric allows for a previously normative concept to be quantified in monetary terms, and in this way, the metric is rendered more suitable for integration into public policy debates. This paper investigates the degree to which risk-based flood insurance premiums are unaffordable in Europe. In addition, this paper compares the outcomes generated by three different definitions of unaffordability in order to investigate the most robust definition. In doing so, the residual income definition was found to be the least sensitive to changes in the threshold. While this paper focuses on Europe, the selected definition can be employed elsewhere in the world and across adaption measures in order to develop a common metric for indicating the potential unaffordability problem.

One paragraph of the manuscript of the paper has been inadvertently omitted in the very final stage of its compilation due to a technical mistake. Since this paragraph discusses the declustering of the used earthquake catalogue and is therefore necessary for the understanding of the seismicity data preprocessing, the authors decided to provide this paragraph in form of a correction. The respective paragraph belongs to chapter 2 of the paper, where it was placed originally, and should be inserted into the published paper before the second to the last paragraph. The omitted text reads as follows:

We present the first high-resolution trace element (Mg/Ca, Sr/Ca, Ba/Ca) record from a stalagmite in southwestern Romania covering the last 3.6 ka, which provides the potential for quantitative climate reconstruction. Precise age control is based on three independent dating methods, in particular for the last 250 yr, where chemical lamina counting is combined with the identification of the 20th century radiocarbon bomb peak and Th-230/U dating. Long-term cave monitoring and model simulations of drip water and speleothem elemental variability indicate that precipitation-related processes are the main drivers of speleothem Mg/Ca ratios. Calibration against instrumental climate data shows a significant anti-correlation of speleothem Mg/Ca ratios with autumn/winter (October to March) precipitation (r = -0.61, p < 0.01), which is statistically robust when considering age uncertainties and auto-correlation. This relationship is used to develop a quantitative reconstruction of autumn/winter precipitation. During the late Holocene, our data suggest a heterogeneous pattern of past regional winter hydroclimate in the Carpathian/Balkan realm, along with intermittent weakening of the dominant influence of North Atlantic forcing. In agreement with other regional paleo-hydrological reconstructions, the observed variability reveals periodically occurring strong NW-SE hydro-climate gradients. We hypothesize, that this pattern is caused by shifts of the eastern edge of the area of influence of the NAO across central eastern Europe due to the confluence of North Atlantic forcing, and other climatic features such as the East Atlantic/Western Russia (EAWR) pattern. (C) 2018 Elsevier B.V. All rights reserved.

Process life cycle assessment (PLCA) is widely used to quantify environmental flows associated with the manufacturing of products and other processes. As PLCA always depends on defining a system boundary, its application involves truncation errors. Different methods of estimating truncation errors are proposed in the literature; most of these are based on artificially constructed system complete counterfactuals. In this article, we review the literature on truncation errors and their estimates and systematically explore factors that influence truncation error estimates. We classify estimation approaches, together with underlying factors influencing estimation results according to where in the estimation procedure they occur. By contrasting different PLCA truncation/error modeling frameworks using the same underlying input-output (I-O) data set and varying cut-off criteria, we show that modeling choices can significantly influence estimates for PLCA truncation errors. In addition, we find that differences in I-O and process inventory databases, such as missing service sector activities, can significantly affect estimates of PLCA truncation errors. Our results expose the challenges related to explicit statements on the magnitude of PLCA truncation errors. They also indicate that increasing the strictness of cut-off criteria in PLCA has only limited influence on the resulting truncation errors. We conclude that applying an additional I-O life cycle assessment or a path exchange hybrid life cycle assessment to identify where significant contributions are located in upstream layers could significantly reduce PLCA truncation errors.

Increasing air temperatures are changing the arctic tundra biome. Permafrost is thawing, snow duration is decreasing, shrub vegetation is proliferating, and boreal wildlife is encroaching. Here we present evidence of the recent range expansion of North American beaver (Castor canadensis) into the Arctic, and consider how this ecosystem engineer might reshape the landscape, biodiversity, and ecosystem processes. We developed a remote sensing approach that maps formation and disappearance of ponds associated with beaver activity. Since 1999, 56 new beaver pond complexes were identified, indicating that beavers are colonizing a predominantly tundra region (18,293km(2)) of northwest Alaska. It is unclear how improved tundra stream habitat, population rebound following overtrapping for furs, or other factors are contributing to beaver range expansion. We discuss rates and likely routes of tundra beaver colonization, as well as effects on permafrost, stream ice regimes, and freshwater and riparian habitat. Beaver ponds and associated hydrologic changes are thawing permafrost. Pond formation increases winter water temperatures in the pond and downstream, likely creating new and more varied aquatic habitat, but specific biological implications are unknown. Beavers create dynamic wetlands and are agents of disturbance that may enhance ecosystem responses to warming in the Arctic.

Rotational motions play a key role in measuring seismic wavefield properties. Using newly developed portable rotational instruments, it is now possible to directly measure rotational motions in a broad frequency range. Here, we investigated the instrumental self-noise and data quality in a huddle test in Fürstenfeldbruck, Germany, in August 2019. We compare the data from six rotational and three translational sensors. We studied the recorded signals using correlation, coherence analysis, and probabilistic power spectral densities. We sorted the coherent noise into five groups with respect to the similarities in frequency content and shape of the signals. These coherent noises were most likely caused by electrical devices, the dehumidifier system in the building, humans, and natural sources such as wind. We calculated self-noise levels through probabilistic power spectral densities and by applying the Sleeman method, a three-sensor method. Our results from both methods indicate that self-noise levels are stable between 0.5 and 40 Hz. Furthermore, we recorded the 29 August 2019 ML 3.4 Dettingen earthquake. The calculated source directions are found to be realistic for all sensors in comparison to the real back azimuth. We conclude that the five tested blueSeis-3A rotational sensors, when compared with respect to coherent noise, self-noise, and source direction, provide reliable and consistent results. Hence, field experiments with single rotational sensors can be undertaken.

Forming as a result of the collision between the Adriatic and European plates, the Alpine orogen exhibits significant lithospheric heterogeneity due to the long history of interplay between these plates, other continental and oceanic blocks in the region, and inherited features from preceeding orogenies. This implies that the thermal and rheological configuration of the lithosphere also varies significantly throughout the region. Lithology and temperature/pressure conditions exert a first order control on rock strength, principally via thermally activated creep deformation and on the distribution at depth of the brittle-ductile transition zone, which can be regarded as the lower bound to the seismogenic zone. Therefore, they influence the spatial distribution of seismicity within a lithospheric plate. In light of this, accurately constrained geophysical models of the heterogeneous Alpine lithospheric configuration, are crucial in describing regional deformation patterns. However, despite the amount of research focussing on the area, different hypotheses still exist regarding the present-day lithospheric state and how it might relate to the present-day seismicity distribution. This dissertaion seeks to constrain the Alpine lithospheric configuration through a fully 3D integrated modelling workflow, that utilises multiple geophysical techniques and integrates from all available data sources. The aim is therefore to shed light on how lithospheric heterogeneity may play a role in influencing the heterogeneous patterns of seismicity distribution observed within the region. This was accomplished through the generation of: (i) 3D seismically constrained, structural and density models of the lithosphere, that were adjusted to match the observed gravity field; (ii) 3D models of the lithospheric steady state thermal field, that were adjusted to match observed wellbore temperatures; and (iii) 3D rheological models of long term lithospheric strength, with the results of each step used as input for the following steps. Results indicate that the highest strength within the crust (~ 1 GPa) and upper mantle (> 2 GPa), are shown to occur at temperatures characteristic for specific phase transitions (more felsic crust: 200 – 400 °C; more mafic crust and upper lithospheric mantle: ~600 °C) with almost all seismicity occurring in these regions. However, inherited lithospheric heterogeneity was found to significantly influence this, with seismicity in the thinner and more mafic Adriatic crust (~22.5 km, 2800 kg m−3, 1.30E-06 W m-3) occuring to higher temperatures (~600 °C) than in the thicker and more felsic European crust (~27.5 km, 2750 kg m−3, 1.3–2.6E-06 W m-3, ~450 °C). Correlation between seismicity in the orogen forelands and lithospheric strength, also show different trends, reflecting their different tectonic settings. As such, events in the plate boundary setting of the southern foreland correlate with the integrated lithospheric strength, occurring mainly in the weaker lithosphere surrounding the strong Adriatic indenter. Events in the intraplate setting of the northern foreland, instead correlate with crustal strength, mainly occurring in the weaker and warmer crust beneath the Upper Rhine Graben. Therefore, not only do the findings presented in this work represent a state of the art understanding of the lithospheric configuration beneath the Alps and their forelands, but also a significant improvement on the features known to significantly influence the occurrence of seismicity within the region. This highlights the importance of considering lithospheric state in regards to explaining observed patterns of deformation.

The variation of Rayleigh ellipticity versus frequency is gaining popularity in site characterization. It becomes a necessary observable to complement dispersion curves when inverting shear wave velocity profiles. Various methods have been proposed so far to extract polarization from ambient vibrations recorded on a single three-component station or with an array of three-component sensors. If only absolute values were recovered 10 yr ago, new array-based techniques were recently proposed with enhanced efficiencies providing also the ellipticity sign. With array processing, higher-order modes are often detected even in the ellipticity domain. We suggest to explore the properties of a high-resolution beamforming where radial and vertical components are explicitly included. If N is the number of three-component sensors, 2N x 2N cross-spectral density matrices are calculated for all presumed directions of propagation. They are built with N radial and N vertical channels. As a first approach, steering vectors are designed to fit with Rayleigh wave properties: the phase shift between radial and vertical components is either -Pi/2 or Pi/2. We show that neglecting the ellipticity tilt due to attenuation has only minor effects on the results. Additionally, we prove analytically that it is possible to retrieve the ellipticity value from the usual maximization of the high-resolution beam power. The method is tested on synthetic data sets and on experimental data. Both are reference sites already analysed by several authors. A detailed comparison with previous results on these cases is provided.

Bedrock landsliding in mountain belts can elevate overall chemical weathering rates through rapid dissolution of exhumed reactive mineral phases in transiently stored deposits. This link between a key process of erosion and the resultant weathering affects the sequestering of carbon dioxide through weathering of silicate minerals and broader links between erosion in active orogens and climate change. Here we address the effect on the carbon cycle of weathering induced by bedrock landsliding in Taiwan and the Western Southern Alps of New Zealand. Using solute chemistry data from samples of seepage from landslide deposits and river discharge from catchments with variable proportions of landsliding, we model the proportion of silicate and carbonate weathering and the balance of sulfuric and carbonic acids that act as weathering agents. We correct for secondary precipitation, geothermal, and cyclic input, to find a closer approximation of the weathering explicitly occurring within landslide deposits. We find highly variable proportions of sulfuric and carbonic acids driving weathering in landslides and stable hillslopes. Despite this variability, the predominance of rapid carbonate weathering within landslides and catchments where mass wasting is prevalent results at best in limited sequestration of carbon dioxide by this process of rapid erosion. In many cases where sulfuric acid is a key weathering agent, a net release of CO2 to the atmosphere occurs. This suggests that a causal link between erosion in mountain belts and climate change through the sequestration of CO2, if it exists, must operate through a process other than chemical weathering driven by landsliding. Plain Language Summary There is a long-standing debate surrounding the link between erosion and climate. It is often suggested that as temperatures increase, rainier and stormier weather could increase erosion of rock; as that rock is exposed, silicate minerals within could break down, which, on long time scales, can remove CO2 from the atmosphere, lowering global temperatures and acting as a negative feedback. Recent studies have shown that landslide deposits are key locations for the link between chemical weathering and physical erosion in some mountain belts. To test how landslides affect the erosion-climate link, we used samples of water seeping through landslides in Taiwan and New Zealand to calculate the amount of carbon dioxide that is either absorbed or released through this chemical reaction. We find that the large amount of freshly exposed rock in Taiwanese landslide deposits contains significant carbonate rock and sulfide minerals; the net result of the weathering of these minerals is a release of carbon dioxide, which inverts the traditional perspective on the role erosion plays in controlling carbon dioxide release. In some mountain belts, it seems that increased erosion and resulting landsliding may act to increase carbon dioxide in the air, opening further questions into the nature of erosional-climatic links.

Models for the predictions of monetary losses from floods mainly blend data deemed to represent a single flood type and region. Moreover, these approaches largely ignore indicators of preparedness and how predictors may vary between regions and events, challenging the transferability of flood loss models. We use a flood loss database of 1812 German flood-affected households to explore how Bayesian multilevel models can estimate normalised flood damage stratified by event, region, or flood process type. Multilevel models acknowledge natural groups in the data and allow each group to learn from others. We obtain posterior estimates that differ between flood types, with credibly varying influences of water depth, contamination, duration, implementation of property-level precautionary measures, insurance, and previous flood experience; these influences overlap across most events or regions, however. We infer that the underlying damaging processes of distinct flood types deserve further attention. Each reported flood loss and affected region involved mixed flood types, likely explaining the uncertainty in the coefficients. Our results emphasise the need to consider flood types as an important step towards applying flood loss models elsewhere. We argue that failing to do so may unduly generalise the model and systematically bias loss estimations from empirical data.

Silicate melts are major components of the Earth’s interior and as such they make an essential contribution in igneous processes, in the dynamics of the solid Earth and the chemical development of the entire Earth. Macroscopic physical and chemical properties such as density, compressibility, viscosity, degree of polymerization etc. are determined by the atomic structure of the melt. Depending on the pressure, but also on the temperature and the chemical composition, silicate melts show different structural properties. These properties are best described by the local coordination environment, i.e. symmetry and number of neighbors (coordination number) of an atom, as well as the distance between the central atom and its neighbors (inter-atomic distance). With increasing pressure and temperature, i.e. with increasing depth in the Earth, the density of the melt increases, which can lead to changes in coordination number and distances. If the coordination number remains the same, the distance usually decreases. If the coordination number increases, the distance can increase. These general trends can, however, vary greatly, which can be attributed in particular to the chemical composition. Due to the fact that natural melts of the deep earth are not accessible to direct investigations, in order to understand their properties under the relevant conditions, extensive experimental and theoretical investigations have been carried out so far. This has often been studied using the example of amorphous samples of the end-members SiO2 and GeO2 , with the latter serving as a structural and chemical analog model to SiO2. Commonly, the experiments were carried out at high pressure and at room temperature. Natural melts are chemically much more complex than the simple end-member SiO2 and GeO2, so that observations made on them may lead to incorrect compression models. Furthermore, the investigations on glasses at room temperature can show potentially strong deviations from the properties of melts under natural thermodynamic conditions. The aim of this thesis was to explain the influence of the composition and the temperature on the structural properties of the melts at high pressures. To understand this, we studied complex alumino-germanate and alumino-silicate glasses. More precisely, we studied synthetic glasses that have a composition like the mineral albite and like a mixture of albite-diopside at the eutectic point. The albite glass is structurally similar to a simplified granitic melt, while the albite-diopside glass simulates a simplified basaltic melt. To study the local coordination environment of the elements, we used X-ray absorption spectroscopy in combination with a diamond anvil cell. Because the diamonds have a high absorbance for X-rays with energies below 10 keV, the direct investigation of the geologically relevant elements such as Si, Al, Ca, Mg etc. with this spectroscopic probe technique in combination with a diamond anvil cell is not possible. Therefore the glasses were doped with Ge and Sr. These elements serve partially or fully as substitutes for important major elements. In this sense, Ge serves as an a substitute for Si and other network formers, while Sr replaces network modifiers such as Ca, Na, Mg etc., as well as other cations with a large ionic radius. In the first step we studied the Ge K-edge in Ge-Albit-glass, NaAlGe3O8, at room temperature up to 131 GPa. This glass has a higher chemical complexity than SiO2 and GeO2, but it is still fully polymerized. The differences in the compression mechanism between this glass and the simple oxides can clearly be attributed to higher chemical complexity. The albite and albite-diopside compositions partially doped with Ge and Sr were probed at room temperature for Ge up to 164 GPa and for Sr up to 42 GPa. While the albite glass is nominally fully polymerized like NaAlGe3O8, the albite-diopside glass is partially depolymerized. The results show that structural changes take place in all three glasses in the first 25 to a maximum of 30 GPa, with both Ge and Sr reaching the maximum coordination number 6 and ∼9, respectively. At higher pressures, only isostructural shrinkage of the coordination polyhedra takes place in the glasses. The most important finding of the high pressure studies on the alumino-silicate and alumino-germanate glasses is that in these complex glasses the polyhedra show a much higher compressibility than what can be observed in the end-members. This is shown in particular by the strong shortening of the Ge-O distances in the amorphous NaAlGe3O8 and albite-diopside glass at pressures above 30 GPa. In addition to the effects of the composition on the compaction process, we investigated the influence of temperature on the structural changes. To do this, we probed the albite-diopside glass, as it is chemically most similar to the melts in the lower mantle. We studied the Ge K edge of the sample with a resistively heated and a laser-heated diamond anvil cell, for a pressure range of up to 48 GPa and a temperature range of up to 5000 K. High temperatures at which the sample is liquid and that are relevant for the Earth mantle, have a significant impact on the structural transformation, with a shift of approx. 30% to significantly lower pressures, compared to the glasses at room temperature and below 1000 K. The results of this thesis represent an important contribution to the understanding of the properties of melts at conditions of the lower mantle. In the context of the discussion about the existence and origin of ultra-dense silicate melts at the core-mantle boundary, these investigations show that the higher density compared to the surrounding material cannot be explained by only structural features, but by a distinct chemical composition. The results also suggest that only very low solubilities of noble gases are to be expected for melts in the lower mantle, so that the structural properties clearly influence the overall budget and transport of noble gases in the Earth’s mantle.

Thermal conductivity (lambda) is an essential physical property of minerals and rocks and fundamental in constraining the thermal field of the lithosphere. In case that adequate samples to measure lambda are not available, it could be indirectly inferred from calculation. One of the most widely applied indirect methods for rocks involve modal mineralogy and porosity as parameters that are incorporated into mathematical mean or mixing models. Robust inferences from these approaches for crystalline rocks were impeded by a small number of studied samples or restriction to certain rock types. We employ this method and examine its applicability to low-porosity plutonic rocks by calculating bulk thermal conductivity lambda(b) for 45 samples covering the entire range from gabbro/diorite to granite. We show that the use of the harmonic-mean model for both rock matrix and porosity provided a good match between lambda(b.meas) and lambda(b.calc) of <10% deviation (2 sigma), with relative and absolute errors amounting to 1.49.7% and 4.44.9%, respectively. The results of our study constitute a big step forward to a robust conclusion on the overall applicability of the harmonic-mean model for inferring lambda(b) of isotropic, low-porosity, mafic to silicic plutonic and metamorphic rocks with an acceptable magnitude of error. Drill cuttings and enclaves form particularly interesting objects for application of this method, as they are poorly suited for direct measurement. Well-derived lambda values for those rocks would permit to calculate heat flow and to model more profoundly the thermal state of the deeper lithosphere.

Knowledge of the present-day crustal in-situ stress field is a key for the understanding of geodynamic processes such as global plate tectonics and earthquakes. It is also essential for the management of geo-reservoirs and underground storage sites for energy and waste. Since 1986, the World Stress Map (WSM) project has systematically compiled the orientation of maximum horizontal stress (S-Hmax). For the 30th anniversary of the project, the WSM database has been updated significantly with 42,870 data records which is double the amount of data in comparison to the database release in 2008. The update focuses on areas with previously sparse data coverage to resolve the stress pattern on different spatial scales. In this paper, we present details of the new WSM database release 2016 and an analysis of global and regional stress pattern. With the higher data density, we can now resolve stress pattern heterogeneities from plate-wide to local scales. In particular, we show two examples of 40 degrees-60 degrees S-Hmax rotations within 70 km. These rotations can be used as proxies to better understand the relative importance of plate boundary forces that control the long wave-length pattern in comparison to regional and local controls of the crustal stress state. In the new WSM project phase IV that started in 2017, we will continue to further refine the information on the S-Hmax orientation and the stress regime. However, we will also focus on the compilation of stress magnitude data as this information is essential for the calibration of geomechanical-numerical models. This enables us to derive a 3-D continuous description of the stress tensor from point-wise and incomplete stress tensor information provided with the WSM database. Such forward models are required for safety aspects of anthropogenic activities in the underground and for a better understanding of tectonic processes such as the earthquake cycle.

Himalayan rivers are frequently hit by catastrophic floods that are caused by the failure of glacial lake and landslide dams; however, the dynamics and long-term impacts of such floods remain poorly understood. We present a comprehensive set of observations that capture the July 2016 glacial lake outburst flood (GLOF) in the Bhotekoshi/Sunkoshi River of Nepal. Seismic records of the flood provide new insights into GLOF mechanics and their ability to mobilize large boulders that otherwise prevent channel erosion. Because of this boulder mobilization, GLOF impacts far exceed those of the annual summer monsoon, and GLOFs may dominate fluvial erosion and channel-hillslope coupling many tens of kilometers downstream of glaciated areas. Long-term valley evolution in these regions may therefore be driven by GLOF frequency and magnitude, rather than by precipitation.

Glacial deposits on the high-altitude, arid southern Central Andean Plateau (CAP), the Puna in northwestern Argentina, document past changes in climate, but the associated geomorphic features have rarely been directly dated. This study provides direct age control of glacial moraine deposits from the central Puna (24 degrees S) at elevations of 3900-5000 m through surface exposure dating with cosmogenic nuclides. Our results show that the most extensive glaciations occurred before 95 ka and an additional major advance occurred between 46 and 39 ka. The latter period is synchronous with the highest lake levels in the nearby Pozuelos basin and the Minchin (Inca Huasi) wet phase on the Altiplano in the northern CAP. None of the dated moraines produced boulder ages corresponding to the Tauca wet phase (24-15 ka). Additionally, the volcanic lithologies of the deposits allow us to establish production ratios at low latitude and high elevation for five different nuclide and mineral systems: Be-10, Ne-21, and Al-26 from quartz (11 or 12 samples) and He-3 and Ne-21 from pyroxene (10 samples). We present production ratios for all combinations of the measured nuclides and cross-calibrated production rates for 21Ne in pyroxene and quartz for the high, (sub-)tropical Andes. The production rates are based on our Be-10-normalized production ratios and a weighted mean of reference 10Be production rates calibrated in the high, tropical Andes (4.02 +/- 0.12 at g(-1) yr(-1)). These are, Ne-21(qtz): 18.1 +/- 1.2 at g(-1) yr(-1) and Ne-21(px): 36.6 +/- 1.8 at g(-1) yr(-1) (En(88-94)) scaled to sea level and high latitude using the Lal/Stone scheme, with 1 sigma uncertainties. As He-3 and Al-26 have been directly calibrated in the tropical Andes, we recommend using those rates. Finally, we compare exposure ages calculated using all measured cosmogenic nuclides from each sample, including 11 feldspar samples measured for Cl-36, and a suite of previously published production rates. (C) 2018 Published by Elsevier B.V.

Understanding the dynamics between the East Asian summer (EASM) and winter monsoon (EAWM) is needed to predict their variability under future global warming scenarios. Here, we investigate the relationship between EASM and EAWM as well as the mechanisms driving their variability during the last 10,000 years by stacking marine and terrestrial (non-speleothem) proxy records from the East Asian realm. This provides a regional and proxy independent signal for both monsoonal systems. The respective signal was subsequently analysed using a linear regression model. We find that the phase relationship between EASM and EAWM is not time-constant and significantly depends on orbital configuration changes. In addition, changes in the Atlantic Meridional Overturning circulation, Arctic sea-ice coverage, El Niño-Southern Oscillation and Sun Spot numbers contributed to millennial scale changes in the EASM and EAWM during the Holocene. We also argue that the bulk signal of monsoonal activity captured by the stacked non-speleothem proxy records supports the previously argued bias of speleothem climatic archives to moisture source changes and/or seasonality.

Thermokarst lakes are prevalent in Arctic coastal lowland regions and sublake permafrost degradation and talik development contributes to greenhouse gas emissions by tapping the large permafrost carbon pool. Whereas lateral thermokarst lake expansion is readily apparent through remote sensing and shoreline measurements, sublake thawed sediment conditions and talik growth are difficult to measure. Here we combine transient electromagnetic surveys with thermal modeling, backed up by measured permafrost properties and radiocarbon ages, to reveal closed-talik geometry associated with a thermokarst lake in continuous permafrost. To improve access to talik geometry data, we conducted surveys along three transient electromagnetic transects perpendicular to lakeshores with different decadal-scale expansion rates of 0.16, 0.38, and 0.58m/year. We modeled thermal development of the talik using boundary conditions based on field data from the lake, surrounding permafrost and a borehole, independent of the transient electromagnetics. A talik depth of 91m was determined from analysis of the transient electromagnetic surveys. Using a lake initiation age of 1400years before present and available subsurface properties the results from thermal modeling of the lake center arrived at a best estimate talk depth of 80m, which is on the same order of magnitude as the results from the transient electromagnetic survey. Our approach has provided a noninvasive estimate of talik geometry suitable for comparable settings throughout circum-Arctic coastal lowland regions.

Reading the sediment record in terms of past climates is challenging since linking climate change to the associated responses of sedimentary systems is not always straightforward. Here we analyze the erosional response of landscapes on the Tibetan Plateau to interglacial climate forcing. Using the theory of dynamical systems on Holocene time series of geochemical proxies, we derive a sedimentary response model that accurately simulates observed proxy variation in three lake records. The model suggests that millennial variations in sediment composition reflect a self-organization of landscapes in response to abrupt climate change between 11.6 and 11.9 ka BP. The self-organization is characterized by oscillations in sediment supply emerging from a feedback between physical and chemical erosion processes, with estimated response times between 3,000 to 18,000 years depending on catchment topography. The implications of our findings emphasize the need for landscape response models to decipher the paleoclimatic code in continental sediment records. Plain Language Summary Lake sediments are an important source of information on past climates. Reading the information is not always straightforward. Complex interactions in landscapes can affect the transmission of climatic signals to the sediment record. However, the exact nature of such complex interactions remains unknown. This study compares sediment deposits of three lakes on the Tibetan Plateau. The deposits are continuous records of landscape responses to climate change during the last 12,000 years. We identified a mathematical model that accurately simulates changes in sediment composition at all sites. The model simulations suggest that an abrupt warming at the end of the last glacial period destabilized the landscapes. This caused fluctuations in the transport of sediments, which persisted for several thousand years. Our findings present evidence for a long-lasting impact of abrupt climate change on fundamental Earth surface processes.

High-latitude treeless ecosystems represent spatially highly heterogeneous landscapes with small net carbon fluxes and a short growing season. Reliable observations and process understanding are critical for projections of the carbon balance of the climate-sensitive tundra. Space-borne remote sensing is the only tool to obtain spatially continuous and temporally resolved information on vegetation greenness and activity in remote circumpolar areas. However, confounding effects from persistent clouds, low sun elevation angles, numerous lakes, widespread surface inundation, and the sparseness of the vegetation render it highly challenging. Here, we conduct an extensive analysis of the timing of peak vegetation productivity as shown by satellite observations of complementary indicators of plant greenness and photosynthesis. We choose to focus on productivity during the peak of the growing season, as it importantly affects the total annual carbon uptake. The suite of indicators are as follows: (1) MODIS-based vegetation indices (VIs) as proxies for the fraction of incident photosynthetically active radiation (PAR) that is absorbed (fPAR), (2) VIs combined with estimates of PAR as a proxy of the total absorbed radiation (APAR), (3) sun-induced chlorophyll fluorescence (SIF) serving as a proxy for photosynthesis, (4) vegetation optical depth (VOD), indicative of total water content and (5) empirically upscaled modelled gross primary productivity (GPP). Averaged over the pan-Arctic we find a clear order of the annual peak as APAR <= GPP < SIF < VIs/VOD. SIF as an indicator of photosynthesis is maximised around the time of highest annual temperatures. The modelled GPP peaks at a similar time to APAR. The time lag of the annual peak between APAR and instantaneous SIF fluxes indicates that the SIF data do contain information on light-use efficiency of tundra vegetation, but further detailed studies are necessary to verify this. Delayed peak greenness compared to peak photosynthesis is consistently found across years and land-cover classes. A particularly late peak of the normalised difference vegetation index (NDVI) in regions with very small seasonality in greenness and a high amount of lakes probably originates from artefacts. Given the very short growing season in circumpolar areas, the average time difference in maximum annual photosynthetic activity and greenness or growth of 3 to 25 days (depending on the data sets chosen) is important and needs to be considered when using satellite observations as drivers in vegetation models.

Geochemical processes such as mineral dissolution and precipitation alter the microstructure of rocks, and thereby affect their hydraulic and mechanical behaviour. Quantifying these property changes and considering them in reservoir simulations is essential for a sustainable utilisation of the geological subsurface. Due to the lack of alternatives, analytical methods and empirical relations are currently applied to estimate evolving hydraulic and mechanical rock properties associated with chemical reactions. However, the predictive capabilities of analytical approaches remain limited, since they assume idealised microstructures, and thus are not able to reflect property evolution for dynamic processes. Hence, aim of the present thesis is to improve the prediction of permeability and stiffness changes resulting from pore space alterations of reservoir sandstones. A detailed representation of rock microstructure, including the morphology and connectivity of pores, is essential to accurately determine physical rock properties. For that purpose, three-dimensional pore-scale models of typical reservoir sandstones, obtained from highly resolved micro-computed tomography (micro-CT), are used to numerically calculate permeability and stiffness. In order to adequately depict characteristic distributions of secondary minerals, the virtual samples are systematically altered and resulting trends among the geometric, hydraulic, and mechanical rock properties are quantified. It is demonstrated that the geochemical reaction regime controls the location of mineral precipitation within the pore space, and thereby crucially affects the permeability evolution. This emphasises the requirement of determining distinctive porosity-permeability relationships by means of digital pore-scale models. By contrast, a substantial impact of spatial alterations patterns on the stiffness evolution of reservoir sandstones are only observed in case of certain microstructures, such as highly porous granular rocks or sandstones comprising framework-supporting cementations. In order to construct synthetic granular samples a process-based approach is proposed including grain deposition and diagenetic cementation. It is demonstrated that the generated samples reliably represent the microstructural complexity of natural sandstones. Thereby, general limitations of imaging techniques can be overcome and various realisations of granular rocks can be flexibly produced. These can be further altered by virtual experiments, offering a fast and cost-effective way to examine the impact of precipitation, dissolution or fracturing on various petrophysical correlations. The presented research work provides methodological principles to quantify trends in permeability and stiffness resulting from geochemical processes. The calculated physical property relations are directly linked to pore-scale alterations, and thus have a higher accuracy than commonly applied analytical approaches. This will considerably improve the predictive capabilities of reservoir models, and is further relevant to assess and reduce potential risks, such as productivity or injectivity losses as well as reservoir compaction or fault reactivation. Hence, the proposed method is of paramount importance for a wide range of natural and engineered subsurface applications, including geothermal energy systems, hydrocarbon reservoirs, CO2 and energy storage as well as hydrothermal deposit exploration.

One common approach to cope with floods is the implementation of structural flood protection measures, such as levees or flood-control reservoirs, which substantially reduce the probability of flooding at the time of implementation. Numerous scholars have problematized this approach. They have shown that increasing the levels of flood protection can attract more settlements and high-value assets in the areas protected by the new measures. Other studies have explored how structural measures can generate a sense of complacency, which can act to reduce preparedness. These paradoxical risk changes have been described as "levee effect", "safe development paradox" or "safety dilemma". In this commentary, we briefly review this phenomenon by critically analysing the intended benefits and unintended effects of structural flood protection, and then we propose an interdisciplinary research agenda to uncover these paradoxical dynamics of risk.

Rheology describes the flow of matter under the influence of stress, and - related to solids- it investigates how solids subjected to stresses deform. As the deformation of the Earth’s outer layers, the lithosphere and the crust, is a major focus of rheological studies, rheology in the geosciences describes how strain evolves in rocks of variable composition and temperature under tectonic stresses. It is here where deformation processes shape the form of ocean basins and mountain belts that ultimately result from the complex interplay between lithospheric plate motion and the susceptibility of rocks to the influence of plate-tectonic forces. A rigorous study of the strength of the lithosphere and deformation phenomena thus requires in-depth studies of the rheological characteristics of the involved materials and the temporal framework of deformation processes. This dissertation aims at analyzing the influence of the physical configuration of the lithosphere on the present-day thermal field and the overall rheological characteristics of the lithosphere to better understand variable expressions in the formation of passive continental margins and the behavior of strike-slip fault zones. The main methodological approach chosen is to estimate the present-day thermal field and the strength of the lithosphere by 3-D numerical modeling. The distribution of rock properties is provided by 3-D structural models, which are used as the basis for the thermal and rheological modeling. The structural models are based on geophysical and geological data integration, additionally constrained by 3-D density modeling. More specifically, to decipher the thermal and rheological characteristics of the lithosphere in both oceanic and continental domains, sedimentary basins in the Sea of Marmara (continental transform setting), the SW African passive margin (old oceanic crust), and the Norwegian passive margin (young oceanic crust) were selected for this study. The Sea of Marmara, in northwestern Turkey, is located where the dextral North Anatolian Fault zone (NAFZ) accommodates the westward escape of the Anatolian Plate toward the Aegean. Geophysical observations indicate that the crust is heterogeneous beneath the Marmara basin, but a detailed characterization of the lateral crustal heterogeneities is presented for the first time in this study. Here, I use different gravity datasets and the general non-uniqueness in potential field modeling, to propose three possible end-member scenarios of crustal configuration. The models suggest that pronounced gravitational anomalies in the basin originate from significant density heterogeneities within the crust. The rheological modeling reveals that associated variations in lithospheric strength control the mechanical segmentation of the NAFZ. Importantly, a strong crust that is mechanically coupled to the upper mantle spatially correlates with aseismic patches where the fault bends and changes its strike in response to the presence of high-density lower crustal bodies. Between the bends, mechanically weaker crustal domains that are decoupled from the mantle are characterized by creep. For the passive margins of SW Africa and Norway, two previously published 3-D conductive and lithospheric-scale thermal models were analyzed. These 3-D models differentiate various sedimentary, crustal, and mantle units and integrate different geophysical data, such as seismic observations and the gravity field. Here, the rheological modeling suggests that the present-day lithospheric strength across the oceanic domain is ultimately affected by the age and past thermal and tectonic processes as well as the depth of the thermal lithosphere-asthenosphere boundary, while the configuration of the crystalline crust dominantly controls the rheological behavior of the lithosphere beneath the continental domains of both passive margins. The thermal and rheological models show that the variations of lithospheric strength are fundamentally influenced by the temperature distribution within the lithosphere. Moreover, as the composition of the lithosphere significantly influences the present-day thermal field, it therefore also affects the rheological characteristics of the lithosphere. Overall my studies add to our understanding of regional tectonic deformation processes and the long-term behavior of sedimentary basins; they confirm other analyses that have pointed out that crustal heterogeneities in the continents result in diverse lithospheric thermal characteristics, which in turn results in higher complexity and variations of rheological behavior compared to oceanic domains with a thinner, more homogeneous crust.

Whistler mode exohiss are the structureless hiss waves observed outside the plasma pause with featured equatorward Poynting flux. An event of the amplification of exohiss as well as chorus waves was recorded by Van Allen Probes during the recovery phase of a weak geomagnetic storm. Amplitudes of both types of the waves showed a significant increase at the regions of electron density enhancements. It is found that the electrons resonant with exohiss and chorus showed moderate pitch angle anisotropies. The ratio of the number of electrons resonating with exohiss to total electron number presented in-phase correlation with density variations, which suggests that exohiss can be amplified due to electron density enhancement in terms of cyclotron instability. The calculation of linear growth rates further supports above conclusion. We suggest that exohiss waves have potential to become more significant due to the background plasma fluctuation.

This study presents the development and a systematic evaluation study of GNSS reflectometry wind speeds. After establishing a wind speed retrieval algorithm, UK TechDemoSat-1 (TDS-1) derived winds, from May 2015 to July 2017, are compared to the Advanced Scatterometer (ASCAT). ERA-Interim wind fields of the European Centre for Medium-range Weather Forecasts (ECMWF) and in situ observation from Tropical Atmosphere Ocean buoy array in the Pacific are taken as reference. One-year averaged TDS-1 global winds demonstrate small differences with ECMWF in a majority of areas as well as discuss under- and overestimations. The pioneering TDS-1 winds demonstrate a root-mean-squared error (RMSE) and bias of 2.77 and -0.33 m/s, which are comparable to the RMSE and bias derived by ASCAT winds, as large as 2.31 and 0.25 m/s, respectively. Using buoys measurements as reference, RMSE and bias of 2.23 and -0.03 m/s for TDS-1 as well as 1.40 and -0.68 m/s for ASCAT are obtained. Utilizing rain microwave-infrared estimates of the Tropical Rainfall Measuring Mission, rain-affected observation of both ASCAT and TDS-1 are collected and evaluated. Although ASCAT winds show a significant performance degradation resulting in an RMSE and bias of 3.16 and 1.03 m/s, respectively, during rain condition, TDS-1 shows a more reliable performance with an RMSE and bias of 2.94 and -0.21 m/s, respectively, which indicates the promising capability of GNSS forward scattering for wind retrievals during rain. A decrease in TDS-1-derived bistatic radar cross sections during rain events, at weak winds, is also demonstrated.

Representation of land-use and hydrologic interactions in respective models has traditionally been problematic. The use of static land-use in most hydrologic models or that of the use of simple hydrologic proxies in land-use change models call for more integrated approaches. The objective of this study is to assess whether dynamic feedback between land-use change and hydrology can (1) improve model performances, and/or (2) produce a more realistic quantification of ecosystem services. To test this, we coupled a land-use change model and a hydrologic mode. First, the land-use change and the hydrologic models were separately developed and calibrated. Then, the two models were dynamically coupled to exchange data at yearly time-steps. The approach is applied to a catchment in South Africa. Performance of coupled models when compared to the uncoupled models were marginal, but the coupled models excelled at the quantification of catchment ecosystem services more robustly.

From north to south, denudation rates from cosmogenic nuclides are similar to 10 t km(-2) yr(-1) at the arid Pan de Aziicar site, similar to 20 t km(2) yr(-1) at the semi-arid site of Santa Gracia, -60 t km(-2) yr(-1) at the Mediterranean climate site of La Campana, and similar to 30 t km(-2) yr(-1) at the humid site of Nahuelbuta. A and B horizons increase in thickness and elemental depletion or enrichment increases from north (similar to 26 degrees S) to south (similar to 38 degrees S) in these horizons. Differences in the degree of chemical weathering, quantified by the chemical depletion fraction (CDF), are significant only between the arid and sparsely vegetated site and the other three sites. Differences in the CDF between the sites, and elemental depletion within the sites are sometimes smaller than the variations induced by the bedrock heterogeneity. Microbial abundances (bacteria and archaea) in saprolite substantially increase from the arid to the semi-arid sites. With this study, we provide a comprehensive dataset characterizing the Critical Zone geochemistry in the Chilean Coastal Cordillera. This dataset confirms climatic controls on weathering and denudation rates and provides prerequisites to quantify the role of biota in future studies.

Modelers can improve a model by addressing the causes for the model errors (data errors and structural errors). This leads to implementing model enhancements (MEs), for example, meteorological data based on more monitoring stations, improved calibration data, and/or modifications in process formulations. However, deciding on which MEs to implement remains a matter of expert knowledge. After implementing multiple MEs, any improvement in model performance is not easily attributed, especially when considering different objectives or aspects of this improvement (e.g., better dynamics vs. reduced bias). We present an approach for comparing the effect of multiple MEs based on real observations and considering multiple objectives (MMEMO). A stepwise selection approach and structured plots help to address the multidimensionality of the problem. Tailored analyses allow a differentiated view on the effect of MEs and their interactions. MMEMO is applied to a case study employing the mesoscale hydro-sedimentological model WASA-SED for the Mediterranean-mountainous Isabena catchment, northeast Spain. The investigated seven MEs show diverse effects: some MEs (e.g., rainfall data) cause improvements for most objectives, while other MEs (e.g., land use data) only affect a few objectives or even decrease model performance. Interaction of MEs was observed for roughly half of the MEs, confirming the need to address them in the analysis. Calibration and increasing the temporal resolution showed by far stronger impact than any of the other MEs. The proposed framework can be adopted in other studies to analyze the effect of MEs and, thus, facilitate the identification and implementation of the most promising MEs for comparable cases.

It is of major interest to estimate the feedback of arctic ecosystems to the global warming we expect in upcoming decades. The speed of this response is driven by the potential of species to migrate, tracking their climate optimum. For this, sessile plants have to produce and disperse seeds to newly available habitats, and pollination of ovules is needed for the seeds to be viable. These two processes are also the vectors that pass genetic information through a population. A restricted exchange among subpopulations might lead to a maladapted population due to diversity losses. Hence, a realistic implementation of these dispersal processes into a simulation model would allow an assessment of the importance of diversity for the migration of plant species in various environments worldwide. To date, dynamic global vegetation models have been optimized for a global application and overestimate the migration of biome shifts in currently warming temperatures. We hypothesize that this is caused by neglecting important fine-scale processes, which are necessary to estimate realistic vegetation trajectories. Recently, we built and parameterized a simulation model LAVESI for larches that dominate the latitudinal treelines in the northernmost areas of Siberia. In this study, we updated the vegetation model by including seed and pollen dispersal driven by wind speed and direction. The seed dispersal is modelled as a ballistic flight, and for the pollination of ovules of seeds produced, we implemented a wind-determined and distance-dependent probability distribution function using a von Mises distribution to select the pollen donor. A local sensitivity analysis of both processes supported the robustness of the model's results to the parameterization, although it highlighted the importance of recruitment and seed dispersal traits for migration rates. This individual-based and spatially explicit implementation of both dispersal processes makes it easily feasible to inherit plant traits and genetic information to assess the impact of migration processes on the genetics. Finally, we suggest how the final model can be applied to substantially help in unveiling the important drivers of migration dynamics and, with this, guide the improvement of recent global vegetation models.

Mechanical and/or chemical removal of material from the subsurface may generate large subsurface cavities, the destabilisation of which can lead to ground collapse and the formation of sinkholes. Numerical simulation of the interaction of cavity growth, host material deformation and overburden collapse is desirable to better understand the sinkhole hazard but is a challenging task due to the involved high strains and material discontinuities. Here, we present 2-D distinct element method numerical simulations of cavity growth and sinkhole development. Firstly, we simulate cavity formation by quasi-static, stepwise removal of material in a single growing zone of an arbitrary geometry and depth. We benchmark this approach against analytical and boundary element method models of a deep void space in a linear elastic material. Secondly, we explore the effects of properties of different uniform materials on cavity stability and sinkhole development. We perform simulated biaxial tests to calibrate macroscopic geotechnical parameters of three model materials representative of those in which sinkholes develop at the Dead Sea shoreline: mud, alluvium and salt. We show that weak materials do not support large cavities, leading to gradual sagging or suffusion-style subsidence. Strong materials support quasi-stable to stable cavities, the overburdens of which may fail suddenly in a caprock or bedrock collapse style. Thirdly, we examine the consequences of layered arrangements of weak and strong materials. We find that these are more susceptible to sinkhole collapse than uniform materials not only due to a lower integrated strength of the overburden but also due to an inhibition of stabilising stress arching. Finally, we compare our model sinkhole geometries to observations at the Ghor Al-Haditha sinkhole site in Jordan. Sinkhole depth / diameter ratios of 0.15 in mud, 0.37 in alluvium and 0.33 in salt are reproduced successfully in the calibrated model materials. The model results suggest that the observed distribution of sinkhole depth / diameter values in each material type may partly reflect sinkhole growth trends.

For the first time, a rain signature in Global Navigation Satellite System Reflectometry (GNSS-R) observations is demonstrated. Based on the argument that the forward quasi-specular scattering relies upon surface gravity waves with lengths larger than several wavelengths of the reflected signal, a commonly made conclusion is that the scatterometric GNSS-R measurements are not sensitive to the surface small-scale roughness generated by raindrops impinging on the ocean surface. On the contrary, this study presents an evidence that the bistatic radar cross section sigma(0) derived from TechDemoSat-1 data is reduced due to rain at weak winds, lower than approximate to 6 m/s. The decrease is as large as approximate to 0.7 dB at the wind speed of 3 m/s due to a precipitation of 0-2 mm/hr. The simulations based on the recently published scattering theory provide a plausible explanation for this phenomenon which potentially enables the GNSS-R technique to detect precipitation over oceans at low winds.

The phi(ev) is calculated from high-resolution discharge and precipitation data for several rain events with a cumulative precipitation P-cum ranging from less than 5mm to more than 80 mm. Because of the high uncertainty of phi(ev) associated with the hydrograph separation method, phi(ev) is calculated with several methods, including graphical methods, digital filters and a tracer-based method. The results indicate that the hydrological response depends on (theta) over bar (ini): during dry conditions phi(ev) is consistently below 0.1, even for events with high and intense precipitation. Above a threshold of (theta) over bar (ini) = 34 vol % phi(ev) can reach values up to 0.99 but there is a high scatter. Some variability can be explained with a weak correlation of phi(ev) with P-cum and rain intensity, but a considerable part of the variability remains unexplained. It is concluded that threshold-based methods can be helpful to prevent overestimation of the hydrological response during dry catchment conditions. The impact of soil moisture on the hydrological response during wet catchment conditions, however, is still insufficiently understood and cannot be generalized based on the present results.

The lateral and vertical temperature distribution in Oman is so far only poorly understood, particularly in the area between Muscat and the Batinah coast, which is the area of this study and which is composed of Cenozoic sediments developed as part of a foreland basin of the Makran Thrust Zone. Temperature logs (T-logs) were run and physical rock properties of the sediments were analyzed to understand the temperature distribution, thermal and hydraulic properties, and heat-transport processes within the sedimentary cover of northern Oman. An advective component is evident in the otherwise conduction-dominated geothermal play system, and is caused by both topography and density driven flow. Calculated temperature gradients (T-gradients) in two wells that represent conductive conditions are 18.7 and 19.5 °C km−1, corresponding to about 70–90 °C at 2000–3000 m depth. This indicates a geothermal potential that can be used for energy intensive applications like cooling or water desalinization. Sedimentation in the foreland basin was initiated after the obduction of the Semail Ophiolite in the late Campanian, and reflects the complex history of alternating periods of transgressive and regressive sequences with erosion of the Oman Mountains. Thermal and hydraulic parameters were analyzed of the basin’s heterogeneous clastic and carbonate sedimentary sequence. Surface heat-flow values of 46.4 and 47.9 mW m−2 were calculated from the T-logs and calculated thermal conductivity values in two wells. The results of this study serve as a starting point for assessing different geothermal applications that may be suitable for northern Oman.

Due to global warming, the problem of assessing water resources and their vulnerability to climate drivers in the Arctic region has become a focus in the recent years. This study is aimed at investigating three lumped hydrological models to predict daily runoff of large-scale Arctic basins in the case of substantial data scarcity. All models were driven only by meteorological forcing reanalysis dataset without any additional information about landscape, soil, or vegetation cover properties of the studied basins. Model parameter regionalization based on transferring the whole parameter set showed good efficiency for predictions in ungauged basins. We run a blind test of the proposed methodology for ensemble runoff predictions on five sub-basins, for which only monthly observations were available, and obtained promising results for current water resources assessment for a broad domain of ungauged basins in the Russian Arctic.