Refine
Year of publication
Document Type
- Article (2623)
- Doctoral Thesis (409)
- Postprint (138)
- Other (75)
- Review (51)
- Preprint (17)
- Monograph/Edited Volume (13)
- Conference Proceeding (12)
- Habilitation Thesis (11)
- Master's Thesis (6)
Language
- English (3357) (remove)
Keywords
- climate change (50)
- Holocene (44)
- erosion (27)
- Himalaya (26)
- permafrost (26)
- Climate change (23)
- remote sensing (23)
- Tibetan Plateau (22)
- Earthquake source observations (21)
- Pollen (21)
Institute
- Institut für Geowissenschaften (3357) (remove)
Analysis and modeling of transient earthquake patterns and their dependence on local stress regimes
(2015)
Investigations in the field of earthquake triggering and associated interactions, which includes aftershock triggering as well as induced seismicity, is important for seismic hazard assessment due to earthquakes destructive power. One of the approaches to study earthquake triggering and their interactions is the use of statistical earthquake models, which are based on knowledge of the basic seismicity properties, in particular, the magnitude distribution and spatiotemporal properties of the triggered events.
In my PhD thesis I focus on some specific aspects of aftershock properties, namely, the relative seismic moment release of the aftershocks with respect to the mainshocks; the spatial correlation between aftershock occurrence and fault deformation; and on the influence of aseismic transients on the aftershock parameter estimation. For the analysis of aftershock sequences I choose a statistical approach, in particular, the well known Epidemic Type Aftershock Sequence (ETAS) model, which accounts for the input of background and triggered seismicity. For my specific purposes, I develop two ETAS model modifications in collaboration with Sebastian Hainzl. By means of this approach, I estimate the statistical aftershock parameters and performed simulations of aftershock sequences as well.
In the case of seismic moment release of aftershocks, I focus on the ratio of cumulative seismic moment release with respect to the mainshocks. Specifically, I investigate the ratio with respect to the focal mechanism of the mainshock and estimate an effective magnitude, which represents the cumulative aftershock energy (similar to Bath's law, which defines the average difference between mainshock and the largest aftershock magnitudes). Furthermore, I compare the observed seismic moment ratios with the results of the ETAS simulations. In particular, I test a restricted ETAS (RETAS) model which is based on results of a clock advanced model and static stress triggering.
To analyze spatial variations of triggering parameters I focus in my second approach on the aftershock occurrence triggered by large mainshocks and the study of the aftershock parameter distribution and their spatial correlation with the coseismic/postseismic slip and interseismic locking. To invert the aftershock parameters I improve the modified ETAS (m-ETAS) model, which is able to take the extension of the mainshock rupture into account. I compare the results obtained by the classical approach with the output of the m-ETAS model.
My third approach is concerned with the temporal clustering of seismicity, which might not only be related to earthquake-earthquake interactions, but also to a time-dependent background rate, potentially biasing the parameter estimations. Thus, my coauthors and I also applied a modification of the ETAS model, which is able to take into account time-dependent background activity. It can be applicable for two different cases: when an aftershock catalog has a temporal incompleteness or when the background seismicity rate changes with time, due to presence of aseismic forces.
An essential part of any research is the testing of the developed models using observational data sets, which are appropriate for the particular study case. Therefore, in the case of seismic moment release I use the global seismicity catalog. For the spatial distribution of triggering parameters I exploit two aftershock sequences of the Mw8.8 2010 Maule (Chile) and Mw 9.0 2011 Tohoku (Japan) mainshocks. In addition, I use published geodetic slip models of different authors. To test our ability to detect aseismic transients my coauthors and I use the data sets from Western Bohemia (Central Europe) and California.
Our results indicate that:
(1) the seismic moment of aftershocks with respect to mainshocks depends on the static stress changes and is maximal for the normal, intermediate for thrust and minimal for strike-slip stress regimes, where the RETAS model shows a good correspondence with the results;
(2) The spatial distribution of aftershock parameters, obtained by the m-ETAS model, shows anomalous values in areas of reactivated crustal fault systems. In addition, the aftershock density is found to be correlated with coseismic slip gradient, afterslip, interseismic coupling and b-values. Aftershock seismic moment is positively correlated with the areas of maximum coseismic slip and interseismically locked areas. These correlations might be related to the stress level or to material properties variations in space;
(3) Ignoring aseismic transient forcing or temporal catalog incompleteness can lead to the significant under- or overestimation of the underlying trigger parameters. In the case when a catalog is complete, this method helps to identify aseismic sources.
Volcanoes are one of the Earth’s most dynamic zones and responsible for many changes in our planet. Volcano seismology aims to provide an understanding of the physical processes in volcanic systems and anticipate the style and timing of eruptions by analyzing the seismic records. Volcanic tremor signals are usually observed in the seismic records before or during volcanic eruptions. Their analysis contributes to evaluate the evolving volcanic activity and potentially predict eruptions. Years of continuous seismic monitoring now provide useful information for operational eruption forecasting. The continuously growing amount of seismic recordings, however, poses a challenge for analysis, information extraction, and interpretation, to support timely decision making during volcanic crises. Furthermore, the complexity of eruption processes and precursory activities makes the analysis challenging.
A challenge in studying seismic signals of volcanic origin is the coexistence of transient signal swarms and long-lasting volcanic tremor signals. Separating transient events from volcanic tremors can, therefore, contribute to improving our understanding of the underlying physical processes. Some similar issues (data reduction, source separation, extraction, and classification) are addressed in the context of music information retrieval (MIR). The signal characteristics of acoustic and seismic recordings comprise a number of similarities. This thesis is going beyond classical signal analysis techniques usually employed in seismology by exploiting similarities of seismic and acoustic signals and building the information retrieval strategy on the expertise developed in the field of MIR.
First, inspired by the idea of harmonic–percussive separation (HPS) in musical signal processing, I have developed a method to extract harmonic volcanic tremor signals and to detect transient events from seismic recordings. This provides a clean tremor signal suitable for tremor investigation along with a characteristic function suitable for earthquake detection. Second, using HPS algorithms, I have developed a noise reduction technique for seismic signals. This method is especially useful for denoising ocean bottom seismometers, which are highly contaminated by noise. The advantage of this method compared to other denoising techniques is that it doesn’t introduce distortion to the broadband earthquake waveforms, which makes it reliable for different applications in passive seismological analysis. Third, to address the challenge of extracting information from high-dimensional data and investigating the complex eruptive phases, I have developed an advanced machine learning model that results in a comprehensive signal processing scheme for volcanic tremors. Using this method seismic signatures of major eruptive phases can be automatically detected. This helps to provide a chronology of the volcanic system. Also, this model is capable to detect weak precursory volcanic tremors prior to the eruption, which could be used as an indicator of imminent eruptive activity. The extracted patterns of seismicity and their temporal variations finally provide an explanation for the transition mechanism between eruptive phases.
Volcanic tremor extraction and earthquake detection using music information retrieval algorithms
(2021)
Volcanic tremor signals are usually observed before or during volcanic eruptions and must be monitored to evaluate the volcanic activity. A challenge in studying seismic signals of volcanic origin is the coexistence of transient signal swarms and long-lasting volcanic tremor signals. Separating transient events from volcanic tremors can, therefore, contrib-ute to improving upon our understanding of the underlying physical processes. Exploiting the idea of harmonic-percussive separation in musical signal processing, we develop a method to extract the harmonic volcanic tremor signals and to detect tran-sient events from seismic recordings. Based on the similarity properties of spectrogram frames in the time-frequency domain, we decompose the signal into two separate spec-trograms representing repeating (harmonic) and nonrepeating (transient) patterns, which correspond to volcanic tremor signals and earthquake signals, respectively. We reconstruct the harmonic tremor signal in the time domain from the complex spectrogram of the repeating pattern by only considering the phase components for the frequency range in which the tremor amplitude spectrum is significantly contribut-ing to the energy of the signal. The reconstructed signal is, therefore, clean tremor signal without transient events. Furthermore, we derive a characteristic function suitable for the detection of tran-sient events (e.g., earthquakes) by integrating amplitudes of the nonrepeating spectro-gram over frequency at each time frame. Considering transient events like earthquakes, 78% of the events are detected for signal-to-noise ratio = 0.1 in our semisynthetic tests. In addition, we compared the number of detected earthquakes using our method for one month of continuous data recorded during the Holuhraun 2014-2015 eruption in Iceland with the bulletin presented in Agustsdottir et al. (2019). Our single station event detection algorithm identified 84% of the bulletin events. Moreover, we detected a total of 12,619 events, which is more than twice the number of the bulletin events.
Records from ocean bottom seismometers (OBSs) are highly contaminated by noise, which is much stronger
compared to data from most land stations, especially on the horizontal components. As a consequence, the high energy of the oceanic noise at frequencies below 1 Hz considerably complicates the analysis of the teleseismic earthquake signals recorded by OBSs.
Previous studies suggested different approaches to remove low-frequency noises from OBS recordings but mainly focused on the vertical component. The records of horizontal components, which are crucial for the application of many methods in passive seismological analysis of body and surface waves, could not be much improved in the teleseismic frequency band. Here we introduce a noise reduction method, which is derived from the harmonic–percussive separation algorithms used in Zali et al. (2021), in order to separate long-lasting narrowband signals from broadband transients in the OBS signal. This leads to significant noise reduction of OBS records on both the vertical and horizontal components and increases the earthquake signal-to-noise ratio (SNR) without distortion of the broadband earthquake waveforms. This is demonstrated through tests with synthetic data. Both SNR and cross-correlation coefficients showed significant improvements for different realistic noise realizations. The application of denoised signals in surface wave analysis and receiver functions is discussed through tests with synthetic and real data.
Records from ocean bottom seismometers (OBSs) are highly contaminated by noise, which is much stronger compared to data from most land stations, especially on the horizontal components. As a consequence, the high energy of the oceanic noise at frequencies below 1 Hz considerably complicates the analysis of the teleseismic earthquake signals recorded by OBSs.
Previous studies suggested different approaches to remove low-frequency noises from OBS recordings but mainly focused on the vertical component. The records of horizontal components, which are crucial for the application of many methods in passive seismological analysis of body and surface waves, could not be much improved in the teleseismic frequency band. Here we introduce a noise reduction method, which is derived from the harmonic–percussive separation algorithms used in Zali et al. (2021), in order to separate long-lasting narrowband signals from broadband transients in the OBS signal. This leads to significant noise reduction of OBS records on both the vertical and horizontal components and increases the earthquake signal-to-noise ratio (SNR) without distortion of the broadband earthquake waveforms. This is demonstrated through tests with synthetic data. Both SNR and cross-correlation coefficients showed significant improvements for different realistic noise realizations. The application of denoised signals in surface wave analysis and receiver functions is discussed through tests with synthetic and real data.
Modern anthropogenic forcing of atmospheric chemistry poses the question of how the Earth System will respond as thousands of gigatons of greenhouse gas are rapidly added to the atmosphere. A similar, albeit nonanthropogenic, situation occurred during the early Paleogene, when catastrophic release of carbon to the atmosphere triggered abrupt increase in global temperatures. The best documented of these events is the Paleocene-Eocene Thermal Maximum (PETM, ~55 Ma) when the magnitude of carbon addition to the oceans and atmosphere was similar to those expected for the future. This event initiated global warming, changes in hydrological cycles, biotic extinction and migrations. A recently proposed hypothesis concerning changes in marine ecosystems suggests that this global warming strongly influenced the shallow-water biosphere, triggering extinctions and turnover in the Larger Foraminifera (LF) community and the demise of corals. The successions from the Adriatic Carbonate Platform (SW Slovenia) represent an ideal location to test the hypothesis of a possible causal link between the PETM and evolution of shallow-water organisms because they record continuous sedimentation from the Late Paleocene to the Early Eocene and are characterized by a rich biota, especially LF, fundamental for detailed biostratigraphic studies. In order to reconstruct paleoenvironmental conditions during deposition, I focused on sedimentological analysis and paleoecological study of benthic assemblages. During the Late Paleocene-earliest Eocene, sedimentation occurred on a shallow-water carbonate ramp system characterized by enhanced nutrient levels. LF represent the common constituent of the benthic assemblages that thrived in this setting throughout the Late Paleocene to the Early Eocene. With detailed biostratigraphic and chemostratigraphic analyses documenting the most complete record to date available for the PETM event in a shallow-water marine environment, I correlated chemostratigraphically for the first time the evolution of LF with the δ¹³C curves. This correlation demonstrated that no major turnover in the LF communities occurred synchronous with the PETM; thus the evolution of LF was mainly controlled by endogenous biotic forces. The study of Late Thanetian metric-sized microbialite-coral mounds which developed in the middle part of the ramp, documented the first Cenozoic occurrence of microbially-cemented mounds. The development of these mounds, with temporary dominance of microbial communities over corals, suggest environmentally-triggered “phase shifts” related to frequent fluctuations of nutrient/turbidity levels during recurrent wet phases which preceding the extreme greenhouse conditions of the PETM. The paleoecological study of the coral community in the microbialites-coral mounds, the study of corals from Early Eocene platform from SW France, and a critical, extensive literature research of Late Paleocene – Early Eocene coral occurrences from the Tethys, the Atlantic, the Caribbean realms suggested that these corals types, even if not forming extensive reefs, are common in the biofacies as small isolated colonies, piles of rubble or small patch-reefs. These corals might have developed ‘alternative’ life strategies to cope with harsh conditions (high/fluctuating nutrients/turbidity, extreme temperatures, perturbation of aragonite saturation state) during the greenhouse times of the early Paleogene, representing a good fossil analogue to modern corals thriving close to their thresholds for survival. These results demonstrate the complexity of the biological responses to extreme conditions, not only in terms of temperature but also nutrient supply, physical disturbance and their temporal variability and oscillating character.
Upper Thanetian microbialite-coral mounds from the Adriatic Carbonate Platform (SW Slovenia) are described herein for the first time, representing an important case study of extensively microbially-cemented boundstones in the Early Paleogene. The mounds are constructed primarily by microbialites associated to small-sized coral colonies, forming metric bioconstructions in a mid-ramp setting. Detailed macroscopic and microscopic studies show that microbes are the major framework builders, playing a prominent role in the stabilization and growth of the mounds, with corals being the second most important component. Microbial carbonates represent up to 70% of the mounds, forming centimetric-thick crusts alternating with coral colonies. The microbial nature of the crusts is demonstrated by their growth form and internal microfabrics, showing accretionary, binding, and encrusting growth fabrics, often with gravity-defying geometries. Thin sections and polished slabs reveal a broad range of mesofabrics, with dense, structureless micrite (leiolite), laminated crusts (stromatolites), and clotted micritic masses (thrombolites). A first layer of micro- encrusters, including leiolites and thrombolites, occurs in cryptic habitats, whereas discontinuous stromatolites encrust the upper surface of corals. A second encrustation, the major mound construction phase, follows and is dominated by thrombolites, encrusting corals and other micro-encrusters. This sequence represents the basic constructional unit horizontally and vertically interlocked, in an irregular pattern, to form the mounds. The processes, which favored the deposition of these microbial carbonates, were mainly related to in situ precipitation, with minor evidences for grain agglutination and trapping processes. Scleractinian corals comprise moderately diversified community of small (centimetric) colonial, massive, platy encrusting, and branching forms. Coral colonies are distributed uniformly throughout the mounds without developing any ecological zonation. These features indicate that coral development remained at the pioneer stage throughout the mound growth. The spatial relationships between corals and microbialites, as well as the characteristics of microbial crusts and coral colonies, indicate a strong ecological competition between corals and microbes. A model for the evolution of the trophic structures during the mound growth is proposed, with changes in the paleoecology of the main bioconstructors triggered by frequent environmental perturbations. Turbidity and nutrient pressure, interpreted here as related to frequent recurrences of wet phases during the warm, humid climate of the Uppermost Thanetian, might have promoted temporary dominance of microbes over corals, causing rapid environmentally- driven "phase shifts" in the dominant biota.
The Paleocene-Eocene thermal maximum represents one of the most rapid and extreme warming events in the Cenozoic. Shallow-water stratigraphic sections from the Adriatic carbonate platform offer a rare opportunity to learn about the nature of Paleocene-Eocene thermal maximum and the effects on shallow-water ecosystems. We use carbon and oxygen isotope stratigraphy, in conjunction with detailed larger benthic foraminiferal biostratigraphy, to establish a high-resolution paleoclimatic record for the Paleocene-Eocene thermal maximum. A prominent negative excursion in delta C-13 curves of bulk-rock (similar to 1 parts per thousand-3 parts per thousand), matrix (similar to 4 parts per thousand), and foraminifera (similar to 6 parts per thousand) is interpreted as the carbon isotope excursion during the Paleocene-Eocene thermal maximum. The strongly C-13-depleted delta(1)d(3)C record of our shallow-marine carbonates compared to open-marine records could result from organic matter oxidation, suggesting intensified weathering, runoff, and organic matter flux.
The Ilerdian larger benthie foraminiferal turnover is documented in detail based on high-resolution correlation with the carbon isotopic excursion. The turnover is described as a two-step process, with the first step (early Ilerdian) marked by a rapid diversification of small alveolinids and nummulitids with weak adult dimorphism, possibly as adaptations to fluctuating Paleocene-Eocene thermal maximum nutrient levels, and a second step (middle Ilerdian) characterized by a further specific diversification, increase of shell size, and well-developed adult dimorphism. Within an evolutionary scheme controlled by long-term biological processes, we argue that high seawater temperatures could have stimulated the early Ilerdian rapid specific diversification. Together, these data help elucidate the effects of global warming and associated feedbacks in shallow-water ecosystems, and by inference, could serve as an assessment analog for future changes.
Today, diverse communities of zooxanthellate corals thrive, but do not build reef, under a wide range of environmental conditions. In these settings they inhabit natural bottom communities, sometimes forming patch-reefs, coral carpets and knobs. Episodes in the fossil record, characterized by limited coral-reef development but widespread occurrence of coral-bearing carbonates, may represent the fossil analogs of these non-reef building, zooxanthellate coral communities. If so, the study of these corals could have valuable implications for paleoenvironmental reconstructions. Here we focus on the evolution of early Paleogene corals as a fossil example of coral communities mainly composed by zooxanthellate corals (or likely zooxanthellate), commonly occurring within carbonate biofacies and with relatively high diversity but with a limited bioconstructional potential as testified by the reduced record of coral reefs. We correlate changes of bioconstructional potential and community compositions of these fossil corals with the main ecological/environmental conditions at that time. The early Paleogene greenhouse climate was characterized by relatively short pulses of warming with the most prominent occurring at the Paleocene-Eocene boundary (PETM event), associated with high weathering rates, nutrient fluxes, and pCO(2) levels. A synthesis of coral occurrences integrated with our data from the Adriatic Carbonate Platform (SW Slovenia) and the Minervois region (SW France), provides evidence for temporal changes in the reef-building capacity of corals associated with a shift in community composition toward forms adapted to tolerate deteriorating sea-water conditions. During the middle Paleocene coral-algal patch reefs and barrier reefs occurred from shallow-water settings, locally with reef-crest structures. A first shift can be traced from middle Paleocene to late Paleocene, with small coral-algal patch reefs and coral-bearing mounds development in shallow to intermediate water depths. In these mounds corals were highly subordinated as bioconstructors to other groups tolerant to higher levels of trophic resources (calcareous red algae, encrusting foraminifera, microbes, and sponges). A second shift occurred at the onset of the early Eocene with a further reduction of coral framework-building capacity. These coral communities mainly formed knobs in shallow-water, turbid settings associated with abundant foraminiferal deposits. We suggest that environmental conditions other than high temperature determined a combination of interrelated stressors that limited the coral-reef construction. A continuous enhancement of sediment load/nutrients combined with geochemical changes of ocean waters likely displaced corals as the main bioconstructors during the late Paleocene-early Eocene times. Nonetheless, these conditions did not affect the capacity of some corals to colonize the substrate, maintain biodiversity, and act as locally important carbonate-sediment producers, suggesting broad environmental tolerance limits of various species of corals. The implications of this study include clues as to how both ancient and modern zooxanthellate corals could respond to changing climate.
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.
The Cretaceous units exposed in the northwestern segment of the Colombian Andes preserve the record of extensional and compressional tectonics prior to the collision with Caribbean oceanic terranes. We integrated field, stratigraphic, sedimentary provenance, whole rock geochemistry, Nd isotopes and U-Pb zircon data to understand the Cretaceous tectonostratigraphic and magmatic record of the Colombian Andes. The results suggest that several sedimentary successions including the Abejorral Fm. were deposited on top of the continental basement in an Early Cretaceous backarc basin (150-100 Ma). Between 120 and 100 Ma, the appearance of basaltic and andesitic magmatism (similar to 115-100 Ma), basin deepening, and seafloor spreading were the result of advanced stages of backarc extension. A change to compressional tectonics took place during the Late Cretaceous (100-80 Ma). During this compressional phase, the extended blocks were reincorporated into the margin, closing the former Early Cretaceous backarc basin. Subsequently, a Late Cretaceous volcanic arc was built on the continental margin: as a result, the volcanic rocks of the Quebradagrande Complex were unconformably deposited on top of the faulted and folded rocks of the Abejorral Fm. Between the Late Cretaceous and the Paleocene (80-60 Ma), an arc-continent collision between the Caribbean oceanic plateau and the South-American continental margin deformed the rocks of the Quebradagrande Complex and shut-down the active volcanic arc. Our results suggest an Early Cretaceous extensional event followed by compressional tectonics prior to the collision with the Caribbean oceanic plateau. (C) 2019 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.
We present new U-Pb LA-ICP-MS data from the Central Andean foreland basins combined with new and published stratigraphic information in order to reconstruct the Miocene fragmentation of the Andean foreland between 26 and 28 degrees S. The disruption of this foreland basin and the subsequent development of elevated intermountain basins have been the focus of several studies. However, the absence of temporal constraints in the Miocene to Pliocene sedimentary record of the low elevation Choromoro and Tucuman foreland basins has presented an obstacle for precise paleogeographic reconstructions. We describe 11 discontinuous stratigraphic sections and use the U-Pb LA-ICP-MS method to date 10 pyroclastic-bearing sediments in order to reconstruct the stratigraphic evolution of the Choromoro and Tucuman basins. We combine our results with published strati graphic and thermochronologic data from adjacent basins to present a refined Miocene paleogeographic model. In a first stage, a continuous Early Miocene foreland lacustrine basin developed, filling up the preexisting Paleogene topography. The second stage is characterized by basin unroofing around similar to 12 Ma; the easily eroded sedimentary cover was removed, leading to the uplift of the underlying basement rocks and the segmentation of the lacustrine system. In the third stage, relief increase took place after similar to 6 Ma due to the low erodibility of the basement blocks; as a result, stable fluvial systems developed. Progressive relief development caused pronounced unconformities in the basins and the development of proximal fluvial-gravitational depositional systems after 3 Ma. This model emphasizes on the relations between tectonics, climate, and erodibility, and their control on the evolution of the depositional systems and relief.
Marked along-strike changes in stratigraphy, mountain belt morphology, basement exhumation, and deformation styles characterize the Andean retroarc; these changes have previously been related to spatiotemporal variations in the subduction angle. We modeled new apatite fission track and apatite (U-Th-Sm)/He data from nine ranges located between 26 degrees S and 28 degrees S. Using new and previously published data, we constructed a Cretaceous to Pliocene paleogeographic model that delineates a four-stage tectonic evolution: extensional tectonics during the Cretaceous (120-75 Ma), the formation of a broken foreland basin between 55 and 30 Ma, reheating due to burial beneath sedimentary rocks (18-13 Ma), and deformation, exhumation, and surface uplift during the Late Miocene and the Pliocene (13-3 Ma). Our model highlights how preexisting upper plate structures control the deformation patterns of broken foreland basins. Because retroarc deformation predates flat-slab subduction, we propose that slab anchoring may have been the precursor of Eocene-Oligocene compression in the Andean retroarc. Our model challenges models which consider broken foreland basins and retroarc deformation in the NW Argentinian Andes to be directly related to Miocene flat subduction.
Marked along-strike changes in stratigraphy, mountain belt morphology, basement exhumation, and deformation styles characterize the Andean retroarc; these changes have previously been related to spatiotemporal variations in the subduction angle. We modeled new apatite fission track and apatite (U-Th-Sm)/He data from nine ranges located between 26 degrees S and 28 degrees S. Using new and previously published data, we constructed a Cretaceous to Pliocene paleogeographic model that delineates a four-stage tectonic evolution: extensional tectonics during the Cretaceous (120-75 Ma), the formation of a broken foreland basin between 55 and 30 Ma, reheating due to burial beneath sedimentary rocks (18-13 Ma), and deformation, exhumation, and surface uplift during the Late Miocene and the Pliocene (13-3 Ma). Our model highlights how preexisting upper plate structures control the deformation patterns of broken foreland basins. Because retroarc deformation predates flat-slab subduction, we propose that slab anchoring may have been the precursor of Eocene-Oligocene compression in the Andean retroarc. Our model challenges models which consider broken foreland basins and retroarc deformation in the NW Argentinian Andes to be directly related to Miocene flat subduction.
Seismological and seismotectonic analysis of the northwestern Argentine Central Andean foreland
(2020)
After a severe M W 5.7 earthquake on October 17, 2015 in El Galpón in the province of Salta NW Argentina, I installed a local seismological network around the estimated epicenter. The network covered an area characterized by inherited Cretaceous normal faults and neotectonic faults with unknown recurrence intervals, some of which may have been reactivated normal faults. The 13 three-component seismic stations recorded data continuously for 15 months.
The 2015 earthquake took place in the Santa Bárbara System of the Andean foreland, at about 17km depth. This region is the easternmost morphostructural region of the central Andes. As a part of the broken foreland, it is bounded to the north by the Subandes fold-and-thrust belt and the Sierras Pampeanas to the south; to the east lies the Chaco-Paraná basin.
A multi-stage morphotectonic evolution with thick-skinned basement uplift and coeval thin-skinned deformation in the intermontane basins is suggested for the study area. The release of stresses associated with the foreland deformation can result in strong earthquakes, as the study area is known for recurrent and historical, destructive earthquakes. The available continuous record reaches back in time, when the strongest event in 1692 (magnitude 7 or intensity IX) destroyed the city of Esteco. Destructive earthquakes and surface deformation are thus a hallmark of this part of the Andean foreland.
With state-of-the-art Python packages (e.g. pyrocko, ObsPy), a semi-automatic approach is followed to analyze the collected continuous data of the seismological network. The resulting 1435 hypocenter locations consist of three different groups: 1.) local crustal earthquakes (nearly half of the events belong to this group), 2.) interplate activity, of regional distance in the slab of the Nazca-plate, and 3.) very deep earthquakes at about 600km depth. My major interest focused on the first event class. Those crustal events are partly aftershock events of the El Galpón earthquake and a second earthquake, in the south of the same fault. Further events can be considered as background seismicity of other faults within the study area. Strikingly, the seismogenic zone encompass the whole crust and propagates brittle deformation down, close to the Moho.
From the collected seismological data, a local seismic velocity model is estimated, using VELEST. After the execution of various stability tests, the robust minimum 1D-velocity model implies guiding values for the composition of the local, subsurface structure of the crust. Afterwards, performing a hypocenter relocation enables the assignment of individual earthquakes to aftershock clusters or extended seismotectonic structures. This allows the mapping of previously unknown seismogenic faults.
Finally, focal mechanisms are modeled for events with acurately located hypocenters, using the newly derived local velocity model. A compressive regime is attested by the majority of focal mechanisms, while the strike direction of the individual seismogenic structures is in agreement with the overall north – south orientation of the Central Andes, its mountain front, and individual mountain ranges in the southern Santa-Bárbara-System.
Millennial-scale palaeoclimate variability has been documented in various terrestrial and marine palaeoclimate proxy records throughout the Northern Hemisphere for the last glacial cycle. Its clear expression and rapid shifts between different states of climate (Greenland Interstadials and Stadials) represents a correlation tool beyond the resolution of e.g. luminescence dating, especially relevant for terrestrial deposits. Usually, comparison of terrestrial proxy datasets and the Greenland ice cores indicates a complex expression of millennial-scale climate variability as recorded in terrestrial geoarchives including loess. Loess is the most widespread terrestrial geoarchive of the Quaternary and especially widespread over Eurasia. However, loess often records a smoothed representation of millennial-scale variability without all fidelity when compared to the Greenland data, this being a relevant limiting feature in integrating loess with other palaeoclimate records. To better understand the loess proxy-response to millennial-scale climate variability, we simulate a proxy signal smoothing by natural processes through application of low-pass filters of delta O-18 data from Greenland, a high-resolution palaeoclimate reference record, alongside speleothem isotope records from the Black Sea-Mediterranean region. We show that low-pass filters represent rather simple models for better constraining the expression of millennial-scale climate variability in low sedimentation environments, and in sediments where proxy-response signals are most likely affected by natural smoothing (by e.g. bioturbation). Interestingly, smoothed datasets from Greenland and the Black Sea-Mediterranean region are most similar in the last similar to 15 ka and between similar to 50-30 ka. Between similar to 30-15 ka, roughly corresponding to the Last Glacial Maximum and the deglaciation, the records show dissimilarities, challenging the construction of robust correlative time-scales in this age range. From our analysis it becomes apparent that patterns of palaeoclimate signals in loess-palaeosol sequences often might be better explained by smoothed Greenland reference data than the original high-resolution Greenland dataset, or other reference data. This opens the possibility to better assess the temporal resolution and palaeoclimate potential of loess-palaeosol sequences in recording supra-regional climate patterns, as well as to securely integrate loess with other chronologically better-resolved palaeoclimate records.
Millennial-scale palaeoclimate variability has been documented in various terrestrial and marine palaeoclimate proxy records throughout the Northern Hemisphere for the last glacial cycle. Its clear expression and rapid shifts between different states of climate (Greenland Interstadials and Stadials) represents a correlation tool beyond the resolution of e.g. luminescence dating, especially relevant for terrestrial deposits. Usually, comparison of terrestrial proxy datasets and the Greenland ice cores indicates a complex expression of millennial-scale climate variability as recorded in terrestrial geoarchives including loess. Loess is the most widespread terrestrial geoarchive of the Quaternary and especially widespread over Eurasia. However, loess often records a smoothed representation of millennial-scale variability without all fidelity when compared to the Greenland data, this being a relevant limiting feature in integrating loess with other palaeoclimate records. To better understand the loess proxy-response to millennial-scale climate variability, we simulate a proxy signal smoothing by natural processes through application of low-pass filters of delta O-18 data from Greenland, a high-resolution palaeoclimate reference record, alongside speleothem isotope records from the Black Sea-Mediterranean region. We show that low-pass filters represent rather simple models for better constraining the expression of millennial-scale climate variability in low sedimentation environments, and in sediments where proxy-response signals are most likely affected by natural smoothing (by e.g. bioturbation). Interestingly, smoothed datasets from Greenland and the Black Sea-Mediterranean region are most similar in the last similar to 15 ka and between similar to 50-30 ka. Between similar to 30-15 ka, roughly corresponding to the Last Glacial Maximum and the deglaciation, the records show dissimilarities, challenging the construction of robust correlative time-scales in this age range. From our analysis it becomes apparent that patterns of palaeoclimate signals in loess-palaeosol sequences often might be better explained by smoothed Greenland reference data than the original high-resolution Greenland dataset, or other reference data. This opens the possibility to better assess the temporal resolution and palaeoclimate potential of loess-palaeosol sequences in recording supra-regional climate patterns, as well as to securely integrate loess with other chronologically better-resolved palaeoclimate records.
According to Dooge (1986) intermediate-scale catchments are systems of organized complexity, being too organized and yet too small to be characterized on a statistical/conceptual basis, but too large and too heterogeneous to be characterized in a deterministic manner. A key requirement for building structurally adequate models precisely for this intermediate scale is a better understanding of how different forms of spatial organization affect storage and release of water and energy. Here, we propose that a combination of the concept of hydrological response units (HRUs) and thermodynamics offers several helpful and partly novel perspectives for gaining this improved understanding. Our key idea is to define functional similarity based on similarity of the terrestrial controls of gradients and resistance terms controlling the land surface energy balance, rainfall runoff transformation, and groundwater storage and release. This might imply that functional similarity with respect to these specific forms of water release emerges at different scales, namely the small field scale, the hillslope, and the catchment scale. We thus propose three different types of "functional units" - specialized HRUs, so to speak - which behave similarly with respect to one specific form of water release and with a characteristic extent equal to one of those three scale levels. We furthermore discuss an experimental strategy based on exemplary learning and replicate experiments to identify and delineate these functional units, and as a promising strategy for characterizing the interplay and organization of water and energy fluxes across scales. We believe the thermodynamic perspective to be well suited to unmask equifinality as inherent in the equations governing water, momentum, and energy fluxes: this is because several combinations of gradients and resistance terms yield the same mass or energy flux and the terrestrial controls of gradients and resistance terms are largely independent. We propose that structurally adequate models at this scale should consequently disentangle driving gradients and resistance terms, because this optionally allow sequifinality to be partly reduced by including available observations, e. g., on driving gradients. Most importantly, the thermodynamic perspective yields an energy-centered perspective on rainfall-runoff transformation and evapotranspiration, including fundamental limits for energy fluxes associated with these processes. This might additionally reduce equifinality and opens up opportunities for testing thermodynamic optimality principles within independent predictions of rainfall-runoff or land surface energy exchange. This is pivotal to finding out whether or not spatial organization in catchments is in accordance with a fundamental organizing principle.
This study presents an application of an innovative sampling strategy to assess soil moisture dynamics in a headwater of the Weisseritz in the German eastern Ore Mountains. A grassland site and a forested site were instrumented with two Spatial TDR clusters (STDR) that consist of 39 and 32 coated TDR probes of 60 cm length. Distributed time series of vertically averaged soil moisture data from both sites/ensembles were analyzed by statistical and geostatistical methods. Spatial variability and the spatial mean at the forested site were larger than at the grassland site. Furthermore, clustering of TDR probes in combination with long-term monitoring allowed identification of average spatial covariance structures at the small field scale for different wetness states. The correlation length of soil water content as well as the sill to nugget ratio at the grassland site increased with increasing average wetness and but, in contrast, were constant at the forested site. As soil properties at both the forested and grassland sites are extremely variable, this suggests that the correlation structure at the forested site is dominated by the pattern of throughfall and interception. We also found a very strong correlation between antecedent soil moisture at the forested site and runoff coefficients of rainfall-runoff events observed at gauge Rehefeld. Antecedent soil moisture at the forest site explains 92% of the variability in the runoff coefficients. By combining these results with a recession analysis we derived a first conceptual model of the dominant runoff mechanisms operating in this catchment. Finally, we employed a physically based hydrological model to shed light on the controls of soil- and plant morphological parameters on soil average soil moisture at the forested site and the grassland site, respectively. A homogeneous soil setup allowed, after fine tuning of plant morphological parameters, most of the time unbiased predictions of the observed average soil conditions observed at both field sites. We conclude that the proposed sampling strategy of clustering TDR probes is suitable to assess unbiased average soil moisture dynamics in critical functional units, in this case the forested site, which is a much better predictor for event scale runoff formation than pre-event discharge. Long term monitoring of such critical landscape elements could maybe yield valuable information for flood warning in headwaters. We thus think that STDR provides a good intersect of the advantages of permanent sampling and spatially highly resolved soil moisture sampling using mobile rods.
Integration of digital elevation models and satellite images to investigate geological processes.
(2006)
In order to better understand the geological boundary conditions for ongoing or past surface processes geologists face two important questions: 1) How can we gain additional knowledge about geological processes by analyzing digital elevation models (DEM) and satellite images and 2) Do these efforts present a viable approach for more efficient research. Here, we will present case studies at a variety of scales and levels of resolution to illustrate how we can substantially complement and enhance classical geological approaches with remote sensing techniques. Commonly, satellite and DEM based studies are being used in a first step of assessing areas of geologic interest. While in the past the analysis of satellite imagery (e.g. Landsat TM) and aerial photographs was carried out to characterize the regional geologic characteristics, particularly structure and lithology, geologists have increasingly ventured into a process-oriented approach. This entails assessing structures and geomorphic features with a concept that includes active tectonics or tectonic activity on time scales relevant to humans. In addition, these efforts involve analyzing and quantifying the processes acting at the surface by integrating different remote sensing and topographic data (e.g. SRTM-DEM, SSM/I, GPS, Landsat 7 ETM, Aster, Ikonos…). A combined structural and geomorphic study in the hyperarid Atacama desert demonstrates the use of satellite and digital elevation data for assessing geological structures formed by long-term (millions of years) feedback mechanisms between erosion and crustal bending (Zeilinger et al., 2005). The medium-term change of landscapes during hundred thousands to millions years in a more humid setting is shown in an example from southern Chile. Based on an analysis of rivers/watersheds combined with landscapes parameterization by using digital elevation models, the geomorphic evolution and change in drainage pattern in the coastal Cordillera can be quantified and put into the context of seismotectonic segmentation of a tectonically active region. This has far-reaching implications for earthquake rupture scenarios and hazard mitigation (K. Rehak, see poster on IMAF Workshop). Two examples illustrate short-term processes on decadal, centennial and millennial time scales: One study uses orogen scale precipitation gradients derived from remotely sensed passive microwave data (Bookhagen et al., 2005a). They demonstrate how debris flows were triggered as a response of slopes to abnormally strong rainfall in the interior parts of the Himalaya during intensified monsoons. The area of the orogen that receives high amounts of precipitation during intensified monsoons also constitutes numerous landslide deposits of up to 1km<sup>3 volume that were generated during intensified monsoon phase at about 27 and 9 ka (Bookhagen et al., 2005b). Another project in the Swiss Alps compared sets of aerial photographs recorded in different years. By calculating high resolution surfaces the mass transport in a landslide could be reconstructed (M. Schwab, Universität Bern). All these examples, although representing only a short and limited selection of projects using remote sense data in geology, have as a common approach the goal to quantify geological processes. With increasing data resolution and new sensors future projects will even enable us to recognize more patterns and / or structures indicative of geological processes in tectonically active areas. This is crucial for the analysis of natural hazards like earthquakes, tsunamis and landslides, as well as those hazards that are related to climatic variability. The integration of remotely sensed data at different spatial and temporal scales with field observations becomes increasingly important. Many of presently highly populated places and increasingly utilized regions are subject to significant environmental pressure and often constitute areas of concentrated economic value. Combined remote sensing and ground-truthing in these regions is particularly important as geologic, seismicity and hydrologic data may be limited here due to the recency of infrastructural development. Monitoring ongoing processes and evaluating the remotely sensed data in terms of recurrence of events will greatly enhance our ability to assess and mitigate natural hazards. <hr> Dokument 1: Foliensatz | Dokument 2: Abstract <hr> Interdisziplinäres Zentrum für Musterdynamik und Angewandte Fernerkundung Workshop vom 9. - 10. Februar 2006
Surface lake sediments, 28 from Hoh Xil, 24 from northeastern China, 99 from Lake Bosten, 31 from Ulungur and 26 from Heihai were collected to determine C-13 and O-18 values. Considering the impact factors, conductivity, alkalinity, pH, TOC, C/N and carbonate-content in the sediments, Cl, P, S, and metal element ratios of Mg/Ca, Sr/Ca, Fe/Mn of bulk sediments as environmental variables enable evaluation of their influences on C-13 and O-18 using principal component analysis (PCA) method. The closure and residence time of lakes can influence the correlation between C-13 and O-18. Lake water will change from fresh to brackish with increasing reduction and eutrophication effects. Mg/Ca in the bulk sediment indicates the characteristic of residence time, Sr/Ca and Fe/Mn infer the salinity of lakes. Carbonate formation processes and types can influence the C-13-O-18 correlation. O-18 will be heavier from Mg-calcite and aragonite formed in a high-salinity water body than calcite formed in freshwater conditions. When carbonate content is less than 30%, there is no relationship with either C-13 or O-18, and also none between C-13 and O-18. More than 30%, carbonate content, however, co-varies highly to C-13 and O-18, and there is also a high correlation between C-13 and O-18. Vegetation conditions and primary productivity of lakes can influence the characteristics of C-13 and O-18, and their co-variance. Total organic matter content (TOC) in the sediments is higher with more terrestrial and submerged plants infilling. In northeastern and northwestern China, when organic matter in the lake sediments comes from endogenous floating organisms and algae, the C-13 value is high. C-13 is in the range of -4%o to 0 parts per thousand when organic matter comes mainly from floating organisms (C/N<6); in the range of -4 parts per thousand to 8 parts per thousand when organic matter comes from diatoms (C/N=6 to 8); and -8 parts per thousand to -4 parts per thousand when organic matter comes from aquatic and terrestrial plants (C/N>8).
A multi-proxy study including analyses of delta C-13(org) for the lake sediment core GN-02 and grain size, TOC. CaCO3 content, delta C-13(carb) and delta O-18(carb) of bulk carbonate, and the mineralogy of the parallel core GN-04 from Gun Nuur was performed to reconstruct the Holocene hydrology and climate on the northern Mongolian Plateau. The chronology was established using 40 C-14 dates of bulk organic matter in addition to nine previously published radiocarbon dates for core GN-02, and further five C-14 dates for the new core GN-04. A lake reservoir effect of 1060 C-14 years was determined as the intercept of the high-resolution GN-02 age-depth model at the modern sediment surface. The size of the reservoir effect is supported by the age of the core-top sample (1200 +/- 40 C-14 years) and the determined difference between a wood-derived radiocarbon age from the GN-02 core base and the age-model inferred age for bulk organic matter at the same stratigraphic level (1000 C-14 years). Low lake level and prevailing aeolian sediment deposition at Gun Nuur under dry conditions were recorded during the earliest Holocene (> 10,800-10,300 cal a BP). Gun Nuur expanded under significantly wetter conditions between 10,300 and 7000 cal a BP. Unstable climate conditions existed in the mid Holocene (7000-2500 cal a BP) and three periods of low lake-levels and significantly drier conditions were recorded between 7000-5700, 4100-3600 and 3000-2500 cal a BP. Intermediate lake levels were inferred for the intervening periods. Around 2500 cal a BP, the climate change and wetter conditions were established again. As a consequence, the lake level of Gun Nuur rose again due to higher effective moisture and the relatively wet present conditions were achieved ca. 1600 cal a BP. Our results suggest that the initial Holocene climate change on the northern Mongolian Plateau was not accompanied by a rapid increase in precipitation as on the Tibetan Plateau. The establishment of wetter conditions in northern Mongolia lagged behind the early Holocene moisture increase on the Tibetan Plateau by ca. 1000 years. Subsiding dry air in the north of the Tibetan Plateau resulted from the strengthened summer monsoon on the Tibetan Plateau during the period of maximum summer insolation and probably inhibited a significant precipitation increase in Mongolia. The significant moisture increase in the Gun Nuur region at ca. 10.3 cal ka BP is probably not related to the northward shift of the present summer monsoon boundary or the moisture delivery from the northern Atlantic through the westerlies. Instead, water from melting snow, ice and frozen ground and the generation of precipitation from the local recycling of moisture are discussed as possible moisture source for the early onset of wetter conditions on the Mongolian Plateau.
The origin of the First Bend of the Yangtze River is key to understanding the birth of the modern Yangtze River. Despite considerable efforts, the timing and mechanism of formation of the First Bend remain highly debated. Inverse river-profile modeling of three tributaries (Chongjiang, Lima, and Gudu) of the Jinsha River, integrated with regional tectonic and geomorphic interpretations, allows the onset of incision at the First Bend to be constrained to 28-20 Ma. The spatio-temporal coincidence of initial river incision and activity of Yulong Thrust Belt in southeastern Tibet highlights thrusting to be fundamental in reshaping the pre-existing stream network at the First Bend. These results enable us to reinterpret a change in sedimentary environment from a braided river to a swamp-like lake in the Jianchuan Basin south of the First Bend, recording the destruction of the hypothesized southwards-flowing paleo-Jinsha and Shuiluo Rivers at ~36-35 Ma by magmatism. During the late Oligoceneearly Miocene, the paleo-Shuiluo River was diverted to the north by focused rock uplift due to thrusting along the Yulong Thrust Belt, which also led to exhumation of the Jianchuan Basin. Diversion of the paleo-Shuiluo River can be explained by capture from a downstream river in the footwall of the Yulong Thrust Belt. Subsequent rapid headward erosion, that was caused by thrusting-induced drop of local base level, is recorded by upstream younging ages for the onset of incision and led to the formation of the First Bend. The combination of new ages for the onset of incision at 28-20 Ma at the First Bend and younger ages upstream indicates northwards expansion of the Jinsha River at a rate of 62 +/- 18 mm/yr. Our results suggest that the origin of the First Bend was likely triggered by thrusting at 28-20 Ma, after which the Yangtze River formed.
The Frasnian-Famennian (F-F) transition of Late Devonian was a critical episode in geological history, recording a major mass extinction event. In this study, we focus on an F-F succession from a deep marine context in Bancheng, southern Guangxi, South China, to investigate coeval changes in pelagic environments of the Paleo-Tethys Ocean. The studied succession is exclusively composed of bedded cherts intercalated with multiple siliceous volcanic ash beds. A SIMS zircon U-Pb Concordia age of 367.8 +/- 2.5 Ma is reported for a tuffaceous layer slightly above the F-F boundary. Geochemical ratios of Al/(Al + Fe + Mn), Ce/Ce*, Y/Ho, and Al, Fe contents in bedded cherts indicate that they are of predominantly biogenic/chemical origin with some terrigenous inputs. Negligible enrichment of redox sensitive elements (Mo, U, V) and low V/Cr ratios (<2) suggest persistently oxic conditions existed in the deep pelagic basin at Bancheng, South China during the F-F transition. These findings call into question the widely held hypothesis that marine anoxia was the primary killing mechanism for the F-F crisis. In contrast, multiple tuffaceous layers throughout the F-F boundary succession indicate frequent volcanic activity, which could have released massive amounts of greenhouse gases into the atmosphere, inducing climate warming. This scenario may have increased continental weathering and riverine fluxes into the ocean, reconciling the increases in Al2O3 content and Al/(Al + Fe + Mn) ratio across the F-F boundary. Documentation of persistently oxic conditions and frequent volcanic activitiy provides new perspectives on the inter-relationship between volcanism, climate, and oceanic redox fluctuation during the F-F biotic crisis.
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.
Pollen influx (number of pollen grains cm−2 year−1) can objectively reflect the dispersal and deposition features of pollen within a certain time and space, and is often used as a basis for the quantitative reconstruction of palaeovegetation; however, little is known about the features and mechanisms of vertical dispersal of pollen. Here we present the results from a 5 year (2006–2010) monitoring program using pollen traps placed at different heights from ground level up to 60 m and surface soil samples in a mixed coniferous and deciduous broad-leaved woodland in the Changbai mountains, northeastern China. The pollen percentages and pollen influx from the traps have very similar characteristics to the highest values for Betula, Fraxinus, Quercus and Pinus, among the tree taxa and Artemisia, Chenopodiaceae and Asteraceae among the herb taxa. Pollen influx values vary significantly with height and show major differences between three distinct layers, above-canopy (≥32 m), within the trunk layer (8 ≤ 32 m) and on the ground (0 m). These differences in pollen influx are explained by differences in (i) the air flows in each of these layers and (ii) the fall speed of pollen of the various taxa. We found that the pollen recorded on the ground surface is a good representation of the major part of the pollen transported in the trunk space of the woodland. Comparison of the pollen influx values with the theoretical, calculated “characteristic pollen source area” (CPSA) of 12 selected taxa indicates that the pollen deposited on the ground surface of the woodland is a fair representation with 85–90 % of the total pollen deposited at a wind speed of 2.4 m s−1 coming from within ca. 1–5 km for Pinus and Quercus, ca. 5–10 km for Ulmus, Tilia, Oleaceae and Betula, ca. 20–40 km for Fraxinus, Poaceae, Chenopodiaceae, Populus and Salix, and ca. 30–60 km for Artemisia; it is also a good representation with 90–98 % of the total pollen deposited coming from within 60 km at a wind speed of 2.4 m s−1, or 100 km at a wind speed: 6 m s−1, for the 12 selected taxa used in the CPSA calculation. Furthermore, comparison with the vegetation map of the area around the sampling site shows that the pollen deposited on the ground represents all plant communities which grow in the study area within 70 km radius of the sampling site. In this study, the pollen percentages obtained from the soil surface samples are significantly biased towards pollen taxa with good preservation due to thick and robust pollen walls. Therefore, if mosses are available instead, soil samples should be avoided for pollen studies, in particular for the study of pollen-vegetation relationships, the estimation of pollen productivities and quantitative reconstruction of past vegetation. The results also indicate that the existing model of pollen dispersal and deposition, Prentice’s model, provides a fair description of the actual pollen dispersal and deposition in this kind of woodland, which suggests that the application of the landscape reconstruction algorithm would be relevant for reconstruction of this type of woodland in the past.
Ostracod distribution and habitat relationships in the Kunlun Mountains, northern Tibetan Plateau
(2013)
Surface sediment samples were collected from the lakes Heihai, Kusai, Haiding Nuur and Yan Hu, and from streams and ponds in the Kunlun Mountains at the northern margin of the Tibetan Plateau to investigate the sub-fossil ostracod (micro-crustacean) fauna of the region. Among 65 collected samples, 46 ostracod shell-rich samples were used to study the relationship between the ostracod distribution and specific conductivity (SC) of the water, which ranged from 0.6 to 53.0 mS cm(-1). A total of eleven ostracod species was identified from this region, with about half of the species restricted to the Tibetan Plateau and its adjacent mountain areas, and the other half representing Holarctic taxa. Tonnacypris cf. estonica and Tonnacypris tonnensis are reported from the Tibetan Plateau for the first time. Leucocythere sp. is the dominant species and Ilyocypris cf. bradyi is also relatively abundant. The other seven species were recorded with limited abundances apparently due to lower SC tolerances. Leucocythere sp. was recorded over the full SC range from 0.6 to 53 mS cm(-1). Eucypris mareotica is a typical brackish and saline water species, which was found at sample sites with high SC (2.8-53.0 mS cm(-1)). In contrast, Leucocythere dorsotuberosa, Candona candida and Eucypris afghanistanensis prefer freshwater to slightly oli-gohaline waters with SC < 1.8 mS cm(-1). The SC optimum and tolerance range for each species were determined and compared to earlier reported data from other regions of Central Asia. The results indicate that species assemblage data from fossil ostracod shells have a large potential to provide information on past SC levels and more general climate-determined moisture conditions. (C) 2013 Elsevier Ltd and INQUA. All rights reserved.
Increasingly available high-frequency data during storm events, when hydrological dynamics most likely activate nitrate storage-flux exchanges, reveal insights into catchment nitrate dynamics. In this study, we explored impacts of seasonality and landscape gradients on nitrate concentration-discharge (C-Q) hysteresis patterns in the Selke catchment, central Germany, which has heterogeneous combinations of meteorological, hydrogeological and land use conditions. Three nested gauging stations established along the main Selke River captured flow and nitrate export dynamics from the uppermost subcatchment (mixed forest and arable land), middle subcatchment (pure steep forest) and lowermost subcatchment (arable and urban land). We collected continuous high-frequency (15-min) discharge and nitrate concentration data from 2012 to 2017 and analyzed the 223 events detected at all three stations. A dominant hysteresis pattern in the uppermost and middle subcatchments was counter-clockwise and combined with an accretion effect, indicating many proximal and mobilized distal nitrate sources. However, 66% of all events at the catchment outlet experienced a dilution effect, possibly due to mechanisms that vary seasonally. During wetting/wet periods (October-March), it was combined mainly with a counter-clockwise pattern due to the dominance of event runoff volume from the uppermost and middle subcatchments. During drying/dry periods (April-September), however, it was combined mainly with a clockwise pattern due to occasional quick surface flows from lowland near-stream urban areas. In addition, the clockwise hysteresis occurred mainly from May-October during mostly drying/dry periods at all three sites, indicating little distal nitrate transport in response to the low terrestrial hydrological connectivity, especially in the lowermost dry and flat sub-catchment. This comprehensive analysis (i.e., clockwise vs. counter-clockwise, accretion vs. dilution) enables in-depth analysis of nitrate export mechanisms during certain periods under different landscape conditions. Specific combination of C-Q relationships could identify target locations for agricultural management actions that decrease nitrate output. Therefore, we strongly encourage long-term multisite and high-frequency monitoring strategies in heterogeneous nested catchment(s), which can help understand process mechanisms, generate data for physical-based water-quality modeling and provide guidance for water and agricultural management.
The Gongjue basin from the eastern Qiangtang terrane is located in the transition region where the regional structural lineation curves from east-west-oriented in Tibet to north-south-oriented in Yunnan. In this study, we sampled the red beds in the basin from the lower Gongjue to upper Ranmugou formations for the first time covering the entire stratigraphic profile. The stratigraphic ages are bracketed within 53-43Ma by new detrital zircon U-Pb ages constraining the maximum deposition age to 52.51.5Ma. Rock magnetic and petrographic studies indicate that detrital magnetite and hematite are the magnetic carriers. Positive reversals and fold tests demonstrate that the characteristic remanent magnetization has a primary origin. The Gongjue and Ranmugou formations yield mean characteristic remanent magnetization directions of D-s/I-s=31.0 degrees/21.3 degrees and D-s/I-s=15.9 degrees/22.0 degrees, respectively. The magnetic inclination of these characteristic remanent magnetizations is significantly shallowed compared to the expected inclination for the locality. However, the elongation/inclination correction method does not provide a meaningful correction, likely because of syn-depositional rotation. Rotations relative to the Eurasian apparent polar wander path occurred in three stages: Stage I, 33.33.4 degrees clockwise rotation during the deposition of the Gongjue and lower Ranmugou formations; Stage II, 26.93.7 degrees counterclockwise rotation during deposition of the lower and middle Ranmugou formation; and Stage III, 17.73.3 degrees clockwise rotation after 43Ma. The complex rotation history recorded in the basin is possibly linked to sinistral shear along the Qiangtang block during India indentation into Asia and the early stage of the extrusion of the northwestern Indochina blocks away from eastern Tibet.
Paleogeographic reconstructions of terranes can greatly benefit from the provenance analysis of sediments. A series of Cenozoic basins provide key sedimentary archives for investigating the growth of the Tibetan Plateau, yet the provenance of the sediments in these basins has never been constrained robustly. Here we report sedimentary petrological and detrital zircon geochronological data from the Paleocene-Eocene Nangqian-Xialaxiu and Gongjue basins. Sandstone detrital modes and zircon morphology suggest that the samples collected in these two basins were sourced from recycled orogen. Detrital zircon geochronology indicates that sediments in the Nangqian-Xialaxiu Basin are characterized by two distinct age populations at 220-280 Ma and 405-445 Ma. In contrast, three predominant age populations of 207-256 Ma, 423-445 Ma, and 1851-1868 Ma, and two subordinate age populations of similar to 50 Ma and similar to 2500 Ma, are recognized in the Gongjue Basin. Comparison with detrital zircon ages from the surrounding terranes suggests that sediments in the Nangqian-Xialaxiu Basin come from the neighboring thrust belts, whereas sediments from the Gongjue Basin are predominantly derived from the distant Songpan-Ganzi Terrane with minor contribution from the surrounding areas. A three-stage Cenozoic evolution of the eastern Tibetan Plateau is proposed. During the Paleocene, the Nangqian-Xialaxiu Basin appeared as a set of small intermontane sub-basins and received plentiful sediments from the neighboring mountain belts; during the Eocene, the Gongjue Basin kept a relatively low altitude and was a depression at the edge of a proto-Plateau; since the Oligocene, the Tibetan Plateau further uplifted and the marginal Gongjue Basin was involved in the Tibetan interior orogeny, indicating the eastward propagation of the Tibetan Plateau.
The Upper Cambrian Lower Qiulitag Group in the Tarim Basin, NW China, is overwhelmingly composed of cyclic dolomites. Based on extensive field investigations and facies analysis from four outcrop sections in the Bachu-Keping area, northwestern Tarim Basin, four main types of facies are recognized: open-marine subtidal, restricted shallow subtidal, intertidal, and supratidal facies, and these are further subdivided into ten lithofacies. In general, these facies are vertically arranged into shallowing-upward, metre-scale cycles. These cycles are commonly composed of a thin basal horizon reflecting abrupt deepening, and a thicker upper succession showing gradual shallowing upwards. Based on the vertical facies arrangements and changes across boundary surfaces, two types of cycle: peritidal and shallow subtidal cycle, are further identified. The peritidal cycles, predominating over the lower-middle Lower Qiulitag Group, commence with shallow subtidal to lower intertidal facies and are capped by inter-supratidal facies. In contrast, the shallow subtidal cycles, dominating the upper Lower Qiulitag Group, are capped by shallow-subtidal facies. Based on vertical lithofacies variations, cycle stacking patterns, and accommodation variations revealed by Fischer plots, six larger-scale third-order depositional sequences (Sq1-Sq6) are recognized. These sequences generally consist of a lower transgressive and an upper regressive systems tract. The transgressive tracts are dominated by thicker-than-average cycles, indicating an overall accommodation increase, whereas the regressive tracts are characterized by thinner-than-average peritidal cycles, indicating an overall accommodation decrease. The sequence boundaries are characterized by transitional zones of stacked thinner-than-average cycles, rather than by a single surface. These sequences can further be grouped into lower-order sequence sets: the lower and upper sequence sets. The lower sequence set, including Sq1-Sq3, is characterized by peritidal facies-dominated sequences and a progressive decrease in accommodation space, indicating a longer-term fall in sea level. In contrast, the upper sequence set (Sq4-Sq6) is characterized by subtidal facies-dominated sequences and a progressive increase in accommodation space, indicating a longer-term rise in sea level.
To provide physically based wind modelling for wind erosion research at regional scale, a 3D computational fluid dynamics (CFD) wind model was developed. The model was programmed in C language based on the Navier-Stokes equations, and it is freely available as open source. Integrated with the spatial analysis and modelling tool (SAMT), the wind model has convenient input preparation and powerful output visualization. To validate the wind model, a series of experiments was conducted in a wind tunnel. A blocking inflow experiment was designed to test the performance of the model on simulation of basic fluid processes. A round obstacle experiment was designed to check if the model could simulate the influences of the obstacle on wind field. Results show that measured and simulated wind fields have high correlations, and the wind model can simulate both the basic processes of the wind and the influences of the obstacle on the wind field. These results show the high reliability of the wind model. A digital elevation model (DEM) of an area (3800 m long and 1700 m wide) in the Xilingele grassland in Inner Mongolia (autonomous region, China) was applied to the model, and a 3D wind field has been successfully generated. The clear implementation of the model and the adequate validation by wind tunnel experiments laid a solid foundation for the prediction and assessment of wind erosion at regional scale.
In order to identify the areas in the Xilingele grassland which are sensitive to wind erosion, a computational fluid dynamics model (CFD-WEM) was used to simulate the wind fields over a region of 37 km(2) which contains different topography and land use types. Previous studies revealed the important influences of topography and land use on wind erosion in the Xilingele grassland. Topography influences wind fields at large scale, and land use influences wind fields near the ground. Two steps were designed to implement the CFD wind simulation, and they were respectively to simulate the influence of topography and surface roughness on the wind. Digital elevation model (DEM) and surface roughness length were the key inputs for the CFD simulation. The wind simulation by CFD-WEM was validated by a wind data set which was measured simultaneously at six positions in the field. Three scenarios with different wind velocities were designed based on observed dust storm events, and wind fields were simulated according to these scenarios to predict the sensitive areas to wind erosion. General assumptions that cropland is the most sensitive area to wind erosion and heavily and moderately grazed grasslands are both sensitive etc. can be refined by the modelling of CFD-WEM. Aided by the results of this study, the land use planning and protection measures against wind erosion can be more efficient. Based on the case study in the Xilingele grassland, a method of regional wind erosion assessment aided by CFD wind simulation is summarized. The essence of this method is a combination of CFD wind simulation and determination of threshold wind velocity for wind erosion. Because of the physically-based simulation and the flexibility of the method, it can be generalised to other regions.
The open source computational fluid dynamics (CFD) wind model (CFD-WEM) for wind erosion research in the Xilingele grassland in Inner Mongolia (autonomous region, China) is compared with two open source CFD models Gerris and OpenFOAM. The evaluation of these models was made according to software technology, implemented methods, handling, accuracy and calculation speed. All models were applied to the same wind tunnel data set. Results show that the simplest CFD-WEM has the highest calculation speed with acceptable accuracy, and the most powerful OpenFOAM produces the simulation with highest accuracy and the lowest calculation speed. Gerris is between CFD-WEM and OpenFOAM. It calculates faster than OpenFOAM, and it is capable to solve different CFD problems. CFD-WEM is the optimal model to be further developed for wind erosion research in Inner Mongolia grassland considering its efficiency and the uncertainties of other input data. However, for other applications using CFD technology, Gerris and OpenFOAM can be good choices. This paper shows the powerful capability of open source CFD software in wind erosion study, and advocates more involvement of open source technology in wind erosion and related ecological researches.
Our study presents the results of a moment tensor inversion of 19 microseismic events with M-L between 2.0 and 3.4, associated with the stimulation operation of an enhanced geothermal reservoir in Basel, Switzerland, in 2006. We adopt a three-step procedure to retrieve point source solution parameters based on full waveform inversion. The inversion is performed by fitting displacement amplitude spectra and displacement seismograms in the first and second step, respectively, assuming a double couple source model and thus obtaining focal solutions for all 19 events. Our results are in agreement with focal mechanisms from a previous study, which employed P wave first-motion polarities from more than 40 stations, whereas our solutions are achieved using full waveform data recorded by less than 10 surface stations. In the last step, a full moment tensor inversion is performed. The results from the moment tensor inversion show an improvement on the waveform fitting compared to the double couple models, which is verified by an F-test. We investigate the stability of the moment tensor solutions by employing different velocity models. The isotropic components of the moment tensor solutions of some events are not negligible, suggesting source volume changes due to fluid injection. Events with significant isotropic components occurred mainly during the stimulation phase and close to the injection well. On the other hand, events that occurred in the post-stimulation phase are predominantly pure shear failure and located further away from the well bore. These spatio-temporal patterns can be explained by the influence of pore pressure variations during and after the hydraulic stimulation at the geothermal site. (C) 2014 Elsevier Ltd. All rights reserved.
We use a data set of 35 surface pollen samples from lake sediments, moss polsters and top soils on the north- eastern Tibetan Plateau to explore the relationship between modern pollen assemblages and contemporary vegetation patterns. The surface pollen transect spanned four vegetation zones--alpine meadow, steppe, steppe desert and desert-- under different climatic/elevational conditions. Relative representation (R (rel)) values and Principal Components Analysis (PCA) were used to determine the relationships between modern pollen and vegetation and regional climate gradients. The results show that the main vegetation zones along the regional and elevational transects can be distinguished by their modern pollen spectra. Relative to Poaceae, a high representation of Artemisia, Nitraria and Chenopodiaceae was found, while Cyperaceae and Gentiana showed values in the middle range, and Ranunculaceae, Asteraceae, Ephedra and Fabaceae had low relative representation values. PCA results indicate a high correlation between the biogeoclimatic zones and annual precipitation and annual temperature and July temperature. The Artemisia/ Chenopodiaceae ratio and the Artemisia/Cyperaceae ratio are useful tools for qualitative and semi-quantitative palaeoenvironmental reconstruction on the north-eastern Tibetan Plateau. Surface lake sediments are found to have different palynomorph spectra from moss cushion and soil samples, reflecting the larger pollen source area in the contemporary vegetation for lakes.
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.
This article offers a reconstruction of the vegetation and climate of the south-western Siberian Baraba forest-steppe area during the last ca. 8000 years. The analysis of palynological data from the sediment core of Lake Bolshie Toroki using quantitative methods has made it possible to reconstruct changes of the dominant types of vegetation and mean July air temperatures. Coniferous forests grew in the vicinity of the lake, and mean July air temperatures were similar to present-day ones between 7.9 and 7.0 kyr BP. The warmest and driest climate occurred at 7.0-5.0 kyr BP. At that time, the region had open steppe landscapes; birch groves began to spread. A cooling trend is seen after 5.5 kyr BP, when forest-steppe began to emerge. Steppe communities started to dominate again after 1.5 kyr BP. Mean July air temperatures lower than now are reconstructed for the period of 1.9-1 kyr BP, and then the temperatures became similar to present-day ones. Comparing the archaeological data on the types of economy of the population which inhabited the Baraba forest-steppe with the data on changes in the natural environment revealed a connection between the gradual transition from hunting and fishing to livestock breeding and the development of forest-steppe landscapes with a decrease in the area covered by forests. The development of the forest-steppe as an ecotonic landscape starting around 5 kyr BP might have contributed to the coexistence of several archaeological cultures with different types of economy on the same territory. (C) 2017 Elsevier Ltd. All rights reserved.
To what extent cities can be made sustainable under the mega-trends of urbanization and climate change remains a matter of unresolved scientific debate. Our inability in answering this question lies partly in the deficient knowledge regarding pivotal humanenvironment interactions. Regarded as the most well documented anthropogenic climate modification, the urban heat island (UHI) effect – the warmth of urban areas relative to the rural hinterland – has raised great public health concerns globally. Worse still, heat waves are being observed and are projected to increase in both frequency and intensity, which further impairs the well-being of urban dwellers. Albeit with a substantial increase in the number of publications on UHI in the recent decades, the diverse urban-rural definitions applied in previous studies have remarkably hampered the general comparability of results achieved. In addition, few studies have attempted to synergize the land use data and thermal remote sensing to systematically assess UHI and its contributing factors.
Given these research gaps, this work presents a general framework to systematically quantify the UHI effect based on an automated algorithm, whereby cities are defined as clusters of maximum spatial continuity on the basis of land use data, with their rural hinterland being defined analogously. By combining land use data with spatially explicit surface skin temperatures from satellites, the surface UHI intensity can be calculated in a consistent and robust manner. This facilitates monitoring, benchmarking, and categorizing UHI intensities for cities across scales. In light of this innovation, the relationship between city size and UHI intensity has been investigated, as well as the contributions of urban form indicators to the UHI intensity.
This work delivers manifold contributions to the understanding of the UHI, which have complemented and advanced a number of previous studies. Firstly, a log-linear relationship between surface UHI intensity and city size has been confirmed among the 5,000 European cities. The relationship can be extended to a log-logistic one, when taking a wider range of small-sized cities into account. Secondly, this work reveals a complex interplay between UHI intensity and urban form. City size is found to have the strongest influence on the UHI intensity, followed by the fractality and the anisometry. However, their relative contributions to the surface UHI intensity depict a pronounced regional heterogeneity, indicating the importance of considering spatial patterns of UHI while implementing UHI adaptation measures.
Lastly, this work presents a novel seasonality of the UHI intensity for individual clusters in the form of hysteresis-like curves, implying a phase shift between the time series of UHI intensity and background temperatures. Combining satellite observation and urban boundary layer simulation, the seasonal variations of UHI are assessed from both screen and skin levels. Taking London as an example, this work ascribes the discrepancies between the seasonality observed at different levels mainly to the peculiarities of surface skin temperatures associated with the incoming solar radiation. In addition, the efforts in classifying cities according to their UHI characteristics highlight the important role of regional climates in determining the UHI.
This work serves as one of the first studies conducted to systematically and statistically scrutinize the UHI. The outcomes of this work are of particular relevance for the overall spatial planning and regulation at meso- and macro levels in order to harness the benefits of rapid urbanization, while proactively minimizing its ensuing thermal stress.
This paper assesses the seasonality of the urban heat island (UHI) effect in the Greater London area (United Kingdom). Combining satellite-based observations and urban boundary layer climate modeling with the UrbClim model, the authors are able to address the seasonality of UHI intensity, on the basis of both land surface temperature (LST) and 2-m air temperature, for four individual times of the day (0130, 1030, 1330, and 2230 local time) and the daily means derived from them. An objective of this paper is to investigate whether the UHI intensities that are based on both quantities exhibit a similar hysteresis-like trajectory that is observed for LST when plotting the UHI intensity against the background temperature. The results show that the UrbClim model can satisfactorily reproduce both the observed urban rural LSTs and 2-m air temperatures as well as their differences and the hysteresis in the surface UHI. The hysteresis-like seasonality is largely absent in both the observed and modeled 2-m air temperatures, however. A sensitivity simulation of the UHI intensity to incoming solar radiation suggests that the hysteresis of the LST can mainly be attributed to the seasonal variation in incoming solar radiation.
We perform a systematic study of all cities in Europe to assess the Urban Heat Island (UHI) intensity by means of remotely sensed land surface temperature data. Defining cities as spatial clusters of urban land cover, we investigate the relationships of the UHI intensity, with the cluster size and the temperature of the surroundings. Our results show that in Europe, the UHI intensity in summer has a strong correlation with the cluster size, which can be well fitted by an empirical sigmoid model. Furthermore, we find a novel seasonality of the UHI intensity for individual clusters in the form of hysteresis-like curves. We characterize the shape and identify apparent regional patterns.
Urban climate is determined by a variety of factors, whose knowledge can help to attenuate heat stress in the context of ongoing urbanization and climate change. We study the influence of city size and urban form on the Urban Heat Island (UHI) phenomenon in Europe and find a complex interplay between UHI intensity and city size, fractality, and anisometry. Due to correlations among these urban factors, interactions in the multi-linear regression need to be taken into account. We find that among the largest 5,000 cities, the UHI intensity increases with the logarithm of the city size and with the fractal dimension, but decreases with the logarithm of the anisometry. Typically, the size has the strongest influence, followed by the compactness, and the smallest is the influence of the degree to which the cities stretch. Accordingly, from the point of view of UHI alleviation, small, disperse, and stretched cities are preferable. However, such recommendations need to be balanced against e.g. positive agglomeration effects of large cities. Therefore, trade-offs must be made regarding local and global aims.
The Indus Molasse records orogenic sedimentation associated with uplift and erosion of the southern margin of Asia in the course of ongoing India-Eurasia collision. Detailed field investigation clarifies the nature and extent of the depositional contact between this molasse and the underlying basement units. We report the first dataset on detrital zircon U-Pb ages, Hf isotopes and apatite U-Pb ages for the autochthonous molasse in the Indus Suture Zone. A latest Oligocene depositional age is proposed on the basis of the youngest detrital zircon U-Pb age peak and is consistent with published biostratigraphic data. Multiple provenance indicators suggest exclusively northerly derivation with no input from India in the lowermost parts of the section. The results provide constraints on the uplift and erosion history of the Ladakh Range following the initial India-Asia collision.
The Puna Plateau, adjacent Eastern Cordillera and the Sierras Pampeanas of the central Andes are largely characterized by thick-skinned, basement-involved deformation. The Puna Plateau hosts similar to N-S trending bedrock ranges bounded by deep-seated reverse faults and sedimentary basins. We contribute to the understanding of thick-skinned dynamics in the Puna Plateau by constraining regional kinematics of the poorly understood southern Puna Plateau through a multidisciplinary approach. On the southeastern plateau, sandstone modal composition and detrital zircon U-Pb and apatite fission-track data from Cenozoic strata indicate basin accumulation during the late Eocene to early Oligocene (similar to 38-28 Ma). Provenance analysis reveals the existence of a regional-scale basin covering the southern Puna Plateau during late Eocene to early Oligocene time (similar to 38-28 Ma) that was sourced from both the western plateau and the eastern plateau margin and had a depocenter located to the west. Petrographic and detrital zircon U-Pb data reveal erosion of proximal western and eastern sources after 12 Ma, in mid-late Miocene time. This indicates that the regional basin was compartmentalized into small-scale depocenters by the growth of basement-cored ranges continuing into the late Miocene (similar to 12-8 Ma). We suggest that the Cenozoic history of the southern Puna Plateau records the formation of a regional basin that was possibly driven by lithospheric flexure during the late Eocene to early Oligocene, before the growth of distributed basement-cored ranges starting as early as the late Oligocene. (C) 2015 Elsevier B.V. All rights reserved.
Barite concretions and bands are widely distributed in black shale-chert horizons in the Yurtus Formation of Lower Cambrian in Aksu area, northwestern Tarim Basin, NW China. They mainly consist of coarse-grained anhedral to euhedral barite crystals with minor dolomites and pyrites. Petrological features indicate these concretions grew from the porewater in unconsolidated sediments at shallow burial below sediment-water interface. The slight deviation of Sr-87/Sr-86 ratios (0.7083 to 0.7090) and significant elevated delta S-34 values (56.8-76.4 parts per thousand CDT) of barite samples with respect to those of the Early Cambrian seawater further support that barite deposits precipitated from the enclosed porewater in sediment column, which evolved from the penecontemporaneous seawater with weak interaction with the host fine-grained siliciclastic sediments and highly-depleted sulfate in response to prolonged strong bacterial sulfate reduction without necessary renewal. The abundant organic matters in the basal Yurtus Formation should have facilitated developing sulfate-depleted methanogenesis zone and sulfate-methane transition zone (SMTZ) slightly after deposition. Therefore, barite deposits in the Yurtus Formation most likely resulted from diagenetic barium cycling and persistently grew from the porewater in the static SMTZ with a low sedimentation rate in the Early Cambrian. In comparison with the distribution of sedimentary barites in geological records, we tentatively proposed that a transition in diagenetic barium cycling and associated mineralization may have occurred from the Precambrian to Cambrian periods; this scenario may be causally linked to the changes in marine ecology (the advent of mesozooplankton and associated faecal pellet) and geochemistry (the increase of seawater sulfate concentration). Thus, the occurrence of diagenetic barite deposits in the Yurtus Formation implies that diagenetic barium cycling and more effective scavenging of barium from CH4- and Ba-rich porewaters within sediments might have become an nonnegligible process in continental margin areas, at least, since the earliest Cambrian, which could have significantly impacted the marine barium cycling. (C) 2015 Elsevier B.V. All rights reserved.
The within-site variability in site response is the randomness in site response at a given site from different earthquakes and is treated as aleatory variability in current seismic hazard/risk analyses.
In this study, we investigate the single-station variability in linear site response at K-NET and KiK-net stations in Japan using a large number of earthquake recordings.
We found that the standard deviation of the horizontal-to-vertical Fourier spectral ratio at individual sites, that is single-station horizontal-to-vertical spectral ratio (HVSR) sigma sigma(HV,s), approximates the within-site variability in site response quantified using surface-to-borehole spectral ratios (for oscillator frequencies higher than the site fundamental frequency) or empirical ground-motion models.
Based on this finding, we then utilize the single-station HVSR sigma as a convenient tool to study the site-response variability at 697 KiK-net and 1169 K-NET sites.
Our results show that at certain frequencies, stiff, rough and shallow sites, as well as small and local events tend to have a higher sigma(HV,s).
However, when being averaged over different sites, the single-station HVSR sigma, that is sigma(HV), increases gradually with decreasing frequency. In the frequency range of 0.25-25 Hz, sigma(HV) is centred at 0.23-0.43 in ln scales (a linear scale factor of 1.26-1.54) with one standard deviation of less than 0.1. sigma(HV) is quite stable across different tectonic regions, and we present a constant, as well as earthquake magnitude- and distance-dependent sigma(HV) models.
In 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.
In this investigation, we examine the uncertainties using the horizontal-to-vertical spectral ratio (HVSR) technique on earthquake recordings to detect site resonant frequencies at 207 KiK-net sites. Our results show that the scenario dependence of response (pseudospectral acceleration) spectral ratio could bias the estimates of resonant frequencies for sites having multiple significant peaks with comparable amplitudes. Thus, the Fourier amplitude spectrum (FAS) should be preferred in computing HVSR. For more than 80% of the investigated sites, the first peak (in the frequency domain) on the average HVSR curve over multiple sites coincides with the highest peak. However, for sites with multiple peaks, the highest peak frequency (f(p)) is less susceptible to the selection criteria of significant peaks and the extent of smoothing to spectrum than the first peak frequency (f(0)). Meanwhile, in comparison to the surface-to-borehole spectral ratio, f(0) tends to underestimate the predominant frequency (at which the largest amplification occurs) more than f(p). In addition, in terms of characterizing linear site response, f(p) shows a better overall performance than f(0). Based on these findings, we thus recommend that seismic network operators provide f(p) on the average HVSRFAS curve as a priority, ideally together with the average HVSRFAS curve in site characterization.
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.
Evaluation of a novel application of earthquake HVSR in site-specific amplification estimation
(2020)
Ground response analyses (GRA) model the vertical propagations of SH waves through flat-layered media (1DSH) and are widely carried out to evaluate local site effects in practice. Horizontal-to-vertical spectral ratio (HVSR) technique is a cost-effective approach to extract certain site-specific information, e.g., site fundamental frequency (f(0)), but HVSR values cannot be directly used to approximate the levels of S-wave amplifications. Motivated by the work of Kawase et al. (2019), we propose a procedure to correct earthquake HVSR amplitudes for direct amplification estimations. The empirical correction compensates HVSR by generic vertical amplification spectra categorized by the vertical fundamental frequency (f(0v)) via kappa-means clustering. In this investigation, we evaluate the effectiveness of the corrected HVSR in approximating observed linear amplifications in comparison with 1DSH modellings. We select a total of 90 KiK-net (Kiban Kyoshin network) surface-downhole sites which are found to have no velocity contrasts below their boreholes and thus of which surface-to-borehole spectral ratios (SBSRs) can be taken as their empirical transfer functions (ETFs). 1DSH-based theoretical transfer functions (TTFs) are computed in the linear domain considering uncertainties in Vs profiles through randomizations. Five goodness-of-fit metrics are adopted to gauge the closeness between observed (ETF) and predicted (i.e., TTF and corrected HVSR) amplifications in both amplitude and spectral shape over frequencies from f(0) to 25 Hz. We find that the empirical correction to HVSR is highly effective and achieves a "good match" in both spectral shape and amplitude at the majority of the 90 KiK-net sites, as opposed to less than one-third for the 1DSH modelling. In addition, the empirical correction does not require a velocity model, which GRAs require, and thus has great potentials in seismic hazard assessments.
In this paper we report a rare and fortunate event of fast magnetosonic (MS, also called equatorial noise) waves modulated by compressional ultralow frequency (ULF) waves measured by Van Allen Probes. The characteristics of MS waves, ULF waves, proton distribution, and their potential correlations are analyzed. The results show that ULF waves can modulate the energetic ring proton distribution and in turn modulate the MS generation. Furthermore, the variation of MS intensities is attributed to not only ULF wave activities but also the variation of background parameters, for example, number density. The results confirm the opinion that MS waves are generated by proton ring distribution and propose a new modulation phenomenon.
Whistler mode exohiss are the structureless hiss waves observed outside the plasma pause with featured equatorward Poynting flux. An event of the amplification of exohiss as well as chorus waves was recorded by Van Allen Probes during the recovery phase of a weak geomagnetic storm. Amplitudes of both types of the waves showed a significant increase at the regions of electron density enhancements. It is found that the electrons resonant with exohiss and chorus showed moderate pitch angle anisotropies. The ratio of the number of electrons resonating with exohiss to total electron number presented in-phase correlation with density variations, which suggests that exohiss can be amplified due to electron density enhancement in terms of cyclotron instability. The calculation of linear growth rates further supports above conclusion. We suggest that exohiss waves have potential to become more significant due to the background plasma fluctuation.
Mechanical behaviors of granite after thermal treatment under loading and unloading conditions
(2021)
Understanding the mechanical behaviors of granite after thermal treatment under loading and unloading conditions is of utmost relevance to deep geothermal energy recovery. In the present study, a series of loading and unloading triaxial compression tests (20, 40 and 60 MPa) on granite specimens after exposure to different temperatures (20, 200, 300, 400, 500 and 600 degrees C) was carried out to quantify the combined effects of thermal treatment and loading/unloading stress conditions on granite strength and deformation. Changes in the microstructure of granite exposed to high temperatures were revealed by optical microscopy. The experimental results indicate that both, thermal treatment and loading/unloading stress conditions, degrade the mechanical behaviors and further decrease the carrying capacity of granite. The gradual degradation of the mechanical characteristics of granite after thermal treatment is mainly associated with the evolution of thermal micro-cracks based on optical microscopy observations. The unloading stress state induces the extension of tension cracks parallel to the axial direction, and thus, the mechanical properties are degraded. Temperatures above 400 degrees C have a more significant influence on the mechanical characteristics of granite than the unloading treatment, whereby 400 degrees C can be treated as a threshold temperature for the delineation of significant deterioration. This study is expected to support feasibility and risk assessments by means of providing data for analytical calculations and numerical simulations on granite exposed to high temperatures during geothermal energy extraction.
Understanding the geologic evolution of Northern Tibetan Plateau with multiple thermochronometers
(2018)
The early onset of deformation following the India-Asia collision, Neogene expanse of uplift, and complex systems that comprise strike-slip faults, thrust faults, and intermontane basins characterize the Cenozoic tectonism of Northern Tibetan Plateau and raise two prominent questions in orogenic geodynamics: 1) What mechanism(s) control(s) the transfer of stress related to the India-Asia collision across the distance of >2000 km; and 2) Why the development of high topography was delayed in the Northern Tibetan Plateau and what does it reveal about how the internal forces and external boundary conditions evolved. To address these two questions, we reconstruct a holistic spatial-temporal deformation history of the Northern Tibetan Plateau by using a range of thermochronometers, with closure temperature spanning from 350 degrees C to-60-70 degrees C. This multi-thermochronometer study reveals three stages of faulting related cooling, in the early Cretaceous, in Paleocene-Eocene and in middle-late Miocene. We observe that Paleocene-Eocene deformation was spatially restricted and mostly occurred on reactivated Cretaceous structures, indicating a control of pre-existing weakness on early Cenozoic deformation. Extensive Neogene deformation contrasts with restricted Paleocene-Eocene deformation and relatively quiescent shortening during the Oligocene-early Miocene, which implies a change in the regional tectonics regime. Global plate reconstructions show that this tectonic reorganization is coeval with an increase in Pacific-Asia plate convergence rates. We argue that this change in regional tectonics is a result of increasing constrictive environment of the eastern plate boundary, which changed the behavior of the Altyn Tagh fault the boundary fault of Northern Tibetan Plateau, causing it to change from feeding slip into structures out of the plateau to feeding slip into structures at plateau margins.
QuestionHow important is the effect of micro-relief and vegetation type on the characteristics of vascular plants and bryophytes in low-centred polygons? LocationSiberian Arctic, Russia. MethodsEight low-centred polygons in northern Siberia were surveyed for vegetation along transects running from the rim to the pond via the rim-pond transition of each polygon and across a vegetation type gradient from open forest to tundra. ResultsThe cover of vascular plants and bryophytes displays no significant differences between the rim and rim-pond transition but is significantly lower in the pond section of the polygons. Alpha-diversity of vascular plants decreases strongly from rim to pond, whereas bryophyte diversity in pond plots is significantly distinct from the rim and the rim-pond transition. There is no clear trend in cover for either plant group along the vegetation type transect and only a weak trend in -diversity. However, both gradients are reflected in the compositional turnover. The applied indicator species analysis identified taxa characteristic of certain environmental conditions. Among others, we found vascular plants primarily characteristic of the rim and bryophyte taxa characteristic of each micro-relief level and vegetation type. ConclusionsThe observed gradual pattern in -diversity and composition of polygonal vegetation suggests that micro-relief is the main driver of changes in the vegetation composition, while vegetation type and the related forest cover change are of subordinate importance for polygonal vegetation patterns along the Siberian tree line.
The spatial and temporal variability of a low-centred polygon on the eastern floodplain area of the lower Anabar River (72.070 degrees N, 113.921 degrees E; northern Yakutia, Siberia) has been investigated using a multi-method approach. The present-day vegetation in each square metre was analysed, revealing a community of Larix, shrubby Betula, and Salix on the polygon rim, a dominance of Carex and Andromeda polifolia in the rim-to-pond transition zone, and a predominantly monospecific Scorpidium scorpioides coverage within the pond. The total organic carbon (TOC) content, TOC/TN (total nitrogen) ratio, grain size, vascular plant macrofossils, moss remains, diatoms, and pollen were analysed for two vertical sections and a sediment core from a transect across the polygon. Radiocarbon dating indicates that the formation of the polygon started at least 1500 yr ago; the general positions of the pond and rim have not changed since that time. Two types of pond vegetation were identified, indicating two contrasting development stages of the polygon. The first was a well-established moss association, dominated by submerged or floating Scorpidium scorpioides and/or Drepanocladus spp. and overgrown by epiphytic diatoms such as Tabellaria flocculosa and Eunotia taxa. This stage coincides temporally with a period in which the polygon was only drained by lateral subsurface water flow, as indicated by mixed grain sizes. A different moss association occurred during times of repeated river flooding (indicated by homogeneous medium-grained sand that probably accumulated during the annual spring snowmelt), characterized by an abundance of Meesia triquetra and a dominance of benthic diatoms (e. g. Navicula vulpina), indicative of a relatively high pH and a high tolerance of disturbance. A comparison of the local polygon vegetation (inferred from moss and macrofossil spectra) with the regional vegetation (inferred from pollen spectra) indicated that the moss association with Scorpidium scorpioides became established during relatively favourable climatic conditions, while the association dominated by Meesia triquetra occurred during periods of harsh climatic conditions. Our study revealed a strong riverine influence (in addition to climatic influences) on polygon development and the type of peat accumulated.
Applying conservation of energy to estimate earthquake frequencies from strain rates and stresses
(2020)
Estimating earthquake occurrence rates from the accumulation rate of seismic moment is an established tool of seismic hazard analysis. We propose an alternative, fault-agnostic approach based on the conservation of energy: the Energy-Conserving Seismicity Framework (ENCOS). Working in energy space has the advantage that the radiated energy is a better predictor of the damage potential of earthquake waves than the seismic moment release. In a region, ENCOS balances the stationary power available to cause earthquakes with the long-term seismic energy release represented by the energy-frequency distribution's first moment. Accumulation and release are connected through the average seismic efficiency, by which we mean the fraction of released energy that is converted into seismic waves. Besides measuring earthquakes in energy, ENCOS differs from moment balance essentially in that the energy accumulation rate depends on the total stress in addition to the strain rate tensor. To validate ENCOS, we exemplarily model the energy-frequency distribution around Southern California. We estimate the energy accumulation rate due to tectonic loading assuming poroelasticity and hydrostasis. Using data from the World Stress Map and assuming the frictional limit to estimate the stress tensor, we obtain a power of 0.8 GW. The uncertainty range, 0.3-2.0GW, originates mainly from the thickness of the seismogenic crust, the friction coefficient on preexisting faults, and models of Global Positioning System (GPS) derived strain rates. Based on a Gutenberg-Richter magnitude-frequency distribution, this power can be distributed over a range of energies consistent with historical earthquake rates and reasonable bounds on the seismic efficiency.
Information on the contemporary in-situ stress state of the earth’s crust is essential for geotechnical applications and physics-based seismic hazard assessment. Yet, stress data records for a data point are incomplete and their availability is usually not dense enough to allow conclusive statements. This demands a thorough examination of the in-situ stress field which is achieved by 3D geomechanicalnumerical models. However, the models spatial resolution is limited and the resulting local stress state is subject to large uncertainties that confine the significance of the findings. In addition, temporal variations of the in-situ stress field are naturally or anthropogenically induced. In my thesis I address these challenges in three manuscripts that investigate (1) the current crustal stress field orientation, (2) the 3D geomechanical-numerical modelling of the in-situ stress state, and (3) the phenomenon of injection induced temporal stress tensor rotations. In the first manuscript I present the first comprehensive stress data compilation of Iceland with 495 data records. Therefore, I analysed image logs from 57 boreholes in Iceland for indicators of the orientation of the maximum horizontal stress component. The study is the first stress survey from different kinds of stress indicators in a geologically very young and tectonically active area of an onshore spreading ridge. It reveals a distinct stress field with a depth independent stress orientation even very close to the spreading centre. In the second manuscript I present a calibrated 3D geomechanical-numerical modelling approach of the in-situ stress state of the Bavarian Molasse Basin that investigates the regional (70x70x10km³) and local (10x10x10km³) stress state. To link these two models I develop a multi-stage modelling approach that provides a reliable and efficient method to derive from the larger scale model initial and boundary conditions for the smaller scale model. Furthermore, I quantify the uncertainties in the models results which are inherent to geomechanical-numerical modelling in general and the multi-stage approach in particular. I show that the significance of the models results is mainly reduced due to the uncertainties in the material properties and the low number of available stress magnitude data records for calibration. In the third manuscript I investigate the phenomenon of injection induced temporal stress tensor rotation and its controlling factors. I conduct a sensitivity study with a 3D generic thermo-hydro-mechanical model. I show that the key control factors for the stress tensor rotation are the permeability as the decisive factor, the injection rate, and the initial differential stress. In particular for enhanced geothermal systems with a low permeability large rotations of the stress tensor are indicated. According to these findings the estimation of the initial differential stress in a reservoir is possible provided the permeability is known and the angle of stress rotation is observed. I propose that the stress tensor rotations can be a key factor in terms of the potential for induced seismicity on pre-existing faults due to the reorientation of the stress field that changes the optimal orientation of faults.
Rotations of the principal stress axes are observed as a result of fluid injection into reservoirs. We use a generic, fully coupled 3-D thermo-hydro-mechanical model to investigate systematically the dependence of this stress rotation on different reservoir properties and injection scenarios. We find that permeability, injection rate, and initial differential stress are the key factors, while other reservoir properties only play a negligible role. In particular, we find that thermal effects do not significantly contribute to stress rotations. For reservoir types with usual differential stress and reservoir treatment the occurrence of significant stress rotations is limited to reservoirs with a permeability of less than approximately 10(-12)m(2). Higher permeability effectively prevents stress rotations to occur. Thus, according to these general findings, the observed principal stress axes rotation can be used as a proxy of the initial differential stress provided that rock permeability and fluid injection rate are known a priori.
Iceland is located on the Mid-Atlantic Ridge which is the plate boundary between the Eurasian and the North American plates. It is one of the few places on earth where an active spreading centre is located onshore but the stress pattern has not been extensively investigated so far. In this paper we present a comprehensive compilation of the orientation of maximum horizontal stress (S-Hmax). In particular we interpret borehole breakouts and drilling induced fractures from borehole image logs in 57 geothermal wells onshore Iceland. The borehole results are combined with other stress indicators including earthquake focal mechanism solutions, geological information and overcoring measurements resulting in a dataset with 495 data records for the S-Hmax orientation. The reliability of each indicator is assessed according to the quality criteria of the World Stress Map project The majority of S-Hmax orientation data records in Iceland is derived from earthquake focal mechanism solutions (35%) and geological fault slip inversions (26%). 20% of the data are borehole related stress indicators. In addition minor shares of S-Hmax orientations are compiled, amongst others, from focal mechanism inversions and the alignment of fissure eruptions. The results show that the S-Hmax orientations derived from different depths and stress indicators are consistent with each other.
The resulting pattern of the present-day stress in Iceland has four distinct subsets of S-Hmax orientations. The S-Hmax orientation is parallel to the rift axes in the vicinity of the active spreading regions. It changes from NE-SW in the South to approximately N-S in central Iceland and NNW-SSE in the North. In the Westfjords which is located far away from the ridge the regional S-Hmax rotates and is parallel to the plate motion. (C) 2016 Elsevier B.V. All rights reserved.
On 27 December 2007, a 1.9 seismic event occurred within a dyke in the deep-level Mponeng Gold Mine, South Africa. From the seismological network of the mine and the one from the Japanese-German Underground Acoustic Emission Research in South Africa (JAGUARS) group, the hypocentral depth (3,509 m), focal mechanism and aftershock location were estimated. Since no mining activity took place in the days before the event, dynamic triggering due to blasting can be ruled out as the cause. To investigate the hypothesis that stress transfer, due to excavation of the gold reef, induced the event, we set up a small-scale high-resolution three-dimensional (3D) geomechanical numerical model. The model consisted of the four different rock units present in the mine: quartzite (footwall), hard lava (hanging wall), conglomerate (gold reef) and diorite (dykes). The numerical solution was computed using a finite-element method with a discretised mesh of approximately elements. The initial stress state of the model is in agreement with in situ data from a neighbouring mine, and the step-wise excavation was simulated by mass removal from the gold reef. The resulting 3D stress tensor and its changes due to mining were analysed based on the Coulomb failure stress changes on the fault plane of the event. The results show that the seismic event was induced regardless of how the Coulomb failure stress changes were calculated and of the uncertainties in the fault plane solution. We also used the model to assess the seismic hazard due to the excavation towards the dyke. The resulting curve of stress changes shows a significant increase in the last in front of the dyke, indicating that small changes in the mining progress towards the dyke have a substantial impact on the stress transfer.
The seismicity of the Kenya rift is characterized by high-frequency low-magnitude events concentrated along the rift axis. Its seismic character is typical for magmatically active continental rifts, where igneous material at a shallow depth causes extensive grid faulting and geothermal activity. Thermal overprinting and dike intrusion prohibit the buildup of large elastic strains, therefore prohibiting the generation of large-magnitude earthquakes. On 6 January 1928, the M-S 6.9 Subukia earthquake occurred on the Laikipia-Marmanet fault, the eastern rift-bounding structure of the central Kenya rift. It is the largest instrumentally recorded seismic event in the Kenya rift, standing in contrast to the current model of the rift's seismic character in which large earthquakes are not anticipated. Furthermore, the proximity of the ruptured fault and the rift axis is intriguing: The rift-bounding structure that ruptured in 1928 remains seismically active, capable of generating large-magnitude earthquakes, even though thermally weakened crust and better oriented structures are present along the rift axis nearby, prohibiting any significant buildup of elastic strain. We excavated the surface rupture of the 1928 Subukia earthquake to find evidence for preceding ground-rupturing earthquakes. We also made a total station survey of the site topography and mapped the site geology. We show that the Laikipia-Marmanet fault was repeatedly activated during the late Quaternary. We found evidence for six ground-rupturing earthquakes, including the 1928 earthquake. The topographic survey around the trench site revealed a degraded fault scarp of approximate to 7.5 m in height, offsetting a small debris slide. Using scarp-diffusion modeling, we estimated an uplift rate of U = 0.09-0.15 mm/yr, constraining the scarp age to 50-85 ka. Assuming an average fault dip of 55 degrees-75 degrees, the preferred uplift rate (0.15 mm/yr) accommodates approximately 10%-20% of the recent rate of extension (0.5 mm/yr) across the Kenya rift.
Questions of identity and provenance of minerals that are parts of masterpieces in museums have become increasingly important in mineralogical and historical studies. Detailed investigations of valuable and unique objects require on-site, nondestructive and noninvasive methods because touching or removing them may cause irreparable damage. A mobile Raman-microprobe has been used to meet these demands for truly in situ mineralogical studies of the large collection of minerals and rocks of the Prussian kings in the Grotto Hall (Grottensaal) of the New Palace (Neues Palais), Park Sanssouci in Potsdam. Minerals on the walls of the Grotto Hall were analyzed to identify them and thereby to complete the data bank of the collection. Fluid and solid inclusions in the interior of a large quartz crystal have been studied to provide evidence of the provenance of the crystal. The fluid inclusions contain aqueous saline solutions, whereas the solid inclusions are needles of anhydrite with a length of about 1.5 mm. The quartz probably originated from an area in the eastern Alps, from the surroundings of Bad Gastein, Austria. This is the first on-site and in situ study of inclusions below the surface of a mineral with a mobile Raman-microprobe outside a laboratory.
Thermobarometrical and mineral-chemical investigations by electron microprobe and LA-ICP-MS on a sillimanite- bearing pegmatoid from the Reinbolt Hills provide important constraints on the P-T-X-age relations of part of East Antarctica during Pan-African tectonism. U-Th-total Pb ages of monazite imply that the pegmatoid of originally Grenvillan age (zircon U-Pb age of ca. 900 Ma) underwent a major, late Pan-African (Cambrian) regional, granulite-facies metamorphism between 500 and 550 Ma. Most of the monazite formed during this event, as result of apatite metasomatism owing to infiltration of high-grade metamorphic fluids. Apatite-biotite and other mineral thermobarometers define the peak metamorphic temperatures and pressures with 850-950 degrees C and 0.8-1.0 GPa. The F-Cl-OH relations in apatite, and biotite, the chemistry of fluid inclusions and the presence of K-feldspar microveins suggest that the metasomatising fluid was a CO2-bearing, diluted KCl brine. The pegmatoid is the first record of monazite-(Ce) formed from fluorapatite that is rich in U (up to 2.6 Wt% UO2) and possesses Th/U ratios <1 (0.09 on average). These chemical signatures are direct reflection of the U and Th concentration patterns in the parental fluorapatite
The 10th edition of the International Congress on the Application of Raman Spectroscopy in Art and Archaeology (RAA2019) was held in Potsdam (Germany) from 3 to 7 September 2019, with eight keynote lectures, 35 oral presentations and 18 Poster Presentations. The number of active participants was 68 delegates from 20 countries among the 236 authors that presented at least one work.
Forests seem to represent low-erosion systems, according to most, but not all, studies of suspended-sediment yield. We surmised that this impression reflects an accidental bias in the selection of monitoring sites towards those with prevailing vertical hydrological flowpaths, rather than a tight causal link between vegetation cover and erosion alone. To evaluate this conjecture, we monitored, over a 2-year period, a 3.3 ha old-growth rainforest catchment prone to frequent and widespread overland flow. We sampled stream flow at two and overland flow at three sites in a nested arrangement on a within-event basis, and monitored the spatial and temporal frequency of overland flow. Suspended-sediment concentrations were modeled with Random Forest and Quantile Regression Forest to be able to estimate the annual yields for the 2 years, which amounted to 1 t ha(-1) and 2 t ha(-1) in a year with below-average and with average precipitation, respectively. These estimates place our monitoring site near the high end of reported suspended-sediment yields and lend credence to the notion that low yields reflect primarily the dominance of vertical flowpaths and not necessarily and exclusively the kind of vegetative cover. Undisturbed forest and surface erosion are certainly no contradiction in terms even in the absence of mass movements.
Saturated hydraulic conductivity (K-s) is an important soil characteristic affecting soil water storage, runoff generation and erosion processes. In some areas where high-intensity rainfall coincides with low K-s values at shallow soil depths, frequent overland flow entails dense drainage networks. Consequently, linear structures such as flowlines alternate with inter-flowline areas. So far, investigations of the spatial variability of K-s mainly relied on isotropic covariance models which are unsuitable to reveal patterns resulting from linear structures. In the present study, we applied two sampling approaches so as to adequately characterize K-s spatial variability in a tropical forest catchment that features a high density of flowlines: A classical nested sampling survey and a purposive sampling strategy adapted to the presence of flowlines. The nested sampling approach revealed the dominance of small-scale variability, which is in line with previous findings. Our purposive sampling, however, detected a strong spatial gradient: surface K-s increased substantially as a function of distance to flowline; 10 m off flowlines, values were similar to the spatial mean of K-s. This deterministic trend can be included as a fixed effect in a linear mixed modeling framework to obtain realistic spatial fields of K-s. In a next step we used probability maps based on those fields and prevailing rainfall intensities to assess the hydrological relevance of the detected pattern. This approach suggests a particularly good agreement between the probability statements of K-s exceedance and observed overland flow occurrence during wet stages of the rainy season.
Numerous studies investigated the influence of abiotic (meteorological conditions) and biotic factors (tree characteristics) on stemflow generation. Although these studies identified the variables that influence stemflow volumes in simply structured forests, the combination of tree characteristics that allows a robust prediction of stemflow volumes in species-rich forests is not well known. Many hydrological applications, however, require at least a rough estimate of stemflow volumes based on the characteristics of a forest stand. The need for robust predictions of stemflow motivated us to investigate the relationships between tree characteristics and stemflow volumes in a species-rich tropical forest located in central Panama. Based on a sampling setup consisting of ten rainfall collectors, 300 throughfall samplers and 60 stemflow collectors and cumulated data comprising 26 rain events, we derive three main findings. Firstly, stemflow represents a minor hydrological component in the studied 1-ha forest patch (1.0% of cumulated rainfall). Secondly, in the studied species-rich forest, single tree characteristics are only weakly related to stemflow volumes. The influence of multiple tree parameters (e.g. crown diameter, presence of large epiphytes and inclination of branches) and the dependencies among these parameters require a multivariate approach to understand the generation of stemflow. Thirdly, predicting stemflow in species-rich forests based on tree parameters is a difficult task. Although our best model can capture the variation in stemflow to some degree, a critical validation reveals that the model cannot provide robust predictions of stemflow. A reanalysis of data from previous studies in species-rich forests corroborates this finding. Based on these results and considering that for most hydrological applications, stemflow is only one parameter among others to estimate, we advocate using the base model, i.e. the mean of the stemflow data, to quantify stemflow volumes for a given study area. Studies in species-rich forests that wish to obtain predictions of stemflow based on tree parameters probably need to conduct a much more extensive sampling than currently implemented by most studies. Copyright (c) 2015 John Wiley & Sons, Ltd.
estimating mean throughfall
(2016)
The selection of an appropriate spatial extent of a sampling plot is one among several important decisions involved in planning a throughfall sampling scheme. In fact, the choice of the extent may determine whether or not a study can adequately characterize the hydrological fluxes of the studied ecosystem. Previous attempts to optimize throughfall sampling schemes focused on the selection of an appropriate sample size, support, and sampling design, while comparatively little attention has been given to the role of the extent. In this contribution, we investigated the influence of the extent on the representativeness of mean throughfall estimates for three forest ecosystems of varying stand structure. Our study is based on virtual sampling of simulated throughfall fields. We derived these fields from throughfall data sampled in a simply structured forest (young tropical forest) and two heterogeneous forests (old tropical forest, unmanaged mixed European beech forest). We then sampled the simulated throughfall fields with three common extents and various sample sizes for a range of events and for accumulated data. Our findings suggest that the size of the study area should be carefully adapted to the complexity of the system under study and to the required temporal resolution of the throughfall data (i.e. event-based versus accumulated). Generally, event-based sampling in complex structured forests (conditions that favor comparatively long autocorrelations in throughfall) requires the largest extents. For event-based sampling, the choice of an appropriate extent can be as important as using an adequate sample size. (C) 2016 Elsevier B.V. All rights reserved.
A wide range of basic and applied problems in water resources research requires high-quality estimates of the spatial mean of throughfall. Many throughfall sampling schemes, however, are not optimally adapted to the system under study. The application of inappropriate sampling schemes may partly reflect the lack of generally applicable guidelines on throughfall sampling strategies. In this study we conducted virtual sampling experiments using simulated fields which are based on empirical throughfall data from three structurally distinct forests (a 12-year old teak plantation, a 5-year old young secondary forest, and a 130-year old secondary forest). In the virtual sampling experiments we assessed the relative error of mean throughfall estimates for 38 different throughfall sampling schemes comprising a variety of funnel- and trough-type collectors and a large range of sample sizes. Moreover, we tested the performance of each scheme for both event-based and accumulated throughfall data. The key findings of our study are threefold. First, as errors of mean throughfall estimates vary as a function of throughfall depth, the decision on which temporal scale (i.e. event-based versus accumulated data) to sample strongly influences the required sampling effort. Second, given a chosen temporal scale throughfall estimates can vary considerably as a function of canopy complexity. Accordingly, throughfall sampling in simply structured forests requires a comparatively modest effort, whereas heterogeneous forests can be extreme in terms of sampling requirements, particularly if the focus is on reliable data of small events. Third, the efficiency of trough-type collectors depends on the spatial structure of throughfall. Strong, long-ranging throughfall patterns decrease the efficiency of troughs substantially. Based on the results of our virtual sampling experiments, which we evaluated by applying two contrasting sampling approaches simultaneously, we derive readily applicable guidelines for throughfall monitoring. (C) 2014 Elsevier B.V. All rights reserved.
Changes in rainfall interception along a secondary forest succession gradient in lowland Panama
(2013)
Secondary forests are rapidly expanding in tropical regions. Yet, despite the importance of understanding the hydrological consequences of land-cover dynamics, the relationship between forest succession and canopy interception is poorly understood. This lack of knowledge is unfortunate because rainfall interception plays an important role in regional water cycles and needs to be quantified for many modeling purposes. To help close this knowledge gap, we designed a throughfall monitoring study along a secondary succession gradient in a tropical forest region of Panama. The investigated gradient comprised 20 forest patches 3 to 130 yr old. We sampled each patch with a minimum of 20 funnel-type throughfall collectors over a continuous 2month period that had nearly 900 mm of rain. During the same period, we acquired forest inventory data and derived several forest structural attributes. We then applied simple and multiple regression models (Bayesian model averaging, BMA) and identified those vegetation parameters that had the strongest influence on the variation of canopy interception. Our analyses yielded three main findings. First, canopy interception changed rapidly during forest succession. After only a decade, throughfall volumes approached levels that are typical for mature forests. Second, a parsimonious (simple linear regression) model based on the ratio of the basal area of small stems to the total basal area outperformed more complex multivariate models (BMA approach). Third, based on complementary forest inventory data, we show that the influence of young secondary forests on interception in realworld fragmented landscapes might be detectable only in regions with a substantial fraction of young forests. Our re-sults suggest that where entire catchments undergo forest regrowth, initial stages of succession may be associated with a substantial decrease of streamflow generation. Our results further highlight the need to study hydrological processes in all forest succession stages, including early ones.
In the humid tropics, continuing high deforestation rates are seen alongside an increasing expansion of secondary forests. In order to understand and model the consequences of these dynamic land-use changes for regional water cycles, the response of soil hydraulic properties to forest disturbance and recovery has to be quantified.At a site in the Brazilian Amazonia, we annually monitored soil infiltrability and saturated hydraulic conductivity (K-s) at 12.5, 20 cm, and 50 cm soil depth after manual forest conversion to pasture (year zero to four after pasture establishment), and during secondary succession after pasture abandonment (year zero to seven after pasture abandonment). We evaluated the hydrological consequences of the detected changes by comparing the soil hydraulic properties with site-specific rainfall intensities and hydrometric observations. Within one year after grazing started, infiltrability and K-s at 12.5 and 20 cm depth decreased by up to one order of magnitude to levels which are typical for 20-year-old pasture. In the three subsequent monitoring years, infiltrability and K-s remained stable. Land use did not impact on subsoil permeability. Whereas infiltrability values are large enough to allow all rainwater to infiltrate even after the conversion, the sudden decline of near-surface K-s is of hydrological relevance as perched water tables and overland flow occur more often on pastures than in forests at our study site. After pasture abandonment and during secondary succession, seven years of recovery did not suffice to significantly increase infiltrability and K-s at 12.5 depth although a slight recovery is obvious. At 20 cm soil depth, we detected a positive linear increase within the seven-year time frame but annual means did not differ significantly. Although more than a doubling of infiltrability and K-s is still required to achieve pre-disturbance levels, which will presumably take more than a decade, the observed slight increases of K-s might already decrease the probability of perched water table generation and overland flow development well before complete recovery.
What is the most appropriate sampling scheme to estimate event-based average throughfall? A satisfactory answer to this seemingly simple question has yet to be found, a failure which we attribute to previous efforts' dependence on empirical studies. Here we try to answer this question by simulating stochastic throughfall fields based on parameters for statistical models of large monitoring data sets. We subsequently sampled these fields with different sampling designs and variable sample supports. We evaluated the performance of a particular sampling scheme with respect to the uncertainty of possible estimated means of throughfall volumes. Even for a relative error limit of 20%, an impractically large number of small, funnel-type collectors would be required to estimate mean throughfall, particularly for small events. While stratification of the target area is not superior to simple random sampling, cluster random sampling involves the risk of being less efficient. A larger sample support, e.g., the use of trough-type collectors, considerably reduces the necessary sample sizes and eliminates the sensitivity of the mean to outliers. Since the gain in time associated with the manual handling of troughs versus funnels depends on the local precipitation regime, the employment of automatically recording clusters of long troughs emerges as the most promising sampling scheme. Even so, a relative error of less than 5% appears out of reach for throughfall under heterogeneous canopies. We therefore suspect a considerable uncertainty of input parameters for interception models derived from measured throughfall, in particular, for those requiring data of small throughfall events.
Soils in various places of the Panama Canal Watershed feature a low saturated hydraulic conductivity (K-s) at shallow depth, which promotes overland-flow generation and associated flashy catchment responses. In undisturbed forests of these areas, overland flow is concentrated in flow lines that extend the channel network and provide hydrological connectivity between hillslopes and streams. To understand the dynamics of overland-flow connectivity, as well as the impact of connectivity on catchment response, we studied an undisturbed headwater catchment by monitoring overland-flow occurrence in all flow lines and discharge, suspended sediment, and total phosphorus at the catchment outlet. We find that connectivity is strongly influenced by seasonal variation in antecedent wetness and can develop even under light rainfall conditions. Connectivity increased rapidly as rainfall frequency increased, eventually leading to full connectivity and surficial drainage of entire hillslopes. Connectivity was nonlinearly related to catchment response. However, additional information on factors such as overland-flow volume would be required to constrain relationships between connectivity, stormflow, and the export of suspended sediment and phosphorus. The effort to monitor those factors would be substantial, so we advocate applying the established links between rain event characteristics, drainage network expansion by flow lines, and catchment response for predictive modeling and catchment classification in forests of the Panama Canal Watershed and in similar regions elsewhere.
Organic matter deposited in ancient, ice-rich permafrost sediments is vulnerable to climate change and may contribute to the future release of greenhouse gases; it is thus important to get a better characterization of the plant organic matter within such sediments. From a Late Quaternary permafrost sediment core from the Buor Khaya Peninsula, we analysed plant-derived sedimentary ancient DNA (sedaDNA) to identify the taxonomic composition of plant organic matter, and undertook palynological analysis to assess the environmental conditions during deposition. Using sedaDNA, we identified 154 taxa and from pollen and non-pollen palynomorphs we identified 83 taxa. In the deposits dated between 54 and 51 kyr BP, sedaDNA records a diverse low-centred polygon plant community including recurring aquatic pond vegetation while from the pollen record we infer terrestrial open-land vegetation with relatively dry environmental conditions at a regional scale. A fluctuating dominance of either terrestrial or swamp and aquatic taxa in both proxies allowed the local hydrological development of the polygon to be traced. In deposits dated between 11.4 and 9.7 kyr BP (13.4-11.1 cal kyr BP), sedaDNA shows a taxonomic turnover to moist shrub tundra and a lower taxonomic richness compared to the older samples. Pollen also records a shrub tundra community, mostly seen as changes in relative proportions of the most dominant taxa, while a decrease in taxonomic richness was less pronounced compared to sedaDNA. Our results show the advantages of using sedaDNA in combination with palynological analyses when macrofossils are rarely preserved. The high resolution of the sedaDNA record provides a detailed picture of the taxonomic composition of plant-derived organic matter throughout the core, and palynological analyses prove valuable by allowing for inferences of regional environmental conditions.
The Fram Strait is an area with a relatively low and irregular distribution of diatom microfossils in surface sediments, and thus microfossil records are scarce, rarely exceed the Holocene, and contain sparse information about past richness and taxonomic composition. These attributes make the Fram Strait an ideal study site to test the utility of sedimentary ancient DNA (sedaDNA) metabarcoding. Amplifying a short, partial rbcL marker from samples of sediment core MSM05/5-712-2 resulted in 95.7% of our sequences being assigned to diatoms across 18 different families, with 38.6% of them being resolved to species and 25.8% to genus level. Independent replicates show a high similarity of PCR products, especially in the oldest samples. Diatom sedaDNA richness is highest in the Late Weichselian and lowest in Mid- and Late Holocene samples. Taxonomic composition is dominated by cold-water and sea-ice-associated diatoms and suggests several reorganisations - after the Last Glacial Maximum, after the Younger Dryas, and after the Early and after the Mid-Holocene. Different sequences assigned to, amongst others, Chaetoceros socialis indicate the detectability of intra-specific diversity using sedaDNA. We detect no clear pattern between our diatom sedaDNA record and the previously published IP25 record of this core, although proportions of pennate diatoms increase with higher IP25 concentrations and proportions of Nitzschia cf. frigida exceeding 2% of the assemblage point towards past sea-ice presence.
The study presents the first evidence of metamorphosed submarine ashes in the Edough Massif, in northeastern Algeria. It occurs below the greenschist-facies Tellian units that represent the thrusted Mesozoic to Eocene passive paleomargin of northern Africa deposited on thinned continental crust. The metamorphic complex consists of tectonically superposed units composed of gneisses (lower unit) and micaschists (upper unit). At the Cap de Garde, these units enclose an "intermediate unit" composed of micaschists and meter-thick layers of marbles, which are sometimes intercalated with amphibolites. The latter occur as discontinuous small lenses and layers. The amphibolites are parallel to the primary bedding of the marbles and the main foliation. Chemical markers and field observations indicate that they are metamorphic equivalents of basic igneous rocks. The lenticular character, low thickness and multiple intercalations with marine sediments and the unusual high lithium concentrations suggest subaqueous near-source basaltic ash-fall deposits in a marine environment. (C) 2014 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.
Determination of the source rocks for the diatexites from the Edough Massif, Annaba, NE Algeria
(2012)
The crystalline Edough Massif is located in the oriental part of the Algerian coastline. It consists of two tectonically superposed units of gneisses, augen-gneisses and migmatitic gneisses in the lower unit and micaschists in the upper unit. The crystalline rocks underwent a low to moderate degree of metamorphism; the gneisses suffered partial melting. They display migmatitic features such as nebulitic structures with contorted leucosome layers and K-feldspar porphyroblasts and thus can be classified as diatexites. The mineralogical composition of these rocks is very homogenous and consists of K-feldspar, micas and quartz. The feldspar-rich, arkosic nature of the outcrop implies a granitic source rock. High K2O/Na2O ratios and high A/CNK > 1.1 indicate an S-type granite source and a peraluminous composition of the protolith respectively. Chondrite normalized REE distribution patterns of the Edough diatexites show gently inclined patterns with a minor negative Eu anomaly (Eu/Eu* = 0.36-0.49), which points to a very slightly differentiated granitic source. The REE pattern and trace element data of the diatexites are similar to those of average Proterozoic upper continental crust, which suggests that they are derived mainly from upper continental crust and were deposited in continental margins.
The sidi Mohamed peridotites (Edough Massif, NE Algeria) - evidence for an upper mantle origin
(2013)
The Hercynian Edough massif is the easternmost crystalline massif of the Algerian coast. It consists of two tectonically superposed units composed of micaschists, gneisses, and peridotite. This study concentrates on the small and isolated Sidi Mohamed peridotite outcrop area (0.03 km(2)). The Sidi Mohamed peridotite is composed mainly of harzburgites (Mg-rich olivine and orthopyroxene as major minerals). The Ni (2051-2920 ppm), Cr (2368-5514 ppm) and MgO (similar to 28-35 wt.%) whole-rock composition and the relative depletion in Nb make these harzburgites comparable to depleted peridotites related to a subduction zone. We suggest that the Sidi Mohamed ultramafic body was derived directly from the upper mantle and tectonically incorporated into the gneiss units of the Edough metamorphic core complex in a subduction environment.
Lava domes are severely hazardous, mound-shaped extrusions of highly viscous lava and commonly erupt at many active stratovolcanoes around the world. Due to gradual growth and flank oversteepening, such lava domes regularly experience partial or full collapses, resulting in destructive and far-reaching pyroclastic density currents. They are also associated with cyclic explosive activity as the complex interplay of cooling, degassing, and solidification of dome lavas regularly causes gas pressurizations on the dome or the underlying volcano conduit. Lava dome extrusions can last from days to decades, further highlighting the need for accurate and reliable monitoring data.
This thesis aims to improve our understanding of lava dome processes and to contribute to the monitoring and prediction of hazards posed by these domes. The recent rise and sophistication of photogrammetric techniques allows for the extraction of observational data in unprecedented detail and creates ideal tools for accomplishing this purpose. Here, I study natural lava dome extrusions as well as laboratory-based analogue models of lava dome extrusions and employ photogrammetric monitoring by Structure-from-Motion (SfM) and Particle-Image-Velocimetry (PIV) techniques. I primarily use aerial photography data obtained by helicopter, airplanes, Unoccupied Aircraft Systems (UAS) or ground-based timelapse cameras. Firstly, by combining a long time-series of overflight data at Volcán de Colima, México, with seismic and satellite radar data, I construct a detailed timeline of lava dome and crater evolution. Using numerical model, the impact of the extrusion on dome morphology and loading stress is further evaluated and an impact on the growth direction is identified, bearing important implications for the location of collapse hazards. Secondly, sequential overflight surveys at the Santiaguito lava dome, Guatemala, reveal surface motion data in high detail. I quantify the growth of the lava dome and the movement of a lava flow, showing complex motions that occur on different timescales and I provide insight into rock properties relevant for hazard assessment inferred purely by photogrammetric processing of remote sensing data. Lastly, I recreate artificial lava dome and spine growth using analogue modelling under controlled conditions, providing new insights into lava extrusion processes and structures as well as the conditions in which they form.
These findings demonstrate the capabilities of photogrammetric data analyses to successfully monitor lava dome growth and evolution while highlighting the advantages of complementary modelling methods to explain the observed phenomena. The results presented herein further bear important new insights and implications for the hazards posed by lava domes.
During eruptive activity of andesitic stratovolcanoes, the extrusion of lava domes, their collapse and intermittent explosions are common volcanic hazards. Many lava domes grow in a preferred direction, in turn affecting the direction of lava flows and pyroclastic density currents. Access to active lava domes is difficult and hazardous, so detailed data characterizing lava dome growth are typically limited, keeping the processes controlling the directionality of extrusions unclear. Here we combine TerraSAR-X satellite radar observations with high-resolution airborne photogrammetry to assess morphological changes, and perform finite element modeling to investigate the impact of loading stress on shallow magma ascent directions associated with lava dome extrusion and crater formation at Volcan de Colima, Mexico. The TerraSAR-X data, acquired in similar to 1-m resolution spotlight mode, enable us to derive a chronology of the eruptive processes from intensity-based time-lapse observations of the general crater and dome evolution. The satellite images are complemented by close-range airborne photos, processed by the Structure-from-Motion workflow. This allows the derivation of high-resolution digital elevation models, providing insight into detailed loading and unloading features. During the observation period from Jan-2013 to Feb-2016, we identify a dominantly W-directed dome growth and lava flow production until Jan-2015. In Feb-2015, following the removal of the active summit dome, the surface crater widened and elongated along a NE-SW axis. Later in May-2015, a new dome grew toward the SW of the crater while a separate vent developed in the NE of the crater, reflecting a change in the direction of magma ascent and possible conduit bifurcation. Finite element models show a significant stress change in agreement with the observed magma ascent direction changes in response to the changing surface loads, both for loading (dome growth) and unloading (crater forming excavation) cases. These results allow insight into shallow dome growth dynamics and the migration of magma ascent in response to changing volcano summit morphology. They further highlight the importance of detailed volcano summit morphology surveillance, as changes in direction or location of dome extrusion may have major implications regarding the directions of potential volcanic hazards, such as pyroclastic density currents generated by dome collapse.
During the period 750-600 Ma ago, prior to the final break-up of the supercontinent Rodinia, the crust of both the North American Craton and Baltica was intruded by significant amounts of rift-related magmas originating from the mantle. In the Proterozoic crust of Southern Norway, the 580 Ma old Fen carbonatite-ultramafic complex is a representative of this type of rocks. In this paper, we report the occurrence of an ultramafic lamprophyre dyke which possibly is linked to the Fen complex, although Ar-40/Ar-39 data from phenocrystic phlogopite from the dyke gave an age of 686 +/- 9 Ma. The lamprophyre dyke was recently discovered in one of the Kongsberg silver mines at Vinoren, Norway. Whole rock geochemistry, geochronological and mineralogical data from the ultramafic lamprophyre dyke are presented aiming to elucidate its origin and possible geodynamic setting. From the whole-rock composition of the Vinoren dyke, the rock could be recognized as transitional between carbonatite and kimberlite-II (orangeite). From its diagnostic mineralogy, the rock is classified as aillikite. The compositions and xenocrystic nature of several of the major and accessory minerals from the Vinoren aillikite are characteristic for diamondiferous rocks (kimberlites/lamproites/UML): Phlogopite with kinoshitalite-rich rims, chromite-spinel-ulvospinel series, Mg- and Mn-rich ilmenites, rutile and lucasite-(Ce). We suggest that the aillikite melt formed during partial melting of a MARID (mica-amphibole-rutile-ilmenite-diopside)-like source under CO2 fluxing. The pre-rifting geodynamic setting of the Vinoren aillikite before the Rodinia supercontinent breakup suggests a relatively thick SCLM (Subcontinental Lithospheric Mantle) during this stage and might indicate a diamond-bearing source for the parental melt. This is in contrast to the about 100 Ma younger Fen complex, which were derived from a thin SCLM.
The Lombok Island is part of the Lesser Sunda Islands (LSI) region – Indonesia, situated along the Sunda-Banda Arcs transition. It lies between zones characterized by the highest intensity geomagnetic anomalies of this region, remarkable as one of the eight most important features provided on the 1st edition of World Digital Magnetic Anomaly Map. The seismicity of this region during the last years is high, while the geological and tectonic structures of this region are still not known in detail. Some local magnetic surveys have been conducted previously during 2004–2005. However, due to the lower accuracy of the used equipment and a limited number of stations, the qualities of the previous measurements are questionable for more interpretations. Thus a more detailed study to better characterize the geomagnetic anomaly -spatially and temporally- over this region and to deeply explore the related regional geology, tectonic and seismicity is needed. The intriguing geomagnetic anomalies over this island region vis-à-vis the socio-cultural situations lead to a study with a special aim to contribute to the assessment of the potential of natural hazards (earthquakes) as well as a new natural resource of energy (geothermal potential).
This study is intended to discuss several crucial questions, including:
i. The real values and the general pattern of magnetic anomalies over the island, as well as their relation to the regional one.
ii. Any temporal changes of regional anomalies over the recent time.
iii. The relationships between the anomalies and the geology and tectonic of this region, especially new insights that can be gained from the geomagnetic observations.
iv. The relationships between the anomalies and the high seismicity of this region, especially some possible links between their variations to the earthquake occurrence.
First, all available geomagnetic data of this region and results of the previous measurements are evaluated. The new geomagnetic surveys carried out in 2006 and 2007/2008 are then presented in detail, followed by the general description of data processing and data quality evaluation. The new results show the general pattern of contiguous negative-positive anomalies, revealing an active arc related subduction region. They agree with earlier results obtained by satellite, aeromagnetic, and marine platforms; and provide a much more detailed picture of the strong anomalies on this island. The temporal characteristics of regional anomalies show a decreasing strength of the dipolar structure, where decreasing of the field intensities is faster than the regional secular variations as defined by the global model (the 10th generation of IGRF). However, some exceptions (increasing of anomalies) have to be noted and further analyzed for several locations.
Thereafter, simultaneous magnetic anomalies and gravity models are generated and interpreted in detail. Three profiles are investigated, providing new insights into the tectonics and geological evolution of the Lombok Island. Geological structure of this island can be divided as two main parts with different consecutive ages: an old part (from late Oligocene to late Miocene) in the South and a younger one (from Pliocene to Holocene) in the North. A new subduction in the back arc region (the Flores Thrust zone) is considered mature and active, showing a tendency of progressive subduction during 2005–2008. Geothermal potential in the northern part of this island can be mapped in more detail using these geomagnetic regional survey data. The earlier estimates of reservoir depth can be confirmed further to a depth of about 800 m. Evaluation of temporal changes of the anomalies gives some possible explanations related to the evolution of the back arc region, large stress accumulations over the LSI region, a specific electrical characteristic of the crust of the Lombok Island region, and a structural discontinuity over this island.
Based on the results, several possible advanced studies involving geomagnetic data and anomaly investigations over the Lombok Island region can be suggested for the future:
i. Monitoring the subduction activity of the back arc region (the Flores Thrust zone) and the accumulated stress over the LSI, that could contribute to middle term hazard assessment with a special attention to the earthquake occurrence in this region. Continuous geomagnetic field measurements from a geomagnetic observatory which can be established in the northern part of the Lombok Island and systematic measurements at several repeat stations can be useful in this regards.
ii. Investigating the specific electrical characteristic (high conductivity) of the crust, that is probably related to some aquifer layers or metal mineralization. It needs other complementary geophysical methods, such as magnetotelluric (MT) or preferably DC resistivity measurements.
iii. Determining the existence of an active structural fault over the Lombok Island, that could be related to long term hazard assessment over the LSI region. This needs an extension of geomagnetic investigations over the neighbouring islands (the Bali Island in the West and the Sumbawa Island in the East; probably also the Sumba and the Flores islands). This seems possible because the regional magnetic lineations might be used to delineate some structural discontinuities, based on the modelling of contrasts in crustal magnetizations.
Throughout the last ~3 million years, the Earth's climate system was characterised by cycles of glacial and interglacial periods. The current warm period, the Holocene, is comparably stable and stands out from this long-term cyclicality. However, since the industrial revolution, the climate has been increasingly affected by a human-induced increase in greenhouse gas concentrations. While instrumental observations are used to describe changes over the past ~200 years, indirect observations via proxy data are the main source of information beyond this instrumental era. These data are indicators of past climatic conditions, stored in palaeoclimate archives around the Earth. The proxy signal is affected by processes independent of the prevailing climatic conditions. In particular, for sedimentary archives such as marine sediments and polar ice sheets, material may be redistributed during or after the initial deposition and subsequent formation of the archive. This leads to noise in the records challenging reliable reconstructions on local or short time scales. This dissertation characterises the initial deposition of the climatic signal and quantifies the resulting archive-internal heterogeneity and its influence on the observed proxy signal to improve the representativity and interpretation of climate reconstructions from marine sediments and ice cores.
To this end, the horizontal and vertical variation in radiocarbon content of a box-core from the South China Sea is investigated. The three-dimensional resolution is used to quantify the true uncertainty in radiocarbon age estimates from planktonic foraminifera with an extensive sampling scheme, including different sample volumes and replicated measurements of batches of small and large numbers of specimen. An assessment on the variability stemming from sediment mixing by benthic organisms reveals strong internal heterogeneity. Hence, sediment mixing leads to substantial time uncertainty of proxy-based reconstructions with error terms two to five times larger than previously assumed.
A second three-dimensional analysis of the upper snowpack provides insights into the heterogeneous signal deposition and imprint in snow and firn. A new study design which combines a structure-from-motion photogrammetry approach with two-dimensional isotopic data is performed at a study site in the accumulation zone of the Greenland Ice Sheet. The photogrammetry method reveals an intermittent character of snowfall, a layer-wise snow deposition with substantial contributions by wind-driven erosion and redistribution to the final spatially variable accumulation and illustrated the evolution of stratigraphic noise at the surface. The isotopic data show the preservation of stratigraphic noise within the upper firn column, leading to a spatially variable climate signal imprint and heterogeneous layer thicknesses. Additional post-depositional modifications due to snow-air exchange are also investigated, but without a conclusive quantification of the contribution to the final isotopic signature.
Finally, this characterisation and quantification of the complex signal formation in marine sediments and polar ice contributes to a better understanding of the signal content in proxy data which is needed to assess the natural climate variability during the Holocene.
Marine sedimentary archives are routinely used to reconstruct past environmental changes. In many cases, bioturbation and sedimentary mixing affect the proxy time-series and the age-depth relationship. While idealized models of bioturbation exist, they usually assume homogeneous mixing, thus that a single sample is representative for the sediment layer it is sampled from.
However, it is largely unknown to which extent this assumption holds for sediments used for paleoclimate reconstructions.
To shed light on
1) the age-depth relationship and its full uncertainty,
2) the magnitude of mixing processes affecting the downcore proxy variations, and
3) the representativity of the discrete sample for the sediment layer, we designed and performed a case study on South China Sea sediment material which was collected using a box corer and which covers the last glacial cycle.
Using the radiocarbon content of foraminiferal tests as a tracer of time, we characterize the spatial age-heterogeneity of sediments in a three-dimensional setup. In total, 118 radiocarbon measurements were performed on defined small- and large-volume bulk samples ( similar to 200 specimens each) to investigate the horizontal heterogeneity of the sediment. Additionally, replicated measurements on small numbers of specimens (10 x 5 specimens) were performed to assess the heterogeneity within a sample volume. Visual assessment of X-ray images and a quantitative assessment of the mixing strength show typical mixing from bioturbation corresponding to around 10 cm mixing depth.
Notably, our 3D radiocarbon distribution reveals that the horizontal heterogeneity (up to 1,250 years), contributing to the age uncertainty, is several times larger than the typically assumed radiocarbon based age-model error (single errors up to 250 years). Furthermore, the assumption of a perfectly bioturbated layer with no mixing underneath is not met.
Our analysis further demonstrates that the age-heterogeneity might be a function of sample size; smaller samples might contain single features from the incomplete mixing and are thus less representative than larger samples.
We provide suggestions for future studies, optimal sampling strategies for quantitative paleoclimate reconstructions and realistic uncertainty in age models, as well as discuss possible implications for the interpretation of paleoclimate records.
Ecologists carry a well-stocked toolbox with a great variety of sampling methods, statistical analyses and modelling tools, and new methods are constantly appearing. Evaluation and optimisation of these methods is crucial to guide methodological choices. Simulating error-free data or taking high-quality data to qualify methods is common practice. Here, we emphasise the methodology of the 'virtual ecologist' (VE) approach where simulated data and observer models are used to mimic real species and how they are 'virtually' observed. This virtual data is then subjected to statistical analyses and modelling, and the results are evaluated against the 'true' simulated data. The VE approach is an intuitive and powerful evaluation framework that allows a quality assessment of sampling protocols, analyses and modelling tools. It works under controlled conditions as well as under consideration of confounding factors such as animal movement and biased observer behaviour. In this review, we promote the approach as a rigorous research tool, and demonstrate its capabilities and practical relevance. We explore past uses of VE in different ecological research fields, where it mainly has been used to test and improve sampling regimes as well as for testing and comparing models, for example species distribution models. We discuss its benefits as well as potential limitations, and provide some practical considerations for designing VE studies. Finally, research fields are identified for which the approach could be useful in the future. We conclude that VE could foster the integration of theoretical and empirical work and stimulate work that goes far beyond sampling methods, leading to new questions, theories, and better mechanistic understanding of ecological systems.
Empirical species distribution models (SDMs) constitute often the tool of choice for the assessment of rapid climate change effects on species vulnerability. Conclusions regarding extinction risks might be misleading, however, because SDMs do not explicitly incorporate dispersal or other demographic processes. Here, we supplement SDMs with a dynamic population model 1) to predict climate-induced range dynamics for black grouse in Switzerland, 2) to compare direct and indirect measures of extinction risks, and 3) to quantify uncertainty in predictions as well as the sources of that uncertainty. To this end, we linked models of habitat suitability to a spatially explicit, individual-based model. In an extensive sensitivity analysis, we quantified uncertainty in various model outputs introduced by different SDM algorithms, by different climate scenarios and by demographic model parameters. Potentially suitable habitats were predicted to shift uphill and eastwards. By the end of the 21st century, abrupt habitat losses were predicted in the western Prealps for some climate scenarios. In contrast, population size and occupied area were primarily controlled by currently negative population growth and gradually declined from the beginning of the century across all climate scenarios and SDM algorithms. However, predictions of population dynamic features were highly variable across simulations. Results indicate that inferring extinction probabilities simply from the quantity of suitable habitat may underestimate extinction risks because this may ignore important interactions between life history traits and available habitat. Also, in dynamic range predictions uncertainty in SDM algorithms and climate scenarios can become secondary to uncertainty in dynamic model components. Our study emphasises the need for principal evaluation tools like sensitivity analysis in order to assess uncertainty and robustness in dynamic range predictions. A more direct benefit of such robustness analysis is an improved mechanistic understanding of dynamic species responses to climate change.
Sub-seasonal thaw slump mass wasting is not consistently energy limited at the landscape scale
(2018)
Predicting future thaw slump activity requires a sound understanding of the atmospheric drivers and geomorphic controls on mass wasting across a range of timescales. On sub-seasonal timescales, sparse measurements indicate that mass wasting at active slumps is often limited by the energy available for melting ground ice, but other factors such as rainfall or the formation of an insulating veneer may also be relevant. To study the sub-seasonal drivers, we derive topographic changes from single-pass radar interferometric data acquired by the TanDEM-X satellites. The estimated elevation changes at 12m resolution complement the commonly observed planimetric retreat rates by providing information on volume losses. Their high vertical precision (around 30 cm), frequent observations (11 days) and large coverage (5000 km(2)) allow us to track mass wasting as drivers such as the available energy change during the summer of 2015 in two study regions. We find that thaw slumps in the Tuktoyaktuk coastlands, Canada, are not energy limited in June, as they undergo limited mass wasting (height loss of around 0 cm day 1) despite the ample available energy, suggesting the widespread presence of early season insulating snow or debris veneer. Later in summer, height losses generally increase (around 3 cm day 1), but they do so in distinct ways. For many slumps, mass wasting tracks the available energy, a temporal pattern that is also observed at coastal yedoma cliffs on the Bykovsky Peninsula, Russia. However, the other two common temporal trajectories are asynchronous with the available energy, as they track strong precipitation events or show a sudden speed-up in late August respectively. The observed temporal patterns are poorly related to slump characteristics like the headwall height. The contrasting temporal behaviour of nearby thaw slumps highlights the importance of complex local and temporally varying controls on mass wasting.
Kijko et al. (2016) present various methods to estimate parameters that are relevant for probabilistic seismic-hazard assessment. One of these parameters, although not the most influential, is the maximum possible earthquake magnitude m(max). I show that the proposed estimation of m(max) is based on an erroneous equation related to a misuse of the estimator in Cooke (1979) and leads to unstable results. So far, reported finite estimations of m(max) arise from data selection, because the estimator in Kijko et al. (2016) diverges with finite probability. This finding is independent of the assumed distribution of earthquake magnitudes. For the specific choice of the doubly truncated Gutenberg-Richter distribution, I illustrate the problems by deriving explicit equations. Finally, I conclude that point estimators are generally not a suitable approach to constrain m(max).
In this study, the spatial and temporal impacts of the Ataturk Dam on agro-meteorological aspects of the Southeastern Anatolia region have been investigated. Change detection and environmental impacts due to water-reserve changes in Ataturk Dam Lake have been determined and evaluated using multi-temporal Landsat satellite imageries and meteorological datasets within a period of 1984-2011. These time series have been evaluated for three time periods. Dam construction period constitutes the first part of the study. Land cover/use changes especially on agricultural fields under the Ataturk Dam Lake and its vicinity have been identified between the periods of 1984-1992. The second period comprises the 10-year period after the completion of filling up the reservoir in 1992. At this period, Landsat and meteorological time-series analyses are examined to assess the impact of the Ataturk Dam Lake on selected irrigated agricultural areas. For the last 9-year period from 2002 to 2011, the relationships between seasonal water-reserve changes and irrigated plains under changing climatic factors primarily driving vegetation activity (monthly, seasonal, and annual fluctuations of rainfall rate, air temperature, humidity) on the watershed have been investigated using a 30-year meteorological time series. The results showed that approximately 368 km(2) of agricultural fields have been affected because of inundation due to the Ataturk Dam Lake. However, irrigated agricultural fields have been increased by 56.3% of the total area (1552 of 2756 km(2)) on Harran Plain within the period of 1984-2011.
Landslides are frequent natural hazards in rugged terrain, when the resisting frictional force of the surface of rupture yields to the gravitational force. These forces are functions of geological and morphological factors, such as angle of internal friction, local slope gradient or curvature, which remain static over hundreds of years; whereas more dynamic triggering events, such as rainfall and earthquakes, compromise the force balance by temporarily reducing resisting forces or adding transient loads. This thesis investigates landslide distribution and orientation due to landslide triggers (e.g. rainfall) at different scales (6-4∙10^5 km^2) and aims to link rainfall movement with the landslide distribution. It additionally explores the local impacts of the extreme rainstorms on landsliding and the role of precursory stability conditions that could be induced by an earlier trigger, such as an earthquake.
Extreme rainfall is a common landslide trigger. Although several studies assessed rainfall intensity and duration to study the distribution of thus triggered landslides, only a few case studies quantified spatial rainfall patterns (i.e. orographic effect). Quantifying the regional trajectories of extreme rainfall could aid predicting landslide prone regions in Japan. To this end, I combined a non-linear correlation metric, namely event synchronization, and radial statistics to assess the general pattern of extreme rainfall tracks over distances of hundreds of kilometers using satellite based rainfall estimates. Results showed that, although the increase in rainfall intensity and duration positively correlates with landslide occurrence, the trajectories of typhoons and frontal storms were insufficient to explain landslide distribution in Japan. Extreme rainfall trajectories inclined northwestwards and were concentrated along some certain locations, such as coastlines of southern Japan, which was unnoticed in the landslide distribution of about 5000 rainfall-triggered landslides. These landslides seemed to respond to the mean annual rainfall rates.
Above mentioned findings suggest further investigation on a more local scale to better understand the mechanistic response of landscape to extreme rainfall in terms of landslides. On May 2016 intense rainfall struck southern Germany triggering high waters and landslides. The highest damage was reported at the Braunsbach, which is located on the tributary-mouth fan formed by the Orlacher Bach. Orlacher Bach is a ~3 km long creek that drains a catchment of about ~6 km^2. I visited this catchment in June 2016 and mapped 48 landslides along the creek. Such high landslide activity was not reported in the nearby catchments within ~3300 km^2, despite similar rainfall intensity and duration based on weather radar estimates. My hypothesis was that several landslides were triggered by rainfall-triggered flash floods that undercut hillslope toes along the Orlacher Bach. I found that morphometric features such as slope and curvature play an important role in landslide distribution on this micro scale study site (<10 km^2). In addition, the high number of landslides along the Orlacher Bach could also be boosted by accumulated damages on hillslopes due karst weathering over longer time scales.
Precursory damages on hillslopes could also be induced by past triggering events that effect landscape evolution, but this interaction is hard to assess independently from the latest trigger. For example, an earthquake might influence the evolution of a landscape decades long, besides its direct impacts, such as landslides that follow the earthquake. Here I studied the consequences of the 2016 Kumamoto Earthquake (MW 7.1) that triggered some 1500 landslides in an area of ~4000 km^2 in central Kyushu, Japan. Topography, i.e. local slope and curvature, both amplified and attenuated seismic waves, thus controlling the failure mechanism of those landslides (e.g. progressive). I found that topography fails in explaining the distribution and the preferred orientation of the landslides after the earthquake; instead the landslides were concentrated around the northeast of the rupture area and faced mostly normal to the rupture plane. This preferred location of the landslides was dominated mainly by the directivity effect of the strike-slip earthquake, which is the propagation of wave energy along the fault in the rupture direction; whereas amplitude variations of the seismic radiation altered the preferred orientation. I suspect that the earthquake directivity and the asymmetry of seismic radiation damaged hillslopes at those preferred locations increasing landslide susceptibility. Hence a future weak triggering event, e.g. scattered rainfall, could further trigger landslides at those damaged hillslopes.