Refine
Year of publication
Document Type
- Article (2597)
- Doctoral Thesis (407)
- Postprint (138)
- Other (74)
- Review (51)
- Preprint (17)
- Monograph/Edited Volume (13)
- Conference Proceeding (12)
- Habilitation Thesis (11)
- Master's Thesis (6)
Language
- English (3328) (remove)
Keywords
- climate change (51)
- Holocene (44)
- erosion (27)
- Himalaya (26)
- permafrost (26)
- Climate change (23)
- remote sensing (23)
- Tibetan Plateau (22)
- Earthquake source observations (21)
- Pollen (20)
Institute
- Institut für Geowissenschaften (3328) (remove)
In this study, we investigated the scale sizes of equatorial plasma irregularities (EPIs) using measurements from the Swarm satellites during its early mission and final constellation phases. We found that with longitudinal separation between Swarm satellites larger than 0.4 degrees, no significant correlation was found any more. This result suggests that EPI structures include plasma density scale sizes less than 44 km in the zonal direction. During the Swarm earlier mission phase, clearly better EPI correlations are obtained in the northern hemisphere, implying more fragmented irregularities in the southern hemisphere where the ambient magnetic field is low. The previously reported inverted-C shell structure of EPIs is generally confirmed by the Swarm observations in the northern hemisphere, but with various tilt angles. From the Swarm spacecrafts with zonal separations of about 150 km, we conclude that larger zonal scale sizes of irregularities exist during the early evening hours (around 1900 LT).
In this study, we investigated the scale sizes of equatorial plasma irregularities (EPIs) using measurements from the Swarm satellites during its early mission and final constellation phases. We found that with longitudinal separation between Swarm satellites larger than 0.4°, no significant correlation was found any more. This result suggests that EPI structures include plasma density scale sizes less than 44 km in the zonal direction. During the Swarm earlier mission phase, clearly better EPI correlations are obtained in the northern hemisphere, implying more fragmented irregularities in the southern hemisphere where the ambient magnetic field is low. The previously reported inverted-C shell structure of EPIs is generally confirmed by the Swarm observations in the northern hemisphere, but with various tilt angles. From the Swarm spacecrafts with zonal separations of about 150 km, we conclude that larger zonal scale sizes of irregularities exist during the early evening hours (around 1900 LT).
By using 3-year global positioning system (GPS) measurements from December 2013 to November 2016, we provide in this study a detailed survey on the climatology of the GPS signal loss of Swarm onboard receivers. Our results show that the GPS signal losses prefer to occur at both low latitudes between +/- 5 and +/- 20 degrees magnetic latitude (MLAT) and high latitudes above 60 degrees MLAT in both hemispheres. These events at all latitudes are observed mainly during equinoxes and December solstice months, while totally absent during June solstice months. At low latitudes the GPS signal losses are caused by the equatorial plasma irregularities shortly after sunset, and at high latitude they are also highly related to the large density gradients associated with ionospheric irregularities. Additionally, the high-latitude events are more often observed in the Southern Hemisphere, occurring mainly at the cusp region and along nightside auroral latitudes. The signal losses mainly happen for those GPS rays with elevation angles less than 20 degrees, and more commonly occur when the line of sight between GPS and Swarm satellites is aligned with the shell structure of plasma irregularities. Our results also confirm that the capability of the Swarm receiver has been improved after the bandwidth of the phase-locked loop (PLL) widened, but the updates cannot radically avoid the interruption in tracking GPS satellites caused by the ionospheric plasma irregularities. Additionally, after the PLL bandwidth increased larger than 0.5 Hz, some unexpected signal losses are observed even at middle latitudes, which are not related to the ionospheric plasma irregularities. Our results suggest that rather than 1.0 Hz, a PLL bandwidth of 0.5 Hz is a more suitable value for the Swarm receiver.
The Relative Pollen Productivities (RPPs) of common steppe species are estimated using Extended R-value (ERV) model based on pollen analysis and vegetation survey of 30 surface soil samples from typical steppe area of northern China. Artemisia, Chenopodiaceae, Poaceae, Cyperaceae, and Asteraceae are the dominant pollen types in pollen assemblages, reflecting the typical steppe communities well. The five dominant pollen types and six common types (Thalictrum, Iridaceae, Potentilla, Ephedra, Brassicaceae, and Ulmus) have strong wind transport abilities; the estimated Relevant Source Area of Pollen (RSAP) is ca. 1000 m when the sediment basin radius is set at 0.5 m. Ulmus, Artemisia, Brassicaceae, Chenopodiaceae, and Thalictrum have relative high RPPs; Poaceae, Cyperaceae, Potentilla, and Ephedra pollen have moderate RPPs; Asteraceae and Iridaceae have low RPPs. The reliability test of RPPs revealed that most of the RPPs are reliable in past vegetation reconstruction. However, the RPPs of Asteraceae and Iridaceae are obviously underestimated, and those of Poaceae, Chenopodiaceae, and Ephedra are either slightly underestimated or slightly overestimated, suggesting that those RPPs should be considered with caution. These RPPs were applied to estimating plant abundances for two fossil pollen spectra (from the Lake Bayanchagan and Lake Haoluku) covering the Holocene in typical steppe area, using the "Regional Estimates of Vegetation Abundance from Large Sites" (REVEALS) model. The RPPs-based vegetation reconstruction revealed that meadow-steppe dominated by Poaceae, Cyperaceae, and Artemisia plants flourished in this area before 6500-5600 cal yr BP, and then was replaced by present typical steppe.
Badlands have long been considered as model landscapes due to their perceived close relationship between form and process. The often intense features of erosion have also attracted many geomorphologists because the associated high rates of erosion appeared to offer the opportunity for studying surface processes and the resulting forms. Recently, the perceived simplicity of badlands has been questioned because the expected relationships between driving forces for erosion and the resulting sediment yield could not be observed. Further, a high variability in erosion and sediment yield has been observed across scales. Finally, denudation based on currently observed erosion rates would have lead to the destruction of most badlands a long time ago. While the perceived simplicity of badlands has sparked a disproportional (compared to the land surface they cover) amount of research, our increasing amount of information has not necessarily increased our understanding of badlands in equal terms. Overall, badlands appear to be more complex than initially assumed. In this paper, we review 40 years of research in the Zin Valley Badlands in Israel to reconcile some of the conflicting results observed there and develop a perspective on the function of badlands as model landscapes. While the data collected in the Zin Valley clearly confirm that spatial and temporal patterns of geomorphic processes and their interaction with topography and surface properties have to be understood, we still conclude that the process of realizing complexity in the "simple" badlands has a model function both for our understanding as well as perspective on all landscape systems.
Coupling land-use change and hydrologic models for quantification of catchment ecosystem services
(2018)
Representation of land-use and hydrologic interactions in respective models has traditionally been problematic. The use of static land-use in most hydrologic models or that of the use of simple hydrologic proxies in land-use change models call for more integrated approaches. The objective of this study is to assess whether dynamic feedback between land-use change and hydrology can (1) improve model performances, and/or (2) produce a more realistic quantification of ecosystem services. To test this, we coupled a land-use change model and a hydrologic mode. First, the land-use change and the hydrologic models were separately developed and calibrated. Then, the two models were dynamically coupled to exchange data at yearly time-steps. The approach is applied to a catchment in South Africa. Performance of coupled models when compared to the uncoupled models were marginal, but the coupled models excelled at the quantification of catchment ecosystem services more robustly.
The equatorial electrojet is an enhanced eastward current in the dayside E region ionosphere flowing along the magnetic equator. The equatorial electrojet is highly variable as it is subject to various forcing mechanisms including atmospheric waves from the lower layers of the atmosphere. There are occasionally times when the intensity of the equatorial electrojet at a fixed longitude shows an oscillatory variation with a period of approximately 6days. We present case studies of such events based on the equatorial electrojet measurements from the CHAMP and Swarm satellites. The spatial and temporal variability of the equatorial electrojet intensity during these events reveals characteristics of a westward propagating wave with zonal wavenumber 1, consistent with the effect of the quasi-6-day planetary wave. Analyses of the geopotential height data from the Aura satellite confirm the presence of the quasi-6-day planetary wave in the lower thermosphere during the events. The amplitude of the quasi-6-day variation in the equatorial electrojet intensity depends on longitude, but no systematic longitudinal dependence is found for different events. During the event of August 2010, quasi-6-day variations are also observed by ground-based magnetometers and a radar in the Peruvian sector. The effect of the quasi-6-day wave accounts for up to +/- 5.9m/s in the equatorial vertical plasma velocity at noon, which is much larger than previously predicted by a numerical model. These results suggest that the quasi-6-day planetary wave is an important source of short-term variability in the equatorial ionosphere.
Understanding the variability of the ionosphere is important for the prediction of space weather and climate. Recent studies have shown that forcing from the lower atmosphere plays a significant role for the short-term (day-to-day) variability of the low-latitude ionosphere. The present study aims to assess the importance of atmospheric forcing for the variability of the daytime equatorial ionospheric electric field on the interannual (year-to-year) time scale. Magnetic field measurements from Huancayo (12.05 degrees S, 75.33 degrees W) are used to augment the equatorial vertical plasma drift velocity (V-Z) measurements from the Jicamarca Unattended Long-term Investigations of the Ionosphere and Atmosphere radar during 2001-2016. V-Z can be regarded as a measure of the zonal electric field. After removing the seasonal variation of similar to 10m/s, midday values of V-Z show an interannual variation of similar to 2m/s with an oscillation period of 2-3years. No evidence of solar cycle influence is found. The Ground-to-topside Atmosphere-Ionosphere model for Aeronomy, which takes into account realistic atmospheric variability below 30km, reproduces the pattern of the observed interannual variation without having to include variable forcing from the magnetosphere. The results indicate that lower atmospheric forcing plays a dominant role for the observed interannual variability of V-Z at 1200 local time.
The Tian Shan range is an inherited intracontinental structure reactivated by the far-field effects of the India-Asia collision. A growing body of thermochronology and magnetostratigraphy datasets shows that the range grew through several tectonic pulses since similar to 25 Ma, however the early Cenozoic history remains poorly constrained. The time-lag between the Eocene India-Asia collision and the Miocene onset of Tian Shan exhumation is particularly enigmatic. This peculiar period is potentially recorded along the southwestern Tian Shan piedmont. There, late Eocene marine deposits of the proto-Paratethys epicontinental sea transition to continental foreland basin sediments of unknown age were recently dated. We provide magnetostratigraphic dating of these continental sediments from the 1700-m-thick Mine section integrated with previously published detrital apatite fission track and U/Pb zircon ages. The most likely correlation to the geomagnetic polarity time scale indicates an age span from 20.8 to 13.3 Ma with a marked increase in accumulation rates at 19-18 Ma. This implies that the entire Oligocene period is missing between the last marine and first continental sediments, as suggested by previous southwestern Tian Shan results. This differs from the southwestern Tarim basin where Eocene marine deposits are continuously overlain by late Eocene-Oligocene continental sediments. This supports a simple evolution model of the western Tarim basin with Eocene-Oligocene foreland basin activation to the south related to northward thrusting of the Kunlun Shan, followed by early Miocene activation of northern foreland basin related to overthrusting of the south Tian Shan. Our data also support southward propagation of the Tian Shan piedmont from 20 to 18 Ma that may relate to motion on the Talas Fergana Fault. The coeval activation of a major right-lateral strike-slip system allowing indentation of the Pamir Salient into the Tarim basin, suggests far-field deformation from the India-Asia collision zone affected the Tian Shan and the Talas Fergana fault by early Miocene. (C) 2015 Elsevier B.V. All rights reserved.
Ground motion with strong-velocity pulses can cause significant damage to buildings and structures at certain periods; hence, knowing the period and velocity amplitude of such pulses is critical for earthquake structural engineering.
However, the physical factors relating the scaling of pulse periods with magnitude are poorly understood.
In this study, we investigate moderate but damaging earthquakes (M-w 6-7) and characterize ground- motion pulses using the method of Shahi and Baker (2014) while considering the potential static-offset effects.
We confirm that the within-event variability of the pulses is large. The identified pulses in this study are mostly from strike-slip-like earthquakes. We further perform simulations using the freq uency-wavenumber algorithm to investigate the causes of the variability of the pulse periods within and between events for moderate strike-slip earthquakes.
We test the effect of fault dips, and the impact of the asperity locations and sizes. The simulations reveal that the asperity properties have a high impact on the pulse periods and amplitudes at nearby stations.
Our results emphasize the importance of asperity characteristics, in addition to earthquake magnitudes for the occurrence and properties of pulses produced by the forward directivity effect.
We finally quantify and discuss within- and between-event variabilities of pulse properties at short distances.
A new view of Ecuador's complex geodynamics has been developed in the course of modeling seismic source zones for probabilistic seismic hazard analysis. This study focuses on two aspects of the plates' interaction at a continental scale: (a) age-related differences in rheology between Farallon and Nazca plates—marked by the Grijalva rifted margin and its inland projection—as they subduct underneath central Ecuador, and (b) the rapidly changing convergence obliquity resulting from the convex shape of the South American northwestern continental margin. Both conditions satisfactorily explain several characteristics of the observed seismicity and of the interseismic coupling. Intermediate-depth seismicity reveals a severe flexure in the Farallon slab as it dips and contorts at depth, originating the El Puyo seismic cluster. The two slabs position and geometry below continental Ecuador also correlate with surface expressions observable in the local and regional geology and tectonics. The interseismic coupling is weak and shallow south of the Grijalva rifted margin and increases northward, with a heterogeneous pattern locally associated to the Carnegie ridge subduction. High convergence obliquity is responsible for the North Andean Block northeastward movement along localized fault systems. The Cosanga and Pallatanga fault segments of the North Andean Block-South American boundary concentrate most of the seismic moment release in continental Ecuador. Other inner block faults located along the western border of the inter-Andean Depression also show a high rate of moderate-size earthquake production. Finally, a total of 19 seismic source zones were modeled in accordance with the proposed geodynamic and neotectonic scheme.
Emerged marine terraces and paleoshorelines along plate margins are prominent geomorphic markers that can be used to quantify the rates and patterns of crustal deformation. The northern margin of the Central Anatolian Plateau has been interpreted as an actively deforming orogenic wedge between the North Anatolian Fault and the Black Sea. Here we use uplifted marine terraces across principal faults on the Sinop Peninsula at the central northern side of the Pontide orogenic wedge to unravel patterns of Quaternary faulting and orogenic wedge behavior. We leveled the present-day elevations of paleoshorelines and dated marine terrace deposits using optically stimulated luminescence (OSL) to determine coastal uplift. The elevations of the paleoshorelines vary between 4 +/- 0.2 and 67 +/- 1.4 m above sea level and OSL ages suggest terrace formation episodes during interglacial periods at ca 125, 190, 400 and 570 ka, corresponding to marine isotopic stages (MIS) 5e, 7a, 11 and 15. Mean apparent vertical displacement rates (without eustatic correction) deduced from these terraces range between 0.02 and 0.18 mm/a, with intermittent faster rates of up to 0.26 mm/a. We obtained higher rates at the eastern and southern parts of the peninsula, toward the hinterland, indicating non-uniform uplift across the different morphotectonic segments of the peninsula. Our data are consistent with active on- and offshore faulting across the Sinop Peninsula. When integrated with regional tectonic observations, the faulting pattern reflects shortening distributed over a broad region of the northern margin of the Central Anatolian Plateau during the Quaternary.
We document Quaternary fluvial incision driven by fault-controlled surface deformation in the inverted intermontane Gökirmak Basin in the Central Pontide mountains along the northern margin of the Central Anatolian Plateau. In-situ-produced Be-10, Ne-21, and Cl-36 concentrations from gravel-covered fluvial terraces and pediment surfaces along the trunk stream of the basin (the Gökirmak River) yield model exposure ages ranging from 71ka to 34645ka and average fluvial incision rates over the past similar to 350ka of 0.280.01mm a(-1). Similarities between river incision rates and coastal uplift rates at the Black Sea coast suggest that regional uplift is responsible for the river incision. Model exposure ages of deformed pediment surfaces along tributaries of the trunk stream range from 605ka to 110 +/- 10ka, demonstrating that the thrust faults responsible for pediment deformation were active after those times and were likely active earlier as well as explaining the topographic relief of the region. Together, our data demonstrate cumulative incision that is linked to active internal shortening and uplift of similar to 0.3mm a(-1) in the Central Pontide orogenic wedge, which may ultimately contribute to the lateral growth of the northern Anatolian Plateau.
Surface uplift at the northern margin of the Central Anatolian Plateau (CAP) is integrally tied to the evolution of the Central Pontides (CP), between the North Anatolian Fault (NAF) and the Black Sea. Our regional morphometric and plate kinematic analyses reveal topographic anomalies, steep channel gradients, and local high relief areas as indicators of ongoing differential surface uplift, which is higher in the western CP compared to the eastern CP and fault-normal components of geodetic slip vectors and the character of tectonic activity of the NAF suggest that stress is accumulated in its broad restraining bend. Seismic reflection and structural field data show evidence for a deep structural detachment horizon responsible for the formation of an actively northward growing orogenic wedge with a positive flower-structure geometry across the CP and the NAF. Taken together, the tectonic, plate kinematic, and geomorphic observations imply that the NAF is the main driving mechanism for wedge tectonics and uplift in the CP. In addition, the NAF Zone defines the boundary between the extensional CAP and the contractional CP. The syntectonic deposits within inverted intermontane basins and deeply incised gorges suggest that the formation of relief, changes in sedimentary dynamics, and > 1 km fluvial incision resulted from accelerated uplift starting in the early Pliocene. The Central Pontides thus provide an example of an accretionary wedge with surface-breaking faults that play a critical role in mountain building processes, sedimentary basin development, and ensuing lateral growth of a continental plateau since the end of the Miocene.
Subduction and collision processes in the Central Andes constrained by converted seismic phases
(2000)
The 26 December 2004 Sumatra-Andaman earthquake of Mw 9.3 triggered a massive tsunami in the Indian Ocean. We here report on observations of the Indian Ocean tsunami at broadband seismic stations located on islands in the area. The tsunami induces long-period (> 1000 s) signals on the horizontal components of the sensor. Frequency-time analysis shows that the long-period signals cannot be due to seismic surface waves, but that it arrives at the expected time of the tsunami. The waveforms are well correlated to tide gauge observations at a location where both observations are available. To explain the signals we favour tilt due to coastal loading but we cannot at the present stage exclude gravitational effects. The density of broadband stations is expected to increase rapidly in the effort of building an earthquake monitoring system. They may unexpectedly become useful tsunami detectors as well
The marine sedimentary record contains unique information about the history of erosion, uplift and climate of the adjacent continent. Inverting this record has been the purpose of many numerical studies. However, limited attention has been given to linking continental erosion to marine sediment transport and deposition in large-scale surface process evolution models. Here we present a new numerical method for marine sediment transport and deposition that is directly coupled to a landscape evolution algorithm solving for the continental fluvial and hillslope erosion equations using implicit and O(N) algorithms. The new method takes into account the sorting of grain sizes (e.g., silt and sand) in the marine domain using a non-linear multiple grain-size diffusion equation and assumes that the sediment flux exported from the continental domain is proportional to the bathymetric slope. Specific transport coefficients and compaction factors are assumed for the two different grain sizes to simulate the stratigraphic architecture. The resulting set of equations is solved using an efficient (O(N) and implicit) algorithm. It can thus be used to invert stratigraphic geometries using a Bayesian approach that requires a large number of simulations. This new method is used to invert the sedimentary geometry of a natural example, the Ogooue Delta (Gabon), over the last similar to 5 Myr. The objective is to unravel the set of erosional histories of the adjacent continental domain compatible with the observed geometry of the offshore delta. For this, we use a Bayesian inversion scheme in which the misfit function is constructed by comparing four geometrical parameters between the natural and the simulated delta: the volume of sediments stored in the delta, the surface slope, the initial and the final shelf lengths. We find that the best-fit values of the transport coefficients for silt in the marine domain are in the range of 300 - 500 m(2)/yr, in agreement with previous studies on offshore diffusion. We also show that, in order to fit the sedimentary geometry, erosion rate on the continental domain must have increased by a factor of 6 to 8 since 5.3 Ma. (C) 2019 Elsevier B.V. All rights reserved.
Southeastern Tibetan Plateau growth revealed by inverse analysis of landscape evolution model
(2022)
The Cenozoic history of the Tibetan Plateau topography is critical for understanding the evolution of the Indian-Eurasian collision, climate, and biodiversity. However, the long-term growth and landscape evolution of the Tibetan Plateau remain ambiguous, it remains unclear if plateau uplift occurred soon after the India-Asia collision in the Paleogene (similar to 50-25 Ma) or later in the Neogene (similar to 20-5 Ma). Here, we reproduce the uplift history of the southeastern Tibetan Plateau using a 2D landscape evolution model, which simultaneously solves fluvial erosion and sediment transport processes in the drainage basins of the Three Rivers region (Yangtze, Mekong, and Salween Rivers). Our model was optimized through a formal inverse analysis with 20,000 forward simulations, which aims to reconcile the transient states of the present-day river profiles. The results, compared to existing paleoelevation and thermochronologic data, suggest initially low elevations (similar to 300-500 m) during the Paleogene, followed by a gradual southeastward propagation of topographic uplift of the plateau margin.
The task of downloading comprehensive datasets of event-based seismic waveforms has been made easier through the development of standardized webservices but is still highly nontrivial because the likelihood of temporary network failures or subtle data errors naturally increases when the amount of requested data is in the order of millions of relatively short segments. This is even more challenging because the typical workflow is not restricted to a single massive download but consists of fetching all possible available input data (e.g., with several repeated download executions) for a processing stage producing any desired user-defined output. Here, we present stream2segment, a highly customizable Python 2+3 package helping the user in the entire workflow of downloading, inspecting, and processing event-based seismic data by means of a relational database management system as archiving storage, which has clear performance and usability advantages, and an integrated processing subroutine requiring a configuration file and a single Python function to produce user-defined output. Stream2segment can also produce diagnostic maps or user-defined plots, which, unlike existing tools, do not require external software dependencies and are not static images but instead are interactive browser-based applications ideally suited for data inspection or annotation tasks and subsequent training of classifiers in foreseen supervised machine-learning applications. Stream2segment has already been used as a data quality tool for datasets within the European Integrated Data Archive and to create a weak-motion database (in the form of a so-called flat file) for the stable continental region of Europe in the context of the European Ground Shaking Intensity Model service, in turn an important building block for seismic hazard studies.
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.