Refine
Year of publication
- 2016 (291) (remove)
Document Type
- Article (221)
- Doctoral Thesis (37)
- Other (14)
- Postprint (12)
- Review (5)
- Habilitation Thesis (1)
- Report (1)
Is part of the Bibliography
- yes (291) (remove)
Keywords
- climate change (5)
- erosion (5)
- water balance (5)
- Pollen (4)
- August 2002 flood (3)
- Body waves (3)
- Central Europe (3)
- Floods Directive (3)
- Himalaya (3)
- Himalayas (3)
Institute
- Institut für Geowissenschaften (291) (remove)
We present a new, seismologically consistent expression for the total area and volume of populations of earthquake-triggered landslides. This model builds on a set of scaling relationships between key parameters, such as landslide spatial density, seismic ground acceleration, fault length, earthquake source depth, and seismic moment. To assess the model we have assembled and normalized a catalog of landslide inventories for 40 shallow, continental earthquakes. Low landscape steepness causes systematic overprediction of the total area and volume of landslides. When this effect is accounted for, the model predicts the total landslide volume of 63% of 40 cases to within a factor 2 of the volume estimated from observations (R-2 = 0.76). The prediction of total landslide area is also sensitive to the landscape steepness, but less so than the total volume, and it appears to be sensitive to controls on the landslide size-frequency distribution, and possibly the shaking duration. Some outliers are likely associated with exceptionally strong rock mass in the epicentral area, while others may be related to seismic source complexities ignored by the model. However, the close match between prediction and estimate for about two thirds of cases in our database suggests that rock mass strength is similar in many cases and that our simple seismic model is often adequate, despite the variety of lithologies and tectonic settings covered. This makes our expression suitable for integration into landscape evolution models and application to the anticipation or rapid assessment of secondary hazards associated with earthquakes.
Bedrock landslides, by excavating deep below saprolite-rock interfaces, create conditions for weathering in which all mineral phases in a lithology are initially unweathered within landslide deposits. As a result, the most labile phases dominate the weathering immediately after mobilisation and during a transient period of depletion. This mode of dissolution can strongly alter the overall output of solutes from catchments and their contribution to global chemical cycles if landslide-derived material is retained in catchments for extended periods after mass wasting.
The South American Andes are frequently exposed to intense rainfall events with varying moisture sources and precipitation-forming processes. In this study, we assess the spatiotemporal characteristics and geographical origins of rainfall over the South American continent. Using high-spatiotemporal resolution satellite data (TRMM 3B42 V7), we define four different types of rainfall events based on their (1) high magnitude, (2) long temporal extent, (3) large spatial extent, and (4) high magnitude, long temporal and large spatial extent combined. In a first step, we analyze the spatiotemporal characteristics of these events over the entire South American continent and integrate their impact for the main Andean hydrologic catchments. Our results indicate that events of type 1 make the overall highest contributions to total seasonal rainfall (up to 50%). However, each consecutive episode of the infrequent events of type 4 still accounts for up to 20% of total seasonal rainfall in the subtropical Argentinean plains. In a second step, we employ complex network theory to unravel possibly non-linear and long-ranged climatic linkages for these four event types on the high-elevation Altiplano-Puna Plateau as well as in the main river catchments along the foothills of the Andes. Our results suggest that one to two particularly large squall lines per season, originating from northern Brazil, indirectly trigger large, long-lasting thunderstorms on the Altiplano Plateau. In general, we observe that extreme rainfall in the catchments north of approximately 20 degrees S typically originates from the Amazon Basin, while extreme rainfall at the eastern Andean foothills south of 20 degrees S and the Puna Plateau originates from southeastern South America.
Non-linear time series analysis of precipitation events using regional climate networks for Germany
(2016)
Synchronous occurrences of heavy rainfall events and the study of their relation in time and space are of large socio-economical relevance, for instance for the agricultural and insurance sectors, but also for the general well-being of the population. In this study, the spatial synchronization structure is analyzed as a regional climate network constructed from precipitation event series. The similarity between event series is determined by the number of synchronous occurrences. We propose a novel standardization of this number that results in synchronization scores which are not biased by the number of events in the respective time series. Additionally, we introduce a new version of the network measure directionality that measures the spatial directionality of weighted links by also taking account of the effects of the spatial embedding of the network. This measure provides an estimate of heavy precipitation isochrones by pointing out directions along which rainfall events synchronize. We propose a climatological interpretation of this measure in terms of propagating fronts or event traces and confirm it for Germany by comparing our results to known atmospheric circulation patterns.
The analytical evaluation of diurnal temperature variation in riverbed sediments provides detailed information on exchange fluxes between rivers and groundwater. The underlying assumption of the stationary, one-dimensional vertical flow field is frequently violated in natural systems where subsurface water flow often has a significant horizontal component. In this paper, we present a new methodology for identifying the geometry of the subsurface flow field using vertical temperature profiles. The statistical analyses are based on model optimisation and selection and are used to evaluate the shape of vertical amplitude ratio profiles. The method was applied to multiple profiles measured around in-stream geomorphological structures in a losing reach of a gravel bed river. The predominant subsurface flow field was systematically categorised in purely vertical and horizontal (hyporheic, parafluvial) components. The results highlight that river groundwater exchange flux at the head, crest and tail of geomorphological structures significantly deviated from the one-dimensional vertical flow, due to a significant horizontal component. The geometry of the subsurface water flow depended on the position around the geomorphological structures and on the river level. The methodology presented in this paper features great potential for characterising the spatial patterns and temporal dynamics of complex subsurface flow geometries by using measured temperature time series in vertical profiles. (C) 2016 Elsevier B.V. All rights reserved.
Im Graduiertenkolleg NatRiskChange der Universität Potsdam und anderen Forschungseinrichtungen werden beobachtete sowie zukünftig mögliche Veränderungen von Naturgefahren untersucht. Teil des strukturierten Doktorandenprogramms sind sogenannte Task-Force-Einsätze, bei denen die Promovierende zeitlich begrenzt ein aktuelles Ereignis auswerten. Im Zuge dieser Aktivität wurde die Sturzflut vom 29.05.2016 in Braunsbach (Baden-Württemberg) untersucht.
In diesem Bericht werden erste Auswertungen zur Einordnung der Niederschläge, zu den hydrologischen und geomorphologischen Prozessen im Einzugsgebiet des Orlacher Bachs sowie zu den verursachten Schäden beleuchtet.
Die Region war Zentrum extremer Regenfälle in der Größenordnung von 100 mm innerhalb von 2 Stunden. Das 6 km² kleine Einzugsgebiet hat eine sehr schnelle Reaktionszeit, zumal bei vorgesättigtem Boden. Im steilen Bachtal haben mehrere kleinere und größere Hangrutschungen über 8000 m³ Geröll, Schutt und Schwemmholz in das Gewässer eingetragen und möglicherweise kurzzeitige Aufstauungen und Durchbrüche verursacht. Neben den großen Wassermengen mit einer Abflussspitze in einer Größenordnung von 100 m³/s hat gerade die Geschiebefracht zu großen Schäden an den Gebäuden entlang des Bachlaufs in Braunsbach geführt.
Carbon and nutrient cycling in kettle hole sediments depending on hydrological dynamics: a review
(2016)
Kettle holes as a specific group of isolated, small lentic freshwater systems (LFS) often are (i) hot spots of biogeochemical cycling and (ii) exposed to frequent sediment desiccation and rewetting. Their ecological functioning is greatly determined by immanent carbon and nutrient transformations. The objective of this review is to elucidate effects of a changing hydrological regime (i.e., dry-wet cycles) on carbon and nutrient cycling in kettle hole sediments. Generally, dry-wet cycles have the potential to increase C and N losses as well as P availability. However, their duration and frequency are important controlling factors regarding direction and intensity of biogeochemical and microbiological responses. To evaluate drought impacts on sediment carbon and nutrient cycling in detail requires the context of the LFS hydrological history. For example, frequent drought events induce physiological adaptation of exposed microbial communities and thus flatten metabolic responses, whereas rare events provoke unbalanced, strong microbial responses. Different potential of microbial resilience to drought stress can irretrievably change microbial communities and functional guilds, gearing cascades of functional responses. Hence, dry-wet events can shift the biogeochemical cycling of organic matter and nutrients to a new equilibrium, thus affecting the dynamic balance between carbon burial and mineralization in kettle holes.
Climatic controls on debris-flow activity and sediment aggradation: The Del Medio fan, NW Argentina
(2016)
In the Central Andes, several studies on alluvial terraces and valley fills have linked sediment aggradation to periods of enhanced sediment supply. However, debate continues over whether tectonic or climatic factors are most important in triggering the enhanced supply. The Del Medio catchment in the Humahuaca Basin (Eastern Cordillera, NW Argentina) is located within a transition zone between subhumid and arid climates and hosts the only active debris-flow fan within this intermontane valley. By combining Be-10 analyses of boulder and sediment samples within the Del Medio catchment, with regional morphometric measurements of nearby catchments, we identify the surface processes responsible for aggradation in the Del Medio fan and their likely triggers. We find that the fan surface has been shaped by debris flows and channel avulsions during the last 400 years. Among potential tectonic, climatic, and autogenic factors that might influence deposition, our analyses point to a combination of several favorable factors that drive aggradation. These are in particular the impact of occasional abundant rainfall on steep slopes in rock types prone to failure, located in a region characterized by relatively low rainfall amounts and limited transport capacity. These characteristics are primarily associated with the climatic transition zone between the humid foreland and the arid orogen interior, which creates an imbalance between sediment supply and sediment transfer. The conditions and processes that drive aggradation in the Del Medio catchment today may provide a modern analog for the conditions and processes that drove aggradation in other nearby tributaries in the past.
Fluvial fill terraces preserve sedimentary archives of landscape responses to climate change, typically over millennial timescales. In the Humahuaca Basin of NW Argentina (Eastern Cordillera, southern Central Andes), our 29 new optically stimulated luminescence ages of late Pleistocene fill terrace sediments demonstrate that the timing of past river aggradation occurred over different intervals on the western and eastern sides of the valley, despite their similar bedrock lithology, mean slopes, and precipitation. In the west, aggradation coincided with periods of increasing precipitation, while in the east, aggradation coincided with decreasing precipitation or more variable conditions. Erosion rates and grain size dependencies in our cosmogenic Be-10 analyses of modern and fill terrace sediments reveal an increased importance of landsliding compared to today on the west side during aggradation, but of similar importance during aggradation on the east side. Differences in the timing of aggradation and the Be-10 data likely result from differences in valley geometry, which causes sediment to be temporarily stored in perched basins on the east side. It appears as if periods of increasing precipitation triggered landslides throughout the region, which induced aggradation in the west, but blockage of the narrow bedrock gorges downstream from the perched basins in the east. As such, basin geometry and fluvial connectivity appear to strongly influence the timing of sediment movement through the system. For larger basins that integrate subbasins with differing geometries or degrees of connectivity (like Humahuaca), sedimentary responses to climate forcing are likely attenuated.
A link between chemical weathering and physical erosion exists at the catchment scale over a wide range of erosion rates(1,2). However, in mountain environments, where erosion rates are highest, weathering may be kinetically limited(3-5) and therefore decoupled from erosion. In active mountain belts, erosion is driven by bedrock landsliding(6) at rates that depend strongly on the occurrence of extreme rainfall or seismicity(7). Although landslides affect only a small proportion of the landscape, bedrock landsliding can promote the collection and slow percolation of surface runoff in highly fragmented rock debris and create favourable conditions for weathering. Here we show from analysis of surface water chemistry in the Southern Alps of New Zealand that weathering in bedrock landslides controls the variability in solute load of these mountain rivers. We find that systematic patterns in surface water chemistry are strongly associated with landslide occurrence at scales from a single hillslope to an entire mountain belt, and that landslides boost weathering rates and river solute loads over decades. We conclude that landslides couple erosion and weathering in fast-eroding uplands and, thus, mountain weathering is a stochastic process that is sensitive to climatic and tectonic controls on mass wasting processes.
This manuscript proposes a method to assess hydrological drought in semi-arid environments under high impoundment rate and applies it to the semi-arid Jaguaribe River basin in Brazil. It analyzes droughts (1) in the largest reservoir systems; (2) in the Upper Basin, considering 4744 reservoirs, 800 wells and almost 18,000 cisterns; and (3) in reservoirs of different sizes during multiyear droughts. Results show that the water demand is constrained in the basin; hydrological and meteorological droughts are often out of phase; there is a negative correlation between storage level and drought severity; and the small systems cannot cope with long-term droughts.
The Dead Sea region has faced substantial environmental challenges in recent decades, including water resource scarcity, similar to 1 m annual decreases in the water level, sinkhole development, ascending-brine freshwater pollution, and seismic disturbance risks. Natural processes are significantly affected by human interference as well as by climate change and tectonic developments over the long term. To get a deep understanding of processes and their interactions, innovative scientific approaches that integrate disciplinary research and education are required. The research project DESERVE (Helmholtz Virtual Institute Dead Sea Research Venue) addresses these challenges in an interdisciplinary approach that includes geophysics, hydrology, and meteorology. The project is implemented by a consortium of scientific institutions in neighboring countries of the Dead Sea (Israel, Jordan, Palestine Territories) and participating German Helmholtz Centres (KIT, GFZ, UFZ). A new monitoring network of meteorological, hydrological, and seismic/geodynamic stations has been established, and extensive field research and numerical simulations have been undertaken. For the first time, innovative measurement and modeling techniques have been applied to the extreme conditions of the Dead Sea and its surroundings. The preliminary results show the potential of these methods. First time ever performed eddy covariance measurements give insight into the governing factors of Dead Sea evaporation. High-resolution bathymetric investigations reveal a strong correlation between submarine springs and neo-tectonic patterns. Based on detailed studies of stratigraphy and borehole information, the extension of the subsurface drainage basin of the Dead Sea is now reliably estimated. Originality has been achieved in monitoring flash floods in an arid basin at its outlet and simultaneously in tributaries, supplemented by spatio-temporal rainfall data. Low-altitude, high resolution photogrammetry, allied to satellite image analysis and to geophysical surveys (e.g. shear-wave reflections) has enabled a more detailed characterization of sinkhole morphology and temporal development and the possible subsurface controls thereon. All the above listed efforts and scientific results take place with the interdisciplinary education of young scientists. They are invited to attend joint thematic workshops and winter schools as well as to participate in field experiments. (C) 2015 The Authors. Published by Elsevier B.V.
Oceanic lithospheric S-wave velocities from the analysis of P-wave polarization at the ocean floor
(2016)
Our knowledge of the absolute S-wave velocities of the oceanic lithosphere is mainly based on global surface wave tomography, local active seismic or compliance measurements using oceanic infragravity waves. The results of tomography give a rather smooth picture of the actual S-wave velocity structure and local measurements have limitations regarding the range of elastic parameters or the geometry of the measurement. Here, we use the P-wave polarization (apparent P-wave incidence angle) of teleseismic events to investigate the S-wave velocity structure of the oceanic crust and the upper tens of kilometres of the mantle beneath single stations. In this study, we present an up to our knowledge new relation of the apparent P-wave incidence angle at the ocean bottom dependent on the half-space S-wave velocity. We analyse the angle in different period ranges at ocean bottom stations (OBSs) to derive apparent S-wave velocity profiles. These profiles are dependent on the S-wave velocity as well as on the thickness of the layers in the subsurface. Consequently, their interpretation results in a set of equally valid models. We analyse the apparent P-wave incidence angles of an OBS data set which was collected in the Eastern Mid Atlantic. We are able to determine reasonable S-wave-velocity-depth models by a three-step quantitative modelling after a manual data quality control, although layer resonance sometimes influences the estimated apparent S-wave velocities. The apparent S-wave velocity profiles are well explained by an oceanic PREM model in which the upper part is replaced by four layers consisting of a water column, a sediment, a crust and a layer representing the uppermost mantle. The obtained sediment has a thickness between 0.3 and 0.9 km with S-wave velocities between 0.7 and 1.4 km s(-1). The estimated total crustal thickness varies between 4 and 10 km with S-wave velocities between 3.5 and 4.3 km s(-1). We find a slight increase of the total crustal thickness from similar to 5 to similar to 8 km towards the South in the direction of a major plate boundary, the Gloria Fault. The observed crustal thickening can be related with the known dominant compression in the vicinity of the fault. Furthermore, the resulting mantle S-wave velocities decrease from values around 5.5 to 4.5 km s(-1) towards the fault. This decrease is probably caused by serpentinization and indicates that the oceanic transform fault affects a broad region in the uppermost mantle. Conclusively, the presented method is useful for the estimation of the local S-wave velocity structure beneath ocean bottom seismic stations. It is easy to implement and consists of two main steps: (1) measurement of apparent P-wave incidence angles in different period ranges for real and synthetic data, and (2) comparison of the determined apparent S-wave velocities for real and synthetic data to estimate S-wave velocity-depth models.
Rapidly uplifting coastlines are frequently associated with convergent tectonic boundaries, like subduction zones, which are repeatedly breached by giant megathrust earthquakes. The coastal relief along tectonically active realms is shaped by the effect of sea-level variations and heterogeneous patterns of permanent tectonic deformation, which are accumulated through several cycles of megathrust earthquakes. However, the correlation between earthquake deformation patterns and the sustained long-term segmentation of forearcs, particularly in Chile, remains poorly understood. Furthermore, the methods used to estimate permanent deformation from geomorphic markers, like marine terraces, have remained qualitative and are based on unrepeatable methods. This contrasts with the increasing resolution of digital elevation models, such as Light Detection and Ranging (LiDAR) and high-resolution bathymetric surveys.
Throughout this thesis I study permanent deformation in a holistic manner: from the methods to assess deformation rates, to the processes involved in its accumulation. My research focuses particularly on two aspects: Developing methodologies to assess permanent deformation using marine terraces, and comparing permanent deformation with seismic cycle deformation patterns under different spatial scales along the M8.8 Maule earthquake (2010) rupture zone. Two methods are developed to determine deformation rates from wave-built and wave-cut terraces respectively. I selected an archetypal example of a wave-built terrace at Santa Maria Island studying its stratigraphy and recognizing sequences of reoccupation events tied with eleven radiocarbon sample ages (14C ages). I developed a method to link patterns of reoccupation with sea-level proxies by iterating relative sea level curves for a range of uplift rates. I find the best fit between relative sea-level and the stratigraphic patterns for an uplift rate of 1.5 +- 0.3 m/ka.
A Graphical User Interface named TerraceM® was developed in Matlab®. This novel software tool determines shoreline angles in wave-cut terraces under different geomorphic scenarios. To validate the methods, I select test sites in areas of available high-resolution LiDAR topography along the Maule earthquake rupture zone and in California, USA. The software allows determining the 3D location of the shoreline angle, which is a proxy for the estimation of permanent deformation rates. The method is based on linear interpolations to define the paleo platform and cliff on swath profiles. The shoreline angle is then located by intersecting these interpolations. The
accuracy and precision of TerraceM® was tested by comparing its results with previous assessments, and through an experiment with students in a computer lab setting at the University
of Potsdam.
I combined the methods developed to analyze wave-built and wave-cut terraces to assess regional patterns of permanent deformation along the (2010) Maule earthquake rupture. Wave-built terraces are tied using 12 Infra Red Stimulated luminescence ages (IRSL ages) and shoreline angles in wave-cut terraces are estimated from 170 aligned swath profiles. The comparison of coseismic slip, interseismic coupling, and permanent deformation, leads to three areas of high permanent uplift, terrace warping, and sharp fault offsets. These three areas correlate with regions of high slip and low coupling, as well as with the spatial limit of at least eight historical megathrust ruptures (M8-9.5). I propose that the zones of upwarping at Arauco and Topocalma reflect changes in frictional properties of the megathrust, which result in discrete boundaries for the propagation of mega earthquakes.
To explore the application of geomorphic markers and quantitative morphology in offshore areas I performed a local study of patterns of permanent deformation inferred from hitherto unrecognized drowned shorelines at the Arauco Bay, at the southern part of the (2010) Maule earthquake rupture zone. A multidisciplinary approach, including morphometry, sedimentology, paleontology, 3D morphoscopy, and a landscape Evolution Model is used to recognize, map, and assess local rates and patterns of permanent deformation in submarine environments. Permanent deformation patterns are then reproduced using elastic models to assess deformation rates of an active submarine splay fault defined as Santa Maria Fault System. The best fit suggests a reverse structure with a slip rate of 3.7 m/ka for the last 30 ka. The register of land level changes during the earthquake cycle at Santa Maria Island suggest that most of the deformation may be accrued through splay fault reactivation during mega earthquakes, like the (2010) Maule event. Considering a recurrence time of 150 to 200 years, as determined from historical and geological observations, slip between 0.3 and 0.7 m per event would be required to account for the 3.7 m/ka millennial slip rate. However, if the SMFS slips only every ~1000 years, representing a few megathrust earthquakes, then a slip of ~3.5 m per event would be required to account for the long- term rate. Such event would be equivalent to a magnitude ~6.7 earthquake capable to generate a local tsunami.
The results of this thesis provide novel and fundamental information regarding the amount of permanent deformation accrued in the crust, and the mechanisms responsible for this accumulation at millennial time-scales along the M8.8 Maule earthquake (2010) rupture zone. Furthermore, the results of this thesis highlight the application of quantitative geomorphology and the use of repeatable methods to determine permanent deformation, improve the accuracy of marine terrace assessments, and estimates of vertical deformation rates in tectonically active coastal areas. This is vital information for adequate coastal-hazard assessments and to anticipate realistic earthquake and tsunami scenarios.
The repeated occurrence of exceptional floods within a few years, such as the Rhine floods in 1993 and 1995 and the Elbe and Danube floods in 2002 and 2013, suggests that floods in Central Europe may be organized in flood-rich and flood-poor periods. This hypothesis is studied by testing the significance of temporal clustering in flood occurrence (peak-over-threshold) time series for 68 catchments across Germany for the period 1932-2005. To assess the robustness of the results, different methods are used: Firstly, the index of dispersion, which quantifies the departure from a homogeneous Poisson process, is investigated. Further, the time-variation of the flood occurrence rate is derived by non-parametric kernel implementation and the significance of clustering is evaluated via parametric and non-parametric tests. Although the methods give consistent overall results, the specific results differ considerably. Hence, we recommend applying different methods when investigating flood clustering. For flood estimation and risk management, it is of relevance to understand whether clustering changes with flood severity and time scale. To this end, clustering is assessed for different thresholds and time scales. It is found that the majority of catchments show temporal clustering at the 5% significance level for low thresholds and time scales of one to a few years. However, clustering decreases substantially with increasing threshold and time scale. We hypothesize that flood clustering in Germany is mainly caused by catchment memory effects along with intra- to inter-annual climate variability, and that decadal climate variability plays a minor role. (C) 2016 Elsevier B.V. All rights reserved.
To understand past flood changes in the Rhine catchment and in particular the role of anthropogenic climate change in extreme flows, an attribution study relying on a proper GCM (general circulation model) downscaling is needed. A downscaling based on conditioning a stochastic weather generator on weather patterns is a promising approach. This approach assumes a strong link between weather patterns and local climate, and sufficient GCM skill in reproducing weather pattern climatology. These presuppositions are unprecedentedly evaluated here using 111 years of daily climate data from 490 stations in the Rhine basin and comprehensively testing the number of classification parameters and GCM weather pattern characteristics. A classification based on a combination of mean sea level pressure, temperature, and humidity from the ERA20C reanalysis of atmospheric fields over central Europe with 40 weather types was found to be the most appropriate for stratifying six local climate variables. The corresponding skill is quite diverse though, ranging from good for radiation to poor for precipitation. Especially for the latter it was apparent that pressure fields alone cannot sufficiently stratify local variability. To test the skill of the latest generation of GCMs from the CMIP5 ensemble in reproducing the frequency, seasonality, and persistence of the derived weather patterns, output from 15 GCMs is evaluated. Most GCMs are able to capture these characteristics well, but some models showed consistent deviations in all three evaluation criteria and should be excluded from further attribution analysis.
During expedition 202 aboard the RV Sonne in 2009, 39 seafloor surface sediment sites were sampled over a wide sector of the North Pacific and adjoining Bering Sea. The data served to infer land-ocean linkages of terrigenous sediment supply in terms of major sources and modes of sediment transport within an over-regional context. This is based on an integrated approach dealing with grain-size analysis, bulk mineralogy and clay mineralogy in combination with statistical data evaluation (end-member modelling of grain-size data, fuzzy cluster analysis of mineralogical data). The findings on clay mineralogy served to update those of earlier work extracted from the literature. Today, two processes of terrigenous sediment supply prevail in the study area: far-distance aeolian sediment supply to the pelagic North Pacific, and hemipelagic sediment dispersal from nearby land sources via ocean currents along the continental margins and island arcs. Aeolian particles show the finest grain sizes (clay and fine silt), whereas hemipelagic sediments have high abundances of coarse silt. Exposed sites on seamounts and the continental slope are partly swept by strong currents, leading to residual enrichment of fine sand. Four sediment sources can be distinguished on the basis of distinct index minerals revealed by statistical data analysis: dust plumes from central Asia (quartz, illite), altered materials from the volcanic regions of Kamchatka and the Aleutian Arc (smectite), detritus from the Alaskan Cordillera (chlorite, hornblende), and fluvial detritus from far-eastern Siberia and the Alaska mainland (quartz, feldspar, illite). These findings confirm those of former studies but considerably expand the geographic range of this suite of proxies as far south as 39A degrees N in the open North Pacific. The present integrated methodological approach proved useful in identifying the major modern processes of terrigenous sediment supply to the study region. This aspect deserves attention in the selection of sediment core sites for future palaeoenvironmental reconstructions related to aeolian and glacial dynamics, as well as the recognition of palaeo-ocean circulation patterns in general.
In the past, floods were basically managed by flood control mechanisms. The focus was set on the reduction of flood hazard. The potential consequences were of minor interest. Nowadays river flooding is increasingly seen from the risk perspective, including possible consequences. Moreover, the large-scale picture of flood risk became increasingly important for disaster management planning, national risk developments and the (re-) insurance industry. Therefore, it is widely accepted that risk-orientated flood management ap-proaches at the basin-scale are needed. However, large-scale flood risk assessment methods for areas of several 10,000 km² are still in early stages. Traditional flood risk assessments are performed reach wise, assuming constant probabilities for the entire reach or basin. This might be helpful on a local basis, but where large-scale patterns are important this approach is of limited use. Assuming a T-year flood (e.g. 100 years) for the entire river network is unrealistic and would lead to an overestimation of flood risk at the large scale. Due to the lack of damage data, additionally, the probability of peak discharge or rainfall is usually used as proxy for damage probability to derive flood risk. With a continuous and long term simulation of the entire flood risk chain, the spatial variability of probabilities could be consider and flood risk could be directly derived from damage data in a consistent way.
The objective of this study is the development and application of a full flood risk chain, appropriate for the large scale and based on long term and continuous simulation. The novel approach of ‘derived flood risk based on continuous simulations’ is introduced, where the synthetic discharge time series is used as input into flood impact models and flood risk is directly derived from the resulting synthetic damage time series.
The bottleneck at this scale is the hydrodynamic simu-lation. To find suitable hydrodynamic approaches for the large-scale a benchmark study with simplified 2D hydrodynamic models was performed. A raster-based approach with inertia formulation and a relatively high resolution of 100 m in combination with a fast 1D channel routing model was chosen.
To investigate the suitability of the continuous simulation of a full flood risk chain for the large scale, all model parts were integrated into a new framework, the Regional Flood Model (RFM). RFM consists of the hydrological model SWIM, a 1D hydrodynamic river network model, a 2D raster based inundation model and the flood loss model FELMOps+r. Subsequently, the model chain was applied to the Elbe catchment, one of the largest catchments in Germany. For the proof-of-concept, a continuous simulation was per-formed for the period of 1990-2003. Results were evaluated / validated as far as possible with available observed data in this period. Although each model part introduced its own uncertainties, results and runtime were generally found to be adequate for the purpose of continuous simulation at the large catchment scale.
Finally, RFM was applied to a meso-scale catchment in the east of Germany to firstly perform a flood risk assessment with the novel approach of ‘derived flood risk assessment based on continuous simulations’. Therefore, RFM was driven by long term synthetic meteorological input data generated by a weather generator. Thereby, a virtual time series of climate data of 100 x 100 years was generated and served as input to RFM providing subsequent 100 x 100 years of spatially consistent river discharge series, inundation patterns and damage values. On this basis, flood risk curves and expected annual damage could be derived directly from damage data, providing a large-scale picture of flood risk. In contrast to traditional flood risk analysis, where homogenous return periods are assumed for the entire basin, the presented approach provides a coherent large-scale picture of flood risk. The spatial variability of occurrence probability is respected. Additionally, data and methods are consistent. Catchment and floodplain processes are repre-sented in a holistic way. Antecedent catchment conditions are implicitly taken into account, as well as physical processes like storage effects, flood attenuation or channel–floodplain interactions and related damage influencing effects. Finally, the simulation of a virtual period of 100 x 100 years and consequently large data set on flood loss events enabled the calculation of flood risk directly from damage distributions. Problems associated with the transfer of probabilities in rainfall or peak runoff to probabilities in damage, as often used in traditional approaches, are bypassed.
RFM and the ‘derived flood risk approach based on continuous simulations’ has the potential to provide flood risk statements for national planning, re-insurance aspects or other questions where spatially consistent, large-scale assessments are required.
Flood loss modeling is an important component for risk analyses and decision support in flood risk management. Commonly, flood loss models describe complex damaging processes by simple, deterministic approaches like depth-damage functions and are associated with large uncertainty. To improve flood loss estimation and to provide quantitative information about the uncertainty associated with loss modeling, a probabilistic, multivariable Bagging decision Tree Flood Loss Estimation MOdel (BT-FLEMO) for residential buildings was developed. The application of BT-FLEMO provides a probability distribution of estimated losses to residential buildings per municipality. BT-FLEMO was applied and validated at the mesoscale in 19 municipalities that were affected during the 2002 flood by the River Mulde in Saxony, Germany. Validation was undertaken on the one hand via a comparison with six deterministic loss models, including both depth-damage functions and multivariable models. On the other hand, the results were compared with official loss data. BT-FLEMO outperforms deterministic, univariable, and multivariable models with regard to model accuracy, although the prediction uncertainty remains high. An important advantage of BT-FLEMO is the quantification of prediction uncertainty. The probability distribution of loss estimates by BT-FLEMO well represents the variation range of loss estimates of the other models in the case study.
Flood generation at the scale of large river basins is triggered by the interaction of the hydrological pre-conditions and the meteorological event conditions at different spatial and temporal scales. This interaction controls diverse flood generating processes and results in floods varying in magnitude and extent, duration as well as socio-economic consequences. For a process-based understanding of the underlying cause-effect relationships, systematic approaches are required. These approaches have to cover the complete causal flood chain, including the flood triggering meteorological event in combination with the hydrological (pre-)conditions in the catchment, runoff generation, flood routing, possible floodplain inundation and finally flood losses.
In this thesis, a comprehensive probabilistic process-based understanding of the causes and effects of floods is advanced. The spatial and temporal dynamics of flood events as well as the geophysical processes involved in the causal flood chain are revealed and the systematic interconnections within the flood chain are deciphered by means of the classification of their associated causes and effects. This is achieved by investigating the role of the hydrological pre-conditions and the meteorological event conditions with respect to flood occurrence, flood processes and flood characteristics as well as their interconnections at the river basin scale.
Broadening the knowledge about flood triggers, which up to now has been limited to linking large-scale meteorological conditions to flood occurrence, the influence of large-scale pre-event hydrological conditions on flood initiation is investigated. Using the Elbe River basin as an example, a classification of soil moisture, a key variable of pre-event conditions, is developed and a probabilistic link between patterns of soil moisture and flood occurrence is established. The soil moisture classification is applied to continuously simulated soil moisture data which is generated using the semi-distributed conceptual rainfall-runoff model SWIM. Applying successively a principal component analysis and a cluster analysis, days of similar soil moisture patterns are identified in the period November 1951 to October 2003.
The investigation of flood triggers is complemented by including meteorological conditions described by a common weather pattern classification that represents the main modes of atmospheric state variability. The newly developed soil moisture classification thereby provides the basis to study the combined impact of hydrological pre-conditions and large-scale meteorological event conditions on flood occurrence at the river basin scale.
A process-based understanding of flood generation and its associated probabilities is attained by classifying observed flood events into process-based flood types such as snowmelt floods or long-rain floods. Subsequently, the flood types are linked to the soil moisture and weather patterns. Further understanding of the processes is gained by modeling of the complete causal flood chain, incorporating a rainfall-runoff model, a 1D/2D hydrodynamic model and a flood loss model. A reshuffling approach based on weather patterns and the month of their occurrence is developed to generate synthetic data fields of meteorological conditions, which drive the model chain, in order to increase the flood sample size. From the large number of simulated flood events, the impact of hydro-meteorological conditions on various flood characteristics is detected through the analysis of conditional cumulative distribution functions and regression trees.
The results show the existence of catchment-scale soil moisture patterns, which comprise of large-scale seasonal wetting and drying components as well as of smaller-scale variations related to spatially heterogeneous catchment processes. Soil moisture patterns frequently occurring before the onset of floods are identified. In winter, floods are initiated by catchment-wide high soil moisture, whereas in summer the flood-initiating soil moisture patterns are diverse and the soil moisture conditions are less stable in time. The combined study of both soil moisture and weather patterns shows that the flood favoring hydro-meteorological patterns as well as their interactions vary seasonally. In the analysis period, 18 % of the weather patterns only result in a flood in the case of preceding soil saturation. The classification of 82 past events into flood types reveals seasonally varying flood processes that can be linked to hydro-meteorological patterns. For instance, the highest flood potential for long-rain floods is associated with a weather pattern that is often detected in the presence of so-called ‘Vb’ cyclones. Rain-on-snow and snowmelt floods are associated with westerly and north-westerly wind directions. The flood characteristics vary among the flood types and can be reproduced by the applied model chain. In total, 5970 events are simulated. They reproduce the observed event characteristics between September 1957 and August 2002 and provide information on flood losses. A regression tree analysis relates the flood processes of the simulated events to the hydro-meteorological (pre-)event conditions and highlights the fact that flood magnitude is primarily controlled by the meteorological event, whereas flood extent is primarily controlled by the soil moisture conditions.
Describing flood occurrence, processes and characteristics as a function of hydro-meteorological patterns, this thesis is part of a paradigm shift towards a process-based understanding of floods. The results highlight that soil moisture patterns as well as weather patterns are not only beneficial to a probabilistic conception of flood initiation but also provide information on the involved flood processes and the resulting flood characteristics.
Variational methods are employed in situations where exact Bayesian inference becomes intractable due to the difficulty in performing certain integrals. Typically, variational methods postulate a tractable posterior and formulate a lower bound on the desired integral to be approximated, e.g. marginal likelihood. The lower bound is then optimised with respect to its free parameters, the so-called variational parameters. However, this is not always possible as for certain integrals it is very challenging (or tedious) to come up with a suitable lower bound. Here, we propose a simple scheme that overcomes some of the awkward cases where the usual variational treatment becomes difficult. The scheme relies on a rewriting of the lower bound on the model log-likelihood. We demonstrate the proposed scheme on a number of synthetic and real examples, as well as on a real geophysical model for which the standard variational approaches are inapplicable.
We present a calibration-set based on modern pollen and satellite-based Advanced Very High Resolution Radiometer (AVHRR) observations of woody cover (including needleleaved, broadleaved and total tree cover) in eastern continental Asia, which shows good performance under cross-validation with the modern analogue technique (all the coefficients of determination between observed and predicted values are greater than 0.65). The calibration-set is used to reconstruct woody cover from a taxonomically harmonized and temporally standardized fossil pollen dataset (including 274 cores) with 500-year resolution over the last 22 kyr. The spatial range of forest has not noticeably changed in eastern continental Asia during the last 22 kyr, although woody cover has, especially at the margin of the eastern Tibetan Plateau and in the forest-steppe transition area of north-central China. Vegetation was sparse during the LGM in the present forested regions, but woody cover increased markedly at the beginning of the Bolling/Allerod period (B/A; ca. 14.5 ka BP) and again at the beginning of the Holocene (ca. 11.5 ka BP), and is related to the enhanced strength of the East Asian Summer Monsoon. Forest flourished in the mid Holocene (ca. 8 ka BP) possibly due to favourable climatic conditions. In contrast, cover was stable in southern China (high cover) and arid central Asia (very low cover) throughout the investigated period. Forest cover increased in the north-eastern part of China during the Holocene. Comparisons of these regional pollen-based results with simulated forest cover from runs of a global climate model (for 9, 6 and 0 ka BP (ECHAM5/JSBACH similar to 1.125 degrees spatial resolution)) reveal many similarities in temporal change. The Holocene woody cover history of eastern continental Asia is different from that of other regions, likely controlled by different climatic variables, i.e. moisture in eastern continental Asia; temperature in northern Eurasia and North America. (C) 2016 Elsevier Ltd. All rights reserved.
Late Quaternary paleoenvironmental records from the Chatanika River valley near Fairbanks (Alaska)
(2016)
Perennially-frozen deposits are considered as excellent paleoenvironmental archives similar to lacustrine, deep marine, and glacier records because of the long-term and good preservation of fossil records under stable permafrost conditions. A permafrost tunnel in the Vault Creek Valley (Chatanika River Valley, near Fairbanks) exposes a sequence of frozen deposits and ground ice that provides a comprehensive set of proxies to reconstruct the late Quaternary environmental history of Interior Alaska. The multi-proxy approach includes different dating techniques (radiocarbon-accelerator mass spectrometry [AMS C-14], optically stimulated luminescence [OSL], thorium/uranium radioisotope disequilibria [Th-230/U]), as well as methods of sedimentology, paleoecology, hydrochemistry, and stable isotope geochemistry of ground ice. The studied sequence consists of 36-m-thick late Quaternary deposits above schistose bedrock. Main portions of the sequence accumulated during the early and middle Wisconsin periods. The lowermost unit A consists of about 9-m-thick ice-bonded fluvial gravels with sand and peat lenses. A late Sangamon (MIS 5a) age of unit A is assumed. Spruce forest with birch, larch, and some shrubby alder dominated the vegetation. High presence of Sphagnum spores and Cyperaceae pollen points to mires in the Vault Creek Valley. The overlying unit B consists of 10-m-thick alternating fluvial gravels, loess-like silt, and sand layers, penetrated by small ice wedges. OSL dates support a stadial early Wisconsin (MIS 4) age of unit B. Pollen and plant macrofossil data point to spruce forests with some birch interspersed with wetlands around the site. The following unit C is composed of 15-m-thick ice-rich loess-like and organic-rich silt with fossil bones and large ice wedges. Unit C formed during the interstadial mid-Wisconsin (MIS 3) and stadial late Wisconsin (MIS 2) as indicated by radiocarbon ages. Post-depositional slope processes significantly deformed both, ground ice and sediments of unit C. Pollen data show that spruce forests and wetlands dominated the area. The macrofossil remains of Picea, Larix, and Alnus incana ssp. tenuifolia also prove the existence of boreal coniferous forests during the mid-Wisconsin interstadial, which were replaced by treeless tundra-steppe vegetation during the late Wisconsin stadial. Unit C is discordantly overlain by the 2-m-thick late Holocene deposits of unit D. The pollen record of unit D indicates boreal forest vegetation similar to the modern one. The permafrost record from the Vault Creek tunnel reflects more than 90 ka of periglacial landscape dynamics triggered by fluvial and eolian accumulation, and formation of ice-wedge polygons and post depositional deformation by slope processes. The record represents a typical Wisconsin valley-bottom facies in Central Alaska. (C) 2016 Elsevier Ltd. All rights reserved.
Waterbodies such as lakes and ponds are abundant in vast Arctic landscapes and strongly affect the thermal state of the surrounding permafrost. In order to gain a better understanding of the impact of small-and medium-sized waterbodies on permafrost and the formation of thermokarst, a land surface model was developed that can represent the vertical and lateral thermal interactions between waterbodies and permafrost. The model was validated using temperature measurements from two typical waterbodies located within the Lena River delta in northern Siberia. Impact simulations were performed under current climate conditions as well as under a moderate and a strong climate-warming scenario. The performed simulations demonstrate that small waterbodies can rise the sediment surface temperature by more than 10 degrees C and accelerate permafrost thaw by a factor of between 4 and 5. Up to 70% of this additional heat flux into the ground was found to be dissipated into the surrounding permafrost by lateral ground heat flux in the case of small, shallow, and isolated waterbodies. Under moderate climate warming, the lateral heat flux was found to reduce permafrost degradation underneath waterbodies by a factor of 2. Under stronger climatic warming, however, the lateral heat flux was too small to prevent rapid permafrost degradation. The lateral heat flux was also found to strongly impede the formation of thermokarst. Despite this stabilizing effect, our simulations have demonstrated that underneath shallow waterbodies (<1 m), thermokarst initiation happens 30 to 40 years earlier than in simulations without preexisting waterbody.
We report the first high-resolution (20-50 years) mid- to late Holocene pollen records from Lake Teletskoye, the largest lake in the Altai Mountains, in south-eastern West Siberia. Generally, the mid- to late Holocene (the last 4250 years) vegetation of the north-eastern Altai, as recorded in two studied sediment cores, is characterised by Siberian pine-spruce-fir forests that are similar to those of the present day. A relatively cool and dry interval with July temperatures lower than those of today occurred between 3.9 and 3.6 ka BP. The widespread distribution of open, steppe-like communities with Artemisia, Chenopodiaceae and Cyperaceae reflects maximum deforestation during this interval. After ca. 3.5 ka BP, the coniferous mountain taiga spread significantly, with maximum woody coverage and taiga biome scores between ca. 2.7 and 1.6 ka BP. This coincides well with the highest July temperature (approximately 1 degrees C higher than today) intervals. A short period of cooling about 13-1.4 ka BP could have been triggered by the increased volcanic activity recorded across the Northern Hemisphere. A new period of cooling started around 1100-1150 CE, with the minimum July temperatures occurring between 1450 and 1800 CE. (C) 2016 Elsevier B.V. All rights reserved.
Freshwater ostracods (Crustacea, Ostracoda) are valuable biological indicators. In Arctic environments, their habitat conditions are barely known and the abundance and diversity of ostracods is documented only in scattered records with incomplete ecological characterization. To determine the taxonomic range of ostracod assemblages and their habitat conditions in polygon ponds in the Indigirka Lowland, north-east Siberia, we collected more than 100 living ostracod individuals per site with a plankton net (mesh size 65 mm) and an exhaustor system from 27 water bodies and studied them in the context of substrate and hydrochemical data. During the summer of 2011, a single pond site and its ostracod population was selected for special study. This first record of the ostracod fauna in the Indigirka Lowland comprises eight species and three additional taxa. Fabaeformiscandona krochini and F. groenlandica were documented for the first time in continental Siberia. Repeated sampling of a low-centre polygon pond yielded insights into the population dynamics of F. pedata. We identified air temperature and precipitation as the main external drivers of water temperatures, water levels, ion concentrations and water stable isotope composition on diurnal and seasonal scales.
Glacial legacies on interglacial vegetation at the Pliocene-Pleistocene transition in NE Asia
(2016)
Broad-scale climate control of vegetation is widely assumed. Vegetation-climate lags are generally thought to have lasted no more than a few centuries. Here our palaeoecological study challenges this concept over glacial–interglacial timescales. Through multivariate analyses of pollen assemblages from Lake El’gygytgyn, Russian Far East and other data we show that interglacial vegetation during the Plio-Pleistocene transition mainly reflects conditions of the preceding glacial instead of contemporary interglacial climate. Vegetation–climate disequilibrium may persist for several millennia, related to the combined effects of permafrost persistence, distant glacial refugia and fire. In contrast, no effects from the preceding interglacial on glacial vegetation are detected. We propose that disequilibrium was stronger during the Plio-Pleistocene transition than during the Mid-Pliocene Warm Period when, in addition to climate, herbivory was important. By analogy to the past, we suggest today’s widespread larch ecosystem on permafrost is not in climate equilibrium. Vegetation-based reconstructions of interglacial climates used to assess atmospheric CO2–temperature relationships may thus yield misleading simulations of past global climate sensitivity.
QuestionHow important is the effect of micro-relief and vegetation type on the characteristics of vascular plants and bryophytes in low-centred polygons? LocationSiberian Arctic, Russia. MethodsEight low-centred polygons in northern Siberia were surveyed for vegetation along transects running from the rim to the pond via the rim-pond transition of each polygon and across a vegetation type gradient from open forest to tundra. ResultsThe cover of vascular plants and bryophytes displays no significant differences between the rim and rim-pond transition but is significantly lower in the pond section of the polygons. Alpha-diversity of vascular plants decreases strongly from rim to pond, whereas bryophyte diversity in pond plots is significantly distinct from the rim and the rim-pond transition. There is no clear trend in cover for either plant group along the vegetation type transect and only a weak trend in -diversity. However, both gradients are reflected in the compositional turnover. The applied indicator species analysis identified taxa characteristic of certain environmental conditions. Among others, we found vascular plants primarily characteristic of the rim and bryophyte taxa characteristic of each micro-relief level and vegetation type. ConclusionsThe observed gradual pattern in -diversity and composition of polygonal vegetation suggests that micro-relief is the main driver of changes in the vegetation composition, while vegetation type and the related forest cover change are of subordinate importance for polygonal vegetation patterns along the Siberian tree line.
Abrupt climate changes and fluctuations over short time scales are superimposed on long-term climate changes. Understanding rapid climate fluctuations at the decadal time scale over the past millennium will enhance our understanding of patterns of climate variability and aid in forecasting climate changes in the future. In this study, climate changes on the southeastern Tibetan Plateau over the past millennium were determined from a 4.82-m-long sediment core from Basomtso Lake. At the centennial time scale, the Medieval Climate Anomaly (MCA), Little Ice Age (LIA) and Current Warm Period (CWP) are distinct in the Basomtso region. Rapid climate fluctuations inferred from five episodes with higher sediment input and likely warmer conditions, as well as seven episodes with lower sediment input and likely colder conditions, were well preserved in our record. These episodes with higher and lower sediment input are characterized by abrupt climate changes and short time durations. Spectral analysis indicates that the climate variations at the centennial scale on the southeastern Tibetan Plateau are influenced by solar activity during the past millennium.
We compare the ability of various site-condition proxies (SCPs) to reduce the aleatory variability of ground motion prediction equations (GMPEs). Three SCPs (measured V-S30, inferred V-S30, local topographic slope) and two accelerometric databases (RESORCE and NGA-West2) are considered. An artificial neural network (ANN) approach including a random-effect procedure is used to derive GMPEs setting the relationship between peak ground acceleration (PGA), peak ground velocity (PGV), pseudo-spectral acceleration [PSA(T)], and explanatory variables (M-w, R-JB, and V-S30 or Slope). The analysis is performed using both discrete site classes and continuous proxy values. All "non-measured" SCPs exhibit a rather poor performance in reducing aleatory variability, compared to the better performance of measured V-S30. A new, fully data-driven GMPE based on the NGA-West2 is then derived, with an aleatory variability value depending on the quality of the SCP. It proves very consistent with previous GMPEs built on the same data set. Measuring V-S30 allows for benefit from an aleatory variability reduction up to 15%.
A new view of Ecuador's complex geodynamics has been developed in the course of modeling seismic source zones for probabilistic seismic hazard analysis. This study focuses on two aspects of the plates' interaction at a continental scale: (a) age-related differences in rheology between Farallon and Nazca plates—marked by the Grijalva rifted margin and its inland projection—as they subduct underneath central Ecuador, and (b) the rapidly changing convergence obliquity resulting from the convex shape of the South American northwestern continental margin. Both conditions satisfactorily explain several characteristics of the observed seismicity and of the interseismic coupling. Intermediate-depth seismicity reveals a severe flexure in the Farallon slab as it dips and contorts at depth, originating the El Puyo seismic cluster. The two slabs position and geometry below continental Ecuador also correlate with surface expressions observable in the local and regional geology and tectonics. The interseismic coupling is weak and shallow south of the Grijalva rifted margin and increases northward, with a heterogeneous pattern locally associated to the Carnegie ridge subduction. High convergence obliquity is responsible for the North Andean Block northeastward movement along localized fault systems. The Cosanga and Pallatanga fault segments of the North Andean Block-South American boundary concentrate most of the seismic moment release in continental Ecuador. Other inner block faults located along the western border of the inter-Andean Depression also show a high rate of moderate-size earthquake production. Finally, a total of 19 seismic source zones were modeled in accordance with the proposed geodynamic and neotectonic scheme.
Event coincidence analysis for quantifying statistical interrelationships between event time series
(2016)
Studying event time series is a powerful approach for analyzing the dynamics of complex dynamical systems in many fields of science. In this paper, we describe the method of event coincidence analysis to provide a framework for quantifying the strength, directionality and time lag of statistical interrelationships between event series. Event coincidence analysis allows to formulate and test null hypotheses on the origin of the observed interrelationships including tests based on Poisson processes or, more generally, stochastic point processes with a prescribed inter-event time distribution and other higher-order properties. Applying the framework to country-level observational data yields evidence that flood events have acted as triggers of epidemic outbreaks globally since the 1950s. Facing projected future changes in the statistics of climatic extreme events, statistical techniques such as event coincidence analysis will be relevant for investigating the impacts of anthropogenic climate change on human societies and ecosystems worldwide.
On 16 September 2015, the M-W = 8.2 Illapel megathrust earthquake ruptured the Central Chilean margin. Combining inversions of displacement measurements and seismic waveforms with high frequency (HF) teleseismic backprojection, we derive a comprehensive description of the rupture, which also predicts deep ocean tsunami wave heights. We further determine moment tensors and obtain accurate depth estimates for the aftershock sequence. The earthquake nucleated near the coast but then propagated to the north and updip, attaining a peak slip of 5-6 m. In contrast, HF seismic radiation is mostly emitted downdip of the region of intense slip and arrests earlier than the long period rupture, indicating smooth slip along the shallow plate interface in the final phase. A superficially similar earthquake in 1943 with a similar aftershock zone had a much shorter source time function, which matches the duration of HF seismic radiation in the recent event, indicating that the 1943 event lacked the shallow slip.
High-resolution topographic data greatly facilitate the remote identification of geomorphic features, furnishing valuable information concerning surface processes and characterization of reference markers for quantifying tectonic deformation. Marine terraces have been used as long baseline geodetic markers of relative past sea-level positions, reflecting the interplay between vertical crustal movements and sea-level oscillations. Uplift rates may be determined from the terrace age and the elevation of its shoreline angle, a geomorphic feature that can be correlated with past sea-levels positions. A precise definition of the shoreline angle in time and space is essential to obtain reliable uplift rates with coherent spatial correlation. To improve our ability to rapidly assess and map shoreline angles at regional and local scales, we have developed TerraceM, a MATLAB (R) graphical user interface that allows the shoreline angle and its associated error to be estimated using high-resolution topography. TerraceM uses topographic swath profiles oriented orthogonally to the terrace riser. Four functions are included to analyze the swath profiles and extract the shoreline angle, from both staircase sequences of multiple terraces and rough coasts characterized by eroded remnants of emerged terrace surfaces. The former are measured by outlining the paleocliffs and paieo-platforms and finding their intersection by extrapolating linear regressions, whereas the latter are assessed by automatically detecting peaks of sea-stack tops and back-projecting them to the modern sea cliff. In the absence of rigorous absolute age determinations of marine terraces, their geomorphic age may be estimated using previously published diffusion models. Postprocessing functions are included to obtain first-order statistics of shoreline-angle elevations and their spatial distribution. TerraceM has the ability to process series of profiles from several sites in an efficient and structured workflow. Results may be exported in Google Earth and ESRI shapefile formats. The precision and accuracy of the method have been estimated from a case study at Santa Cruz, California, by comparing TerraceM results with published field measurements. The repeatability was evaluated using multiple measurements made by inexperienced users. TerraceM will improve the efficiency and precision of estimating shoreline-angle elevations in wave-cut terraces in both marine and lacustrine environments.
Although the climate development over the Holocene in the Northern Hemisphere is well known, palaeolimnological climate reconstructions reveal spatiotemporal variability in northern Eurasia. Here we present a multi-proxy study from north-eastern Siberia combining sediment geochemistry, and diatom and pollen data from lake-sediment cores covering the last 38,000 cal. years. Our results show major changes in pyrite content and fragilarioid diatom species distributions, indicating prolonged seasonal lake-ice cover between similar to 13,500 and similar to 8900 cal. years BP and possibly during the 8200 cal. years BP cold event. A pollen-based climate reconstruction generated a mean July temperature of 17.8 degrees C during the Holocene Thermal Maximum (HTM) between similar to 8900 and similar to 4500 cal. years BP. Naviculoid diatoms appear in the late Holocene indicating a shortening of the seasonal ice cover that continues today. Our results reveal a strong correlation between the applied terrestrial and aquatic indicators and natural seasonal climate dynamics in the Holocene. Planktonic diatoms show a strong response to changes in the lake ecosystem due to recent climate warming in the Anthropocene. We assess other palaeolimnological studies to infer the spatiotemporal pattern of the HTM and affirm that the timing of its onset, a difference of up to 3000 years from north to south, can be well explained by climatic teleconnections. The westerlies brought cold air to this part of Siberia until the Laurentide ice sheet vanished 7000 years ago. The apparent delayed ending of the HTM in the central Siberian record can be ascribed to the exceedance of ecological thresholds trailing behind increases in winter temperatures and decreases in contrast in insolation between seasons during the mid to late Holocene as well as lacking differentiation between summer and winter trends in paleolimnological reconstructions. (C) 2015 Elsevier Ltd. All rights reserved.