Refine
Has Fulltext
- no (150) (remove)
Year of publication
- 2020 (150) (remove)
Document Type
- Article (140)
- Doctoral Thesis (4)
- Other (3)
- Monograph/Edited Volume (2)
- Review (1)
Is part of the Bibliography
- yes (150)
Keywords
- climate change (5)
- Himalaya (3)
- model (3)
- models (3)
- precipitation (3)
- Arctic (2)
- CO2 (2)
- China (2)
- Chinese loess (2)
- Europe (2)
Institute
- Institut für Geowissenschaften (150) (remove)
This study, based on new and high quality in situ observations, quantifies for the first time, the individual contributions of light-absorbing aerosols (black carbon (BC), brown carbon (BrC) and dust) to aerosol absorption over the Indo-Gangetic Plain (IGP) and the Himalayan foothill region, a relatively poorly studied region with several sensitive ecosystems of global importance, as well as highly vulnerable populations. The annual and seasonal average single scattering albedo (SSA) over Kathmandu is the lowest of all the locations. The SSA over Kathmandu is < 0.89 during all seasons, which confirms the dominance of light-absorbing carbonaceous aerosols from local and regional sources over Kathmandu. It is observed here that the SSA decreases with increasing elevation, confirming the dominance of light absorbing carbonaceous aerosols at higher elevations. In contrast, the SSA over the IGP does not exhibit a pronounced spatial variation. BC dominates (>= 75%) the aerosol absorption over the IGP and the Himalayan foothills throughout the year. Higher BC concentration at elevated locations in the Himalayas leads to lower SSA at elevated locations in the Himalayas. The contribution of dust to aerosol absorption is higher throughout the year over the IGP than over the Himalayan foothills. The aerosol absorption over South Asia is very high, exceeding available observations over East Asia, and also exceeds previous model estimates. This quantification will be valuable as observational constraints to help improve regional simulations of climate change, impacts on the glaciers and the hydrological cycle, and will help to direct the focus towards BC as the main contributor to aerosol-induced warming in the region.
Receiver function approaches have proven to be valuable for the investigation of crustal and upper mantle discontinuities whose sharp changes in seismic velocities cause wave conversions. While the crustal and mantle transition zone discontinuities are largely understood, the X-discontinuity at 250-350 km depth is still an object of controversial debate. The origin and global distribution of this structure with a velocity jump of 1.5-4.8% for compressional and shear waves is still unexplained. Although the crustal and mantle transition zone discontinuities beneath SW Morocco and surroundings have been investigated, only a few studies observed the X-discontinuity and place the depth at 260-370 km beneath the region of western Morocco. In order to better locate and characterize the X-discontinuity beneath southwest Morocco, we create P-wave receiver functions using data recorded by the Morocco-Munster array and detect the X-discontinuity at apparent depths of 285-350 km. In the western part of our study region we find apparent depths of similar to 310-340 km. The eastern part of the study area appears more complex: we locate two velocity jumps at apparent depths of around 285-295 km and 330-350 km in the northeast, and in the southeast we find a discontinuity at apparent depths of 340-350 km. Due to the large depth range and the twofold appearance of the X-discontinuity, we suggest that two different phase transitions cause the X-discontinuity beneath SW Morocco. The velocity contrasts at larger depths likely point to the coesite-stishovite phase transition occurring in deep eclogitic pools. The shallower depths can be explained by the transition from orthoenstatite to high-pressure clinoenstatite which requires the reaction between eclogite and peridotite to form orthopyroxene-rich peridotite. This reaction is likely related to previously proposed small-scale mantle upwellings beneath SW Morocco. Since both phase transitions require eclogite occurrence, the location of the X-discontinuity in this region can be used to indicate the location of recycled oceanic crust.
Sediment Transit Time and Floodplain Storage Dynamics in Alluvial Rivers Revealed by Meteoric 10Be
(2020)
Quantifying the time scales of sediment transport and storage through river systems is fundamental for understanding weathering processes, biogeochemical cycling, and improving watershed management, but measuring sediment transit time is challenging. Here we provide the first systematic test of measuring cosmogenic meteoric Beryllium-10 (10Bem) in the sediment load of a large alluvial river to quantify sediment transit times. We take advantage of a natural experiment in the Rio Bermejo, a lowland alluvial river traversing the east Andean foreland basin in northern Argentina. This river has no tributaries along its trunk channel for nearly 1,300 km downstream from the mountain front. We sampled suspended sediment depth profiles along the channel and measured the concentrations of 10Bem in the chemically extracted grain coatings. We calculated depth-integrated 10Bem concentrations using sediment flux data and found that 10Bem concentrations increase 230% from upstream to downstream, indicating a mean total sediment transit time of 8.4 ± 2.2 kyr. Bulk sediment budget-based estimates of channel belt and fan storage times suggest that the 10Bem tracer records mixing of old and young sediment reservoirs. On a reach scale, 10Bem transit times are shorter where the channel is braided and superelevated above the floodplain, and longer where the channel is incised and meandering, suggesting that transit time is controlled by channel morphodynamics. This is the first systematic application of 10Bem as a sediment transit time tracer and highlights the method's potential for inferring sediment routing and storage dynamics in large river systems.
Holocene temperature proxy records are commonly used in quantitative synthesis and model-data comparisons. However, comparing correlations between time series from records collected in proximity to one another with the expected correlations based on climate model simulations indicates either regional or noisy climate signals in Holocene temperature proxy records. In this study, we evaluate the consistency of spatial correlations present in Holocene proxy records with those found in data from the Last Glacial Maximum (LGM). Specifically, we predict correlations expected in LGM proxy records if the only difference to Holocene correlations would be due to more time uncertainty and more climate variability in the LGM. We compare this simple prediction to the actual correlation structure in the LGM proxy records. We found that time series data of ice-core stable isotope records and planktonic foraminifera Mg/Ca ratios were consistent between the Holocene and LGM periods, while time series of Uk'37 proxy records were not as we found no correlation between nearby LGM records. Our results support the finding of highly regional or noisy marine proxy records in the compilation analysed here and suggest the need for further studies on the role of climate proxies and the processes of climate signal recording and preservation.
Groningen is the largest onshore gas field under production in Europe. The pressure depletion of the gas field started in 1963. In 1991, the first induced micro-earthquakes have been located at reservoir level with increasing rates in the following decades. Most of these events are of magnitude less than 2.0 and cannot be felt. However, maximum observed magnitudes continuously increased over the years until the largest, significant event with ML=3.6 was recorded in 2014, which finally led to the decision to reduce the production. This causal sequence displays the crucial role of understanding and modeling the relation between production and induced seismicity for economic planing and hazard assessment. Here we test whether the induced seismicity related to gas exploration can be modeled by the statistical response of fault networks with rate-and-state-dependent frictional behavior. We use the long and complete local seismic catalog and additionally detailed information on production-induced changes at the reservoir level to test different seismicity models. Both the changes of the fluid pressure and of the reservoir compaction are tested as input to approximate the Coulomb stress changes. We find that the rate-and-state model with a constant tectonic background seismicity rate can reproduce the observed long delay of the seismicity onset. In contrast, so-called Coulomb failure models with instantaneous earthquake nucleation need to assume that all faults are initially far from a critical state of stress to explain the delay. Our rate-and-state model based on the fluid pore pressure fits the spatiotemporal pattern of the seismicity best, where the fit further improves by taking the fault density and orientation into account. Despite its simplicity with only three free parameters, the rate-and-state model can reproduce the main statistical features of the observed activity.
The structural evolution of calderas in rifts helps to characterize the spatiotemporal relationships between magmatism, long wavelength crustal deformation and the formation of tectonic deformation zones along the rift axis. We document the structural characteristics of the c. 36 ka old Menengai Caldera located within a young zone of extension in the central Kenya Rift. Field mapping and high-resolution digital surface models show that NNE-striking Holocene normal faults perpendicular to the regional ESE-WNWextension direction dominate the interior sectors of the rift. Inside the caldera, these structures are overprinted by post-collapse doming and faulting of the magmatic centre, resulting in obliquely slipping normal faults bounding a resurgence horst. Radiocarbon dating of faulted units as young as 5 ka cal BP and the palaeo-shorelines of a lake formed during the African Humid Period in the Nakuru Basin indicate that volcanism and fault activity inside and in the vicinity of Menengai must have been sustained during the Holocene. Our analysis confirms that the caldera is located at the centre of an extending rift segment and suggests that other magmatic centres and young zones of faulting along the volcano-tectonic axis of the Kenya Rift may constitute nucleation points of faulting that ultimately foster future continental break-up.
The southern Central Andes (SCA) (between 27 degrees S and 40 degrees S) is bordered to the west by the convergent margin between the continental South American Plate and the oceanic Nazca Plate. The subduction angle along this margin is variable, as is the deformation of the upper plate. Between 33 degrees S and 35 degrees S, the subduction angle of the Nazca plate increases from sub-horizontal (< 5 degrees) in the north to relatively steep (similar to 30 degrees) in the south. The SCA contain inherited lithological and structural heterogeneities within the crust that have been reactivated and overprinted since the onset of subduction and associated Cenozoic deformation within the Andean orogen. The distribution of the deformation within the SCA has often been attributed to the variations in the subduction angle and the reactivation of these inherited heterogeneities. However, the possible influence that the thickness and composition of the continental crust have had on both short-term and long-term deformation of the SCA is yet to be thoroughly investigated. For our investigations, we have derived density distributions and thicknesses for various layers that make up the lithosphere and evaluated their relationships with tectonic events that occurred over the history of the Andean orogeny and, in particular, investigated the short- and long-term nature of the present-day deformation processes. We established a 3D model of lithosphere beneath the orogen and its foreland (29 degrees S-39 degrees S) that is consistent with currently available geological and geophysical data, including the gravity data. The modelled crustal configuration and density distribution reveal spatial relationships with different tectonic domains: the crystalline crust in the orogen (the magmatic arc and the main orogenic wedge) is thicker (similar to 55 km) and less dense (similar to 2900 kg/m(3)) than in the forearc (similar to 35 km, similar to 2975 kg/m(3)) and foreland (similar to 30 km, similar to 3000 kg/m(3)). Crustal thickening in the orogen probably occurred as a result of stacking of low-density domains, while density and thickness variations beneath the forearc and foreland most likely reflect differences in the tectonic evolution of each area following crustal accretion. No clear spatial relationship exists between the density distribution within the lithosphere and previously proposed boundaries of crustal terranes accreted during the early Paleozoic. Areas with ongoing deformation show a spatial correlation with those areas that have the highest topographic gradients and where there are abrupt changes in the average crustal-density contrast. This suggests that the short-term deformation within the interior of the Andean orogen and its foreland is fundamentally influenced by the crustal composition and the relative thickness of different crustal layers. A thicker, denser, and potentially stronger lithosphere beneath the northern part of the SCA foreland is interpreted to have favoured a strong coupling between the Nazca and South American plates, facilitating the development of a sub-horizontal slab.
The chapter illustrates how cities can be understood from a system–theory perspective as complex social systems. It argues that the classical and linear intervention methods are often no (longer) suitable for the complex structure, temporal dynamics and multifaceted processuality of urban development. It offers a systemic and systems theory-inspired method as an alternative approach, which allows for extended possibilities that are more appropriate for dealing with urban development processes. The method was developed on the basis of practical experience and theoretical insights. The approach should demonstrate for local decision-makers potential areas of activity for organising urban changes through co-design.
Permafrost is warming in the northern high latitudes, inducing highly dynamic thaw-related permafrost disturbances across the terrestrial Arctic. Monitoring and tracking of permafrost disturbances is important as they impact surrounding landscapes, ecosystems and infrastructure. Remote sensing provides the means to detect, map, and quantify these changes homogeneously across large regions and time scales. Existing Landsat-based algorithms assess different types of disturbances with similar spatiotemporal requirements. However, Landsat-based analyses are restricted in northern high latitudes due to the long repeat interval and frequent clouds, in particular at Arctic coastal sites. We therefore propose to combine Landsat and Sentinel-2 data for enhanced data coverage and present a combined annual mosaic workflow, expanding currently available algorithms, such as LandTrendr, to achieve more reliable time series analysis. We exemplary test the workflow for twelve sites across the northern high latitudes in Siberia. We assessed the number of images and cloud-free pixels, the spatial mosaic coverage and the mosaic quality with spectral comparisons. The number of available images increased steadily from 1999 to 2019 but especially from 2016 onward with the addition of Sentinel-2 images. Consequently, we have an increased number of cloud-free pixels even under challenging environmental conditions, which then serve as the input to the mosaicking process. In a comparison of annual mosaics, the Landsat+Sentinel-2 mosaics always fully covered the study areas (99.9–100 %), while Landsat-only mosaics contained data-gaps in the same years, only reaching coverage percentages of 27.2 %, 58.1 %, and 69.7 % for Sobo Sise, East Taymyr, and Kurungnakh in 2017, respectively. The spectral comparison of Landsat image, Sentinel-2 image, and Landsat+Sentinel-2 mosaic showed high correlation between the input images and mosaic bands (e.g., for Kurungnakh 0.91–0.97 between Landsat and Landsat+Sentinel-2 mosaic and 0.92–0.98 between Sentinel-2 and Landsat+Sentinel-2 mosaic) across all twelve study sites, testifying good quality mosaic results. Our results show that especially the results for northern, coastal areas was substantially improved with the Landsat+Sentinel-2 mosaics. By combining Landsat and Sentinel-2 data we accomplished to create reliably high spatial resolution input mosaics for time series analyses. Our approach allows to apply a high temporal continuous time series analysis to northern high latitude permafrost regions for the first time, overcoming substantial data gaps, and assess permafrost disturbance dynamics on an annual scale across large regions with algorithms such as LandTrendr by deriving the location, timing and progression of permafrost thaw disturbances
Two different approaches are used to assess the impacts associated with natural hazards and climate change in cities. A bottom-up approach uses high resolution data on constituent assets within the urban area. In contrast, a top-down approach uses less detailed information but is consequently more readily transferable. Here, we compare damage curves generated by each approach for coastal flooding in London. To compare them, we fit a log-logistic regression with three parameters to the calculated damage curves. We find that the functions are remarkably similar in their shape, albeit with different inflection points and a maximum damage that differs by 13%-25%. If rescaled, the curves agree almost exactly, which enables damage assessment to be undertaken following the calculation of the three parameters.
Controversy over the plate tectonic affinity and evolution of the Saxon granulites in a two- or multi-plate setting during inter- or intracontinental collision makes the Saxon Granulite Massif a key area for the understanding of the Palaeozoic Variscan orogeny. The massif is a large dome structure in which tectonic slivers of metapelite and metaophiolite units occur along a shear zone separating a diapir-like body of high-Pgranulite below from low-Pmetasedimentary rocks above. Each of the upper structural units records a different metamorphic evolution until its assembly with the exhuming granulite body. New age and petrologic data suggest that the metaophiolites developed from early Cambrian protoliths during high-Pamphibolite facies metamorphism in the mid- to late-Devonian and thermal overprinting by the exhuming hot granulite body in the early Carboniferous. A correlation of new Ar-Ar biotite ages with publishedP-T-tdata for the granulites implies that exhumation and cooling of the granulite body occurred at average rates of similar to 8 mm/year and similar to 80 degrees C/Ma, with a drop in exhumation rate from similar to 20 to similar to 2.5 mm/year and a slight rise in cooling rate between early and late stages of exhumation. A time lag ofc. 2 Ma between cooling through the closure temperatures for argon diffusion in hornblende and biotite indicates a cooling rate of 90 degrees C/Ma when all units had assembled into the massif. A two-plate model of the Variscan orogeny in which the above evolution is related to a short-lived intra-Gondwana subduction zone conflicts with the oceanic affinity of the metaophiolites and the timescale ofc. 50 Ma for the metamorphism. Alternative models focusing on the internal Variscan belt assume distinctly different material paths through the lower or upper crust for strikingly similar granulite massifs. An earlier proposed model of bilateral subduction below the internal Variscan belt may solve this problem.
Subduction of Neo-Tethys oceanic lithosphere beneath the Iranian plate during the Mesozoic formed several igneous bodies of ultramafic to intermediate and felsic composition. Intrusion of these magmas into a regional metamorphic sequence (the Sanandaj-Sirjan Zone) caused partial melting and formation of migmatites with meta-pelitic protoliths. The Alvand complex (west Iran) is a unique area comprising migmatites of both mafic and pelitic protoliths. In this area, the gabbroic rocks contain veins of leucosome at their contact with pyroxenite and olivine gabbro. These leucosomes are geochemically and mineralogically different from leucosomes of the meta-pelitic migmatites and clearly show properties of I-type granites. Microscopic observations and whole rock compositions of the mafic migmatite leucosomes show that migmatization occurred through partial melting of biotite, hornblende and plagioclase. Thermobarometric calculations indicate 800 degrees C and 3.7 kbar for partial melting, although phase diagram modeling demonstrates that the presence of water could decrease the solidus temperature by about 40 degrees C. Our results suggest an asthenospheric magma upwelling as the source of heat for partial melting of the gabbroic rock during subduction of Neo-Tethys oceanic crust under the western edge of the Iranian plate. The present study also reveals relationships between migmatization and formation of S- and I -type granites in the area.
Flood loss data collection and modeling are not standardized, and previous work has indicated that losses from different flood types (e.g., riverine and groundwater) may follow different driving forces. However, different flood types may occur within a single flood event, which is known as a compound flood event. Therefore, we aimed to identify statistical similarities between loss-driving factors across flood types and test whether the corresponding losses should be modeled separately. In this study, we used empirical data from 4,418 respondents from four survey campaigns studying households in Germany that experienced flooding. These surveys sought to investigate several features of the impact process (hazard, socioeconomic, preparedness, and building characteristics, as well as flood type). While the level of most of these features differed across flood type subsamples (e.g., degree of preparedness), they did so in a nonregular pattern. A variable selection process indicates that besides hazard and building characteristics, information on property-level preparedness was also selected as a relevant predictor of the loss ratio. These variables represent information, which is rarely adopted in loss modeling. Models shall be refined with further data collection and other statistical methods. To save costs, data collection efforts should be steered toward the most relevant predictors to enhance data availability and increase the statistical power of results. Understanding that losses from different flood types are driven by different factors is a crucial step toward targeted data collection and model development and will finally clarify conditions that allow us to transfer loss models in space and time. <br /> Key Points <br /> Survey data of flood-affected households show different concurrent flood types, undermining the use of a single-flood-type loss model Thirteen variables addressing flood hazard, the building, and property level preparedness are significant predictors of the building loss ratio Flood type-specific models show varying significance across the predictor variables, indicating a hindrance to model transferability
A silica undersaturated alkali-olivine basanitic magma intruded the late Paleocene/early Eocene Jafnayn Formation near Muscat. Geochemical analyses indicate that a significant amount of host rock (limestone) was assimilated into the magma. We dated the basanite as 42.7 +/- 1.0 Ma (2 sigma error; late Lutetian), using the whole rock Ar-40/Ar-39 step-wise heating technique. Intrusion occurred in the hanging wall of a major regional extensional shear zone (Frontal Range Fault, FRF) bounding the northern margin of two domes within the Oman Mountains (Jabal Akhdar and Saih Hatat domes). Two shear intervals along the FRF have been documented. The first interval lasted immediately after emplacement of the Semail Ophiolite (latest Cretaceous-early Eocene) while the second and poorly constrained interval was assumed to have occurred during the Oligocene.
The proximity of the basanite to the FRF suggests that magma used extensional faults for the upper part of its ascent path. Reactivated Permian rift faults of the Pangaea rift or other preexisting faults may have been used for the lower ascent part.
We conclude that the basanite intrusion coincided with the onset of the second deformation interval along the FRF, because (1) the position of the basanite is near a dextral releasing bend, associated with the second shear interval, (2) the overlap of our Ar-40/Ar-39 age with the cooling curves for rocks from the nearby Jabal Akhdar Dome, and (3) the basanite postdates the first FRF deformation episode by > 10 Ma. Thus, the second interval along the FRF had started already during the late Lutetian and probably lasted into the Miocene.
Drainage divide networks
(2020)
Drainage divides are organized into tree-like networks that may record information about drainage divide mobility. However, views diverge about how to best assess divide mobility. Here, we apply a new approach of automatically extracting and ordering drainage divide networks from digital elevation models to results from landscape evolution model experiments. We compared landscapes perturbed by strike-slip faulting and spatiotemporal variations in erodibility to a reference model to assess which topographic metrics (hillslope relief, flow distance, and chi) are diagnostic of divide mobility. Results show that divide segments that are a minimum distance of similar to 5 km from river confluences strive to attain constant values of hillslope relief and flow distance to the nearest stream. Disruptions of such patterns can be related to mobile divides that are lower than stable divides, closer to streams, and often asymmetric in shape. In general, we observe that drainage divides high up in the network, i.e., at great distances from river confluences, are more susceptible to disruptions than divides closer to these confluences and are thus more likely to record disturbance for a longer time period. We found that across-divide differences in hillslope relief proved more useful for assessing divide migration than other tested metrics. However, even stable drainage divide networks exhibit across-divide differences in any of the studied topographic metrics. Finally, we propose a new metric to quantify the connectivity of divide junctions.
Der zentralasiatische Naturraum, wie er sich uns heute präsentiert, ist das Ergebnis eines Zusammenwirkens vieler verschiedener Faktoren über Jahrmillionen hinweg. Im aktuellen Kontext des Klimawandels zeigt sich jedoch, wie stark sich Stoffflüsse auch kurzfristig ändern und dabei das Gesicht der Landschaft verwandeln können. Die Gobi-Wüste in der Inneren Mongolei (China), als Teil der gleichnamigen Trockenregionen Nordwestchinas, ist aufgrund der Ausgestaltung ihrer landschaftsprägenden Elemente sowie ihrer Landschaftsdynamik, im Zusammenhang mit der Lage zum Tibet-Plateau, in den Fokus der klimageschichtlichen Grundlagenforschung gerückt. Als großes Langzeitarchiv unterschiedlichster fluvialer, lakustriner und äolischer Sedimente stellt sie eine bedeutende Lokalität zur Rekonstruktion von lokalen und regionalen Stoffflüssen dar.. Andererseits ist die Gobi-Wüste zugleich auch eine bedeutende Quelle für den überregionalen Staubtransport, da sie aufgrund der klimatischen Bedingungen insbesondere der Erosion durch Ausblasung preisgegeben wird. Vor diesem Hintergrund erfolgten zwischen 2011 und 2014, im Rahmen des BMBF-Verbundprogramms WTZ Zentralasien – Monsundynamik & Geoökosysteme (Förderkennzeichen 03G0814), mehrere deutsch-chinesische Expeditionen in das Ejina-Becken (Innere Mongolei) und das Qilian Shan-Vorland. Im Zuge dieser Expeditionen wurden für eine Bestimmung potenzieller Sedimentquellen erstmals zahlreiche Oberflächenproben aus dem gesamten Einzugsgebiet des Heihe (schwarzer Fluss) gesammelt. Zudem wurden mit zwei Bohrungen im inneren des Ejina-Beckens, ergänzende Sedimentbohrkerne zum bestehenden Bohrkern D100 (siehe Wünnemann (2005)) abgeteuft, um weit reichende, ergänzende Informationen zur Landschaftsgeschichte und zum überregionalen Sedimenttransfer zu erhalten. Gegenstand und Ziel der vorliegenden Doktorarbeit ist die sedimentologisch-mineralogische Charakterisierung des Untersuchungsgebietes in Bezug auf potenzielle Sedimentquellen und Stoffflüsse des Ejina-Beckens sowie die Rekonstruktion der Ablagerungsgeschichte eines dort erbohrten, 19m langen Sedimentbohrkerns (GN100). Schwerpunkt ist hierbei die Klärung der Sedimentherkunft innerhalb des Bohrkerns sowie die Ausweisung von Herkunftssignalen und möglichen Sedimentquellen bzw. Sedimenttransportpfaden. Die methodische Herangehensweise basiert auf einem Multi-Proxy-Ansatz zur Charakterisierung der klastischen Sedimentfazies anhand von Geländebeobachtungen, lithologisch-granulometrischen und mineralogisch-geochemischen Analysen sowie statistischen Verfahren. Für die mineralogischen Untersuchungen der Sedimente wurde eine neue, rasterelektronenmikroskopische Methode zur automatisierten Partikelanalyse genutzt und den traditionellen Methoden gegenübergestellt. Die synoptische Betrachtung der granulometrischen, geochemischen und mineralogischen Befunde der Oberflächensedimente ergibt für das Untersuchungsgebiet ein logisches Kaskadenmodell mit immer wiederkehrenden Prozessbereichen und ähnlichen Prozesssignalen. Die umfangreichen granulometrischen Analysen deuten dabei auf abnehmende Korngrößen mit zunehmender Entfernung vom Qilian Shan hin und ermöglichen die Identifizierung von vier texturellen Signalen: den fluvialen Sanden, den Dünensanden, den Stillwassersedimenten und Stäuben. Diese Ergebnisse können als Interpretationsgrundlage für die Korngrößenanalysen des Bohrkerns genutzt werden. Somit ist es möglich, die Ablagerungsgeschichte der Bohrkernsedimente zu rekonstruieren und in Verbindung mit eigenen und literaturbasierten Datierungen in einen Gesamtkontext einzuhängen. Für das Untersuchungsgebiet werden somit vier Ablagerungsphasen ausgewiesen, die bis in die Zeit des letzten glazialen Maximums (LGM) zurückreichen. Während dieser Ablagerungsphasen kam es im Zuge unterschiedlicher Aktivitäts- und Stabilitätsphasen zu einer kontinuierlichen Progradation und Überprägung des Schwemmfächers. Eine besonders aktive Phase kann zwischen 8 ka und 4 ka BP festgestellt werden, während der es aufgrund zunehmender fluvialer Aktivitäten zu einer deutlich verstärkten Schwemmfächerdynamik gekommen zu sein scheint. In den Abschnitten davor und danach waren es vor allem äolische Prozesse, die zu einer Überprägung des Schwemmfächers geführt haben. Hinsichtlich der mineralogischen Herkunftssignale gibt es eine große Variabilität. Dies spiegelt die enorme Heterogenität der Geologie des Untersuchungsgebietes wider, wodurch die räumlichen Signale nicht sehr stark ausgeprägt sind. Dennoch, können für das Einzugsgebiet drei größere Bereiche deklariert werden, die als Herkunftsgebiet in Frage kommen. Das östliche Qilian Shan Vorland zeichnet sich dabei durch deutlich höhere Chloritgehalte als primäre Quelle für die Sedimente im Ejina-Becken aus. Sie unterscheiden sich insbesondere durch stark divergierende Chloritgehalte in der Tonmineral- und Gesamtmineralfraktion, was das östliche Qilian Shan Vorland als primäre Quelle für die Sedimente im Ejina-Becken auszeichnet. Dies steht in Zusammenhang mit den Grünschiefern, Ophioliten und Serpentiniten in diesem Bereich. Geochemisch deutet vor allem das Cr/Rb-Verhältnis eine große Variabilität innerhalb des Einzugsgebietes an. Auch hier ist es das östliche Vorland, welches aufgrund seines hohen Anteils an mafischen Gesteinen reich an Chromiten und Spinellen ist und sich somit vom restlichen Untersuchungsgebiet abhebt. Die zeitliche aber auch die generelle Variabilität der Sedimentherkunft lässt sich in den Bohrkernsedimenten nicht so deutlich nachzeichnen. Die mineralogisch-sedimentologischen Eigenschaften der erbohrten klastischen Sedimente zeugen zwar von zwischenzeitlichen Änderungen bei der Sedimentherkunft, diese sind jedoch nicht so deutlich ausgeprägt, wie es die Quellsignale in den Oberflächensedimenten vermuten lassen. Ein Grund dafür scheint die starke Vermischung unterschiedlichster Sedimente während des Transportes zu sein. Die Kombination der Korngrößenergebnisse mit den Befunden der Gesamt- und Schwermineralogie deuten darauf hin, dass es zwischenzeitlich eine Phase mit überwiegend äolischen Prozessen gegeben hat, die mit einem Sedimenteintrag aus dem westlichen Bei Shan in Verbindung stehen. Neben der Zunahme ultrastabiler Schwerminerale wie Zirkon und Granat und der Abnahme opaker Schwerminerale, weisen vor allem die heutigen Verhältnisse darauf hin. Der Vergleich der traditionellen Schwermineralanalyse mit der Computer-Controlled-Scanning-Electron-Microscopy (kurz: CCSEM), die eine automatisierte Partikelauswertung der Proben ermöglicht, zeigt den deutlichen Vorteil der modernen Analysemethode. Neben einem zeitlichen Vorteil, den man durch die automatisierte Abarbeitung der vorbereiteten Proben erlangen kann, steht vor allem die deutlich größere statistische Signifikanz des Ergebnisses im Vordergrund. Zudem können mit dieser Methode auch chemische Varietäten einiger Schwerminerale bestimmt werden, die eine noch feinere Klassifizierung und sicherere Aussagen zu einer möglichen Sedimentherkunft ermöglichen. Damit ergeben sich außerdem verbesserte Aussagen zu Zusammensetzungen und Entstehungsprozessen der abgelagerten Sedimente. Die Studie verdeutlicht, dass die Sedimentherkunft innerhalb des Untersuchungsgebietes sowie die ablaufenden Prozesse zum Teil stark von lokalen Gegebenheiten abhängen. Die Heterogenität der Geologie und die Größe des Einzugsgebietes sowie die daraus resultierende Komplexität der Sedimentgenese, machen exakte Zuordnungen zu klar definierten Sedimentquellen sehr schwer. Dennoch zeigen die Ergebnisse, dass die Sedimentzufuhr in das Ejina-Becken in erster Linie durch fluviale klastische Sedimente des Heihe aus dem Qilian Shan erfolgt sein muss. Die Untersuchungsergebnisse zeigen jedoch ebenso die Notwendigkeit einer ergänzenden Bearbeitung angrenzender Untersuchungsgebiete, wie beispielsweise den Gobi-Altai im Norden oder den Beishan im Westen, sowie die Verdichtung der Oberflächenbeprobung zur feineren Auflösung von lokalen Sedimentquellen.
Machine learning (ML) algorithms are being increasingly used in Earth and Environmental modeling studies owing to the ever-increasing availability of diverse data sets and computational resources as well as advancement in ML algorithms. Despite advances in their predictive accuracy, the usefulness of ML algorithms for inference remains elusive. In this study, we employ two popular ML algorithms, artificial neural networks and random forest, to analyze a large data set of flood events across Germany with the goals to analyze their predictive accuracy and their usability to provide insights to hydrologic system functioning. The results of the ML algorithms are contrasted against a parametric approach based on multiple linear regression. For analysis, we employ a model-agnostic framework named Permuted Feature Importance to derive the influence of models' predictors. This allows us to compare the results of different algorithms for the first time in the context of hydrology. Our main findings are that (1) the ML models achieve higher prediction accuracy than linear regression, (2) the results reflect basic hydrological principles, but (3) further inference is hindered by the heterogeneity of results across algorithms. Thus, we conclude that the problem of equifinality as known from classical hydrological modeling also exists for ML and severely hampers its potential for inference. To account for the observed problems, we propose that when employing ML for inference, this should be made by using multiple algorithms and multiple methods, of which the latter should be embedded in a cross-validation routine.
Flood loss modeling is a central component of flood risk analysis. Conventionally, this involves univariable and deterministic stage-damage functions. Recent advancements in the field promote the use of multivariable and probabilistic loss models, which consider variables beyond inundation depth and account for prediction uncertainty. Although companies contribute significantly to total loss figures, novel modeling approaches for companies are lacking. Scarce data and the heterogeneity among companies impede the development of company flood loss models. We present three multivariable flood loss models for companies from the manufacturing, commercial, financial, and service sector that intrinsically quantify prediction uncertainty. Based on object-level loss data (n = 1,306), we comparatively evaluate the predictive capacity of Bayesian networks, Bayesian regression, and random forest in relation to deterministic and probabilistic stage-damage functions, serving as benchmarks. The company loss data stem from four postevent surveys in Germany between 2002 and 2013 and include information on flood intensity, company characteristics, emergency response, private precaution, and resulting loss to building, equipment, and goods and stock. We find that the multivariable probabilistic models successfully identify and reproduce essential relationships of flood damage processes in the data. The assessment of model skill focuses on the precision of the probabilistic predictions and reveals that the candidate models outperform the stage-damage functions, while differences among the proposed models are negligible. Although the combination of multivariable and probabilistic loss estimation improves predictive accuracy over the entire data set, wide predictive distributions stress the necessity for the quantification of uncertainty.
New Zealand's Alpine Fault is a large, platebounding strike-slip fault, which ruptures in large (M-w > 8) earthquakes. We conducted field and laboratory analyses of fault rocks to assess its fault zone architecture. Results reveal that the Alpine Fault Zone has a complex geometry, comprising an anastomosing network of multiple slip planes that have accommodated different amounts of displacement. This contrasts with the previous perception of the Alpine Fault Zone, which assumes a single principal slip zone accommodated all displacement. This interpretation is supported by results of drilling projects and geophysical investigations. Furthermore, observations presented here show that the young, largely unconsolidated sediments that constitute the footwall at shallow depths have a significant influence on fault gouge rheological properties and structure.
Many of the volcanic plateau margins of the eastern, formerly glaciated, foreland of the Patagonian Andes are undermined by giant landslides (>= 10(8) m(3)). One cluster of such landslides extends along the margin of the Meseta del Lago Buenos Aires (MLBA) plateau that is formed mainly by Neogene-Quaternary basalts. The dry climate is at odds with numerous >2-km long earthflows nested within older and larger compound landslides. We present a hydrological analysis, a detailed geomorphic map, interpretations of exposed landslide interiors, and radiocarbon dating of the El Mirador landslide, which is one of the largest and morphologically most representative landslide. We find that the presence of lakes on top of the plateau, causing low infiltration rates, correlates negatively with the abundance of earthflows on compound landslides along the plateau margins. Field outcrops show that the pattern of compound landslides and earthflows is likely controlled by groundwater seepage at the contact between the basalts and underlying soft Miocene molasse. Numerous peat bogs store water and sediment and are more abundant in earthflow-affected areas than in their contributing catchment areas. <br /> Radiocarbon dates indicate that these earthflows displaced metre-thick layers of peat in the late Holocene (<2.5 ka). We conclude that earthflows of the MLBA plateau might be promising proxies of past hydroclimatic conditions in the Patagonian foreland, if strong earthquakes or gradual crustal stress changes due to glacioisostatic rebound can be ruled out.
The Coulomb failure stress (CFS) criterion is the most commonly used method for predicting spatial distributions of aftershocks following large earthquakes. However, large uncertainties are always associated with the calculation of Coulomb stress change. The uncertainties mainly arise due to nonunique slip inversions and unknown receiver faults; especially for the latter, results are highly dependent on the choice of the assumed receiver mechanism. Based on binary tests (aftershocks yes/no), recent studies suggest that alternative stress quantities, a distance-slip probabilistic model as well as deep neural network (DNN) approaches, all are superior to CFS with predefined receiver mechanism. To challenge this conclusion, which might have large implications, we use 289 slip inversions from SRCMOD database to calculate more realistic CFS values for a layered half-space and variable receiver mechanisms. We also analyze the effect of the magnitude cutoff, grid size variation, and aftershock duration to verify the use of receiver operating characteristic (ROC) analysis for the ranking of stress metrics. The observations suggest that introducing a layered half-space does not improve the stress maps and ROC curves. However, results significantly improve for larger aftershocks and shorter time periods but without changing the ranking. We also go beyond binary testing and apply alternative statistics to test the ability to estimate aftershock numbers, which confirm that simple stress metrics perform better than the classic Coulomb failure stress calculations and are also better than the distance-slip probabilistic model.
The 2.7-2.9 Ma Ertsberg East Skarn System (EESS) is a world-class Cu-Au skarn that formed within and adjacent to an intrusion within a paleodepth of 0.5 km and > 2.5 km. Its economic mineralisation developed by sustained reaction of magmatic fluid with contact metamorphosed siliciclastic and carbonate rocks at the margin of the adjacent Ertsberg quartz monzodiorite intrusion. Based on high-resolution mineral mapping, chemical analysis and thermodynamic calculations, the multistage formation processes of the exoskarn components of the EESS are examined in the context of changing pressure, temperature, fluid composition and fluid phase. We show that contact metamorphism of dolomitic sediments occurred at 51 +/- 5 MPa, between 700 degrees C and 800 degrees C and in the presence of a H2O-CO2-fluid containing similar to 10 to similar to 70 mol% CO2. This prograde metamorphism formed a forsterite + diopside + calcite + phlogopite + spinel assemblage. Such forsterite-dominated skarns account for similar to 55 vol% of the EESS exoskarns. Rare pargasite (previously unrecognized in this deposit) formed locally in the metamorphosed carbonate sequence where the protolith was composed of supratidal evaporites with dolomitic carbonate and interlayered calc-silicate rocks. The subsequent flux of a lower pressure magmatic gas containing SO2(g) caused sulphate metasomatism. This high temperature gas alteration of the metamorphic assemblage also caused skarn Cu-Fe-sulphide mineralisation. The influx of a SO2 gas through fracture permeability occurred at a temperature between similar to 600 and 700 degrees C and caused calcite to be replaced by anhydrite, with the coupled release of H2S(g). This in-situ release of H2S(g) scavenged trace Cu from the gas phase to deposit Cu-Fe-sulphides, which make the economic value of the distinct. We demonstrate that the formation of metal sulphides within forsterite skarns of the Ertsberg East Skarn System required a minimum flux of similar to 1,050 Mt SO2(g) and show that volcanic degassing may have occurred over a time span of similar to 3,900 years. As the system waned, the ambient fluid resulted in partial retrograde serpentinization of olivine and diopside without carbonation, and at temperatures sufficiently high to preserve anhydrite.
High-pressure experiments were performed to investigate the effectiveness, rate and mechanism of carbonation of serpentinites by a carbon-saturated COH fluid at 1.5-2.5 GPa and 375-700 degrees C. This allows a better understanding of the fate and redistribution of slab-derived carbonic fluids when they react with the partially hydrated mantle within and above the subducting slab under pressure and temperature conditions corresponding to the forearc mantle. Interactions between carbon-saturated CO2-H2O-CH4 fluids and serpentinite were investigated using natural serpentinite cylinders with natural grain sizes and shapes in piston-cylinder experiments. The volatile composition of post-run fluids was quantified by gas chromatography. Solid phases were examined by Raman spectroscopy, electron microscopy and laser ablation inductively coupled plasma mass spectrometry. Textures, porosity and phase abundances of recovered rock cores were visualized and quantified by three-dimensional, high-resolution computed tomography. We find that carbonation of serpentinites is efficient at sequestering CO2 from the interacting fluid into newly formed magnesite. Time-series experiments demonstrate that carbonation is completed within similar to 96 h at 2 GPa and 600 degrees C. With decreasing CO2, aq antigorite is replaced first by magnesite + quartz followed by magnesite + talc + chlorite in distinct, metasomatic fronts. Above antigorite stability magnesite + enstatite + talc + chlorite occur additionally. The formation of fluid-permeable reaction zones enhances the reaction rate and efficiency of carbonation. Carbonation probably occurs via an interface-coupled replacement process, whereby interconnected porosity is present within reaction zones after the experiment. Consequently, carbonation of serpentinites is self-promoting and efficient even if fluid flow is channelized into veins. We conclude that significant amounts of carbonates may accumulate, over time, in the hydrated forearc mantle.
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.
High Mountain Asia (HMA) is dependent upon both the amount and timing of snow and glacier meltwater. Previous model studies and coarse resolution (0.25° × 0.25°, ∼25 km × 25 km) passive microwave assessments of trends in the volume and timing of snowfall, snowmelt, and glacier melt in HMA have identified key spatial and seasonal heterogeneities in the response of snow to changes in regional climate. Here we use recently developed, continuous, internally consistent, and high-resolution passive microwave data (3.125 km × 3.125 km, 1987–2016) from the special sensor microwave imager instrument family to refine and extend previous estimates of changes in the snow regime of HMA. We find an overall decline in snow volume across HMA; however, there exist spatially contiguous regions of increasing snow volume—particularly during the winter season in the Pamir, Karakoram, Hindu Kush, and Kunlun Shan. Detailed analysis of changes in snow-volume trends through time reveal a large step change from negative trends during the period 1987–1997, to much more positive trends across large regions of HMA during the periods 1997–2007 and 2007–2016. We also find that changes in high percentile monthly snow-water volume exhibit steeper trends than changes in low percentile snow-water volume, which suggests a reduction in the frequency of high snow-water volumes in much of HMA. Regions with positive snow-water storage trends generally correspond to regions of positive glacier mass balances.
The magnetic equator in the Brazilian region has moved over 1,100 km northward since 1957, passing the geomagnetic observatory Tatuoca (TTB), in northern Brazil, around 2013. We recovered and processed TTB hourly mean values of the geomagnetic field horizontal (H) component from 1957 until 2019, allowing the investigation of long-term changes in the daily variation due to the influence of secular variation, solar activity, season, and lunar phase. The H day-to-day variability and the occurrence of the counter electrojet at TTB were also investigated. Until the 1990s, ionospheric solar quiet currents dominated the quiet-time daily variation at TTB. After 2000, the magnitude of the daily variation became appreciably greater due to the equatorial electrojet (EEJ) contribution. The H seasonal and day-to-day variability increased as the magnetic equator approached, but their amplitudes normalized to the average daily variation remained at similar levels. Meanwhile, the amplitude of the lunar variation, normalized in the same way, increased from 5% to 12%. Within the EEJ region, the occurrence rate of the morning counter electrojet (MCEJ) increased with proximity to the magnetic equator, while the afternoon counter electrojet (ACEJ) did not. EEJ currents derived from CHAMP and Swarm satellite data revealed that the MCEJ rate varies with magnetic latitude within the EEJ region while the ACEJ rate is largely constant. Simulations with the Thermosphere-Ionosphere-Electrodynamics General Circulation Model based on different geomagnetic main field configurations suggest that long-term changes in the geomagnetic daily variation at TTB can be attributed to the main field secular variation.
Accelerograms are the primary source for characterizing strong ground motion. It is therefore of paramount interest to have high-quality recordings free from any nonphysical contamination. Frequently, accelerograms are affected by baseline jumps and drifts, either related to the instrument and/or a major earthquake. In this work, I propose a correction method for these undesired baseline drifts based on segmented linear least squares. The algorithm operates on the integrated waveforms and combines all three instrument components to estimate a model that modifies the baseline to be at zero continuously. The procedure consists of two steps: first a suite of models with variable numbers of discontinuities is derived for all three instrument components. During this process, the number of discontinuities is reduced in a parsimonious way, for example, two very close discontinuities are merged into a single one. In the second step, the optimal model is selected on the basis of the Bayesian information criterion. I exemplify the application on synthetic waveforms with known discontinuities and on observed waveforms from a unified strong-motion database of the Japan Meteorological Agency (JMA) and the National Research Institute for Earth Science and Disaster Prevention (NIED, Japan) networks for the major events of the 2016 Kumamoto earthquakes. After the baseline jump correction, the waveforms are furthermore corrected for displacement according to Wang et al.(2011). The resulting displacements are comparable to the Interferometric Synthetic Aperture Radar-derived displacement estimates for the Kumamoto earthquake sequence.
Earthquakes often rupture across more than one fault segment. If such rupture segmentation occurs on a significant scale, a simple point-source or one-fault model may not represent the rupture process well. As a consequence earthquake characteristics inferred, based on one-source assumptions, may become systematically wrong. This might have effects on follow-up analyses, for example regional stress field inversions and seismic hazard assessments. While rupture segmentation is evident for most M-w > 7 earthquakes, also smaller ones with 5.5 < M-w < 7 can be segmented. We investigate the sensitivity of globally available data sets to rupture segmentation and their resolution to reliably estimate the mechanisms in presence of segmentation. We focus on the sensitivity of InSAR (Interferometric Synthetic Aperture Radar) data in the static near-field and seismic waveforms in the far-field of the rupture and carry out non-linear and Bayesian optimizations of single-source and two-sources kinematic models (double-couple point sources and finite, rectangular sources) using InSAR and teleseismic waveforms separately. Our case studies comprises of four M-w 6-7 earthquakes: the 2009 L'Aquila and 2016 Amatrice (Italy) and the 2005 and 2008 Zhongba (Tibet) earthquakes. We contrast the data misfits of different source complexity by using the Akaike informational criterion (AIC). We find that the AIC method is well suited for data-driven inferences on significant rupture segmentation for the given data sets. This is based on our observation that an AIC-stated significant improvement of data fit for two-segment models over one-segment models correlates with significantly different mechanisms of the two source segments and their average compared to the single-segment mechanism. We attribute these modelled differences to a sufficient sensitivity of the data to resolve rupture segmentation. Our results show that near-field data are generally more sensitive to rupture segmentation of shallow earthquakes than far-field data but that also teleseismic data can resolve rupture segmentation in the studied magnitude range. We further conclude that a significant difference in the modelled source mechanisms for different segmentations shows that an appropriate choice of model segmentation matters for a robust estimation of source mechanisms. It reduces systematic biases and trade-off and thereby improves the knowledge on the rupture. Our study presents a strategy and method to detect significant rupture segmentation such that an appropriate model complexity can be used in the source mechanism inference. A similar, systematic investigation of earthquakes in the range of M-w 5.5-7 could provide important hazard-relevant statistics on rupture segmentation. In these cases single-source models introduce a systematic bias. Consideration of rupture segmentation therefore matters for a robust estimation of source mechanisms of the studied earthquakes.
AM(W)= 5.1 earthquake on January 21st, 2016 marked the beginning of a significant seismic sequence in the southern Alboran Sea, culminating in aM(W)= 6.3 earthquake on January 25th, and continuing with further moderate magnitude earthquakes until March. We use data from 35 seismic broadband stations in Spain, Morocco and Portugal to relocate the seismicity, estimate seismic moment tensors, and isolate regional apparent source time functions for the main earthquake. Relocation and regional moment tensor inversion consistently yield very shallow depths for the majority of events. We obtain 50 moment tensors for the sequence, showing a mixture of strike-slip faulting for the foreshock and the main event and reverse faulting for the major aftershocks. The leading role of reverse focal mechanisms among the aftershocks may be explained by the geometry of the fault network. The mainshock nucleates at a bend along the left-lateral Al-Idrisi fault, introducing local transpression within the transtensional Alboran Basin. The shallow depths of the 2016 Alboran Sea earthquakes may favor slip-partitioning on the involved faults. Apparent source durations for the main event suggest a similar to 21 km long, asymmetric rupture that propagates primarily toward NE into the restraining fault segment, with fast rupture speed of similar to 3.0 km/s. Consistently, the inversion for laterally variable fault displacement situates the main slip in the restraining segment. The partitioning into strike-slip rupture and dip-slip aftershocks confirms a non-optimal orientation of this segment, and suggests that the 2016 event settled a slip deficit from previous ruptures that could not propagate into the stronger restraining segment.
Filling in for Governments?
(2020)
The 2015 Paris Agreement on climate change abandons the Kyoto Protocol's paradigm of binding emissions targets and relies instead on countries' voluntary contributions. However, the Paris Agreement encourages not only governments but also sub-national governments, corporations and civil society to contribute to reaching ambitious climate goals. In a transition from the regulated architecture of the Kyoto Protocol to the open system of the Paris Agreement, the Agreement seeks to integrate non-state actors into the treaty-based climate regime. In 2014 the secretariat of the United Nations Framework Convention on Climate Change Peru and France created the Non-State Actor Zone for Climate Action (and launched the Global Climate Action portal). In December 2019, this portal recorded more than twenty thousand climate-commitments of private and public non-state entities, making the non-state venues of international climate meetings decisively more exciting than the formal negotiation space. This level engagement and governments' response to it raises a flurry of questions in relation to the evolving nature of the climate regime and climate change governance, including the role of private actors as standard setters and the lack of accountability mechanisms for non-state actions. This paper takes these developments as occasion to discuss the changing role of private actors in the climate regime.
The origin of Asian monsoons
(2020)
The Cenozoic inception and development of the Asian monsoons remain unclear and have generated much debate, as several hypotheses regarding circulation patterns at work in Asia during the Eocene have been proposed in the few last decades. These include (a) the existence of modern-like monsoons since the early Eocene; (b) that of a weak South Asian monsoon (SAM) and little to no East Asian monsoon (EAM); or (c) a prevalence of the Intertropical Convergence Zone (ITCZ) migrations, also referred to as Indonesian-Australian monsoon (I-AM). As SAM and EAM are supposed to have been triggered or enhanced primarily by Asian palaeogeographic changes, their possible inception in the very dynamic Eocene palaeogeographic context remains an open question, both in the modelling and field-based communities. We investigate here Eocene Asian climate conditions using the IPSL-CM5A2 (Sepulchre et al., 2019) earth system model and revised palaeogeographies. Our Eocene climate simulation yields atmospheric circulation patterns in Asia substantially different from modern conditions. A large high-pressure area is simulated over the Tethys ocean, which generates intense low tropospheric winds blowing southward along the western flank of the proto-Himalayan-Tibetan plateau (HTP) system. This low-level wind system blocks, to latitudes lower than 10 degrees N, the migration of humid and warm air masses coming from the Indian Ocean. This strongly contrasts with the modern SAM, during which equatorial air masses reach a latitude of 20-25 degrees N over India and southeastern China. Another specific feature of our Eocene simulation is the widespread subsidence taking place over northern India in the midtroposphere (around 5000 m), preventing deep convective updraught that would transport water vapour up to the condensation level. Both processes lead to the onset of a broad arid region located over northern India and over the HTP. More humid regions of high seasonality in precipitation encircle this arid area, due to the prevalence of the Intertropical Convergence Zone (ITCZ) migrations (or Indonesian-Australian monsoon, I-AM) rather than monsoons. Although the existence of this central arid region may partly result from the specifics of our simulation (model dependence and palaeogeographic uncertainties) and has yet to be confirmed by proxy records, most of the observational evidence for Eocene monsoons are located in the highly seasonal transition zone between the arid area and the more humid surroundings. We thus suggest that a zonal arid climate prevailed over Asia before the initiation of monsoons that most likely occurred following Eocene palaeogeographic changes. Our results also show that precipitation seasonality should be used with caution to infer the presence of a monsoonal circulation and that the collection of new data in this arid area is of paramount importance to allow the debate to move forward.
Sediment archives in the terrestrial and marine realm are regularly analyzed to infer changes in climate, tectonic, or anthropogenic boundary conditions of the past. However, contradictory observations have been made regarding whether short period events are faithfully preserved in stratigraphic archives; for instance, in marine sediments offshore large river systems. On the one hand, short period events are hypothesized to be non-detectable in the signature of terrestrially derived sediments due to buffering during sediment transport along large river systems. On the other hand, several studies have detected signals of short period events in marine records offshore large river systems. We propose that this apparent discrepancy is related to the lack of a differentiation between different types of signals and the lack of distinction between river response times and signal propagation times. In this review, we (1) expand the definition of the term ‘signal’ and group signals in sub-categories related to hydraulic grain size characteristics, (2) clarify the different types of ‘times’ and suggest a precise and consistent terminology for future use, and (3) compile and discuss factors influencing the times of signal transfer along sediment routing systems and how those times vary with hydraulic grain size characteristics. Unraveling different types of signals and distinctive time periods related to signal propagation addresses the discrepancies mentioned above and allows a more comprehensive exploration of event preservation in stratigraphy – a prerequisite for reliable environmental reconstructions from terrestrially derived sedimentary records.
Barite scales in geothermal installations are a highly unwanted effect of circulating deep saline fluids. They build up in the reservoir if supersaturated fluids are re-injected, leading to irreversible loss of injectivity. A model is presented for calculating the total expected barite precipitation. To determine the related injectivity decline over time, the spatial precipitation distribution in the subsurface near the injection well is assessed by modelling barite growth kinetics in a radially diverging Darcy flow domain. Flow and reservoir properties as well as fluid chemistry are chosen to represent reservoirs subject to geothermal exploration located in the North German Basin (NGB) and the Upper Rhine Graben (URG) in Germany. Fluids encountered at similar depths are hotter in the URG, while they are more saline in the NGB. The associated scaling amount normalised to flow rate is similar for both regions. The predicted injectivity decline after 10 years, on the other hand, is far greater for the NGB (64%) compared to the URG (24%), due to the temperature- and salinity-dependent precipitation rate. The systems in the NGB are at higher risk. Finally, a lightweight score is developed for approximating the injectivity loss using the Damkohler number, flow rate and total barite scaling potential. This formula can be easily applied to geothermal installations without running complex reactive transport simulations.
In near- surface geophysics, ground-based mapping surveys are routinely used in a variety of applications including those from archaeology, civil engineering, hydrology, and soil science. The resulting geophysical anomaly maps of, for example, magnetic or electrical parameters are usually interpreted to laterally delineate subsurface structures such as those related to the remains of past human activities, subsurface utilities and other installations, hydrological properties, or different soil types. To ease the interpretation of such data sets, we have developed a multiscale processing, analysis, and visualization strategy. Our approach relies on a discrete redundant wavelet transform (RWT) implemented using cubic-spline filters and the a trous algorithm, which allows to efficiently compute a multiscale decomposition of 2D data using a series of 1D convolutions. The basic idea of the approach is presented using a synthetic test image, whereas our archaeogeophysical case study from northeast Germany demonstrates its potential to analyze and process rather typical geophysical anomaly maps including magnetic and topographic data. Our vertical-gradient magnetic data show amplitude variations over several orders of magnitude, complex anomaly patterns at various spatial scales, and typical noise patterns, whereas our topographic data show a distinct hill structure superimposed by a microtopographic stripe pattern and random noise. Our results demonstrate that the RWT approach is capable to successfully separate these components and that selected wavelet planes can be scaled and combined so that the reconstructed images allow for a detailed, multiscale structural interpretation also using integrated visualizations of magnetic and topographic data. Because our analysis approach is straightforward to implement without laborious parameter testing and tuning, computationally efficient, and easily adaptable to other geophysical data sets, we believe that it can help to rapidly analyze and interpret different geophysical mapping data collected to address a variety of near-surface applications from engineering practice and research.
Neodymium isotopic composition (epsilon Nd) has enjoyed widespread use as a palaeotracer, principally because it behaves quasi-conservatively in the modern ocean. However, recent bottom water epsilon Nd reconstructions from the eastern North Atlantic are difficult to interpret under assumptions of conservative behaviour. The observation that this apparent departure from conservative behaviour increases with enhanced ice-rafted debris (IRD) fluxes has resulted in the suggestion that IRD leads to the overprinting of bottom water epsilon Nd through reversible scavenging. In this study, a simple water column model successfully reproduces epsilon Nd reconstructions from the eastern North Atlantic at the Last Glacial Maximum and Heinrich Stadial 1, and demonstrates that the changes in scavenging intensity required for good model-data fit is in good agreement with changes in the observed IRD flux. Although uncertainties in model parameters preclude a more definitive conclusion, the results indicate that the suggestion of IRD as a source of non-conservative behaviour in the epsilon Nd tracer is reasonable and that further research into the fundamental chemistry underlying the marine neodymium cycle is necessary to increase confidence in assumptions of conservative epsilon Nd behaviour in the past.
Multiple P-T-d-t paths reveal the evolution of the final Nuna assembly in northeast Australia
(2020)
The final assembly of the Mesoproterozoic supercontinent Nuna was marked by the collision of Laurentia and Australia at 1.60 Ga, which is recorded in the Georgetown Inlier of NE Australia. Here, we decipher the metamorphic evolution of this final Nuna collisional event using petrostructural analysis, major and trace element compositions of key minerals, thermodynamic modelling, and multi-method geochronology. The Georgetown Inlier is characterised by deformed and metamorphosed 1.70-1.62 Ga sedimentary and mafic rocks, which were intruded byc. 1.56 Ga old S-type granites. Garnet Lu-Hf and monazite U-Pb isotopic analyses distinguish two major metamorphic events (M1 atc. 1.60 Ga and M2 atc. 1.55 Ga), which allows at least two composite fabrics to be identified at the regional scale-c. 1.60 Ga S1 (consisting in fabrics S1a and S1b) andc. 1.55 Ga S2 (including fabrics S2a and S2b). Also, three tectono-metamorphic domains are distinguished: (a) the western domain, with S1 defined by low-P(LP) greenschist facies assemblages; (b) the central domain, where S1 fabric is preserved as medium-P(MP) amphibolite facies relicts, and locally as inclusion trails in garnet wrapped by the regionally dominant low-Pamphibolite facies S2 fabric; and (c) the eastern domain dominated by upper amphibolite to granulite facies S2 foliation. In the central domain, 1.60 GaMP-medium-T(MT) metamorphism (M1) developed within the staurolite-garnet stability field, with conditions ranging from 530-550 degrees C at 6-7 kbar (garnet cores) to 620-650 degrees C at 8-9 kbar (garnet rims), and it is associated with S1 fabric. The onset of 1.55 GaLP-high-T(HT) metamorphism (M2) is marked by replacement of staurolite by andalusite (M2a/D2a), which was subsequently pseudomorphed by sillimanite (M2b/D2b) where granite and migmatite are abundant.P-Tconditions ranged from 600 to 680 degrees C and 4-6 kbar for the M2b sillimanite stage. 1.60 Ga garnet relicts within the S2 foliation highlight the progressive obliteration of the S1 fabric by regional S2 in the central zone during peak M2 metamorphism. In the eastern migmatitic complex, partial melting of paragneiss and amphibolite occurred syn- to post-S2, at 730-770 degrees C and 6-8 kbar, and at 750-790 degrees C and 6 kbar, respectively. The pressure-temperature-deformation-time paths reconstructed for the Georgetown Inlier suggest ac. 1.60 Ga M1/D1 event recorded under greenschist facies conditions in the western domain and under medium-Pand medium-Tconditions in the central domain. This event was followed by the regional 1.56-1.54 Ga low-Pand high-Tphase (M2/D2), extensively recorded in the central and eastern domains. Decompression between these two metamorphic events is ascribed to an episode of exhumation. The two-stage evolution supports the previous hypothesis that the Georgetown Inlier preserves continental collisional and subsequent thermal perturbation associated with granite emplacement.
The steady increase of ground-motion data not only allows new possibilities but also comes with new challenges in the development of ground-motion models (GMMs). Data classification techniques (e.g., cluster analysis) do not only produce deterministic classifications but also probabilistic classifications (e.g., probabilities for each datum to belong to a given class or cluster). One challenge is the integration of such continuous classification in regressions for GMM development such as the widely used mixed-effects model. We address this issue by introducing an extension of the mixed-effects model to incorporate data weighting. The parameter estimation of the mixed-effects model, that is, fixed-effects coefficients of the GMMs and the random-effects variances, are based on the weighted likelihood function, which also provides analytic uncertainty estimates. The data weighting permits for earthquake classification beyond the classical, expert-driven, binary classification based, for example, on event depth, distance to trench, style of faulting, and fault dip angle. We apply Angular Classification with Expectation-maximization, an algorithm to identify clusters of nodal planes from focal mechanisms to differentiate between, for example, interface- and intraslab-type events. Classification is continuous, that is, no event belongs completely to one class, which is taken into account in the ground-motion modeling. The theoretical framework described in this article allows for a fully automatic calibration of ground-motion models using large databases with automated classification and processing of earthquake and ground-motion data. As an example, we developed a GMM on the basis of the GMM by Montalva et al. (2017) with data from the strong-motion flat file of Bastias and Montalva (2016) with similar to 2400 records from 319 events in the Chilean subduction zone. Our GMM with the data-driven classification is comparable to the expert-classification-based model. Furthermore, the model shows temporal variations of the between-event residuals before and after large earthquakes in the region.
Variation of deuterium excess in surface waters across a 5000-m elevation gradient in eastern Nepal
(2020)
The strong elevation gradient of the Himalaya allows for investigation of altitude and orographic impacts on surface water delta O-18 and delta D stable isotope values. This study differentiates the time- and altitude-variable contributions of source waters to the Arun River in eastern Nepal. It provides isotope data along a 5000-m gradient collected from tributaries as well as groundwater, snow, and glacial-sourced surface waters and time-series data from April to October 2016. We find nonlinear trends in delta O-18 and delta D lapse rates with high-elevation lapse rates (4000-6000 masl) 5-7 times more negative than low-elevation lapse rates (1000-3000 masl). A distinct seasonal signal in delta O-18 and delta D lapse rates indicates time-variable source-water contributions from glacial and snow meltwater as well as precipitation transitions between the Indian Summer Monsoon and Winter Westerly Disturbances. Deuterium excess correlates with the extent of snowpack and tracks melt events during the Indian Summer Monsoon season. Our analysis identifies the influence of snow and glacial melt waters on river composition during low-flow conditions before the monsoon (April/May 2016) followed by a 5-week transition to the Indian Summer Monsoon-sourced rainfall around mid-June 2016. In the post-monsoon season, we find continued influence from glacial melt waters as well as ISM-sourced groundwater.
The change of the mechanical properties of granular materials with pressure is an important topic associated with many industrial applications. In this paper we investigate the influence of hydrostatic pressure (P-e) on the effective bulk compressibility (C-eff) of a granular material by applying two modified theoretical approaches that are based on contact mechanics and micromechanics, respectively. For a granular material composed of rough grains, an extended contact model is developed to elucidate the effect of roughness of grain surfaces on bulk compressibility. At relatively low pressures, the model predicts that the decrease of bulk compressibility with pressure may be described by a power law with an exponent of -1/2 (i.e., C-eff proportional to P-e(1/2) ), but deviates at intermediate pressures. At elevated pressures beyond full contact, bulk compressibility remains almost unchanged, which may be roughly evaluated by continuum contact mechanics. As an alternative explanation of pressure-dependent bulk compressibility, we suggest a micromechanical model that accounts for effects of different types of pore space present in granular materials. Narrow and compliant inter-granular cracks are approximated by three-dimensional oblate spheroidal cracks with rough surfaces, whereas the equant and stiff pores surrounded by three and four neighboring grains are modeled as tubular pores with cross sections of three and four cusp-like corners, respectively. In this model, bulk compressibility is strongly reduced with increasing pressure by progressive closure of rough-walled cracks. At pressures exceeding crack closure pressure, deformation of the remaining equant pores is largely insensitive to pressure, with almost no further change in bulk compressibility. To validate these models, we performed hydrostatic compression tests on Bentheim sandstone (a granular rock consisting of quartz with high porosity) under a wide range of pressure. The relation between observed microstructures and measured pressure-dependent bulk compressibility is well explained by both suggested models.
Understanding the relation between injection-induced seismic moment release and operational parameters is crucial for early identification of possible seismic hazards associated with fluid-injection projects. We conducted laboratory fluid-injection experiments on permeable sandstone samples containing a critically stressed fault at different fluid pressurization rates. The observed fluid-induced fault deformation is dominantly aseismic. Fluid-induced stick-slip and fault creep reveal that total seismic moment release of acoustic emission (AE) events is related to total injected volume, independent of respective fault slip behavior. Seismic moment release rate of AE scales with measured fault slip velocity. For injection-induced fault slip in a homogeneous pressurized region, released moment shows a linear scaling with injected volume for stable slip (steady slip and fault creep), while we find a cubic relation for dynamic slip. Our results highlight that monitoring evolution of seismic moment release with injected volume in some cases may assist in discriminating between stable slip and unstable runaway ruptures.
A ground motion logic tree for seismic hazard analysis in the stable cratonic region of Europe
(2020)
Regions of low seismicity present a particular challenge for probabilistic seismic hazard analysis when identifying suitable ground motion models (GMMs) and quantifying their epistemic uncertainty. The 2020 European Seismic Hazard Model adopts a scaled backbone approach to characterise this uncertainty for shallow seismicity in Europe, incorporating region-to-region source and attenuation variability based on European strong motion data. This approach, however, may not be suited to stable cratonic region of northeastern Europe (encompassing Finland, Sweden and the Baltic countries), where exploration of various global geophysical datasets reveals that its crustal properties are distinctly different from the rest of Europe, and are instead more closely represented by those of the Central and Eastern United States. Building upon the suite of models developed by the recent NGA East project, we construct a new scaled backbone ground motion model and calibrate its corresponding epistemic uncertainties. The resulting logic tree is shown to provide comparable hazard outcomes to the epistemic uncertainty modelling strategy adopted for the Eastern United States, despite the different approaches taken. Comparison with previous GMM selections for northeastern Europe, however, highlights key differences in short period accelerations resulting from new assumptions regarding the characteristics of the reference rock and its influence on site amplification.
Probabilistic assessment of seismic hazard and risk over a geographical region presents the modeler with challenges in the characterization of the site amplification that are not present in site-specific assessment. Using site-to-site residuals from a ground motion model fit to observations from the Japanese KiK-net database, correlations between measured local amplifications and mappable proxies such as topographic slope and geology are explored. These are used subsequently to develop empirical models describing amplification as a direct function of slope, conditional upon geological period. These correlations also demonstrate the limitations of inferring 30-m shearwave velocity from slope and applying them directly into ground motion models. Instead, they illustrate the feasibility of deriving spectral acceleration amplification factors directly from sets of observed records, which are calibrated to parameters that can be mapped uniformly on a regional scale. The result is a geologically calibrated amplification model that can be incorporated into national and regional seismic hazard and risk assessment, ensuring that the corresponding total aleatory variability reflects the predictive capability of the mapped site proxy.
High-resolution surface velocities and strain for Anatolia from Sentinel-1 InSAR and GNSS data
(2020)
Measurements of present-day surface deformation are essential for the assessment of long-term seismic hazard. The European Space Agency's Sentinel-1 satellites enable global, high-resolution observation of crustal motion from Interferometric Synthetic Aperture Radar (InSAR). We have developed automated InSAR processing systems that exploit the first similar to 5 years of Sentinel-1 data to measure surface motions for the similar to 800,000-km(2) Anatolian region. Our new 3-D velocity and strain rate fields illuminate deformation patterns dominated by westward motion of Anatolia relative to Eurasia, localized strain accumulation along the North and East Anatolian Faults, and rapid vertical signals associated with anthropogenic activities and to a lesser extent extension across the grabens of western Anatolia. We show that automatically processed Sentinel-1 InSAR data can characterize details of the velocity and strain rate fields with high resolution and accuracy over large regions. These results are important for assessing the relationship between strain accumulation and release in earthquakes. <br /> Plain Language Summary Satellite-based measurements of small rates of motion of the Earth's surface made at high spatial resolutions and over large areas are important for many geophysical applications including improving earthquake hazard models. We take advantage of recent advances in geodetic techniques in order to measure surface velocities and tectonic strain accumulation across the Anatolia region, including the highly seismogenic and often deadly North Anatolian Fault. We show that by combining Sentinel-1 Interferometric Synthetic Aperture Radar (InSAR) data with Global Navigation Satellite System (GNSS) measurements we can enhance our view of surface deformation associated with active tectonics, the earthquake cycle, and anthropogenic processes.
Geochemical processes change the microstructure of rocks and thereby affect their physical behaviour at the macro scale. A micro-computer tomography (micro-CT) scan of a typical reservoir sandstone is used to numerically examine the impact of three spatial alteration patterns on pore morphology, permeability and elastic moduli by correlating precipitation with the local flow velocity magnitude. The results demonstrate that the location of mineral growth strongly affects the permeability decrease with variations by up to four orders in magnitude. Precipitation in regions of high flow velocities is characterised by a predominant clogging of pore throats and a drastic permeability reduction, which can be roughly described by the power law relation with an exponent of 20. A continuous alteration of the pore structure by uniform mineral growth reduces the permeability comparable to the power law with an exponent of four or the Kozeny-Carman relation. Preferential precipitation in regions of low flow velocities predominantly affects smaller throats and pores with a minor impact on the flow regime, where the permeability decrease is considerably below that calculated by the power law with an exponent of two. Despite their complete distinctive impact on hydraulics, the spatial precipitation patterns only slightly affect the increase in elastic rock properties with differences by up to 6.3% between the investigated scenarios. Hence, an adequate characterisation of the spatial precipitation pattern is crucial to quantify changes in hydraulic rock properties, whereas the present study shows that its impact on elastic rock parameters is limited. The calculated relations between porosity and permeability, as well as elastic moduli can be applied for upscaling micro-scale findings to reservoir-scale models to improve their predictive capabilities, what is of paramount importance for a sustainable utilisation of the geological subsurface.
Relative pollen productivity (RPP) estimates are fractionate values, often in relation to Poaceae, that allow vegetation cover to be estimated from pollen counts with the help of models. RPP estimates are especially used in the scientific community in Europe and China, with a few studies in North America. Here we present a comprehensive compilation of available northern hemispheric RPP studies and their results arising from 51 publications with 60 sites and 131 taxa. This compilation allows scientists to identify data gaps in need of further RPP analyses but can also aid them in finding an RPP set for their study region. We also present a taxonomically harmonised, unified RPP dataset for the Northern Hemisphere and subsets for North America (including Greenland), Europe (including arctic Russia), and China, which we generated from the available studies. The unified dataset gives the mean RPP for 55 harmonised taxa as well as fall speeds, which are necessary to reconstruct vegetation cover from pollen counts and RPP values. Data are openly available at https://doi.org/10.1594/PANGAEA.922661 (Wieczorek and Herzschuh, 2020).
Different upper tail indicators exist to characterize heavy tail phenomena, but no comparative study has been carried out so far. We evaluate the shape parameter (GEV), obesity index, Gini index and upper tail ratio (UTR) against a novel benchmark of tail heaviness - the surprise factor. Sensitivity analyses to sample size and changes in scale-to-location ratio are carried out in bootstrap experiments. The UTR replicates the surprise factor best but is most uncertain and only comparable between records of similar length. For samples with symmetric Lorenz curves, shape parameter, obesity and Gini indices provide consistent indications. For asymmetric Lorenz curves, however, the first two tend to overestimate, whereas Gini index tends to underestimate tail heaviness. We suggest the use of a combination of shape parameter, obesity and Gini index to characterize tail heaviness. These indicators should be supported with calculation of the Lorenz asymmetry coefficients and interpreted with caution.
In recent years, nature-based solutions are receiving increasing attention in the field of disaster risk reduction and climate change adaptation as inclusive, no regret approaches. Ecosystem-based adaptation (EbA) can mitigate the impacts of climate change, build resilience and tackle environmental degradation thereby supporting the targets set by the 2030 Agenda, the Paris Agreement and the Sendai Framework. Despite these benefits, EbA is still rarely implemented in practice. To better understand the barriers to implementation, this research examines policy-makers' perceptions of EbA, using an extended version of Protection Motivation Theory as an analytical framework. Through semi-structured interviews with policy-makers at regional and provincial level in Central Vietnam, it was found that EbA is generally considered a promising response option, mainly due to its multiple ecosystem-service benefits. The demand for EbA measures was largely driven by the perceived consequences of natural hazards and climate change. Insufficient perceived response efficacy and time-lags in effectiveness for disaster risk reduction were identified as key impediments for implementation. Pilot projects and capacity building on EbA are important means to overcome these perceptual barriers. This paper contributes to bridging the knowledge-gap on political decision-making regarding EbA and can, thereby, promote its mainstreaming into policy plans.
Marked along-strike changes in stratigraphy, mountain belt morphology, basement exhumation, and deformation styles characterize the Andean retroarc; these changes have previously been related to spatiotemporal variations in the subduction angle. We modeled new apatite fission track and apatite (U-Th-Sm)/He data from nine ranges located between 26 degrees S and 28 degrees S. Using new and previously published data, we constructed a Cretaceous to Pliocene paleogeographic model that delineates a four-stage tectonic evolution: extensional tectonics during the Cretaceous (120-75 Ma), the formation of a broken foreland basin between 55 and 30 Ma, reheating due to burial beneath sedimentary rocks (18-13 Ma), and deformation, exhumation, and surface uplift during the Late Miocene and the Pliocene (13-3 Ma). Our model highlights how preexisting upper plate structures control the deformation patterns of broken foreland basins. Because retroarc deformation predates flat-slab subduction, we propose that slab anchoring may have been the precursor of Eocene-Oligocene compression in the Andean retroarc. Our model challenges models which consider broken foreland basins and retroarc deformation in the NW Argentinian Andes to be directly related to Miocene flat subduction.
Millennial-scale palaeoclimate variability has been documented in various terrestrial and marine palaeoclimate proxy records throughout the Northern Hemisphere for the last glacial cycle. Its clear expression and rapid shifts between different states of climate (Greenland Interstadials and Stadials) represents a correlation tool beyond the resolution of e.g. luminescence dating, especially relevant for terrestrial deposits. Usually, comparison of terrestrial proxy datasets and the Greenland ice cores indicates a complex expression of millennial-scale climate variability as recorded in terrestrial geoarchives including loess. Loess is the most widespread terrestrial geoarchive of the Quaternary and especially widespread over Eurasia. However, loess often records a smoothed representation of millennial-scale variability without all fidelity when compared to the Greenland data, this being a relevant limiting feature in integrating loess with other palaeoclimate records. To better understand the loess proxy-response to millennial-scale climate variability, we simulate a proxy signal smoothing by natural processes through application of low-pass filters of delta O-18 data from Greenland, a high-resolution palaeoclimate reference record, alongside speleothem isotope records from the Black Sea-Mediterranean region. We show that low-pass filters represent rather simple models for better constraining the expression of millennial-scale climate variability in low sedimentation environments, and in sediments where proxy-response signals are most likely affected by natural smoothing (by e.g. bioturbation). Interestingly, smoothed datasets from Greenland and the Black Sea-Mediterranean region are most similar in the last similar to 15 ka and between similar to 50-30 ka. Between similar to 30-15 ka, roughly corresponding to the Last Glacial Maximum and the deglaciation, the records show dissimilarities, challenging the construction of robust correlative time-scales in this age range. From our analysis it becomes apparent that patterns of palaeoclimate signals in loess-palaeosol sequences often might be better explained by smoothed Greenland reference data than the original high-resolution Greenland dataset, or other reference data. This opens the possibility to better assess the temporal resolution and palaeoclimate potential of loess-palaeosol sequences in recording supra-regional climate patterns, as well as to securely integrate loess with other chronologically better-resolved palaeoclimate records.