Refine
Year of publication
Document Type
- Article (2545) (remove)
Language
- English (2545) (remove)
Keywords
- Holocene (35)
- climate change (28)
- Earthquake source observations (19)
- Climate change (17)
- Himalaya (17)
- Pollen (17)
- Tibetan Plateau (17)
- erosion (17)
- permafrost (17)
- Seismicity and tectonics (15)
Institute
- Institut für Geowissenschaften (2545) (remove)
Southeastern Tibetan Plateau growth revealed by inverse analysis of landscape evolution model
(2022)
The Cenozoic history of the Tibetan Plateau topography is critical for understanding the evolution of the Indian-Eurasian collision, climate, and biodiversity. However, the long-term growth and landscape evolution of the Tibetan Plateau remain ambiguous, it remains unclear if plateau uplift occurred soon after the India-Asia collision in the Paleogene (similar to 50-25 Ma) or later in the Neogene (similar to 20-5 Ma). Here, we reproduce the uplift history of the southeastern Tibetan Plateau using a 2D landscape evolution model, which simultaneously solves fluvial erosion and sediment transport processes in the drainage basins of the Three Rivers region (Yangtze, Mekong, and Salween Rivers). Our model was optimized through a formal inverse analysis with 20,000 forward simulations, which aims to reconcile the transient states of the present-day river profiles. The results, compared to existing paleoelevation and thermochronologic data, suggest initially low elevations (similar to 300-500 m) during the Paleogene, followed by a gradual southeastward propagation of topographic uplift of the plateau margin.
Fast-localized electron loss, resulting from interactions with electromagnetic ion cyclotron (EMIC) waves, can produce deepening minima in phase space density (PSD) radial profiles. Here, we perform a statistical analysis of local PSD minima to quantify how readily these are associated with radiation belt depletions. The statistics of PSD minima observed over a year are compared to the Versatile Electron Radiation Belts (VERB) simulations, both including and excluding EMIC waves. The observed minima distribution can only be achieved in the simulation including EMIC waves, indicating their importance in the dynamics of the radiation belts. By analyzing electron flux depletions in conjunction with the observed PSD minima, we show that, in the heart of the outer radiation belt (L* < 5), on average, 53% of multi-MeV electron depletions are associated with PSD minima, demonstrating that fast localized loss by interactions with EMIC waves are a common and crucial process for ultra-relativistic electron populations.
The investigation of stresses, faults, structure and seismic hazards requires a good understanding and mapping of earthquake rupture and slip. Constraining the finite source of earthquakes from seismic and geodetic waveforms is challenging because the directional effects of the rupture itself are small and dynamic numerical solutions often include a large number of free parameters. The computational effort is large and therefore difficult to use in an exploratory forward modelling or inversion approach. Here, we use a simplified self-similar fracture model with only a few parameters, where the propagation of the fracture front is decoupled from the calculation of the slip. The approximative method is flexible and computationally efficient. We discuss the strengths and limitations of the model with real-case examples of well-studied earthquakes. These include the M-w 8.3 2015 Illapel, Chile, megathrust earthquake at the plate interface of a subduction zone and examples of continental intraplate strike-slip earthquakes like the M-w 7.1 2016 Kumamoto, Japan, multisegment variable slip event or the M-w 7.5 2018 Palu, Indonesia, supershear earthquake. Despite the simplicity of the model, a large number of observational features ranging from different rupture-front isochrones and slip distributions to directional waveform effects or high slip patches are easy to model. The temporal evolution of slip rate and rise time are derived from the incremental growth of the rupture and the stress drop without imposing other constraints. The new model is fast and implemented in the open-source Python seismology toolbox Pyrocko, ready to study the physics of rupture and to be used in finite source inversions.
Pressure induced structural changes in silicate melts have a great impact on their physico-chemical properties and hence on their behaviour in the deep Earth's interior. In order to gain a deeper understanding we have studied the densification mechanism in multicomponent aluminosilicate glasses (albitic and albit-diopside composition) by means of extended X-ray absorption fine structure spectroscopy coupled to a diamond anvil cell up to 164 GPa. We have monitored the structural modifications from the network-former Ge as well as the network-modifier Sr. Notably, we tracked the evolution of Ge-O and Sr-O bond lengths (RGe-O, RSr-O) and their coordination number with pressure. We show that RGe-O increases strongly up to about 32 GPa, whereas RSr-O increases only slightly up to similar to 26 GPa. We assign these extensions to the increase of the coordination number from 4 to 6 (Ge) and from similar to 6 to at least 9 (Sr). Upon further compression RGe-O and RSr-O exhibit a continuous decrease to the highest probed pressure. These bond contractions, notably of RGe-O, that are continuous and exceed the one observed in pure SiO2 and GeO2, reflect a higher structural flexibility of multi-component glasses compared to those simple systems. Particularly, the high fraction of non-bridging oxygen atoms due to the presence of Na, Sr, Ca, Mg in the studied glasses, favours the simple compression of the highly-coordinated polyhedra of Si and Ge at pressure greater than 30 GPa. This is in strong contrast to pure oxides where cation polyhedral distortions govern the densification mechanism of the glass. The results of this study demonstrate that low field-strength alkali and alkaline earth cations, ubiquitous in deep Earth's melts, have a profound influence on the densification mechanism of glasses. Our results provide important constrains for interpreting the observed low velocity anomalies at the Earth's core-mantle boundary that have been, beyond others, referred to the presence of high-density melts. The hypothesis that non-buoyant melts at the Earth's core-mantle boundary can be formed by peculiar structural transformations in melts leading to higher coordination numbers compared to their crystalline equivalents is not supported from the present observations. The present results rather suggest that if velocity anomalies are to be explained by melts, these likely have considerable differences in chemical composition to the surrounding crystalline phase assemblage.
Marine sedimentary archives are routinely used to reconstruct past environmental changes. In many cases, bioturbation and sedimentary mixing affect the proxy time-series and the age-depth relationship. While idealized models of bioturbation exist, they usually assume homogeneous mixing, thus that a single sample is representative for the sediment layer it is sampled from.
However, it is largely unknown to which extent this assumption holds for sediments used for paleoclimate reconstructions.
To shed light on
1) the age-depth relationship and its full uncertainty,
2) the magnitude of mixing processes affecting the downcore proxy variations, and
3) the representativity of the discrete sample for the sediment layer, we designed and performed a case study on South China Sea sediment material which was collected using a box corer and which covers the last glacial cycle.
Using the radiocarbon content of foraminiferal tests as a tracer of time, we characterize the spatial age-heterogeneity of sediments in a three-dimensional setup. In total, 118 radiocarbon measurements were performed on defined small- and large-volume bulk samples ( similar to 200 specimens each) to investigate the horizontal heterogeneity of the sediment. Additionally, replicated measurements on small numbers of specimens (10 x 5 specimens) were performed to assess the heterogeneity within a sample volume. Visual assessment of X-ray images and a quantitative assessment of the mixing strength show typical mixing from bioturbation corresponding to around 10 cm mixing depth.
Notably, our 3D radiocarbon distribution reveals that the horizontal heterogeneity (up to 1,250 years), contributing to the age uncertainty, is several times larger than the typically assumed radiocarbon based age-model error (single errors up to 250 years). Furthermore, the assumption of a perfectly bioturbated layer with no mixing underneath is not met.
Our analysis further demonstrates that the age-heterogeneity might be a function of sample size; smaller samples might contain single features from the incomplete mixing and are thus less representative than larger samples.
We provide suggestions for future studies, optimal sampling strategies for quantitative paleoclimate reconstructions and realistic uncertainty in age models, as well as discuss possible implications for the interpretation of paleoclimate records.
The within-site variability in site response is the randomness in site response at a given site from different earthquakes and is treated as aleatory variability in current seismic hazard/risk analyses.
In this study, we investigate the single-station variability in linear site response at K-NET and KiK-net stations in Japan using a large number of earthquake recordings.
We found that the standard deviation of the horizontal-to-vertical Fourier spectral ratio at individual sites, that is single-station horizontal-to-vertical spectral ratio (HVSR) sigma sigma(HV,s), approximates the within-site variability in site response quantified using surface-to-borehole spectral ratios (for oscillator frequencies higher than the site fundamental frequency) or empirical ground-motion models.
Based on this finding, we then utilize the single-station HVSR sigma as a convenient tool to study the site-response variability at 697 KiK-net and 1169 K-NET sites.
Our results show that at certain frequencies, stiff, rough and shallow sites, as well as small and local events tend to have a higher sigma(HV,s).
However, when being averaged over different sites, the single-station HVSR sigma, that is sigma(HV), increases gradually with decreasing frequency. In the frequency range of 0.25-25 Hz, sigma(HV) is centred at 0.23-0.43 in ln scales (a linear scale factor of 1.26-1.54) with one standard deviation of less than 0.1. sigma(HV) is quite stable across different tectonic regions, and we present a constant, as well as earthquake magnitude- and distance-dependent sigma(HV) models.
In this study we present a novel method for the automatic detection of minerals and elements using hyperspectral transmittance imaging microscopy measurements of complete thin sections (HyperTIM).
This is accomplished by using a hyperspectral camera system that operates in the visible and near-infrared (VNIR) range with a specifically designed sample holder, scanning setup, and a microscope lens.
We utilize this method on a monazite ore thin section from Steenkampskraal (South Africa), which we analyzed for the rare earth element (REE)-bearing mineral monazite ((Ce,Nd,La)PO4), with high concentrations of Nd. The transmittance analyses with the hyperspectral VNIR camera can be used to identify REE minerals and Nd in thin sections.
We propose a three-point band depth index, the Nd feature depth index (NdFD), and its related product the Nd band depth index (NdBDI), which enables automatic mineral detection and classification for the Nd-bearing monazites in thin sections. In combination with the average concentration of the relative Nd content, it permits a destruction-free, total concentration calculation for Nd across the entire thin section.
Groundwater recharge (GWR) is one of the most challenging water fluxes to estimate, as it relies on observed data that are often limited in many developing countries.
This study developed an innovative water budget method using satellite products for estimating the spatially distributed GWR at monthly and annual scales in tropical wet sedimentary regions despite cloudy conditions.
The distinctive features proposed in this study include the capacity to address 1) evapotranspiration estimations in tropical wet regions frequently overlaid by substantial cloud cover; and 2) seasonal root-zone water storage estimations in sedimentary regions prone to monthly variations.
The method also utilises satellite-based information of the precipitation and surface runoff. The GWR was estimated and validated for the hydrologically contrasting years 2016 and 2017 over a tropical wet sedimentary region located in North-eastern Brazil, which has substantial potential for groundwater abstraction.
This study showed that applying a cloud-cleaning procedure based on monthly compositions of biophysical data enables the production of a reasonable proxy for evapotranspiration able to estimate groundwater by the water budget method.
The resulting GWR rates were 219 (2016) and 302 (2017) mm yr(-1), showing good correlations (CC = 0.68 to 0.83) and slight underestimations (PBIAS =-13 to-9%) when compared with the referenced estimates obtained by the water table fluctuation method for 23 monitoring wells. Sensitivity analysis shows that water storage changes account for +19% to-22% of our monthly evaluation.
The satellite-based approach consistently demonstrated that the consideration of cloud-cleaned evapotranspiration and root-zone soil water storage changes are essential for a proper estimation of spatially distributed GWR in tropical wet sedimentary regions because of their weather seasonality and cloudy conditions.
Carbonate minerals are common in both marine and lacustrine records, and are frequently used for paleoenvironmental reconstructions. The sedimentary sequence of the endorheic Dead Sea and its precursors contain aragonite laminae that provide a detailed sedimentary archive of climatic, hydrologic, limnologic and environmental conditions since the Pleistocene. However, the interpretation of these archives requires a detailed understanding of the constraints and mechanisms affecting CaCO3 precipitation, which are still debated. The implications of aragonite precipitation in the Dead Sea and in its late Pleistocene predecessor (Lake Lisan) were investigated in this study by mixing natural and synthetic brines with a synthetic bicarbonate solution that mimics flash-floods composition, with and without the addition of extracellular polymeric substances (EPS). Aragonite precipitation was monitored, and precipitation rates and carbonate yields were calculated and are discussed with respect to modern aquatic environments. The experimental insights on aragonite precipitation are then integrated with microfacies analyses in order to reconstruct and constrain prevailing limnogeological processes and their hydroclimatic drivers under low (interglacial) and high (glacial) lake level stands. Aragonite precipitation took place within days to several weeks after the mixing of the brines with a synthetic bicarbonate solution. Incubation time was proportional to bicarbonate concentration, and precipitation rates were partially influenced by ionic strength. Additionally, extracellular polymeric substances inhibited aragonite precipitation for several months. As for the lake's water budget, our calculations suggest that the precipitation of a typical aragonite lamina (0.5 mm thick) during high lake stand requires unreasonable freshwater inflow from either surface or subsurface sources. This discrepancy can be resolved by considering one or a combination of the following scenarios; (1) discontinuous aragonite deposition over parts of the lake floor; (2) supply of additional carbonate flux (or fluxes) to the lake from aeolian dust and the remobilization and dissolution of dust deposits at the watershed; (3) carbonate production via oxidation of organic carbon by sulfate-reducing bacteria. Altogether, it is suggested that aragonite laminae thickness cannot be directly interpreted for quantitatively reconstructing the hydrological balance for the entire lake, they may still prove valuable for identifying inherent hydroclimatic periodicities at a single site.
Improving nitrogen (N) status in European water bodies is a pressing issue. N levels depend not only on current but also past N inputs to the landscape, that have accumulated through time in legacy stores (e.g., soil, groundwater).
Catchment-scale N models, that are commonly used to investigate in-stream N levels, rarely examine the magnitude and dynamics of legacy components.
This study aims to gain a better understanding of the long-term fate of the N inputs and its uncertainties, using a legacy-driven N model (ELEMeNT) in Germany's largest national river basin (Weser; 38,450 km(2)) over the period 1960-2015.
We estimate the nine model parameters based on a progressive constraining strategy, to assess the value of different observational data sets.
We demonstrate that beyond in-stream N loading, soil N content and in-stream N concentration allow to reduce the equifinality in model parameterizations.
We find that more than 50% of the N surplus denitrifies (1480-2210 kg ha(-1)) and the stream export amounts to around 18% (410-640 kg ha(-1)), leaving behind as much as around 230-780 kg ha(-1) of N in the (soil) source zone and 10-105 kg ha(-1) in the subsurface.
A sensitivity analysis reveals the importance of different factors affecting the residual uncertainties in simulated N legacies, namely hydrologic travel time, denitrification rates, a coefficient characterizing the protection of organic N in source zone and N surplus input.
Our study calls for proper consideration of uncertainties in N legacy characterization, and discusses possible avenues to further reduce the equifinality in water quality modeling.
Drainage-divide migration, controlled by rock-uplift and rainfall patterns, may play a major role in the geomorphic evolution of mountain ranges.
However, divide-migration rates over geologic timescales have only been estimated by theoretical studies and remain empirically poorly constrained.
Geomorphological evidence suggests that the Sierra de Aconquija, on the eastern side of the southern Central Andes, northwest Argentina, is undergoing active westward drainage-divide migration. The mountain range has been subjected to steep rock trajectories and pronounced orographic rainfall for the last several million years, presenting an ideal setting for using low-temperature thermochronometric data to explore its topographic evolution.
We perform three-dimensional thermal-kinematic modeling of previously published thermochronometric data spanning the windward and leeward sides of the range to explore the most likely structural and topographic evolution of the range.
We find that the data can be explained by scenarios involving drainage-divide migration alone, or by scenarios that also involve changes in the structures that have accommodated deformation through time.
By combining new Be-10-derived catchment-average denudation rates with geomorphic constraints on probable fault activity, we conclude that the evolution of the range was likely dominated by west-vergent faulting on a high-angle reverse fault underlying the range, together with westward drainage-divide migration at a rate of several km per million years.
Our findings place new constraints on the magnitudes and rates of drainage-divide migration in real landscapes, quantify the effects of orographic rainfall and erosion on the topographic evolution of a mountain range, and highlight the importance of considering drainage-divide migration when interpreting thermochronometer age patterns.
Sedimentary ancient DNA-based studies have been used to probe centuries of climate and environmental changes and how they affected cyanobacterial assemblages in temperate lakes. Due to cyanobacteria containing potential bloom-forming and toxin-producing taxa, their approximate reconstruction from sediments is crucial, especially in lakes lacking long-term monitoring data. To extend the resolution of sediment record interpretation, we used high-throughput sequencing, amplicon sequence variant (ASV) analysis, and quantitative PCR to compare pelagic cyanobacterial composition to that in sediment traps (collected monthly) and surface sediments in Lake Tiefer See. Cyanobacterial composition, species richness, and evenness was not significantly different among the pelagic depths, sediment traps and surface sediments (p > 0.05), indicating that the cyanobacteria in the sediments reflected the cyanobacterial assemblage in the water column. However, total cyanobacterial abundances (qPCR) decreased from the metalimnion down the water column. The aggregate-forming (Aphanizomenon) and colony-forming taxa (Snowella) showed pronounced sedimentation. In contrast, Planktothrix was only very poorly represented in sediment traps (meta- and hypolimnion) and surface sediments, despite its highest relative abundance at the thermocline (10 m water depth) during periods of lake stratification (May-October). We conclude that this skewed representation in taxonomic abundances reflects taphonomic processes, which should be considered in future DNA-based paleolimnological investigations.
Megathrust earthquakes impose changes of differential stress and pore pressure in the lithosphere-asthenosphere system that are transiently relaxed during the postseismic period primarily due to afterslip, viscoelastic and poroelastic processes.
Especially during the early postseismic phase, however, the relative contribution of these processes to the observed surface deformation is unclear.
To investigate this, we use geodetic data collected in the first 48 days following the 2010 Maule earthquake and a poro-viscoelastic forward model combined with an afterslip inversion.
This model approach fits the geodetic data 14% better than a pure elastic model. Particularly near the region of maximum coseismic slip, the predicted surface poroelastic uplift pattern explains well the observations.
If poroelasticity is neglected, the spatial afterslip distribution is locally altered by up to +/- 40%.
Moreover, we find that shallow crustal aftershocks mostly occur in regions of increased postseismic pore-pressure changes, indicating that both processes might be mechanically coupled.
On their way from inland to the ocean, flowing water bodies, their constituents and their biotic communities are ex-posed to complex transport and transformation processes. However, detailed process knowledge as revealed by La-grangian measurements adjusted to travel time is rare in large rivers, in particular at hydrological extremes. To fill this gap, we investigated autotrophic processes, heterotrophic carbon utilization, and micropollutant concentrations applying a Lagrangian sampling design in a 600 km section of the River Elbe (Germany) at historically low discharge. Under base flow conditions, we expect the maximum intensity of instream processes and of point source impacts. Phy-toplankton biomass and photosynthesis increased from upstream to downstream sites but maximum chlorophyll con-centration was lower than at mean discharge. Concentrations of dissolved macronutrients decreased to almost complete phosphate depletion and low nitrate values. The longitudinal increase of bacterial abundance and production was less pronounced than in wetter years and bacterial community composition changed downstream. Molecular analyses revealed a longitudinal increase of many DOM components due to microbial production, whereas saturated lipid-like DOM, unsaturated aromatics and polyphenols, and some CHOS surfactants declined. In decomposition exper-iments, DOM components with high O/C ratios and high masses decreased whereas those with low O/C ratios, low masses, and high nitrogen content increased at all sites. Radiocarbon age analyses showed that DOC was relatively old (890-1870 years B.P.), whereas the mineralized fraction was much younger suggesting predominant oxidation of algal lysis products and exudates particularly at downstream sites. Micropollutants determining toxicity for algae (terbuthylazine, terbutryn, isoproturon and lenacil), hexachlorocyclohexanes and DDTs showed higher concentrations from the middle towards the downstream part but calculated toxicity was not negatively correlated to phytoplankton. Overall, autotrophic and heterotrophic process rates and micropollutant concentrations increased from up-to down-stream reaches, but their magnitudes were not distinctly different to conditions at medium discharges.
LegacyPollen 1.0
(2022)
Here we describe the LegacyPollen 1.0, a dataset of 2831 fossil pollen records with metadata, a harmonized taxonomy, and standardized chronologies.
A total of 1032 records originate from North America, 1075 from Europe, 488 from Asia, 150 from Latin America, 54 from Africa, and 32 from the Indo-Pacific.
The pollen data cover the late Quaternary (mostly the Holocene). The original 10 110 pollen taxa names (including variations in the notations) were harmonized to 1002 terrestrial taxa (including Cyperaceae), with woody taxa and major herbaceous taxa harmonized to genus level and other herbaceous taxa to family level.
The dataset is valuable for synthesis studies of, for example, taxa areal changes, vegetation dynamics, human impacts (e.g., deforestation), and climate change at global or continental scales.
The harmonized pollen and metadata as well as the harmonization table are available from PANGAEA (https://doi.org/10.1594/PANGAEA.929773; Herzschuh et al., 2021). R code for the harmonization is provided at Zenodo (https://doi.org/10.5281/zenodo.5910972; Herzschuh et al., 2022) so that datasets at a customized harmonization level can be easily established.
Boreal forests cover over half of the global permafrost area and protect underlying permafrost. Boreal forest development, therefore, has an impact on permafrost evolution, especially under a warming climate.
Forest disturbances and changing climate conditions cause vegetation shifts and potentially destabilize the carbon stored within the vegetation and permafrost. Disturbed permafrost-forest ecosystems can develop into a dry or swampy bush- or grasslands, shift toward broadleaf- or evergreen needleleaf-dominated forests, or recover to the pre-disturbance state.
An increase in the number and intensity of fires, as well as intensified logging activities, could lead to a partial or complete ecosystem and permafrost degradation. We study the impact of forest disturbances (logging, surface, and canopy fires) on the thermal and hydrological permafrost conditions and ecosystem resilience.
We use a dynamic multilayer canopy-permafrost model to simulate different scenarios at a study site in eastern Siberia. We implement expected mortality, defoliation, and ground surface changes and analyze the interplay between forest recovery and permafrost. We find that forest loss induces soil drying of up to 44%, leading to lower active layer thicknesses and abrupt or steady decline of a larch forest, depending on disturbance intensity.
Only after surface fires, the most common disturbances, inducing low mortality rates, forests can recover and overpass pre-disturbance leaf area index values. We find that the trajectory of larch forests after surface fires is dependent on the precipitation conditions in the years after the disturbance. Dryer years can drastically change the direction of the larch forest development within the studied period.
Measuring the variability of incoming neutrons locally would be usefull for the cosmic-ray neutron sensing (CRNS) method. As the measurement of high energy neutrons is not so easy, alternative particles can be considered for such purpose. Among them, muons are particles created from the same cascade of primary cosmic-ray fluxes that generate neutrons at the ground. In addition, they can be easily detected by small and relatively inexpensive detectors. For these reasons they could provide a suitable local alternative to incoming corrections based on remote neutron monitor data. The reported measurements demonstrated that muon detection system can detect incoming cosmic-ray variations locally. Furthermore the precision of this measurement technique is considered adequate for many CRNS applications.
A large landslide (frozen debris avalanche) occurred at Assapaat on the south coast of the Nuussuaq Peninsula in Central West Greenland on June 13, 2021, at 04:04 local time. We present a compilation of available data from field observations, photos, remote sensing, and seismic monitoring to describe the event. Analysis of these data in combination with an analysis of pre- and post-failure digital elevation models results in the first description of this type of landslide. The frozen debris avalanche initiated as a 6.9 * 10(6) m(3) failure of permafrozen talus slope and underlying colluvium and till at 600-880 m elevation. It entrained a large volume of permafrozen colluvium along its 2.4 km path in two subsequent entrainment phases accumulating a total volume between 18.3 * 10(6) and 25.9 * 10(6) m(3). About 3.9 * 10(6) m(3) is estimated to have entered the Vaigat strait; however, no tsunami was reported, or is evident in the field. This is probably because the second stage of entrainment along with a flattening of slope angle reduced the mobility of the frozen debris avalanche. We hypothesise that the initial talus slope failure is dynamically conditioned by warming of the ice matrix that binds the permafrozen talus slope. When the slope ice temperature rises to a critical level, its shear resistance is reduced, resulting in an unstable talus slope prone to failure. Likewise, we attribute the large-scale entrainment to increasing slope temperature and take the frozen debris avalanche as a strong sign that the permafrost in this region is increasingly at a critical state. Global warming is enhanced in the Arctic and frequent landslide events in the past decade in Western Greenland let us hypothesise that continued warming will lead to an increase in the frequency and magnitude of these types of landslides. Essential data for critical arctic slopes such as precipitation, snowmelt, and ground and surface temperature are still missing to further test this hypothesis. It is thus strongly required that research funds are made available to better predict the change of landslide threat in the Arctic.
A novel idea for an optimal time delay state space reconstruction from uni- and multivariate time series is presented. The entire embedding process is considered as a game, in which each move corresponds to an embedding cycle and is subject to an evaluation through an objective function. This way the embedding procedure can be modeled as a tree, in which each leaf holds a specific value of the objective function. By using a Monte Carlo ansatz, the proposed algorithm populates the tree with many leafs by computing different possible embedding paths and the final embedding is chosen as that particular path, which ends at the leaf with the lowest achieved value of the objective function. The method aims to prevent getting stuck in a local minimum of the objective function and can be used in a modular way, enabling practitioners to choose a statistic for possible delays in each embedding cycle as well as a suitable objective function themselves. The proposed method guarantees the optimization of the chosen objective function over the parameter space of the delay embedding as long as the tree is sampled sufficiently. As a proof of concept, we demonstrate the superiority of the proposed method over the classical time delay embedding methods using a variety of application examples. We compare recurrence plot-based statistics inferred from reconstructions of a Lorenz-96 system and highlight an improved forecast accuracy for map-like model data as well as for palaeoclimate isotope time series. Finally, we utilize state space reconstruction for the detection of causality and its strength between observables of a gas turbine type thermoacoustic combustor.