Refine
Year of publication
- 2019 (276) (remove)
Document Type
- Article (216)
- Doctoral Thesis (24)
- Postprint (15)
- Other (9)
- Review (7)
- Habilitation Thesis (3)
- Conference Proceeding (1)
- Master's Thesis (1)
Is part of the Bibliography
- yes (276)
Keywords
- climate (5)
- remote sensing (4)
- Andes (3)
- Climate (3)
- Holocene (3)
- geochronology (3)
- melt inclusions (3)
- paleoclimate (3)
- photogrammetry (3)
- silicon (3)
Institute
- Institut für Geowissenschaften (276) (remove)
Glacial-interglacial variations in CO2 and methane in polar ice cores have been attributed, in part, to changes in global wetland extent, but the wetland distribution before the Last Glacial Maximum (LGM, 21 ka to 18 ka) remains virtually unknown. We present a study of global peatland extent and carbon (C) stocks through the last glacial cycle (130 ka to present) using a newly compiled database of 1,063 detailed stratigraphic records of peat deposits buried by mineral sediments, as well as a global peatland model. Quantitative agreement between modeling and observations shows extensive peat accumulation before the LGM in northern latitudes (> 40 degrees N), particularly during warmer periods including the last interglacial (130 ka to 116 ka, MIS 5e) and the interstadial (57 ka to 29 ka, MIS 3). During cooling periods of glacial advance and permafrost formation, the burial of northern peatlands by glaciers and mineral sediments decreased active peatland extent, thickness, and modeled C stocks by 70 to 90% from warmer times. Tropical peatland extent and C stocks show little temporal variation throughout the study period. While the increased burial of northern peats was correlated with cooling periods, the burial of tropical peat was predominately driven by changes in sea level and regional hydrology. Peat burial by mineral sediments represents a mechanism for long-term terrestrial C storage in the Earth system. These results show that northern peatlands accumulate significant C stocks during warmer times, indicating their potential for C sequestration during the warming Anthropocene.
This study aims to identify the best-performing site characterization proxy alternative and complementary to the conventional 30 m average shear-wave velocity V-S30, as well as the optimal combination of proxies in characterizing linear site response. Investigated proxies include T-0 (site fundamental period obtained from earthquake horizontal-to-vertical spectral ratios), V-Sz (measured average shear-wave velocities to depth z, z = 5, 10, 20 and 30 m), Z(0.8) and Z(1.0) (measured site depths to layers having shear-wave velocity 0.8 and 1.0 km/s, respectively), as well as Z(x-infer) (inferred site depths from a regional velocity model, x = 0.8 and 1.0, 1.5 and 2.5 km/s). To evaluate the performance of a site proxy or a combination, a total of 1840 surface-borehole recordings is selected from KiK-net database. Site amplifications are derived using surface-to-borehole response-, Fourier- and cross-spectral ratio techniques and then are compared across approaches. Next, the efficacies of 7 single-proxies and 11 proxy-pairs are quantified based on the site-to-site standard deviation of amplification residuals of observation about prediction using the proxy or the pair. Our results show that T-0 is the best-performing single-proxy among T-0, Z(0.8), Z(1.0) and V-Sz. Meanwhile, T-0 is also the best-performing proxy among T-0, Z(0.8), Z(1.0) and Z(x-infer) complementary to V-S30 in accounting for the residual amplification after V-S30-correction. Besides, T-0 alone can capture most of the site effects and should be utilized as the primary site indicator. Though (T-0, V-S30) is the best-performing proxy pair among (V-S30, T-0), (V-S30, Z(0.8)), (V-S30, Z(1.0)), (V-S30, Z(x-infer)) and (T-0, V-Sz), it is only slightly better than (T-0, V-S20). Considering both efficacy and engineering utility, the combination of T-0 (primary) and V-S20 (secondary) is recommended. Further study is needed to test the performances of various proxies on sites in deep sedimentary basins.
The geochemical composition of oceanic basalts provides us with a window into the distribution of geochemical elements within the Earth’s mantle in space and time. In conjunction with a throughout knowledge on how the different elements behave e.g. during melt formation and evolution or on their partition behaviour between e.g. minerals and melts this information has been transformed into various models on how oceanic crust is formed along plume influenced or normal mid-ocean ridge segments, how oceanic crust evolves in response to seawater, on subduction recycling of oceanic crust and so forth. The work presented in this habilitation was aimed at refining existing models, putting further constraints on some of the major open questions in this field of research while at the same time trying to increase our knowledge on the behaviour of noble gases as a tracer for melt formation and evolution processes. In the line of this work the author and her co-workers were able to answer one of the major questions concerning the formation of oceanic crust along plume-influenced ridges – in which physical state does the plume material enter the ridge? Based on submarine volcanic glass He, Ne and Ar data, the author and her co-workers have shown that the interaction of mantle plumes with mid-ocean ridges occurs in the physical form of melts. In addition, the author and her co-workers have also put further constraints on one of the major questions concerning the formation of oceanic crust along normal mid-ocean ridges – namely how is the mid-ocean ridge system effectively cooled to form the lower oceanic crust? Based on Ne and Ar data in combination with Cl/K ratios of basaltic glass from the Mid-Atlantic ridge and estimates of crystallisation pressures they have shown, that seawater penetration reaches lower crustal levels close to the Moho, indicating that hydrothermal circulation might be an effective cooling mechanism even for the deep parts of the oceanic crust. Considering subduction recycling, the heterogeneity of the Earth’s mantle and mantle dynamic processes the key question is on which temporal and spatial scales is the Earth’s mantle geochemically heterogeneous? In the line of this work the author along with her co-workers have shown based on Cl/K ratios in conjunction with the Sr, Nd, and Pb isotopes of the OIBs representing the type localities for the different mantle endmembers that the quantity of Cl recycled into the mantle via subduction is not uniform and that neither the HIMU nor the EM1 and EM2 mantle components can be considered as distinct mantle endmembers. In addition, we have shown, based on He, Ne and Ar isotope and trace-element data from the Foundation hotspot that the near ridge seamounts of the Foundation seamount chain formed by the Foundation hotspot erupt lavas with a trace-element signature clearly characteristic of oceanic gabbro which indicates the existence of recycled, virtually unchanged lower oceanic crust in the plume source. This is a clear sign of the inefficiency of the stirring mechanism existing at mantle depth. Similar features are seen in other near-axis hotspot magmas around the world. Based on He, Sr, Nd, Pb and O isotopes and trace elements in primitive mafic dykes from the Etendeka flood basalts, NW Namibia the author along with her co-workers have shown that deep, less degassed mantle material carried up by a mantle plume contributed significantly to the flood basalt magmatism. The Etendeka flood basalts are part of the South Atlantic LIP, which is associated with the breakup of Gondwana, the formation of the Paraná-Etendeka flood basalts and the Walvis Ridge - Tristan da Cunha hotspot track. Thus reinforcing the lately often-challenged concept of mantle plumes and the role of mantle plumes in the formation of large igneous provinces. Studying the behaviour of noble gases during melt formation and evolution the author along with her co-workers has shown that He can be considerable more susceptible to changes during melt formation and evolution resulting not only in a complete decoupling of He isotopes from e.g. Ne or Pb isotopes but also in a complete loss of the primary mantle isotope signal. They have also shown that this decoupling occurs mainly during the melt formation processes requiring He to be more compatible during mantle melting than Ne. In addition, the author along with her co workers were able to show that incorporation of atmospheric noble gases into igneous rocks is in general a two-step process: (1) magma contamination by assimilation of altered oceanic crust results in the entrainment of air-equilibrated seawater noble gases; (2) atmospheric noble gases are adsorbed onto grain surfaces during sample preparation. This implies, considering the ubiquitous presence of the contamination signal, that magma contamination by assimilation of a seawater-sourced component is an integral part of mid-ocean ridge basalt evolution.
Visible-Light-Mediated Photodynamic Water Disinfection @ Bimetallic-Doped Hybrid Clay Nanocomposites
(2019)
This study reports a new class of photocatalytic hybrid clay nanocomposites prepared from low-cost sources (kaolinite clay and Carica papaya seeds) doped with Zn and Cu salts via a solvothermal process. X-ray diffraction analysis suggests that Cu-doping and Cu/Zn-doping introduce new phases into the crystalline structure of Kaolinite clay, which is linked to the reduced band gap of kaolinite from typically between 4.9 and 8.2 eV to 2.69 eV for Cu-doped and 1.5 eV for Cu/Zn hybrid clay nanocomposites (Nisar, J.; Arhammar, C.; Jamstorp, E.; Ahuja, R. Phys. Rev. B 2011, 84, 075120). In the presence of solar light irradiation, Cu- and Cu/Zn-doped nanocomposites facilitate the electron hole pair separation. This promotes the generation of singlet oxygen which in turn improves the water disinfection efficiencies of these novel nanocomposite materials. The nanocomposite materials were further characterized using high-resolution scanning electron microscopy, fluorimetry, therrnogravimetric analysis, and Raman spectroscopy. The breakthrough times of the nanocomposites for a fixed bed mode of disinfection of water contaminated with 2.32 x 10(7) cfu/mL E. coli ATCC 25922 under solar light irradiation are 25 h for Zn-doped, 30 h for Cu-doped, and 35 h for Cu/Zn-doped nanocomposites. In the presence of multidrug and multimetal resistant strains of E. coli, the breakthrough time decreases significantly. Zn-only doped nanocomposites are not photocatalytically active. In the absence of light, the nanocomposites are still effective in decontaminating water, although less efficient than under solar light irradiation. Electrostatic interaction, metal toxicity, and release of singlet oxygen (only in the Cu-doped and Cu/Zn-doped nanocomposites) are the three disinfection mechanisms by which these nanocomposites disinfect water. A regrowth study indicates the absence of any living E. coli cells in treated water even after 4 days. These data and the long hydraulic times (under gravity) exhibited by these nanocomposites during photodisinfection of water indicate an unusually high potential of these nanocomposites as efficient, affordable, and sustainable point-of-use systems for the disinfection of water in developing countries.
The sedimentary record of the Dead Sea provides an exceptional high-resolution archive of past climate changes in the drought-sensitive eastern Mediterranean-Levant, a key region for the development of humankind at the boundary of global climate belts. Moreover, it is the only deep hypersaline lake known to have deposited long sequences of finely laminated, annually deposited sediments (i.e. varves) of varied compositions, including aragonite, gypsum, halite and clastic sediments. Vast efforts have been made over the years to decipher the environmental information stored in these evaporitic-clastic sequences spanning from the Pleistocene Lake Amora to the Holocene Dead Sea. A general characterisation of sediment facies has been derived from exposed sediment sections, as well as from shallow- and deep-water sediment cores. During high lake stands and episodes of positive water budget, mostly during glacial times, alternating aragonite and detritus laminae (‘aad’ facies) were accumulated, whereas during low lake stands and droughts, prevailing during interglacials, laminated detritus (‘ld’ facies) and laminated halite (‘lh’ facies) dominate the sequence. In this paper, we (i) review the three types of laminated sediments of the Dead Sea sedimentary record (‘aad’, ‘ld’ and ‘lh’ facies), (ii) discuss their modes of formation, deposition and accumulation, and their interpretation as varves, and (iii) illustrate how Dead Sea varves are utilized for palaeoclimate reconstructions and for establishing floating chronologies.
Small-scale variations in mineral chemistry, textures, and platinum group element (PGE) mineralization were investigated in the Lower and Middle Group chromitite layers LG6, LG6a, MG1, MG2, and MG2 II from vertical drill core profiles at the Thaba mine in the northwestern limb of the Bushveld Complex. We present detailed geochemical profiles of chromite composition and chromite crystal size distribution curves to shed light on the processes of chromite accumulation and textural modification as well as mineralization. Multiple samples within each layer were assayed for PGE concentrations, and the respective platinum group mineral association was determined by mineral liberation analysis (MLA). There is strong evidence for postcumulus changes in the chromitites. The crystal size distribution curves suggest that the primary chromite texture was coarsened by a combination of adcumulus growth and textural equilibration, while compaction of the crystal mush played only a minor role. Mineral compositions were also modified by postcumulus processes, but because of the very high modal amount of chromite and its local preservation in orthopyroxene oikocrysts, that phase retained much primary information. Vertical variations of chromite composition within chromitite layers and from one layer to another do not support the idea of chromite accumulation from crystal-rich slurries or crystal settling from a large magma chamber. Instead, we favor a successive buildup of chromitite layers by repeated injections of relatively thin layers of chromite-saturated magmas, with in situ crystallization occurring at the crystal mush-magma interface. The adcumulus growth of chromite grains to form massive chromitite required addition of Cr to the layers, which we attribute to downward percolation from the overlying magma. The PGE concentrations are elevated in all chromitite layers compared to adjacent silicate rocks and show a systematic increase upward from LG6 (avg 807 ppb Ir + Ru + Rh + Pt + Pd + Au) to MG2 II (avg 2,062 ppb). There are also significant internal variations in all layers, with enrichments at hanging and/or footwalls. The enriched nature of chromitites in PGEs compared to host pyroxenites is a general feature, independent of the layer thickness. The MLA results distinguish two principal groups of PGE mineral associations: the LG6, LG6, and MG1 are dominated by the malanite series, laurite, and PGE sulfarsenides, while the MG2 and MG2 II layers are characterized by laurite and PGE sulfides as well as Pt-Fe-Sn and PGE-Sb-Bi-Pb alloys. Differences in the PGE associations are attributed to postcumulus alteration of the MG2 and MG2 II layer, while the chromitites below, particularly LG6 and LG6a, contain a more pristine association.
Variability of the Cold Season Climate in Central Asia. Part II: Hydroclimatic Predictability
(2019)
Central Asia (CA) is subjected to a large variability of precipitation. This study presents a statistical model, relating precipitation anomalies in three subregions of CA in the cold season (November-March) with various predictors in the preceding October. Promising forecast skill is achieved for two subregions covering 1) Uzbekistan, Turkmenistan, Kyrgyzstan, Tajikistan, and southern Kazakhstan and 2) Iran, Afghanistan, and Pakistan. ENSO in October is identified as the major predictor. Eurasian snow cover and the quasi-biennial oscillation further improve the forecast performance. To understand the physical mechanisms, an analysis of teleconnections between these predictors and the wintertime circulation over CA is conducted. The correlation analysis of predictors and large-scale circulation indices suggests a seasonal persistence of tropical circulation modes and a dynamical forcing of the westerly circulation by snow cover variations over Eurasia. An EOF analysis of pressure and humidity patterns allows separating the circulation variability over CA into westerly and tropical modes and confirms that the identified predictors affect the respective circulation characteristics. Based on the previously established weather type classification for CA, the predictors are investigated with regard to their effect on the regional circulation. The results suggest a modification of the Hadley cell due to ENSO variations, with enhanced moisture supply from the Arabian Gulf during El Nino. They further indicate an influence of Eurasian snow cover on the wintertime Arctic Oscillation (AO) and Northern Hemispheric Rossby wave tracks. Positive anomalies favor weather types associated with dry conditions, while negative anomalies promote the formation of a quasi-stationary trough over CA, which typically occurs during positive AO conditions.
The interactions between atmosphere and steep topography in the eastern south–central Andes result in complex relations with inhomogenous rainfall distributions. The atmospheric conditions leading to deep convection and extreme rainfall and their spatial patterns—both at the valley and mountain-belt scales—are not well understood. In this study, we aim to identify the dominant atmospheric conditions and their spatial variability by analyzing the convective available potential energy (CAPE) and dew-point temperature (Td). We explain the crucial effect of temperature on extreme rainfall generation along the steep climatic and topographic gradients in the NW Argentine Andes stretching from the low-elevation eastern foreland to the high-elevation central Andean Plateau in the west. Our analysis relies on version 2.0 of the ECMWF’s (European Centre for Medium-RangeWeather Forecasts) Re-Analysis (ERA-interim) data and TRMM (Tropical Rainfall Measuring Mission) data. We make the following key observations: First, we observe distinctive gradients along and across strike of the Andes in dew-point temperature and CAPE that both control rainfall distributions. Second, we identify a nonlinear correlation between rainfall and a combination of dew-point temperature and CAPE through a multivariable regression analysis. The correlation changes in space along the climatic and topographic gradients and helps to explain controlling factors for extreme-rainfall generation. Third, we observe more contribution (or higher importance) of Td in the tropical low-elevation foreland and intermediate-elevation areas as compared to the high-elevation central Andean Plateau for 90th percentile rainfall. In contrast, we observe a higher contribution of CAPE in the intermediate-elevation area between low and high elevation, especially in the transition zone between the tropical and subtropical areas for the 90th percentile rainfall. Fourth, we find that the parameters of the multivariable regression using CAPE and Td can explain rainfall with higher statistical significance for the 90th percentile compared to lower rainfall percentiles. Based on our results, the spatial pattern of rainfall-extreme events during the past ∼16 years can be described by a combination of dew-point temperature and CAPE in the south–central Andes.
The interactions between atmosphere and steep topography in the eastern south–central Andes result in complex relations with inhomogenous rainfall distributions. The atmospheric conditions leading to deep convection and extreme rainfall and their spatial patterns—both at the valley and mountain-belt scales—are not well understood. In this study, we aim to identify the dominant atmospheric conditions and their spatial variability by analyzing the convective available potential energy (CAPE) and dew-point temperature (Td). We explain the crucial effect of temperature on extreme rainfall generation along the steep climatic and topographic gradients in the NW Argentine Andes stretching from the low-elevation eastern foreland to the high-elevation central Andean Plateau in the west. Our analysis relies on version 2.0 of the ECMWF’s (European Centre for Medium-RangeWeather Forecasts) Re-Analysis (ERA-interim) data and TRMM (Tropical Rainfall Measuring Mission) data. We make the following key observations: First, we observe distinctive gradients along and across strike of the Andes in dew-point temperature and CAPE that both control rainfall distributions. Second, we identify a nonlinear correlation between rainfall and a combination of dew-point temperature and CAPE through a multivariable regression analysis. The correlation changes in space along the climatic and topographic gradients and helps to explain controlling factors for extreme-rainfall generation. Third, we observe more contribution (or higher importance) of Td in the tropical low-elevation foreland and intermediate-elevation areas as compared to the high-elevation central Andean Plateau for 90th percentile rainfall. In contrast, we observe a higher contribution of CAPE in the intermediate-elevation area between low and high elevation, especially in the transition zone between the tropical and subtropical areas for the 90th percentile rainfall. Fourth, we find that the parameters of the multivariable regression using CAPE and Td can explain rainfall with higher statistical significance for the 90th percentile compared to lower rainfall percentiles. Based on our results, the spatial pattern of rainfall-extreme events during the past ∼16 years can be described by a combination of dew-point temperature and CAPE in the south–central Andes.
Shrinking glaciers in the Hindu Kush-Karakoram-Himalaya-Nyainqentanglha (HKKHN) region have formed several thousand moraine-dammed glacial lakes(1-3), some of these having grown rapidly in past decades(3,4). This growth may promote more frequent and potentially destructive glacial lake outburst floods (GLOFs)(5-7). Testing this hypothesis, however, is confounded by incomplete databases of the few reliable, though selective, case studies. Here we present a consistent Himalayan GLOF inventory derived automatically from all available Landsat imagery since the late 1980s. We more than double the known GLOF count and identify the southern Himalayas as a hotspot region, compared to the more rarely affected Hindu Kush-Karakoram ranges. Nevertheless, the average annual frequency of 1.3 GLOFs has no credible posterior trend despite reported increases in glacial lake areas in most of the HKKHN3,8, so that GLOF activity per unit lake area has decreased since the late 1980s. We conclude that learning more about the frequency and magnitude of outburst triggers, rather than focusing solely on rapidly growing glacial lakes, might improve the appraisal of GLOF hazards.
Skarn deposits are found on every continents and were formed at different times from Precambrian to Tertiary. Typically, the formation of a skarn is induced by a granitic intrusion in carbonates-rich sedimentary rocks. During contact metamorphism, fluids derived from the granite interact with the sedimentary host rocks, which results in the formation of calc-silicate minerals at the expense of carbonates. Those newly formed minerals generally develop in a metamorphic zoned aureole with garnet in the proximal and pyroxene in the distal zone. Ore elements contained in magmatic fluids are precipitated due to the change in fluid composition. The temperature decrease of the entire system, due to the cooling of magmatic fluids and the entering of meteoric water, allows retrogression of some prograde minerals.
The Hämmerlein skarn deposit has a multi-stage history with a skarn formation during regional metamorphism and a retrogression of primary skarn minerals during the granitic intrusion. Tin was mobilized during both events. The 340 Ma old tin-bearing skarn minerals show that tin was present in sediments before the granite intrusion, and that the first Sn enrichment occurred during the skarn formation by regional metamorphism fluids. In a second step at ca. 320 Ma, tin-bearing fluids were produced with the intrusion of the Eibenstock granite. Tin, which has been added by the granite and remobilized from skarn calc-silicates, precipitated as cassiterite.
Compared to clay or marl, the skarn is enriched in Sn, W, In, Zn, and Cu. These metals have been supplied during both regional metamorphism and granite emplacement. In addition, the several isotopic and chemical data of skarn samples show that the granite selectively added elements such as Sn, and that there was no visible granitic contribution to the sedimentary signature of the skarn
The example of Hämmerlein shows that it is possible to form a tin-rich skarn without associated granite when tin has already been transported from tin-bearing sediments during regional metamorphism by aqueous metamorphic fluids. These skarns are economically not interesting if tin is only contained in the skarn minerals. Later alteration of the skarn (the heat and fluid source is not necessarily a granite), however, can lead to the formation of secondary cassiterite (SnO2), with which the skarn can become economically highly interesting.
We analyze trends in compound flooding resulting from high coastal water levels (HCWLs) and peak river discharge over northwestern Europe during 1901-2014. Compound peak discharge associated with 37 stream gauges with at least 70 years of record availability near the North and Baltic Sea coasts is used. Compound flooding is assessed using a newly developed index, compound hazard ratio, that compares the severity of river flooding associated with HCWL with the at-site, T-year (a flood with 1/T chance of being exceeded in any given year) fluvial peak discharge. Our findings suggest a spatially coherent pattern in the dependence between HCWL and river peaks and in compound flood magnitudes and frequency. For higher return levels, we find upward trends in compound hazard ratio frequency at midlatitudes (gauges from 47 degrees N to 60 degrees N) and downward trends along the high latitude (>60 degrees N) regions of northwestern Europe. Plain Language Summary Compound floods in delta areas, that is, the co-occurrence of high coastal water levels (HCWLs) and high river discharge, are a particular challenge for disaster management. Such events are caused by two distinct mechanisms: (1) HCWLs may affect river flows and water levels by backwater effects or by reversing the seaward flow of rivers, particularly in regions with elevation less than 10 m in northwestern Europe. (2) The correlation between HCWL and river flow peaks may also stem from a common meteorological driver. Severe storm periods may be associated with high winds leading to storm surges, and at the same time with high precipitation followed by inland flooding. Understanding the historical trends in compound flooding, owing to changes in relative sea levels, in river flooding and in the dependence between these two drivers, is essential for projecting future changes and disaster management. The risk assessment frameworks are often limited to assessing flood risk from a single driver only. We present a new approach to assess compound flood severity resulting from extreme coastal water level and peak river discharge. We find upward trends in compound flooding for midlatitude regions and downward trends for high latitudes in northwestern Europe.
Different tectonic episodes from Late Triassic to recent times in the eastern Binalud Mountains have resulted from convergence and transpression between the Turan and Central Iran plates. Heterogeneous deformation and variable portions of pure and simple shear, demonstrated by finite strain and vorticity analysis in the Mashhad metamorphic rocks, indicate strain partitioning during the first tectonic episode. Modern strain partitioning is characterized by reverse and strike-slip faulting along the Neyshabur fault system and Shandiz fault zone in the southern and northern flanks of the eastern Binalud, respectively. Time-transgressive regional deformation migrated from the hinterland of the belt into the foreland basin, from northeast to the southwest of the mountains. Different generations of deformation resulted in obliteration of the subduction-related accretionary wedge, and growth of an orogenic wedge resulted from collision between the Central Iran and Turan plates.
Ring current electrons (1–100 keV) have received significant attention in recent decades, but many questions regarding their major transport and loss mechanisms remain open. In this study, we use the four‐dimensional Versatile Electron Radiation Belt code to model the enhancement of phase space density that occurred during the 17 March 2013 storm. Our model includes global convection, radial diffusion, and scattering into the Earth's atmosphere driven by whistler‐mode hiss and chorus waves. We study the sensitivity of the model to the boundary conditions, global electric field, the electric field associated with subauroral polarization streams, electron loss rates, and radial diffusion coefficients. The results of the code are almost insensitive to the model parameters above 4.5 RERE, which indicates that the general dynamics of the electrons between 4.5 RE and the geostationary orbit can be explained by global convection. We found that the major discrepancies between the model and data can stem from the inaccurate electric field model and uncertainties in lifetimes. We show that additional mechanisms that are responsible for radial transport are required to explain the dynamics of ≥40‐keV electrons, and the inclusion of the radial diffusion rates that are typically assumed in radiation belt studies leads to a better agreement with the data. The overall effect of subauroral polarization streams on the electron phase space density profiles seems to be smaller than the uncertainties in other input parameters. This study is an initial step toward understanding the dynamics of these particles inside the geostationary orbit.
Two principal groups of processes shape mass fluxes from and into a soil: vertical profile development and lateral soil redistribution. Periods having predominantly progressive soil forming processes (soil profile development) alternate with periods having predominantly regressive processes (erosion). As a result, short‐term soil redistribution – years to decades – can differ substantially from long‐term soil redistribution; i.e. centuries to millennia. However, the quantification of these processes is difficult and consequently their rates are poorly understood. To assess the competing roles of erosion and deposition we determined short‐ and long‐term soil redistribution rates in a formerly glaciated area of the Uckermark, northeast Germany. We compared short‐term erosion or accumulation rates using plutonium‐239 and ‐240 (239+240Pu) and long‐term rates using both in situ and meteoric cosmogenic beryllium‐10 (10Be). Three characteristic process domains have been analysed in detail: a flat landscape position having no erosion/deposition, an erosion‐dominated mid‐slope, and a deposition‐dominated lower‐slope site. We show that the short‐term mass erosion and accumulation rates are about one order of magnitude higher than long‐term redistribution rates. Both, in situ and meteoric 10Be provide comparable results. Depth functions, and therefore not only an average value of the topsoil, give the most meaningful rates. The long‐term soil redistribution rates were in the range of −2.1 t ha‐1 yr‐1 (erosion) and +0.26 t ha‐1 yr‐1 (accumulation) whereas the short‐term erosion rates indicated strong erosion of up to 25 t ha‐1 yr‐1 and accumulation of 7.6 t ha‐1 yr‐1. Our multi‐isotope method identifies periods of erosion and deposition, confirming the ‘time‐split approach’ of distinct different phases (progressive/regressive) in soil evolution. With such an approach, temporally‐changing processes can be disentangled, which allows the identification of both the dimensions of and the increase in soil erosion due to human influence
We investigate the thermal and structural evolution of asymmetric rifted margin using numerical modeling and geological observations derived from the Western Pyrenees. Our numerical model provides a self-consistent physical evolution of the top basement heat flow during asymmetric rifting. The model shows a pronounced thermal asymmetry that is caused by migration of the rift center toward the upper plate. The same process creates a diachronism for the record of maximum heat flow and maximum temperatures (T-max) in basal rift sequences. The Mauleon-Arzacq basin (W-Pyrenees) corresponds to a former mid-Cretaceous asymmetric hyperextended rift basin. New vitrinite reflectance data in addition to existing data sets from this basin reveal an asymmetry in the distribution of peak heat (T-max) with respect to the rift shoulders, where highest values are located at the former upper- to lower-plate transition. This data set from the Arzacq-Mauleon field study confirms for the first time the thermal asymmetry predicted by numerical models. Numerical modeling results also suggest that complexities in synrift thermal architecture could arise when hanging-wall-derived extensional allochthons and related T-max become part of the lower plate and are transported away from the upper- to lower-plate transition. This study emphasizes the limitations of the common approach to integrate punctual thermal data from pre-rift to synrift sedimentary sequences in order to describe the rift-related thermal evolution and paleothermal gradients at the scale of a rift basin or a rifted margin.
The thermal structure of subduction zones exerts a major influence on deep-seated mechanical and chemical processes controlling arc magmatism, seismicity, and global element cycles. Accretionary complexes exposed inland may comprise tectonic blocks with contrasting pressure-temperature (P-T) histories, making it possible to investigate the dynamics and thermal evolution of former subduction interfaces. With this aim, we present new Lu-Hf geochronological results for mafic rocks of the Halilbagi Complex (Anatolia) that evolved along different thermal gradients. Samples include a lawsonite-epidote blueschist, a lawsonite-epidote eclogite, and an epidote eclogite (all with counter-clockwise P-T paths), a prograde lawsonite blueschist with a "hairpin"-type P-T path, and a garnet amphibolite from the overlying sub-ophiolitic metamorphic sole. Equilibrium phase diagrams suggest that the garnet amphibolite formed at similar to 0.6-0.7 GPa and 800-850 degrees C, whereas the prograde lawsonite blueschist records burial from 2.1 GPa and 420 degrees C to 2.6 GPa and 520 degrees C. Well-defined Lu-Hf isochrons were obtained for the epidote eclogite (92.38 +/- 0.22 Ma) and the lawsonite-epidote blueschist (90.19 +/- 0.54 Ma), suggesting rapid garnet growth. The lawsonite-epidote eclogite (87.30 +/- 0.39 Ma) and the prograde lawsonite blueschist (ca. 86 Ma) are younger, whereas the garnet amphibolite (104.5 +/- 3.5 Ma) is older. Our data reveal a consistent trend of progressively decreasing geothermal gradient from granulite-facies conditions at similar to 104 Ma to the epidote-eclogite facies around 92 Ma, and the lawsonite blueschist-facies between 90 Ma and 86 Ma. Three Lu-Hf garnet dates (between 92 Ma and 87 Ma) weighted toward the growth of post-peak rims (as indicated by Lu distribution in garnet) suggest that the HP/LT rocks were exhumed continuously and not episodically. We infer that HP/LT metamorphic rocks within the Halilbagi Complex were subjected to continuous return flow, with "warm" rocks being exhumed during the tectonic burial of "cold" ones. Our results, combined with regional geological constraints, allow us to speculate that subduction started at a transform fault near a mid-oceanic spreading centre. Following its formation, this ancient subduction interface evolved thermally over more than 15 Myr, most likely as a result of heat dissipation rather than crustal underplating. (C) 2018, China University of Geosciences (Beijing) and Peking University. Production and hosting by Elsevier B.V.
The foreland of the Andes in South America is characterised by distinct along strike changes in surface deformational styles. These styles are classified into two end-members, the thin-skinned and the thick-skinned style. The superficial expression of thin-skinned deformation is a succession of narrowly spaced hills and valleys, that form laterally continuous ranges on the foreland facing side of the orogen. Each of the hills is defined by a reverse fault that roots in a basal décollement surface within the sedimentary cover, and acted as thrusting ramp to stack the sedimentary pile. Thick-skinned deformation is morphologically characterised by spatially disparate, basement-cored mountain ranges. These mountain ranges are uplifted along reactivated high-angle crustal-scale discontinuities, such as suture zones between different tectonic terranes.
Amongst proposed causes for the observed variation are variations in the dip angle of the Nazca plate, variation in sediment thickness, lithospheric thickening, volcanism or compositional differences. The proposed mechanisms are predominantly based on geological observations or numerical thermomechanical modelling, but there has been no attempt to understand the mechanisms from a point of data-integrative 3D modelling. The aim of this dissertation is therefore to understand how lithospheric structure controls the deformational behaviour. The integration of independent data into a consistent model of the lithosphere allows to obtain additional evidence that helps to understand the causes for the different deformational styles. Northern Argentina encompasses the transition from the thin-skinned fold-and-thrust belt in Bolivia, to the thick-skinned Sierras Pampeanas province, which makes this area a well suited location for such a study. The general workflow followed in this study first involves data-constrained structural- and density-modelling in order to obtain a model of the study area. This model was then used to predict the steady-state thermal field, which was then used to assess the present-day rheological state in northern Argentina.
The structural configuration of the lithosphere in northern Argentina was determined by means of data-integrative, 3D density modelling verified by Bouguer gravity. The model delineates the first-order density contrasts in the lithosphere in the uppermost 200 km, and discriminates bodies for the sediments, the crystalline crust, the lithospheric mantle and the subducting Nazca plate. To obtain the intra-crustal density structure, an automated inversion approach was developed and applied to a starting structural model that assumed a homogeneously dense crust. The resulting final structural model indicates that the crustal structure can be represented by an upper crust with a density of 2800 kg/m³, and a lower crust of 3100 kg/m³. The Transbrazilian Lineament, which separates the Pampia terrane from the Río de la Plata craton, is expressed as a zone of low average crustal densities.
In an excursion, we demonstrate in another study, that the gravity inversion method developed to obtain intra-crustal density structures, is also applicable to obtain density variations in the uppermost lithospheric mantle. Densities in such sub-crustal depths are difficult to constrain from seismic tomographic models due to smearing of crustal velocities. With the application to the uppermost lithospheric mantle in the north Atlantic, we demonstrate in Tan et al. (2018) that lateral density trends of at least 125\,km width are robustly recovered by the inversion method, thereby providing an important tool for the delineation of subcrustal density trends.
Due to the genetic link between subduction, orogenesis and retroarc foreland basins the question rises whether the steady-state assumption is valid in such a dynamic setting. To answer this question, I analysed (i) the impact of subduction on the conductive thermal field of the overlying continental plate, (ii) the differences between the transient and steady-state thermal fields of a geodynamic coupled model. Both studies indicate that the assumption of a thermal steady-state is applicable in most parts of the study area. Within the orogenic wedge, where the assumption cannot be applied, I estimated the transient thermal field based on the results of the conducted analyses.
Accordingly, the structural model that had been obtained in the first step, could be used to obtain a 3D conductive steady-state thermal field. The rheological assessment based on this thermal field indicates that the lithosphere of the thin-skinned Subandean ranges is characterised by a relatively strong crust and a weak mantle. Contrarily, the adjacent foreland basin consists of a fully coupled, very strong lithosphere. Thus, shortening in northern Argentina can only be accommodated within the weak lithosphere of the orogen and the Subandean ranges. The analysis suggests that the décollements of the fold-and-thrust belt are the shallow continuation of shear zones that reside in the ductile sections of the orogenic crust. Furthermore, the localisation of the faults that provide strain transfer between the deeper ductile crust and the shallower décollement is strongly influenced by crustal weak zones such as foliation. In contrast to the northern foreland, the lithosphere of the thick-skinned Sierras Pampeanas is fully coupled and characterised by a strong crust and mantle. The high overall strength prevents the generation of crustal-scale faults by tectonic stresses. Even inherited crustal-scale discontinuities, such as sutures, cannot sufficiently reduce the strength of the lithosphere in order to be reactivated. Therefore, magmatism that had been identified to be a precursor of basement uplift in the Sierras Pampeanas, is the key factor that leads to the broken foreland of this province. Due to thermal weakening, and potentially lubrication of the inherited discontinuities, the lithosphere is locally weakened such that tectonic stresses can uplift the basement blocks. This hypothesis explains both the spatially disparate character of the broken foreland, as well as the observed temporal delay between volcanism and basement block uplift.
This dissertation provides for the first time a data-driven 3D model that is consistent with geophysical data and geological observations, and that is able to causally link the thermo-rheological structure of the lithosphere to the observed variation of surface deformation styles in the retroarc foreland of northern Argentina.
The habilitation thesis presented here includes results from several studies dealing with fluid-rock interactions and rock deformation processes in active fault zones. The focus in all of these studies is on the influence of clay minerals on the geochemical and the hydro-mechanical behavior of the fault rocks. The research was conducted on rock cores and cuttings from four scientific drilling projects at the San Andreas Fault (USA), the Nankai Trough subduction zone and the Japan Trench subduction zone (Japan), as well as the Alpine Fault in New Zealand. These ICDP (International Continental Scientific Drilling Program) and IODP (International Ocean Discovery Program) funded projects were all conducted with the aim to monitor and better understand earthquakes.
Chapter 1 contains a short introduction to the topic with basic principles and objectives regarding the research approach. Chapter 2 describes the state of the art in clay mineral and fault zone science, gives a short description of the individual drilling projects and their locations on which the research was based, and summarizes the most important analytical methods used. Chapter 3 comprises ten peer-reviewed publications that are connected thematically and methodologically. The papers were published in the years 2006-2015, and additional related publications including myself as co-author are given in the literature list. The ten publications address different questions concerning the formation of clay minerals and processes of fluid-rock interaction in active fault zones. Six papers contain results from the SAFOD drilling project, USA (San Andreas Fault Observatory at Depth), with the main focus on fluid-rock interaction processes in fault rocks and the formation and location of clay minerals. Three publications report on research from the NanTroSEIZE drilling project (Nankai Trough Seismogenic Zone Experiment) and the JFAST drilling project (Japan Trench Fast Drilling Project). Both projects are situated in Japan. Here, the swelling behavior of smectite clay minerals in relation to changing environmental conditions (e.g. temperature and/or humidity) was investigated. The last publication included here concerns a study from the DFDP project (Deep Fault Drilling Project) in New Zealand, where I investigated the deformation of clay minerals on the context of the hydro-mechanical behavior of the fault zone rocks. I was first author in nine of the publications and in charge of the project preparation, measurements and data analyses, and the completion of the manuscript. As co-author on the other publication I was responsible for electronmicroscopy analyses (SEM and TEM) and their interpretation.
The key results from the publications in Chapter 3 are discussed in Chapter 4 with additional considerations from more recent papers. Following the major theses in Chapter 5, Chapter 6 highlights a future research project in clay mineralogy research at the GFZ. An appendix includes more detailed descriptions of the laboratory equipment and lists of all publications, conference contributions and teaching courses and modules.
Partial melting is a first order process for the chemical differentiation of the crust (Vielzeuf et al., 1990). Redistribution of chemical elements during melt generation crucially influences the composition of the lower and upper crust and provides a mechanism to concentrate and transport chemical elements that may also be of economic interest. Understanding of the diverse processes and their controlling factors is therefore not only of scientific interest but also of high economic importance to cover the demand for rare metals.
The redistribution of major and trace elements during partial melting represents a central step for the understanding how granite-bound mineralization develops (Hedenquist and Lowenstern, 1994). The partial melt generation and mobilization of ore elements (e.g. Sn, W, Nb, Ta) into the melt depends on the composition of the sedimentary source and melting conditions. Distinct source rocks have different compositions reflecting their deposition and alteration histories. This specific chemical “memory” results in different mineral assemblages and melting reactions for different protolith compositions during prograde metamorphism (Brown and Fyfe, 1970; Thompson, 1982; Vielzeuf and Holloway, 1988). These factors do not only exert an important influence on the distribution of chemical elements during melt generation, they also influence the volume of melt that is produced, extraction of the melt from its source, and its ascent through the crust (Le Breton and Thompson, 1988). On a larger scale, protolith distribution and chemical alteration (weathering), prograde metamorphism with partial melting, melt extraction, and granite emplacement are ultimately depending on a (plate-)tectonic control (Romer and Kroner, 2016). Comprehension of the individual stages and their interaction is crucial in understanding how granite-related mineralization forms, thereby allowing estimation of the mineralization potential of certain areas. Partial melting also influences the isotope systematics of melt and restite. Radiogenic and stable isotopes of magmatic rocks are commonly used to trace back the source of intrusions or to quantify mixing of magmas from different sources with distinct isotopic signatures (DePaolo and Wasserburg, 1979; Lesher, 1990; Chappell, 1996). These applications are based on the fundamental requirement that the isotopic signature in the melt reflects that of the bulk source from which it is derived. Different minerals in a protolith may have isotopic compositions of radiogenic isotopes that deviate from their whole rock signature (Ayres and Harris, 1997; Knesel and Davidson, 2002). In particular, old minerals with a distinct parent-to-daughter (P/D) ratio are expected to have a specific radiogenic isotope signature. As the partial melting reaction only involves selective phases in a protolith, the isotopic signature of the melt reflects that of the minerals involved in the melting reaction and, therefore, should be different from the bulk source signature. Similar considerations hold true for stable isotopes.
When dealing with issues that are of high societal relevance, Earth sciences still face a lack of acceptance, which is partly rooted in insufficient communication strategies on the individual and local community level. To increase the efficiency of communication routines, science has to transform its outreach concepts to become more aware of individual needs and demands. The “encoding/decoding” concept as well as critical intercultural communication studies can offer pivotal approaches for this transformation.
The thesis comprises three experimental studies, which were carried out to unravel the short- as well as the long-term mechanical properties of shale rocks. Short-term mechanical properties such as compressive strength and Young’s modulus were taken from recorded stress-strain curves of constant strain rate tests. Long-term mechanical properties are represented by the time– dependent creep behavior of shales. This was obtained from constant stress experiments, where the test duration ranged from a couple minutes up to two weeks. A profound knowledge of the mechanical behavior of shales is crucial to reliably estimate the potential of a shale reservoir for an economical and sustainable extraction of hydrocarbons (HC). In addition, healing of clay-rich forming cap rocks involving creep and compaction is important for underground storage of carbon dioxide and nuclear waste.
Chapter 1 introduces general aspects of the research topic at hand and highlights the motivation for conducting this study. At present, a shift from energy recovered from conventional resources e.g., coal towards energy provided by renewable resources such as wind or water is a big challenge. Gas recovered from unconventional reservoirs (shale plays) is considered a potential bridge technology.
In Chapter 2, short-term mechanical properties of two European mature shale rocks are presented, which were determined from constant strain rate experiments performed at ambient and in situ deformation conditions (confining pressure, pc ≤ 100 MPa, temperature, T ≤ 125 °C, representing pc, T - conditions at < 4 km depth) using a Paterson– type gas deformation apparatus. The investigated shales were mainly from drill core material of Posidonia (Germany) shale and weathered material of Bowland (United Kingdom) shale. The results are compared with mechanical properties of North American shales. Triaxial compression tests performed perpendicular to bedding revealed semibrittle deformation behavior of Posidonia shale with pronounced inelastic deformation. This is in contrast to Bowland shale samples that deformed brittle and displayed predominantly elastic deformation. The static Young’s modulus, E, and triaxial compressive strength, σTCS, determined from recorded stress-strain curves strongly depended on the applied confining pressure and sample composition, whereas the influence of temperature and strain rate on E and σTCS was minor. Shales with larger amounts of weak minerals (clay, mica, total organic carbon) yielded decreasing E and σTCS. This may be related to a shift from deformation supported by a load-bearing framework of hard phases (e.g., quartz) towards deformation of interconnected weak minerals, particularly for higher fractions of about 25 – 30 vol% weak phases. Comparing mechanical properties determined at reservoir conditions with mechanical data applying effective medium theories revealed that E and σTCS of Posidonia and Bowland shale are close to the lower (Reuss) bound. Brittleness B is often quoted as a measure indicating the response of a shale formation to stimulation and economic production. The brittleness, B, of Posidonia and Bowland shale, estimated from E, is in good agreement with the experimental results. This correlation may be useful to predict B from sonic logs, from which the (dynamic) Young’s modulus can be retrieved.
Chapter 3 presents a study of the long-term creep properties of an immature Posidonia shale. Constant stress experiments (σ = const.) were performed at elevated confining pressures (pc = 50 – 200 MPa) and temperatures (T = 50 – 200 °C) to simulate reservoir pc, T - conditions. The Posidonia shale samples were acquired from a quarry in South Germany. At stresses below ≈ 84 % compressive strength of Posidonia shale, at high temperature and low confining pressure, samples showed pronounced transient (primary) creep with high deformation rates in the semibrittle regime. Sample deformation was mainly accommodated by creep of weak sample constituents and pore space reduction. An empirical power law relation between strain and time, which also accounts for the influence of pc, T and σ on creep strain was formulated to describe the primary creep phase. Extrapolation of the results to a creep period of several years, which is the typical time interval for a large production decline, suggest that fracture closure is unlikely at low stresses. At high stresses as expected for example at the contact between the fracture surfaces and proppants added during stimulation measures, subcritical crack growth may lead to secondary and tertiary creep. An empirical power law is suggested to describe secondary creep of shale rocks as a function of stress, pressure and temperature. The predicted closure rates agree with typical production decline curves recorded during the extraction of hydrocarbons. At the investigated conditions, the creep behavior of Posidonia shale was found to correlate with brittleness, calculated from sample composition.
In Chapter 4 the creep properties of mature Posidonia and Bowland shales are presented. The observed long-term creep behavior is compared to the short-term behavior determined in Chapter 2. Creep experiments were performed at simulated reservoir conditions of pc = 50 – 115 MPa and T = 75 – 150 °C. Similar to the mechanical response of immature Posidonia shale samples investigated in Chapter 3, creep strain rates of mature Bowland and Posidonia shales were enhanced with increasing stress and temperature and decreasing confining pressures. Depending on applied deformation conditions, samples displayed either only a primary (decelerating) or in addition also a secondary (quasi-steady state) and subsequently a tertiary (accelerating) creep phase before failure. At the same deformation conditions, creep strain of Posidonia shale, which is rich in weak constituents, is tremendously higher than of quartz-rich Bowland shale. Typically, primary creep strain is again mostly accommodated by deformation of weak minerals and local pore space reduction. At the onset of tertiary creep most of the deformation was accommodated by micro crack growth. A power law was used to characterize the primary creep phase of Posidonia and Bowland shale. Primary creep strain of shale rocks is inversely correlated to triaxial compressive strength and brittleness, as described in Chapter 2.
Chapter 5 provides a synthesis of the experimental findings and summarizes the major results of the studies presented in Chapters 2 – 4 and potential applications in the Exploration & Production industry.
Chapter 6 gives a brief outlook on potential future experimental research that would help to further improve our understanding of processes leading to fracture closure involving proppant embedment in unconventional shale gas reservoirs. Such insights may allow to improve stimulation techniques aimed at maintaining economical extraction of hydrocarbons over several years.
With the onset of the global food crisis, the discussion about the use and misuse of agricultural market interventions regained academic attention. As a result of economies of scale, centralized policy implementation at the regional level has the potential to reduce the budgetary costs of policies. Borrowing from the literature on international unions and international policy coordination, we develop a conceptual framework to analyze when regional policy implementation makes sense. This is the case whenever spill-overs from centralization are large and policy preferences, driven by country-specific characteristics, are homogeneous. Subsequently, we examine the advantageousness of centralized policy implementation for the West African region regarding the most common food security policies. We show that centralization of trade policies and emergency food reserves is beneficial, while buffer stocks, safety net policies, and producer support policies should be implemented at the national level.
The Atlantic Meridional Overturning Circulation (AMOC) is likely the most well-known system of ocean currents on Earth, redistributing heat, nutrients and carbon over a large part of the Earth’s surface and affecting global climate as a result. Due to enhanced freshwater fluxes into the subpolar North Atlantic as a response to global warming, the AMOC is expected, and may have already started, to weaken and these changes will likely have global impacts. It is therefore of considerable relevance to improve our understanding of past and future AMOC changes. My thesis tries to answer some of the open questions in this field by giving strong evidence that the AMOC has already weakened over the last century, by narrowing future projections of this slowdown and
by studying the impacts on global surface warming.
While there have been various studies trying to reconstruct the strength of the overturning circulation in the past, often based on model simulations in combination with observations (Jackson et al., 2016, Kanzow et al., 2010) or proxies (Frajka-Williams, 2015, Latif et al., 2006), the results so far, due to lack of direct measurements, have been inconclusive. In the first paper I build on previous work that links the anomalously low sea surface temperatures (SSTs) in the North Atlantic with the reduced meridional heat transport due to a weaker AMOC. Using the output of a high-resolution global climate model, I derive a characteristic spatial and seasonal SST fingerprint of an AMOC slowdown and an improved SST-based AMOC index. The same fingerprint is seen in
the observational SSTs since the late 19th Century, giving strong evidence that since then the AMOC has slowed down. In addition, the reconstruction of the historical overturning strength with the new AMOC index agrees well with and extends the results of earlier studies as well as the direct measurements from the RAPID project and shows a strong decline of the AMOC by about 15% (3±1 Sv) since the mid-20th Century (Caesar et al., 2018).
The reconstruction of the historical overturning strength with the AMOC index enables us to weight future AMOC projections based on their skill in modeling the historical AMOC as described in the second paper of this thesis (Olson et al., 2018). Using Bayesian model averaging we considerably narrow the projections of the CMIP5 ensemble to a decrease of -4.0 Sv and -6.8 Sv between the years 1960-1999 and 2060-2099 for the RCP4.5 and RCP8.5 emission scenarios, respectively. These values fit to, yet are at the lower end of, previously published estimates.
In the third paper I examine how the AMOC slowdown affects the global mean surface temperature (GMST) with a focus on how it will change the ocean heat uptake (OHC). Accounting for the effect of changes in the radiative forcing on the GMST, I test how AMOC variations correlate with the residual part of surface temperature changes in the past. I find that the correlation is positive which fits the understanding that the deep-water formation that is important in driving the AMOC cools the deep ocean and therefore warms the surface (Caesar et al., 2019). The future weakening of the overturning circulation could therefore delay global surface warming.
Due to nonlinear behavior and scale specific changes it can be difficult to study the dominant processes and modes that drive climate variability. In the fourth paper we develop and test a new technique based on the wavelet multiscale correlation (WMC) similarity measure to study climate variability on different temporal and spatial scales (Agarwal et al., 2018). In a fifth contribution to my thesis this method is applied to the observed sea surface temperatures. The results reconfirm well-known relations between SST anomalies such as the El Niño-Southern Oscillation (ENSO) and the Pacific Decadal Oscillation (PDO) on inter-annual and decadal timescales, respectively. They
furthermore give new insights into the characteristics and origins of long-range teleconnections, for example, that the teleconnection between ENSO and Indian Ocean dipole exist mainly between the northern part of the ENSO tongue and the equatorial Indian Ocean, and provides therefore valuable knowledge about the regions that are necessary to include when modeling regional climate variability at a certain scale (Agarwal et al., 2019).
In summary, my PhD thesis investigates past and future AMOC variability and its effects on global mean surface temperature by utilizing a combination of observational sea surface data and the output of historical and future climate model simulations from both the high-resolution CM2.6 model as well as the CMIP5 ensemble. It further includes the development and validation of a new method to study climate variability, that, applied to the observed sea surface temperatures, gives new insight about teleconnections in the Earth System. My findings provide evidence that the AMOC has already slowed down, will continue to do so in the future, and will impact the global mean temperature. Further impacts of an AMOC slowdown may include increased sea-level rise at the U.S. east coast (Ezer, 2015), heat extremes in Europe (Duchez et al., 2016) and increased storm activity in the North Atlantic region (Jackson et al., 2015), all of which have significant socio-economic implications.
The Cumbres Calchaquies Range forms part of the Famatinian metamorphic basement of the Eastern Sierras Pampeanas. The sedimentary protoliths of the metamorphic sequence were deposited in a marine basin alongside the western margin of Gondwana during the Neoproterozoic. New petrologic, geochemical and thermobarometric data give insight into the evolution of the sedimentary basin, its sediment source area, its later metamorphic overprint and its regional relationship to other parts of the Famatinian basement. The metamorphic series studied here consists of banded schists and gneisses and rare calcsilcate-rocks and migmatites that have been reworked by mid-to deep-crustal metamorphic and tectonic processes. The bulk rock compositions indicate shale, wacke, marl and litharenitic protoliths. The metamorphosed elastic sediments have major and trace element compositions indicating a continental granitoid-dominated source area with low sediment recycling. Low SiO2/Al2O3 ratios suggest a relatively low maturity of the sedimentary protoliths. Therefore, the Cumbres Calchaquies section represents a sequence of turbidity currents with progressive shallowing of the depositional environment, as indicated by quartz- and carbonate-rich sediments. The overall data are consistent with the geodynamic environment of a basin adjacent to a continental magmatic arc as the most probable scenario. Whereas the sedimentary protoliths of the metamorphic basement in the Sierra de Ancasti and Sierra de Aconquija, located ca 100-300 km south of the study area are interpreted as originating in an evolving back-arc basin, our results from the Cumbres Calchaquies region indicate a sedimentary source in a felsic continental arc with no significant influx of basic rocks. The Famatinian metamorphic evolution of the Cumbres Calchaquies rocks is of typical Barrow-type, culminating in partial melting of the metasediments. Conventional thermobarometry combined with thermodynamic models (pseudosections) reveal a prograde evolution reaching peak conditions of ca 665 degrees C/6.1 Kbar. This implies a geothermal gradient of ca 35 degrees C/km, which is slightly higher than the average for continental crust and suggests a period of crustal thinning, as known from back-arc basins, or additional heat supply by voluminous intrusions.
Although the Paris Agreement arguably made some progress, interest in supplementary approaches to climate change co-operation persist. This article examines the conditions under which a climate club might emerge and grow. Using agent-based simulations, it shows that even with less than a handful of major actors as initial members, a club can eventually reduce global emissions effectively. To succeed, a club must be initiated by the ‘right’ constellation of enthusiastic actors, offer sufficiently large incentives for reluctant countries and be reasonably unconstrained by conflicts between members over issues beyond climate change. A climate club is particularly likely to persist and grow if initiated by the United States and the European Union. The combination of club-good benefits and conditional commitments can produce broad participation under many conditions.
Climate science is highly interdisciplinary by nature, so understanding interactions between Earth processes inherently warrants the use of analytical software that can operate across the disciplines of Earth science. Toward this end, we present the Climate Data Toolbox for MATLAB, which contains more than 100 functions that span the major climate-related disciplines of Earth science. The toolbox enables streamlined, entirely scriptable workflows that are intuitive to write and easy to share. Included are functions to evaluate uncertainty, perform matrix operations, calculate climate indices, and generate common data displays. Documentation is presented pedagogically, with thorough explanations of how each function works and tutorials showing how the toolbox can be used to replicate results of published studies. As a well-tested, well-documented platform for interdisciplinary collaborations, the Climate Data Toolbox for MATLAB aims to reduce time spent writing low-level code, let researchers focus on physics rather than coding and encourage more efficacious code sharing. Plain Language Summary This article describes a collection of computer code that has recently been released to help scientists analyze many types of Earth science data. The code in this toolbox makes it easy to investigate things like global warming, El Nino, or other major climate-related processes such as how winds affect ocean circulation. Although the toolbox was designed to be used by expert climate scientists, its instruction manual is well written, and beginners may be able to learn a great deal about coding and Earth science, simply by following along with the provided examples. The toolbox is intended to help scientists save time, help them ensure their analysis is accurate, and make it easy for other scientists to repeat the results of previous studies.
The early Cretaceous Koegel Fontein intrusive complex is situated near the Atlantic coast in South Africa, about 350 km northwest of Cape Town. The complex comprises felsic units of granite and syenite with compositionally related dykes, and a single intrusive plug of diorite. Existing zircon U-Pb ages of 144 +/- 2 Ma for the syenite and 133.9 +/- 1.3 Ma for the granite suggest that the emplacement of the complex took place over a period of about 10 My. This study provides additional and independent ages of the Koegel Fontein complex by Ar-40/Ar-39 dating to confirm the onset and duration of magmatism and better define the sequence of igneous units that comprise it. New laser step-heating Ar-40/Ar-3(9) ages on plagioclase and biotite from the main intrusive units in the complex are presented here, including samples previously dated by U-Pb dating. The Ar-40/Ar-39 ages for the granite and syenite units (131.1 +/- 0.9 Ma and 143.3 +/- 0.9, respectively) are in good agreement with the zircon U-Pb ages. Other units not previously dated include the Rooivleitjie alkaline granite (150.7 +/- 0.6 Ma), two quartz-porphyry dykes (143.0 +/- 0.9 and 139.4 +/- 1.7 Ma) and the Zout Rivier diorite plug (133.0 +/- 1.0 Ma). The new results confirm an early onset of magmatism at Koegel Fontein relative to that of the Etendeka Province some 1000 km to the north, which is consistent with the regional south-to-north propagation of South Atlantic rifting. The youngest Ar-40/Ar-3(9) ages at Koegel Fontein (134 to 131 Ma, Rietpoort Granite and 133 Ma, Zout Rivier diorite) correspond to the age of the first magnetic seafloor-spreading anomaly offshore, and we suggest that the longevity of Koegel Fontein magmatism relates to a superposition of pre-drift magmatism onshore and spreading-related magmatism as continental separation began.
An essential, respected, and critical aspect of the modern practice of science and scientific publishing is peer review. The process of peer review facilitates best practices in scientific conduct and communication, ensuring that manuscripts published are as accurate, valuable, and clearly communicated. The over 216 papers published in Tectonics in 2018 benefit from the time, effort, and expertise of our reviewers who have provided thoughtfully considered advice on each manuscript. This role is critical to advancing our understanding of the evolution of the continents and their margins, as these reviews lead to even clearer and higher-quality papers. In 2018, the over 443 papers submitted to Tectonics were the beneficiaries of more than 1,010 reviews provided by 668 members of the tectonics community and related disciplines. To everyone who has volunteered their time and intellect to peer reviewing, thank you for helping Tectonics and all other AGU Publications provide the best science possible.
In the Next Generation Attenuation West2 (NGA-West2) project, a 3D subsurface structure model (Japan Seismic Hazard Information Station [J-SHIS]) was queried to establish depths to 1.0 and 2.5 km/s velocity isosurfaces for sites without depth measurement in Japan. In this article, we evaluate the depth parameters in the J-SHIS velocity model by comparing them with their corresponding site-specific depth measurements derived from selected KiK-net velocity profiles. The comparison indicates that the J-SHIS model underestimates site depths at shallow sites and overestimates depths at deep sites. Similar issues were also identified in the southern California basin model. Our results also show that these underestimations and over-estimations have a potentially significant impact on ground-motion prediction using NGA-West2 ground-motion models (GMMs). Site resonant period may be considered as an alternative to depth parameter in the site term of a GMM.
The morphology of marine and lacustrine terraces has been largely used to measure past sea- and lake-level positions and estimate vertical deformation in a wealth of studies focused on climate and tectonic processes. To obtain accurate morphometric assessments of terrace morphology we present TerraceM-2, an improved version of our MatlabR (R) graphic-user interface that provides new methodologies for morphometric analyses as well as landscape evolution and fault-dislocation modeling. The new version includes novel routines to map the elevation and spatial distribution of terraces, to model their formation and evolution, and to estimate fault-slip rates from terrace deformation patterns. TerraceM-2 has significantly improves its processing speed and mapping capabilities, and includes separate functions for developing customized workflows beyond the graphic-user interface. We illustrate these new mapping and modeling capabilities with three examples: mapping lacustrine shorelines in the Dead Sea to estimate deformation across the Dead Sea Fault, landscape evolution modeling to estimate a history of uplift rates in southern Peru, and dislocation modeling of deformed marine terraces in California. These examples also illustrate the need to use topographic data of different resolutions. The new modeling and mapping routines of TerraceM-2 highlight the advantages of an integrated joint mapping and modeling approach to improve the efficiency and precision of coastal terrace metrics in both marine and lacustrine environments.
River bank filtration (RBF) is considered to efficiently remove nitrate and trace organic micropollutants (OMP) from polluted surface waters. This is essential for maintaining good groundwater quality and providing high quality drinking water. Predicting the fate of OMP during RBF is difficult as the biogeochemical factors controlling the removal efficiency are not fully understood. To determine in-situ removal efficiency and degradation rates of nitrate and OMP indicator substances we conducted a field study in a RBF system during a period of one and a half years incorporating temporally and spatially varying redox conditions and temperature changes typically occurring in temperate climates. RBF was analyzed by means of mixing ratios between infiltrated river water and groundwater as well as average residence times of surface water towards the individual groundwater observation wells. These results were used to calculate temperature dependent first order degradation rates of redox sensitive species and several OMP. Five out of ten investigated OMP were completely removed along RBF pathways. We demonstrate that degradation rates of several OMP during bank filtration were controlled by redox conditions and temperature whereby temperature itself also had a significant influence on the extent of the most reactive oxic zone. The seasonal variations in temperature alone could explain a considerable percentage of the variance in dissolved oxygen (34%), nitrate (81%) as well as the OMPs diclofenac (44%) and sulfamethoxazole (76%). Estimated in-situ degradation rates roughly varied within one order of magnitude for temperature changes between 5 degrees C and 20 degrees C. This study highlights that temporal variability in temperature and redox zonation is a significant factor for migration and degradation of nitrate and several OMPs. (C) 2019 Elsevier Ltd. All rights reserved.
The northward indentation of the Pamir salient into the Tarim basin at the western syntaxis of the India-Asia collision zone is the focus of controversial models linking lithospheric to surface and atmospheric processes. Here we report on tectonic events recorded in the most complete and best-dated sedimentary sequences from the western Tarim basin flanking the eastern Pamir (the Aertashi section), based on sedimentologic, provenance, and magnetostratigraphic analyses. Increased tectonic subsidence and a shift from marine to continental fluvio-deltaic deposition at 41Ma indicate that far-field deformation from the south started to affect the Tarim region. A sediment accumulation hiatus from 24.3 to 21.6Ma followed by deposition of proximal conglomerates is linked to fault propagation into the Tarim basin. From 21.6 to 15.0Ma, increasing accumulation rates of fining upward clastics is interpreted as the expression of a major dextral transtensional system linking the Kunlun to the Tian Shan ahead of the northward Pamir indentation. At 15.0Ma, the appearance of North Pamir-sourced conglomerates followed at 11Ma by Central Pamir-sourced volcanics coincides with a shift to E-W compression, clockwise vertical-axis rotations and the onset of growth strata associated with the activation of the local east vergent Qimugen thrust wedge. Together, this enables us to interpret that Pamir indentation into Tarim had started by 24.3Ma, reached the study location by 15.0Ma and had passed it by 11Ma, providing kinematic constraints on proposed tectonic models involving intracontinental subduction and delamination.
Underground coal gasification (UCG) enables utilization of coal reserves, currently not economically exploitable due to complex geological boundary conditions. Hereby, UCG produces a high-calorific synthesis gas that can be used for generation of electricity, fuels, and chemical feedstock. The present study aims to identify economically-competitive, site-specific end-use options for onshore- and offshore-produced UCG synthesis gas, taking into account the capture and storage (CCS) and/or utilization (CCU) of produced CO2. Modeling results show that boundary conditions favoring electricity, methanol, and ammonia production expose low costs for air separation, low compression power requirements, and appropriate shares of H-2/N-2. Hereby, a gasification agent ratio of more than 30% oxygen by volume is not favorable from the economic and CO2 mitigation viewpoints. Compared to the costs of an offshore platform with its technical equipment, offshore drilling costs are marginal. Thus, uncertainties related to parameters influenced by drilling costs are negligible. In summary, techno-economic process modeling results reveal that air-blown gasification scenarios are the most cost-effective ones, while offshore UCG-CCS/CCU scenarios are up to 1.7 times more expensive than the related onshore processes. Hereby, all investigated onshore scenarios except from ammonia production under the assumed worst-case conditions are competitive on the European market.
The Kp index is a measure of the midlatitude global geomagnetic activity and represents short-term magnetic variations driven by solar wind plasma and interplanetary magnetic field. The Kp index is one of the most widely used indicators for space weather alerts and serves as input to various models, such as for the thermosphere and the radiation belts. It is therefore crucial to predict the Kp index accurately. Previous work in this area has mostly employed artificial neural networks to nowcast Kp, based their inferences on the recent history of Kp and on solar wind measurements at L1. In this study, we systematically test how different machine learning techniques perform on the task of nowcasting and forecasting Kp for prediction horizons of up to 12 hr. Additionally, we investigate different methods of machine learning and information theory for selecting the optimal inputs to a predictive model. We illustrate how these methods can be applied to select the most important inputs to a predictive model of Kp and to significantly reduce input dimensionality. We compare our best performing models based on a reduced set of optimal inputs with the existing models of Kp, using different test intervals, and show how this selection can affect model performance.
The heavy rainfall events in recent years have caused great damage, which has increased the public awareness of the topic of heavy rainfall. For this reason, this article discusses how a systematic integration of heavy rainfall within the framework of the European Floods Directive would be possible and reasonable. For this purpose, a matrix covering possible synergies and barriers was created for all steps of the directive, which were then examined in 15 semi-structured interviews with representatives from specialized administration, the private sector and academia. Although there are some synergies, the additional effort required, especially regarding the identification of the risk areas and the higher level of detail required for risk modeling, would be so high that the European Floods Directive cannot be deemed to be an appropriate framework for heavy rainfall risk management. Nevertheless, there is a need for action, e.g. in the field of self-protection, improved risk communication to the population, combined with increased public and interagency cooperation.
The fall into the Oligocene icehouse is marked by a steady decline in global temperature with punctuated cooling at the Eocene-Oligocene transition, both of which are well documented in the marine realm. However, the chronology and mechanisms of cooling on land remain unclear. Here, we use clumped isotope thermometry on northeastern Tibetan continental carbonates to reconstruct a detailed Paleogene surface temperature record for the Asian continental interior, and correlate this to an enhanced pollen data set. Our results show two successive dramatic (>9 degrees C) temperature drops, at 37 Ma and at 33.5 Ma. These large-magnitude decreases in continental temperatures can only be explained by a combination of both regional cooling and shifts of the rainy season to cooler months, which we interpret to reflect a decline of monsoonal intensity. Our results suggest that the response of Asian surface temperatures and monsoonal rainfall to the steady decline of atmospheric CO2 and global temperature through the late Eocene was nonlinear and occurred in two steps separated by a period of climatic instability. Our results support the onset of the Antarctic Circumpolar Current coeval to the Oligocene isotope event 1 (Oi-1) glaciation at 33.5 Ma, reshaping the distribution of surface heat worldwide; however, the origin of the 37 Ma cooling event remains less clear.
We address the question of whether all large-magnitude earthquakes produce an erosion peak in the subaerial components of fluvial catchments. We evaluate the sediment flux response to the Maule earthquake in the Chilean Andes (Mw 8.8) using daily suspended sediment records from 31 river gauges. The catchments cover drainage areas of 350 to around 10,000 km(2), including a wide range of topographic slopes and vegetation cover of the Andean western flank. We compare the 3- to 8-year postseismic record of sediment flux to each of the following preseismic periods: (1) all preseismic data, (2) a 3-year period prior to the seismic event, and (3) the driest preseismic periods, as drought conditions prevailed in the postseismic period. Following the earthquake, no increases in suspended sediment flux were observed for moderate to high percentiles of the streamflow distribution (mean, median, and >= 75th percentile). However, more than half of the examined stations showed increased sediment flux during baseflow. By using a Random Forest approach, we evaluate the contributions of seismic intensities, peak ground accelerations, co-seismic landslides, hydroclimatic conditions, topography, lithology, and land cover to explain the observed changes in suspended sediment concentration and fluxes. We find that the best predictors are hillslope gradient, low-vegetation cover, and changes in streamflow discharge. This finding suggests a combined first-order control of topography, land cover, and hydrology on the catchment-wide erosion response. We infer a reduced sediment connectivity due to the postseismic drought, which increased the residence time of sediment detached and remobilized following the Maule earthquake.
Plate tectonics is among the most important geological processes on Earth, but its emergence and evolution remain unclear. Here we extrapolate models of present-day plate tectonics to the past and propose that since about three billion years ago the rise of continents and the accumulation of sediments at continental edges and in trenches has provided lubrication for the stabilization of subduction and has been crucial in the development of plate tectonics on Earth. We conclude that the two largest surface erosion and subduction lubrication events occurred after the Palaeoproterozoic Huronian global glaciations (2.45 to 2.2 billion years ago), leading to the formation of the Columbia supercontinent, and after the Neoproterozoic ‘snowball’ Earth glaciations (0.75 to 0.63 billion years ago). The snowball Earth event followed the ‘boring billion’—a period of reduced plate tectonic activity about 1.75 to 0.75 billion years ago that was probably caused by a shortfall of sediments in trenches—and it kick-started the modern episode of active plate tectonics.
During lower sea levels in glacial periods, deep permafrost formed on large continental shelf areas of the Arctic Ocean. Subsequent sea level rise and coastal erosion created subsea permafrost, which generally degrades after inundation under the influence of a complex suite of marine, near-shore processes. Global warming is especially pronounced in the Arctic, and will increase the transition to and the degradation of subsea permafrost, with implications for atmospheric climate forcing, offshore infrastructure, and aquatic ecosystems.
This thesis combines new geophysical, borehole observational and modelling approaches to enhance our understanding of subsea permafrost dynamics. Three specific areas for advancement were identified: (I) sparsity of observational data, (II) lacking implementation of salt infiltration mechanisms in models, and (III) poor understanding of the regional differences in key driving parameters. This study tested the combination of spectral ratios of the ambient vibration seismic wavefield, together with estimated shear wave velocity from seismic interferometry analysis, for estimating the thickness of the unfrozen sediment overlying the ice-bonded permafrost offshore. Mesoscale numerical calculations (10^1 to 10^2 m, thousands of years) were employed to develop and solve the coupled heat diffusion and salt transport equations including phase change effects. Model soil parameters were constrained by borehole data, and the impact of a variety of influences during the transgression was tested in modelling studies. In addition, two inversion schemes (particle swarm optimization and a least-square method) were used to reconstruct temperature histories for the past 200-300 years in the Laptev Sea region in Siberia from two permafrost borehole temperature records. These data were evaluated against larger scale reconstructions from the region.
It was found (I) that peaks in spectral ratios modelled for three-layer, one-dimensional systems corresponded with thaw depths. Around Muostakh Island in the central Laptev Sea seismic receivers were deployed on the seabed. Derived depths of the ice-bonded permafrost table were between 3.7-20.7 m ± 15 %, increasing with distance from the coast. (II) Temperatures modelled during the transition to subsea permafrost resembled isothermal conditions after about 2000 years of inundation at Cape Mamontov Klyk, consistent with observations from offshore boreholes. Stratigraphic scenarios showed that salt distribution and infiltration had a large impact on the ice saturation in the sediments. Three key factors were identified that, when changed, shifted the modelled permafrost thaw depth most strongly: bottom water temperatures, shoreline retreat rate and initial temperature before inundation. Salt transport based on diffusion and contribution from arbitrary density-driven mechanisms only accounted for about 50 % of observed thaw depths at offshore sites hundreds to thousands of years after inundation. This bias was found consistently at all three sites in the Laptev Sea region. (III) In the temperature reconstructions, distinct differences in the local temperature histories between the western Laptev Sea and the Lena Delta sites were recognized, such as a transition to warmer temperatures a century later in the western Laptev Sea as well as a peak in warming three decades later. The local permafrost surface temperature history at Sardakh Island in the Lena Delta was reminiscent of the circum-Arctic regional average trends. However, Mamontov Klyk in the western Laptev Sea was consistent to Arctic trends only in the most recent decade and was more similar to northern hemispheric mean trends. Both sites were consistent with a rapid synoptic recent warming.
In conclusion, the consistency between modelled response, expected permafrost distribution, and observational data suggests that the passive seismic method is promising for the determination of the thickness of unfrozen sediment on the continental Arctic shelf. The quantified gap between currently modelled and observed thaw depths means that the impact of degradation on climate forcing, ecosystems, and infrastructure is larger than current models predict. This discrepancy suggests the importance of further mechanisms of salt penetration and thaw that have not been considered – either pre-inundation or post-inundation, or both. In addition, any meaningful modelling of subsea permafrost would have to constrain the identified key factors and their regional differences well. The shallow permafrost boreholes provide missing well-resolved short-scale temperature information in the coastal permafrost tundra of the Arctic. As local differences from circum-Arctic reconstructions, such as later warming and higher warming magnitude, were shown to exist in this region, these results provide a basis for local surface temperature record parameterization of climate and, in particular, permafrost models. The results of this work bring us one step further to understanding the full picture of the transition from terrestrial to subsea permafrost.
Carbonate-rich silicate and carbonate melts play a crucial role in deep Earth magmatic processes and their melt structure is a key parameter, as it controls physical and chemical properties. Carbonate-rich melts can be strongly enriched in geochemically important trace elements. The structural incorporation mechanisms of these elements are difficult to study because such melts generally cannot be quenched to glasses, which are usually employed for structural investigations. This thesis investigates the influence of CO2 on the local environments of trace elements contained in silicate glasses with variable CO2 concentrations as well as in silicate and carbonate melts. The compositions studied include sodium-rich peralkaline silicate melts and glasses and carbonate melts similar to those occurring naturally at Oldoinyo Lengai volcano, Tanzania.
The local environments of the three elements yttrium (Y), lanthanum (La) and strontium (Sr) were investigated in synthesized glasses and melts using X-ray absorption fine structure (XAFS) spectroscopy. Especially extended X-ray absorption fine structure spectroscopy (EXAFS) provides element specific information on local structure, such as bond lengths, coordination numbers and the degree of disorder. To cope with the enhanced structural disorder present in glasses and melts, EXAFS analysis was based on fitting approaches using an asymmetric distribution function as well as a correlation model according to bond valence theory. Firstly, silicate glasses quenched from high pressure/temperature melts with up to 7.6 wt % CO2 were investigated. In strongly and extremely peralkaline glasses the local structure of Y is unaffected by the CO2 content (with oxygen bond lengths of ~ 2.29 Å). Contrary, the bond lengths for Sr-O and La-O increase with increasing CO2 content in the strongly peralkaline glasses from ~ 2.53 to ~ 2.57 Å and from ~ 2.52 to ~ 2.54 Å, respectively, while they remain constant in extremely peralkaline glasses (at ~ 2.55 Å and 2.54 Å, respectively). Furthermore, silicate and unquenchable carbonate melts were investigated in-situ at high pressure/temperature conditions (2.2 to 2.6 GPa, 1200 to 1500 °C) using a Paris-Edinburgh press. A novel design of the pressure medium assembly for this press was developed, which features increased mechanical stability as well as enhanced transmittance at relevant energies to allow for low content element EXAFS in transmission. Compared to glasses the bond lengths of Y-O, La-O and Sr-O are elongated by up to + 3 % in the melt and exhibit higher asymmetric pair distributions. For all investigated silicate melt compositions Y-O bond lengths were found constant at ~ 2.37 Å, while in the carbonate melt the Y-O length increases slightly to 2.41 Å. The La-O bond lengths in turn, increase systematically over the whole silicate – carbonate melt joint from 2.55 to 2.60 Å. Sr-O bond lengths in melts increase from ~ 2.60 to 2.64 Å from pure silicate to silicate-bearing carbonate composition with constant elevated bond length within the carbonate region.
For comparison and deeper insight, glass and melt structures of Y and Sr bearing sodium-rich silicate to carbonate compositions were simulated in an explorative ab initio molecular dynamics (MD) study. The simulations confirm observed patterns of CO2-dependent local changes around Y and Sr and additionally provide further insights into detailed incorporation mechanisms of the trace elements and CO2. Principle findings include that in sodium-rich silicate compositions carbon either is mainly incorporated as a free carbonate-group or shares one oxygen with a network former (Si or [4]Al) to form a non-bridging carbonate. Of minor importance are bridging carbonates between two network formers. Here, a clear preference for two [4]Al as adjacent network formers occurs, compared to what a statistical distribution would suggest. In C-bearing silicate melts minor amounts of molecular CO2 are present, which is almost totally dissolved as carbonate in the quenched glasses.
The combination of experiment and simulation provides extraordinary insights into glass and melt structures. The new data is interpreted on the basis of bond valence theory and is used to deduce potential mechanisms for structural incorporation of investigated elements, which allow for prediction on their partitioning behavior in natural melts. Furthermore, it provides unique insights into the dissolution mechanisms of CO2 in silicate melts and into the carbonate melt structure. For the latter, a structural model is suggested, which is based on planar CO3-groups linking 7- to 9-fold cation polyhedra, in accordance to structural units as found in the Na-Ca carbonate nyerereite. Ultimately, the outcome of this study contributes to rationalize the unique physical properties and geological phenomena related to carbonated silicate-carbonate melts.
Fold and thrust belts are characteristic features of collisional orogen that grow laterally through time by deforming the upper crust in response to stresses caused by convergence. The deformation propagation in the upper crust is accommodated by shortening along major folds and thrusts. The formation of these structures is influenced by the mechanical strength of décollements, basement architecture, presence of preexisting structures and taper of the wedge. These factors control not only the sequence of deformation but also cause differences in the structural style.
The Himalayan fold and thrust belt exhibits significant differences in the structural style from east to west. The external zone of the Himalayan fold and thrust belt, also called the Subhimalaya, has been extensively studied to understand the temporal development and differences in the structural style in Bhutan, Nepal and India; however, the Subhimalaya in Pakistan remains poorly studied. The Kohat and Potwar fold and thrust belts (herein called Kohat and Potwar) represent the Subhimalaya in Pakistan. The Main Boundary Thrust (MBT) marks the northern boundary of both Kohat and Potwar, showing that these belts are genetically linked to foreland-vergent deformation within the Himalayan orogen, despite the pronounced contrast in structural style. This contrast becomes more pronounced toward south, where the active strike-slip Kalabagh Fault Zone links with the Kohat and Potwar range fronts, known as the Surghar Range and the Salt Range, respectively. The Surghar and Salt Ranges developed above the Surghar Thrust (SGT) and Main Frontal Thrust (MFT). In order to understand the structural style and spatiotemporal development of the major structures in Kohat and Potwar, I have used structural modeling and low temperature thermochronolgy methods in this study. The structural modeling is based on construction of balanced cross-sections by integrating surface geology, seismic reflection profiles and well data. In order to constrain the timing and magnitude of exhumation, I used apatite (U-Th-Sm)/He (AHe) and apatite fission track (AFT) dating. The results obtained from both methods are combined to document the Paleozoic to Recent history of Kohat and Potwar.
The results of this research suggest two major events in the deformation history. The first major deformation event is related to Late Paleozoic rifting associated with the development of the Neo-Tethys Ocean. The second major deformation event is related to the Late Miocene to Pliocene development of the Himalayan fold and thrust belt in the Kohat and Potwar. The Late Paleozoic rifting is deciphered by inverse thermal modelling of detrital AFT and AHe ages from the Salt Range. The process of rifting in this area created normal faulting that resulted in the exhumation/erosion of Early to Middle Paleozoic strata, forming a major unconformity between Cambrian and Permian strata that is exposed today in the Salt Range. The normal faults formed in Late Paleozoic time played an important role in localizing the Miocene-Pliocene deformation in this area. The combination of structural reconstructions and thermochronologic data suggest that deformation initiated at 15±2 Ma on the SGT ramp in the southern part of Kohat. The early movement on the SGT accreted the foreland into the Kohat deforming wedge, forming the range front. The development of the MBT at 12±2 Ma formed the northern boundary of Kohat and Potwar. Deformation propagated south of the MBT in the Kohat on double décollements and in the Potwar on a single basal décollement. The double décollement in the Kohat adopted an active roof-thrust deformation style that resulted in the disharmonic structural style in the upper and lower parts of the stratigraphic section. Incremental shortening resulted in the development of duplexes in the subsurface between two décollements and imbrication above the roof thrust. Tectonic thickening caused by duplexes resulted in cooling and exhumation above the roof thrust by removal of a thick sequence of molasse strata. The structural modelling shows that the ramps on which duplexes formed in Kohat continue as tip lines of fault propagation folds in the Potwar. The absence of a double décollement in the Potwar resulted in the preservation of a thick sequence of molasse strata there. The temporal data suggest that deformation propagated in-sequence from ~ 8 to 3 Ma in the northern part of Kohat and Potwar; however, internal deformation in the Kohat was more intense, probably required for maintaining a critical taper after a significant load was removed above the upper décollement. In the southern part of Potwar, a steeper basement slope (β≥3°) and the presence of salt at the base of the stratigraphic section allowed for the complete preservation of the stratigraphic wedge, showcased by very little internal deformation. Activation of the MFT at ~4 Ma allowed the Salt Range to become the range front of the Potwar. The removal of a large amount of molasse strata above the MFT ramp enhanced the role of salt in shaping the structural style of the Salt Range and Kalabagh Fault Zone. Salt accumulation and migration resulted in the formation of normal faults in both areas. Salt migration in the Kalabagh fault zone has triggered out-of-sequence movement on ramps in the Kohat.
The amount of shortening calculated between the MBT and the SGT in Kohat is 75±5 km and between the MBT and the MFT in Potwar is 65±5 km. A comparable amount of shortening is accommodated in the Kohat and Potwar despite their different widths: 70 km Kohat and 150 km Potwar. In summary, this research suggests that deformation switched between different structures during the last ~15 Ma through different modes of fault propagation, resulting in different structural styles and the out-of-sequence development of Kohat and Potwar.
The task of downloading comprehensive datasets of event-based seismic waveforms has been made easier through the development of standardized webservices but is still highly nontrivial because the likelihood of temporary network failures or subtle data errors naturally increases when the amount of requested data is in the order of millions of relatively short segments. This is even more challenging because the typical workflow is not restricted to a single massive download but consists of fetching all possible available input data (e.g., with several repeated download executions) for a processing stage producing any desired user-defined output. Here, we present stream2segment, a highly customizable Python 2+3 package helping the user in the entire workflow of downloading, inspecting, and processing event-based seismic data by means of a relational database management system as archiving storage, which has clear performance and usability advantages, and an integrated processing subroutine requiring a configuration file and a single Python function to produce user-defined output. Stream2segment can also produce diagnostic maps or user-defined plots, which, unlike existing tools, do not require external software dependencies and are not static images but instead are interactive browser-based applications ideally suited for data inspection or annotation tasks and subsequent training of classifiers in foreseen supervised machine-learning applications. Stream2segment has already been used as a data quality tool for datasets within the European Integrated Data Archive and to create a weak-motion database (in the form of a so-called flat file) for the stable continental region of Europe in the context of the European Ground Shaking Intensity Model service, in turn an important building block for seismic hazard studies.
Localization processes in the viscous lower crust generate ductile shear zones over a broad range of scales affecting long‐term lithosphere deformation and the mechanical response of faults during the seismic cycle. Here we use centimeter‐scale numerical models in order to gain detailed insight into the processes involved in strain localization and rheological weakening in viscously deforming rocks. Our 2‐D Cartesian models are benchmarked to high‐temperature and high‐pressure torsion experiments on Carrara marble samples containing a single weak Solnhofen limestone inclusion. The models successfully reproduce bulk stress‐strain transients and final strain distributions observed in the experiments by applying a simple, first‐order softening law that mimics rheological weakening. We find that local stress concentrations forming at the inclusion tips initiate strain localization inside the host matrix. At the tip of the propagating shear zone, weakening occurs within a process zone, which expands with time from the inclusion tips toward the matrix. Rheological weakening is a precondition for shear zone localization, and the width of this shear zone is found to be controlled by the degree of softening. Introducing a second softening step at elevated strain, a high strain layer develops inside the localized shear zone, analogous to the formation of ultramylonite bands in mylonites. These results elucidate the transient evolution of stress and strain rate during inception and maturation of ductile shear zones.
Lake sediments are increasingly explored as reliable paleoflood archives. In addition to established flood proxies including detrital layer thickness, chemical composition, and grain size, we explore stable oxygen and carbon isotope data as paleoflood proxies for lakes in catchments with carbonate bedrock geology. In a case study from Lake Mondsee (Austria), we integrate high-resolution sediment trapping at a proximal and a distal location and stable isotope analyses of varved lake sediments to investigate flood-triggered detrital sediment flux. First, we demonstrate a relation between runoff, detrital sediment flux, and isotope values in the sediment trap record covering the period 2011-2013 CE including 22 events with daily (hourly) peak runoff ranging from 10 (24) m(3) s(-1) to 79 (110) m(3) s(-1). The three- to ten-fold lower flood-triggered detrital sediment deposition in the distal trap is well reflected by attenuated peaks in the stable isotope values of trapped sediments. Next, we show that all nine flood-triggered detrital layers deposited in a sediment record from 1988 to 2013 have elevated isotope values compared with endogenic calcite. In addition, even two runoff events that did not cause the deposition of visible detrital layers are distinguished by higher isotope values. Empirical thresholds in the isotope data allow estimation of magnitudes of the majority of floods, although in some cases flood magnitudes are overestimated because local effects can result in too-high isotope values. Hence we present a proof of concept for stable isotopes as reliable tool for reconstructing flood frequency and, although with some limitations, even for flood magnitudes.
Lake sediments are increasingly explored as reliable paleoflood archives. In addition to established flood proxies including detrital layer thickness, chemical composition, and grain size, we explore stable oxygen and carbon isotope data as paleoflood proxies for lakes in catchments with carbonate bedrock geology. In a case study from Lake Mondsee (Austria), we integrate high-resolution sediment trapping at a proximal and a distal location and stable isotope analyses of varved lake sediments to investigate flood-triggered detrital sediment flux. First, we demonstrate a relation between runoff, detrital sediment flux, and isotope values in the sediment trap record covering the period 2011-2013 CE including 22 events with daily (hourly) peak runoff ranging from 10 (24) m(3) s(-1) to 79 (110) m(3) s(-1). The three- to ten-fold lower flood-triggered detrital sediment deposition in the distal trap is well reflected by attenuated peaks in the stable isotope values of trapped sediments. Next, we show that all nine flood-triggered detrital layers deposited in a sediment record from 1988 to 2013 have elevated isotope values compared with endogenic calcite. In addition, even two runoff events that did not cause the deposition of visible detrital layers are distinguished by higher isotope values. Empirical thresholds in the isotope data allow estimation of magnitudes of the majority of floods, although in some cases flood magnitudes are overestimated because local effects can result in too-high isotope values. Hence we present a proof of concept for stable isotopes as reliable tool for reconstructing flood frequency and, although with some limitations, even for flood magnitudes.
Using optical absorption and Raman spectroscopic measurements, in conjunction with the first-principles calculations, a pressure-induced high-spin (HS)-to-low-spin (LS) state electronic transition of Fe2+ (M2-octahedral site) was resolved around 76-80GPa in a natural triphylite-lithiophilite sample with chemical composition (LiFe0.7082+Mn0.292PO4)-Li-M1-Fe-M2 (theoretical composition (LiFe0.52+Mn0.5PO4)-Li-M1-Fe-M2). The optical absorption spectra at ambient conditions consist of a broad doublet band with two constituents (1) (similar to 9330cm(-1)) and (2) (similar to 7110cm(-1)), resulting from the electronic spin-allowed transition (T2gEg)-T-5-E-5 of octahedral (HSFe2+)-Fe-M2. Both (1) and (2) bands shift non-linearly with pressure to higher energies up to similar to 55GPa. In the optical absorption spectrum measured at similar to 81GPa, the aforementioned HS-related bands disappear, whereas a new broadband with an intensity maximum close to 16,360cm(-1) appears, superimposed on the tail of the high-energy ligand-to-metal O2-Fe2+ charge-transfer absorption edge. We assign this new band to the electronic spin-allowed dd-transition (1)A(1g)(1)T(1g) of LS Fe2+ in octahedral coordination. The high-pressure Raman spectra evidence the Fe2+ HS-to-LS transition mainly from the abrupt shift of the P-O symmetric stretching modes to lower frequencies at similar to 76GPa, the highest pressure achieved in the Raman spectroscopic experiments. Calculations indicated that the presence of Mn-M2(2+) simply shifts the isostructural HS-to-LS transition to higher pressures compared to the triphylite Fe-M2(2+) end-member, in qualitative agreement with our experimental observations.
Rapidly changing climate in the Northern Hemisphere and associated socio-economic impacts require reliable understanding of lake systems as important freshwater resources and sensitive sentinels of environmental change. To better understand time-series data in lake sediment cores, it is necessary to gain information on within-lake spatial variabilities of environmental indicator data. Therefore, we retrieved a set of 38 samples from the sediment surface along spatial habitat gradients in the boreal, deep, and yet pristine Lake Bolshoe Toko in southern Yakutia, Russia. Our methods comprise laboratory analyses of the sediments for multiple proxy parameters, including diatom and chironomid taxonomy, oxygen isotopes from diatom silica, grain-size distributions, elemental compositions (XRF), organic carbon content, and mineralogy (XRD). We analysed the lake water for cations, anions, and isotopes. Our results show that the diatom assemblages are strongly influenced by water depth and dominated by planktonic species, i.e. Pliocaenicus bolshetokoensis. Species richness and diversity are higher in the northern part of the lake basin, associated with the availability of benthic, i.e. periphytic, niches in shallower waters. delta O-18(diatom) values are higher in the deeper south-western part of the lake, probably related to water temperature differences. The highest amount of the chironomid taxa underrepresented in the training set used for palaeoclimate inference was found close to the Utuk River and at southern littoral and profundal sites. Abiotic sediment components are not symmetrically distributed in the lake basin, but vary along restricted areas of differential environmental forcing. Grain size and organic matter are mainly controlled by both river input and water depth. Mineral (XRD) data distributions are influenced by the methamorphic lithology of the Stanovoy mountain range, while elements (XRF) are intermingled due to catchment and diagenetic differences. We conclude that the lake represents a valuable archive for multiproxy environmental reconstruction based on diatoms (including oxygen isotopes), chironomids, and sediment-geochemical parameters. Our analyses suggest multiple coring locations preferably at intermediate depth in the northern basin and the deep part in the central basin, to account for representative bioindicator distributions and higher temporal resolution, respectively.
Identifying the provenance signature and geodynamic setting on which sedimentary basins at convergent margins grow is challenging since they result from coupled erosional and tectonic processes, which shape the evolution of source areas and the stress regime. The Early Cretaceous evolution of the northern Andes of Colombia is characterized by extensional tectonics and the subsequent formation of a marginal basin. The Abejorral Formation and coeval volcano-sedimentary rocks are exposed along the western flank and axis of the Central Cordillera. They comprise an Early Cretaceous transgressive sequence initially accumulated in fluvial deltaic environments, which switched towards a deep-marine setting, and are interpreted as the infilling record of a marginal back-arc basin. Available provenance data suggest that Permian-Triassic metamorphic and less abundant Jurassic magmatic rocks forming the basement of the Central Cordillera sourced the Abejorral Formation. New detailed volcanic and metamorphic lithics analyses, conventional and varietal study of heavy minerals, detrital rutile mineral chemistry, allowed us to document changes in the source areas defined by the progressive appearance of both higher-grade and more distal low-grade metamorphic sources, which switched from pelitic to dominantly mafic in composition. Crystallochemical indexes of clay minerals of fine-grained rocks of the Abejorral Formation suggest that samples located close to the Romeral Fault System show characteristics of low-medium P-T low-grade metamorphism, whereas rocks located farther to the northeast preserve primary diagenetic features, which suggest a high heat-flow accumulation setting. We interpret that the Abejorral Formation records the progressive unroofing of the Central Cordillera basement that was being rapidly exhumed, as well as the incorporation of distal subduction-related metamorphic complexes to the west in response either to the widening of extensional front or the reactivation of fault structures on the oceanward margin of the basin. Although the deformational record of the Abejorral Formation would have resulted from over-imposed episodes, our new geochronological constraints suggest that this sedimentary sequence must have been deformed before the Paleocene due to the presence of arc-related intrusive non-deformed magmatic rocks with a crystallization age of ca. 60 Ma.
Sociocultural valuation (SCV) of ecosystem services (ES) discloses the principles, importance or preferences expressed by people towards nature. Although ES research has increasingly addressed sociocultural values in past years, little effort has been made to systematically review the components of sociocultural valuation applications for different decision contexts (i.e. awareness raising, accounting, priority setting, litigation and instrument design). In this analysis, we investigate the characteristics of 48 different sociocultural valuation applications—characterised by unique combinations of decision context, methods, data collection formats and participants—across ten European case studies. Our findings show that raising awareness for the sociocultural value of ES by capturing people’s perspective and establishing the status quo, was found the most frequent decision context in case studies, followed by priority setting and instrument development. Accounting and litigation issues were not addressed in any of the applications. We reveal that applications for particular decision contexts are methodologically similar, and that decision contexts determine the choice of methods, data collection formats and participants involved. Therefore, we conclude that understanding the decision context is a critical first step to designing and carrying out fit-for-purpose sociocultural valuation of ES in operational ecosystem management.
Small-aperture array as a tool to monitor fluid injection- and extraction-induced microseismicity
(2019)
The monitoring of microseismicity during temporary human activities such as fluid injections for hydrofracturing, hydrothermal stimulations or wastewater disposal is a difficult task. The seismic stations often cannot be installed on hard rock, and at quiet places, noise is strongly increased during the operation itself and the installation of sensors in deep wells is costly and often not feasible. The combination of small-aperture seismic arrays with shallow borehole sensors offers a solution. We tested this monitoring approach at two different sites: (1) accompanying a fracking experiment in sedimentary shale at 4km depth and (2) above a gas field under depletion. The small-aperture arrays were planned according to theoretical wavenumber studies combined with simulations considering the local noise conditions. We compared array recordings with recordings available from shallow borehole sensors and give examples of detection and location performance. Although the high-frequency noise on the 50-m-deep borehole sensors was smaller compared to the surface noise before the injection experiment, the signals were highly contaminated during injection by the pumping activities. Therefore, a set of three small-aperture arrays at different azimuths was more suited to detect small events, since noise recorded on these arrays is uncorrelated with each other. Further, we developed recommendations for the adaptation of the monitoring concept to other sites experiencing induced seismicity.
We analyze a large transient strainmeter signal recorded at 62.5 m depth along the southern shore of the eastern Sea of Marmara region in northwestern Turkey. This region represents a passage of stress transfer from the Izmit rupture to the Marmara seismic gap. The strain signal was recorded at the Esenkoy site by one of the ICDP-GONAF (International Continental Drilling Programme - Geophysical Observatory at the North Anatolian Fault) strainmeters on the Armutlu peninsula with a maximum amplitude of 5 microstrain and lasting about 50 days. The onset of the strain signal coincided with the origin time of a M-w 4.4 earthquake offshore Yalova, which occurred as part of a seismic sequence including eight M-w >= 3.5 earthquakes. The Mw 4.4 event occurred at a distance of about 30 km from Esenkoy on June 25th 2016 representing the largest earthquake in this region since 2008. Before the event, the maximum horizontal strain was subparallel to the regional maximum horizontal stress derived from stress inversion of local seismicity. During the strain transient, we observe a clockwise rotation in the local horizontal strain field of about 20 degrees. The strain signal does not correlate with known environmental parameters such as annual changes of sea level, rainfall or temperature. The strain signal could indicate local slow slip on the Cinarcik fault and thus a transfer of stress to the eastern Marmara seismic gap.
The Raman spectra of pure N-2, CO2, and CH4 were analyzed over the range 10 to 500 bars and from -160 degrees C to 200 degrees C (N-2), 22 degrees C to 350 degrees C (CO2), and -100 degrees C to 450 degrees C (CH4). At constant temperature, Raman peak position, including the more intense CO2 peak (nu+), decreases (shifts to lower wave number) with increasing pressure for all three gases over the entire pressure and temperature (PT) range studied. At constant pressure, the peak position for CO2 and CH4 increases (shifts to higher wave number) with increasing temperature over the entire PT range studied. In contrast, N-2 first shows an increase in peak position with increasing temperature at constant pressure, followed by a decrease in peak position with increasing temperature. The inflection temperature at which the trend reverses for N-2 is located between 0 degrees C and 50 degrees C at pressures above similar to 50 bars and is pressure dependent. Below similar to 50 bars, the inflection temperature was observed as low as -120 degrees C. The shifts in Raman peak positions with PT are related to relative density changes, which reflect changes in intermolecular attraction and repulsion. A conceptual model relating the Raman spectral properties of N-2, CO2, and CH4 to relative density (volume) changes and attractive and repulsive forces is presented here. Additionally, reduced temperature-dependent densimeters and barometers are presented for each pure component over the respective PT ranges. The Raman spectral behavior of the pure gases as a function of temperature and pressure is assessed to provide a framework for understanding the behavior of each component in multicomponent N-2-CO2-CH4 gas systems in a future study.
An understanding of the depositional environment and paleogeography of the Siwalik foreland basin are crucial in interpreting the basin configuration, sediment transport pathways and its evolutionary history. This study examines the sedimentology of the Siwalik succession of the Kameng River valley, Arunachal Himalaya, northeastern India. The facies characteristics of the fine-grained, well-sorted sediments of the Dafla Formation and its complex, polymodal paleocurrent pattern in this section, reveals deposition in a variety of open marine to deltaic environment. The overlying Subansiri Formation, characterized by coarse-grained, thick, multistoried sandstone, and showing more consistent SW-ward paleocurrent, indicate deposition from a large, axial braided river system. The proposed redefinition of the boundary between the Lower Siwalik Dafia and the Middle Siwalik Subansiri formations implies their transition at around 7.5 Ma, instead of 10.5 Ma, suggested earlier. The revised age of the transition is consistent with the age of arrival of the Transhimalayan sediments at 7 Ma and also denotes the time of marine to fluvial transition in this area. Presence of marine sediments in the Kameng section, with similar records further west, indicates the existence of an extensive seaway in the eastern Himalaya during the lower Siwalik time. The extant paleodrainage reconstructions have been recast on the basis of new data on the sedimentology and paleocurrent from this section. It is inferred that the changing sea level, uplifting Shillong Plateau and drainage evolution in the eastern Himalayan foreland during the middle Miocene time controlled the marine to fluvial transition in the basin.
Cosmic-ray neutron sensing (CRNS) is a promising non-invasive technique to estimate snow water equivalent (SWE) over large areas. In contrast to preliminary studies focusing on shallow snow conditions (SWE <130 mm), more recently the method was shown experimentally to be sensitive also to deeper snowpacks providing the basis for its use at mountain experimental sites. However, hysteretic neutron response has been observed for complex snow cover including patchy snow-free areas. In the present study we aimed to understand and support the experimental findings using a comprehensive neutron modeling approach. Several simulations have been set up in order to disentangle the effect on the signal of different land surface characteristics and to reproduce multiple observations during periods of snow melt and accumulation. To represent the actual land surface heterogeneity and the complex snow cover, the model used data from terrestrial laser scanning. The results show that the model was able to accurately reproduce the CRNS signal and particularly the hysteresis effect during accumulation and melting periods. Moreover, the sensor footprint was found to be anisotropic and affected by the spatial distribution of liquid water and snow as well as by the topography of the nearby mountains. Under fully snow-covered conditions the CRNS is able to accurately estimate SWE without prior knowledge about snow density profiles or other spatial anomalies. These results provide new insights into the characteristics of the detected neutron signal in complex terrain and support the use of CRNS for long-term snow monitoring in high elevated mountain environments.
Selenite pseudomorphs
(2019)
The large, shallow earthquakes at Northridge, California (1994), Chi-Chi, Taiwan (1999), and Wenchuan, China (2008), each triggered thousands of landslides. We have determined the position of these landslides along hillslopes, normalizing for statistical bias. The landslide patterns have a co-seismic signature, with clustering at ridge crests and slope toes. A cross-check against rainfall-induced landslide inventories seems to confirm that crest clustering is specific to seismic triggering as observed in previous studies. In our three study areas, the seismic ground motion parameters and lithologic and topographic features used do not seem to exert a primary control on the observed patterns of landslide clustering. However, we show that at the scale of the epicentral area, crest and toe clustering occur in areas with specific geological features. Toe clustering of seismically induced landslides tends to occur along regional major faults. Crest clustering is concentrated at sites where the lithology along hillslopes is approximately uniform, or made of alternating soft and hard strata, and without strong overprint of geological structures. Although earthquake-induced landslides locate higher on hillslopes in a statistically significant way, geological features strongly modulate the landslide position along the hillslopes. As a result the observation of landslide clustering on topographic ridges cannot be used as a definite indicator of the topographic amplification of ground shaking.
The removal, redistribution, and transient storage of sediments in tectonically active mountain belts is thought to exert a first-order control on shallow crustal stresses, fault activity, and hence on the spatiotemporal pattern of regional deformation processes. Accordingly, sediment loading and unloading cycles in intermontane sedimentary basins may inhibit or promote intrabasinal faulting, respectively, but unambiguous evidence for this potential link has been elusive so far. Here we combine 2D numerical experiments that simulate contractional deformation in a broken-foreland setting (i.e., a foreland where shortening is diachronously absorbed by spatially disparate, reverse faults uplifting basement blocks) with field data from intermontane basins in the NW Argentine Andes. Our modeling results suggest that thicker sedimentary fills (>0.7-1.0 km) may suppress basinal faulting processes, while thinner fills (<0.7 km) tend to delay faulting. Conversely, the removal of sedimentary loads via fluvial incision and basin excavation promotes renewed intrabasinal faulting. These results help to better understand the tectono-sedimentary history of intermontane basins that straddle the eastern border of the Andean Plateau in northwestern Argentina. For example, the Santa Maria and the Humahuaca basins record intrabasinal deformation during or after sediment unloading, while the Quebrada del Toro Basin reflects the suppression of intrabasinal faulting due to loading by coarse conglomerates. We conclude that sedimentary loading and unloading cycles may exert a fundamental control on spatiotemporal deformation patterns in intermontane basins of tectonically active broken forelands. (C) 2018 Elsevier B.V. All rights reserved.
In this study, a low-cost unmanned aerial vehicle was used to obtain multi-spectral high-resolution imagery (1.4 cmpx(-1)) from2 microcatchments (3.3 ha) with burned Mediterranean shrubland and pine forests. This imagery was used to calculate the blue normalized differential vegetation index and to generate digital elevation models for estimating the sediment connectivity index. Both indices enabled an integrated approach for deciphering how hydrological and sediment connectivity interact with vegetation as well as soil conservation structures. The application of spatial analysis improves our understanding of the feedback between biological and geomorphological processes. Local spatial data analysis established a significant link between local geomorphological and biological factors, enabling a precise identification of homogeneous areas at micro-catchment scale and the minimal size of vegetation units reacting to geomorphology as natural groups at plot-scale where management strategies and efforts should be applied. Establishing this local relationship between sediment connectivity and vegetation patterns through new and interdisciplinary methodologies represents a new strategy for the assessment of ecosystem dynamics and management.
The semiarid northeast of Brazil is one of the most densely populated dryland regions in the world and recurrently affected by severe droughts. Thus, reliable seasonal forecasts of streamflow and reservoir storage are of high value for water managers. Such forecasts can be generated by applying either hydrological models representing underlying processes or statistical relationships exploiting correlations among meteorological and hydrological variables. This work evaluates and compares the performances of seasonal reservoir storage forecasts derived by a process-based hydrological model and a statistical approach.
Driven by observations, both models achieve similar simulation accuracies. In a hindcast experiment, however, the accuracy of estimating regional reservoir storages was considerably lower using the process-based hydrological model, whereas the resolution and reliability of drought event predictions were similar by both approaches. Further investigations regarding the deficiencies of the process-based model revealed a significant influence of antecedent wetness conditions and a higher sensitivity of model prediction performance to rainfall forecast quality.
Within the scope of this study, the statistical model proved to be the more straightforward approach for predictions of reservoir level and drought events at regionally and monthly aggregated scales. However, for forecasts at finer scales of space and time or for the investigation of underlying processes, the costly initialisation and application of a process-based model can be worthwhile. Furthermore, the application of innovative data products, such as remote sensing data, and operational model correction methods, like data assimilation, may allow for an enhanced exploitation of the advanced capabilities of process-based hydrological models.
The semiarid northeast of Brazil is one of the most densely populated dryland regions in the world and recurrently affected by severe droughts. Thus, reliable seasonal forecasts of streamflow and reservoir storage are of high value for water managers. Such forecasts can be generated by applying either hydrological models representing underlying processes or statistical relationships exploiting correlations among meteorological and hydrological variables. This work evaluates and compares the performances of seasonal reservoir storage forecasts derived by a process-based hydrological model and a statistical approach.
Driven by observations, both models achieve similar simulation accuracies. In a hindcast experiment, however, the accuracy of estimating regional reservoir storages was considerably lower using the process-based hydrological model, whereas the resolution and reliability of drought event predictions were similar by both approaches. Further investigations regarding the deficiencies of the process-based model revealed a significant influence of antecedent wetness conditions and a higher sensitivity of model prediction performance to rainfall forecast quality.
Within the scope of this study, the statistical model proved to be the more straightforward approach for predictions of reservoir level and drought events at regionally and monthly aggregated scales. However, for forecasts at finer scales of space and time or for the investigation of underlying processes, the costly initialisation and application of a process-based model can be worthwhile. Furthermore, the application of innovative data products, such as remote sensing data, and operational model correction methods, like data assimilation, may allow for an enhanced exploitation of the advanced capabilities of process-based hydrological models.
Hydrometeorological hazards caused losses of approximately 110 billion U.S. Dollars in 2016 worldwide. Current damage estimations do not consider the uncertainties in a comprehensive way, and they are not consistent between spatial scales. Aggregated land use data are used at larger spatial scales, although detailed exposure data at the object level, such as openstreetmap.org, is becoming increasingly available across the globe.We present a probabilistic approach for object-based damage estimation which represents uncertainties and is fully scalable in space. The approach is applied and validated to company damage from the flood of 2013 in Germany. Damage estimates are more accurate compared to damage models using land use data, and the estimation works reliably at all spatial scales. Therefore, it can as well be used for pre-event analysis and risk assessments. This method takes hydrometeorological damage estimation and risk assessments to the next level, making damage estimates and their uncertainties fully scalable in space, from object to country level, and enabling the exploitation of new exposure data.
Hydrometeorological hazards caused losses of approximately 110 billion U.S. Dollars in 2016 worldwide. Current damage estimations do not consider the uncertainties in a comprehensive way, and they are not consistent between spatial scales. Aggregated land use data are used at larger spatial scales, although detailed exposure data at the object level, such as openstreetmap.org, is becoming increasingly available across the globe.We present a probabilistic approach for object-based damage estimation which represents uncertainties and is fully scalable in space. The approach is applied and validated to company damage from the flood of 2013 in Germany. Damage estimates are more accurate compared to damage models using land use data, and the estimation works reliably at all spatial scales. Therefore, it can as well be used for pre-event analysis and risk assessments. This method takes hydrometeorological damage estimation and risk assessments to the next level, making damage estimates and their uncertainties fully scalable in space, from object to country level, and enabling the exploitation of new exposure data.
Niche-based species distribution models (SDMs) have become an essential tool in conservation and restoration planning. Given the current threats to freshwater biodiversity, it is of fundamental importance to address scale effects on the performance of niche-based SDMs of freshwater species’ distributions. The scale effects are addressed here in the context of hierarchical catchment ordering, considered as counterpart to coarsening grain-size by increasing grid-cell size. We combine fish occurrence data from the Danube River Basin, the hierarchical catchment ordering and multiple environmental factors representing topographic, climatic and anthropogenic effects to model fish occurrence probability across multiple scales. We focus on 1st to 5th order catchments. The spatial scale (hierarchical catchment order) only marginally influences the mean performance of SDMs, however the uncertainty of the estimates increases with scale. Key predictors and their relative importance are scale and species dependent. Our findings have useful implications for choosing proper species dependent spatial scales for river rehabilitation measures, and for conservation planning in areas where fine grain species data are unavailable.
Interannual variations in ecosystem primary productivity are dominated by water availability. Until recently, characterizing the photosynthetic response of different ecosystems to soil moisture anomalies was hampered by observational limitations. Here, we use a number of satellite-based proxies for productivity, including spectral indices, sun-induced chlorophyll fluorescence, and data-driven estimates of gross primary production, to reevaluate the relationship between terrestrial photosynthesis and water. In contrast to nonwoody vegetation, we find a resilience of forested ecosystems to reduced soil moisture. Sun-induced chlorophyll fluorescence and data-driven gross primary production indicate an increase in photosynthesis as a result of the accompanying higher amounts of light and temperature despite lowered light-use-efficiency. Conversely, remote sensing indicators of greenness reach their detection limit and largely remain stable. Our study thus highlights the differential responses of ecosystems along a tree cover gradient and illustrates the importance of differentiating photosynthesis indicators from those of greenness for the monitoring and understanding of ecosystems. Plain Language Summary The capacity of vegetation to thrive and to sequester carbon depends on how much water they can have access to. In this work, we evaluate how different types of satellite observations can describe the response of vegetation to changes in soil moisture over the entire planet. The first source of observation measures only the greenness of the land surface, the second measures light that is emitted by pigments in plants which are photosynthetically active (chlorophyll fluorescence), and the third are simulations of gross carbon uptake derived from machine learning techniques. For periods of water shortage all three indicate a reduction of growth in ecosystems with few trees. However, in cold boreal forests, when soil moisture is particularly low, we still detect an increase in photosynthesis due to higher light and temperature conditions, but this is not reflected in the greenness indicator. This work illustrates how lack of water is not necessarily harmful for catching carbon through photosynthesis, but to monitor this effect, we need remote sensing indicators that measure more than just how green the plants are, and fluorescence is likely a good candidate.
The surface deformation associated with the 2010 M-w 8.8 Maule earthquake in Chile was recorded in great detail before, during and after the event. The high data quality of the continuous GPS (cGPS) observations has facilitated a number of studies that model the postseismic deformation signal with a combination of relocking, afterslip and viscoelastic relaxation using linear rheology for the upper mantle. Here, we investigate the impact of using linear Maxwell or power-law rheology with a 2D geomechanical-numerical model to better understand the relative importance of the different processes that control the postseismic deformation signal. Our model results reveal that, in particular, the modeled cumulative vertical postseismic deformation pattern in the near field (< 300 km from the trench) is very sensitive to the location of maximum afterslip and choice of rheology. In the model with power-law rheology, the afterslip maximum is located at 20-35 km rather than > 50 km depth as suggested in previous studies. The explanation for this difference is that in the model with power-law rheology the relaxation of coseismically imposed differential stresses occurs mainly in the lower crust. However, even though the model with power-law rheology probably has more potential to explain the vertical postseismic signal in the near field, the uncertainty of the applied temperature field is substantial, and this needs further investigations and improvements.
The Sentinel Application Platform (SNAP) architecture facilitates Earth Observation data processing. In this work, we present results from a new Snow Processor for SNAP. We also describe physical principles behind the developed snow property retrieval technique based on the analysis of Ocean and Land Colour Instrument (OLCI) onboard Sentinel-3A/B measurements over clean and polluted snow fields. Using OLCI spectral reflectance measurements in the range 400-1020 nm, we derived important snow properties such as spectral and broadband albedo, snow specific surface area, snow extent and grain size on a spatial grid of 300 m. The algorithm also incorporated cloud screening and atmospheric correction procedures over snow surfaces. We present validation results using ground measurements from Antarctica, the Greenland ice sheet and the French Alps. We find the spectral albedo retrieved with accuracy of better than 3% on average, making our retrievals sufficient for a variety of applications. Broadband albedo is retrieved with the average accuracy of about 5% over snow. Therefore, the uncertainties of satellite retrievals are close to experimental errors of ground measurements. The retrieved surface grain size shows good agreement with ground observations. Snow specific surface area observations are also consistent with our OLCI retrievals. We present snow albedo and grain size mapping over the inland ice sheet of Greenland for areas including dry snow, melted/melting snow and impurity rich bare ice. The algorithm can be applied to OLCI Sentinel-3 measurements providing an opportunity for creation of long-term snow property records essential for climate monitoring and data assimilation studies-especially in the Arctic region, where we face rapid environmental changes including reduction of snow/ice extent and, therefore, planetary albedo.
Numerous preflight investigations were necessary prior to the exposure experiment BIOMEX on the International Space Station to test the basic potential of selected microorganisms to resist or even to be active under Mars-like conditions. In this study, methanogenic archaea, which are anaerobic chemolithotrophic microorganisms whose lifestyle would allow metabolism under the conditions on early and recent Mars, were analyzed. Some strains from Siberian permafrost environments have shown a particular resistance. In this investigation, we analyzed the response of three permafrost strains (Methanosarcina soligelidi SMA-21, Candidatus Methanosarcina SMA-17, Candidatus Methanobacterium SMA-27) and two related strains from non-permafrost environments (Methanosarcina mazei, Methanosarcina barkeri) to desiccation conditions (-80 degrees C for 315 days, martian regolith analog simulants S-MRS and P-MRS, a 128-day period of simulated Mars-like atmosphere). Exposure of the different methanogenic strains to increasing concentrations of magnesium perchlorate allowed for the study of their metabolic shutdown in a Mars-relevant perchlorate environment. Survival and metabolic recovery were analyzed by quantitative PCR, gas chromatography, and a new DNA-extraction method from viable cells embedded in S-MRS and P-MRS. All strains survived the two Mars-like desiccating scenarios and recovered to different extents. The permafrost strain SMA-27 showed an increased methanogenic activity by at least 10-fold after deep-freezing conditions. The methanogenic rates of all strains did not decrease significantly after 128 days S-MRS exposure, except for SMA-27, which decreased 10-fold. The activity of strains SMA-17 and SMA-27 decreased after 16 and 60 days P-MRS exposure. Non-permafrost strains showed constant survival and methane production when exposed to both desiccating scenarios. All strains showed unaltered methane production when exposed to the perchlorate concentration reported at the Phoenix landing site (2.4 mM) or even higher concentrations. We conclude that methanogens from (non-)permafrost environments are suitable candidates for potential life in the martian subsurface and therefore are worthy of study after space exposure experiments that approach Mars-like surface conditions.
Terrestrial gravimetry is increasingly used to monitor mass transport processes in geophysics boosted by the ongoing technological development of instruments. Resolving a particular phenomenon of interest, however, requires a set of gravity corrections of which the uncertainties have not been addressed up to now. In this study, we quantify the time domain uncertainty of tide, global atmospheric, large-scale hydrological, and nontidal ocean loading corrections. The uncertainty is assessed by comparing the majority of available global models for a suite of sites worldwide. The average uncertainty expressed as root-mean-square error equals 5.1nm/s(2), discounting local hydrology or air pressure. The correction-induced uncertainty of gravity changes over various time periods of interest ranges from 0.6nm/s(2) for hours up to a maximum of 6.7nm/s(2) for 6months. The corrections are shown to be significant and should be applied for most geophysical applications of terrestrial gravimetry. From a statistical point of view, however, resolving subtle gravity effects in the order of few nanometers per square second is challenged by the uncertainty of the corrections. Plain Language Summary Many scientists are exploring ways to benefit from gravity measurements in fields of high societal relevance such as monitoring of volcanoes or measuring the amount of water in underground. Any application of such new methods, however, requires careful preparation of the gravity measurements. The intention of the preparation process is to ensure that the measurements do not contain information about processes that are not of interest. For that reason, the influence of atmosphere, ocean, tides, and hydrology needs to be reduced from the gravity. In this study, we investigate how this reduction process influences the quality of the measurement. We found that the precision degrades especially owing to the hydrology. The ocean plays an important role at sites close to the coast and the atmosphere at sites located in mountains. The overall errors of the reductions may complicate a reliable use of gravity measurements in certain studies focusing on very small signals. Nevertheless, the precision of gravity reductions alone does not obstruct a meaningful use of gravity measurements in most research fields. Details specifying the reduction precision are provided in this study allowing scientist dealing with gravity measurements to decide if their signal of interest can be reliably resolved.
The application of electrical resistivity tomography to peatlands supports conventional coring by providing data on the current condition of peatlands, including data on stratigraphy, peat properties and thickness of organic deposits. Data on the current condition of drained peatlands are particularly required to improve estimates of carbon storage as well as losses and emissions from agriculturally used peatlands. However, most of the studies focusing on electrical resistivity tomography surveys have been conducted on natural peatlands with higher groundwater levels. Peatlands drained for agriculture have not often been studied using geophysical techniques. Drained sites are characterized by low groundwater levels and high groundwater fluctuations during the year, which lead to varying levels of water saturation. To validate better electrical resistivity tomography surveys of drained peatlands, the aim of this laboratory study is to investigate the influence of varying water saturation levels on electrical conductivity (reciprocal of resistivity) for a variety of peat and gyttja types, as well as for different degrees of peat decomposition. Results show that different levels of water saturation strongly influence bulk electrical conductivity. Distinct differences in this relationship exist between peat and gyttja substrates and between different degrees of peat decomposition. Peat shows an exponential relationship for all degrees of decomposition, whereas gyttja, in particular organic-rich gyttja, is characterized by a rather unimodal relationship. The slopes for the relationship between electrical conductivity and water content are steeper at high degrees of decomposition than for peat of low degrees of decomposition. These results have direct implications for field electrical resistivity tomography surveys. In drained peatlands that are strongly susceptible to drying, electrical resistivity tomography surveys have a high potential to monitor the actual field water content. In addition, at comparable water saturations, high or low degrees of decomposition can be inferred from electrical conductivity.
Nature-based solutions (NBS) have recently received attention due to their potential ability to sustainably reduce hydro-meteorological risks, providing co-benefits for both ecosystems and affected people. Therefore, pioneering research has dedicated efforts to optimize the design of NBS, to evaluate their wider co-benefits and to understand promoting and/or hampering governance conditions for the uptake of NBS. In this article, we aim to complement this research by conducting a comprehensive literature review of factors shaping people’s perceptions of NBS as a means to reduce hydro-meteorological risks. Based on 102 studies, we identified six topics shaping the current discussion in this field of research: (1) valuation of the co-benefits (including those related to ecosystems and society); (2) evaluation of risk reduction efficacy; (3) stakeholder participation; (4) socio-economic and location-specific conditions; (5) environmental attitude, and (6) uncertainty. Our analysis reveals that concerned empirical insights are diverse and even contradictory, they vary in the depth of the insights generated and are often not comparable for a lack of a sound theoretical-methodological grounding. We, therefore, propose a conceptual model outlining avenues for future research by indicating potential inter-linkages between constructs underlying perceptions of NBS to hydro-meteorological risks.
Colombia's agriculture, forestry and other land use sector accounts for nearly half of its total greenhouse gas (GHG) emissions. The importance of smallholder deforestation is comparatively high in relation to its regional counterparts, and livestock agriculture represents the largest driver of primary forest depletion. Silvopastoral systems (SPSs) are presented as agroecological solutions that synergistically enhance livestock productivity, improve local farmers' livelihoods and hold the potential to reduce pressure on forest conversion. The department of Caquetá represents Colombia's most important deforestation hotspot. Targeting smallholder livestock farms through survey data, in this work we investigate the GHG mitigation potential of implementing SPSs for smallholder farms in this region. Specifically, we assess whether the carbon sequestration taking place in the soil and biomass of SPSs is sufficient to offset the per-hectare increase in livestock GHG emissions resulting from higher stocking rates. To address these questions we use data on livestock population characteristics and historic land cover changes reported from a survey covering 158 farms and model the carbon sequestration occurring in three different scenarios of progressively-increased SPS complexity using the CO2 fix model. We find that, even with moderate tree planting densities, the implementation of SPSs can reduce GHG emissions by 2.6 Mg CO2e ha−1 yr−1 in relation to current practices, while increasing agriculture productivity and contributing to the restoration of severely degraded landscapes.
The habilitation deals with the numerical analysis of the recurrence properties of geological and climatic processes. The recurrence of states of dynamical processes can be analysed with recurrence plots and various recurrence quantification options. In the present work, the meaning of the structures and information contained in recurrence plots are examined and described. New developments have led to extensions that can be used to describe the recurring patterns in both space and time. Other important developments include recurrence plot-based approaches to identify abrupt changes in the system's dynamics, to detect and investigate external influences on the dynamics of a system, the couplings between different systems, as well as a combination of recurrence plots with the methodology of complex networks. Typical problems in geoscientific data analysis, such as irregular sampling and uncertainties, are tackled by specific modifications and additions. The development of a significance test allows the statistical evaluation of quantitative recurrence analysis, especially for the identification of dynamical transitions. Finally, an overview of typical pitfalls that can occur when applying recurrence-based methods is given and guidelines on how to avoid such pitfalls are discussed. In addition to the methodological aspects, the application potential especially for geoscientific research questions is discussed, such as the identification and analysis of transitions in past climates, the study of the influence of external factors to ecological or climatic systems, or the analysis of landuse dynamics based on remote sensing data.
Results of a paleolimnological investigation of a well-dated lake sediment section from Shikotan Island (Southern Kurils) showed that from ca 8.0 to 5.8 cal ka BP a warm and humid period corresponding to middle Holocene optimum took place. Cooling thereafter corresponds to Neoglacial. A reconstructed from ca 0.9 to ca 0.58 cal ka BP warm period can correspond to a Medieval Warm Period. Cooling after 0.58cal ka BP can be correlated with the LIA. Marine regression stages were identified at ca 6.2-5.9, 5.5-5.1 and 1.07-0.36 cal ka BP. The general chronology of major climatic events of Holocene in the island is in accordance with the climate records from the North Pacific region. Revealed spatial differences in timing and magnitude of the Late Holocene climatic episodes (LIA, MWP) in the region needs further investigations.
Reconstructing rates and patterns of colluvial soil redistribution in agrarian (hummocky) landscapes
(2019)
Humans have triggered or accelerated erosion processes since prehistoric times through agricultural practices. Optically stimulated luminescence (OSL) is widely used to quantify phases and rates of the corresponding landscape change, by measuring the last moment of daylight exposure of sediments. However, natural and anthropogenic mixing processes, such as bioturbation and tillage, complicate the use of OSL as grains of different depositional ages become mixed, and grains become exposed to light even long after the depositional event of interest. Instead, OSL determines the stabilization age, indicating when sediments were buried below the active mixing zone. These stabilization ages can cause systematic underestimation when calculating deposition rates. Our focus is on colluvial deposition in a kettle hole in the Uckermark region, northeastern Germany. We took 32 samples from five locations in the colluvium filling the kettle hole to study both spatial and temporal patterns in colluviation. We combined OSL dating with advanced age modelling to determine the stabilization age of colluvial sediments. These ages were combined with an archaeological reconstruction of historical ploughing depths to derive the levels of the soil surface at the moment of stabilization; the deposition depths, which were then used to calculate unbiased deposition rates. We identified two phases of colluvial deposition. The oldest deposits (similar to 5 ka) were located at the fringe of the kettle hole and accumulated relatively slowly, whereas the youngest deposits (<0.3 ka) rapidly filled the central kettle hole with rates of two orders of magnitude higher. We suggest that the latter phase is related to artificial drainage, facilitating accessibility in the central depression for agricultural practices. Our results show the need for numerical dating techniques that take archaeological and soil-geomorphological information into account to identify spatiotemporal patterns of landscape change, and to correctly interpret landscape dynamics in anthropogenically influenced hilly landscapes. (c) 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.
RAA2019
(2019)
These abstracts result from the 10th International Congress on the Application of Raman Spectroscopy in Art and Archaeology held 03.09. – 07.09.2019 in Potsdam (Germany).
The RAA is an established biennial international conference series. Since the beginning in 2001, the RAA conferences promote Raman Spectroscopy and play an important role in increasing the field of its applications in art history, history, archaeology, palaeontology, conservation and restoration, museology, degradation of cultural heritage, archaeometry, etc. Furthermore, the development of new instrumentation, especially for non-invasive measurements, receives great attention.
The Congress covers all topics of Raman spectroscopic applications in art and archaeology and focuses on the following themes:
• Material characterization and degradation processes
• Conservation issues affecting cultural heritage
• Raman spectroscopy of biological and organic materials
• Surface enhanced Raman spectroscopy
• Chemometrics in Raman spectroscopy
• Development of Raman techniques
• New Raman instrumentation and applications in cultural heritage objects investigations
• Raman spectroscopy in paleontology, paleoenvironment and archaeology
Oxidation of particulate organic carbon (POC) during fluvial transit releases CO2 to the atmosphere and can influence global climate. Field data show large POC oxidation fluxes in lowland rivers; however, it is unclear if POC losses occur predominantly during in-river transport, where POC is in continual motion within an aerated environment, or during transient storage in floodplains, which may be anoxic. Determination of the locus of POC oxidation in lowland rivers is needed to develop process-based models to predict POC losses, constrain carbon budgets, and unravel links between climate and erosion. However, sediment exchange between rivers and floodplains makes differentiating POC oxidation during in-river transport from oxidation during floodplain storage difficult. Here, we isolated inriver POC oxidation using flume experiments transporting petrogenic and biospheric POC without floodplain storage. Our experiments showed solid phase POC losses of 0%-10% over similar to 10(3) km of fluvial transport, compared to similar to 7% to >50% losses observed in rivers over similar distances. The production of dissolved organic carbon (DOC) and dissolved rhenium (a proxy for petrogenic POC oxidation) was consistent with small POC lasses, and replicate experiments in static water tanks gave similar results. Our results show that fluvial sediment transport, particle abrasion, and turbulent mixing have a minimal role on POC oxidation, and they suggest that POC losses may accrue primarily in floodplain storage.
We study optimal and equilibrium sizes of cities in a city system model with pollution. Pollution is a function of population size. If pollution is local or per-capita pollution increases with population, equilibrium cities are too large under symmetry; with asymmetric cities, the largest cities are too large and the smallest too small. When pollution is global and per-capita pollution declines with city size, cities may be too small under symmetry; with asymmetric cities, the largest cities are too small and the smallest too large if the marginal damage of pollution is large enough. We calibrate the model to US cities and find that the largest cities may be undersized by 3-4%.
Pollen-based quantitative land-cover reconstruction for northern Asia covering the last 40 ka cal BP
(2019)
We collected the available relative pollen productivity estimates (PPEs) for 27 major pollen taxa from Eurasia and applied them to estimate plant abundances during the last 40 ka cal BP (calibrated thousand years before present) using pollen counts from 203 fossil pollen records in northern Asia (north of 40 degrees N). These pollen records were organized into 42 site groups and regional mean plant abundances calculated using the REVEALS (Regional Estimates of Vegetation Abundance from Large Sites) model. Time-series clustering, constrained hierarchical clustering, and detrended canonical correspondence analysis were performed to investigate the regional pattern, time, and strength of vegetation changes, respectively. Reconstructed regional plant functional type (PFT) components for each site group are generally consistent with modern vegetation in that vegetation changes within the regions are characterized by minor changes in the abundance of PFTs rather than by an increase in new PFTs, particularly during the Holocene. We argue that pollen-based REVEALS estimates of plant abundances should be a more reliable reflection of the vegetation as pollen may overestimate the turnover, particularly when a high pollen producer invades areas dominated by low pollen producers. Comparisons with vegetation-independent climate records show that climate change is the primary factor driving land-cover changes at broad spatial and temporal scales. Vegetation changes in certain regions or periods, however, could not be explained by direct climate change, e.g. inland Siberia, where a sharp increase in evergreen conifer tree abundance occurred at ca. 7-8 ka cal BP despite an unchanging climate, potentially reflecting their response to complex climate-permafrost-fire-vegetation interactions and thus a possible long-term lagged climate response.
The intermontane Humahuaca Basin in the Eastern Cordillera of the northwest Argentine Andes lies leeward of an orographic barrier to easterly derived moisture. An average of >2000 mm/yr of rainfall along the eastern flanks of the barrier contrasts with <200 mm/yr in the orogen interior. Paleoenvironmental reconstructions suggest that the basin became disconnected from the foreland during the Miocene-Pliocene by the growth of fault-bounded mountain ranges. Fossil records, sedimentology, and stable isotope data imply that rerouting of the fluvial network by 4.2 Ma and reduced rainfall by ca. 3 Ma were consequences of that range uplift. Here, we present cosmogenic nuclide-derived (Be-10) paleodenudation rates from 6 to 2 Ma fluvial deposits collected from the Humahuaca Basin. Despite increased tectonic activity, our Be-10 data show a tenfold decrease in denudation rates at ca. 3 Ma, documenting a link between uplift-induced semiarid conditions and decreasing hillslope denudation rates. This new data set thus demonstrates the influence of hydrological change on spatiotemporal denudation patterns in tectonically active mountain areas.
Sinkholes and depressions are typical landforms of karst regions. They pose a considerable natural hazard to infrastructure, agriculture, economy and human life in affected areas worldwide. The physio-chemical processes of sinkholes and depression formation are manifold, ranging from dissolution and material erosion in the subsurface to mechanical subsidence/failure of the overburden. This thesis addresses the mechanisms leading to the development of sinkholes and depressions by using complementary methods: remote sensing, distinct element modelling and near-surface geophysics.
In the first part, detailed information about the (hydro)-geological background, ground structures, morphologies and spatio-temporal development of sinkholes and depressions at a very active karst area at the Dead Sea are derived from satellite image analysis, photogrammetry and geologic field surveys. There, clusters of an increasing number of sinkholes have been developing since the 1980s within large-scale depressions and are distributed over different kinds of surface materials: clayey mud, sandy-gravel alluvium and lacustrine evaporites (salt). The morphology of sinkholes differs depending in which material they form: Sinkholes in sandy-gravel alluvium and salt are generally deeper and narrower than sinkholes in the interbedded evaporite and mud deposits. From repeated aerial surveys, collapse precursory features like small-scale subsidence, individual holes and cracks are identified in all materials. The analysis sheds light on the ongoing hazardous subsidence process, which is driven by the base-level fall of the Dead Sea and by the dynamic formation of subsurface water channels.
In the second part of this thesis, a novel, 2D distinct element geomechanical modelling approach with the software PFC2D-V5 to simulating individual and multiple cavity growth and sinkhole and large-scale depression development is presented. The approach involves a stepwise material removal technique in void spaces of arbitrarily shaped geometries and is benchmarked by analytical and boundary element method solutions for circular cavities. Simulated compression and tension tests are used to calibrate model parameters with bulk rock properties for the materials of the field site. The simulations show that cavity and sinkhole evolution is controlled by material strength of both overburden and cavity host material, the depth and relative speed of the cavity growth and the developed stress pattern in the subsurface. Major findings are: (1) A progressively deepening differential subrosion with variable growth speed yields a more fragmented stress pattern with stress interaction between the cavities. It favours multiple sinkhole collapses and nesting within large-scale depressions. (2) Low-strength materials do not support large cavities in the material removal zone, and subsidence is mainly characterised by gradual sagging into the material removal zone with synclinal bending. (3) High-strength materials support large cavity formation, leading to sinkhole formation by sudden collapse of the overburden. (4) Large-scale depression formation happens either by coalescence of collapsing holes, block-wise brittle failure, or gradual sagging and lateral widening.
The distinct element based approach is compared to results from remote sensing and geophysics at the field site. The numerical simulation outcomes are generally in good agreement with derived morphometrics, documented surface and subsurface structures as well as seismic velocities. Complementary findings on the subrosion process are provided from electric and seismic measurements in the area.
Based on the novel combination of methods in this thesis, a generic model of karst landform evolution with focus on sinkhole and depression formation is developed. A deepening subrosion system related to preferential flow paths evolves and creates void spaces and subsurface conduits. This subsequently leads to hazardous subsidence, and the formation of sinkholes within large-scale depressions. Finally, a monitoring system for shallow natural hazard phenomena consisting of geodetic and geophysical observations is proposed for similarly affected areas.
A well-preserved sequence, by Archean standards, of mantle-derived metabasalts and metakomatiites forms large parts of the lower Onverwacht Group of the Barberton Greenstone Belt (South Africa). To elucidate the origin of mafic and ultramafic rocks from this 3.55 to 3.45 Ga sequence, we present a comprehensive geochemical dataset including major and trace elements as well as Lu-Hf and Sm-Nd isotope compositions for a variety of metavolcanic rocks. These include metabasalts of the amphibolite-facies Sandspruit and Theespruit Formations as well as metabasalts and metakomatiites of the lower greenschist-facies Komati Formation. Based on their incompatible trace element patterns, the basalts of the Sandspruit and Theespruit Formations can be subdivided into a light rare earth element (LREE) depleted group, a LREE-undepleted group, and a LREE-enriched group. Positive epsilon Hf-(t) and epsilon Nd-(t) values of ca. +3 to +4 and 0 to +2, respectively, together with depletions in Th and La-CN/Yb-CN indicate derivation of the LREE-depleted basalts from a depleted mantle source. However, chondritic epsilon Hf-(t) and epsilon Nd-(t) values combined with positive Th and La-CN/Yb-CN of the LREE-enriched samples indicate a contribution from older granitoid crust in the petrogenesis of these samples. Trace element patterns of komatiites and basalts of the Komati Formation are generally flat relative to primitive mantle with slight depletions in heavy rare earth elements and Th and overall positive epsilon Hf-(t) of + 2.5 +/- 3.5 (2 s.d.) and epsilon Nd-(t) of + 0.5 +/- 2.2 (2 s. d.). The coherence in trace element characteristics suggests a common magmatic origin for basalts and komatiites. This study reveals that the two lavas were derived from the same mantle plume, i. e. komatiites were formed by high degrees of melting of a depleted mantle source containing residual garnet and the basalts were formed by moderate degrees of partial melting in shallower regions of the mantle. Based on the current dataset, combined with published data, we propose a geodynamic model for the oldest units of the Barberton Greenstone Belt that describes the development from a submerged continental setting (for the Sandspruit and Theespruit Formations) to a submarine plateau setting (for the Komati Formation) as a consequence of continental rifting.
Main group pallasite meteorites are samples of a single early magmatic planetesimal, dominated by metal and olivine but containing accessory chromite, sulfide, phosphide, phosphates, and rare phosphoran olivine. They represent mixtures of core and mantle materials, but the environment of formation is poorly understood, with a quiescent core-mantle boundary, violent core-mantle mixture, or surface mixture all recently suggested. Here, we review main group pallasite data sets and petrologic characteristics, and present new observations on the low-MnO pallasite Brahin that contains abundant fragmental olivine, but also rounded and angular olivine and potential evidence of sulfide-phosphide liquid immiscibility. A reassessment of the literature shows that low-MnO and high-FeO subgroups preferentially host rounded olivine and low-temperature P2O5-rich phases such as the Mg-phosphate farringtonite and phosphoran olivine. These phases form after metal and silicate reservoirs back-react during decreasing temperature after initial separation, resulting in oxidation of phosphorus and chromium. Farringtonite and phosphoran olivine have not been found in the common subgroup PMG, which are mechanical mixtures of olivine, chromite with moderate Al2O3 contents, primitive solid metal, and evolved liquid metal. Lower concentrations of Mn in olivine of the low-MnO PMG subgroup, and high concentrations of Mn in low-Al2O3 chromites, trace the development and escape of sulfide-rich melt in pallasites and the partially chalcophile behavior for Mn in this environment. Pallasites with rounded olivine indicate that the core-mantle boundary of their planetesimal may not be a simple interface but rather a volume in which interactions between metal, silicate, and other components occur.
We measure valence-to-core x-ray emission spectra of compressed crystalline GeO₂ up to 56 GPa and of amorphous GeO₂ up to 100 GPa. In a novel approach, we extract the Ge coordination number and mean Ge-O distances from the emission energy and the intensity of the Kβ'' emission line. The spectra of high-pressure polymorphs are calculated using the Bethe-Salpeter equation. Trends observed in the experimental and calculated spectra are found to match only when utilizing an octahedral model. The results reveal persistent octahedral Ge coordination with increasing distortion, similar to the compaction mechanism in the sequence of octahedrally coordinated crystalline GeO₂ high-pressure polymorphs.
We measure valence-to-core x-ray emission spectra of compressed crystalline GeO₂ up to 56 GPa and of amorphous GeO₂ up to 100 GPa. In a novel approach, we extract the Ge coordination number and mean Ge-O distances from the emission energy and the intensity of the Kβ'' emission line. The spectra of high-pressure polymorphs are calculated using the Bethe-Salpeter equation. Trends observed in the experimental and calculated spectra are found to match only when utilizing an octahedral model. The results reveal persistent octahedral Ge coordination with increasing distortion, similar to the compaction mechanism in the sequence of octahedrally coordinated crystalline GeO₂ high-pressure polymorphs.
James Ross Island (JRI) offers the exceptional opportunity to study microbial-driven pedogenesis without the influence of vascular plants or faunal activities (e.g., penguin rookeries). In this study, two soil profiles from JRI (one at Santa Martha Cove - SMC, and another at Brandy Bay BB) were investigated, in order to gain information about the initial state of soil formation and its interplay with prokaryotic activity, by combining pedological, geochemical and microbiological methods. The soil profiles are similar with respect to topographic position and parent material but are spatially separated by an orographic barrier and therefore represent windward and leeward locations towards the mainly southwesterly winds. These different positions result in differences in electric conductivity of the soils caused by additional input of bases by sea spray at the windward site and opposing trends in the depth functions of soil pH and electric conductivity. Both soils are classified as Cryosols, dominated by bacterial taxa such as Actinobacteria, Proteobacteria, Acidobacteria, Gemmatimonadetes and Chloroflexi. A shift in the dominant taxa was observed below 20 cm in both soils as well as an increased abundance of multiple operational taxonomic units (OTUs) related to potential chemolithoautotrophic Acidiferrobacteraceae. This shift is coupled by a change in microstructure. While single/pellicular grain microstructure (SMC) and platy microstructure (BB) are dominant above 20 cm, lenticular microstructure is dominant below 20 cm in both soils. The change in microstructure is caused by frequent freeze-thaw cycles and a relative high water content, and it goes along with a development of the pore spacing and is accompanied by a change in nutrient content. Multivariate statistics revealed the influence of soil parameters such as chloride, sulfate, calcium and organic carbon contents, grain size distribution and pedogenic oxide ratios on the overall microbial community structure and explained 49.9% of its variation. The correlation of the pedogenic oxide ratios with the compositional distribution of microorganisms as well as the relative abundance certain microorganisms such as potentially chemolithotrophic Acidiferrobacteraceae-related OTUs could hint at an interplay between soil-forming processes and microorganisms.
Past climate and continentality inferred from ice wedges at Batagay Highlands, interior Yakutia
(2019)
Ice wedges in the Yana Highlands of interior Yakutia - the most continental region of the Northern Hemisphere - were investigated to elucidate changes in winter climate and continentality that have taken place since the Middle Pleistocene. The Batagay megaslump exposes ice wedges and composite wedges that were sampled from three cryostratigraphic units: the lower ice complex of likely pre-Marine Isotope Stage (MIS) 6 age, the upper ice complex (Yedoma) and the upper sand unit (both MIS 3 to 2). A terrace of the nearby Adycha River provides a Late Holocene (MIS 1) ice wedge that serves as a modern reference for interpretation. The stable-isotope composition of ice wedges in the MIS 3 upper ice complex at Batagay is more depleted (mean delta O-18 about -35 parts per thousand) than those from 17 other ice-wedge study sites across coastal and central Yakutia. This observation points to lower winter temperatures and therefore higher continentality in the Yana Highlands during MIS 3. Likewise, more depleted isotope values are found in Holocene wedge ice (mean delta O-18 about -29 parts per thousand) compared to other sites in Yakutia. Ice-wedge isotopic signatures of the lower ice complex mean delta O-18 about -33 parts per thousand) and of the MIS 3-2 upper sand unit (mean delta O-18 from about -33 parts per thousand to -30 parts per thousand) are less distinctive regionally. The latter unit preserves traces of fast formation in rapidly accumulating sand sheets and of post-depositional isotopic fractionation.
Rapid local adaptation frequently occurs during the spread of invading species. It remains unclear, however, how consistent, and therefore potentially predictable, such patterns of local adaptation are. One approach to this question is to measure patterns of local differentiation in functional traits and plasticity levels in invasive species in multiple regions. Finding consistent patterns of local differentiation in replicate regions suggests that these patterns are adaptive. Further, this outcome indicates that the invading species likely responds predictably to selection along environmental gradients, even though standing genetic variation is likely to have been reduced during introduction. We studied local differentiation in the invasive annual plant Erodium cicutarium in two invaded regions, California and Chile. We collected seeds from across strong gradients in precipitation and temperature in Mediterranean-climate parts of the two regions (10 populations per region). We grew seeds from maternal families from these populations through two generations and exposed the second generation to contrasting levels of water and nutrient availability. We measured growth, flowering time and leaf functional traits across these treatments to obtain trait means and plasticity measures. We found strong differentiation among populations in all traits. Plants from drier environments flowered earlier, were less plastic in flowering time and reached greater size in all treatments. Correlations among traits within regions suggested a coordinated evolutionary response along environmental gradients associated with growing season length. There was little divergence in traits and trait intercorrelations between regions, but strongly parallel divergence in traits within regions. Similar, statistically consistent patterns of local trait differentiation across two regions suggest that local adaptation to environmental gradients has aided the spread of this invasive species, and that the formation of ecotypes in newly invaded environments has been relatively consistent and predictable.
Interactions and feedbacks between tectonics, climate, and upper plate architecture control basin geometry, relief, and depositional systems. The Andes is part of a longlived continental margin characterized by multiple tectonic cycles which have strongly modified the Andean upper plate architecture. In the Andean retroarc, spatiotemporal variations in the structure of the upper plate and tectonic regimes have resulted in marked along-strike variations in basin geometry, stratigraphy, deformational style, and mountain belt morphology. These along-strike variations include high-elevation plateaus (Altiplano and Puna) associated with a thin-skin fold-and-thrust-belt and thick-skin deformation in broken foreland basins such as the Santa Barbara system and the Sierras Pampeanas. At the confluence of the Puna Plateau, the Santa Barbara system and the Sierras Pampeanas, major along-strike changes in upper plate architecture, mountain belt morphology, basement exhumation, and deformation style can be recognized. I have used a source to sink approach to unravel the spatiotemporal tectonic evolution of the Andean retroarc between 26 and 28°S. I obtained a large low-temperature thermochronology data set from basement units which includes apatite fission track, apatite U-Th-Sm/He, and zircon U-Th/He (ZHe) cooling ages. Stratigraphic descriptions of Miocene units were temporally constrained by U-Pb LA-ICP-MS zircon ages from interbedded pyroclastic material.
Modeled ZHe ages suggest that the basement of the study area was exhumed during the Famatinian orogeny (550-450 Ma), followed by a period of relative tectonic quiescence during the Paleozoic and the Triassic. The basement experienced horst exhumation during the Cretaceous development of the Salta rift. After initial exhumation, deposition of thick Cretaceous syn-rift strata caused reheating of several basement blocks within the Santa Barbara system. During the Eocene-Oligocene, the Andean compressional setting was responsible for the exhumation of several disconnected basement blocks. These exhumed blocks were separated by areas of low relief, in which humid climate and low erosion rates facilitated the development of etchplains on the crystalline basement. The exhumed basement blocks formed an Eocene to Oligocene broken foreland basin in the back-bulge depozone of the Andean foreland. During the Early Miocene, foreland basin strata filled up the preexisting Paleogene topography. The basement blocks in lower relief positions were reheated; associated geothermal gradients were higher than 25°C/km. Miocene volcanism was responsible for lateral variations on the amount of reheating along the Campo-Arenal basin. Around 12 Ma, a new deformational phase modified the drainage network and fragmented the lacustrine system. As deformation and rock uplift continued, the easily eroded sedimentary cover was efficiently removed and reworked by an ephemeral fluvial system, preventing the development of significant relief. After ~6 Ma, the low erodibility of the basement blocks which began to be exposed caused relief increase, leading to the development of stable fluvial systems. Progressive relief development modified atmospheric circulation, creating a rainfall gradient. After 3 Ma, orographic rainfall and high relief lead to the development of proximal fluvial-gravitational depositional systems in the surrounding basins.
Paleogene evolution and demise of the proto-Paratethys Sea in Central Asia (Tarim and Tajik basins)
(2019)
The proto-Paratethys Sea covered a vast area extending from the Mediterranean Tethys to the Tarim Basin in western China during Cretaceous and early Paleogene. Climate modelling and proxy studies suggest that Asian aridification has been governed by westerly moisture modulated by fluctuations of the proto-Paratethys Sea. Transgressive and regressive episodes of the proto-Paratethys Sea have been previously recognized but their timing, extent and depositional environments remain poorly constrained. This hampers understanding of their driving mechanisms (tectonic and/or eustatic) and their contribution to Asian aridification. Here, we present a new chronostratigraphic framework based on biostratigraphy and magnetostratigraphy as well as a detailed palaeoenvironmental analysis for the Paleogene proto-Paratethys Sea incursions in the Tajik and Tarim basins. This enables us to identify the major drivers of marine fluctuations and their potential consequences on Asian aridification. A major regional restriction event, marked by the exceptionally thick (<= 400 m) shelf evaporites is assigned a Danian-Selandian age (ca. 63-59 Ma) in the Aertashi Formation. This is followed by the largest recorded proto-Paratethys Sea incursion with a transgression estimated as early Thanetian (ca. 59-57 Ma) and a regression within the Ypresian (ca. 53-52 Ma), both within the Qimugen Formation. The transgression of the next incursion in the Kalatar and Wulagen formations is now constrained as early Lutetian (ca. 47-46 Ma), whereas its regression in the Bashibulake Formation is constrained as late Lutetian (ca. 41 Ma) and is associated with a drastic increase in both tectonic subsidence and basin infilling. The age of the final and least pronounced sea incursion restricted to the westernmost margin of the Tarim Basin is assigned as Bartonian-Priabonian (ca. 39.7-36.7 Ma). We interpret the long-term westward retreat of the proto-Paratethys Sea starting at ca. 41 Ma to be associated with far-field tectonic effects of the Indo-Asia collision and Pamir/Tibetan plateau uplift. Short-term eustatic sea level transgressions are superimposed on this long-term regression and seem coeval with the transgression events in the other northern Peri-Tethyan sedimentary provinces for the 1st and 2nd sea incursions. However, the 3rd sea incursion is interpreted as related to tectonism. The transgressive and regressive intervals of the proto-Paratethys Sea correlate well with the reported humid and arid phases, respectively in the Qaidam and Xining basins, thus demonstrating the role of the proto-Paratethys Sea as an important moisture source for the Asian interior and its regression as a contributor to Asian aridification.
A number of sedimentary provenance studies have been undertaken in order to determine whether the palaeo-Red River was once a river of continental proportions into which the upper reaches of the Yangtze, Salween, Mekong, Irrawaddy, and Yarlung drained. We have assessed the evidence that the Yarlung originally flowed into the palaeo-Red river, and then sequentially into the Irrawaddy and Brahmaputra, connecting to the latter first via the Lohit and then the Siang. For this river system, we have integrated our new data from the Paleogene-Recent Irrawaddy drainage basin (detrital zircon U-Pb with Hf and fission track, rutile U-Pb, mica Ar-Ar, bulk rock Sr-Nd, and petrography) with previously published data, to produce a palaeodrainage model that is consistent with all datasets. In our model, the Yarlung never flowed into the Irrawaddy drainage: during the Paleogene, the Yarlung suture zone was an internally drained basin, and from Neogene times onwards the Yarlung drained into the Brahmaputra in the Bengal Basin. The Central Myanmar Basin, through which the Irrawaddy River flows today, received predominantly locally-derived detritus until the Middle Eocene, the Irrawaddy initiated as a through-going river draining the Mogok Metamorphic Belt and Bomi-Chayu granites to the north sometime in the Late Eocene to Early Oligocene, and the river was dominated by a stable MMB-dominated drainage throughout the Neogene to present day. Existing evidence does not support any connection between the Yarlung and the Red River in the past, but there is a paucity of suitable palaeo-Red River deposits with which to make a robust comparison. We argue that this limitation also precludes a robust assessment of a palaeo-connection between the Yangtze/ Salween/Mekong and the Red River; it is difficult to unequivocally interpret the recorded provenance changes as the result of specific drainage reorganisations. We highlight the palaeo-Red River deposits of the Hanoi Basin as a potential location for future research focus in view of the near-complete Cenozoic record of palaeo-Red River deposits at this location. A majority of previous studies consider that if a major continental-scale drainage ever existed at all, it fragmented early in the Cenozoic. Such a viewpoint would agree with the growing body of evidence from palaeoaltitude studies that large parts of SE Tibet were uplifted by this period. This then leads towards the intriguing question as to the mechanisms which caused the major period of river incision in the Miocene in this region.
The progress of science is tied to the standardization of measurements, instruments, and data. This is especially true in the Big Data age, where analyzing large data volumes critically hinges on the data being standardized. Accordingly, the lack of community-sanctioned data standards in paleoclimatology has largely precluded the benefits of Big Data advances in the field. Building upon recent efforts to standardize the format and terminology of paleoclimate data, this article describes the Paleoclimate Community reporTing Standard (PaCTS), a crowdsourced reporting standard for such data. PaCTS captures which information should be included when reporting paleoclimate data, with the goal of maximizing the reuse value of paleoclimate data sets, particularly for synthesis work and comparison to climate model simulations. Initiated by the LinkedEarth project, the process to elicit a reporting standard involved an international workshop in 2016, various forms of digital community engagement over the next few years, and grassroots working groups. Participants in this process identified important properties across paleoclimate archives, in addition to the reporting of uncertainties and chronologies; they also identified archive-specific properties and distinguished reporting standards for new versus legacy data sets. This work shows that at least 135 respondents overwhelmingly support a drastic increase in the amount of metadata accompanying paleoclimate data sets. Since such goals are at odds with present practices, we discuss a transparent path toward implementing or revising these recommendations in the near future, using both bottom-up and top-down approaches.