Refine
Year of publication
- 2022 (164) (remove)
Document Type
- Article (132)
- Doctoral Thesis (19)
- Postprint (9)
- Monograph/Edited Volume (2)
- Part of a Book (1)
- Review (1)
Is part of the Bibliography
- yes (164)
Keywords
- climate change (7)
- permafrost (7)
- Andes (3)
- diffusion (3)
- machine learning (3)
- ocean color remote sensing (3)
- radiation belts (3)
- Arctic ocean (2)
- Argentina (2)
- Central Andes (2)
Institute
- Institut für Geowissenschaften (164) (remove)
Ground motion with strong-velocity pulses can cause significant damage to buildings and structures at certain periods; hence, knowing the period and velocity amplitude of such pulses is critical for earthquake structural engineering.
However, the physical factors relating the scaling of pulse periods with magnitude are poorly understood.
In this study, we investigate moderate but damaging earthquakes (M-w 6-7) and characterize ground- motion pulses using the method of Shahi and Baker (2014) while considering the potential static-offset effects.
We confirm that the within-event variability of the pulses is large. The identified pulses in this study are mostly from strike-slip-like earthquakes. We further perform simulations using the freq uency-wavenumber algorithm to investigate the causes of the variability of the pulse periods within and between events for moderate strike-slip earthquakes.
We test the effect of fault dips, and the impact of the asperity locations and sizes. The simulations reveal that the asperity properties have a high impact on the pulse periods and amplitudes at nearby stations.
Our results emphasize the importance of asperity characteristics, in addition to earthquake magnitudes for the occurrence and properties of pulses produced by the forward directivity effect.
We finally quantify and discuss within- and between-event variabilities of pulse properties at short distances.
Among the multitude of geomorphological processes, aeolian shaping processes are of special character, Pedogenic dust is one of the most important sources of atmospheric aerosols and therefore regarded as a key player for atmospheric processes. Soil dust emissions, being complex in composition and properties, influence atmospheric processes and air quality and has impacts on other ecosystems. In this because even though their immediate impact can be considered low (exceptions exist), their constant and large-scale force makes them a powerful player in the earth system. dissertation, we unravel a novel scientific understanding of this complex system based on a holistic dataset acquired during a series of field experiments on arable land in La Pampa, Argentina. The field experiments as well as the generated data provide information about topography, various soil parameters, the atmospheric dynamics in the very lower atmosphere (4m height) as well as measurements regarding aeolian particle movement across a wide range of particle size classes between 0.2μm up to the coarse sand.
The investigations focus on three topics: (a) the effects of low-scale landscape structures on aeolian transport processes of the coarse particle fraction, (b) the horizontal and vertical fluxes of the very fine particles and (c) the impact of wind gusts on particle emissions.
Among other considerations presented in this thesis, it could in particular be shown, that even though the small-scale topology does have a clear impact on erosion and deposition patterns, also physical soil parameters need to be taken into account for a robust statistical modelling of the latter. Furthermore, specifically the vertical fluxes of particulate matter have different characteristics for the particle size classes. Finally, a novel statistical measure was introduced to quantify the impact of wind gusts on the particle uptake and its application on the provided data set. The aforementioned measure shows significantly increased particle concentrations during points in time defined as gust event.
With its holistic approach, this thesis further contributes to the fundamental understanding of how atmosphere and pedosphere are intertwined and affect each other.
Wetlands are dynamic ecosystems that require continuous monitoring and assessment of degradation status to design strategies for their sustainable management. While hydrology provides the primary functional control for the wetland ecosystem, the loss of landscape connectivity influences wetland degradation in a major way as it leads to fragmentation. This article aims to integrate hydrogeomorphic and ecological concepts for the assessment of degradation status and its causal factors for a large wetland in the western Ganga plains, India, the Haiderpur, using a wetlandscape approach. We have used a remote-sensing-based approach, which offers a powerful tool for assessing and linking cross-scale structures, functions, and controls in a wetlandscape. The Haiderpur, a Ramsar site since December 2021, is an artificial wetland located on the right bank of the Ganga River wherein the inflows are controlled by a barrage constructed on the Ganga River apart from smaller tributaries flowing in from the north. A novel aspect of this work is the integration of river dynamics and its connectivity to the wetlandscape to understand the spatiotemporal variability in the waterspread area in the wetland. In this work, we have developed an integrated wetlandscape assessment approach by evaluating wetland's geomorphic and hydrological connectivity status for the period 1993-2019 (25 years) across three different spatial scales - regional, catchment, and wetland. We have highlighted the ecological implications of connectivity and patch dynamics for developing sustainable wetland management plans.
Humankind and their environment need to be protected from the harmful effects of spent nuclear fuel, and therefore disposal in deep geological formations is favoured worldwide. Suitability of potential host rocks is evaluated, among others, by the retention capacity with respect to radionuclides. Safety assessments are based on the quantification of radionuclide migration lengths with numerical simulations as experiments cannot cover the required temporal (1 Ma) and spatial scales (>100 m).
Aim of the present thesis is to assess the migration of uranium, a geochemically complex radionuclide, in the potential host rock Opalinus Clay. Radionuclide migration in clay formations is governed by diffusion due to their low permeability and retarded by sorption. Both processes highly depend on pore water geochemistry and mineralogy that vary between different facies. Diffusion is quantified with the single-component (SC) approach using one diffusion coefficient for all species and the process-based multi-component (MC) option. With this, each species is assigned its own diffusion coefficient and the interaction with the diffuse double layer is taken into account. Sorption is integrated via a bottom-up approach using mechanistic surface complexation models and cation exchange. Therefore, reactive transport simulations are conducted with the geochemical code PHREEQC to quantify uranium migration, i.e. diffusion and sorption, as a function of mineralogical and geochemical heterogeneities on the host rock scale.
Sorption processes are facies dependent. Migration lengths vary between the Opalinus Clay facies by up to 10 m. Thereby, the geochemistry of the pore water, in particular the partial pressure of carbon dioxide (pCO2), is more decisive for the sorption capacity than the amount of clay minerals. Nevertheless, higher clay mineral quantities compensate geochemical variations. Consequently, sorption processes must be quantified as a function of pore water geochemistry in contact with the mineral assemblage.
Uranium diffusion in the Opalinus Clay is facies independent. Speciation is dominated by aqueous ternary complexes of U(VI) with calcium and carbonate. Differences in the migration lengths between SC and MC diffusion are with +/-5 m negligible. Further, the application of the MC approach highly depends on the quality and availability of the underlying data. Therefore, diffusion processes can be adequately quantified with the SC approach using experimentally determined diffusion coefficients.
The hydrogeological system governs pore water geochemistry within the formation rather than the mineralogy. Diffusive exchange with the adjacent aquifers established geochemical gradients over geological time scales that can enhance migration by up to 25 m. Consequently, uranium sorption processes must be quantified following the identified priority: pCO2 > hydrogeology > mineralogy.
The presented research provides a workflow and orientation for other potential disposal sites with similar pore water geochemistry due to the identified mechanisms and dependencies. With a maximum migration length of 70 m, the retention capacity of the Opalinus Clay with respect to uranium is sufficient to fulfill the German legal minimum requirement of a thickness of at least 100 m.
Groundwater is critical in supporting current and future reliable water supply throughout Africa. Although continental maps of groundwater storage and recharge have been developed, we currently lack a clear understanding on how the controls on groundwater recharge vary across the entire continent. Reviewing the existing literature, we synthesize information on reported groundwater recharge controls in Africa. We find that 15 out of 22 of these controls can be characterised using global datasets. We develop 11 descriptors of climatic, topographic, vegetation, soil and geologic properties using global datasets, to characterise groundwater recharge controls in Africa. These descriptors cluster Africa into 15 Recharge Landscape Units for which we expect recharge controls to be similar. Over 80% of the continents land area is organized by just nine of these units. We also find that aggregating the Units by similarity into four broader Recharge Landscapes (Desert, Dryland, Wet tropical and Wet tropical forest) provides a suitable level of landscape organisation to explain differences in ground-based long-term mean annual recharge and recharge ratio (annual recharge / annual precipitation) estimates. Furthermore, wetter Recharge Landscapes are more efficient in converting rainfall to recharge than drier Recharge Landscapes as well as having higher annual recharge rates. In Dryland Recharge Landscapes, we found that annual recharge rates largely varied according to mean annual precipitation, whereas recharge ratio estimates increase with increasing monthly variability in P-PET. However, we were unable to explain why ground based estimates of recharge signatures vary across other Recharge Landscapes, in which there are fewer ground based recharge estimates, using global datasets alone. Even in dryland regions, there is still considerable unexplained variability in the estimates of annual recharge and recharge ratio, stressing the limitations of global datasets for investigating ground-based information.
Statistical distributions of flood peak discharge often show heavy tail behavior, that is, extreme floods are more likely to occur than would be predicted by commonly used distributions that have exponential asymptotic behavior.
This heavy tail behavior may surprise flood managers and citizens, as human intuition tends to expect light tail behavior, and the heaviness of the tails is very difficult to predict, which may lead to unnecessarily high flood damage.
Despite its high importance, the literature on the heavy tail behavior of flood distributions is rather fragmented.
In this review, we provide a coherent overview of the processes causing heavy flood tails and the implications for science and practice.
Specifically, we propose nine hypotheses on the mechanisms causing heavy tails in flood peak distributions related to processes in the atmosphere, the catchment, and the river system.
We then discuss to which extent the current knowledge supports or contradicts these hypotheses.
We also discuss the statistical conditions for the emergence of heavy tail behavior based on derived distribution theory and relate them to the hypotheses and flood generation mechanisms.
We review the degree to which the heaviness of the tails can be predicted from process knowledge and data. Finally, we recommend further research toward testing the hypotheses and improving the prediction of heavy tails.
Thousands of glacier lakes have been forming behind natural dams in high mountains following glacier retreat since the early 20th century. Some of these lakes abruptly released pulses of water and sediment with disastrous downstream consequences. Yet it remains unclear whether the reported rise of these glacier lake outburst floods (GLOFs) has been fueled by a warming atmosphere and enhanced meltwater production, or simply a growing research effort. Here we estimate trends and biases in GLOF reporting based on the largest global catalog of 1,997 dated glacier-related floods in six major mountain ranges from 1901 to 2017. We find that the positive trend in the number of reported GLOFs has decayed distinctly after a break in the 1970s, coinciding with independently detected trend changes in annual air temperatures and in the annual number of field-based glacier surveys (a proxy of scientific reporting). We observe that GLOF reports and glacier surveys decelerated, while temperature rise accelerated in the past five decades. Enhanced warming alone can thus hardly explain the annual number of reported GLOFs, suggesting that temperature-driven glacier lake formation, growth, and failure are weakly coupled, or that outbursts have been overlooked. Indeed, our analysis emphasizes a distinct geographic and temporal bias in GLOF reporting, and we project that between two to four out of five GLOFs on average might have gone unnoticed in the early to mid-20th century. We recommend that such biases should be considered, or better corrected for, when attributing the frequency of reported GLOFs to atmospheric warming.
On 7 January 2020, an M-w 6.4 earthquake occurred in the northeastern Caribbean, a few kilometers offshore of the island of Puerto Rico. It was the mainshock of a complex seismic sequence, characterized by a large number of energetic earthquakes illuminating an east-west elongated area along the southwestern coast of Puerto Rico. Deformation fields constrained by Interferometric Synthetic Aperture Radar and Global Navigation Satellite System data indicate that the coseismic movements affected only the western part of the island. To assess the mainshock's source fault parameters, we combined the geodetically derived coseismic deformation with teleseismic waveforms using Bayesian inference. The results indicate a roughly east-west oriented fault, dipping northward and accommodating similar to 1.4 m of transtensional motion. Besides, the determined location and orientation parameters suggest an offshore continuation of the recently mapped North Boqueron Bay-Punta Montalva fault in southwest Puerto Rico. This highlights the existence of unmapped faults with moderate-to-large earthquake potential within the Puerto Rico region.
We investigate whether the distribution of maximum seasonal streamflow is significantly affected by catchment or climate state of the season/month ahead. We fit the Generalized Extreme Value (GEV) distribution to extreme seasonal streamflow for around 600 stations across Europe by conditioning the GEV location and scale parameters on 14 indices, which represent the season-ahead climate or catchment state. The comparison of these climate-informed models with the classical GEV distribution, with time-constant parameters, suggests that there is a substantial potential for seasonal forecasting of flood probabilities. The potential varies between seasons and regions. Overall, the season-ahead catchment wetness shows the highest potential, although climate indices based on large-scale atmospheric circulation, sea surface temperature or sea ice concentration also show some skill for certain regions and seasons. Spatially coherent patterns and a substantial fraction of climate-informed models are promising signs towards early alerts to increase flood preparedness already a season ahead.
Efforts have been made in the past to enhance building exposure models on a regional scale with increasing spatial resolutions by integrating different data sources. This work follows a similar path and focuses on the downscaling of the existing SARA exposure model that was proposed for the residential building stock of the communes of Valparaiso and Vina del Mar (Chile). Although this model allowed great progress in harmonising building classes and characterising their differential physical vulnerabilities, it is now outdated, and in any case, it is spatially aggregated over large administrative units. Hence, to more accurately consider the impact of future earthquakes on these cities, it is necessary to employ more reliable exposure models. For such a purpose, we propose updating this existing model through a Bayesian approach by integrating ancillary data that has been made increasingly available from Volunteering Geo-Information (VGI) activities. Its spatial representation is also optimised in higher resolution aggregation units that avoid the inconvenience of having incomplete building-by-building footprints. A worst-case earthquake scenario is presented to calculate direct economic losses and highlight the degree of uncertainty imposed by exposure models in comparison with other parameters used to generate the seismic ground motions within a sensitivity analysis. This example study shows the great potential of using increasingly available VGI to update worldwide building exposure models as well as its importance in scenario-based seismic risk assessment.
It is widely recognized that collisional mountain belt topography is generated by crustal thickening and lowered by river bedrock erosion, linking climate and tectonics(1-4). However, whether surface processes or lithospheric strength control mountain belt height, shape and longevity remains uncertain. Additionally, how to reconcile high erosion rates in some active orogens with long-term survival of mountain belts for hundreds of millions of years remains enigmatic. Here we investigate mountain belt growth and decay using a new coupled surface process(5,6) and mantle-scale tectonic model(7). End-member models and the new non-dimensional Beaumont number, Bm, quantify how surface processes and tectonics control the topographic evolution of mountain belts, and enable the definition of three end-member types of growing orogens: type 1, non-steady state, strength controlled (Bm > 0.5); type 2, flux steady state(8), strength controlled (Bm approximate to 0.4-0.5); and type 3, flux steady state, erosion controlled (Bm < 0.4). Our results indicate that tectonics dominate in Himalaya-Tibet and the Central Andes (both type 1), efficient surface processes balance high convergence rates in Taiwan (probably type 2) and surface processes dominate in the Southern Alps of New Zealand (type 3). Orogenic decay is determined by erosional efficiency and can be subdivided into two phases with variable isostatic rebound characteristics and associated timescales. The results presented here provide a unified framework explaining how surface processes and lithospheric strength control the height, shape, and longevity of mountain belts.
Seismology, like many scientific fields, e.g., music information retrieval and speech signal pro- cessing, is experiencing exponential growth in the amount of data acquired by modern seismo- logical networks. In this thesis, I take advantage of the opportunities offered by "big data" and by the methods developed in the areas of music information retrieval and machine learning to predict better the ground motion generated by earthquakes and to study the properties of the surface layers of the Earth. In order to better predict seismic ground motions, I propose two approaches based on unsupervised deep learning methods, an autoencoder network and Generative Adversarial Networks. The autoencoder technique explores a massive amount of ground motion data, evaluates the required parameters, and generates synthetic ground motion data in the Fourier amplitude spectra (FAS) domain. This method is tested on two synthetic datasets and one real dataset. The application on the real dataset shows that the substantial information contained within the FAS data can be encoded to a four to the five-dimensional manifold. Consequently, only a few independent parameters are required for efficient ground motion prediction. I also propose a method based on Conditional Generative Adversarial Networks (CGAN) for simulating ground motion records in the time-frequency and time domains. CGAN generates the time-frequency domains based on the parameters: magnitude, distance, and shear wave velocities to 30 m depth (VS30). After generating the amplitude of the time-frequency domains using the CGAN model, instead of classical conventional methods that assume the amplitude spectra with a random phase spectrum, the phase of the time-frequency domains is recovered by minimizing the observed and reconstructed spectrograms. In the second part of this dissertation, I propose two methods for the monitoring and characterization of near-surface materials and site effect analyses. I implement an autocorrelation function and an interferometry method to monitor the velocity changes of near-surface materials resulting from the Kumamoto earthquake sequence (Japan, 2016). The observed seismic velocity changes during the strong shaking are due to the non-linear response of the near-surface materials. The results show that the velocity changes lasted for about two months after the Kumamoto mainshock. Furthermore, I used the velocity changes to evaluate the in-situ strain-stress relationship. I also propose a method for assessing the site proxy "VS30" using non-invasive analysis. In the proposed method, a dispersion curve of surface waves is inverted to estimate the shear wave velocity of the subsurface. This method is based on the Dix-like linear operators, which relate the shear wave velocity to the phase velocity. The proposed method is fast, efficient, and stable. All of the methods presented in this work can be used for processing "big data" in seismology and for the analysis of weak and strong ground motion data, to predict ground shaking, and to analyze site responses by considering potential time dependencies and nonlinearities.
The 100 km wide Merida Andes extend from the Colombian/Venezuelan border to the Coastal Cordillera. The mountain chain and its associated major strike-slip fault systems in western Venezuela formed due to oblique convergence of the Caribbean with the South American Plates and the north-eastwards expulsion of the North Andean Block. Due to the limited knowledge of lithospheric structures related to the formation of the Merida Andes research projects have been developed to illuminate this zone with deep geophysical data. In this study, we present three-dimensional inversion of broadband magnetotelluric data, collected along a 240 km long profile crossing the Merida Andes and the Maracaibo and Barinas-Apure foreland basins. The distribution of the stations limits resolution of the model to off-profile features. Combining 3D inversion of synthetic data sets derived from 3D modelling with 3D inversion of measured data, we could derive a 10 to 15 km wide corridor with good lateral resolution to develop hypotheses about the origin of deep-reaching anomalies of high electrical conductivity. The Merida Andes appear generally as electrically resistive structures, separated by anomalies associated with the most important fault systems of the region, the Bocono and Valera faults. Sensitivity tests suggest that the Valera Fault reaches to depths of up to 12 km and the Bocono Fault to more than 35 km depth. Both structures are connected to a sizeable conductor located east of the profile at 12-15 km depth. We propose that the high conductivity associated with this off-profile conductor may be related to the detachment of the Trujillo Block. We also identified a conductive zone that correlates spatially with the location of a gravity low, possibly representing a SE tilt of the Maracaibo Triangular Block under the mountain chain to great depths (>30 km). The relevance of these tectonic blocks in our models at crustal depths seems to be consistent with proposed theories that describe the geodynamics of western Venezuela as dominated by floating blocks or orogens. Our results stress the importance of the Trujillo Block for the current tectonic evolution of western Venezuela and confirm the relevance of the Bocono Fault carrying deformation to the lower crust and upper mantle. The Barinas-Apure and the Maracaibo sedimentary basins are imaged as electrically conductive with depths of 4 to 5 km and 5 to 10 km, respectively. The Barinas-Apure basin is imaged as a simple 1D structure, in contrast to the Maracaibo Basin, where a series of conductive and resistive bodies could be related to active deformation causing the juxtaposition of older geological formations and younger basin sediments.
The Andes are a ~7000 km long N-S trending mountain range developed along the South American western continental margin. Driven by the subduction of the oceanic Nazca plate beneath the continental South American plate, the formation of the northern and central parts of the orogen is a type case for a non-collisional orogeny. In the southern Central Andes (SCA, 29°S-39°S), the oceanic plate changes the subduction angle between 33°S and 35°S from almost horizontal (< 5° dip) in the north to a steeper angle (~30° dip) in the south. This sector of the Andes also displays remarkable along- and across- strike variations of the tectonic deformation patterns. These include a systematic decrease of topographic elevation, of crustal shortening and foreland and orogenic width, as well as an alternation of the foreland deformation style between thick-skinned and thin-skinned recorded along- and across the strike of the subduction zone. Moreover, the SCA are a very seismically active region. The continental plate is characterized by a relatively shallow seismicity (< 30 km depth) which is mainly focussed at the transition from the orogen to the lowland areas of the foreland and the forearc; in contrast, deeper seismicity occurs below the interiors of the northern foreland. Additionally, frequent seismicity is also recorded in the shallow parts of the oceanic plate and in a sector of the flat slab segment between 31°S and 33°S. The observed spatial heterogeneity in tectonic and seismic deformation in the SCA has been attributed to multiple causes, including variations in sediment thickness, the presence of inherited structures and changes in the subduction angle of the oceanic slab. However, there is no study that inquired the relationship between the long-term rheological configuration of the SCA and the spatial deformation patterns. Moreover, the effects of the density and thickness configuration of the continental plate and of variations in the slab dip angle in the rheological state of the lithosphere have been not thoroughly investigated yet. Since rheology depends on composition, pressure and temperature, a detailed characterization of the compositional, structural and thermal fields of the lithosphere is needed. Therefore, by using multiple geophysical approaches and data sources, I constructed the following 3D models of the SCA lithosphere: (i) a seismically-constrained structural and density model that was tested against the gravity field; (ii) a thermal model integrating the conversion of mantle shear-wave velocities to temperature with steady-state conductive calculations in the uppermost lithosphere (< 50 km depth), validated by temperature and heat-flow measurements; and (iii) a rheological model of the long-term lithospheric strength using as input the previously-generated models.
The results of this dissertation indicate that the present-day thermal and rheological fields of the SCA are controlled by different mechanisms at different depths. At shallow depths (< 50 km), the thermomechanical field is modulated by the heterogeneous composition of the continental lithosphere. The overprint of the oceanic slab is detectable where the oceanic plate is shallow (< 85 km depth) and the radiogenic crust is thin, resulting in overall lower temperatures and higher strength compared to regions where the slab is steep and the radiogenic crust is thick. At depths > 50 km, largest temperatures variations occur where the descending slab is detected, which implies that the deep thermal field is mainly affected by the slab dip geometry.
The outcomes of this thesis suggests that long-term thermomechanical state of the lithosphere influences the spatial distribution of seismic deformation. Most of the seismicity within the continental plate occurs above the modelled transition from brittle to ductile conditions. Additionally, there is a spatial correlation between the location of these events and the transition from the mechanically strong domains of the forearc and foreland to the weak domain of the orogen. In contrast, seismicity within the oceanic plate is also detected where long-term ductile conditions are expected. I therefore analysed the possible influence of additional mechanisms triggering these earthquakes, including the compaction of sediments in the subduction interface and dehydration reactions in the slab. To that aim, I carried out a qualitative analysis of the state of hydration in the mantle using the ratio between compressional- and shear-wave velocity (vp/vs ratio) from a previous seismic tomography. The results from this analysis indicate that the majority of the seismicity spatially correlates with hydrated areas of the slab and overlying continental mantle, with the exception of the cluster within the flat slab segment. In this region, earthquakes are likely triggered by flexural processes where the slab changes from a flat to a steep subduction angle.
First-order variations in the observed tectonic patterns also seem to be influenced by the thermomechanical configuration of the lithosphere. The mechanically strong domains of the forearc and foreland, due to their resistance to deformation, display smaller amounts of shortening than the relatively weak orogenic domain. In addition, the structural and thermomechanical characteristics modelled in this dissertation confirm previous analyses from geodynamic models pointing to the control of the observed heterogeneities in the orogen and foreland deformation style. These characteristics include the lithospheric and crustal thickness, the presence of weak sediments and the variations in gravitational potential energy.
Specific conditions occur in the cold and strong northern foreland, which is characterized by active seismicity and thick-skinned structures, although the modelled crustal strength exceeds the typical values of externally-applied tectonic stresses. The additional mechanisms that could explain the strain localization in a region that should resist deformation are: (i) increased tectonic forces coming from the steepening of the slab and (ii) enhanced weakening along inherited structures from pre-Andean deformation events. Finally, the thermomechanical conditions of this sector of the foreland could be a key factor influencing the preservation of the flat subduction angle at these latitudes of the SCA.
The stabilizing properties of mineral-organic carbon (OC) interactions have been studied in many soil environments (temperate soils, podzol lateritic soils, and paddy soils). Recently, interest in their role in permafrost regions is increasing as permafrost was identified as a hotspot of change. In thawing ice-rich permafrost regions, such as the Yedoma domain, 327-466 Gt of frozen OC is buried in deep sediments. Interactions between minerals and OC are important because OC is located very near the mineral matrix. Mineral surfaces and elements could mitigate recent and future greenhouse gas emissions through physical and/or physicochemical protection of OC. The dynamic changes in redox and pH conditions associated with thermokarst lake formation and drainage trigger metal-oxide dissolution and precipitation, likely influencing OC stabilization and microbial mineralization. However, the influence of thermokarst processes on mineral-OC interactions remains poorly constrained. In this study, we aim to characterize Fe, Mn, Al, and Ca minerals and their potential protective role for OC. Total and selective extractions were used to assess the crystalline and amorphous oxides or complexed metal pools as well as the organic acids found within these pools. We analyzed four sediment cores from an ice-rich permafrost area in Central Yakutia, which were drilled (i) in undisturbed Yedoma uplands, (ii) beneath a recent lake formed within Yedoma deposits, (iii) in a drained thermokarst lake basin, and (iv) beneath a mature thermokarst lake from the early Holocene period. We find a decrease in the amount of reactive Fe, Mn, Al, and Ca in the deposits on lake formation (promoting reduction reactions), and this was largely balanced by an increase in the amount of reactive metals in the deposits on lake drainage (promoting oxidation reactions). We demonstrate an increase in the metal to C molar ratio on thermokarst process, which may indicate an increase in metal-C bindings and could provide a higher protective role against microbial mineralization of organic matter. Finally, we find that an increase in mineral-OC interactions corresponded to a decrease in CO2 and CH4 gas emissions on thermokarst process. Mineral-OC interactions could mitigate greenhouse gas production from permafrost thaw as soon as lake drainage occurs.
Boreal forests cover over half of the global permafrost area and protect underlying permafrost. Boreal forest development, therefore, has an impact on permafrost evolution, especially under a warming climate.
Forest disturbances and changing climate conditions cause vegetation shifts and potentially destabilize the carbon stored within the vegetation and permafrost. Disturbed permafrost-forest ecosystems can develop into a dry or swampy bush- or grasslands, shift toward broadleaf- or evergreen needleleaf-dominated forests, or recover to the pre-disturbance state.
An increase in the number and intensity of fires, as well as intensified logging activities, could lead to a partial or complete ecosystem and permafrost degradation. We study the impact of forest disturbances (logging, surface, and canopy fires) on the thermal and hydrological permafrost conditions and ecosystem resilience.
We use a dynamic multilayer canopy-permafrost model to simulate different scenarios at a study site in eastern Siberia. We implement expected mortality, defoliation, and ground surface changes and analyze the interplay between forest recovery and permafrost. We find that forest loss induces soil drying of up to 44%, leading to lower active layer thicknesses and abrupt or steady decline of a larch forest, depending on disturbance intensity.
Only after surface fires, the most common disturbances, inducing low mortality rates, forests can recover and overpass pre-disturbance leaf area index values. We find that the trajectory of larch forests after surface fires is dependent on the precipitation conditions in the years after the disturbance. Dryer years can drastically change the direction of the larch forest development within the studied period.
Reconstructing thermal histories in thrust belts is commonly used to infer the age and rates of thrusting and hence the driving mechanisms of orogenesis.
In areas where ancient basins have been incorporated into the orogenic wedge, a quantitative reconstruction of the thermal history helps distinguish among potential mechanisms responsible for heating events.
We present such a reconstruction for the Ischigualasto-Villa Union basin in the western Pampean Ranges of Argentina, where Triassic rifting and late Cretaceous-Cenozoic retroarc foreland basin development has been widely documented, including Miocene flat-slab subduction.
We report results of organic and inorganic thermal indicators acquired along three stratigraphic sections, including vitrinite reflectance and X-ray diffractometry in claystones and new thermochronological [(apatite fission-track and apatite and zircon [U-Th]/He)] analyses.
Despite up to 5 km-thick Cenozoic overburden and unlike previously thought, the thermal peak in the basin is not due to Cenozoic burial but occurred in the Triassic, associated with a high heat flow of up to 90 mWm(-2) and <2 km of burial, which heated the base of the Triassic strata to similar to 160 degrees C. Following exhumation, attested by the development of an unconformity between the Triassic and Late-Cretaceous-Cenozoic sequences, Cenozoic re-burial increased the temperature to similar to 110 degrees C at the base of the Triassic section and only similar to 50 degrees C 7 km upsection, suggesting a dramatic decrease in the thermal gradient.
The onset of Cenozoic cooling occurred at similar to 10(-8) Ma, concomitant with sediment accumulation and thus preceding the latest Miocene onset of thrusting that has been independently documented by stratigraphic-cross-cutting relationships.
We argue that the onset of cooling is associated with lithospheric refrigeration following establishment of flat-slab subduction, leading to the eastward displacement of the asthenospheric wedge beneath the South American plate.
Our study places time and temperature constraints on flat-slab cooling that calls for a careful interpretation of exhumation signals in thrustbelts inferred from thermochronology only.
The spatial pattern of extreme precipitation from 40 years of gauge data in the central Himalaya
(2022)
The topography of the Himalaya exerts a substantial control on the spatial distribution of monsoonal rainfall, which is a vital water source for the regional economy and population. But the occurrence of short-lived and high-intensity precipitation results in socio-economic losses. This study relies on 40 years of daily data from 204 ground stations in Nepal to derive extreme precipitation thresholds, amounts, and days at the 95th percentile. We additionally determine the precipitation magnitude-frequency relation. We observe that extreme precipitation amounts follow an almost uniform band parallel to topographic contour lines in the southern Himalaya mountains in central and eastern Nepal but not in western Nepal. The relationship of extreme precipitation indices with topographic relief shows that extreme precipitation thresholds decrease with increasing elevation, but extreme precipitation days increase in higher elevation areas. Furthermore, stations above 1 km elevation exhibit a power-law relation in the rainfall magnitude-frequency framework. Stations at higher elevations generally have lower values of power-law exponents than low elevation areas. This suggests a fundamentally different behaviour of the rainfall distribution and an increased occurrence of extreme rainfall storms in the high elevation areas of Nepal.
Agricultural production worldwide has been increasing in the last decades at a very fast pace and with it the waste generation. Livestock activities are one of the largest producers of residues in the agricultural sector and contribute greatly to climate change. The present chapter gives an introduction and an in-depth analysis of the waste management of livestock for the conversion in a circular agriculture and economy based on research and experience in the sector conducted in the last decades. The conversion of animal waste into energy generation is an opportunity for farmers to obtain additional economic benefits, while contributing to the environment by preventing the release of GHGs into the atmosphere. The use of animal waste for energy generation through anaerobic digestion is a progressive technique and is being widely accepted in Europe, where Germany is the leading country in the use of biogas plants for energy production among others in the European Union. Economically speaking, the livestock industry faces the challenge of converting its production into a clean and more profitable production. The goal of this chapter is to analyze the economic benefit as well as the environmental contribution and future challenges of the use of livestock waste in the biorefineries sector from different perspectives, based on an intensive literature review. This review is accompanied by a geospatial analysis component, mapping biogas reactor hotspots and clusters in Germany, by means of methods of spatial statistics as analysis methods as kernel density estimations (KDE) and K-means clustering, based on volunteer geographic data. The applied methods easily can be transferred to other regions and allow a quick macroscopic overview over existing biogas reactors; furthermore, an identification of cluster and hotspots with a high biogas potential, that in a subsequent step can be analyzed in depth in larger scales.
The Walker Circulation (WC) is an east-west trending band of atmospheric circulation cells along the equator and the predominant controller of heat and moisture transport in the tropics. Its variability is closely linked to the sea-surface temperature (SST) changes across the Pacific, the Indian and the Atlantic Oceans and can have pronounced effects on the humidity regimes of the adjacent continents. In recent years, the evolution of the WC during the Plioand Pleistocene epochs has been intensely studied in the context of the effectiveness of the tropics in modulating global climate change (e.g., the intensification of Northern Hemisphere glaciation). However, the onset of the modern WC pattern as well as its global impact during the Plioand Pleistocene is controversially assessed in the literature. For its onset, previous studies have suggested dates ranging between 2.4 and 0.8 million years ago (Myr), while its argued impact ranges from crucially influencing the increase of Northern Hemisphere ice sheet growth by channelling heat and moisture from the tropics into the high latitudes to having no effect on global ice volume changes. In order to achieve a comprehensive understanding of the spatiotemporal evolution of the WC during this time frame, we statistically analysed 30 globally distributed SST records covering the low and high latitudes between 3.5 and 1.5 Myr, encompassing the Late Pliocene to Early Pleistocene. We utilized a statistical change-point regression model to determine significant change points in the SST evolution of the (sub)-tropics and high latitudes that potentially relate to changes in the WC. We find that the WC experienced a multifaceted evolution between the Late Pliocene and the Early Pleistocene with significant transitional steps at-2.7 and-2.1 Ma. Our results suggest after the Late Pliocene, a pre-modern WC set in, which was characterized by a progressively strengthened Pacific Walker Cell alongside a weakened Indian Ocean Walker Cell. This change was potentially triggered by the constriction of the Indonesian seaway, an important transmitter between the Pacific and Indian Ocean. The ensuing mode of the WC intensified until-2.1 Myr, when SST values around the global scale signalled a progressive strengthening of the Indian Walker Cell in phase with the progressive strengthening of the Pacific and Atlantic Cells. Our findings indicate that a shift from a pre-modern to a modern-like WC potentially only occurred during the mid-Pleistocene.