Refine
Year of publication
Document Type
- Article (2695)
- Doctoral Thesis (504)
- Postprint (139)
- Other (72)
- Review (51)
- Monograph/Edited Volume (34)
- Preprint (17)
- Conference Proceeding (12)
- Habilitation Thesis (12)
- Master's Thesis (6)
Keywords
- climate change (49)
- Holocene (44)
- erosion (28)
- Himalaya (26)
- remote sensing (25)
- permafrost (23)
- Climate change (21)
- Tibetan Plateau (21)
- Pollen (20)
- climate (20)
Institute
- Institut für Geowissenschaften (3549) (remove)
Frequent and intense rainfall events demand innovative techniques to better predict the extreme rainfall dynamics. This task requires essentially the assessment of the basic types of atmospheric processes that trigger extreme rainfall, and then to examine the differences between those processes, which may help to identify key patterns to improve predictive algorithms. We employ tools from network theory to compare the spatial features of extreme rainfall over the Japanese archipelago and surrounding areas caused by two atmospheric processes: the Baiu front, which occurs mainly in June and July (JJ), and the tropical storms from August to November (ASON). We infer from complex networks of satellite-derived rainfall data, which are based on the nonlinear correlation measure of event synchronization. We compare the spatial scales involved in both systems and identify different regions which receive rainfall due to the large spatial scale of the Baiu and tropical storm systems. We observed that the spatial scales involved in the Baiu driven rainfall extremes, including the synoptic processes behind the frontal development, are larger than tropical storms, which even have long tracks during extratropical transitions. We further delineate regions of coherent rainfall during the two seasons based on network communities, identifying the horizontal (east-west) rainfall bands during JJ over the Japanese archipelago, while during ASON these bands align with the island arc of Japan.
The polar and subtropical jet streams are strong upper-level winds with a crucial influence on weather throughout the Northern Hemisphere midlatitudes. In particular, the polar jet is located between cold arctic air to the north and warmer subtropical air to the south. Strongly meandering states therefore often lead to extreme surface weather. The parameter values of the detection scheme are optimized using simulated annealing and a skill function that accounts for the zonal-mean jet stream position (Rikus, 2015). After the successful optimization process, we apply our scheme to reanalysis data covering 1979-2015 and calculate seasonal-mean probabilistic maps and trends in wind strength and position of jet streams. We present longitudinally defined probability distributions of the positions for both jets for all on the Northern Hemisphere seasons. This shows that winter is characterized by two well-separated jets over Europe and Asia (ca. 20 degrees W to 140 degrees E). In contrast, summer normally has a single merged jet over the western hemisphere but can have both merged and separated jet states in the eastern hemisphere.
Flow accumulation algorithms estimate the steady state of flow on real or modeled topographic surfaces and are crucial for hydrological and geomorphological assessments, including delineation of river networks, drainage basins, and sediment transport processes. Existing flow accumulation algorithms are typically designed to compute flows on regular grids and are not directly applicable to arbitrarily sampled topographic data such as lidar point clouds. In this study we present a random sampling scheme that generates homogeneous point densities, in combination with a novel flow path tracing approach-the Facet-Flow Network (FFN)-that estimates flow accumulation in terms of specific catchment area (SCA) on triangulated surfaces. The random sampling minimizes biases due to spatial sampling and the FFN allows for direct flow estimation from point clouds. We validate our approach on a Gaussian hill surface and study the convergence of its SCA compared to the analytical solution. Here, our algorithm outperforms the multiple flow direction algorithm, which is optimized for divergent surfaces. We also compute the SCA of a 6-km(2)-steep, vegetated catchment on Santa Cruz Island, California, based on airborne lidar point-cloud data. Point-cloud-based SCA values estimated by our method compare well with those estimated by the D-infinity or multiple flow direction algorithm on gridded data. The advantage of computing SCA from point clouds becomes relevant especially for divergent topography and for small drainage areas: These are depicted with much more detail due to the higher sampling density of point clouds.
The complementary advantages of GPS and seismic measurements are well recognized in seismotectonic monitoring studies. Therefore, integrated processing of the two data streams has been proposed recently in an attempt to obtain accurate and reliable information of surface displacements associated with earthquakes. A hitherto still critical issue in the integrated processing is real-time detection and precise estimation of the transient baseline error in the seismic records. Here, we report on a new approach by introducing the seismic acceleration corrected by baseline errors into the state equation system. The correction is performed and regularly updated in short epochs (with increments which may be as short as seconds), so that station position, velocity, and acceleration can be constrained very tightly and baseline error can be estimated as a random-walk process. With the adapted state equation system, our study highlights the use of a new approach developed for integrated processing of GPS and seismic data by means of sequential least-squares adjustment. The efficiency of our approach is demonstrated and validated using simulated, experimental, and real datasets. The latter were collected at collocated GPS and seismic stations around the 4 April 2010, E1 Mayor-Cucapah earthquake (Mw, 7.2). The results have shown that baseline errors of the strong-motion sensors are corrected precisely and high-precision seismic displacements are real-timely obtained by the new approach.
With increasing amount of strong motion data, Ground Motion Prediction Equation (GMPE) developers are able to quantify empirical site amplification functions (delta S2S(s)) from GMPE residuals, for use in site-specific Probabilistic Seismic Hazard Assessment. In this study, we first derive a GMPE for 5% damped Pseudo Spectral Acceleration (g) of Active Shallow Crustal earthquakes in Japan with 3.4 <= M-w <= 7.3 and 0 <= R-JB <= 600km. Using k-mean spectral clustering technique, we then classify our estimated delta S2S(s)(T = 0.01 - 2s) of 588 wellcharacterized sites, into 8 site clusters with distinct mean site amplification functions, and within-cluster site-tosite variability similar to 50% smaller than the overall dataset variability (phi(S2S)). Following an evaluation of existing schemes, we propose a revised data-driven site classification characterized by kernel density distributions of V-s30, V-s10, H-800, and predominant period (T-G) of the site clusters.
Quantitative estimates of sea-level rise in the Mediterranean Basin become increasingly accurate thanks to detailed satellite monitoring. However, such measuring campaigns cover several years to decades, while longer-term sea-level records are rare for the Mediterranean. We used a data archeological approach to reanalyze monthly mean sea-level data of the Antalya-I (1935–1977) tide gauge to fill this gap. We checked the accuracy and reliability of these data before merging them with the more recent records of the Antalya-II (1985–2009) tide gauge, accounting for an eight-year hiatus. We obtain a composite time series of monthly and annual mean sea levels spanning some 75 years, providing the longest record for the eastern Mediterranean Basin, and thus an essential tool for studying the region's recent sea-level trends. We estimate a relative mean sea-level rise of 2.2 ± 0.5 mm/year between 1935 and 2008, with an annual variability (expressed here as the standard deviation of the residuals, σresiduals = 41.4 mm) above that at the closest tide gauges (e.g., Thessaloniki, Greece, σresiduals = 29.0 mm). Relative sea-level rise accelerated to 6.0 ± 1.5 mm/year at Antalya-II; we attribute roughly half of this rate (~3.6 mm/year) to tectonic crustal motion and anthropogenic land subsidence. Our study highlights the value of data archeology for recovering and integrating historic tide gauge data for long-term sea-level and climate studies.
Sulfur is an important component in volcanic gases at the Earth surface but also present in the deep Earth in hydrothermal or magmatic fluids. Little is known about the evolution of such fluids during ascent in the crust. A new optical cell was developed for in situ Raman spectroscopic investigations on fluids allowing abrupt or continuous changes of pressure up to 200 MPa at temperatures up to 750 degrees C. The concept is based on a flexible gold bellow, which separates the sample fluid from the pressure medium water. To avoid reactions between aggressive fluids and the pressure cell, steel components in contact with the fluid are shielded by gold foil. The cell was tested to study redox reactions in fluids using aqueous ammonium sulfate solutions as a model system. During heating at constant pressure of 130 MPa, sulfate ions transform first to HSO4- ions and then to molecular units such as H2SO4. Variation of pressure shows that the stability of sulfate species relies on fluid density, i.e., highly charged species are stable only in high-density fluids. Partial decomposition of ammonium was evident above 550 degrees C by the occurrence of a nitrogen peak in the Raman spectra. Reduced sulfur species were observed above 700 degrees C by Raman signals near 2590 cm(-1) assigned to HS- and H2S. No clear evidence for the formation of sulfur dioxide was found in contrary to previous studies on aqueous H2SO4, suggesting very reducing conditions in our experiments. Fluid-mineral interaction was studied by inserting into the cell a small, semi-open capsule filled with a mixture of pyrite and pyrrhotite. Oxidation of the sample assembly was evident by transformation of pyrite to pyrrhotite. As a consequence, sulfide species were observed in the fluid already at temperatures of similar to 600 degrees C.
The age models of fluvio-lacustrine sedimentary sequences are often subject of discussions in paleoclimate research. The techniques employed to build an age model are very diverse, ranging from visual or intuitive estimation of the age-depth relationship over linear or spline interpolations between age control points to sophisticated Bayesian techniques also taking into account the most likely deposition times of the type of sediment within the sequence. All these methods, however, fail in detecting abrupt variations in sedimentation rates, including the possibility of episodes of no deposition (hiatus), which is the strength of the method presented in this work. The new technique simply compares the deposition time of equally thick sediment slices from the differences of subsequent radiometric age dates and the unit deposition times of the various sediment types. The percentage overlap of the distributions of these two sources of information, together with the evidence from the sedimentary record, helps to build an age model of complex sequences including abrupt variations in the rate of deposition including one or many hiatuses. (C) 2014 Elsevier B.V. All rights reserved.
A new view of Ecuador's complex geodynamics has been developed in the course of modeling seismic source zones for probabilistic seismic hazard analysis. This study focuses on two aspects of the plates' interaction at a continental scale: (a) age-related differences in rheology between Farallon and Nazca plates—marked by the Grijalva rifted margin and its inland projection—as they subduct underneath central Ecuador, and (b) the rapidly changing convergence obliquity resulting from the convex shape of the South American northwestern continental margin. Both conditions satisfactorily explain several characteristics of the observed seismicity and of the interseismic coupling. Intermediate-depth seismicity reveals a severe flexure in the Farallon slab as it dips and contorts at depth, originating the El Puyo seismic cluster. The two slabs position and geometry below continental Ecuador also correlate with surface expressions observable in the local and regional geology and tectonics. The interseismic coupling is weak and shallow south of the Grijalva rifted margin and increases northward, with a heterogeneous pattern locally associated to the Carnegie ridge subduction. High convergence obliquity is responsible for the North Andean Block northeastward movement along localized fault systems. The Cosanga and Pallatanga fault segments of the North Andean Block-South American boundary concentrate most of the seismic moment release in continental Ecuador. Other inner block faults located along the western border of the inter-Andean Depression also show a high rate of moderate-size earthquake production. Finally, a total of 19 seismic source zones were modeled in accordance with the proposed geodynamic and neotectonic scheme.
In the past, floods were basically managed by flood control mechanisms. The focus was set on the reduction of flood hazard. The potential consequences were of minor interest. Nowadays river flooding is increasingly seen from the risk perspective, including possible consequences. Moreover, the large-scale picture of flood risk became increasingly important for disaster management planning, national risk developments and the (re-) insurance industry. Therefore, it is widely accepted that risk-orientated flood management ap-proaches at the basin-scale are needed. However, large-scale flood risk assessment methods for areas of several 10,000 km² are still in early stages. Traditional flood risk assessments are performed reach wise, assuming constant probabilities for the entire reach or basin. This might be helpful on a local basis, but where large-scale patterns are important this approach is of limited use. Assuming a T-year flood (e.g. 100 years) for the entire river network is unrealistic and would lead to an overestimation of flood risk at the large scale. Due to the lack of damage data, additionally, the probability of peak discharge or rainfall is usually used as proxy for damage probability to derive flood risk. With a continuous and long term simulation of the entire flood risk chain, the spatial variability of probabilities could be consider and flood risk could be directly derived from damage data in a consistent way.
The objective of this study is the development and application of a full flood risk chain, appropriate for the large scale and based on long term and continuous simulation. The novel approach of ‘derived flood risk based on continuous simulations’ is introduced, where the synthetic discharge time series is used as input into flood impact models and flood risk is directly derived from the resulting synthetic damage time series.
The bottleneck at this scale is the hydrodynamic simu-lation. To find suitable hydrodynamic approaches for the large-scale a benchmark study with simplified 2D hydrodynamic models was performed. A raster-based approach with inertia formulation and a relatively high resolution of 100 m in combination with a fast 1D channel routing model was chosen.
To investigate the suitability of the continuous simulation of a full flood risk chain for the large scale, all model parts were integrated into a new framework, the Regional Flood Model (RFM). RFM consists of the hydrological model SWIM, a 1D hydrodynamic river network model, a 2D raster based inundation model and the flood loss model FELMOps+r. Subsequently, the model chain was applied to the Elbe catchment, one of the largest catchments in Germany. For the proof-of-concept, a continuous simulation was per-formed for the period of 1990-2003. Results were evaluated / validated as far as possible with available observed data in this period. Although each model part introduced its own uncertainties, results and runtime were generally found to be adequate for the purpose of continuous simulation at the large catchment scale.
Finally, RFM was applied to a meso-scale catchment in the east of Germany to firstly perform a flood risk assessment with the novel approach of ‘derived flood risk assessment based on continuous simulations’. Therefore, RFM was driven by long term synthetic meteorological input data generated by a weather generator. Thereby, a virtual time series of climate data of 100 x 100 years was generated and served as input to RFM providing subsequent 100 x 100 years of spatially consistent river discharge series, inundation patterns and damage values. On this basis, flood risk curves and expected annual damage could be derived directly from damage data, providing a large-scale picture of flood risk. In contrast to traditional flood risk analysis, where homogenous return periods are assumed for the entire basin, the presented approach provides a coherent large-scale picture of flood risk. The spatial variability of occurrence probability is respected. Additionally, data and methods are consistent. Catchment and floodplain processes are repre-sented in a holistic way. Antecedent catchment conditions are implicitly taken into account, as well as physical processes like storage effects, flood attenuation or channel–floodplain interactions and related damage influencing effects. Finally, the simulation of a virtual period of 100 x 100 years and consequently large data set on flood loss events enabled the calculation of flood risk directly from damage distributions. Problems associated with the transfer of probabilities in rainfall or peak runoff to probabilities in damage, as often used in traditional approaches, are bypassed.
RFM and the ‘derived flood risk approach based on continuous simulations’ has the potential to provide flood risk statements for national planning, re-insurance aspects or other questions where spatially consistent, large-scale assessments are required.
Nitraria is a halophytic taxon (i.e., adapted to saline environments) that belongs to the plant family Nitrariaceae and is distributed from the Mediterranean, across Asia into the south-eastern tip of Australia. This taxon is thought to have originated in Asia during the Paleogene (66-23 Ma), alongside the proto-Paratethys epicontinental sea. The evolutionary history of Nitraria might hold important clues on the links between climatic and biotic evolution but limited taxonomic documentation of this taxon has thus far hindered this line of research. Here we investigate if the pollen morphology and the chemical composition of the pollen wall are informative of the evolutionary history of Nitraria and could explain if origination along the proto-Paratethys and dispersal to the Tibetan Plateau was simultaneous or a secondary process. To answer these questions, we applied a novel approach consisting of a combination of Fourier Transform Infrared spectroscopy (FTIR), to determine the chemical composition of the pollen wall, and pollen morphological analyses using Light Microscopy (LM) and Scanning Electron Microscopy (SEM). We analysed our data using ordinations (principal components analysis and non-metric multidimensional scaling), and directly mapped it on the Nitrariaceae phylogeny to produce a phylomorphospace and a phylochemospace. Our LM, SEM and FTIR analyses show clear morphological and chemical differences between the sister groups Peganum and Nitraria. Differences in the morphological and chemical characteristics of highland species (Nitraria schoberi, N. sphaerocarpa, N. sibirica and N. tangutorum) and lowland species (Nitraria billardierei and N. retusa) are very subtle, with phylogenetic history appearing to be a more important control on Nitraria pollen than local environmental conditions. Our approach shows a compelling consistency between the chemical and morphological characteristics of the eight studied Nitrariaceae species, and these traits are in agreement with the phylogenetic tree. Taken together, this demonstrates how novel methods for studying fossil pollen can facilitate the evolutionary investigation of living and extinct taxa, and the environments they represent.
Stream restoration aims at an enhancement of ecological habitats, an increase of water retention within a landscape and sometimes even at an improvement of biogeochemical functions of lotic ecosystems. For the latter, good exchange between groundwater and stream water is often considered to be of major importance. In this study hydraulic connectivity between river and aquifer was investigated for a four years period, covering the restoration of an old oxbow after the second year. The oxbow became reconnected to the stream and the clogging layer in the oxbow was excavated. We expected increasing hydraulic connectivity between oxbow and aquifer after restoration of the stream, and decreasing hydraulic connectivity for the former shortcut due to increased clogging. To test that hypothesis, the spatial and temporal characteristics of the coupled groundwater-stream water system before and after the restoration were analysed by principal component analyses of time series of groundwater heads and stream water levels. The first component depicted between 53% and 70% of the total variance in the dataset for the different years. It captured the propagation of the pressure signal induced by stream water level fluctuations throughout the adjacent aquifer. Thus it could be used as a measure of hydraulic connectivity between stream and aquifer. During the first year, the impact of stream water level fluctuations decreased with distance from the regulated river (shortcut), whereas the hydraulic connection of the oxbow to the adjacent aquifer was very low. After restoration of the stream we observed a slight but not significant increase of hydraulic connectivity in the oxbow in the second year after restoration, but no change for the former shortcut. There is some evidence that the pattern of hydraulic connectivity at the study site is by far more determined by the natural heterogeneity of hydraulic conductivities of the floodplain sediments and the initial construction of the shortcut rather than by the clogging layer in the oxbow. (C) 2015 The Authors. Published by Elsevier B.V.
In soils and sediments there is a strong coupling between local biogeochemical processes and the distribution of water, electron acceptors, acids and nutrients. Both sides are closely related and affect each other from small scale to larger scales. Soil structures such as aggregates, roots, layers or macropores enhance the patchiness of these distributions. At the same time it is difficult to access the spatial distribution and temporal dynamics of these parameter. Noninvasive imaging techniques with high spatial and temporal resolution overcome these limitations. And new non-invasive techniques are needed to study the dynamic interaction of plant roots with the surrounding soil, but also the complex physical and chemical processes in structured soils. In this study we developed an efficient non-destructive in-situ method to determine biogeochemical parameters relevant to plant roots growing in soil. This is a quantitative fluorescence imaging method suitable for visualizing the spatial and temporal pH changes around roots. We adapted the fluorescence imaging set-up and coupled it with neutron radiography to study simultaneously root growth, oxygen depletion by respiration activity and root water uptake. The combined set up was subsequently applied to a structured soil system to map the patchy structure of oxic and anoxic zones induced by a chemical oxygen consumption reaction for spatially varying water contents. Moreover, results from a similar fluorescence imaging technique for nitrate detection were complemented by a numerical modeling study where we used imaging data, aiming to simulate biodegradation under anaerobic, nitrate reducing conditions.
The spread of shrubs in Namibian savannas raises questions about the resilience of these ecosystems to global change. This makes it necessary to understand the past dynamics of the vegetation, since there is no consensus on whether shrub encroachment is a new phenomenon, nor on its main drivers. However, a lack of long-term vegetation datasets for the region and the scarcity of suitable palaeoecological archives, makes reconstructing past vegetation and land cover of the savannas a challenge.
To help meet this challenge, this study addresses three main research questions: 1) is pollen analysis a suitable tool to reflect the vegetation change associated with shrub encroachment in savanna environments? 2) Does the current encroached landscape correspond to an alternative stable state of savanna vegetation? 3) To what extent do pollen-based quantitative vegetation reconstructions reflect changes in past land cover?
The research focuses on north-central Namibia, where despite being the region most affected by shrub invasion, particularly since the 21st century, little is known about the dynamics of this phenomenon.
Field-based vegetation data were compared with modern pollen data to assess their correspondence in terms of composition and diversity along precipitation and grazing intensity gradients. In addition, two sediment cores from Lake Otjikoto were analysed to reveal changes in vegetation composition that have occurred in the region over the past 170 years and their possible drivers. For this, a multiproxy approach (fossil pollen, sedimentary ancient DNA (sedaDNA), biomarkers, compound specific carbon (δ13C) and deuterium (δD) isotopes, bulk carbon isotopes (δ13Corg), grain size, geochemical properties) was applied at high taxonomic and temporal resolution. REVEALS modelling of the fossil pollen record from Lake Otjikoto was run to quantitatively reconstruct past vegetation cover. For this, we first made pollen productivity estimates (PPE) of the most relevant savanna taxa in the region using the extended R-value model and two pollen dispersal options (Gaussian plume model and Lagrangian stochastic model). The REVEALS-based vegetation reconstruction was then validated using remote sensing-based regional vegetation data.
The results show that modern pollen reflects the composition of the vegetation well, but diversity less well. Interestingly, precipitation and grazing explain a significant amount of the compositional change in the pollen and vegetation spectra. The multiproxy record shows that a state change from open Combretum woodland to encroached Terminalia shrubland can occur over a century, and that the transition between states spans around 80 years and is characterized by a unique vegetation composition. This transition is supported by gradual environmental changes induced by management (i.e. broad-scale logging for the mining industry, selective grazing and reduced fire activity associated with intensified farming) and related land-use change. Derived environmental changes (i.e. reduced soil moisture, reduced grass cover, changes in species composition and competitiveness, reduced fire intensity) may have affected the resilience of Combretum open woodlands, making them more susceptible to change to an encroached state by stochastic events such as consecutive years of precipitation and drought, and by high concentrations of pCO2. We assume that the resulting encroached state was further stabilized by feedback mechanisms that favour the establishment and competitiveness of woody vegetation.
The REVEALS-based quantitative estimates of plant taxa indicate the predominance of a semi-open landscape throughout the 20th century and a reduction in grass cover below 50% since the 21st century associated with the spread of encroacher woody taxa. Cover estimates show a close match with regional vegetation data, providing support for the vegetation dynamics inferred from multiproxy analyses. Reasonable PPEs were made for all woody taxa, but not for Poaceae.
In conclusion, pollen analysis is a suitable tool to reconstruct past vegetation dynamics in savannas. However, because pollen cannot identify grasses beyond family level, a multiproxy approach, particularly the use of sedaDNA, is required. I was able to separate stable encroached states from mere woodland phases, and could identify drivers and speculate about related feedbacks. In addition, the REVEALS-based quantitative vegetation reconstruction clearly reflects the magnitude of the changes in the vegetation cover that occurred during the last 130 years, despite the limitations of some PPEs.
This research provides new insights into pollen-vegetation relationships in savannas and highlights the importance of multiproxy approaches when reconstructing past vegetation dynamics in semi-arid environments. It also provides the first time series with sufficient taxonomic resolution to show changes in vegetation composition during shrub encroachment, as well as the first quantitative reconstruction of past land cover in the region. These results help to identify the different stages in savanna dynamics and can be used to calibrate predictive models of vegetation change, which are highly relevant to land management.
Early Ilerdian (Early Eocene, Shallow Benthic Zones 5 and 6) carbonate systems of the Pyrenees shelf were deposited after a time of severe climatic ('Paleocene-Eocene Thermal Maximum, PETM') and phylogenetic ('Larger Foraminifer Turnover') changes. They reflect the radiation of nummulitid, alveolinid, and orbitolitid larger foraminifera after remarkable biotic changes at the end of the Paleocene, and announce their subsequent flourishing in the Middle Eocene. A paleoenvironmental model for tropical carbonate environments of this particular time interval is provided herein. During the Early Ilerdian, the inner and middle ramp deposits from Minerve, Campo and Serraduy revealed the end-member of a tropical carbonate factory with carbonate production dominated by the end-members of biotically (photo-autotrophic skeletal) controlled and biotically induced carbonate precipitation. Inner platform environments are dominated by alveolinids and in part by orbitolitids, middle platform environments are dominated by nummulitids. Corals are present, but they do not form reefs, which is a typical feature for the Eocene. Nummulite shoal complexes, which are well-known from the Middle Eocene are also absent during the studied Early Ilerdian interval, which may reflect the early evolutionary stage of this group
A partially non-ergodic ground-motion prediction equation is estimated for Europe and the Middle East. Therefore, a hierarchical model is presented that accounts for regional differences. For this purpose, the scaling of ground-motion intensity measures is assumed to be similar, but not identical in different regions. This is achieved by assuming a hierarchical model, where some coefficients are treated as random variables which are sampled from an underlying global distribution. The coefficients are estimated by Bayesian inference. This allows one to estimate the epistemic uncertainty in the coefficients, and consequently in model predictions, in a rigorous way. The model is estimated based on peak ground acceleration data from nine different European/Middle Eastern regions. There are large differences in the amount of earthquakes and records in the different regions. However, due to the hierarchical nature of the model, regions with only few data points borrow strength from other regions with more data. This makes it possible to estimate a separate set of coefficients for all regions. Different regionalized models are compared, for which different coefficients are assumed to be regionally dependent. Results show that regionalizing the coefficients for magnitude and distance scaling leads to better performance of the models. The models for all regions are physically sound, even if only very few earthquakes comprise one region.
Extra-tropical circulation systems impede poleward moisture advection by the Indian Summer Monsoon. In this context, the Himalayan range is believed to insulate the south Asian circulation from extra-tropical influences and to delineate the northern extent of the Indian Summer Monsoon in central Asia. Paleoclimatic evidence, however, suggests increased moisture availability in the Early Holocene north of the Himalayan range which is attributed to an intensification of the Indian Summer Monsoon. Nevertheless, mechanisms leading to a surpassing of the Himalayan range and the northern maximum extent of summer monsoonal influence remain unknown. Here we show that the Kunlun barrier on the northern Tibetan Plateau [similar to 36 degrees N] delimits Indian Summer Monsoon precipitation during the Holocene. The presence of the barrier relocates the insulation effect 1,000 km further north, allowing a continental low intensity branch of the Indian Summer Monsoon which is persistent throughout the Holocene. Precipitation intensities at its northern extent seem to be driven by differentiated solar heating of the Northern Hemisphere indicating dependency on energy-gradients rather than absolute radiation intensities. The identified spatial constraints of monsoonal precipitation will facilitate the prediction of future monsoonal precipitation patterns in Central Asia under varying climatic conditions.
This study aims to establish, evaluate, and apply a modern pollen-climate transfer function from the transition zone between arctic tundra and light-needled taiga in Arctic Siberia. Lacustrine samples (n = 96) from the northern Siberian lowlands of Yakutia were collected along four north-to-south transects crossing the arctic forest line. Samples span a broad temperature and precipitation gradient (mean July temperature, T-July: 7.5-18.7 degrees C; mean annual precipitation, P-ann: 114-315 mm/yr). Redundancy analyses are used to examine the relationship between the modern pollen signal and corresponding vegetation types and climate. Performance of transfer functions for T-July and P-ann were cross-validated and tested for spatial autocorrelation effects. The root mean square errors of prediction are 1.67 degrees C for T-July and 40 mm/yr for P-ann. A climate reconstruction based on fossil pollen spectra from a Siberian Arctic lake sediment core spanning the Holocene yielded cold conditions for the Late Glacial (1-2 degrees C below present T-July). Warm and moist conditions were reconstructed for the early to mid Holocene (2 degrees C higher T-July than present), and climate conditions similar to modern ones were reconstructed for the last 4000 years. In conclusion, our modern pollen data set fills the gap of existing regional calibration sets with regard to the underrepresented Siberian tundra-taiga transition zone. The Holocene climate reconstruction indicates that the temperature deviation from modern values was only moderate despite the assumed Arctic sensitivity to present climate change.
Residential assets, comprising buildings and household contents, are a major source of direct flood losses. Existing damage models are mostly deterministic and limited to particular countries or flood types. Here, we compile building-level losses from Germany, Italy and the Netherlands covering a wide range of fluvial and pluvial flood events. Utilizing a Bayesian network (BN) for continuous variables, we find that relative losses (i.e. loss relative to exposure) to building structure and its contents could be estimated with five variables: water depth, flow velocity, event return period, building usable floor space area and regional disposable income per capita. The model's ability to predict flood losses is validated for the 11 flood events contained in the sample. Predictions for the German and Italian fluvial floods were better than for pluvial floods or the 1993 Meuse river flood. Further, a case study of a 2010 coastal flood in France is used to test the BN model's performance for a type of flood not included in the survey dataset. Overall, the BN model achieved better results than any of 10 alternative damage models for reproducing average losses for the 2010 flood. An additional case study of a 2013 fluvial flood has also shown good performance of the model. The study shows that data from many flood events can be combined to derive most important factors driving flood losses across regions and time, and that resulting damage models could be applied in an open data framework.
Extreme weather events are likely to occur more often under climate change and the resulting effects on ecosystems could lead to a further acceleration of climate change. But not all extreme weather events lead to extreme ecosystem response. Here, we focus on hazardous ecosystem behaviour and identify coinciding weather conditions. We use a simple probabilistic risk assessment based on time series of ecosystem behaviour and climate conditions. Given the risk assessment terminology, vulnerability and risk for the previously defined hazard are estimated on the basis of observed hazardous ecosystem behaviour.
We apply this approach to extreme responses of terrestrial ecosystems to drought, defining the hazard as a negative net biome productivity over a 12-month period. We show an application for two selected sites using data for 1981-2010 and then apply the method to the pan-European scale for the same period, based on numerical modelling results (LPJmL for ecosystem behaviour; ERA-Interim data for climate).
Our site-specific results demonstrate the applicability of the proposed method, using the SPEI to describe the climate condition. The site in Spain provides an example of vulnerability to drought because the expected value of the SPEI is 0.4 lower for hazardous than for non-hazardous ecosystem behaviour. In northern Germany, on the contrary, the site is not vulnerable to drought because the SPEI expectation values imply wetter conditions in the hazard case than in the non-hazard case.
At the pan-European scale, ecosystem vulnerability to drought is calculated in the Mediterranean and temperate region, whereas Scandinavian ecosystems are vulnerable under conditions without water shortages. These first model- based applications indicate the conceptual advantages of the proposed method by focusing on the identification of critical weather conditions for which we observe hazardous ecosystem behaviour in the analysed data set. Application of the method to empirical time series and to future climate would be important next steps to test the approach.
Extreme weather events are likely to occur more often under climate change and the resulting effects on ecosystems could lead to a further acceleration of climate change. But not all extreme weather events lead to extreme ecosystem response. Here, we focus on hazardous ecosystem behaviour and identify coinciding weather conditions. We use a simple probabilistic risk assessment based on time series of ecosystem behaviour and climate conditions. Given the risk assessment terminology, vulnerability and risk for the previously defined hazard are estimated on the basis of observed hazardous ecosystem behaviour.
We apply this approach to extreme responses of terrestrial ecosystems to drought, defining the hazard as a negative net biome productivity over a 12-month period. We show an application for two selected sites using data for 1981-2010 and then apply the method to the pan-European scale for the same period, based on numerical modelling results (LPJmL for ecosystem behaviour; ERA-Interim data for climate).
Our site-specific results demonstrate the applicability of the proposed method, using the SPEI to describe the climate condition. The site in Spain provides an example of vulnerability to drought because the expected value of the SPEI is 0.4 lower for hazardous than for non-hazardous ecosystem behaviour. In northern Germany, on the contrary, the site is not vulnerable to drought because the SPEI expectation values imply wetter conditions in the hazard case than in the non-hazard case.
At the pan-European scale, ecosystem vulnerability to drought is calculated in the Mediterranean and temperate region, whereas Scandinavian ecosystems are vulnerable under conditions without water shortages. These first model- based applications indicate the conceptual advantages of the proposed method by focusing on the identification of critical weather conditions for which we observe hazardous ecosystem behaviour in the analysed data set. Application of the method to empirical time series and to future climate would be important next steps to test the approach.
Extracellular DNA (eDNA) is a ubiquitous biological compound in aquatic sediment and soil. Previous studies suggested that eDNA plays an important role in biogeochemical element cycling, horizontal gene transfer and stabilization of biofilm structures. Previous methods for eDNA extraction were either not suitable for oligotrophic sediments or only allowed quantification but no genetic analyses. Our procedure is based on cell detachment and eDNA liberation from sediment particles by sequential washing with an alkaline sodium phosphate buffer followed by a separation of cells and eDNA. The separated eDNA is then bound onto silica particles and purified, whereas the intracellular DNA from the separated cells is extracted using a commercial kit. The method provides extra- and intracellular DNA of high purity that is suitable for downstream applications like PCR. Extracellular DNA was extracted from organic-rich shallow sediment of the Baltic Sea, glacially influenced sediment of the Barents Sea and from the oligotrophic South Pacific Gyre. The eDNA concentration in these samples varied from 23 to 626 ng g(-1) wet weight sediment. A number of experiments were performed to verify each processing step. Although extraction efficiency is higher than other published methods, it is not fully quantitative. (C) 2014 Elsevier B.V. All rights reserved.
The computation of such synthetic GFs is computationally and operationally demanding. As a consequence, the onthe-fly recalculation of synthetic GFs in each iteration of an optimisation is time-consuming and impractical. Therefore, the pre-calculation and efficient storage of synthetic GFs on a dense grid of source to receiver combinations enables the efficient lookup and utilisation of GFs in time-critical scenarios. We present a Python-based framework and toolkit - Pyrocko-GF - that enables the pre-calculation of synthetic GF stores, which are independent of their numerical calculation method and GF transfer function. The framework aids in the creation of such GF stores by interfacing a suite of established numerical forward modelling codes in seismology (computational back ends). So far, interfaces to back ends for layered Earth model cases have been provided; however, the architecture of Pyrocko-GF is designed to cover back ends for other geometries (e.g. full 3-D heterogeneous media) and other physical quantities (e.g. gravity, pressure, tilt). Therefore, Pyrocko-GF defines an extensible GF storage format suitable for a wide range of GF types, especially handling elasticity and wave propagation problems. The framework assists with visualisations, quality control, and the exchange of GF stores, which is supported through an online platform that provides many pre-calculated GF stores for local, regional, and global studies. The Pyrocko-GF toolkit comes with a well-documented application programming interface (API) for the Python programming language to efficiently facilitate forward modelling of geophysical processes, e.g. synthetic waveforms or static displacements for a wide range of source models.
Python is at the forefront of scientific computation for seismologists and therefore should be introduced to students interested in becoming seismologists. On its own, Python is open source and well designed with extensive libraries. However, Python code can also be executed, visualized, and communicated to others with "Jupyter Notebooks". Thus, Jupyter Notebooks are ideal for teaching students Python and scientific computation. In this article, we designed an openly available Python library and collection of Jupyter Notebooks based on defined scientific computation learning goals for seismology students. The Notebooks cover topics from an introduction to Python to organizing data, earthquake catalog statistics, linear regression, and making maps. Our Python library and collection of Jupyter Notebooks are meant to be used as course materials for an upper-division data analysis course in an Earth Science Department, and the materials were tested in a Probabilistic Seismic Hazard course. However, seismologists or anyone else who is interested in Python for data analysis and map making can use these materials.
Investigation of transient soil moisture profiles yields valuable information of near- surface processes. A recently developed reconstruction algorithm based on the telegraph equation allows the inverse estimation of soil moisture profiles along coated, three rod TDR probes. Laboratory experiments were carried out to prove the results of the inversion and to understand the influence of probe rod deformation and solid objects close to the probe in heterogonous media. Differences in rod geometry can lead to serious misinterpretations in the soil moisture profile but have small influence on the average soil moisture along the probe. Solids in the integration volume have almost no effect on average soil moisture but result in locally slightly decreased moisture values. Inverted profiles obtained in a loamy soil with a clay content of about 16% were in good agreement with independent measurements.
Investigation of transient soil moisture profiles yields valuable information of near- surface processes. A recently developed reconstruction algorithm based on the telegraph equation allows the inverse estimation of soil moisture profiles along coated, three rod TDR probes. Laboratory experiments were carried out to prove the results of the inversion and to understand the influence of probe rod deformation and solid objects close to the probe in heterogeneous media. Differences in rod geometry can lead to serious misinterpretations in the soil moisture profile, but have small influence on the average soil moisture along the probe. Solids in the integration volume have almost no effect on average soil moisture, but result in locally slightly decreased moisture values. Inverted profiles obtained in a loamy soil with a clay content of about 16% were in good agreement with independent measurements.
Dynamic earthquake rupture modeling provides information on the rupture physics as the rupture velocity, frictions or tractions acting during the rupture process. Nevertheless, as often based on spatial gridded preset geometries, dynamic modeling is depending on many free parameters leading to both a high non-uniqueness of the results and large computation times. That decreases the possibilities of full Bayesian error analysis.
To assess the named problems we developed the quasi-dynamic rupture model which is presented in this work. It combines the kinematic Eikonal rupture model with a boundary element method for quasi-static slip calculation.
The orientation of the modeled rupture plane is defined by a previously performed moment tensor inversion. The simultanously inverted scalar seismic moment allows an estimation of the extension of the rupture. The modeled rupture plane is discretized by a set of rectangular boundary elements. For each boundary element an applied traction vector is defined as the boundary value.
For insights in the dynamic rupture behaviour the rupture front propagation is calculated for incremental time steps based on the 2D Eikonal equation. The needed location-dependent rupture velocity field is assumed to scale linearly with a layered shear wave velocity field.
At each time all boundary elements enclosed within the rupture front are used to calculate the quasi-static slip distribution. Neither friction nor stress propagation are considered. Therefore the algorithm is assumed to be “quasi-static”. A series of the resulting quasi-static slip snapshots can be used as a quasi-dynamic model of the rupture process.
As many a priori information is used from the earth model (shear wave velocity and elastic parameters) and the moment tensor inversion (rupture extension and orientation) our model is depending on few free parameters as the traction field, the linear factor between rupture and shear wave velocity and the nucleation point and time. Hence stable and fast modeling results are obtained as proven from the comparison to different infinite and finite static crack solutions.
First dynamic applications show promissing results. The location-dependent rise time is automatically derived by the model. Different simple kinematic models as the slip-pulse or the penny-shaped crack model can be reproduced as well as their corresponding slip rate functions. A source time function (STF) approximation calculated from the cumulative sum of moment rates of each boundary element gives results similar to theoretical and empirical known STFs.
The model was also applied to the 2015 Illapel earthquake. Using a simple rectangular rupture geometry and a 2-layered traction regime yields good estimates of both the rupture front propagation and the slip patterns which are comparable to literature results. The STF approximation shows a good fit with previously published STFs.
The quasi-dynamic rupture model is hence able to fastly calculate reproducable slip results. That allows to test full Bayesian error analysis in the future. Further work on a full seismic source inversion or even a traction field inversion can also extend the scope of our model.
Sea-surface temperatures (SSTs) based on radiolarian assemblage changes are estimated for the last 160 kyr, from a sediment core (Y9) recovered from Pukaki Saddle, northeast of Campbell Plateau. Site Y9 lies beneath Subantarctic Surface Water (SAW) immediately to the north of the Subantarctic Front (SAF), which in this region is bathymetrically constrained by the edges of Campbell Plateau and defines the northern boundary of the Antarctic Circumpolar Current (ACC). Radiolarian assemblages are characterised by an exceptionally high abundance of the Antarctic to subantarctic species Antarctissa spp. (up to 68%), especially during glacial intervals. SST estimates are derived using Factor Analysis and the Modern Analog Technique. Both methods capture the glacial-interglacial (G-I) pattern. The SST reconstructions show the changing relative influence of distinct water masses during the past G-I cycle, with major temperature variations of the order of 7-9 degrees C at glacial Terminations. Glacials (marine isotope stages (MIS) 6 and 2) are associated with particularly cool SSTs that are indicative of a more vigorous SAF/ACC and an enhancement of the inflow through Pukaki Saddle and/or frequent development of cold-core eddies at the SAF. By contrast, the influence of warmer waters and relaxation of the ACC during interglacials can be inferred from temperatures slightly warmer (e.g., mid-Holocene) and/or comparable to present day (e.g., MIS 5e). During these intervals, relatively warmer temperatures most likely indicate a higher warmcore eddy activity due to a strengthened Subtropical Front and/or a weakened inflow of cool water through Pukaki Saddle and/or an increased stratification in the Campbell Plateau region. Furthermore, the SST record is characterised by an abrupt warming at ca. 10 kyr (i.e., Termination l), the occurrence of a reversal at Termination I, and a warming event at the end of MIS 4, coinciding with the A4 event in the Byrd ice core. These characteristics, together with the pronounced G-I cycle shown by the SST estimates, suggest that Site Y9 is influenced by major oceanographic changes in the SW Pacific and responds to thermal changes at high southern latitudes. (C) 2014 Elsevier B.V. All rights reserved.
We propose a RAndom Interacting Network (RAIN) model to study the interactions between a pair of complex networks. The model involves two major steps: (i) the selection of a pair of nodes, one from each network, based on intra-network node-based characteristics, and (ii) the placement of a link between selected nodes based on the similarity of their relative importance in their respective networks. Node selection is based on a selection fitness function and node linkage is based on a linkage probability defined on the linkage scores of nodes. The model allows us to relate within-network characteristics to between-network structure. We apply the model to the interaction between the USA and Schengen airline transportation networks (ATNs). Our results indicate that two mechanisms: degree-based preferential node selection and degree-assortative link placement are necessary to replicate the observed inter-network degree distributions as well as the observed inter-network assortativity. The RAIN model offers the possibility to test multiple hypotheses regarding the mechanisms underlying network interactions. It can also incorporate complex interaction topologies. Furthermore, the framework of the RAIN model is general and can be potentially adapted to various real-world complex systems.
We study the source properties of the 2005 Kashmir earthquake and its aftershocks to unravel the seismotectonics of the NW Himalayan syntaxis. The mainshock and larger aftershocks have been simultaneously relocated using phase data. We use back-projection of high-frequency energy from multiple teleseismic arrays to model the spatio-temporal evolution of the mainshock rupture. Our analysis reveal a bilateral rupture, which initially propagated SE and then NW of the epicenter, with an average rupture velocity of similar to 2 km s(-1). The area of maximum energy release is parallel to and bound by the surface rupture. Incorporating rupture propagation and velocity, we model the mainshock as a line source using P- and SH-waveform inversion. Our result confirms that the mainshock occurred on a NE dipping (similar to 35 degrees) fault plane, with centroid depth of similar to 10 km. Integrated source time function show that majority of the energy was released in the first similar to 20 s, and was confined above the hypocenter. From waveform inverted fault dimension and seismic moment, we argue that the mainshock had an additional similar to 25 km blind rupture beyond the NW Himalayan syntaxis. Combining this with findings from previous studies, we conjecture that the blind rupture propagated NW of the syntaxis underneath a weak detachment overlain by infra-Cambrian salt layer, and terminated in a wedge thrust. All moderate-to-large aftershocks, NW of the mainshock rupture, are concentrated at the edge of the blind rupture termination. Source modeling of these aftershocks reveal thrust mechanism with centroid depths of 2-10 km, and fault planes oriented subparallel to the mainshock rupture. To study the influence of mainshock rupture on aftershock occurrence, we compute Coulomb failure stress on aftershock faults. All these aftershocks lie in the positive Coulomb stress change region. This suggest that the aftershocks have been triggered by either co-seismic or post-seismic slip on the mainshock fault.
We report on a receiver function study of the crust and upper mantle within DESERT, a multidisciplinary geophysical project to study the lithosphere across the Dead Sea Transform (DST). A temporary seismic network was operated on both sides of the DST between 2000 April and 2001 June. The depth of the Moho increases smoothly from about 30 to 34-38 km towards the east across the DST, with significant north-south variations east of the DST. These Moho depth estimates from receiver functions are consistent with results from steep-and wide-angle controlled-source techniques. Steep-angle reflections and receiver functions reveal an additional discontinuity in the lower crust, but only east of the DST. This leads to the conclusion that the internal crustal structure east and west of the DST is different. The P to S converted phases from both discontinuities at 410 and 660 km are delayed by 2 s with respect to the IASP91 global reference model. This would indicate that the transition zone is consistent with the global average, but the upper mantle above 410 km is 3-4 per cent slower than the standard earth model
The Seismic Hazard Inferred from Tectonics based on the Global Strain Rate Map (SHIFT_GSRM) earthquake forecast was designed to provide high-resolution estimates of global shallow seismicity to be used in seismic hazard assessment. This model combines geodetic strain rates with global earthquake parameters to characterize long-term rates of seismic moment and earthquake activity. Although SHIFT_GSRM properly computes seismicity rates in seismically active continental regions, it underestimates earthquake rates in subduction zones by an average factor of approximately 3. We present a complementary method to SHIFT_GSRM to more accurately forecast earthquake rates in 37 subduction segments, based on the conservation of moment principle and the use of regional interface seismicity parameters, such as subduction dip angles, corner magnitudes, and coupled seismogenic thicknesses. In seven progressive steps, we find that SHIFT_GSRM earthquake-rate underpredictions are mainly due to the utilization of a global probability function of seismic moment release that poorly captures the great variability among subduction megathrust interfaces. Retrospective test results show that the forecast is consistent with the observations during the 1 January 1977 to 31 December 2014 period. Moreover, successful pseudoprospective evaluations for the 1 January 2015 to 31 December 2018 period demonstrate the power of the regionalized earthquake model to properly estimate subduction-zone seismicity.
Damage due to floods has increased during the last few decades, and further increases are expected in several regions due to climate change and growing vulnerability. To address the projected increase in flood risk, a combination of structural and non-structural flood risk mitigation measures is considered as a promising adaptation strategy. Such a combination takes into account that flood defence systems may fail, and prepares for unexpected crisis situations via land-use planning and private damage reduction, e.g. via building precautionary measures, and disaster response. However, knowledge about damage-reducing measures is scarce and often fragmented since based on case studies. For instance, it is believed that private precautionary measures, like shielding with water shutters or building fortification, are especially effective in areas with frequent flood events and low flood water levels. However, some of these measures showed a significant damage-reducing effect also during the extreme flood event in 2002 in Germany. This review analyses potentials of land-use planning and private flood precautionary measures as components of adaptation strategies for global change. Focus is on their implementation, their damage-reducing effects and their potential contribution to address projected changes in flood risk, particularly in developed countries.
A review of source models to further the understanding of the seismicity of the Groningen field
(2022)
The occurrence of felt earthquakes due to gas production in Groningen has initiated numerous studies and model attempts to understand and quantify induced seismicity in this region. The whole bandwidth of available models spans the range from fully deterministic models to purely empirical and stochastic models. In this article, we summarise the most important model approaches, describing their main achievements and limitations. In addition, we discuss remaining open questions and potential future directions of development.
Like almost all fields of science, hydrology has benefited to a large extent from the tremendous improvements in scientific instruments that are able to collect long-time data series and an increase in available computational power and storage capabilities over the last decades. Many model applications and statistical analyses (e.g., extreme value analysis) are based on these time series. Consequently, the quality and the completeness of these time series are essential. Preprocessing of raw data sets by filling data gaps is thus a necessary procedure. Several interpolation techniques with different complexity are available ranging from rather simple to extremely challenging approaches. In this paper, various imputation methods available to the hydrological researchers are reviewed with regard to their suitability for filling gaps in the context of solving hydrological questions. The methodological approaches include arithmetic mean imputation, principal component analysis, regression-based methods and multiple imputation methods. In particular, autoregressive conditional heteroscedasticity (ARCH) models which originate from finance and econometrics will be discussed regarding their applicability to data series characterized by non-constant volatility and heteroscedasticity in hydrological contexts. The review shows that methodological advances driven by other fields of research bear relevance for a more intensive use of these methods in hydrology. Up to now, the hydrological community has paid little attention to the imputation ability of time series models in general and ARCH models in particular.
The Canadian database on contemporary crustal stress has not been revised systematically in the past two decades. Here we present the results of our new compilation that contains 514 new data records for the orientation data of maximum compressive horizontal stress and 188 data records that were re-assessed. In total the Canadian stress database has now 1667 data records, which is an increase of about 45%. From these data, a new Canadian Stress map as well as one for the Province of Alberta is presented.
To analyse the stress pattern, we use the quasi median on the circle as a smoothing algorithm that generates a smoothed stress map of the maximum compressive horizontal stress orientation on a regular grid. The newly introduced quasi interquartile range on the circle estimates the spreading of the data and is used as a measure for the wave-length of the stress pattern. The result of the hybrid wavelength analysis confirms that long spatial wavelength stress patterns (>= 1000 km) exist in large areas in Canada. The observed stress pattern is transmitted through the intra-plate regions.
The results reveal that shorter spatial wave length variation of the maximum compressive horizontal stress orientation of less than 200 km, prevails particularly in south-eastern and western Canada. Regional stress sources such as density contrasts, active fault systems, crustal structures, etc. might have a significant impact in these regions. In contrast to these variations, the observed stress pattern in the Alberta Basin is very homogeneous and mainly controlled by plate boundary forces and body forces. The influence of curvature of the Rocky Mountains salient in southern Alberta is minimal. The present-day horizontal stress orientations determined herein have important implications for the production of hydrocarbons and geothermal energy in the Alberta Basin. (C) 2014 Elsevier B.V. All rights reserved.
Cosmic-Ray neutron sensing (CRS) is a unique approach to measure soil moisture at field scale filling the gap of current methodologies. However, CRS signal is affected by all the hydrogen pools on the land surface and understanding their relative importance plays an important role for the application of the method e.g., validation of remote sensing products and data assimilation. In this study, a soil moisture scaling approach is proposed to estimate directly the correct CRS soil moisture based on the soil moisture profile measured at least in one position within the field. The approach has the advantage to avoid the need to introduce one correction for each hydrogen contribution and to estimate indirectly all the related time-varying hydrogen pools. Based on the data collected in three crop seasons, the scaling approach shows its ability to identify and to quantify the seasonal biomass water equivalent. Additionally, the analysis conducted at sub-daily time resolution is able to quantify the daily vertical redistribution of the water biomass and the rainfall interception, showing promising applications of the CRS method also for these types of measurements. Overall, the study underlines how not only soil moisture but all the specific hydrological processes in the soil-plant-atmosphere continuum should be considered for a proper evaluation of the CRS signal. For this scope, the scaling approach reveals to be a simple and pragmatic analysis that can be easily extended to other experimental sites. (C) 2015 Elsevier B.V. All rights reserved.
Changes in the steepness of river profiles or abrupt vertical steps (i.e. waterfalls) are thought to be indicative of changes in erosion rates, lithology or other factors that affect landscape evolution. These changes are referred to as knickpoints or knickzones and are pervasive in bedrock river systems. Such features are thought to reveal information about landscape evolution and patterns of erosion, and therefore their locations are often reported in the geomorphic literature. It is imperative that studies reporting knickpoints and knickzones use a reproducible method of quantifying their locations, as their number and spatial distribution play an important role in interpreting tectonically active landscapes. In this contribution we introduce a reproducible knickpoint and knickzone extraction algorithm that uses river profiles transformed by integrating drainage area along channel length (the so-called integral or chi method). The profile is then statistically segmented and the differing slopes and step changes in the elevations of these segments are used to identify knickpoints, knickzones and their relative magnitudes. The output locations of identified knickpoints and knickzones compare favourably with human mapping: we test the method on Santa Cruz Island, CA, using previously reported knickzones and also test the method against a new dataset from the Quadrilatero Ferrifero in Brazil. The algorithm allows for the extraction of varying knickpoint morphologies, including stepped, positive slope-break (concave upward) and negative slope-break knickpoints. We identify parameters that most affect the resulting knickpoint and knickzone locations and provide guidance for both usage and outputs of the method to produce reproducible knickpoint datasets.
Volcanic eruptions are often preceded by seismic activity that can be used to quantify the volcanic activity. In order to allow consistent inference of the volcanic activity state from the observed seismicity patterns, objective and time-invariant classification results achievable by automatic systems should be preferred. Most automatic classification approaches need a large preclassified data set for training the system. However, in case of a volcanic crisis, we are often confronted with a lack of training data due to insufficient prior observations. In the worst case (e. g., volcanic crisis related reconfiguration of stations), there are even no prior observations available. Finally, due to the imminent crisis there might be no time for the time-consuming process of preparing a training data set. For this reason, we have developed a novel seismic-event spotting technique in order to be less dependent on previously acquired data bases and classification schemes. We are using a learning-while-recording approach based on a minimum number of reference waveforms, thus allowing for the build-up of a classification scheme as early as interesting events have been identified. First, short-term wave-field parameters (here, polarization and spectral attributes) are extracted from a continuous seismic data stream. The sequence of multidimensional feature vectors is then used to identify a fixed number of clusters in the feature space. Based on this general description of the overall wave field by a mixture of multivariate Gaussians, we are able to learn particular event classifiers (here, hidden Markov models) from a single waveform example. To show the capabilities of this new approach we apply the algorithm to a data set recorded at Soufriere Hills volcano, Montserrat. Supported by very high classification rates, we conclude that the suggested approach provides a valuable tool for volcano monitoring systems.
Pseudotachylyte veins frequently associated with mylonites and ultramylonites occur within migmatitic paragneisses, metamonzodiorites, as well as felsic and mafic granulites at the base of the section of the Hercynian lower crust exposed in Calabria (Southern Italy). The crustal section is tectonically superposed on lower grade units. Ultramylonites and pseudotachylytes are particularly well developed in migmatitic paragneisses, whereas sparse fault-related pseudotachylytes and thin mylonite/ultramylonite bands occur in granulite-facies rocks. The presence of sillimanite and clinopyroxene in ultramylonites and mylonites indicates that relatively high-temperature conditions preceded the formation of pseudotachylytes. We have analysed pseudotachylytes from different rock types to ascertain their deep crustal origin and to better understand the relationships between brittle and ductile processes during deformation of the deeper crust. Different protoliths were selected to test how lithology controls pseudotachylyte composition and textures. In migmatites and felsic granulites, euhedral or cauliflower-shaped garnets directly crystallized from pseudotachylyte melts of near andesitic composition. This indicates that pseudotachylytes originated at deep crustal conditions (> 0.75 GPa). In mafic protoliths, quenched needle-to-feather-shaped high-alumina orthopyroxene occurs in contact with newly crystallized plagioclase. The pyroxene crystallizes in garnet-free and garnet-bearing veins. The simultaneous growth of orthopyroxene and plagioclase as well as almandine, suggests lower crustal origin, with pressures in excess of 0.85 GPa. The existence of melts of different composition in the same vein indicates the stepwise, non-equilibrium conditions of frictional melting. Melt formed and intruded into pre-existing anisotropies. In mafic granulites, brittle faulting is localized in a previously formed thin high-temperature mylonite bands. migmatitic gneisses are deformed into ultramylonite domains characterized by s-c fabric. Small grain size and fluids lowered the effective stress on the c planes favouring a seismic event and the consequent melt generation. Microstructures and ductile deformation of pseudotachylytes suggest continuous ductile flow punctuated by episodes of high-strain rate, leading to seismic events and melting.
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.
Processes driving the production, transformation and transport of methane (CH4 / in wetland ecosystems are highly complex. We present a simple calculation algorithm to separate open-water CH4 fluxes measured with automatic chambers into diffusion-and ebullition-derived components. This helps to reveal underlying dynamics, to identify potential environmental drivers and, thus, to calculate reliable CH4 emission estimates. The flux separation is based on identification of ebullition-related sudden concentration changes during single measurements. Therefore, a variable ebullition filter is applied, using the lower and upper quartile and the interquartile range (IQR). Automation of data processing is achieved by using an established R script, adjusted for the purpose of CH4 flux calculation. The algorithm was validated by performing a laboratory experiment and tested using flux measurement data (July to September 2013) from a former fen grassland site, which converted into a shallow lake as a result of rewetting. Ebullition and diffusion contributed equally (46 and 55 %) to total CH4 emissions, which is comparable to ratios given in the literature. Moreover, the separation algorithm revealed a concealed shift in the diurnal trend of diffusive fluxes throughout the measurement period. The water temperature gradient was identified as one of the major drivers of diffusive CH4 emissions, whereas no significant driver was found in the case of erratic CH4 ebullition events.
Other than commonly assumed in seismology, the phase velocity of Rayleigh waves is not necessarily a single-valued function of frequency. In fact, a single Rayleigh mode can exist with three different values of phase velocity at one frequency. We demonstrate this for the first higher mode on a realistic shallow seismic structure of a homogeneous layer of unconsolidated sediments on top of a half-space of solid rock (LOH). In the case of LOH a significant contrast to the half-space is required to produce the phenomenon. In a simpler structure of a homogeneous layer with fixed (rigid) bottom (LFB) the phenomenon exists for values of Poisson's ratio between 0.19 and 0.5 and is most pronounced for P-wave velocity being three times S-wave velocity (Poisson's ratio of 0.4375). A pavement-like structure (PAV) of two layers on top of a half-space produces the multivaluedness for the fundamental mode. Programs for the computation of synthetic dispersion curves are prone to trouble in such cases. Many of them use mode-follower algorithms which loose track of the dispersion curve and miss the multivalued section. We show results for well established programs. Their inability to properly handle these cases might be one reason why the phenomenon of multivaluedness went unnoticed in seismological Rayleigh wave research for so long. For the very same reason methods of dispersion analysis must fail if they imply wave number k(l)(omega) for the lth Rayleigh mode to be a single-valued function of frequency.. This applies in particular to deconvolution methods like phase-matched filters. We demonstrate that a slant-stack analysis fails in the multivalued section, while a Fourier-Bessel transformation captures the complete Rayleigh-wave signal. Waves of finite bandwidth in the multivalued section propagate with positive group-velocity and negative phase-velocity. Their eigenfunctions appear conventional and contain no conspicuous feature.
Mineral resource exploration and mining is an essential part of today's high-tech industry. Elements such as rare-earth elements (REEs) and copper are, therefore, in high demand. Modern exploration techniques from multiple platforms (e.g., spaceborne and airborne), to detect and map the spectral characteristics of the materials of interest, require spectral libraries as an essential reference. They include field and laboratory spectral information in combination with geochemical analyses for validation. Here, we present a collection of REE- and copper-related hyperspectral spectra with associated geochemical information. The libraries contain reflectance spectra from rare-earth element oxides, REE-bearing minerals, copper-bearing minerals and mine surface samples from the Apliki copper-gold-pyrite mine in the Republic of Cyprus. The samples were measured with the HySpex imaging spectrometers in the visible and near infrared (VNIR) and shortwave infrared (SWIR) range (400-2500 nm). The geochemical validation of each sample is provided with the reflectance spectra. The spectral libraries are openly available to assist future mineral mapping campaigns and laboratory spectroscopic analyses. The spectral libraries and corresponding geochemistry are published via GFZ Data Services with the following DOIs: https://doi.org/10.5880/GFZ.1.4.2019.004 (13 REE-bearing minerals and 16 oxide powders, Koerting et al., 2019a), https://doi.org/10.5880/GFZ.1.4.2019.003 (20 copper-bearing minerals, Koellner et al., 2019), and https://doi.org/10.5880/GFZ.1.4.2019.005 (37 copper-bearing surface material samples from the Apliki coppergold-pyrite mine in Cyprus, Koerting et al., 2019b). All spectral libraries are united and comparable by the internally consistent method of hyperspectral data acquisition in the laboratory.
A SSHAC Level 3 Probabilistic Seismic Hazard Analysis for a New-Build Nuclear Site in South Africa
(2015)
A probabilistic seismic hazard analysis has been conducted for a potential nuclear power plant site on the coast of South Africa, a country of low-to-moderate seismicity. The hazard study was conducted as a SSHAC Level 3 process, the first application of this approach outside North America. Extensive geological investigations identified five fault sources with a non-zero probability of being seismogenic. Five area sources were defined for distributed seismicity, the least active being the host zone for which the low recurrence rates for earthquakes were substantiated through investigations of historical seismicity. Empirical ground-motion prediction equations were adjusted to a horizon within the bedrock at the site using kappa values inferred from weak-motion analyses. These adjusted models were then scaled to create new equations capturing the range of epistemic uncertainty in this region with no strong motion recordings. Surface motions were obtained by convolving the bedrock motions with site amplification functions calculated using measured shear-wave velocity profiles.
During the last years Urban Green Infrastructure (UGI) has evolved as a research focus across Europe. UGI can be understood as a multifunctional network of different urban green spaces and elements contributing to urban benefits. Urban agriculture has gained increasing research interest in this context. While a strong focus has been made on functions and benefits of small scale activities, the question is still open, whether these findings can be up-scaled and transferred to the farmland scale. Furthermore, multifunctionality of urban and peri-urban agriculture is rarely being considered in the landscape context. This research aims to address these gaps and harnesses the question if agricultural landscapes – which in many European metropolitan regions provide significant spatial potential – can contribute to UGI as multifunctional green spaces. This work considers multifunctionality qualitatively based on stakeholder opinion, using a participatory research approach. This study provides new insights in peri-urban farmland potentials for UGI development, resulting into a strategy framework. Furthermore, it reflects on the role of the stakeholder involvement for `multifunctionality planning´. It suggests that it helps to define meaningful bundles of intertwined functions that interact on different scales, helping to deal with non-linearity of multiple functions and to better manage them simultaneously.
Solar wind observations show that geomagnetic storms are mainly driven by interplanetary coronal mass ejections (ICMEs) and corotating or stream interaction regions (C/SIRs). We present a binary classifier that assigns one of these drivers to 7,546 storms between 1930 and 2015 using ground‐based geomagnetic field observations only. The input data consists of the long‐term stable Hourly Magnetospheric Currents index alongside the corresponding midlatitude geomagnetic observatory time series. This data set provides comprehensive information on the global storm time magnetic disturbance field, particularly its spatial variability, over eight solar cycles. For the first time, we use this information statistically with regard to an automated storm driver identification. Our supervised classification model significantly outperforms unskilled baseline models (78% accuracy with 26[19]% misidentified interplanetary coronal mass ejections [corotating or stream interaction regions]) and delivers plausible driver occurrences with regard to storm intensity and solar cycle phase. Our results can readily be used to advance related studies fundamental to space weather research, for example, studies connecting galactic cosmic ray modulation and geomagnetic disturbances. They are fully reproducible by means of the underlying open‐source software (Pick, 2019, http://doi.org/10.5880/GFZ.2.3.2019.003)
Solar wind observations show that geomagnetic storms are mainly driven by interplanetary coronal mass ejections (ICMEs) and corotating or stream interaction regions (C/SIRs). We present a binary classifier that assigns one of these drivers to 7,546 storms between 1930 and 2015 using ground‐based geomagnetic field observations only. The input data consists of the long‐term stable Hourly Magnetospheric Currents index alongside the corresponding midlatitude geomagnetic observatory time series. This data set provides comprehensive information on the global storm time magnetic disturbance field, particularly its spatial variability, over eight solar cycles. For the first time, we use this information statistically with regard to an automated storm driver identification. Our supervised classification model significantly outperforms unskilled baseline models (78% accuracy with 26[19]% misidentified interplanetary coronal mass ejections [corotating or stream interaction regions]) and delivers plausible driver occurrences with regard to storm intensity and solar cycle phase. Our results can readily be used to advance related studies fundamental to space weather research, for example, studies connecting galactic cosmic ray modulation and geomagnetic disturbances. They are fully reproducible by means of the underlying open‐source software (Pick, 2019, http://doi.org/10.5880/GFZ.2.3.2019.003)
Response spectra are of fundamental importance in earthquake engineering and represent a standard measure in seismic design for the assessment of structural performance. However, unlike Fourier spectral amplitudes, the relationship of response spectral amplitudes to seismological source, path, and site characteristics is not immediately obvious and might even be considered counterintuitive for high oscillator frequencies. The understanding of this relationship is nevertheless important for seismic-hazard analysis. The purpose of the present study is to comprehensively characterize the variation of response spectral amplitudes due to perturbations of the causative seismological parameters. This is done by calculating the absolute parameter sensitivities (sensitivity coefficients) defined as the partial derivatives of the model output with respect to its input parameters. To derive sensitivities, we apply algorithmic differentiation (AD). This powerful approach is extensively used for sensitivity analysis of complex models in meteorology or aerodynamics. To the best of our knowledge, AD has not been explored yet in the seismic-hazard context. Within the present study, AD was successfully implemented for a proven and extensively applied simulation program for response spectra (Stochastic Method SIMulation [SMSIM]) using the TAPENADE AD tool. We assess the effects and importance of input parameter perturbations on the shape of response spectra for different regional stochastic models in a quantitative way. Additionally, we perform sensitivity analysis regarding adjustment issues of groundmotion prediction equations.
A submerged pine forest from the early Holocene in the Mecklenburg Lake District, northern Germany
(2018)
For the first time, evidence of a submerged pine forest from the early Holocene can be documented in a central European lake. Subaquatic tree stumps were discovered in Lake Giesenschlagsee at a depth of between 2 and 5m using scuba divers, side-scan sonar and a remotely operated vehicle. Several erect stumps, anchored to the ground by roots, represent an insitu record of this former forest. Botanical determination revealed the stumps to be Scots pine (Pinus sylvestris) with an individual tree age of about 80years. The trees could not be dated by means of dendrochronology, as they are older than the regional reference chronology for pine. Radiocarbon ages from the wood range from 10880 +/- 210 to 10370 +/- 130cal. a BP, which is equivalent to the mid-Preboreal to early Boreal biozones. The trees are rooted in sedge peat, which can be dated to this period as well, using pollen stratigraphical analysis. Tilting of the peat bed by 4m indicates subsidence of the ground due to local dead ice melting, causing the trees to become submerged and preserved for millennia. Together with recently detected Lateglacial insitu tree occurrences in nearby lakes, the submerged pine forest at Giesenschlagsee represents a new and highly promising type of geo-bio-archive for the wider region. Comparable insitu pine remnants occur at some terrestrial (buried setting) and marine (submerged setting) sites in northern central Europe and beyond, but they partly differ in age. In general, the insitu pine finds document shifts of the zonal boreal forest ecosystem during the late Quaternary.
Ambitious climate policies, as well as economic development, education, technological progress and less resource-intensive lifestyles, are crucial elements for progress towards the UN Sustainable Development Goals (SDGs). However, using an integrated modelling framework covering 56 indicators or proxies across all 17 SDGs, we show that they are insufficient to reach the targets. An additional sustainable development package, including international climate finance, progressive redistribution of carbon pricing revenues, sufficient and healthy nutrition and improved access to modern energy, enables a more comprehensive sustainable development pathway. We quantify climate and SDG outcomes, showing that these interventions substantially boost progress towards many aspects of the UN Agenda 2030 and simultaneously facilitate reaching ambitious climate targets. Nonetheless, several important gaps remain; for example, with respect to the eradication of extreme poverty (180 million people remaining in 2030). These gaps can be closed by 2050 for many SDGs while also respecting the 1.5 °C target and several other planetary boundaries.
Current evidence suggests that many of the major events in hominin evolution occurred in East Africa. Hence, over the past two decades, there has been intensive work undertaken to understand African palaeoclimate and tectonics in order to put together a coherent picture of how the environment of Africa has varied over the past 10 Myr. A new consensus is emerging that suggests the unusual geology and climate of East Africa created a complex, environmentally very variable setting. This new understanding of East African climate has led to the pulsed climate variability hypothesis that suggests the long-term drying trend in East Africa was punctuated by episodes of short alternating periods of extreme humidity and aridity which may have driven hominin speciation, encephalization and dispersals out of Africa. This hypothesis is unique as it provides a conceptual framework within which other evolutionary theories can be examined: first, at macro-scale comparing phylogenetic gradualism and punctuated equilibrium; second, at a more focused level of human evolution comparing allopatric speciation, aridity hypothesis, turnover pulse hypothesis, variability selection hypothesis, Red Queen hypothesis and sympatric speciation based on sexual selection. It is proposed that each one of these mechanisms may have been acting on hominins during these short periods of climate variability, which then produce a range of different traits that led to the emergence of new species. In the case of Homo erectus (sensu lato), it is not just brain size that changes but life history (shortened inter-birth intervals, delayed development), body size and dimorphism, shoulder morphology to allow thrown projectiles, adaptation to long-distance running, ecological flexibility and social behaviour. The future of evolutionary research should be to create evidence-based meta-narratives, which encompass multiple mechanisms that select for different traits leading ultimately to speciation.
High-pressure is a key feature of deep subsurface environments. High partial pressure of dissolved gasses plays an important role in microbial metabolism, because thermodynamic feasibility of many reactions depends on the concentration of reactants. For gases, this is controlled by their partial pressure, which can exceed 1 MPa at in situ conditions. Therefore, high hydrostatic pressure alone is not sufficient to recreate true deep subsurface in situ conditions, but the partial pressure of dissolved gasses has to be controlled as well. We developed an incubation system that allows for incubations at hydrostatic pressure up to 60 MPa, temperatures up to 120 degrees C, and at high gas partial pressure. The composition and partial pressure of gasses can be manipulated during the experiment. To keep costs low, the system is mainly made from off-the-shelf components with only very few custommade parts. A flexible and inert PVDF (polyvinylidene fluoride) incubator sleeve, which is almost impermeable for gases, holds the sample and separates it from the pressure fluid. The flexibility of the incubator sleeve allows for sub-sampling of the medium without loss of pressure. Experiments can be run in both static and flow-through mode. The incubation system described here is usable for versatile purposes, not only the incubation of microorganisms and determination of growth rates, but also for chemical degradation or extraction experiments under high gas saturation, e.g., fluid-gas-rock-interactions in relation to carbon dioxide sequestration. As an application of the system we extracted organic compounds from sub-bituminous coal using H2O as well as a H2O-CO2 mixture at elevated temperature (90 degrees C) and pressure (5 MPa). Subsamples were taken at different time points during the incubation and analyzed by ion chromatography. Furthermore we demonstrated the applicability of the system for studies of microbial activity, using samples from the Isis mud volcano. We could detect an increase in sulfate reduction rate upon the addition of methane to the sample.
Most South Asian countries have challenges in ensuring water, energy, and food (WEF) security, which are often interacting positively or negatively. To address these challenges, the nexus approach provides a framework to identify the interactions of the WEF sectors as an integrated system. However, most nexus studies only qualitatively discuss the interactions between these sectors. This study conducts a systematic analysis of the WEF security nexus in South Asia by using open data sources at the country scale. We analyze interactions between the WEF sectors statistically, defining positive and negative correlations between the WEF security indicators as synergies and trade-offs, respectively. By creating networks of the synergies and trade-offs, we further identify most positively and negatively influencing indicators in the WEF security nexus. We observe a larger share of trade-offs than synergies within the water and energy sectors and a larger share of synergies than trade-offs among the WEF sectors for South Asia. However, these observations vary across the South Asian countries. Our analysis highlights that strategies on promoting sustainable energy and discouraging fossil fuel use could have overall positive effects on the WEF security nexus in the countries. This study provides evidence for considering the WEF security nexus as an integrated system rather than just a combination of three different sectors or securities.
Sustainable development goals (SDGs) have set the 2030 agenda to transform our world by tackling multiple challenges humankind is facing to ensure well-being, economic prosperity, and environmental protection. In contrast to conventional development agendas focusing on a restricted set of dimensions, the SDGs provide a holistic and multidimensional view on development. Hence, interactions among the SDGs may cause diverging results. To analyze the SDG interactions we systematize the identification of synergies and trade-offs using official SDG indicator data for 227 countries. A significant positive correlation between a pair of SDG indicators is classified as a synergy while a significant negative correlation is classified as a trade-off. We rank synergies and trade-offs between SDGs pairs on global and country scales in order to identify the most frequent SDG interactions. For a given SDG, positive correlations between indicator pairs were found to outweigh the negative ones in most countries. Among SDGs the positive and negative correlations between indicator pairs allowed for the identification of particular global patterns. SDG 1 (No poverty) has synergetic relationship with most of the other goals, whereas SDG 12 (Responsible consumption and production) is the goal most commonly associated with trade-offs. The attainment of the SDG agenda will greatly depend on whether the identified synergies among the goals can be leveraged. In addition, the highlighted trade-offs, which constitute obstacles in achieving the SDGs, need to be negotiated and made structurally nonobstructive by deeper changes in the current strategies.
A tale of shifting relations
(2021)
Understanding the dynamics between the East Asian summer (EASM) and winter monsoon (EAWM) is needed to predict their variability under future global warming scenarios. Here, we investigate the relationship between EASM and EAWM as well as the mechanisms driving their variability during the last 10,000 years by stacking marine and terrestrial (non-speleothem) proxy records from the East Asian realm. This provides a regional and proxy independent signal for both monsoonal systems. The respective signal was subsequently analysed using a linear regression model. We find that the phase relationship between EASM and EAWM is not time-constant and significantly depends on orbital configuration changes. In addition, changes in the Atlantic Meridional Overturning circulation, Arctic sea-ice coverage, El Niño-Southern Oscillation and Sun Spot numbers contributed to millennial scale changes in the EASM and EAWM during the Holocene. We also argue that the bulk signal of monsoonal activity captured by the stacked non-speleothem proxy records supports the previously argued bias of speleothem climatic archives to moisture source changes and/or seasonality.
A tale of shifting relations
(2021)
Understanding the dynamics between the East Asian summer (EASM) and winter monsoon (EAWM) is needed to predict their variability under future global warming scenarios. Here, we investigate the relationship between EASM and EAWM as well as the mechanisms driving their variability during the last 10,000 years by stacking marine and terrestrial (non-speleothem) proxy records from the East Asian realm. This provides a regional and proxy independent signal for both monsoonal systems. The respective signal was subsequently analysed using a linear regression model. We find that the phase relationship between EASM and EAWM is not time-constant and significantly depends on orbital configuration changes. In addition, changes in the Atlantic Meridional Overturning circulation, Arctic sea-ice coverage, El Niño-Southern Oscillation and Sun Spot numbers contributed to millennial scale changes in the EASM and EAWM during the Holocene. We also argue that the bulk signal of monsoonal activity captured by the stacked non-speleothem proxy records supports the previously argued bias of speleothem climatic archives to moisture source changes and/or seasonality.
Water research is introduced from the combined perspectives of natural and social science and cases of citizen and stakeholder coproduction of knowledge. Using the overarching notion of transdisciplinarity, we examine how interdisciplinary and participatory water research has taken place and could be developed further. It becomes apparent that water knowledge is produced widely within society, across certified disciplinary experts and noncertified expert stakeholders and citizens. However, understanding and management interventions may remain partial, or even conflicting, as much research across and between traditional disciplines has failed to integrate disciplinary paradigms due to philosophical, methodological, and communication barriers. We argue for more agonistic relationships that challenge both certified and noncertified knowledge productively. These should include examination of how water research itself embeds and is embedded in social context and performs political work. While case studies of the cultural and political economy of water knowledge exist, we need more empirical evidence on how exactly culture, politics, and economics have shaped this knowledge and how and at what junctures this could have turned out differently. We may thus channel the coproductionist critique productively to bring perspectives, alternative knowledges, and implications into water politics where they were not previously considered; in an attempt to counter potential lock-in to particular water policies and technologies that may be inequitable, unsustainable, or unacceptable. While engaging explicitly with politics, transdisciplinary water research should remain attentive to closing down moments in the research process, such as framings, path-dependencies, vested interests, researchers’ positionalities, power, and scale.
A transparent and data-driven global tectonic regionalization model for seismic hazard assessment
(2018)
A key concept that is common to many assumptions inherent within seismic hazard assessment is that of tectonic similarity. This recognizes that certain regions of the globe may display similar geophysical characteristics, such as in the attenuation of seismic waves, the magnitude scaling properties of seismogenic sources or the seismic coupling of the lithosphere. Previous attempts at tectonic regionalization, particularly within a seismic hazard assessment context, have often been based on expert judgements; in most of these cases, the process for delineating tectonic regions is neither reproducible nor consistent from location to location. In this work, the regionalization process is implemented in a scheme that is reproducible, comprehensible from a geophysical rationale, and revisable when new relevant data are published. A spatial classification-scheme is developed based on fuzzy logic, enabling the quantification of concepts that are approximate rather than precise. Using the proposed methodology, we obtain a transparent and data-driven global tectonic regionalization model for seismic hazard applications as well as the subjective probabilities (e.g. degree of being active/degree of being cratonic) that indicate the degree to which a site belongs in a tectonic category.
The central European Bohemian Massif has undergone over two centuries of scientific investigation which has made it a pivotal area for the development and testing of modern geological theories. The discovery of melt inclusions in high-grade rocks, either crystallized as nanogranitoids or as glassy inclusions, prompted the re-evaluation of the area with an ‘inclusionist’ eye. Melt inclusions have been identified in a wide range of rocks, including felsic/perpotassic granulites, migmatites, eclogites and garnet clinopyroxenites, all the result of melting events albeit over a wide range of pressure/temperature conditions (800–1000°C/0.5–5 GPa). This contribution provides an overview of such inclusions and discusses the qualitative and quantitative constraints they provide for melting processes, and the nature of melts and fluids involved in these processes. In particular, data on trace-element signatures of melt inclusions trapped at mantle depths are presented and discussed. Moreover, experimental re-homogenization of nanogranitoids provided microstructural criteria allowing assessment of the conditions at which melt and host are mutually stable during melting. Overall this work aims to provide guidelines and suggestions for petrologists wishing to explore the fascinating field of melt inclusions in metamorphic terranes worldwide, based on the newest discoveries from the still-enigmatic Bohemian Massif.
We analyze the spatiotemporal evolution of seismicity during a sequence of moderate (an M-w 4.7 foreshock and M-w 5.8 mainshock) earthquakes occurring in September 2019 at the transition between a creeping and a locked segment of the North Anatolian fault in the central Sea of Marmara, northwest Turkey. To investigate in detail the seismicity evolution, we apply a matched-filter technique to continuous waveforms, thus reducing the magnitude threshold for detection. Sequences of foreshocks preceding the two largest events are clearly seen, exhibiting two different behaviors: a long-term activation of the seismicity along the entire fault segment and a short-term concentration around the epicenters of the large events. We suggest a two-scale preparation phase, with aseismic slip preparing the mainshock final rupture a few days before, and a cascade mechanism leading to the nucleation of the mainshock. Thus, our study shows a combination of seismic and aseismic slip during the foreshock sequence changing the strength of the fault, bringing it closer to failure.
We present a new algorithm for solving the common problem of flow trapped in closed depressions within digital elevation models, as encountered in many applications relying on flow routing. Unlike other approaches (e.g., the Priority-Flood depression filling algorithm), this solution is based on the explicit computation of the flow paths both within and across the depressions through the construction of a graph connecting together all adjacent drainage basins. Although this represents many operations, a linear time complexity can be reached for the whole computation, making it very efficient. Compared to the most optimized solutions proposed so far, we show that this algorithm of flow path enforcement yields the best performance when used in landscape evolution models. In addition to its efficiency, our proposed method also has the advantage of letting the user choose among different strategies of flow path enforcement within the depressions (i.e., filling vs. carving). Furthermore, the computed graph of basins is a generic structure that has the potential to be reused for solving other problems as well, such as the simulation of erosion. This sequential algorithm may be helpful for those who need to, e.g., process digital elevation models of moderate size on single computers or run batches of simulations as part of an inference study.
The first 1400-year floating varve chronology for north-eastern Germany covering the late Allered to the early Holocene has been established by microscopic varve counts from the Rehwiese palaeolake sediment record. The Laacher See Tephra (LST), at the base of the studied interval, forms the tephrochronological anchor point. The fine laminations were examined using a combination of micro-facies and mu XRF analyses and are typical of calcite varves, which in this case provide mainly a warm season signal. Two varve types with different sub-layer structures have been distinguished: (I) complex varves consisting of up to four seasonal sub-layers formed during the Allered and early Holocene periods, and, (II) simple two sub-layer type varves only occurring during the Younger Dryas. The precision of the chronology has been improved by varve-to-varve comparison of two independently analyzed sediment profiles based on well-defined micro-marker layers. This has enabled both (1) the precise location of single missing varies in one of the sediment profiles, and, (2) the verification of varve interpolation in disturbed varve intervals in the parallel core. Inter-annual and decadal-scale variability in sediment deposition processes were traced by multi-proxy data series including seasonal layer thickness, high-resolution element scans and total organic and inorganic carbon data at a five-varve resolution. These data support the idea of a two-phase Younger Dryas, with the first interval (12,675-12,275 varve years BP) characterised by a still significant but gradually decreasing warm-season calcite precipitation and a second phase (12,275-11,690 varve years BP) with only weak calcite precipitation. Detailed correlation of these two phases with the Meerfelder Maar record based on the LST isochrone and independent varve counts provides clues about regional differences and seasonal aspects of YD climate change along a transect from a location proximal to the North Atlantic in the west to a more continental site in the east
The interplay between topography and Indian summer monsoon circulation profoundly controls precipitation distribution, sediment transport, and river discharge along the Southern Himalayan Mountain Front (SHF). The Higher Himalayas form a major orographic barrier that separates humid sectors to the south and and regions to the north. During the Indian summer monsoon, vortices transport moisture from the Bay of Bengal, swirl along the SHF to the northwest, and cause heavy rainfall when colliding with the mountain front. In the eastern and central parts of the Himalaya, precipitation measurements derived from passive microwave analysis (SSM/I) show a strong gradient, with high values at medium elevations and extensive penetration of moisture along major river valleys into the orogen. The end of the monsoonal conveyer belt is near the Sutlej Valley in the NW Himalaya, where precipitation is lower and rainfall maxima move to lower elevations. This region thus comprises a climatic transition zone that is very sensitive to changes in Indian summer monsoon strength. To constrain magnitude, temporal, and spatial distribution of precipitation, we analyzed high-resolution passive microwave data from the last decade and identified an abnormal monsoon year (AMY) in 2002. During the 2002 AMY, violent rainstorms conquered orographic barriers and penetrated far into otherwise and regions in the northwest Himalaya at elevations in excess of 3 km asl. While precipitation in these regions was significantly increased and triggered extensive erosional processes (i.e., debris flows) on sparsely vegetated, steep hillslopes, mean rainfall along the low to medium elevations was not significantly greater in magnitude. This shift may thus play an important role in the overall sediment flux toward the Himalayan foreland. Using extended precipitation and sediment flux records for the last century, we show that these events have a decadal recurrence interval during the present-day monsoon circulation. Hence, episodically occurring AMYs control geomorphic processes primarily in the high-elevation and sectors of the orogen, while annual recurring monsoonal rainfall distribution dominates erosion in the low- to medium- elevation parts along the SHF. (C) 2004 Elsevier B.V. All rights reserved
Different lake systems might reflect different climate elements of climate changes, while the responses of lake systems are also divers, and are not completely understood so far. Therefore, a comparison of lakes in different climate zones, during the high-amplitude and abrupt climate fluctuations of the Last Glacial to Holocene transition provides an exceptional opportunity to investigate distinct natural lake system responses to different abrupt climate changes. The aim of this doctoral thesis was to reconstruct climatic and environmental fluctuations down to (sub-) annual resolution from two different lake systems during the Last Glacial-Interglacial transition (~17 and 11 ka). Lake Gościąż, situated in the temperate central Poland, developed in the Allerød after recession of the Last Glacial ice sheets. The Dead Sea is located in the Levant (eastern Mediterranean) within a steep gradient from sub-humid to hyper-arid climate, and formed in the mid-Miocene. Despite their differences in sedimentation processes, both lakes form annual laminations (varves), which are crucial for studies of abrupt climate fluctuations. This doctoral thesis was carried out within the DFG project PALEX-II (Paleohydrology and Extreme Floods from the Dead Sea ICDP Core) that investigates extreme hydro-meteorological events in the ICDP core in relation to climate changes, and ICLEA (Virtual Institute of Integrated Climate and Landscape Evolution Analyses) that intends to better the understanding of climate dynamics and landscape evolutions in north-central Europe since the Last Glacial. Further, it contributes to the Helmholtz Climate Initiative REKLIM (Regional Climate Change and Humans) Research Theme 3 “Extreme events across temporal and spatial scales” that investigates extreme events using climate data, paleo-records and model-based simulations. The three main aims were to (1) establish robust chronologies of the lakes, (2) investigate how major and abrupt climate changes affect the lake systems, and (3) to compare the responses of the two varved lakes to these hemispheric-scale climate changes.
Robust chronologies are a prerequisite for high-resolved climate and environmental reconstructions, as well as for archive comparisons. Thus, addressing the first aim, the novel chronology of Lake Gościąż was established by microscopic varve counting and Bayesian age-depth modelling in Bacon for a non-varved section, and was corroborated by independent age constrains from 137Cs activity concentration measurements, AMS radiocarbon dating and pollen analysis. The varve chronology reaches from the late Allerød until AD 2015, revealing more Holocene varves than a previous study of Lake Gościąż suggested. Varve formation throughout the complete Younger Dryas (YD) even allowed the identification of annually- to decadal-resolved leads and lags in proxy responses at the YD transitions.
The lateglacial chronology of the Dead Sea (DS) was thus far mainly based on radiocarbon and U/Th-dating. In the unique ICDP core from the deep lake centre, continuous search for cryptotephra has been carried out in lateglacial sediments between two prominent gypsum deposits – the Upper and Additional Gypsum Units (UGU and AGU, respectively). Two cryptotephras were identified with glass analyses that correlate with tephra deposits from the Süphan and Nemrut volcanoes indicating that the AGU is ~1000 years younger than previously assumed, shifting it into the YD, and the underlying varved interval into the Bølling/Allerød, contradicting previous assumptions.
Using microfacies analyses, stable isotopes and temperature reconstructions, the second aim was achieved at Lake Gościąż. The YD lake system was dynamic, characterized by higher aquatic bioproductivity, more re-suspended material and less anoxia than during the Allerød and Early Holocene, mainly influenced by stronger water circulation and catchment erosion due to stronger westerly winds and less lake sheltering. Cooling at the YD onset was ~100 years longer than the final warming, while environmental proxies lagged the onset of cooling by ~90 years, but occurred contemporaneously during the termination of the YD. Chironomid-based temperature reconstructions support recent studies indicating mild YD summer temperatures. Such a comparison of annually-resolved proxy responses to both abrupt YD transitions is rare, because most European lake archives do not preserve varves during the YD.
To accomplish the second aim at the DS, microfacies analyses were performed between the UGU (~17 ka) and Holocene onset (~11 ka) in shallow- (Masada) and deep-water (ICDP core) environments. This time interval is marked by a huge but fluctuating lake level drop and therefore the complete transition into the Holocene is only recorded in the deep-basin ICDP core. In this thesis, this transition was investigated for the first time continuously and in detail. The final two pronounced lake level drops recorded by deposition of the UGU and AGU, were interrupted by one millennium of relative depositional stability and a positive water budget as recorded by aragonite varve deposition interrupted by only a few event layers. Further, intercalation of aragonite varves between the gypsum beds of the UGU and AGU shows that these generally dry intervals were also marked by decadal- to centennial-long rises in lake level. While continuous aragonite varves indicate decadal-long stable phases, the occurrence of thicker and more frequent event layers suggests general more instability during the gypsum units. These results suggest a pattern of complex and variable hydroclimate at different time scales during the Lateglacial at the DS.
The third aim was accomplished based on the individual studies above that jointly provide an integrated picture of different lake responses to different climate elements of hemispheric-scale abrupt climate changes during the Last Glacial-Interglacial transition. In general, climatically-driven facies changes are more dramatic in the DS than at Lake Gościąż. Further, Lake Gościąż is characterized by continuous varve formation nearly throughout the complete profile, whereas the DS record is widely characterized by extreme event layers, hampering the establishment of a continuous varve chronology. The lateglacial sedimentation in Lake Gościąż is mainly influenced by westerly winds and minor by changes in catchment vegetation, whereas the DS is primarily influenced by changes in winter precipitation, which are caused by temperature variations in the Mediterranean. Interestingly, sedimentation in both archives is more stable during the Bølling/Allerød and more dynamic during the YD, even when sedimentation processes are different.
In summary, this doctoral thesis presents seasonally-resolved records from two lake archives during the Lateglacial (ca 17-11 ka) to investigate the impact of abrupt climate changes in different lake systems. New age constrains from the identification of volcanic glass shards in the lateglacial sediments of the DS allowed the first lithology-based interpretation of the YD in the DS record and its comparison to Lake Gościąż. This highlights the importance of the construction of a robust chronology, and provides a first step for synchronization of the DS with other eastern Mediterranean archives. Further, climate reconstructions from the lake sediments showed variability on different time scales in the different archives, i.e. decadal- to millennial fluctuations in the lateglacial DS, and even annual variations and sub-decadal leads and lags in proxy responses during the rapid YD transitions in Lake Gościąż. This showed the importance of a comparison of different lake archives to better understand the regional and local impacts of hemispheric-scale climate variability. An unprecedented example is demonstrated here of how different lake systems show different lake responses and also react to different climate elements of abrupt climate changes. This further highlights the importance of the understanding of the respective lake system for climate reconstructions.
Abrupt or gradual?
(2018)
We used a change point analysis on a late Pleistocene-Holocene lake-sediment record from the Chew Bahir basin in the southern Ethiopian Rift to determine the amplitude and duration of past climate transitions. The most dramatic changes occurred over 240 yr (from similar to 15,700 to 15,460 yr) during the onset of the African Humid Period (AHP), and over 990 yr (from similar to 4875 to 3885 yr) during its protracted termination. The AHP was interrupted by a distinct dry period coinciding with the high-latitude Younger Dryas stadial, which had an abrupt onset (less than similar to 100 yr) at similar to 13,260 yr and lasted until similar to 11,730 yr. Wet-dry-wet transitions prior to the AHP may reflect the high-latitude Dansgaard-Oeschger cycles, as indicated by cross-correlation of the potassium record with the NorthGRIP ice core record between similar to 45-20 ka. These findings may contribute to the debates regarding the amplitude, and duration and mechanisms of past climate transitions, and their possible influence on the development of early modern human cultures.
Identifying abrupt transitions is a key question in various disciplines. Existing transition detection methods, however, do not rigorously account for time series uncertainties, often neglecting them altogether or assuming them to be independent and qualitatively similar. Here, we introduce a novel approach suited to handle uncertainties by representing the time series as a time-ordered sequence of probability density functions. We show how to detect abrupt transitions in such a sequence using the community structure of networks representing probabilities of recurrence. Using our approach, we detect transitions in global stock indices related to well-known periods of politico-economic volatility. We further uncover transitions in the El Niño-Southern Oscillation which coincide with periods of phase locking with the Pacific Decadal Oscillation. Finally, we provide for the first time an ‘uncertainty-aware’ framework which validates the hypothesis that ice-rafting events in the North Atlantic during the Holocene were synchronous with a weakened Asian summer monsoon.
Environmental parameters such as rainfall, temperature and relative humidity can affect the composition of higher plant leaf wax. The abundance and distribution of leaf wax biomarkers, such as long chain n-alkanes, in sedimentary archives have therefore been proposed as proxies reflecting climate change. However, a robust palaeoclimatic interpretation requires a thorough understanding of how environmental changes affect leaf wax n-alkane distributions in living plants. We have analysed the concentration and chain length distribution of leaf wax n-alkanes in Acacia and Eucalyptus species along a 1500 km climatic gradient in northern Australia that ranges from subtropical to arid. We show that aridity affected the concentration and distribution of n-alkanes for plants in both genera. For both Acacia and Eucalyptus n-alkane concentration increased by a factor of ten to the dry centre of Australia, reflecting the purpose of the wax in preventing water loss from the leaf. Furthermore, Acacian-alkanes decreased in average chain length (ACL) towards the arid centre of Australia, whereas Eucalyptus ACL increased under arid conditions. Our observations demonstrate that n-alkane concentration and distribution in leaf wax are sensitive to hydroclimatic conditions. These parameters could therefore potentially be employed in palaeorecords to estimate past environmental change. However, our finding of a distinct response of n-alkane ACL values to hydrological changes in different taxa also implies that the often assumed increase in ACL under drier conditions is not a robust feature for all plant species and genera and as such additional information about the prevalent vegetation are required when ACL values are used as a palaeoclimate proxy.
Accessory minerals of the Caledonian Rumburk granite are investigated to gain insight into its magmatic and post-magmatic evolution history. Recent geothermometers calibrated for trace elements in rutile (Zr), zircon (Ti), and quartz (Ti) were used to determine mineral-formation temperatures, which are compared with T data obtained from melt and fluid-inclusion Studies on quartz. Improved electron-microprobe analytical conditions allowed distinguishing several generations of rutile. Submicron-sized rutile needles included in quartz crystallized at around 739 +/- 13 degrees C and, thus, are evidently magmatic. Simultaneous crystallization of the high-T rutile and quartz is the favoured concept compared with an exsolution model for the needles. Th-U-total Pb dating of xenotime-(Y) by electron microprobe yielded a bimodal age distribution of 494 +/- 8 Ma (2 sigma; n = 44) and 311 +/- 8 Ma (2 sigma; n = 48), which is missing in monazite-(Ce). The older age correlates with the early Ordovician granite emplacement age Suggested by earlier isotopic Studies. The younger Carboniferous age also may be geologically reasonable, because the granite experienced a minor tectonothermal overprint during the Variscan orogenesis. However, whether this event has caused the resetting of the isotopic system in the xenotime is uncertain. This also holds for the age of the partial breakdown of monazite and xenotime into reaction coronas composed of fluorapatite, allanite-(Ce), epidote +/- clinozoisite. This alteration assemblage was likely produced already during autometasomatic reworking of the solidifying magma in Ordovician time, but it cannot be excluded that it relates to a Carboniferous fluid imprint connected with late-Variscan processes.
Vermetid reefs and rocky shores are hot spots of biodiversity, often referred to as the subtropical equivalent of coral reefs. The development of the ecosystem depends on the activity of several reef builders, including red crustose coralline algae (CCA) such as Neogoniolithon brassica-florida. Despite its importance, little is known about Neogoniolithon sp. acclimation to rapid changes in light intensity and corresponding photosynthetic activity. To overcome the large spatial variability in the light field (due to location and the porous nature of the rocks) we grew Neogoniolithon sp. on glass slides and characterized its photosynthetic performance in response to various light intensities by following O-2 exchange and fluorescence parameters. This was also performed on rock-inhabiting thalli collected from the east Mediterranean basin. Generally, maximal photosynthetic rate was reached when Neogoniolithon sp. thalli grown under low illumination (such as in protected niches where the light intensity can be as low as 1% of surface illumination) were examined. When exposed to light intensities higher than those experienced during growth, Neogoniolithon sp. activates adaptive/protective mechanisms such as state transition and nonphotochemical fluorescence quenching and increases the dark respiration thereafter. We find that the Fv/Fm parameter (variable/maximal fluorescence) is not suitable to assess photosynthetic performance in Neogoniolithon sp. and propose using instead an alternative parameter recently developed. Our findings help to clarify why Neogoniolithon sp. is usually observed in shaded niches along the reef surfaces.
Centroid moment tensor (CMT) parameters can be estimated from seismic waveforms. Since these data indirectly observe the deformation process, CMTs are inferred as solutions to inverse problems which are generally underdetermined and require significant assumptions, including assumptions about data noise. Broadly speaking, we consider noise to include both theory and measurement errors, where theory errors are due to assumptions in the inverse problem and measurement errors are caused by the measurement process. While data errors are routinely included in parameter estimation for full CMTs, less attention has been paid to theory errors related to velocity-model uncertainties and how these affect the resulting moment-tensor (MT) uncertainties. Therefore, rigorous uncertainty quantification for CMTs may require theory-error estimation which becomes a problem of specifying noise models. Various noise models have been proposed, and these rely on several assumptions. All approaches quantify theory errors by estimating the covariance matrix of data residuals. However, this estimation can be based on explicit modelling, empirical estimation and/or ignore or include covariances. We quantitatively compare several approaches by presenting parameter and uncertainty estimates in nonlinear full CMT estimation for several simulated data sets and regional field data of the M-1 4.4, 2015 June 13 Fox Creek, Canada, event. While our main focus is at regional distances, the tested approaches are general and implemented for arbitrary source model choice. These include known or unknown centroid locations, full MTs, deviatoric MTs and double-couple MTs. We demonstrate that velocity-model uncertainties can profoundly affect parameter estimation and that their inclusion leads to more realistic parameter uncertainty quantification. However, not all approaches perform equally well. Including theory errors by estimating non-stationary (non-Toeplitz) error covariance matrices via iterative schemes during Monte Carlo sampling performs best and is computationally most efficient. In general, including velocity-model uncertainties is most important in cases where velocity structure is poorly known.
We herein reappraise the pressure-temperature (PT) evolution of the high-pressure and low-temperature (HP-LT) Tavsanli zone (western Turkey) in order to (i) better characterize rock units exhumed along a cooling subduction interface, from birth to steady state and (ii) constrain exhumation and detachment dynamics, as well as mechanical coupling between plates. Based on PT estimates and field observations three oceanic complexes are recognized between the HP-LT continental margin and the obducted ophiolite, with PT estimates ranging from incipient metamorphism to blueschist-fades conditions. PT conditions for the continental unit are reappraised to 24 kbar and similar to 500 degrees C on the basis of pseudosection modelling and Raman spectroscopy on carbonaceous material. A tentative reconstruction of the subduction zone evolution is proposed using available radiometric and palaeogeographic data and recent thermomechanical modelling. Both PT conditions and field observations point out to the slicing of km-sized units at different preferred depths along the subduction interface, thus providing constraints on the dynamics of accretion and underplating. In particular, the comparison of PT estimates for the Tavsanli zone and for other broadly similar fossil subduction settings (i.e., Oman, Corsica, New Caledonia, Franciscan, Schistes Lustres) suggests that units are detached preferentially from the slab at specific depths of 30-40 km (i.e., downdip of the seismogenic zone) and similar to 80 km. We propose that these depths are controlled by major changes in mechanical coupling along the plate interface, whereas exhumation through time would rather be controlled by large-scale geodynamic boundary conditions. (C) 2015 Elsevier B.V. All rights reserved.
Landscape evolution models (LEMs) allow the study of earth surface responses to changing climatic and tectonic forcings. While much effort has been devoted to the development of LEMs that simulate a wide range of processes, the numerical accuracy of these models has received less attention. Most LEMs use first-order accurate numerical methods that suffer from substantial numerical diffusion. Numerical diffusion particularly affects the solution of the advection equation and thus the simulation of retreating landforms such as cliffs and river knickpoints. This has potential consequences for the integrated response of the simulated landscape. Here we test a higher-order flux-limiting finite volume method that is total variation diminishing (TVD-FVM) to solve the partial differential equations of river incision and tectonic displacement. We show that using the TVD-FVM to simulate river incision significantly influences the evolution of simulated landscapes and the spatial and temporal variability of catchment-wide erosion rates. Furthermore, a two-dimensional TVD-FVM accurately simulates the evolution of landscapes affected by lateral tectonic displacement, a process whose simulation was hitherto largely limited to LEMs with flexible spatial discretization. We implement the scheme in TTLEM (TopoToolbox Landscape Evolution Model), a spatially explicit, raster-based LEM for the study of fluvially eroding landscapes in TopoToolbox 2.
The Lower Rhine Graben (Central Europe) is a prime example of a seismically active low-strain rift zone characterized by pronounced anthropogenic and climatic overprint of structures, and long recurrence intervals of large earthquakes. These factors render the identification of active faults and surface ruptures difficult. We investigated two fault scarps in the Lower Rhine Graben, to decipher their structural character, offset and potential seismogenic origin. Both scarps were modified by anthropogenic activity. The Hemmerich site lies c. 20 km SW of Cologne, along the Erft Fault. The Untermaubach site lies SW of Duren, where the Schafberg Fault projects into the Rur River valley. At the Hemmerich site, geomorphic and geophysical data, as well as exploratory coring reveal evidence of repeated normal faulting. Geophysical analysis and palaeoseismological excavation at the Untermaubach site reveal a complex fault zone in Holocene gravels characterized by subtle gravel deformation. Differentiation of tectonic and fluvial features was only possible with trenching, because fault structures and grain sizes of the sediments were below the resolution of the geophysical data. Despite these issues, our investigation demonstrates that valuable insight into past earthquakes and seismogenic deformation in a low-strain environment can be revealed using a multidisciplinary approach.
Convexities in the longitudinal profiles of actively incising rivers are typically considered to represent the morphologic signal of a transient response to external perturbations in tectonic or climatic forcing. Distinguishing such knickzones from those that may be anchored to the channel network by spatial variations in rock uplift, however, can be challenging. Here, we combine stream profile analysis, Be-10 watershed-averaged erosion rates, and numerical modeling of stream profile evolution to evaluate whether knickzones in the Abukuma massif of northeast Japan represent a temporal or spatial change in rock uplift rate in relation to forearc shortening. Knickzones in channels that drain the eastern flank of the Abukuma massif are characterized by breaks in slope-area scaling and separate low-gradient, alluvial upper-channel segments from high-gradient, deeply-incised lower channel segments. Average erosion rates inferred from Be-10 concentrations in modern sediment below knickzones exceed erosion rates above knickzones by 20-50%. Although profile convexities could be interpreted as a transient response to an increase in rock uplift rate associated with slip on the range-bounding fault, geologic constraints on the initiation of fault slip and the magnitude of displacement cannot be reconciled with a recent, spatially uniform increase in slip rate. Rather, we find that knickzone position, stream profile gradients, and basin averaged erosion rates are best explained by a relatively abrupt spatial increase in uplift rate localized above a flat-ramp transition in the fault system. These analyses highlight the importance of considering spatially non-uniform uplift in the interpretation of stream profile evolution and demonstrate that the adjustment of river profiles to fault displacement can provide constraints on fault geometry in actively eroding landscapes. (C) 2013 Elsevier B.V. All rights reserved.
Terrestrial mud volcanoes (TMVs) represent geochemically diverse habitats with varying sulfur sources and yet sulfur cycling in these environments remains largely unexplored. Here we characterized the sulfur-metabolizing microorganisms and activity in four TMVs in Azerbaijan. A combination of geochemical analyses, biological rate measurements and molecular diversity surveys (targeting metabolic genes aprA and dsrA and SSU ribosomal RNA) supported the presence of active sulfur-oxidizing and sulfate-reducing guilds in all four TMVs across a range of physiochemical conditions, with diversity of these guilds being unique to each TMV. The TMVs varied in potential sulfate reduction rates (SRR) by up to four orders of magnitude with highest SRR observed in sediments where in situ sulfate concentrations were highest. Maximum temperatures at which SRR were measured was 60 degrees C in two TMVs. Corresponding with these trends in SRR, members of the potentially thermophilic, spore-forming, Desulfotomaculum were detected in these TMVs by targeted 16S rRNA analysis. Additional sulfate-reducing bacterial lineages included members of the Desulfobacteraceae and Desulfobulbaceae detected by aprA and dsrA analyses and likely contributing to the mesophilic SRR measured. Phylotypes affiliated with sulfide-oxidizing Gamma- and Betaproteobacteria were abundant in aprA libraries from low sulfate TMVs, while the highest sulfate TMV harboured 16S rRNA phylotypes associated with sulfur-oxidizing Epsilonproteobacteria. Altogether, the biogeochemical and microbiological data indicate these unique terrestrial habitats support diverse active sulfur-cycling microorganisms reflecting the in situ geochemical environment.
We applied the geomorphic indices (hypsometry and stream length gradient) to evaluate the differential uplift of the central and southern Longmenshan, a mountain range characterized by rapid erosion, strong tectonic uplift, and devastating seismic hazards. The results of the geomorphic analysis indicate that the Beichuan-Yingxiu fault and the Shuangshi-Dachuan fault act as major tectonic boundaries separating areas experiencing rapid uplift from slow uplift. The results of the geomorphic analysis also suggest that the Beichuan-Yingxiu fault is the most active fault with the largest relative uplift rates compared to the rest of the faults in the Longmenshan fault system. We compared reflected relative uplift rates based on the hypsometry and stream length gradient indices with geological/geodetic absolute rates. Along-strike and across-strike variations in the hypsometry and stream length gradient correlate with the spatial patterns derived from the apatite fission track exhumation rates, the leveling-derived uplift rate, and coseismic vertical displacements during the 2008 Wenchuan earthquake. These data defined multiple fault relationships in a complex thrust zone and provided geomorphic evidence to evaluate the potential seismic hazards of the southern Longmenshan range.
Magmatic continental rifts often constitute the earliest stage of nascent plate boundaries. These extensional tectonic provinces are characterized by ubiquitous normal faulting and volcanic activity; the spatial pattern, the geometry, and the age of these normal faults can help to unravel the spatiotemporal relationships between extensional deformation, magmatism, and long-wavelength crustal deformation of continental rift provinces. This study focuses on the active faulting in the Kenya Rift of the Cenozoic East African Rift System (EARS) with a focus on the mid-Pleistocene to the present-day.
To examine the early stages of continental break-up in the EARS, this thesis presents a time-averaged minimum extension rate for the inner graben of the Northern Kenya Rift (NKR) for the last 0.5 m.y. Using the TanDEM-X digital elevation model, fault-scarp geometries and associated throws are determined across the volcano-tectonic axis of the inner graben of the NKR. By integrating existing geochronology of faulted units with new ⁴⁰Ar/³⁹Ar radioisotopic dates, time-averaged extension rates are calculated. This study reveals that in the inner graben of the NKR, the long-term extension rate based on mid-Pleistocene to recent brittle deformation has minimum values of 1.0 to 1.6 mm yr⁻¹, locally with values up to 2.0 mm yr⁻¹. In light of virtually inactive border faults of the NKR, we show that extension is focused in the region of the active volcano-tectonic axis in the inner graben, thus highlighting the maturing of continental rifting in the NKR.
The phenomenon of focused extension is further investigated with a structural analysis of the youngest volcanic manifestations of the Kenya Rift, their relationship with extensional structures, and their overprint by Holocene faulting. In this context I analyzed the fault characteristics at the ~36 ka old Menengai Caldera and adjacent areas in the Central Kenya Rift using detailed field mapping and a structure-from-motion-based DEM generated from UAV data. In general, the Holocene intra-rift normal faults are dip-slip faults which strike NNE and thus reflect the present-day tectonic stress field; however, inside Menengai caldera persistent magmatic activity and magmatic resurgence overprints these young structures significantly. The caldera is located at the center of an actively extending rift segment and this and the other volcanic edifices of the Kenya Rift may constitute nucleation points of faulting an magmatic extensional processes that ultimately lead into a future stage of magma-assisted rifting.
When viewed at the scale of the entire Kenya Rift the protracted normal faulting in this region compartmentalizes the larger rift depressions, and influences the sedimentology and the hydrology of the intra-rift basins at a scale of less than 100 km. In the present day, most of the fault-bounded sub-basins of the Kenya Rift are hydrologically isolated due to this combination of faulting and magmatic activity that has generated efficient hydrological barriers that maintain these basins as semi-independent geomorphic entities. This isolation, however, was overcome during wetter climatic conditions during the past when the basins were transiently connected. I therefore also investigated the hydrological connectivity of the rift basins during the African Humid Period of the early Holocene, when climate was wetter. With the help of DEM analysis, lake-highstand indicators, radiocarbon dating, and a review of the fossil record, two lake-river-cascades could be identified: one directed southward, and one directed northward. Both cascades connected presently isolated rift basins during the early Holocene via spillovers of lakes and incised river gorges. This hydrological connection fostered the dispersal of aquatic faunas along the rift, and in addition, the water divide between the two river systems represented the only terrestrial dispersal corridor across the Kenya Rift. The reconstruction explains isolated distributions of Nilotic fish species in Kenya Rift lakes and of Guineo-Congolian mammal species in forests east of the Kenya Rift. On longer timescales, repeated episodes of connectivity and isolation must have occurred. To address this problem I participated in research to analyze a sediment drill core from the Koora basin of the Southern Kenya Rift, which provides a paleo-environmental record of the last 1 Ma. Based on this record it can be concluded that at ~400 ka relatively stable environmental conditions were disrupted by tectonic, hydrological, and ecological changes, resulting in increasingly large and frequent fluctuations in water availability, grassland communities, and woody plant cover. The major environmental shifts reflected in the drill core data coincide with phases where volcano-tectonic activity affected the basin. This thesis therefore shows how protracted extensional tectonic processes and the resulting geomorphologic conditions can affect the hydrology, the paleo-environment and the biodiversity of extensional zones in Kenya and elsewhere.
We performed leaching tests at elevated temperatures and pressures with an Alum black shale from Bomholm, Denmark and a Posidonia black shale from Lower Saxony, Germany. The Alum shale is a carbonate free black shale with pyrite and barite, containing 74.4 mu g/g U. The Posidonia shales is a calcareous shale with pyrite but without detectable amounts of barite containing 3.6 mu g/g U. Pyrite oxidized during the tests forming sulfuric acid which lowered the pH on values between 2 and 3 of the extraction fluid from the Alum shale favoring a release of U from the Alum shale to the fluid during the short-term and in the beginning of the long-term experiments. The activity concentration of U-238 is as high as 23.9 mBq/ml in the fluid for those experiments. The release of U and Th into the fluid is almost independent of pressure. The amount of uranium in the European shales is similar to that of the Marcellus Shale in the United States but the daughter product of U-238, the Ra-226 activity concentrations in the experimentally derived leachates from the European shales are quite low in comparison to that found in industrially derived flowback fluids from the Marcellus shale. This difference could mainly be due to missing Cl in the reaction fluid used in our experiments and a lower fluid to solid ratio in the industrial plays than in the experiments due to subsequent fracking and minute cracks from which Ra can easily be released.
Flood risk management in Germany follows an integrative approach in which both private households and businesses can make an important contribution to reducing flood damage by implementing property-level adaptation measures. While the flood adaptation behavior of private households has already been widely researched, comparatively less attention has been paid to the adaptation strategies of businesses. However, their ability to cope with flood risk plays an important role in the social and economic development of a flood-prone region. Therefore, using quantitative survey data, this study aims to identify different strategies and adaptation drivers of 557 businesses damaged by a riverine flood in 2013 and 104 businesses damaged by pluvial or flash floods between 2014 and 2017. Our results indicate that a low perceived self-efficacy may be an important factor that can reduce the motivation of businesses to adapt to flood risk. Furthermore, property-owners tended to act more proactively than tenants. In addition, high experience with previous flood events and low perceived response costs could strengthen proactive adaptation behavior. These findings should be considered in business-tailored risk communication.
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.
The electrical resistivity tomography (ERT) method is widely used to investigate geological, geotechnical, and hydrogeological problems in inland and aquatic environments (i.e., lakes, rivers, and seas). The objective of the ERT method is to obtain reliable resistivity models of the subsurface that can be interpreted in terms of the subsurface structure and petrophysical properties. The reliability of the resulting resistivity models depends not only on the quality of the acquired data, but also on the employed inversion strategy. Inversion of ERT data results in multiple solutions that explain the measured data equally well. Typical inversion approaches rely on different deterministic (local) strategies that consider different smoothing and damping strategies to stabilize the inversion. However, such strategies suffer from the trade-off of smearing possible sharp subsurface interfaces separating layers with resistivity contrasts of up to several orders of magnitude. When prior information (e.g., from outcrops, boreholes, or other geophysical surveys) suggests sharp resistivity variations, it might be advantageous to adapt the parameterization and inversion strategies to obtain more stable and geologically reliable model solutions. Adaptations to traditional local inversions, for example, by using different structural and/or geostatistical constraints, may help to retrieve sharper model solutions. In addition, layer-based model parameterization in combination with local or global inversion approaches can be used to obtain models with sharp boundaries.
In this thesis, I study three typical layered near-surface environments in which prior information is used to adapt 2D inversion strategies to favor layered model solutions. In cooperation with the coauthors of Chapters 2-4, I consider two general strategies. Our first approach uses a layer-based model parameterization and a well-established global inversion strategy to generate ensembles of model solutions and assess uncertainties related to the non-uniqueness of the inverse problem. We apply this method to invert ERT data sets collected in an inland coastal area of northern France (Chapter~2) and offshore of two Arctic regions (Chapter~3). Our second approach consists of using geostatistical regularizations with different correlation lengths. We apply this strategy to a more complex subsurface scenario on a local intermountain alluvial fan in southwestern Germany (Chapter~4). Overall, our inversion approaches allow us to obtain resistivity models that agree with the general geological understanding of the studied field sites. These strategies are rather general and can be applied to various geological environments where a layered subsurface structure is expected. The flexibility of our strategies allows adaptations to invert other kinds of geophysical data sets such as seismic refraction or electromagnetic induction methods, and could be considered for joint inversion approaches.
The real-time recognition and precise correction of baseline shifts in strong-motion records is a critical issue for GPS and accelerometer combined processing. This paper proposes a method to adaptively recognize and correct baseline shifts in strong-motion records by utilizing GPS measurements using two phases Kalman filter. By defining four kinds of learning statistics and criteria, the time series of estimated baseline shifts can be divided into four time intervals: initialization, static, transient and permanent. During the time interval in which the transient baseline shift is recognized, the dynamic noise of the Kalman filter system and the length of the baseline shifts estimation window are adaptively adjusted to yield a robust integration solution. The validations from an experimental and real datasets show that acceleration baseline shifts can be precisely recognized and corrected, thus, the combined system adaptively adjusted the estimation strategy to get a more robust solution. (C) 2014 COSPAR. Published by Elsevier Ltd. All rights reserved.
The Cluster mission has produced a large data set of electron flux measurements in the Earth's magnetosphere since its launch in late 2000. Electron fluxes are measured using Research with Adaptive Particle Imaging Detector (RAPID)/Imaging Electron Spectrometer (IES) detector as a function of energy, pitch angle, spacecraft position, and time. However, no adiabatic invariants have been calculated for Cluster so far. In this paper we present a step-by-step guide to calculations of adiabatic invariants and conversion of the electron flux to phase space density (PSD) in these coordinates. The electron flux is measured in two RAPID/IES energy channels providing pitch angle distribution at energies 39.2-50.5 and 68.1-94.5 keV in nominal mode since 2004. A fitting method allows to expand the conversion of the differential fluxes to the range from 40 to 150 keV. Best data coverage for phase space density in adiabatic invariant coordinates can be obtained for values of second adiabatic invariant, K, similar to 10(2), and values of the first adiabatic invariant mu in the range approximate to 5-20 MeV/G. Furthermore, we describe the production of a new data product "LSTAR," equivalent to the third adiabatic invariant, available through the Cluster Science Archive for years 2001-2018 with 1-min resolution. The produced data set adds to the availability of observations in Earth's radiation belts region and can be used for long-term statistical purposes.
The Cluster mission has produced a large data set of electron flux measurements in the Earth's magnetosphere since its launch in late 2000. Electron fluxes are measured using Research with Adaptive Particle Imaging Detector (RAPID)/Imaging Electron Spectrometer (IES) detector as a function of energy, pitch angle, spacecraft position, and time. However, no adiabatic invariants have been calculated for Cluster so far. In this paper we present a step-by-step guide to calculations of adiabatic invariants and conversion of the electron flux to phase space density (PSD) in these coordinates. The electron flux is measured in two RAPID/IES energy channels providing pitch angle distribution at energies 39.2-50.5 and 68.1-94.5 keV in nominal mode since 2004. A fitting method allows to expand the conversion of the differential fluxes to the range from 40 to 150 keV. Best data coverage for phase space density in adiabatic invariant coordinates can be obtained for values of second adiabatic invariant, K, similar to 10(2), and values of the first adiabatic invariant mu in the range approximate to 5-20 MeV/G. Furthermore, we describe the production of a new data product "LSTAR," equivalent to the third adiabatic invariant, available through the Cluster Science Archive for years 2001-2018 with 1-min resolution. The produced data set adds to the availability of observations in Earth's radiation belts region and can be used for long-term statistical purposes.
Advance in geocomputation
(2014)
Advances in hydrogravimetry
(2023)
The interest of the hydrological community in the gravimetric method has steadily increased within the last decade. This is reflected by numerous studies from many different groups with a broad range of approaches and foci. Many of those are traditionally rather hydrology-oriented groups who recognized gravimetry as a potential added value for their hydrological investigations. While this resulted in a variety of interesting and useful findings, contributing to extend the respective knowledge and confirming the methodological potential, on the other hand, many interesting and unresolved questions emerged.
This thesis manifests efforts, analyses and solutions carried out in this regard. Addressing and evaluating many of those unresolved questions, the research contributes to advancing hydrogravimetry, the combination of gravimetric and hydrological methods, in showing how gravimeters are a highly useful tool for applied hydrological field research.
In the first part of the thesis, traditional setups of stationary terrestrial superconducting gravimeters are addressed. They are commonly installed within a dedicated building, the impermeable structure of which shields the underlying soil from natural exchange of water masses (infiltration, evapotranspiration, groundwater recharge). As gravimeters are most sensitive to mass changes directly beneath the meter, this could impede their suitability for local hydrological process investigations, especially for near-surface water storage changes (WSC). By studying temporal local hydrological dynamics at a dedicated site equipped with traditional hydrological measurement devices, both below and next to the building, the impact of these absent natural dynamics on the gravity observations were quantified. A comprehensive analysis with both a data-based and model-based approach led to the development of an alternative method for dealing with this limitation. Based on determinable parameters, this approach can be transferred to a broad range of measurement sites where gravimeters are deployed in similar structures. Furthermore, the extensive considerations on this topic enabled a more profound understanding of this so called umbrella effect.
The second part of the thesis is a pilot study about the field deployment of a superconducting gravimeter. A newly developed field enclosure for this gravimeter was tested in an outdoor installation adjacent to the building used to investigate the umbrella effect. Analyzing and comparing the gravity observations from both indoor and outdoor gravimeters showed performance with respect to noise and stable environmental conditions was equivalent while the sensitivity to near-surface WSC was highly increased for the field deployed instrument. Furthermore it was demonstrated that the latter setup showed gravity changes independent of the depth where mass changes occurred, given their sufficiently wide horizontal extent. As a consequence, the field setup suits monitoring of WSC for both short and longer time periods much better. Based on a coupled data-modeling approach, its gravity time series was successfully used to infer and quantify local water budget components (evapotranspiration, lateral subsurface discharge) on the daily to annual time scale.
The third part of the thesis applies data from a gravimeter field deployment for applied hydrological process investigations. To this end, again at the same site, a sprinkling experiment was conducted in a 15 x 15 m area around the gravimeter. A simple hydro-gravimetric model was developed for calculating the gravity response resulting from water redistribution in the subsurface. It was found that, from a theoretical point of view, different subsurface water distribution processes (macro pore flow, preferential flow, wetting front advancement, bypass flow and perched water table rise) lead to a characteristic shape of their resulting gravity response curve. Although by using this approach it was possible to identify a dominating subsurface water distribution process for this site, some clear limitations stood out. Despite the advantage for field installations that gravimetry is a non-invasive and integral method, the problem of non-uniqueness could only be overcome by additional measurements (soil moisture, electric resistivity tomography) within a joint evaluation. Furthermore, the simple hydrological model was efficient for theoretical considerations but lacked the capability to resolve some heterogeneous spatial structures of water distribution up to a needed scale. Nevertheless, this unique setup for plot to small scale hydrological process research underlines the high potential of gravimetery and the benefit of a field deployment.
The fourth and last part is dedicated to the evaluation of potential uncertainties arising from the processing of gravity observations. The gravimeter senses all mass variations in an integral way, with the gravitational attraction being directly proportional to the magnitude of the change and inversely proportional to the square of the distance of the change. Consequently, all gravity effects (for example, tides, atmosphere, non-tidal ocean loading, polar motion, global hydrology and local hydrology) are included in an aggregated manner. To isolate the signal components of interest for a particular investigation, all non-desired effects have to be removed from the observations. This process is called reduction. The large-scale effects (tides, atmosphere, non-tidal ocean loading and global hydrology) cannot be measured directly and global model data is used to describe and quantify each effect. Within the reduction process, model errors and uncertainties propagate into the residual, the result of the reduction. The focus of this part of the thesis is quantifying the resulting, propagated uncertainty for each individual correction. Different superconducting gravimeter installations were evaluated with respect to their topography, distance to the ocean and the climate regime. Furthermore, different time periods of aggregated gravity observation data were assessed, ranging from 1 hour up to 12 months. It was found that uncertainties were highest for a frequency of 6 months and smallest for hourly frequencies. Distance to the ocean influences the uncertainty of the non-tidal ocean loading component, while geographical latitude affects uncertainties of the global hydrological component. It is important to highlight that the resulting correction-induced uncertainties in the residual have the potential to mask the signal of interest, depending on the signal magnitude and its frequency. These findings can be used to assess the value of gravity data across a range of applications and geographic settings.
In an overarching synthesis all results and findings are discussed with a general focus on their added value for bringing hydrogravimetric field research to a new level. The conceptual and applied methodological benefits for hydrological studies are highlighted. Within an outlook for future setups and study designs, it was once again shown what enormous potential is offered by gravimeters as hydrological field tools.
Forest structure is a crucial component in the assessment of whether a forest is likely to act as a carbon sink under changing climate. Detailed 3D structural information about the tundra–taiga ecotone of Siberia is mostly missing and still underrepresented in current research due to the remoteness and restricted accessibility. Field based, high-resolution remote sensing can provide important knowledge for the understanding of vegetation properties and dynamics. In this study, we test the applicability of consumer-grade Unmanned Aerial Vehicles (UAVs) for rapid calculation of stand metrics in treeline forests. We reconstructed high-resolution photogrammetric point clouds and derived canopy height models for 10 study sites from NE Chukotka and SW Yakutia. Subsequently, we detected individual tree tops using a variable-window size local maximum filter and applied a marker-controlled watershed segmentation for the delineation of tree crowns. With this, we successfully detected 67.1% of the validation individuals. Simple linear regressions of observed and detected metrics show a better correlation (R2) and lower relative root mean square percentage error (RMSE%) for tree heights (mean R2 = 0.77, mean RMSE% = 18.46%) than for crown diameters (mean R2 = 0.46, mean RMSE% = 24.9%). The comparison between detected and observed tree height distributions revealed that our tree detection method was unable to representatively identify trees <2 m. Our results show that plot sizes for vegetation surveys in the tundra–taiga ecotone should be adapted to the forest structure and have a radius of >15–20 m to capture homogeneous and representative forest stands. Additionally, we identify sources of omission and commission errors and give recommendations for their mitigation. In summary, the efficiency of the used method depends on the complexity of the forest’s stand structure.