Refine
Year of publication
- 2017 (298) (remove)
Document Type
- Article (229)
- Doctoral Thesis (28)
- Postprint (20)
- Other (11)
- Review (5)
- Part of a Book (2)
- Bachelor Thesis (1)
- Conference Proceeding (1)
- Master's Thesis (1)
Keywords
- Holocene (9)
- climate change (5)
- Arctic (3)
- Chile (3)
- Climate change (3)
- Erosion (3)
- Palaeoclimate (3)
- Tibetan Plateau (3)
- change detection (3)
- earthquake (3)
Institute
- Institut für Geowissenschaften (298) (remove)
Anthropogenically amplified erosion leads to increased fine-grained sediment input into the fluvial system in the 15.000 km2 Kharaa River catchment in northern Mongolia and constitutes a major stressing factor for the aquatic ecosystem. This study uniquely combines the application of intensive monitoring, source fingerprinting and catchment modelling techniques to allow for the comparison of the credibility and accuracy of each single method. High-resolution discharge data were used in combination with daily suspended solid measurements to calculate the suspended sediment budget and compare it with estimations of the sediment budget model SedNet. The comparison of both techniques showed that the development of an overall sediment budget with SedNet was possible, yielding results in the same order of magnitude (20.3 kt a- 1 and 16.2 kt a- 1).
Radionuclide sediment tracing, using Be-7, Cs-137 and Pb-210 was applied to differentiate sediment sources for particles < 10μm from hillslope and riverbank erosion and showed that riverbank erosion generates 74.5% of the suspended sediment load, whereas surface erosion contributes 21.7% and gully erosion only 3.8%. The contribution of the single subcatchments of the Kharaa to the suspended sediment load was assessed based on their variation in geochemical composition (e.g. in Ti, Sn, Mo, Mn, As, Sr, B, U, Ca and Sb). These variations were used for sediment source discrimination with geochemical composite fingerprints based on Genetic Algorithm driven Discriminant Function Analysis, the Kruskal–Wallis H-test and Principal Component Analysis. The contributions of the individual sub-catchment varied from 6.4% to 36.2%, generally showing higher contributions from the sub-catchments in the middle, rather than the upstream portions of the study area.
The results indicate that river bank erosion generated by existing grazing practices of livestock is the main cause for elevated fine sediment input. Actions towards the protection of the headwaters and the stabilization of the river banks within the middle reaches were identified as the highest priority. Deforestation and by lodging and forest fires should be prevented to avoid increased hillslope erosion in the mountainous areas. Mining activities are of minor importance for the overall catchment sediment load but can constitute locally important point sources for particular heavy metals in the fluvial system.
The characteristics of an aboveground cosmic-ray neutron sensor (CRNS) are evaluated for monitoring a mountain snowpack in the Austrian Alps from March 2014 to June 2016. Neutron counts were compared to continuous point-scale snow depth (SD) and snow-water-equivalent (SWE) measurements from an automatic weather station with a maximum SWE of 600 mm (April 2014). Several spatially distributed Terrestrial Laser Scanning (TLS)-based SD and SWE maps were additionally used. A strong nonlinear correlation is found for both SD and SWE. The representative footprint of the CRNS is in the range of 230-270 m. In contrast to previous studies suggesting signal saturation at around 100 mm of SWE, no complete signal saturation was observed. These results imply that CRNS could be transferred into an unprecedented method for continuous detection of spatially averaged SD and SWE for alpine snowpacks, though with sensitivity decreasing with increasing SWE. While initially different functions were found for accumulation and melting season conditions, this could be resolved by accounting for a limited measurement depth. This depth limit is in the range of 200 mm of SWE for dense snowpacks with high liquid water contents and associated snow density values around 450 kg m(-3) and above. In contrast to prior studies with shallow snowpacks, interannual transferability of the results is very high regardless of presnowfall soil moisture conditions. This underlines the unexpectedly high potential of CRNS to close the gap between point-scale measurements, hydrological models, and remote sensing of the cryosphere in alpine terrain.
Both Alpine and Mediterranean areas are considered sensitive to so-called global change, considered as the combination of climate and land use changes. All panels on climate evolution predict future scenarios of increasing frequency and magnitude of floods which are likely to lead to huge geomorphic adjustments of river channels so major metamorphosis of fluvial systems is expected as a result of global change. Such pressures are likely to give rise to major ecological and economic changes and challenges that governments need to address as a matter of priority. Changes in river flow regimes associated with global change are therefore ushering in a new era, where there is a critical need to evaluate hydro-geomorphological hazards from headwaters to lowland areas (flooding can be not just a problem related to being under the water). A key question is how our understanding of these hazards associated with global change can be improved; improvement has to come from integrated research which includes the climatological and physical conditions that could influence the hydrology and sediment generation and hence the conveyance of water and sediments (including the river’s capacity, i.e. amount of sediment, and competence, i.e. channel deformation) and the vulnerabilities and economic repercussions of changing hydrological hazards (including the evaluation of the hydro-geomorphological risks too).
Within this framework, the purpose of this international symposium is to bring together researchers from several disciplines as hydrology, fluvial geomorphology, hydraulic engineering, environmental science, geography, economy (and any other related discipline) to discuss the effects of global change over the river system in relation with floods. The symposium is organized by means of invited talks given by prominent experts, oral lectures, poster sessions and discussion sessions for each individual topic; it will try to improve our understanding of how rivers are likely to evolve as a result of global change and hence address the associated hazards of that fluvial environmental change concerning flooding.
Four main topics are going to be addressed:
- Modelling global change (i.e. climate and land-use) at relevant spatial (regional, local) and temporal (from the long-term to the single-event) scales.
- Measuring and modelling river floods from the hydrological, sediment transport (both suspended and bedload) and channel morphology points of view at different spatial (from the catchment to the reach) and temporal (from the long-term to the single-event) scales.
- Evaluation and assessment of current and future river flooding hazards and risks in a global change perspective.
- Catchment management to face river floods in a changing world.
We are very pleased to welcome you to Potsdam. We hope you will enjoy your participation at the International Symposium on the Effects of Global Change on Floods, Fluvial Geomorphology and Related Hazards in Mountainous Rivers and have an exciting and profitable experience. Finally, we would like to thank all speakers, participants, supporters, and sponsors for their contributions that for sure will make of this event a very remarkable and fruitful meeting. We acknowledge the valuable support of the European Commission (Marie Curie Intra-European Fellowship, Project ‘‘Floodhazards’’, PIEF-GA-2013-622468, Seventh EU Framework Programme) and the Deutschen Forschungsgemeinschaft (Research Training Group “Natural Hazards and Risks in a Changing World” (NatRiskChange; GRK 2043/1) as the symposium would not have been possible without their help. Without your cooperation, this symposium would not be either possible or successful.
Climate or land use?
(2017)
This study intends to contribute to the ongoing discussion on whether land use and land cover changes (LULC) or climate trends have the major influence on the observed increase of flood magnitudes in the Sahel. A simulation-based approach is used for attributing the observed trends to the postulated drivers. For this purpose, the ecohydrological model SWIM (Soil and Water Integrated Model) with a new, dynamic LULC module was set up for the Sahelian part of the Niger River until Niamey, including the main tributaries Sirba and Goroul. The model was driven with observed, reanalyzed climate and LULC data for the years 1950–2009. In order to quantify the shares of influence, one simulation was carried out with constant land cover as of 1950, and one including LULC. As quantitative measure, the gradients of the simulated trends were compared to the observed trend. The modeling studies showed that for the Sirba River only the simulation which included LULC was able to reproduce the observed trend. The simulation without LULC showed a positive trend for flood magnitudes, but underestimated the trend significantly. For the Goroul River and the local flood of the Niger River at Niamey, the simulations were only partly able to reproduce the observed trend. In conclusion, the new LULC module enabled some first quantitative insights into the relative influence of LULC and climatic changes. For the Sirba catchment, the results imply that LULC and climatic changes contribute in roughly equal shares to the observed increase in flooding. For the other parts of the subcatchment, the results are less clear but show, that climatic changes and LULC are drivers for the flood increase; however their shares cannot be quantified. Based on these modeling results, we argue for a two-pillar adaptation strategy to reduce current and future flood risk: Flood mitigation for reducing LULC-induced flood increase, and flood adaptation for a general reduction of flood vulnerability.
The characteristics of a landscape pose essential factors for hydrological processes. Therefore, an adequate representation of the landscape of a catchment in hydrological models is vital. However, many of such models exist differing, amongst others, in spatial concept and discretisation. The latter constitutes an essential pre-processing step, for which many different algorithms along with numerous software implementations exist. In that context, existing solutions are often model specific, commercial, or depend on commercial back-end software, and allow only a limited or no workflow automation at all.
Consequently, a new package for the scientific software and scripting environment R, called lumpR, was developed. lumpR employs an algorithm for hillslope-based landscape discretisation directed to large-scale application via a hierarchical multi-scale approach. The package addresses existing limitations as it is free and open source, easily extendible to other hydrological models, and the workflow can be fully automated. Moreover, it is user-friendly as the direct coupling to a GIS allows for immediate visual inspection and manual adjustment. Sufficient control is furthermore retained via parameter specification and the option to include expert knowledge. Conversely, completely automatic operation also allows for extensive analysis of aspects related to landscape discretisation.
In a case study, the application of the package is presented. A sensitivity analysis of the most important discretisation parameters demonstrates its efficient workflow automation. Considering multiple streamflow metrics, the employed model proved reasonably robust to the discretisation parameters. However, parameters determining the sizes of subbasins and hillslopes proved to be more important than the others, including the number of representative hillslopes, the number of attributes employed for the lumping algorithm, and the number of sub-discretisations of the representative hillslopes.
The Limpopo Basin in southern Africa is prone to droughts which affect the livelihood of millions of people in South Africa, Botswana, Zimbabwe and Mozambique. Seasonal drought early warning is thus vital for the whole region. In this study, the predictability of hydrological droughts during the main runoff period from December to May is assessed using statistical approaches. Three methods (multiple linear models, artificial neural networks, random forest regression trees) are compared in terms of their ability to forecast streamflow with up to 12 months of lead time. The following four main findings result from the study. 1. There are stations in the basin at which standardised streamflow is predictable with lead times up to 12 months. The results show high inter-station differences of forecast skill but reach a coefficient of determination as high as 0.73 (cross validated). 2. A large range of potential predictors is considered in this study, comprising well-established climate indices, customised teleconnection indices derived from sea surface temperatures and antecedent streamflow as a proxy of catchment conditions. El Nino and customised indices, representing sea surface temperature in the Atlantic and Indian oceans, prove to be important teleconnection predictors for the region. Antecedent streamflow is a strong predictor in small catchments (with median 42% explained variance), whereas teleconnections exert a stronger influence in large catchments. 3. Multiple linear models show the best forecast skill in this study and the greatest robustness compared to artificial neural networks and random forest regression trees, despite their capabilities to represent nonlinear relationships. 4. Employed in early warning, the models can be used to forecast a specific drought level. Even if the coefficient of determination is low, the forecast models have a skill better than a climatological forecast, which is shown by analysis of receiver operating characteristics (ROCs). Seasonal statistical forecasts in the Limpopo show promising results, and thus it is recommended to employ them as complementary to existing forecasts in order to strengthen preparedness for droughts.
Flood events can be expressed by a variety of characteristics such as flood magnitude and extent, event duration or incurred loss. Flood estimation and management may benefit from understanding how the different flood characteristics relate to the hydrological catchment conditions preceding the event and to the meteorological conditions throughout the event. In this study, we therefore propose a methodology to investigate the hydro-meteorological controls on different flood characteristics, based on the simulation of the complete flood risk chain from the flood triggering precipitation event, through runoff generation in the catchment, flood routing and possible inundation in the river system and floodplains to flood loss. Conditional cumulative distribution functions and regression tree analysis delineate the seasonal varying flood processes and indicate that the effect of the hydrological pre-conditions, i.e. soil moisture patterns, and of the meteorological conditions, i.e. weather patterns, depends on the considered flood characteristic. The methodology is exemplified for the Elbe catchment. In this catchment, the length of the build-up period, the event duration and the number of gauges undergoing at least a 10-year flood are governed by weather patterns. The affected length and the number of gauges undergoing at least a 2-year flood are however governed by soil moisture patterns. In case of flood severity and loss, the controlling factor is less pronounced. Severity is slightly governed by soil moisture patterns whereas loss is slightly governed by weather patterns. The study highlights that flood magnitude and extent arise from different flood generation processes and concludes that soil moisture patterns as well as weather patterns are not only beneficial to inform on possible flood occurrence but also on the involved flood processes and resulting flood characteristics.
Long-term estimates of the dissolution of silicate rock are generally derived from a range of isotopic proxies, such as the radiogenic strontium isotope ratio (Sr-87/Sr-86), which are preserved in sediment archives. For these systems to fairly represent silicate weathering, the changes in isotopic ratios in terrestrial surface waters should correspond to changes in the overall silicate dissolution. This assumes that the silicate mineral phases that act as sources of a given isotope dissolve at a rate that is proportional to the overall silicate weathering. Bedrock landsliding exhumes large quantities of fresh rock for weathering in transient storage, and rapid weathering in these deposits is controlled primarily by dissolution of the most reactive phases. In this study, we test the hypothesis that preferential weathering of these labile minerals can decouple the dissolution of strontium sources from the actual silicate weathering rates in the rapidly eroding Western Southern Alps (WSA) of New Zealand. We find that rapid dissolution of relatively radiogenic calcite and biotite in landslides leads to high local fluxes in strontium with isotopic ratios that offer no clear discrimination between sources. These higher fluxes of radiogenic strontium are in contrast to silicate weathering rates in landslides that are not systematically elevated. On a mountain belt scale, radiogenic strontium fluxes are not coupled to volumes of recent landslides in large (>100 km(2)) catchments, but silicate weathering fluxes are. Such decoupling is likely due first to the broad variability in the strontium content of carbonate minerals, and second to the combination of radiogenic strontium released from both biotite and carbonate in recent landslides. This study supports previous work suggesting the limited utility of strontium isotopes as a system to study silicate weathering in the WSA. Crucially however, in settings where bedrock landsliding is a dominant erosive process there is potential for both random and systematic bias in isotope proxies if the most reactive phases exposed for dissolution by landslides disproportionately contribute to the proxy of choice. This clearly suggests that the isotopic composition of marine Sr is a proxy for periods of rapid mountain uplift and erosion rather than for the associated enhanced silicate weathering. (C) 2017 Elsevier Ltd. All rights reserved.
Erosion and tectonic uplift are widely thought to be coupled through feedbacks involving orographic precipitation, relief development, and crustal weakening. In many orogenic systems, it can be difficult to distinguish whether true feedbacks exist, or whether observed features are a consequence of tectonic forcing. To help elucidate these interactions, we examine Gongga Shan, a 7556 m peak on the eastern margin of the Tibetan Plateau where cosmogenic Be-10 basin-wide erosion rates reach >5 mm/yr, defining a region of localized rapid erosion associated with a restraining bend in the left-lateral Xianshuihe Fault. Erosion rates are consistent with topography, thermochronometry, and geodetic data, suggesting a stable pattern of uplift and exhumation over at least the past 2-3 My. Transpression along the Xianshuihe Fault, orographically enhanced precipitation, thermally weakened crust, and substantial local relief all developed independently in the Gongga region and existed there prior to the uplift of Gongga Shan. However, only where all of these conditions are present do the observed topographic and erosional extremes exist, and their relative timing indicates that these conditions are not a consequence of rapid uplift. We conclude that their collocation at 3-4 Ma set into motion a series of feedbacks between erosion and uplift that has resulted in the exceptionally high topography and rapid erosion rates observed today.
Mountain channels can be strongly coupled with adjacent hillslopes, exchanging both mass and energy. However, hypotheses of the underlying cause and effect relations are based on indirect observations that do not resolve the mechanics of channel-hillslope coupling at the process scale. Here we present direct observational data of a coupled channel-hillslope system in the catchment area of the Erlenbach, a mountain stream in Switzerland. A slow-moving landslide flanking this alpine stream failed after a flood had eroded an alluvial step in the channel at its base, representing evidence for an upsystem link in channel-hillslope coupling. Progressive accumulation of landslide debris in the channel resulted in a renewed step, stabilizing the hillslope and restoring the channel step in a downsystem link. Thus, upsystem and downsystem coupling mechanisms are joined in a negative feedback cycle. In this cycle, debuttressing and rebuttressing due to channel bed erosion and alluviation are the dominant controls on hillslope stability. Based on an order of magnitude estimate it is plausible that the observed feedback mechanism is a relevant process in the production of coarse (>2 mm) sediment in the Erlenbach.
Permafrost is a distinct feature of the terrestrial Arctic and is vulnerable to climate warming. Permafrost degrades in different ways, including deepening of a seasonally unfrozen surface and localized but rapid development of deep thaw features. Pleistocene ice-rich permafrost with syngenetic ice-wedges, termed Yedoma deposits, are widespread in Siberia, Alaska, and Yukon, Canada and may be especially prone to rapid-thaw processes. Freeze-locked organic matter in such deposits can be re-mobilized on short time-scales and contribute to a carbon-cycle climate feedback. Here we synthesize the characteristics and vulnerability of Yedoma deposits by synthesizing studies on the Yedoma origin and the associated organic carbon pool. We suggest that Yedoma deposits accumulated under periglacial weathering, transport, and deposition dynamics in non-glaciated regions during the late Pleistocene until the beginning of late glacial warming. The deposits formed due to a combination of aeolian, colluvial, nival, and alluvial deposition and simultaneous ground ice accumulation. We found up to 130 gigatons organic carbon in Yedoma, parts of which are well-preserved and available for fast decomposition after thaw. Based on incubation experiments, up to 10% of the Yedoma carbon is considered especially decomposable and may be released upon thaw. The substantial amount of ground ice in Yedoma makes it highly vulnerable to disturbances such as thermokarst and thermo-erosion processes. Mobilization of permafrost carbon is expected to increase under future climate warming. Our synthesis results underline the need of accounting for Yedoma carbon stocks in next generation Earth-System-Models for a more complete representation of the permafrost-carbon feedback.
Presentation and Analysis of a Worldwide Database of Earthquake-Induced Landslide Inventories
(2017)
Earthquake-induced landslide (EQIL) inventories are essential tools to extend our knowledge of the relationship between earthquakes and the landslides they can trigger. Regrettably, such inventories are difficult to generate and therefore scarce, and the available ones differ in terms of their quality and level of completeness. Moreover, access to existing EQIL inventories is currently difficult because there is no centralized database. To address these issues, we compiled EQIL inventories from around the globe based on an extensive literature study. The database contains information on 363 landslide-triggering earthquakes and includes 66 digital landslide inventories. To make these data openly available, we created a repository to host the digital inventories that we have permission to redistribute through the U.S. Geological Survey ScienceBase platform. It can grow over time as more authors contribute their inventories. We analyze the distribution of EQIL events by time period and location, more specifically breaking down the distribution by continent, country, and mountain region. Additionally, we analyze frequency distributions of EQIL characteristics, such as the approximate area affected by landslides, total number of landslides, maximum distance from fault rupture zone, and distance from epicenter when the fault plane location is unknown. For the available digital EQIL inventories, we examine the underlying characteristics of landslide size, topographic slope, roughness, local relief, distance to streams, peak ground acceleration, peak ground velocity, and Modified Mercalli Intensity. Also, we present an evaluation system to help users assess the suitability of the available inventories for different types of EQIL studies and model development.
Prediction of the area affected by earthquake-induced landsliding based on seismological parameters
(2017)
We present an analytical, seismologically consistent expression for the surface area of the region within which most landslides triggered by an earthquake are located (landslide distribution area). This expression is based on scaling laws relating seismic moment, source depth, and focal mechanism with ground shaking and fault rupture length and assumes a globally constant threshold of acceleration for onset of systematic mass wasting. The seismological assumptions are identical to those recently used to propose a seismologically consistent expression for the total volume and area of landslides triggered by an earthquake. To test the accuracy of the model we gathered geophysical information and estimates of the landslide distribution area for 83 earthquakes. To reduce uncertainties and inconsistencies in the estimation of the landslide distribution area, we propose an objective definition based on the shortest distance from the seismic wave emission line containing 95% of the total landslide area. Without any empirical calibration the model explains 56% of the variance in our dataset, and predicts 35 to 49 out of 83 cases within a factor of 2, depending on how we account for uncertainties on the seismic source depth. For most cases with comprehensive landslide inventories we show that our prediction compares well with the smallest region around the fault containing 95% of the total landslide area. Aspects ignored by the model that could explain the residuals include local variations of the threshold of acceleration and processes modulating the surface ground shaking, such as the distribution of seismic energy release on the fault plane, the dynamic stress drop, and rupture directivity. Nevertheless, its simplicity and first-order accuracy suggest that the model can yield plausible and useful estimates of the landslide distribution area in near-real time, with earthquake parameters issued by standard detection routines.
Soil properties show high heterogeneity at different spatial scales and their correct characterization remains a crucial challenge over large areas. The aim of the study is to quantify the impact of different types of uncertainties that arise from the unresolved soil spatial variability on simulated hydrological states and fluxes. Three perturbation methods are presented for the characterization of uncertainties in soil properties. The methods are applied on the soil map of the upper Neckar catchment (Germany), as an example. The uncertainties are propagated through the distributed mesoscale hydrological model (mHM) to assess the impact on the simulated states and fluxes. The model outputs are analysed by aggregating the results at different spatial and temporal scales. These results show that the impact of the different uncertainties introduced in the original soil map is equivalent when the simulated model outputs are analysed at the model grid resolution (i.e. 500 m). However, several differences are identified by aggregating states and fluxes at different spatial scales (by subcatchments of different sizes or coarsening the grid resolution). Streamflow is only sensitive to the perturbation of long spatial structures while distributed states and fluxes (e.g. soil moisture and groundwater recharge) are only sensitive to the local noise introduced to the original soil properties. A clear identification of the temporal and spatial scale for which finer-resolution soil information is (or is not) relevant is unlikely to be universal. However, the comparison of the impacts on the different hydrological components can be used to prioritize the model improvements in specific applications, either by collecting new measurements or by calibration and data assimilation approaches. In conclusion, the study underlines the importance of a correct characterization of uncertainty in soil properties. With that, soil maps with additional information regarding the unresolved soil spatial variability would provide strong support to hydrological modelling applications.
Seismic-hazard assessment is of great importance within the field of engineering seismology. Nowadays, it is common practice to define future seismic demands using probabilistic seismic-hazard analysis (PSHA). Often it is neither obvious nor transparent how PSHA responds to changes in its inputs. In addition, PSHA relies on many uncertain inputs. Sensitivity analysis (SA) is concerned with the assessment and quantification of how changes in the model inputs affect the model response and how input uncertainties influence the distribution of the model response. Sensitivity studies are challenging primarily for computational reasons; hence, the development of efficient methods is of major importance. Powerful local (deterministic) methods widely used in other fields can make SA feasible, even for complex models with a large number of inputs; for example, automatic/algorithmic differentiation (AD)-based adjoint methods. Recently developed derivative-based global sensitivity measures can combine the advantages of such local SA methods with efficient sampling strategies facilitating quantitative global sensitivity analysis (GSA) for complex models. In our study, we propose and implement exactly this combination. It allows an upper bounding of the sensitivities involved in PSHA globally and, therefore, an identification of the noninfluential and the most important uncertain inputs. To the best of our knowledge, it is the first time that derivative-based GSA measures are combined with AD in practice. In addition, we show that first-order uncertainty propagation using the delta method can give satisfactory approximations of global sensitivity measures and allow a rough characterization of the model output distribution in the case of PSHA. An illustrative example is shown for the suggested derivative-based GSA of a PSHA that uses stochastic ground-motion simulations.
The aim of this paper is to investigate the ability of various site-condition proxies (SCPs) to reduce ground-motion aleatory variability and evaluate how SCPs capture nonlinearity site effects. The SCPs used here are time-averaged shear-wave velocity in the top 30 m (V-S30), the topographical slope (slope), the fundamental resonance frequency (f(0)) and the depth beyond which V-s exceeds 800 m/s (H800). We considered first the performance of each SCP taken alone and then the combined performance of the 6 SCP pairs [V-S30-f(0)], [V-S30-H-800], [f(0)-slope], [H-800-slope], [V-S30-slope] and [f(0)-H-800]. This analysis is performed using a neural network approach including a random effect applied on a KiK-net subset for derivation of ground-motion prediction equations setting the relationship between various ground-motion parameters such as peak ground acceleration, peak ground velocity and pseudo-spectral acceleration PSA (T), and Mw, RJB, focal depth and SCPs. While the choice of SCP is found to have almost no impact on the median groundmotion prediction, it does impact the level of aleatory uncertainty. VS30 is found to perform the best of single proxies at short periods (T < 0.6 s), while f(0) and H-800 perform better at longer periods; considering SCP pairs leads to significant improvements, with particular emphasis on [V-S30-H-800] and [f(0)-slope] pairs. The results also indicate significant nonlinearity on the site terms for soft sites and that the most relevant loading parameter for characterising nonlinear site response is the "stiff" spectral ordinate at the considered period.
In hydrological models, parameters are used to represent the time-invariant characteristics of catchments and to capture different aspects of hydrological response. Hence, model parameters need to be identified based on their role in controlling the hydrological behaviour. For the identification of meaningful parameter values, multiple and complementary performance criteria are used that compare modelled and measured discharge time series. The reliability of the identification of hydrologically meaningful model parameter values depends on how distinctly a model parameter can be assigned to one of the performance criteria.& para;& para;To investigate this, we introduce the new concept of connective strength between model parameters and performance criteria. The connective strength assesses the intensity in the interrelationship between model parameters and performance criteria in a bijective way. In our analysis of connective strength, model simulations are carried out based on a latin hypercube sampling. Ten performance criteria including Nash-Sutcliffe efficiency (NSE), Kling-Gupta efficiency (KGE) and its three components (alpha, beta and r) as well as RSR (the ratio of the root mean square error to the standard deviation) for different segments of the flow duration curve (FDC) are calculated.& para;& para;With a joint analysis of two regression tree (RT) approaches, we derive how a model parameter is connected to different performance criteria. At first, RTs are constructed using each performance criterion as the target variable to detect the most relevant model parameters for each performance criterion. Secondly, RTs are constructed using each parameter as the target variable to detect which performance criteria are impacted by changes in the values of one distinct model parameter. Based on this, appropriate performance criteria are identified for each model parameter.& para;& para;In this study, a high bijective connective strength between model parameters and performance criteria is found for low- and mid-flow conditions. Moreover, the RT analyses emphasise the benefit of an individual analysis of the three components of KGE and of the FDC segments. Furthermore, the RT analyses highlight under which conditions these performance criteria provide insights into precise parameter identification. Our results show that separate performance criteria are required to identify dominant parameters on low- and mid-flow conditions, whilst the number of required performance criteria for high flows increases with increasing process complexity in the catchment. Overall, the analysis of the connective strength between model parameters and performance criteria using RTs contribute to a more realistic handling of parameters and performance criteria in hydrological modelling.
Reliable flood risk analyses, including the estimation of damage, are an important prerequisite for efficient risk management. However, not much is known about flood damage processes affecting companies. Thus, we conduct a flood damage assessment of companies in Germany with regard to two aspects. First, we identify relevant damage-influencing variables. Second, we assess the prediction performance of the developed damage models with respect to the gain by using an increasing amount of training data and a sector-specific evaluation of the data. Random forests are trained with data from two postevent surveys after flood events occurring in the years 2002 and 2013. For a sector-specific consideration, the data set is split into four subsets corresponding to the manufacturing, commercial, financial, and service sectors. Further, separate models are derived for three different company assets: buildings, equipment, and goods and stock. Calculated variable importance values reveal different variable sets relevant for the damage estimation, indicating significant differences in the damage process for various company sectors and assets. With an increasing number of data used to build the models, prediction errors decrease. Yet the effect is rather small and seems to saturate for a data set size of several hundred observations. In contrast, the prediction improvement achieved by a sector-specific consideration is more distinct, especially for damage to equipment and goods and stock. Consequently, sector-specific data acquisition and a consideration of sector-specific company characteristics in future flood damage assessments is expected to improve the model performance more than a mere increase in data.
Voluminous magmatism during the South Atlantic opening has been considered as a classical example for plume related continental breakup. We present a study of the crustal structure around Walvis Ridge, near the intersection with the African margin. Two wide-angle seismic profiles were acquired. One is oriented NNW–SSE, following the continent–ocean transition and crossing Walvis Ridge. A second amphibious profile runs NW–SE from the Angola Basin into continental Namibia. At the continent–ocean boundary (COB) the mafic crust beneath Walvis Ridge is up to 33 km thick, with a pronounced high-velocity lower crustal body. Towards the south there is a smooth transition to 20–25 km thick crust underlying the COB in the Walvis Basin, with a similar velocity structure, indicating a gabbroic lower crust with associated cumulates at the base. The northern boundary of Walvis Ridge towards the Angola Basin shows a sudden change to oceanic crust only 4–6 km thick, coincident with the projection of the Florianopolis Fracture Zone, one of the most prominent tectonic features of the South Atlantic ocean basin. In the amphibious profile the COB is defined by a sharp transition from oceanic to rifted continental crust, with a magmatic overprint landward of the intersection of Walvis Ridge with the Namibian margin. The continental crust beneath the Congo Craton is 40 km thick, shoaling to 35 km further SE. The velocity models show that massive high-velocity gabbroic intrusives are restricted to a narrow zone directly underneath Walvis Ridge and the COB in the south. This distribution of rift-related magmatism is not easily reconciled with models of continental breakup following the establishment of a large, axially symmetric plume in the Earth's mantle. Rift-related lithospheric stretching and associated transform faulting play an overriding role in locating magmatism, dividing the margin in a magma-dominated southern and an essentially amagmatic northern segment.
Vast areas of the terrestrial Subarctic and Arctic are underlain by permafrost. Landscape evolution is therefore largely controlled by climate-driven periglacial processes. The response of the frozen ground to late Quaternary warm and cold stages is preserved in permafrost sequences, and deducible by multi-proxy palaeoenvironmental approaches. Here, we analyse radiocarbon-dated mid-Wisconsin Interstadial and Holocene lacustrine deposits preserved in the Kit-1 pingo permafrost sequence combined with water and surface sediment samples from nine modern water bodies on Seward Peninsula (NW Alaska) to reconstruct thermokarst dynamics and determine major abiotic factors that controlled the aquatic ecosystem variability. Our methods comprise taxonomical diatom analyses as well as Detrended Correspondence Analysis (DCA) and Redundancy Analysis (RDA). Our results show, that the fossil diatom record reflects thermokarst lake succession since about 42 C-14 kyr BP. Different thermolcarst lake stages during the mid-Wisconsin Interstadial, the late Wisconsin and the early Holocene are mirrored by changes in diatom abundance, diversity, and ecology. We interpret the taxonomical changes in the fossil diatom assemblages in combination with both modern diatom data from surrounding ponds and existing micropalaeontological, sedimentological and mineralogical data from the pingo sequence. A diatom based quantitative reconstruction of lake water pH indicates changing lake environments during mid-Wisconsin to early Holocene stages. Mineralogical analyses indicate presence of tephra fallout and its impact on fossil diatom communities. Our comparison of modern and fossil diatom communities shows the highest floristic similarity of modern polygon ponds to the corresponding initial (shallow water) development stages of thermolcarst lakes. We conclude, that mid-Wisconsin thermokarst processes in the study area could establish during relatively warm interstadial climate conditions accompanied by increased precipitation due to approaching coasts, while still high continentality and hence high seasonal temperature gradients led to warm summers in the central part of Beringia. (C) 2017 Elsevier B.V. All rights reserved.
We investigated a well-dated sediment section of a palaeolake situated in the coastal zone of Shikotan Island (Lesser Kurils) for organic sediment-geochemistry and biotic components (diatoms, chironomids, pollen) in order to provide a reconstruction of the palaeoenvironmental changes and palaeo-events (tsunamis, sea-level fluctuations and landslides) in Holocene. During the ca 8000 years of sedimentation the changes in organic sediment-geochemistry and in composition of the diatoms and chironomids as well as the shifts in composition of terrestrial vegetation suggest that the period until ca 5800 cal yr BP was characterized by a warm and humid climate (corresponds to middle Holocene optimum) with climate cooling thereafter. A warm period reconstructed from ca 900 to at least ca 580 cal yr BP corresponds to a transition to a Nara-Heian-Kamakura warm stage and can be correlated to a Medieval Warm Period. After 580 cal yr PB, the lake gradually dried out and climatic signals could not be obtained from the declining lacustrine biological communities, but the increasing role of spruce and disappearance of the oak from the vegetation give evidences of the climate cooling that can be correlated with the LIA. The marine regression stages at the investigated site are identified for ca 6200-5900 (at the end of the middle Holocene transgression), ca 5500-5100 (Middle Jomon regression or Kemigawa regression), and ca 1070-360 cal yr BP (at the end of Heian transgression). The lithological structure of sediments and the diatom compositions give evidences for the multiple tsunami events of different strengths in the Island. Most remarkable of them can be dated at around ca 7000, 6460, 5750, 4800, 950 cal yr BP. The new results help to understand the Holocene environmental history of the Southern Kurils as a part of the Kuril-Kamchatka and Aleutian Marginal Sea-Island Arc Systems in the North-Western Pacific region.
Tree stands in the boreal treeline ecotone are, in addition to climate change, impacted by disturbances such as fire, water-related disturbances and logging. We aim to understand how these disturbances affect growth, age structure, and spatial patterns of larch stands in the north-eastern Siberian treeline ecotone (lower Kolyma River region), an insufficiently researched region. Stand structure of Larix cajanderi Mayr was studied at seven sites impacted by disturbances. Maximum tree age ranged from 44 to 300 years. Young to medium-aged stands had, independent of disturbance type, the highest stand densities with over 4000 larch trees per ha. These sites also had the highest growth rates for tree height and stem diameter. Overall lowest stand densities were found in a polygonal field at the northern end of the study area, with larches growing in distinct " tree islands". At all sites, saplings are significantly clustered. Differences in fire severity led to contrasting stand structures with respect to tree, recruit, and overall stand densities. While a low severity fire resulted in low-density stands with high proportions of small and young larches, high severity fires resulted in high-density stands with high proportions of big trees. At waterdisturbed sites, stand structure varied between waterlogged and drained sites and latitude. These mixed effects of climate and disturbance make it difficult to predict future stand characteristics and the treeline position.
This article presents a review of the current data on the level of paleolimnological knowledge about lakes in the Russian part of the northern Eurasia. The results of investigation of the northwestern European part of Russia as the best paleolimnologically studied sector of the Russian north is presented in detail. The conditions of lacustrine sedimentation at the boundary between the Late Pleistocene and Holocene and the role of different external factors in formation of their chemical composition, including active volcanic activity and possible large meteorite impacts, are also discussed. The results of major paleoclimatic and paleoecological reconstructions in northern Siberia are presented. Particular attention is given to the databases of abiotic and biotic parameters of lake ecosystems as an important basis for quantitative reconstructions of climatic and ecological changes in the Late Pleistocene and Holocene. Keywords: paleolimnology, lakes, bottom sediments, northern.
Adaptation to flood risk
(2017)
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.
Trends in precipitation over Germany and the Rhine basin related to changes in weather patterns
(2017)
Precipitation as the central meteorological feature for agriculture, water security, and human well-being amongst others, has gained special attention ever since. Lack of precipitation may have devastating effects such as crop failure and water scarcity. Abundance of precipitation, on the other hand, may as well result in hazardous events such as flooding and again crop failure. Thus, great effort has been spent on tracking changes in precipitation and relating them to underlying processes. Particularly in the face of global warming and given the link between temperature and atmospheric water holding capacity, research is needed to understand the effect of climate change on precipitation.
The present work aims at understanding past changes in precipitation and other meteorological variables. Trends were detected for various time periods and related to associated changes in large-scale atmospheric circulation. The results derived in this thesis may be used as the foundation for attributing changes in floods to climate change. Assumptions needed for the downscaling of large-scale circulation model output to local climate stations are tested and verified here.
In a first step, changes in precipitation over Germany were detected, focussing not only on precipitation totals, but also on properties of the statistical distribution, transition probabilities as a measure for wet/dry spells, and extreme precipitation events.
Shifting the spatial focus to the Rhine catchment as one of the major water lifelines of Europe and the largest river basin in Germany, detected trends in precipitation and other meteorological variables were analysed in relation to states of an ``optimal'' weather pattern classification. The weather pattern classification was developed seeking the best skill in explaining the variance of local climate variables.
The last question addressed whether observed changes in local climate variables are attributable to changes in the frequency of weather patterns or rather to changes within the patterns itself. A common assumption for a downscaling approach using weather patterns and a stochastic weather generator is that climate change is expressed only as a changed occurrence of patterns with the pattern properties remaining constant. This assumption was validated and the ability of the latest generation of general circulation models to reproduce the weather patterns was evaluated.
% Paper 1
Precipitation changes in Germany in the period 1951-2006 can be summarised briefly as negative in summer and positive in all other seasons. Different precipitation characteristics confirm the trends in total precipitation: while winter mean and extreme precipitation have increased, wet spells tend to be longer as well (expressed as increased probability for a wet day followed by another wet day). For summer the opposite was observed: reduced total precipitation, supported by decreasing mean and extreme precipitation and reflected in an increasing length of dry spells.
Apart from this general summary for the whole of Germany, the spatial distribution within the country is much more differentiated. Increases in winter precipitation are most pronounced in the north-west and south-east of Germany, while precipitation increases are highest in the west for spring and in the south for autumn. Decreasing summer precipitation was observed in most regions of Germany, with particular focus on the south and west.
The seasonal picture, however, was again differently represented in the contributing months, e.g.\ increasing autumn precipitation in the south of Germany is formed by strong trends in the south-west in October and in the south-east in November. These results emphasise the high spatial and temporal organisation of precipitation changes.
% Paper 2
The next step towards attributing precipitation trends to changes in large-scale atmospheric patterns was the derivation of a weather pattern classification that sufficiently stratifies the local climate variables under investigation. Focussing on temperature, radiation, and humidity in addition to precipitation, a classification based on mean sea level pressure, near-surface temperature, and specific humidity was found to have the best skill in explaining the variance of the local variables. A rather high number of 40 patterns was selected, allowing typical pressure patterns being assigned to specific seasons by the associated temperature patterns. While the skill in explaining precipitation variance is rather low, better skill was achieved for radiation and, of course, temperature.
Most of the recent GCMs from the CMIP5 ensemble were found to reproduce these weather patterns sufficiently well in terms of frequency, seasonality, and persistence.
% Paper 3
Finally, the weather patterns were analysed for trends in pattern frequency, seasonality, persistence, and trends in pattern-specific precipitation and temperature. To overcome uncertainties in trend detection resulting from the selected time period, all possible periods in 1901-2010 with a minimum length of 31 years were considered. Thus, the assumption of a constant link between patterns and local weather was tested rigorously. This assumption was found to hold true only partly. While changes in temperature are mainly attributable to changes in pattern frequency, for precipitation a substantial amount of change was detected within individual patterns.
Magnitude and even sign of trends depend highly on the selected time period. The frequency of certain patterns is related to the long-term variability of large-scale circulation modes.
Changes in precipitation were found to be heterogeneous not only in space, but also in time - statements on trends are only valid for the specific time period under investigation. While some part of the trends can be attributed to changes in the large-scale circulation, distinct changes were found within single weather patterns as well.
The results emphasise the need to analyse multiple periods for thorough trend detection wherever possible and add some note of caution to the application of downscaling approaches based on weather patterns, as they might misinterpret the effect of climate change due to neglecting within-type trends.
The Cotentin Peninsula (Normandy, France) displays sequences of marine terraces and rasas, the latter being wide Late Cenozoic coastal erosion surfaces, that are typical of Western European coasts in Portugal, Spain, France and southern England. Remote sensing imagery and field mapping enabled reappraisal of the Cotentin coastal sequences. From bottom to top, the N Cotentin sequence includes four previously recognized Pleistocene marine terraces (T1 to T4) at elevations <40 m as well as four higher and older rasas (R1 to R4) reaching 200 +/- 5 m in elevation. Low-standing marine terraces are not observed in the central part of the Peninsula and a limited number of terraces are described to the south. The high-standing rasas are widespread all over the peninsula. Such strandline distributions reveal major changes during the Late Cenozoic. Progressive uplift of an irregular sea-floor led to subaerial exposure of bathymetric highs that were carved into rocky platforms, rasas and marine terraces. Eventually, five main islands coalesced and connected to the mainland to the south to form the Cotentin Peninsula. On the basis of previous dating of the last interglacial maximum terrace (i.e. Marine Isotopic Stage, MIS 5e), sequential morphostratigraphy and modelling, we have reappraised uplift rates and derived: (i) mean Upper Pleistocene (i.e. since MIS 5e similar to 122 +/- 6 ka, i.e. kilo annum) apparent uplift rates of 0.04 +/- 0.01 mm/yr, (ii) mean Middle Pleistocene eustasy-corrected uplift rates of 0.09 +/- 0.03 mm/yr, and (iii) low mean Pleistocene uplift rates of 0.01 mm/yr. Extrapolations of these slow rates combined with geological evidence implies that the formation of the sequences from the Cotentin Peninsula occurred between 3 Ma (Pliocene) and 15 Ma (Miocene), which cannot be narrowed down further without additional research. Along the coasts of Western Europe, sequences of marine terraces and rasas are widespread (169 preserve the MIS Se benchmark). In Spain, Portugal, S England and other parts of western France, the sequences morphostratigraphy is very similar to that of Cotentin. The onset of such Western European sequences occurred during the Miocene (e.g. Spain) or Pliocene (e.g. Portugal). We interpret this Neogene-Quaternary coastal uplift as a symptom of the increasing lithospheric compression that accompanies Cenozoic orogenies. (C) 2017 Elsevier B.V. All rights reserved.
Pollen-based quantitative reconstructions of past climate variables is a standard palaeoclimatic approach. Despite knowing that the spatial extent of the calibration-set affects the reconstruction result, guidance is lacking as to how to determine a suitable spatial extent of the pollen-climate calibration-set. In this study, past mean annual precipitation (P-ann) during the Holocene (since 11.5 cal ka BP) is reconstructed repeatedly for pollen records from Qinghai Lake (36.7 degrees N, 100.5 degrees E; north-east Tibetan Plateau), Gonghai Lake (38.9 degrees N, 112.2 degrees E; north China) and Sihailongwan Lake (42.3 degrees N, 126.6 degrees E; north-east China) using calibration-sets of varying spatial extents extracted from the modern pollen dataset of China and Mongolia (2559 sampling sites and 168 pollen taxa in total). Results indicate that the spatial extent of the calibration-set has a strong impact on model performance, analogue quality and reconstruction diagnostics (absolute value, range, trend, optimum). Generally, these effects are stronger with the modern analogue technique (MAT) than with weighted averaging partial least squares (WA-PLS). With respect to fossil spectra from northern China, the spatial extent of calibration-sets should be restricted to radii between ca. 1000 and 1500 km because small-scale calibration-sets (<800 km radius) will likely fail to include enough spatial variation in the modern pollen assemblages to reflect the temporal range shifts during the Holocene, while too broad a scale calibration-set (>1500 km radius) will include taxa with very different pollen-climate relationships. (C) 2017 Elsevier Ltd. All rights reserved.
The Indian Summer Monsoon (ISM) is one of the most important climate systems, whose variability and driving mechanisms are of broad interest for academic and societal communities. Here, we present a well-dated high-resolution pollen analysis from a 4.82-m long sediment core taken from Basomtso, in the southeastern Tibetan Plateau (TP), which depicts the regional climate changes of the past millennium. Our results show that subalpine coniferous forest was dominant around Basomtso from ca. 867 to ca. 750 cal. yr BP, indicating a warm and semi-humid climate. The timberline in the study area significantly decreased from ca. 750 to ca.100 cal. yr BP, and a cold climate, corresponding to the Little Ice Age (LIA) prevailed. Since ca. 100 cal. yr BP, the vegetation type changed to forest-meadow with rising temperatures and moisture. Ordination analysis reveals that the migration of vegetation was dominated by regional temperatures and then by moisture. Further comparisons between the Basomtso pollen record and the regional temperature reconstructions underscore the relevance of the Basomtso record from the southeastern TP for regional and global climatologies. Our pollen based moisture reconstruction demonstrates the strong multicentennial-scale link to ISM variability, providing solid evidence for the increase of monsoonal strengths over the past four centuries. Spectral analysis indicates the potential influence of solar forcing. However, a closer relationship has been observed between multicentennial ISM variations and Indian Ocean sea surface temperature anomalies (SSTs), suggesting that the variations in monsoonal precipitation over the southeastern TP are probably driven by the Indian Ocean Dipole on the multicentennial scale. (C) 2017 Elsevier Ltd. All rights reserved.
Early agriculture can be detected in palaeovegetation records, but quantification of the relative importance of climate and land use in influencing regional vegetation composition since the onset of agriculture is a topic that is rarely addressed. We present a novel approach that combines pollen-based REVEALS estimates of plant cover with climate, anthropogenic land-cover and dynamic vegetation modelling results. This is used to quantify the relative impacts of land use and climate on Holocene vegetation at a sub-continental scale, i.e. northern and western Europe north of the Alps. We use redundancy analysis and variation partitioning to quantify the percentage of variation in vegetation composition explained by the climate and land-use variables, and Monte Carlo permutation tests to assess the statistical significance of each variable. We further use a similarity index to combine pollen based REVEALS estimates with climate-driven dynamic vegetation modelling results. The overall results indicate that climate is the major driver of vegetation when the Holocene is considered as a whole and at the sub-continental scale, although land use is important regionally. Four critical phases of land-use effects on vegetation are identified. The first phase (from 7000 to 6500 BP) corresponds to the early impacts on vegetation of farming and Neolithic forest clearance and to the dominance of climate as a driver of vegetation change. During the second phase (from 4500 to 4000 BP), land use becomes a major control of vegetation. Climate is still the principal driver, although its influence decreases gradually. The third phase (from 2000 to 1500 BP) is characterised by the continued role of climate on vegetation as a consequence of late-Holocene climate shifts and specific climate events that influence vegetation as well as land use. The last phase (from 500 to 350 BP) shows an acceleration of vegetation changes, in particular during the last century, caused by new farming practices and forestry in response to population growth and industrialization. This is a unique signature of anthropogenic impact within the Holocene but European vegetation remains climatically sensitive and thus may continue to respond to ongoing climate change. (C) 2017 Elsevier Ltd. All rights reserved.
Climatic and limnological changes at Lake Karakul (Tajikistan) during the last similar to 29 cal ka
(2017)
We present results of analyses on a sediment core from Lake Karakul, located in the eastern Pamir Mountains, Tajikistan. The core spans the last similar to 29 cal ka. We investigated and assessed processes internal and external to the lake to infer changes in past moisture availability. Among the variables used to infer lake-external processes, high values of grain-size end-member (EM) 3 (wide grain-size distribution that reflects fluvial input) and high Sr/Rb and Zr/Rb ratios (coinciding with coarse grain sizes), are indicative of moister conditions. High values in EM1, EM2 (peaks of small grain sizes that reflect long-distance dust transport or fine, glacially derived clastic input) and TiO2 (terrigenous input) are thought to reflect greater influence of dry air masses, most likely of Westerly origin. High input of dust from distant sources, beginning before the Last Glacial Maximum (LGM) and continuing to the late glacial, reflects the influence of dry Westerlies, whereas peaks in fluvial input suggest increased moisture availability. The early to early-middle Holocene is characterised by coarse mean grain sizes, indicating constant, high fluvial input and moister conditions in the region. A steady increase in terrigenous dust and a decrease in fluvial input from 6.6 cal ka BP onwards points to the Westerlies as the predominant atmospheric circulation through to present, and marks a return to drier and even arid conditions in the area. Proxies for productivity (TOC, TOC/TN, TOCBr), redox potential (Fe/Mn) and changes in the endogenic carbonate precipitation (TIC, delta(18) OCarb) indicate changes within the lake. Low productivity characterised the lake from the late Pleistocene until 6.6 cal ka BP, and increased rapidly afterwards. Lake level remained low until the LGM, but water depth increased to a maximum during the late glacial and remained high into the early Holocene. Subsequently, the water level decreased to its present stage. Today the lake system is mainly climatically controlled, but the depositional regime is also driven by internal limnogeological processes.
Arctic and alpine treelines worldwide differ in their reactions to climate change. A northward advance of or densification within the treeline ecotone will likely influence climate-vegetation feedback mechanisms. In our study, which was conducted in the Taimyr Depression in the North Siberian Lowlands, w present a combined field-and model-based approach helping us to better understand the population processes involved in the responses of the whole treeline ecotone, spanning from closed forest to single-tree tundra, to climate warming. Using information on stand structure, tree age, and seed quality and quantity from seven sites, we investigate effects of intra-specific competition and seed availability on the specific impact of recent climate warming on larch stands. Field data show that tree density is highest in the forest-tundra, and average tree size decreases from closed forest to single-tree tundra. Age-structure analyses indicate that the trees in the closed forest and forest-tundra have been present for at least similar to 240 yr. At all sites except the most southerly ones, past establishment is positively correlated with regional temperature increase. In the single-tree tundra, however, a change in growth form from krummholz to erect trees, beginning similar to 130 yr ago, rather than establishment date has been recorded. Seed mass decreases from south to north, while seed quantity increases. Simulations with LAVESI (Larix Vegetation Simulator) further suggest that relative density changes strongly in response to a warming signal in the forest-tundra while intra-specific competition limits densification in the closed forest and seed limitation hinders densification in the single-tree tundra. We find striking differences in strength and timing of responses to recent climate warming. While forest-tundra stands recently densified, recruitment is almost non-existent at the southern and northern end of the ecotone due to autecological processes. Palaeo-treelines may therefore be inappropriate to infer past temperature changes at a fine scale. Moreover, a lagged treeline response to past warming will, via feedback mechanisms, influence climate change in the future.
The concept of a Global Monsoon (GM) has been proposed based on modern precipitation observations, but its application over a wide range of temporal scales is still under debate. Here, we present a synthesis of 268 continental paleo-moisture records collected from monsoonal systems in the Eastern Hemisphere, including the East Asian Monsoon (EAsM), the Indian Monsoon (IM), the East African Monsoon (EAfM), and the Australian Monsoon (AuM) covering the last 18,000 years. The overall pattern of late Glacial to Holocene moisture change is consistent with those inferred from ice cores and marine records. With respect to the last 10,000 years (10 ka), i.e. a period that has high spatial coverage, a Fuzzy c-Means clustering analysis of the moisture index records together with "Xie-Beni" index reveals four clusters of our data set. The paleoclimatic meaning of each cluster is interpreted considering the temporal evolution and spatial distribution patterns. The major trend in the tropical AuM, EAfM, and IM regions is a gradual decrease in moisture conditions since the early Holocene. Moisture changes in the EAsM regions show maximum index values between 8 and 6 ka. However, records located in nearby subtropical areas, i.e. in regions not influenced by the intertropical convergence zone, show an opposite trend compared to the tropical monsoon regions (AuM, EAfM and IM), i.e. a gradual increase. Analyses of modern meteorological data reveal the same spatial patterns as in the paleoclimate records such that, in times of overall monsoon strengthening, lower precipitation rates are observed in the nearby subtropical areas. We explain this pattern as the effect of a strong monsoon circulation suppressing air uplift in nearby subtropical areas, and hence hindering precipitation. By analogy to the modern system, this would mean that during the early Holocene strong monsoon period, the intensified ascending airflows within the monsoon domains led to relatively weaker ascending or even descending airflows in the adjacent subtropical regions, resulting in a precipitation deficit compared to the late Holocene. Our conceptual model therefore integrates regionally contrasting moisture changes into the Global Monsoon hypothesis. (C) 2017 Elsevier Ltd. All rights reserved.
Retrogressive thaw slumps (RTSs) are among the most active landforms in the Arctic; their number has increased significantly over the past decades. While processes initiating discrete RTSs are well identified, the major terrain controls on the development of coastal RTSs at a regional scale are not yet defined. Our research reveals the main geomorphic factors that determine the development of RTSs along a 238km segment of the Yukon Coast, Canada. We (1) show the current extent of RTSs, (2) ascertain the factors controlling their activity and initiation, and (3) explain the spatial differences in the density and areal coverage of RTSs. We mapped and classified 287 RTSs using high-resolution satellite images acquired in 2011. We highlighted the main terrain controls over their development using univariate regression trees model. Coastal geomorphology influenced both the activity and initiation of RTSs: active RTSs and RTSs initiated after 1972 occurred primarily on terrains with slope angles greater than 3.9 degrees and 5.9 degrees, respectively. The density and areal coverage of RTSs were constrained by the volume and thickness of massive ice bodies. Differences in rates of coastal change along the coast did not affect the model. We infer that rates of coastal change averaged over a 39year period are unable to reflect the complex relationship between RTSs and coastline dynamics. We emphasize the need for large-scale studies of RTSs to evaluate their impact on the ecosystem and to measure their contribution to the global carbon budget. Plain Language Summary Retrogressive thaw slumps, henceforth slumps are a type of landslides that occur when permafrost thaws. Slumps are active landforms: they develop quickly and extend over several hectares. Satellite imagery allows to map such slumps over large areas. Our research shows where slumps develop along a 238 km segment of the Yukon Coast in Canada and explains which environments are most suitable for slump occurrence. We found that active and newly developed slumps were triggered where coastal slopes were greater than 3.9 degrees and 5.9 degrees, respectively. We explain that coastal erosion influences the development of slumps by modifying coastal slopes. We found that the highest density of slumps as well as the largest slumps occurred on terrains with high amounts of ice bodies in the ground. This study provides tools to better identify areas in the Arctic that are prone to slump development.
Radiocarbon and optically stimulated luminescence dating of sediments from Lake Karakul, Tajikistan
(2017)
Lake Karakul in the eastern Pamirs is a large and closed-basin lake in a partly glaciated catchment. Two parallel sediment cores were collected from 12 m water depth. The cores were correlated using XRF analysis and dated using radiocarbon and OSL techniques. The age results of the two dating methods are generally in agreement. The correlated composite core of 12.26 m length represents continuous accumulation of sediments in the lake basin since 31 ka. The lake reservoir effect (LRE) remained relatively constant over this period. High sediment accumulation rates (SedARs) were recorded before 23 ka and after 6.5 ka. The relatively close position of the coring location near the eastern shore of the lake implies that high SedARs resulted from low lake levels. Thus, high SedARs and lower lake levels before 23 ka probably reflect cold and dry climate conditions that inhibited the arrival of moist air at high elevation in the eastern Pamirs. Low lake levels after 6.5 ka were probably caused by declining temperatures after the warmer early Holocene, which had caused a reduction in water resources stored as snow, ice and frozen ground in the catchment. Low SedARs during 23-6.5 ka suggest increased lake levels in Lake Karakul. A short-lived increase of SedARs at 15 ka probably corresponds to the rapid melting of glaciers in the Karakul catchment during the Greenland Interstadial le, shortly after glaciers in the catchment had reached their maximum extents. The sediment cores from Lake Karakul represent an important climate archive with robust chronology for the last glacial interglacial cycle from Central Asia. (C) 2017 Elsevier B.V. All rights reserved.
We analyzed chlorophyll-a and Colored Dissolved Organic Matter (CDOM) dynamics from field measurements and assessed the potential of multispectral satellite data for retrieving water-quality parameters in three small surface reservoirs in the Brazilian semiarid region. More specifically, this work is comprised of: (i) analysis of Chl-a and trophic dynamics; (ii) characterization of CDOM; (iii) estimation of Chl-a and CDOM from OLI/Landsat-8 and RapidEye imagery. The monitoring lasted 20 months within a multi-year drought, which contributed to water-quality deterioration. Chl-a and trophic state analysis showed a highly eutrophic status for the perennial reservoir during the entire study period, while the non-perennial reservoirs ranged from oligotrophic to eutrophic, with changes associated with the first events of the rainy season. CDOM characterization suggests that the perennial reservoir is mostly influenced by autochthonous sources, while allochthonous sources dominate the non-perennial ones. Spectral-group classification assigned the perennial reservoir as a CDOM-moderate and highly eutrophic reservoir, whereas the non-perennial ones were assigned as CDOM-rich and oligotrophic-dystrophic reservoirs. The remote sensing initiative was partially successful: the Chl-a was best modelled using RapidEye for the perennial one; whereas CDOM performed best with Landsat-8 for non-perennial reservoirs. This investigation showed potential for retrieving water quality parameters in dry areas with small reservoirs.
The large variety of atmospheric circulation systems affecting the eastern Asian climate is reflected by the complex Asian vegetation distribution. Particularly in the transition zones of these circulation systems, vegetation is supposed to be very sensitive to climate change. Since proxy records are scarce, hitherto a mechanistic understanding of the past spatio-temporal climate-vegetation relationship is lacking. To assess the Holocene vegetation change and to obtain an ensemble of potential mid-Holocene biome distributions for eastern Asia, we forced the diagnostic biome model BIOME4 with climate anomalies of different transient Holocene climate simulations performed in coupled atmosphere-ocean(-vegetation) models. The simulated biome changes are compared with pollen-based biome records for different key regions. In all simulations, substantial biome shifts during the last 6000 years are confined to the high northern latitudes and the monsoon-westerly wind transition zone, but the temporal evolution and amplitude of change strongly depend on the climate forcing. Large parts of the southern tundra are replaced by taiga during the mid-Holocene due to a warmer growing season and the boreal treeline in northern Asia is shifted northward by approx. 4 degrees in the ensemble mean, ranging from 1.5 to 6 degrees in the individual simulations, respectively. This simulated treeline shift is in agreement with pollen-based reconstructions from northern Siberia. The desert fraction in the transition zone is reduced by 21% during the mid-Holocene compared to pre-industrial due to enhanced precipitation. The desert-steppe margin is shifted westward by 5 degrees (1-9 degrees in the individual simulations). The forest biomes are expanded north-westward by 2 degrees, ranging from 0 to 4 degrees in the single simulations. These results corroborate pollen-based reconstructions indicating an extended forest area in north-central China during the mid-Holocene. According to the model, the forest-to-non-forest and steppe-to-desert changes in the climate transition zones are spatially not uniform and not linear since the mid-Holocene.
Temporal and spatial stability of the vegetation climate relationship is a basic ecological assumption for pollen-based quantitative inferences of past climate change and for predicting future vegetation. We explore this assumption for the Holocene in eastern continental Asia (China, Mongolia). Boosted regression trees (BRT) between fossil pollen taxa percentages (Abies, Artemisia, Betula, Chenopodiaceae, Cyperaceae, Ephedra, Picea, Pinus, Poaceae and Quercus) and climate model outputs of mean annual precipitation (P-ann) and mean temperature of the warmest month (Mt(wa)) for 9 and 6 ka (ka = thousand years before present) were set up and results compared to those obtained from relating modern pollen to modern climate. Overall, our results reveal only slight temporal differences in the pollen climate relationships. Our analyses suggest that the importance of P-ann compared with Mt(wa) for taxa distribution is higher today than it was at 6 ka and 9 ka. In particular, the relevance of P-ann for Picea and Pinus increases and has become the main determinant. This change in the climate tree pollen relationship parallels a widespread tree pollen decrease in north-central China and the eastern Tibetan Plateau. We assume that this is at least partly related to vegetation climate disequilibrium originating from human impact. Increased atmospheric CO2 concentration may have permitted the expansion of moisture-loving herb taxa (Cyperaceae and Poaceae) during the late Holocene into arid/semi-arid areas. We furthermore find that the pollen climate relationship between north-central China and the eastern Tibetan Plateau is generally similar, but that regional differences are larger than temporal differences. In summary, vegetation climate relationships in China are generally stable in space and time, and pollen-based climate reconstructions can be applied to the Holocene. Regional differences imply the calibration-set should be restricted spatially.
We have analyzed the recently developed pan-European strong motion database, RESORCE-2012: spectral parameters, such as stress drop (stress parameter, Delta sigma), anelastic attenuation (Q), near surface attenuation (kappa(0)) and site amplification have been estimated from observed strong motion recordings. The selected dataset exhibits a bilinear distance-dependent Q model with average kappa(0) value 0.0308 s. Strong regional variations in inelastic attenuation were also observed: frequency-independent Q(0) of 1462 and 601 were estimated for Turkish and Italian data respectively. Due to the strong coupling between Q and kappa(0), the regional variations in Q have strong impact on the estimation of near surface attenuation kappa(0). kappa(0) was estimated as 0.0457 and 0.0261 s for Turkey and Italy respectively. Furthermore, a detailed analysis of the variability in estimated kappa(0) revealed significant within-station variability. The linear site amplification factors were constrained from residual analysis at each station and site-class type. Using the regional Q(0) model and a site-class specific kappa(0), seismic moments (M-0) and source corner frequencies f (c) were estimated from the site corrected empirical Fourier spectra. Delta sigma did not exhibit magnitude dependence. The median Delta sigma value was obtained as 5.75 and 5.65 MPa from inverted and database magnitudes respectively. A comparison of response spectra from the stochastic model (derived herein) with that from (regional) ground motion prediction equations (GMPEs) suggests that the presented seismological parameters can be used to represent the corresponding seismological attributes of the regional GMPEs in a host-to-target adjustment framework. The analysis presented herein can be considered as an update of that undertaken for the previous Euro-Mediterranean strong motion database presented by Edwards and Fah (Geophys J Int 194(2):1190-1202, 2013a).
Receiver functions (RF) have been used for several decades to study structures beneath seismic stations. Although most available stations are deployed on shore, the number of ocean bottom station (OBS) experiments has increased in recent years. Almost all OBSs have to deal with higher noise levels and a limited deployment time (approximate to 1year), resulting in a small number of usable records of teleseismic earthquakes. Here we use OBSs deployed as midaperture array in the deep ocean (4.5-5.5km water depth) of the eastern mid-Atlantic. We use evaluation criteria for OBS data and beamforming to enhance the quality of the RFs. Although some stations show reverberations caused by sedimentary cover, we are able to identify the Moho signal, indicating a normal thickness (5-8km) of oceanic crust. Observations at single stations with thin sediments (300-400m) indicate that a probable sharp lithosphere-asthenosphere boundary (LAB) might exist at a depth of approximate to 70-80km which is in line with LAB depth estimates for similar lithospheric ages in the Pacific. The mantle discontinuities at approximate to 410km and approximate to 660km are clearly identifiable. Their delay times are in agreement with PREM. Overall the usage of beam-formed earthquake recordings for OBS RF analysis is an excellent way to increase the signal quality and the number of usable events.
Two nuclear explosions were carried out by the Democratic People’s Republic of North Korea in January and September 2016. Epicenters were located close to those of the 2006, 2009, and 2013 previous explosions. We perform a seismological analysis of the 2016 events combining the analysis of full waveforms at regional distances and seismic array beams at teleseismic distances. We estimate the most relevant source parameters, such as source depth, moment release, and full moment tensor (MT). The best MT solution can be decomposed into an isotropic source, directly related with the explosion and an additional deviatoric term, likely due to near‐source interactions with topographic and/or underground facilities features. We additionally perform an accurate resolution test to assess source parameters uncertainties and trade‐offs. This analysis sheds light on source parameters inconsistencies among studies on previous shallow explosive sources. The resolution of the true MT is hindered by strong source parameters trade‐offs, so that a broad range of well‐fitting MT solutions can be found, spanning from a dominant positive isotropic term to a dominant negative vertical compensated linear vector dipole. The true mechanism can be discriminated by additionally modeling first‐motion polarities at seismic arrays at teleseismic distances. A comparative assessment of the 2016 explosion with earlier nuclear tests documents similar vertical waveforms but a significant increase of amplitude for the 2016 explosions, which proves that the 9 September 2016 was the largest nuclear explosion ever performed in North Korea with a magnitude Mw 4.9 and a shallow depth of less than 2 km, although there are no proofs of a fusion explosion. Modeling transversal component waveforms suggests variable size and orientation of the double‐couple components of the 2009, 2013, and 2016 sources.
The increasing numbers of recordings at individual sites allows quantification of empirical linear site-response adjustment factors (delta S2S(s)) from the ground motion prediction equation (GMPE) residuals. The delta S2S(s) are then used to linearly scale the ergodic GMPE predictions to obtain site-specific ground motion predictions in a partially non-ergodic Probabilistic Seismic Hazard Assessment (PSHA). To address key statistical and conceptual issues in the current practice, we introduce a novel empirical region-and site-specific PSHA methodology wherein, (1) site-to-site variability (phi(S2S)) is first estimated as a random-variance in a mixed-effects GMPE regression, (2) delta S2S(s) at new sites with strong motion are estimated using the a priori phi(S2S), and (3) the GMPE site-specific single-site aleatory variability sigma(ss,s) is replaced with a generic site-corrected aleatory variability sigma(0). Comparison of region- and site-specific hazard curves from our method against the traditional ergodic estimates at 225 sites in Europe and Middle East shows an approximate 50% difference in predicted ground motions over a range of hazard levels-a strong motivation to increase seismological monitoring of critical facilities and enrich regional ground motion data sets.
Site-Corrected Magnitude- and Region-Dependent Correlations of Horizontal Peak Spectral Amplitudes
(2017)
Empirical correlations of horizontal peak spectral amplitudes (PSA) are modeled using the total-residuals obtained in a ground motion prediction equation (GMPE) regression. Recent GMPEs moved toward partially non-ergodic region-and site-specific predictions, while the residual correlation models remained largely ergodic. Using mixed-effects regression, we decompose the total-residuals of a pan-European GMPE into between-event, between-site, and event-and-site corrected residuals to investigate the ergodicity in empirical PSA correlations. We first observed that the between-event correlations are magnitude-dependent, partially due to the differences in source spectra, and influence of stress-drop parameter on small and large events. Next, removing the between-site residuals from within-event residuals yields the event-and-site corrected residuals which are found to be region-dependent, possibly due to the regional differences in distance-decay of short period PSAs. Using our site-corrected magnitude- and region-dependent correlations, and the between-site residuals as empirical site-specific ground motion adjustments, we compute partially non-ergodic conditional mean spectra at four well-recorded sites in Europe and Middle Eastern regions.
Most of the deformation associated with the seismic cycle in subduction zones occurs offshore and has been therefore difficult to quantify with direct observations at millennial timescales. Here we study millennial deformation associated with an active splay-fault system in the Arauco Bay area off south central Chile. We describe hitherto unrecognized drowned shorelines using high-resolution multibeam bathymetry, geomorphic, sedimentologic, and paleontologic observations and quantify uplift rates using a Landscape Evolution Model. Along a margin-normal profile, uplift rates are 1.3m/ka near the edge of the continental shelf, 1.5m/ka at the emerged Santa Maria Island, -0.1m/ka at the center of the Arauco Bay, and 0.3m/ka in the mainland. The bathymetry images a complex pattern of folds and faults representing the surface expression of the crustal-scale Santa Maria splay-fault system. We modeled surface deformation using two different structural scenarios: deep-reaching normal faults and deep-reaching reverse faults with shallow extensional structures. Our preferred model comprises a blind reverse fault extending from 3km depth down to the plate interface at 16km that slips at a rate between 3.0 and 3.7m/ka. If all the splay-fault slip occurs during every great megathrust earthquake, with a recurrence of similar to 150-200years, the fault would slip similar to 0.5m per event, equivalent to a magnitude similar to 6.4 earthquake. However, if the splay-fault slips only with a megathrust earthquake every similar to 1000years, the fault would slip similar to 3.7m per event, equivalent to a magnitude similar to 7.5 earthquake.
A sediment core from Khatanga-12 Lake (Taimyr Peninsula, Krasnoyarsk krai) has been studied. The 131.5-cm-long core covers ca. 7100 years of sedimentation. Chironomid analysis, a qualitative reconstruction of the paleoenvironment in the region, and a quantitative reconstruction of variations of the mean July air temperature and in the water depth of the lake have been performed using Northern Russia chironomid-inferred mean July temperature models (Nazarova et al., 2008, 2011, 2015). Khatanga-12 Lake was formed during the Middle Holocene warming as a result of thermokarst processes. The development of the lake ecosystem at different stages of its development was influenced by climatic and cryolithogenic factors. The Middle Holocene warming, which occurred around 7100-6250 cal. years BP, activated thermokarst processes and resulted in the formation of the lake basin. Later, between 6250 and 4500 cal. years BP, a period of cooling took place, as is proved by chironomid analysis. The bottom sediments of the lake during this period were formed by erosion processes on the lake shores. The reconstructed conditions were close to the modern after 2500 cal. years BP.
Aim: Across the planet, grass-dominated biomes are experiencing shrub encroachment driven by atmospheric CO2 enrichment and land-use change. By altering resource structure and availability, shrub encroachment may have important impacts on vertebrate communities. We sought to determine the magnitude and variability of these effects across climatic gradients, continents, and taxa, and to learn whether shrub thinning restores the structure of vertebrate communities. Location: Worldwide. Time period: Contemporary. Major taxa studied: Terrestrial vertebrates. Methods: We estimated relationships between percentage shrub cover and the structure of terrestrial vertebrate communities (species richness, Shannon diversity and community abundance) in experimentally thinned and unmanipulated shrub-encroached grass-dominated biomes using systematic review and meta-analyses of 43 studies published from 1978 to 2016. We modelled the effects of continent, biome, mean annual precipitation, net primary productivity and the normalized difference vegetation index (NDVI) on the relationship between shrub cover and vertebrate community structure. Results: Species richness, Shannon diversity and total abundance had no consistent relationship with shrub encroachment and experimental thinning did not reverse encroachment effects on vertebrate communities. However, some effects of shrub encroachment on vertebrate communities differed with net primary productivity, amongst vertebrate groups, and across continents. Encroachment had negative effects on vertebrate diversity at low net primary productivity. Mammalian and herpetofaunal diversity decreased with shrub encroachment. Shrub encroachment also had negative effects on species richness and total abundance in Africa but positive effects in North America. Main conclusions: Biodiversity conservation and mitigation efforts responding to shrub encroachment should focus on low-productivity locations, on mammals and herpetofauna, and in Africa. However, targeted research in neglected regions such as central Asia and India will be needed to fill important gaps in our knowledge of shrub encroachment effects on vertebrates. Additionally, our findings provide an impetus for determining the mechanisms associated with changes in vertebrate diversity and abundance in shrub-encroached grass-dominated biomes.