Refine
Has Fulltext
- no (237)
Year of publication
- 2018 (237) (remove)
Document Type
- Article (237) (remove)
Language
- English (237) (remove)
Is part of the Bibliography
- yes (237)
Keywords
- Himalaya (5)
- Germany (4)
- climate change (4)
- permafrost (4)
- Palynology (3)
- Quaternary (3)
- landslides (3)
- Arctic (2)
- Australia (2)
- Chile (2)
Institute
- Institut für Geowissenschaften (237) (remove)
Quantitative estimates of sea-level rise in the Mediterranean Basin become increasingly accurate thanks to detailed satellite monitoring. However, such measuring campaigns cover several years to decades, while longer-term sea-level records are rare for the Mediterranean. We used a data archeological approach to reanalyze monthly mean sea-level data of the Antalya-I (1935–1977) tide gauge to fill this gap. We checked the accuracy and reliability of these data before merging them with the more recent records of the Antalya-II (1985–2009) tide gauge, accounting for an eight-year hiatus. We obtain a composite time series of monthly and annual mean sea levels spanning some 75 years, providing the longest record for the eastern Mediterranean Basin, and thus an essential tool for studying the region's recent sea-level trends. We estimate a relative mean sea-level rise of 2.2 ± 0.5 mm/year between 1935 and 2008, with an annual variability (expressed here as the standard deviation of the residuals, σresiduals = 41.4 mm) above that at the closest tide gauges (e.g., Thessaloniki, Greece, σresiduals = 29.0 mm). Relative sea-level rise accelerated to 6.0 ± 1.5 mm/year at Antalya-II; we attribute roughly half of this rate (~3.6 mm/year) to tectonic crustal motion and anthropogenic land subsidence. Our study highlights the value of data archeology for recovering and integrating historic tide gauge data for long-term sea-level and climate studies.
In this study, the spatial and temporal impacts of the Ataturk Dam on agro-meteorological aspects of the Southeastern Anatolia region have been investigated. Change detection and environmental impacts due to water-reserve changes in Ataturk Dam Lake have been determined and evaluated using multi-temporal Landsat satellite imageries and meteorological datasets within a period of 1984-2011. These time series have been evaluated for three time periods. Dam construction period constitutes the first part of the study. Land cover/use changes especially on agricultural fields under the Ataturk Dam Lake and its vicinity have been identified between the periods of 1984-1992. The second period comprises the 10-year period after the completion of filling up the reservoir in 1992. At this period, Landsat and meteorological time-series analyses are examined to assess the impact of the Ataturk Dam Lake on selected irrigated agricultural areas. For the last 9-year period from 2002 to 2011, the relationships between seasonal water-reserve changes and irrigated plains under changing climatic factors primarily driving vegetation activity (monthly, seasonal, and annual fluctuations of rainfall rate, air temperature, humidity) on the watershed have been investigated using a 30-year meteorological time series. The results showed that approximately 368 km(2) of agricultural fields have been affected because of inundation due to the Ataturk Dam Lake. However, irrigated agricultural fields have been increased by 56.3% of the total area (1552 of 2756 km(2)) on Harran Plain within the period of 1984-2011.
Sub-seasonal thaw slump mass wasting is not consistently energy limited at the landscape scale
(2018)
Predicting future thaw slump activity requires a sound understanding of the atmospheric drivers and geomorphic controls on mass wasting across a range of timescales. On sub-seasonal timescales, sparse measurements indicate that mass wasting at active slumps is often limited by the energy available for melting ground ice, but other factors such as rainfall or the formation of an insulating veneer may also be relevant. To study the sub-seasonal drivers, we derive topographic changes from single-pass radar interferometric data acquired by the TanDEM-X satellites. The estimated elevation changes at 12m resolution complement the commonly observed planimetric retreat rates by providing information on volume losses. Their high vertical precision (around 30 cm), frequent observations (11 days) and large coverage (5000 km(2)) allow us to track mass wasting as drivers such as the available energy change during the summer of 2015 in two study regions. We find that thaw slumps in the Tuktoyaktuk coastlands, Canada, are not energy limited in June, as they undergo limited mass wasting (height loss of around 0 cm day 1) despite the ample available energy, suggesting the widespread presence of early season insulating snow or debris veneer. Later in summer, height losses generally increase (around 3 cm day 1), but they do so in distinct ways. For many slumps, mass wasting tracks the available energy, a temporal pattern that is also observed at coastal yedoma cliffs on the Bykovsky Peninsula, Russia. However, the other two common temporal trajectories are asynchronous with the available energy, as they track strong precipitation events or show a sudden speed-up in late August respectively. The observed temporal patterns are poorly related to slump characteristics like the headwall height. The contrasting temporal behaviour of nearby thaw slumps highlights the importance of complex local and temporally varying controls on mass wasting.
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.
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.
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.
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.
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 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.
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.
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 central Andes
(2018)
The Central Andes and the Atacama Desert represent a unique geological, climatic, and magmatic setting on our planet. It is the only place on Earth where subduction of an oceanic plate below an active continental margin has led to an extensive mountain chain and an orogenic plateau that is second in size only to the Tibetan Plateau, which resulted from continental collision. In this article, we introduce the history of the Central Andes and the evolution of its landscape. We also discuss links between tectonic forces, magmatism, and the extreme hyperarid climate of this land that, in turn, has led to rich deposits of precious ores and minerals.
Dumortierite was synthesized in piston-cylinder experiments at 2.5-4.0 GPa, 650-700 degrees C in the Al2O3 -B2O3-SiO2-H2O (ABSH) system. Electron-microprobe (EMP) analyses reveal significant boron-excess (up to 0.26 B-[4] per formula unit, pfu) and silicon-deficiency relative to the ideal anhydrous dumortierite stoichiometry Al7BSi3O18 . The EMP data in conjunction with results from single-crystal Raman spectroscopy and powder X-ray diffraction provide evidence that silicon at the tetrahedral site is replaced by excess boron via the substitution Si-[4] <--> B-[4] + H. The Raman spectrum of synthetic dumortierite in the frequency region 2000 4000 cm(-1) comprises eight bands, of which six are located at frequencies below 3400 cm(-1). This points to strong hydrogen bonding, most likely O2-H center dot center dot center dot O7 and O7-H center dot center dot center dot O2, arising from a high number of octahedral vacancies at the All site and substitution of trivalent Al3+ and B3+ for Si4+ at Si1 and Si2 sites, causing decreasing acceptor-donor distances and lower incident valence at the acceptor oxygen. Contrary to the synthetic high-pressure ABSH-dumortierite, magnesiodumortierite from the Dora-Maira Massif, which is assumed to have formed at similar conditions (2.5-3.0 GPa, 700 degrees C), does not show any B-excess. Tourmaline shows an analogous behaviour in that magnesium-rich (e.g., dravitic) tourmaline formed at high pressure shows no or only minor amounts of tetrahedral boron, whereas natural aluminum-rich tourmaline and synthetic olenitic tourmaline formed at high pressures can incorporate significant amounts of tetrahedral boron. Two mechanisms might account for this discrepancy: (i) Structural avoidance of Mg-[6]-(OR3+)-R-[4] configurations in magnesiodumortierite due to charge deficieny at the oxygens O2 and O7 and strong local distortion of M1 due to decreased O2-O7 bond length, and/or (ii) decreasing fluid mobility of boron in Al-rich systems at high pressures.
Nitraria is a halophytic taxon (i.e., adapted to saline environments) that belongs to the plant family Nitrariaceae and is distributed from the Mediterranean, across Asia into the south-eastern tip of Australia. This taxon is thought to have originated in Asia during the Paleogene (66-23 Ma), alongside the proto-Paratethys epicontinental sea. The evolutionary history of Nitraria might hold important clues on the links between climatic and biotic evolution but limited taxonomic documentation of this taxon has thus far hindered this line of research. Here we investigate if the pollen morphology and the chemical composition of the pollen wall are informative of the evolutionary history of Nitraria and could explain if origination along the proto-Paratethys and dispersal to the Tibetan Plateau was simultaneous or a secondary process. To answer these questions, we applied a novel approach consisting of a combination of Fourier Transform Infrared spectroscopy (FTIR), to determine the chemical composition of the pollen wall, and pollen morphological analyses using Light Microscopy (LM) and Scanning Electron Microscopy (SEM). We analysed our data using ordinations (principal components analysis and non-metric multidimensional scaling), and directly mapped it on the Nitrariaceae phylogeny to produce a phylomorphospace and a phylochemospace. Our LM, SEM and FTIR analyses show clear morphological and chemical differences between the sister groups Peganum and Nitraria. Differences in the morphological and chemical characteristics of highland species (Nitraria schoberi, N. sphaerocarpa, N. sibirica and N. tangutorum) and lowland species (Nitraria billardierei and N. retusa) are very subtle, with phylogenetic history appearing to be a more important control on Nitraria pollen than local environmental conditions. Our approach shows a compelling consistency between the chemical and morphological characteristics of the eight studied Nitrariaceae species, and these traits are in agreement with the phylogenetic tree. Taken together, this demonstrates how novel methods for studying fossil pollen can facilitate the evolutionary investigation of living and extinct taxa, and the environments they represent.
Proglacial environments are ideal for studying the development of soils through the changes of rocks exposed by glacier retreat to weathering and microbial processes. Carbon (C) and nitrogen (N) contents as well as soil pH and soil elemental compositions are thought to be dominant factors structuring the bacterial, archaeal and fungal communities in the early stages of soil ecosystem formation. However, the functional linkages between C and N contents, soil composition and microbial community structures remain poorly understood. Here, we describe a multivariate analysis of geochemical properties and associated microbial community structures between a moraine and a glaciofluvial outwash in the proglacial area of a High Arctic glacier (Longyearbreen, Svalbard). Our results reveal distinct differences in developmental stages and heterogeneity between the moraine and the glaciofluvial outwash. We observed significant relationships between C and N contents, delta C-13(org) and delta N-15 isotopic ratios, weathering and microbial abundance and community structures. We suggest that the observed differences in microbial and geochemical parameters between the moraine and the glaciofluvial outwash are primarily a result of geomorphological variations of the proglacial terrain.
The impact of reservoir heterogeneities on High-Temperature aquifer thermal energy storage systems
(2018)
We conducted a geoscientific feasibility study for the development of a high-temperature thermal aquifer energy storage system (HT-ATES) outside the capital of Muscat, northern Oman. The aquifer storage is part of a solar geothermal cooling project for the sustainable and continuous cooling of office buildings. The main concept is that excess solar energy will be stored in the subsurface through hot water injection and subsequently utilised as auxiliary energy source during peak demand times. The characterisation of aquifer heterogeneities is thus essential to predict subsurface thermal heat plume development and recovery efficiency of the storage system. We considered two aquifer systems as potential storage horizons, (i) a clastic-dominated alluvial fan system where individual channel systems in combination with diagenetic alterations constitute the main heterogeneities and (ii) a carbonate-dominated system represented by a homogenous layer-cake architecture. The feasibility study included a multidisciplinary approach from initial field work, geocellular reservoir modelling to finite element fluid flow and thermal modelling. Our results show that for the HT-ATES system, with a high frequency of injection and production cycles, heat loss mainly occurs due to heterogeneities in the permeability field of the aquifer in combination with buoyancy driven vertical fluid flow. An impermeable cap-rock is needed to keep the heat plume in place. Conductive heat loss is a minor issue. Highly complex heat plume geometries are apparent in the clastic channel system and ATES well planning is challenging due to the complex and interconnected high permeable channels. The carbonate sequence shows uniform plume geometries due to the layer cake architecture of the system and is tentatively more suitable for ATES development. Based on our findings we propose the general concept of HT-ATES traps, incorporating and building on expertise and knowledge from petroleum and reservoir geology regarding reservoir rocks and suitable trap&seal geometries. The concept can be used as guideline for future high-temperature aquifer storage exploration and development.
The Red Indian Line (RIL) in central Newfoundland is the suture, where the main tract of the Iapetus Ocean was closed at similar to 452 Ma during accretion of the peri-Gondwanan Victoria arc with the composite active Laurentian margin. The protracted deformation history of this soft collision started at similar to 471 Ma with accretion of oceanic terranes to the active composite Laurentian margin. After Iapetus closure both colliding active margins were progressively deformed and metamorphosed during Silurian and Devonian (Salinic, Acadian and Neoacadian orogenic cycles). Peak conditions of the very low- to medium-grade, heterogeneously distributed metamorphism were determined by pseudosection techniques within the range of 2-7 kbar, 230-450 degrees C during increase of the metamorphic field gradient from similar to 12 degrees C/km to similar to 32 degrees C/km over time. Multiple metamorphic crystallisation stages were dated by white mica Ar-40/Ar-39 spot and plateau ages, additional Rb-Sr mineral isochrons involving white mica and one U/Pb age of titanite. All resulting ages between 439 +/- 4 Ma and 356 +/- 16 Ma postdate the closure of Iapetus. Results differ along two transects: The oldest ages of 443-421 Ma (Salinic orogenic cycle) were observed along the northern transect through the RIL zone with minimal younger overprint. Hence low temperature, intermediate to high pressure conditions (4.0-7.0 kbar, 230-340 degrees C) achieved during Taconic-Salinic underthrusting are well preserved. During Acadian dextral transpression the Taconic-Salinic structural wedge was tilted subvertically. In contrast, rocks along the southern transect through the RIL zone mainly show Acadian ages of 408-390 Ma with local preservation of older ages. Acadian deformation occurred under low temperature/low pressure conditions (similar to 250-450 degrees C, 2.5-4.6 kbar). Also Silurian terrestrial cover rocks were buried under these conditions. Acadian-Neoacadian deformation (393-340 Ma) becomes younger towards the northwest and progressively localized in transcurrent fault zones. This final foreland deformation at shallow crustal level established the Acadian/Neoacadian orogenic front in central Newfoundland slightly northwest of the RIL.
An overview is given on the current state of X-ray absorption measurements on silicate melts and glasses. The challenges, limitations, and achievements of analyzing X-ray absorption spectra measured in liquids to determine structural properties of major and minor elements in magmas are described, with particular focus on describing non-Gaussian pair distribution functions in highly disordered glasses and melts, measured at in situ conditions. This includes a discussion on the progress of combining experiments with data from molecular dynamics simulations. For the measurements at conditions of the deep Earth, various experimental approaches and necessities are discussed and two examples are described in more detail. Finally, the achievements and prospects are presented for measuring X-ray absorption spectra indirectly by X-ray Raman scattering.
We performed leaching tests at elevated temperatures and pressures with an Alum black shale from Bomholm, Denmark and a Posidonia black shale from Lower Saxony, Germany. The Alum shale is a carbonate free black shale with pyrite and barite, containing 74.4 mu g/g U. The Posidonia shales is a calcareous shale with pyrite but without detectable amounts of barite containing 3.6 mu g/g U. Pyrite oxidized during the tests forming sulfuric acid which lowered the pH on values between 2 and 3 of the extraction fluid from the Alum shale favoring a release of U from the Alum shale to the fluid during the short-term and in the beginning of the long-term experiments. The activity concentration of U-238 is as high as 23.9 mBq/ml in the fluid for those experiments. The release of U and Th into the fluid is almost independent of pressure. The amount of uranium in the European shales is similar to that of the Marcellus Shale in the United States but the daughter product of U-238, the Ra-226 activity concentrations in the experimentally derived leachates from the European shales are quite low in comparison to that found in industrially derived flowback fluids from the Marcellus shale. This difference could mainly be due to missing Cl in the reaction fluid used in our experiments and a lower fluid to solid ratio in the industrial plays than in the experiments due to subsequent fracking and minute cracks from which Ra can easily be released.
Quantifying rock weakening due to decreasing calcite mineral content by numerical simulations
(2018)
The quantification of changes in geomechanical properties due to chemical reactions is of paramount importance for geological subsurface utilisation, since mineral dissolution generally reduces rock stiffness. In the present study, the effective elastic moduli of two digital rock samples, the Fontainebleau and Bentheim sandstones, are numerically determined based on micro-CT images. Reduction in rock stiffness due to the dissolution of 10% calcite cement by volume out of the pore network is quantified for three synthetic spatial calcite distributions (coating, partial filling and random) using representative sub-cubes derived from the digital rock samples. Due to the reduced calcite content, bulk and shear moduli decrease by 34% and 38% in maximum, respectively. Total porosity is clearly the dominant parameter, while spatial calcite distribution has a minor impact, except for a randomly chosen cement distribution within the pore network. Moreover, applying an initial stiffness reduced by 47% for the calcite cement results only in a slightly weaker mechanical behaviour. Using the quantitative approach introduced here substantially improves the accuracy of predictions in elastic rock properties compared to general analytical methods, and further enables quantification of uncertainties related to spatial variations in porosity and mineral distribution.
As an effort to reduce parameter uncertainties in constructing recurrence plots, and in particular to avoid potential artefacts, this paper presents a technique to derive artefact-safe region of parameter sets. This technique exploits both deterministic (incl. chaos) and stochastic signal characteristics of recurrence quantification (i.e. diagonal structures). It is useful when the evaluated signal is known to be deterministic. This study focuses on the recurrence plot generated from the reconstructed phase space in order to represent many real application scenarios when not all variables to describe a system are available (data scarcity). The technique involves random shuffling of the original signal to destroy its original deterministic characteristics. Its purpose is to evaluate whether the determinism values of the original and the shuffled signal remain closely together, and therefore suggesting that the recurrence plot might comprise artefacts. The use of such determinism-sensitive region shall be accompanied by standard embedding optimization approaches, e.g. using indices like false nearest neighbor and mutual information, to result in a more reliable recurrence plot parameterization.
One main challenge in constructing a reliable recurrence plot (RP) and, hence, its quantification [recurrence quantification analysis (RQA)] of a continuous dynamical system is the induced noise that is commonly found in observation time series. This induced noise is known to cause disrupted and deviated diagonal lines despite the known deterministic features and, hence, biases the diagonal line based RQA measures and can lead to misleading conclusions. Although discontinuous lines can be further connected by increasing the recurrence threshold, such an approach triggers thick lines in the plot. However, thick lines also influence the RQA measures by artificially increasing the number of diagonals and the length of vertical lines [e.g., Determinism (DET) and Laminarity (LAM) become artificially higher]. To take on this challenge, an extended RQA approach for accounting disrupted and deviated diagonal lines is proposed. The approach uses the concept of a sliding diagonal window with minimal window size that tolerates the mentioned deviated lines and also considers a specified minimal lag between points as connected. This is meant to derive a similar determinism indicator for noisy signal where conventional RQA fails to capture. Additionally, an extended local minima approach to construct RP is also proposed to further reduce artificial block structures and vertical lines that potentially increase the associated RQA like LAM. The methodology and applicability of the extended local minima approach and DET equivalent measure are presented and discussed, respectively.
X-ray emission and X-ray Raman scattering spectroscopy are powerful tools to investigate the local electronic and atomic structure of high and low Z elements in situ. Notably, these methods can be applied for in situ spectroscopy at high pressure and high temperature using resistively or laser-heated diamond anvil cells in order to achieve thermodynamic conditions which appear in the Earth's interior. We present a setup for combined X-ray emission and X-ray Raman scattering studies at beamline P01 of PETRA III using a portable wavelength-dispersive von Hamos spectrometer together with the permanently installed multiple-analyzer Johann-type spectrometer. The capabilities of this setup are exemplified by investigating the iron spin crossover of siderite FeCO3 up to 49.3 GPa by measuring the Fe M2,3-edge and the Fe Kβ1,3 emission line simultaneously. With this setup, the Fe valence-to-core emission can be detected together with the Kβ1,3 emission line providing complementary information on the sample's electronic structure. By implementing a laser-heating device, we demonstrate the strength of using a von Hamos type spectrometer for spin state mapping at extreme conditions. Finally, we give different examples of low Z elements' absorption edges relevant for application in geoscience that are accessible with the Johann-type XRS spectrometer. With this setup new insights into the spin transition and compression mechanisms of Earth's mantle materials can be obtained of importance for comprehension of the macroscopic physical and chemical properties of the Earth's interior.
The variation of Rayleigh ellipticity versus frequency is gaining popularity in site characterization. It becomes a necessary observable to complement dispersion curves when inverting shear wave velocity profiles. Various methods have been proposed so far to extract polarization from ambient vibrations recorded on a single three-component station or with an array of three-component sensors. If only absolute values were recovered 10 yr ago, new array-based techniques were recently proposed with enhanced efficiencies providing also the ellipticity sign. With array processing, higher-order modes are often detected even in the ellipticity domain. We suggest to explore the properties of a high-resolution beamforming where radial and vertical components are explicitly included. If N is the number of three-component sensors, 2N x 2N cross-spectral density matrices are calculated for all presumed directions of propagation. They are built with N radial and N vertical channels. As a first approach, steering vectors are designed to fit with Rayleigh wave properties: the phase shift between radial and vertical components is either -Pi/2 or Pi/2. We show that neglecting the ellipticity tilt due to attenuation has only minor effects on the results. Additionally, we prove analytically that it is possible to retrieve the ellipticity value from the usual maximization of the high-resolution beam power. The method is tested on synthetic data sets and on experimental data. Both are reference sites already analysed by several authors. A detailed comparison with previous results on these cases is provided.
We present the first high-resolution trace element (Mg/Ca, Sr/Ca, Ba/Ca) record from a stalagmite in southwestern Romania covering the last 3.6 ka, which provides the potential for quantitative climate reconstruction. Precise age control is based on three independent dating methods, in particular for the last 250 yr, where chemical lamina counting is combined with the identification of the 20th century radiocarbon bomb peak and Th-230/U dating. Long-term cave monitoring and model simulations of drip water and speleothem elemental variability indicate that precipitation-related processes are the main drivers of speleothem Mg/Ca ratios. Calibration against instrumental climate data shows a significant anti-correlation of speleothem Mg/Ca ratios with autumn/winter (October to March) precipitation (r = -0.61, p < 0.01), which is statistically robust when considering age uncertainties and auto-correlation. This relationship is used to develop a quantitative reconstruction of autumn/winter precipitation. During the late Holocene, our data suggest a heterogeneous pattern of past regional winter hydroclimate in the Carpathian/Balkan realm, along with intermittent weakening of the dominant influence of North Atlantic forcing. In agreement with other regional paleo-hydrological reconstructions, the observed variability reveals periodically occurring strong NW-SE hydro-climate gradients. We hypothesize, that this pattern is caused by shifts of the eastern edge of the area of influence of the NAO across central eastern Europe due to the confluence of North Atlantic forcing, and other climatic features such as the East Atlantic/Western Russia (EAWR) pattern. (C) 2018 Elsevier B.V. All rights reserved.
Process life cycle assessment (PLCA) is widely used to quantify environmental flows associated with the manufacturing of products and other processes. As PLCA always depends on defining a system boundary, its application involves truncation errors. Different methods of estimating truncation errors are proposed in the literature; most of these are based on artificially constructed system complete counterfactuals. In this article, we review the literature on truncation errors and their estimates and systematically explore factors that influence truncation error estimates. We classify estimation approaches, together with underlying factors influencing estimation results according to where in the estimation procedure they occur. By contrasting different PLCA truncation/error modeling frameworks using the same underlying input-output (I-O) data set and varying cut-off criteria, we show that modeling choices can significantly influence estimates for PLCA truncation errors. In addition, we find that differences in I-O and process inventory databases, such as missing service sector activities, can significantly affect estimates of PLCA truncation errors. Our results expose the challenges related to explicit statements on the magnitude of PLCA truncation errors. They also indicate that increasing the strictness of cut-off criteria in PLCA has only limited influence on the resulting truncation errors. We conclude that applying an additional I-O life cycle assessment or a path exchange hybrid life cycle assessment to identify where significant contributions are located in upstream layers could significantly reduce PLCA truncation errors.
In this study, we developed an autodetection technique for the equatorial plasma depletions (EPDs) and their occurrence and depletion amplitudes based on in situ electron density measurements gathered by Swarm A satellite. For the first time, comparisons are made among the detected EPDs and their amplitudes with the loss of Global Positioning System (GPS) signal of receivers onboard Swarm A, and the Swarm Level-2 product, Ionospheric Bubble Index (IBI). It has been found that the highest rate of EPD occurrence takes place generally between 2200 and 0000 magnetic local time (MLT), in agreement with the IBI. However, the largest amplitudes of EPD are detected earlier at about 1900-2100 MLT. This coincides with the moment of higher background electron density and the largest occurrence of GPS signal loss. From a longitudinal perspective, the higher depletion amplitude is always witnessed in spatial bins with higher background electron density. At most longitudes, the occurrence rate of postmidnight EPDs is reduced compared to premidnight ones; while more postmidnight EPDs are observed at African longitudes. CHAMP observations confirm this point regardless of high or low solar activity condition. Further by comparing with previous studies and the plasma vertical drift velocity from ROCSAT-1, we suggest that while the F region vertical plasma drift plays a key role in dominating the occurrence of EPDs during premidnight hours, the postmidnight EPDs are the combined results from the continuing of former EPDs and newborn EPDs, especially during June solstice. And these newborn EPDs during postmidnight hours seem to be less related to the plasma vertical drift.
Many hydrological models have been calibrated and validated using hydrographs alone. Because streamflow integrates water fluxes in space, many distributed hydrological models tend to have multiple feasible descriptions of hydrological processes. This equifinality usually leads to substantial prediction uncertainty. In this study, additional constraintsnamely, the spatial patterns of long-term average evapotranspiration (ET), shallow groundwater level, and land cover changewere used to investigate the reduction of equifinality and prediction uncertainty in the Soil and Water Assessment Tool (SWAT) in the Wami River basin in Tanzania. The additional constraints were used in the set-up, parameter emulation and calibration of the SWAT model termed an improved hydrological model (IHM). The IHM was then compared with a classical hydrological model (CHM) that was also developed using the SWAT model but without additional constraints. In the calibration, the CHM used only the hydrograph, but the IHM used the hydrograph and the spatial pattern of long-term average ET as an additional constraint. The IHM produced a single, unique behavioural simulation, whereas the CHM produced many behavioural simulations that resulted in prediction uncertainty. The performance of the IHM with respect to the hydrograph was more consistent than that of the CHM, and the former clearly captured the mean behaviour of ET in the river basin. Therefore, we conclude that additional constraints substantially reduce equifinality and prediction uncertainty in a distributed hydrological model.
High-latitude treeless ecosystems represent spatially highly heterogeneous landscapes with small net carbon fluxes and a short growing season. Reliable observations and process understanding are critical for projections of the carbon balance of the climate-sensitive tundra. Space-borne remote sensing is the only tool to obtain spatially continuous and temporally resolved information on vegetation greenness and activity in remote circumpolar areas. However, confounding effects from persistent clouds, low sun elevation angles, numerous lakes, widespread surface inundation, and the sparseness of the vegetation render it highly challenging. Here, we conduct an extensive analysis of the timing of peak vegetation productivity as shown by satellite observations of complementary indicators of plant greenness and photosynthesis. We choose to focus on productivity during the peak of the growing season, as it importantly affects the total annual carbon uptake. The suite of indicators are as follows: (1) MODIS-based vegetation indices (VIs) as proxies for the fraction of incident photosynthetically active radiation (PAR) that is absorbed (fPAR), (2) VIs combined with estimates of PAR as a proxy of the total absorbed radiation (APAR), (3) sun-induced chlorophyll fluorescence (SIF) serving as a proxy for photosynthesis, (4) vegetation optical depth (VOD), indicative of total water content and (5) empirically upscaled modelled gross primary productivity (GPP). Averaged over the pan-Arctic we find a clear order of the annual peak as APAR <= GPP < SIF < VIs/VOD. SIF as an indicator of photosynthesis is maximised around the time of highest annual temperatures. The modelled GPP peaks at a similar time to APAR. The time lag of the annual peak between APAR and instantaneous SIF fluxes indicates that the SIF data do contain information on light-use efficiency of tundra vegetation, but further detailed studies are necessary to verify this. Delayed peak greenness compared to peak photosynthesis is consistently found across years and land-cover classes. A particularly late peak of the normalised difference vegetation index (NDVI) in regions with very small seasonality in greenness and a high amount of lakes probably originates from artefacts. Given the very short growing season in circumpolar areas, the average time difference in maximum annual photosynthetic activity and greenness or growth of 3 to 25 days (depending on the data sets chosen) is important and needs to be considered when using satellite observations as drivers in vegetation models.
This paper investigates the transferability of calibrated HBV model parameters under stable and contrasting conditions in terms of flood seasonality and flood generating processes (FGP) in five Norwegian catchments with mixed snowmelt/rainfall regimes. We apply a series of generalized (differential) split-sample tests using a 6-year moving window over (i) the entire runoff observation periods, and (ii) two subsets of runoff observations distinguished by the seasonal occurrence of annual maximum floods during either spring or autumn. The results indicate a general model performance loss due to the transfer of calibrated parameters to independent validation periods of -5 to -17%, on average. However, there is no indication that contrasting flood seasonality exacerbates performance losses, which contradicts the assumption that optimized parameter sets for snowmelt-dominated floods (during spring) perform particularly poorly on validation periods with rainfall-dominated floods (during autumn) and vice versa.
The selection of earthquake focal mechanisms (FMs) for stress tensor inversion (STI) is commonly done on a spatial basis, that is, hypocentres. However, this selection approach may include data that are undesired, for example, by mixing events that are caused by different stress tensors when for the STI a single stress tensor is assumed. Due to the significant increase of FM data in the past decades, objective data-driven data selection is feasible, allowing more refined FM catalogues that avoid these issues and provide data weights for the STI routines. We present the application of angular classification with expectation-maximization (ACE) as a tool for data selection. ACE identifies clusters of FM without a priori information. The identified clusters can be used for the classification of the style-of-faulting and as weights of the FM data. We demonstrate that ACE effectively selects data that can be associated with a single stress tensor. Two application examples are given for weighted STI from South America. We use the resulting clusters and weights as a priori information for an STI for these regions and show that uncertainties of the stress tensor estimates are reduced significantly.
Flood damage estimation is a core task in flood risk assessments and requires reliable flood loss models. Identifying the driving factors of flood loss at residential buildings and gaining insight into their relations is important to improve our understanding of flood damage processes. For that purpose, we learn probabilistic graphical models, which capture and illustrate (in-)dependencies between the considered variables. The models are learned based on postevent surveys with flood-affected residents after six flood events, which occurred in Germany between 2002 and 2013. Besides the sustained building damage, the survey data contain information about flooding parameters, early warning and emergency measures, property-level mitigation measures and preparedness, socioeconomic characteristics of the household, and building characteristics. The analysis considers the entire data set with a total of 4,468 cases as well as subsets of the data set partitioned into single flood events and flood types: river floods, levee breaches, surface water flooding, and groundwater floods, to reveal differences in the damaging processes. The learned networks suggest that the flood loss ratio of residential buildings is directly influenced by hydrological and hydraulic aspects as well as by building characteristics and property-level mitigation measures. The study demonstrates also that for different flood events and process types the building damage is influenced by varying factors. This suggests that flood damage models need to be capable of reproducing these differences for spatial and temporal model transfers.
Several thousands of moraine-dammed and supraglacial lakes spread over the Hindu Kush Himalayan (HKH) region, and some have grown rapidly in past decades due to glacier retreat. The sudden emptying of these lakes releases large volumes of water and sediment in destructive glacial lake outburst floods (GLOFs), one of the most publicised natural hazards to the rapidly growing Himalayan population. Despite the growing number and size of glacial lakes, the frequency of documented GLOFs is remarkably constant. We explore this possible reporting bias and offer a new processing chain for establishing a more complete Himalayan GLOF inventory. We make use of the full seasonal archive of Landsat images between 1988 and 2016, and track automatically where GLOFs left shrinking water bodies, and tails of sediment at high elevations. We trained a Random Forest classifier to generate fuzzy land cover maps for 2491 images, achieving overall accuracies of 91%. We developed a likelihood-based change point technique to estimate the timing of GLOFs at the pixel scale. Our method objectively detected ten out of eleven documented GLOFs, and another ten lakes that gave rise to previously unreported GLOFs. We thus nearly doubled the existing GLOF record for a study area covering similar to 10% of the HKH region. Remaining challenges for automatically detecting GLOFs include image insufficiently accurate co-registration, misclassifications in the land cover maps and image noise from clouds, shadows or ice. Yet our processing chain is robust and has the potential for being applied on the greater HKH and mountain ranges elsewhere, opening the door for objectively expanding the knowledge base on GLOF activity over the past three decades.
The ability of water to transport and transform soil materials is one of the main drivers of soil and landscape development. In turn, soil and landscape properties determine how water is distributed in soil landscapes. Understanding the complex dynamics of this co-evolution of soils, landscapes and the hydrological system is fundamental in adapting land management to changes in climate. Soil-Landscape Evolution Models (SLEMs) are used to simulate the development and evolution of soils and landscapes. However, many hydrologic processes, such as preferential flow and subsurface lateral flow, are currently absent in these models. This limits the applicability of SLEMs to improve our understanding of feedbacks in the hydro-pedo-geomorphological system. Implementation of these hydrologic processes in SLEMs faces several complications related to calculation demands, limited methods for linking pedogenic and hydrologic processes, and limited data on quantification of changes in the hydrological system over time. In this contribution, we first briefly review processes and feedbacks in soil-landscape-hydrological systems. Next, we elaborate on the development required to include these processes in SLEMs. We discuss the state-of-the-art knowledge, identify complications, give partial solutions and suggest important future development. The main requirements for incorporating hydrologic processes in SLEMs are: (1) designing a model framework that can deal with varying timescales for different sets of processes, (2) developing and implementing methods for simulating pedogenesis as a function of water flow, (3) improving and implementing knowledge on the evolution and dynamics of soil hydraulic properties over different timescales, and (4) improving the database on temporal changes and dynamics of flow paths.
The phi(ev) is calculated from high-resolution discharge and precipitation data for several rain events with a cumulative precipitation P-cum ranging from less than 5mm to more than 80 mm. Because of the high uncertainty of phi(ev) associated with the hydrograph separation method, phi(ev) is calculated with several methods, including graphical methods, digital filters and a tracer-based method. The results indicate that the hydrological response depends on (theta) over bar (ini): during dry conditions phi(ev) is consistently below 0.1, even for events with high and intense precipitation. Above a threshold of (theta) over bar (ini) = 34 vol % phi(ev) can reach values up to 0.99 but there is a high scatter. Some variability can be explained with a weak correlation of phi(ev) with P-cum and rain intensity, but a considerable part of the variability remains unexplained. It is concluded that threshold-based methods can be helpful to prevent overestimation of the hydrological response during dry catchment conditions. The impact of soil moisture on the hydrological response during wet catchment conditions, however, is still insufficiently understood and cannot be generalized based on the present results.
Effects of upper mantle heterogeneities on the lithospheric stress field and dynamic topography
(2018)
The orientation and tectonic regime of the observed crustal/lithospheric stress field contribute to our knowledge of different deformation processes occurring within the Earth's crust and lithosphere. In this study, we analyze the influence of the thermal and density structure of the upper mantle on the lithospheric stress field and topography. We use a 3-D lithosphere–asthenosphere numerical model with power-law rheology, coupled to a spectral mantle flow code at 300 km depth. Our results are validated against the World Stress Map 2016 (WSM2016) and the observation-based residual topography. We derive the upper mantle thermal structure from either a heat flow model combined with a seafloor age model (TM1) or a global S-wave velocity model (TM2). We show that lateral density heterogeneities in the upper 300 km have a limited influence on the modeled horizontal stress field as opposed to the resulting dynamic topography that appears more sensitive to such heterogeneities. The modeled stress field directions, using only the mantle heterogeneities below 300 km, are not perturbed much when the effects of lithosphere and crust above 300 km are added. In contrast, modeled stress magnitudes and dynamic topography are to a greater extent controlled by the upper mantle density structure. After correction for the chemical depletion of continents, the TM2 model leads to a much better fit with the observed residual topography giving a good correlation of 0.51 in continents, but this correction leads to no significant improvement of the fit between the WSM2016 and the resulting lithosphere stresses. In continental regions with abundant heat flow data, TM1 results in relatively small angular misfits. For example, in western Europe the misfit between the modeled and observation-based stress is 18.3°. Our findings emphasize that the relative contributions coming from shallow and deep mantle dynamic forces are quite different for the lithospheric stress field and dynamic topography.
Lithospheric plates move over the low-viscosity asthenosphere balancing several forces, which generate plate motions. We use a global 3-D lithosphere-asthenosphere model (SLIM3D) with visco-elasto-plastic rheology coupled to a spectral model of mantle flow at 300 km depth to quantify the influence of intra-plate friction and asthenospheric viscosity on plate velocities. We account for the brittle-ductile deformation at plate boundaries (yield stress) using a plate boundary friction coefficient to predict the present-day plate motion and net rotation of the lithospheric plates. Previous modeling studies have suggested that small friction coefficients (mu < 0.1, yield stress similar to 100 MPa) can lead to plate tectonics in models of mantle convection. Here we show that in order to match the observed present-day plate motion and net rotation, the frictional parameter must be less than 0.05. We obtain a good fit with the magnitude and orientation of the observed plate velocities (NUVEL-1A) in a no-net-rotation (NNR) reference frame with mu < 0.05 and a minimum asthenosphere viscosity of similar to 5 . 10(19) Pas to 10(20) Pas. Our estimates of net rotation (NR) of the lith-osphere suggest that amplitudes similar to 0.1-0.2 (degrees/Ma), similar to most observation-based estimates, can be obtained with asthenosphere viscosity cutoff values of similar to 10(19) Pas to 5 . 10(19) Pas and friction coefficients mu < 0.05.
Modern mobile devices (i.e. smartphones and tablet computers) are widespread, everyday tools, which are equipped with a variety of sensors including three-axis magnetometers. Here, we investigate the feasibility and the potential of using such mobile devices to mimic geophysical experiments in the classroom in a table-top setup. We focus on magnetic surveying and present a basic setup of a table-top experiment for collecting three-component magnetic data across well-defined source bodies and structures. Our results demonstrate that the quality of the recorded data is sufficient to address a number of important basic concepts in the magnetic method. The shown examples cover the analysis of magnetic data recorded across different kinds of dipole sources, thus illustrating the complexity of magnetic anomalies. In addition, we analyze the horizontal resolution capabilities using a pair of dipole sources placed at different horizontal distances to each other. Furthermore, we demonstrate that magnetic data recorded with a mobile device can even be used to introduce filtering, transformation, and inversion approaches as they are typically used when processing magnetic data sets recorded for real-world field applications. Thus, we conclude that such table-top experiments represent an easy-to-implement experimental procedure (as student exercise or classroom demonstration) and can provide first hands-on experience in the basic principles of magnetic surveying including the fundamentals of data acquisition, analysis and processing, as well as data evaluation and interpretation.
Flooding is an imminent natural hazard threatening most river deltas, e.g. the Mekong Delta. An appropriate flood management is thus required for a sustainable development of the often densely populated regions. Recently, the traditional event-based hazard control shifted towards a risk management approach in many regions, driven by intensive research leading to new legal regulation on flood management. However, a large-scale flood risk assessment does not exist for the Mekong Delta. Particularly, flood risk to paddy rice cultivation, the most important economic activity in the delta, has not been performed yet. Therefore, the present study was developed to provide the very first insight into delta-scale flood damages and risks to rice cultivation. The flood hazard was quantified by probabilistic flood hazard maps of the whole delta using a bivariate extreme value statistics, synthetic flood hydrographs, and a large-scale hydraulic model. The flood risk to paddy rice was then quantified considering cropping calendars, rice phenology, and harvest times based on a time series of enhanced vegetation index (EVI) derived from MODIS satellite data, and a published rice flood damage function. The proposed concept provided flood risk maps to paddy rice for the Mekong Delta in terms of expected annual damage. The presented concept can be used as a blueprint for regions facing similar problems due to its generic approach. Furthermore, the changes in flood risk to paddy rice caused by changes in land use currently under discussion in the Mekong Delta were estimated. Two land-use scenarios either intensifying or reducing rice cropping were considered, and the changes in risk were presented in spatially explicit flood risk maps. The basic risk maps could serve as guidance for the authorities to develop spatially explicit flood management and mitigation plans for the delta. The land-use change risk maps could further be used for adaptive risk management plans and as a basis for a cost-benefit of the discussed land-use change scenarios. Additionally, the damage and risks maps may support the recently initiated agricultural insurance programme in Vietnam.
Abrupt or gradual?
(2018)
We used a change point analysis on a late Pleistocene-Holocene lake-sediment record from the Chew Bahir basin in the southern Ethiopian Rift to determine the amplitude and duration of past climate transitions. The most dramatic changes occurred over 240 yr (from similar to 15,700 to 15,460 yr) during the onset of the African Humid Period (AHP), and over 990 yr (from similar to 4875 to 3885 yr) during its protracted termination. The AHP was interrupted by a distinct dry period coinciding with the high-latitude Younger Dryas stadial, which had an abrupt onset (less than similar to 100 yr) at similar to 13,260 yr and lasted until similar to 11,730 yr. Wet-dry-wet transitions prior to the AHP may reflect the high-latitude Dansgaard-Oeschger cycles, as indicated by cross-correlation of the potassium record with the NorthGRIP ice core record between similar to 45-20 ka. These findings may contribute to the debates regarding the amplitude, and duration and mechanisms of past climate transitions, and their possible influence on the development of early modern human cultures.
Recent studies indicate that the tropical belt has been expanding during recent decades, which can significantly influence precipitation in subtropical climates. Often the location of the tropical border is identified using the Hadley cell edge (HCE) or the subtropical jet stream (STJ), but most studies concentrated on the zonal-mean state, thereby missing regional impacts. Here we detect longitudinal-resolved trends in STJ cores and HCEs over 1979-2016 in both hemispheres at a higher spatial and temporal resolution than previous studies. Besides pronounced regional trend differences in both sign and magnitude, we show that winter HCE and STJ variability increased in the Mediterranean region and decreased over the American and Asian continents. Rainfall variability in these regions changed likewise, and most of those changes can be explained by changes in HCE/STJ variability. This highlights the importance of understanding future tropical belt changes both regionally and in terms of variability. Plain Language Summary We applied a new network-based method to detect motion of the tropical climate border with longitudinal resolution. Depending on the longitudinal position, there are differences in both direction and magnitude of the border motion. In addition, we demonstrate that the rainfall variability is increasing in the Mediterranean region and decreasing over the American and Asian continents, which can be explained by the variability of the tropical belt location.
Terrestrial cosmogenic nuclide (TCN) concentrations in fluvial sediment, from which denudation rates are commonly inferred, can be affected by hillslope processes. TCN concentrations in gravel and sand may differ if localized, deep-excavation processes (e.g. landslides, debris flows) affect the contributing catchment, whereas the TCN concentrations of sand and gravel tend to be more similar when diffusional processes like soil creep and sheetwash are dominant. To date, however, no study has systematically compared TCN concentrations in different detrital grain-size fractions with a detailed inventory of hillslope processes from the entire catchment. Here we compare concentrations of the TCN Be-10 in 20 detrital sand samples from the Quebrada del Toro (southern Central Andes, Argentina) to a hillslope-process inventory from each contributing catchment. Our comparison reveals a shift from low-slope gullying and scree production in slowly denuding, low-slope areas to steep-slope gullying and landsliding in fast-denuding, steep areas. To investigate whether the nature of hillslope processes (locally excavating or more uniformly denuding) may be reflected in a comparison of the Be-10 concentrations of sand and gravel, we define the normalized sand-gravel index (NSGI) as the Be-10-concentration difference between sand and gravel divided by their summed concentrations. We find a positive, linear relationship between the NSGI and median slope, such that our NSGI values broadly reflect the shift in hillslope processes from low-slope gullying and scree production to steep-slope gullying and landsliding. Higher NSGI values characterize regions affected by steep-slope gullying or landsliding. We relate the large scatter in the relationship, which is exhibited particularly in low-slope areas, to reduced hillslope-channel connectivity and associated transient sediment storage within those catchments. While high NSGI values in well-connected catchments are a reliable signal of deep-excavation processes, hillslope excavation processes may not be reliably recorded by NSGI values where sediment experiences transient storage. (c) 2018 John Wiley & Sons, Ltd.
Recent global warming is pronounced in high-latitude regions (e.g. northern Asia), and will cause the vegetation to change. Future vegetation trends (e.g. the "arctic greening") will feed back into atmospheric circulation and the global climate system. Understanding the nature and causes of past vegetation changes is important for predicting the composition and distribution of future vegetation communities. Fossil pollen records from 468 sites in northern and eastern Asia were biomised at selected times between 40 cal ka bp and today. Biomes were also simulated using a climate-driven biome model and results from the two approaches compared in order to help understand the mechanisms behind the observed vegetation changes. The consistent biome results inferred by both approaches reveal that long-term and broad-scale vegetation patterns reflect global- to hemispheric-scale climate changes. Forest biomes increase around the beginning of the late deglaciation, become more widespread during the early and middle Holocene, and decrease in the late Holocene in fringe areas of the Asian Summer Monsoon. At the southern and southwestern margins of the taiga, forest increases in the early Holocene and shows notable species succession, which may have been caused by winter warming at ca. 7 cal ka bp. At the northeastern taiga margin (central Yakutia and northeastern Siberia), shrub expansion during the last deglaciation appears to prevent the permafrost from thawing and hinders the northward expansion of evergreen needle-leaved species until ca. 7 cal ka bp. The vegetation-climate disequilibrium during the early Holocene in the taiga-tundra transition zone suggests that projected climate warming will not cause a northward expansion of evergreen needle-leaved species.
Along the NE Pamir margin, flights of late Quaternary fluvial terraces span actively deforming fault-related folds. We present detailed results on two terraces dated using optically stimulated luminescence (OSL) and cosmogenic radionuclide Be-10 (CRN) techniques. Quartz OSL dating of two different grain sizes (4-11 mu m and 90-180 mu m) revealed the fine-grain quartz fraction may overestimate the terrace ages by up to a factor of ten. Two-mm, small-aliquot, coarse-grain quartz OSL ages, calculated using the minimum age model, yielded stratigraphically consistent ages within error and dated times of terrace deposition to similar to 9 and similar to 16 ka. We speculate that, in this arid environment, fine-grain samples can be transported and deposited in single, turbid, and (sometimes) night-time floods that prevent thorough bleaching and, thereby, can lead to relatively large residual OSL signals. In contrast, sand in the fluvial system is likely to have a much longer residence time during transport, thereby providing greater opportunities for thorough bleaching. CRN Be-10 depth profiles date the timing of terrace abandonment to similar to 8 and similar to 14 ka: ages that generally agree with the coarse-grain quartz OSL ages. Our new terrace age of similar to 13-14 ka is broadly consistent with other terraces in the region that indicate terrace deposition and subsequent abandonment occurred primarily during glacial-interglacial transitions, thereby suggesting a climatic control on the formation of these terraces on the margins of the Tarim Basin. Furthermore, tectonic shortening rates calculated from these deformed terraces range from similar to 1.2 to similar to 4.6 mm/a and, when combined with shortening rates from other structures in the region, illuminate the late Quaternary basinward migration of deformation to faults and folds along the Pamir-Tian Shan collisional interface.
Tundra be dammed
(2018)
Increasing air temperatures are changing the arctic tundra biome. Permafrost is thawing, snow duration is decreasing, shrub vegetation is proliferating, and boreal wildlife is encroaching. Here we present evidence of the recent range expansion of North American beaver (Castor canadensis) into the Arctic, and consider how this ecosystem engineer might reshape the landscape, biodiversity, and ecosystem processes. We developed a remote sensing approach that maps formation and disappearance of ponds associated with beaver activity. Since 1999, 56 new beaver pond complexes were identified, indicating that beavers are colonizing a predominantly tundra region (18,293km(2)) of northwest Alaska. It is unclear how improved tundra stream habitat, population rebound following overtrapping for furs, or other factors are contributing to beaver range expansion. We discuss rates and likely routes of tundra beaver colonization, as well as effects on permafrost, stream ice regimes, and freshwater and riparian habitat. Beaver ponds and associated hydrologic changes are thawing permafrost. Pond formation increases winter water temperatures in the pond and downstream, likely creating new and more varied aquatic habitat, but specific biological implications are unknown. Beavers create dynamic wetlands and are agents of disturbance that may enhance ecosystem responses to warming in the Arctic.
Although hydrologic models provide hypothesis testing of complex dynamics occurring at catchments, fresh-water quality modeling is still incipient at many subtropical headwaters. In Brazil, a few modeling studies assess freshwater nutrients, limiting policies on hydrologic ecosystem services. This paper aims to compare freshwater quality scenarios under different land-use and land-cover (LULC) change, one of them related to ecosystem-based adaptation (EbA), in Brazilian headwaters. Using the spatially semi-distributed Soil and Water Assessment Tool (SWAT) model, nitrate, total phosphorous (TP) and sediment were modeled in catchments ranging from 7.2 to 1037 km(2). These head-waters were eligible areas of the Brazilian payment for ecosystem services (PES) projects in the Cantareira water supply system, which had supplied water to 9 million people in the Sao Paulo metropolitan region (SPMR). We considered SWAT modeling of three LULC scenarios: (i) recent past scenario (S1), with historical LULC in 1990; (ii) current land-use scenario (S2), with LULC for the period 2010-2015 with field validation; and (iii) future land-use scenario with PES (S2 + EbA). This latter scenario proposed forest cover restoration through EbA following the river basin plan by 2035. These three LULC scenarios were tested with a selected record of rainfall and evapotranspiration observed in 2006-2014, with the occurrence of extreme droughts. To assess hydrologic services, we proposed the hydrologic service index (HSI), as a new composite metric comparing water pollution levels (WPL) for reference catchments, related to the grey water footprint (greyWF) and water yield. On the one hand, water quality simulations allowed for the regionalization of greyWF at spatial scales under LULC scenarios. According to the critical threshold, HSI identified areas as less or more sustainable catchments. On the other hand, conservation practices simulated through the S2 + EbA scenario envisaged not only additional and viable best management practices (BMP), but also preventive decision-making at the headwaters of water supply systems.
Mountain rivers respond to strong earthquakes by rapidly aggrading to accommodate excess sediment delivered by co-seismic landslides. Detailed sediment budgets indicate that rivers need several years to decades to recover from seismic disturbances, depending on how recovery is defined. We examine three principal proxies of river recovery after earthquake-induced sediment pulses around Pokhara, Nepal's second largest city. Freshly exhumed cohorts of floodplain trees in growth position indicate rapid and pulsed sedimentation that formed a fan covering 150 km2 in a Lesser Himalayan basin with tens of metres of debris between the 11th and 15th centuries AD. Radiocarbon dates of buried trees are consistent with those of nearby valley deposits linked to major medieval earthquakes, such that we can estimate average rates of re-incision since. We combine high-resolution digital elevation data, geodetic field surveys, aerial photos, and dated tree trunks to reconstruct geomorphic marker surfaces. The volumes of sediment relative to these surfaces require average net sediment yields of up to 4200 t km–2 yr–1 for the 650 years since the last inferred earthquake-triggered sediment pulse. The lithological composition of channel bedload differs from that of local bedrock, confirming that rivers are still mostly evacuating medieval valley fills, locally incising at rates of up to 0.2 m yr–1. Pronounced knickpoints and epigenetic gorges at tributary junctions further illustrate the protracted fluvial response; only the distal portions of the earthquake-derived sediment wedges have been cut to near their base. Our results challenge the notion that mountain rivers recover speedily from earthquakes within years to decades. The valley fills around Pokhara show that even highly erosive Himalayan rivers may need more than several centuries to adjust to catastrophic perturbations. Our results motivate some rethinking of post-seismic hazard appraisals and infrastructural planning in active mountain regions.
Major Countries
(2018)
Management intensity modifies soil properties, e.g., organic carbon (C-org) concentrations and soil pH with potential feedbacks on plant diversity. These changes might influence microbial P concentrations (P-mic) in soil representing an important component of the Pcycle. Our objectives were to elucidate whether abiotic and biotic variables controlling P-mic concentrations in soil are the same for forests and grasslands, and to assess the effect of region and management on P-mic concentrations in forest and grassland soils as mediated by the controlling variables. In three regions of Germany, Schwabische Alb, Hanich-Dun, and Schorfheide-Chorin, we studied forest and grassland plots (each n=150) differing in plant diversity and land-use intensity. In contrast to controls of microbial biomass carbon (C-mic), P-mic was strongly influenced by soil pH, which in turn affected phosphorus (P) availability and thus microbial Puptake in forest and grassland soils. Furthermore, P-mic concentrations in forest and grassland soils increased with increasing plant diversity. Using structural equation models, we could show that soil C-org is the profound driver of plant diversity effects on P-mic in grasslands. For both forest and grassland, we found regional differences in P-mic attributable to differing environmental conditions (pH, soil moisture). Forest management and tree species showed no effect on P-mic due to a lack of effects on controlling variables (e.g., C-org). We also did not find management effects in grassland soils which might be caused by either compensation of differently directed effects across sites or by legacy effects of former fertilization constraining the relevance of actual practices. We conclude that variables controlling P-mic or C-mic in soil differ in part and that regional differences in controlling variables are more important for P-mic in soil than those induced by management.
Alkaline igneous rocks are relatively rare in settings of tectonic convergence and little is known about their petrogenesis in these settings. This study aims to contribute to a better understanding of the formation of alkaline igneous rocks by an investigation of the Tezhsar volcano-intrusive alkaline ring complex (TAC) in the Armenian Lesser Caucasus, which is located between the converging Eurasian and Arabian plates. We present new petrological, geochemical and Sr Nd isotope data for the TAC to constrain magma genesis and magma source characteristics. Moreover, we provide a new Ar-40/Ar-39 age of 41.0 +/- 0.5 Ma on amphibole from a nepheline syenite that is integrated into the regional context of ongoing regional convergence and widespread magmatism. The TAC is spatially concentric and measures similar to 10 km in diameter representing the relatively shallow plumbing system of a major stratovolcano juxtaposed by ring faulting with its extrusive products. The plutonic units comprise syenites and nepheline syenites, whereas the extrusive units are dominated by trachytic-phonolitic rocks. The characteristic feature of the TAC is the development of pseudomorphs after leucite in all types of the volcanic, subvolcanic and intrusive alkaline rocks. Whole-rock major element data show a metaluminous (Alkalinity Index = 0-0.1), alkalic and silica-undersaturated (Feldspathoid Silica-Saturation Index <0) character of the TAC. The general trace element enrichment and strong fractionation of REEs (La-N/Yb-N up to 70) indicate a relatively enriched magma source and small degrees of partial melting. All TAC rocks show a negative Nb Ta anomalies typical of subduction zone settings. The initial 87Sr/85Sr ratios (0.704-0.705) and positive sNd values (+3 to +5) indicate an isotopically depleted upper mantle and lack of significant crustal influence, which in turn suggests the TAC magma has formed via differentiation from lithospheric mantle melts.
The equatorial electrojet occasionally reverses during morning and afternoon hours, leading to periods of westward current in the ionospheric E region that are known as counter electrojet (CEJ) events. We present the first analysis of CEJ climatology and CEJ dependence on solar flux and lunar phase for the Brazilian sector, based on an extensive ground-based data set for the years 2008 to 2017 from the geomagnetic observatory Tatuoca (1.2 degrees S, 48.5 degrees W), and we compare it to the results found for Huancayo (12.0 degrees S, 75.3 degrees W) observatory in the Peruvian sector. We found a predominance of morning CEJ events for both sectors. The afternoon CEJ occurrence rate in the Brazilian sector is twice as high as in the Peruvian sector. The afternoon CEJ occurrence rate strongly depends on season, with maximum rates occurring during the northern-hemisphere summer for the Brazilian sector and during the northern-hemisphere winter for the Peruvian sector. Significant discrepancies between the two sectors are also found for morning CEJ rates during the northern-hemisphere summer. These longitudinal differences are in agreement with a CEJ climatology derived from contemporary Swarm satellite data and can be attributed in part to the well-known longitudinal wave-4 structure in the background equatorial electrojet strength that results from nonmigrating solar tides and stationary planetary waves. Simulations with the Thermosphere-Ionosphere-Electrodynamics General Circulation Model show that the remaining longitudinal variability in CEJ during northern summer can be explained by the effect of migrating tides in the presence of the varying geomagnetic field in the South Atlantic Anomaly.
Snow meltwaters account for most of the yearly water budgets of many catchments in High Mountain Asia (HMA). We examine trends in snow water equivalent (SWE) using passive microwave data (1987 to 2009). We find an overall decrease in SWE in HMA, despite regions of increased SWE in the Pamir, Kunlun Shan, Eastern Himalaya, and Eastern Tien Shan. Although the average decline in annual SWE across HMA (contributing area, 2641 x 10(3) km(2)) is low (average, -0.3%), annual SWE losses conceal distinct seasonal and spatial heterogeneities across the study region. For example, the Tien Shan has seen both strong increases in winter SWE and sharp declines in spring and summer SWE. In the majority of catchments, the most negative SWE trends are found in mid-elevation zones, which often correspond to the regions of highest snow-water storage and are somewhat distinct from glaciated areas. Negative changes in SWE storage in these mid-elevation zones have strong implications for downstream water availability.
The landscape of the semiarid Pampa in central Argentina is characterized by late Pleistocene aeolian deposits, covering large plains with sporadic dune structures. Since the current land use changed from extensive livestock production within the Caldenal forest ecosystem to arable land, the wind erosion risk increased distinctly. We measured wind erosion and deposition patterns at the plot scale and investigated the spatial variability of the erosion processes. The wind-induced mass-transport was measured with 18 Modified Wilson and Cooke samplers (MWAC), installed on a 1.44 ha large field in a 20 x 40 m grid. Physical and chemical soil properties from the upper soil as well as a digital elevation model were recorded in a 20 x 20 m grid. In a 5-month measuring campaign data from seven storms with three different wind directions was obtained. Results show very heterogeneous patterns of erosion and deposition for each storm and indicate favoured erosion on windward and deposits on leeward terrain positions. Furthermore, a multiple regression model was build, explaining up to 70% of the spatial variance of erosion by just using four predictors: topsoil thickness, relative elevation, soil organic carbon content and slope direction. Our findings suggest a structure-process-structure complex where the landscape structure determines the effects of recent wind erosion processes which again slowly influence the structure, leading to a gradual increase of soil heterogeneity.
The lower atmospheric forcing effects on the ionosphere are particularly evident during extreme meteorological events known as sudden stratospheric warmings (SSWs). During SSWs, the polar stratosphere and ionosphere, two distant atmospheric regions, are coupled through the SSW-induced modulation of atmospheric migrating and nonmigrating tides. The changes in the migrating semidiurnal solar and lunar tides are the major source of ionospheric variabilities during SSWs. In this study, we use 55 years of ground-magnetometer observations to investigate the composite characteristics of the lunar tide of the equatorial electrojet (EEJ) during SSWs. These long-term observations allow us to capture the EEJ lunar tidal response to the SSWs in a statistical sense. Further, we examine the influence of solar flux conditions and the phases of quasi-biennial oscillation (QBO) on the lunar tide and find that the QBO phases and solar flux conditions modulate the EEJ lunar tidal response during SSWs in a similar way as they modulate the wintertime Arctic polar vortex. This work provides first evidence of modulation of the EEJ lunar tide due to QBO. Plain Language Summary This study focuses on the vertical coupling between the polar stratosphere and equatorial ionosphere during sudden stratospheric warmings (SSWs). Extreme meteorological events such as SSWs induce variabilities in the ionosphere by modulating the atmospheric migrating and nonmigrating tides, and these variabilities can be comparable to a moderate geomagnetic storm. Observations and modeling studies have found that the changes in the migrating semidiurnal solar and lunar tides are a major source of ionospheric variabilities during SSWs. The equatorial electrojet (EEJ) is a narrow ribbon of current flowing over the dip equator in the ionosphere and is particularly sensitive to tidal changes. Long-term ground-magnetometer recordings have been used in this study to estimate the variations induced in EEJ during SSWs due to the lunar semidiurnal tide in a statistical sense. The wintertime Arctic polar vortex and the occurrence of SSWs are modulated by solar flux conditions and the phases of quasi-biennial oscillation. In this work, we find the first evidence of lunar tidal modulation of EEJ due to quasi-biennial oscillation during SSWs. Our findings will be useful in providing improved predictions of ionospheric variations due to SSWs. The aeronomy community will be the most impacted by this paper.
Originally developed for use in the petroleum industry, Rock-Eval pyrolysis is a technique commonly applied to lake sediments to infer paleoenvironmental reconstructions. The standard Rock-Eval parameters provide information on the amount of total organic and inorganic carbon (TOC and MinC, respectively), and are usually interpreted as proxies for the source (aquatic or terrestrial) of the primary production of organic matter (Hydrogen Index vs Oxygen Index). Although this method usually provides valuable evidence, the common presence of siderite in tropical lake sediments can alter the primary signal of the sedimentary organic matter (SOM). Indeed, the CO2 and CO released by the pyrolysis of siderite are integral to the calculation of the SOM-related standard Rock-Eval parameters. In this study, we analyze sediments from a core collected in the Lake Barombi (southwest Cameroon) and describe the impact of siderite on standard Rock-Eval parameters. We propose a workflow that allows standard Rock-Eval parameters to be corrected, based on the analysis of thermograms. The proposed corrections provide siderite-effect-free parameters, accurately reflecting the changes in sedimentary organic matter composition. (C) 2018 Elsevier Ltd. All rights reserved.
The lateral and vertical temperature distribution in Oman is so far only poorly understood, particularly in the area between Muscat and the Batinah coast, which is the area of this study and which is composed of Cenozoic sediments developed as part of a foreland basin of the Makran Thrust Zone. Temperature logs (T-logs) were run and physical rock properties of the sediments were analyzed to understand the temperature distribution, thermal and hydraulic properties, and heat-transport processes within the sedimentary cover of northern Oman. An advective component is evident in the otherwise conduction-dominated geothermal play system, and is caused by both topography and density driven flow. Calculated temperature gradients (T-gradients) in two wells that represent conductive conditions are 18.7 and 19.5 °C km−1, corresponding to about 70–90 °C at 2000–3000 m depth. This indicates a geothermal potential that can be used for energy intensive applications like cooling or water desalinization. Sedimentation in the foreland basin was initiated after the obduction of the Semail Ophiolite in the late Campanian, and reflects the complex history of alternating periods of transgressive and regressive sequences with erosion of the Oman Mountains. Thermal and hydraulic parameters were analyzed of the basin’s heterogeneous clastic and carbonate sedimentary sequence. Surface heat-flow values of 46.4 and 47.9 mW m−2 were calculated from the T-logs and calculated thermal conductivity values in two wells. The results of this study serve as a starting point for assessing different geothermal applications that may be suitable for northern Oman.
Landslide hazard motivates the need for a deeper understanding of the events that occur before, during, and after catastrophic slope failures. Due to the destructive nature of such events, in situ observation is often difficult or impossible. Here, we use data from a network of 58 seismic stations to characterise a large landslide at the Askja caldera, Iceland, on 21 July 2014. High data quality and extensive network coverage allow us to analyse both long- and short-period signals associated with the landslide, and thereby obtain information about its triggering, initiation, timing, and propagation. At long periods, a landslide force history inversion shows that the Askja landslide was a single, large event starting at the SE corner of the caldera lake at 23:24:05 UTC and propagating to the NW in the following 2 min The bulk sliding mass was 7-16 x 10(10) kg, equivalent to a collapsed volume of 35-80 x 10(6) m(3). The sliding mass was displaced downslope by 1260 +/- 250 m. At short periods, a seismic tremor was observed for 30 min before the landslide. The tremor is approximately harmonic with a fundamental frequency of 2.3 Hz and shows time-dependent changes of its frequency content. We attribute the seismic tremor to stick-slip motion along the landslide failure plane. Accelerating motion leading up to the catastrophic slope failure culminated in an aseismic quiescent period for 2 min before the landslide. We propose that precursory seismic signals may be useful in landslide early-warning systems. The 8 h after the main landslide failure are characterised by smaller slope failures originating from the destabilised caldera wall decaying in frequency and magnitude. We introduce the term "afterslides" for this subsequent, declining slope activity after a large landslide.
Plain Language Summary The 2015 Gorkha earthquake in Nepal caused severe losses in the hydropower sector. The country temporarily lost similar to 20% of its hydropower capacity, and >30 hydropower projects were damaged. The projects hit hardest were those that were affected by earthquake-triggered landslides. We show that these projects are located along very steep rivers with towering sidewalls that are prone to become unstable during strong seismic ground shaking. A statistical classification based on a topographic metric that expresses river steepness and earthquake ground acceleration is able to approximately predict hydropower damage during future earthquakes, based on successful testing of past cases. Thus, our model enables us to estimate earthquake damages to hydropower projects in other parts of the Himalayas. We find that >10% of the Himalayan drainage network may be unsuitable for hydropower infrastructure given high probabilities of high earthquake damages.
The Alpine Fault is capable of generating large (MW > 8) earthquakes and is the main geohazard on South Island, NZ, and late in its 250–291-year seismic cycle. To minimize its hazard potential, it is indispensable to identify and understand the processes influencing the geomechanical behavior and strength-evolution of the fault. High-resolution microstructural, mineralogical and geochemical analyses of the Alpine Fault's core demonstrate wall rock fragmentation, assisted by mineral dissolution, and cementation resulting in the formation of a fine-grained principal slip zone (PSZ). A complex network of anastomosing and mutually cross-cutting calcite veins implies that faulting occurred during episodes of dilation, slip and sealing. Fluid-assisted dilatancy leads to a significant volume increase accommodated by vein formation in the fault core. Undeformed euhedral chlorite crystals and calcite veins that have cut footwall gravels demonstrate that these processes occurred very close to the Earth's surface. Microstructural evidence indicates that cataclastic processes dominate the deformation and we suggest that powder lubrication and grain rolling, particularly influenced by abundant nanoparticles, play a key role in the fault core's velocity-weakening behavior rather than frictional sliding. This is further supported by the absence of smectite, which is reasonable given recently measured geothermal gradients of more than 120 °C km−1 and the impermeable nature of the PSZ, which both limit the growth of this phase and restrict its stability to shallow depths. Our observations demonstrate that high-temperature fluids can influence authigenic mineral formation and thus control the fault's geomechanical behavior and the cyclic evolution of its strength.
Intercomparison of cosmic-ray neutron sensors and water balance monitoring in an urban environment
(2018)
Sensor-to-sensor variability is a source of error common to all geoscientific instruments that needs to be assessed before comparative and applied research can be performed with multiple sensors. Consistency among sensor systems is especially critical when subtle features of the surrounding terrain are to be identified. Cosmic-ray neutron sensors (CRNSs) are a recent technology used to monitor hectometre-scale environmental water storages, for which a rigorous comparison study of numerous co-located sensors has not yet been performed. In this work, nine stationary CRNS probes of type "CRS1000" were installed in relative proximity on a grass patch surrounded by trees, buildings, and sealed areas. While the dynamics of the neutron count rates were found to be similar, offsets of a few percent from the absolute average neutron count rates were found. Technical adjustments of the individual detection parameters brought all instruments into good agreement. Furthermore, we found a critical integration time of 6 h above which all sensors showed consistent dynamics in the data and their RMSE fell below 1% of gravimetric water content. The residual differences between the nine signals indicated local effects of the complex urban terrain on the scale of several metres. Mobile CRNS measurements and spatial simulations with the URANOS neutron transport code in the surrounding area (25 ha) have revealed substantial sub-footprint heterogeneity to which CRNS detectors are sensitive despite their large averaging volume. The sealed and constantly dry structures in the footprint furthermore damped the dynamics of the CRNS-derived soil moisture. We developed strategies to correct for the sealed-area effect based on theoretical insights about the spatial sensitivity of the sensor. This procedure not only led to reliable soil moisture estimation during dry-out periods, it further revealed a strong signal of intercepted water that emerged over the sealed surfaces during rain events. The presented arrangement offered a unique opportunity to demonstrate the CRNS performance in complex terrain, and the results indicated great potential for further applications in urban climate research.
Following the widespread assumption that a majority of ubiquitous marine microplastic particles originate from land-based sources, recent studies identify rivers as important pathways for microplastic particles (MPP) to the oceans. Yet a detailed understanding of the underlying processes and dominant sources is difficult to obtain with the existing accurate but extremely time-consuming methods available for the identification of MPP. Thus in the presented study, a novel approach applying short-wave infrared imaging spectroscopy for the quick and semi-automated identification of MPP is applied in combination with a multitemporal survey concept. Volume-reduced surface water samples were taken from transects at ten points along a major watercourse running through the South of Berlin, Germany, on six dates. After laboratory treatment, the samples were filtered onto glass fiber filters, scanned with an imaging spectrometer and analyzed by image processing. The presented method allows to count MPP, classify the plastic types and determine particle sizes. At the present stage of development particles larger than 450 m in diameter can be identified and a visual validation showed that the results are reliable after a subsequent visual final check of certain typical error types. Therefore, the method has the potential to accelerate microplastic identification by complementing FTIR and Raman microspectroscopy. Technical advancements (e.g. new lens) will allow lower detection limits and a higher grade of automatization in the near future. The resulting microplastic concentrations in the water samples are discussed in a spatio-temporal context with respect to the influence (i) of urban areas, (ii) of effluents of three major Berlin wastewater treatment plants discharging into the canal and (iii) of precipitation events. Microplastic concentrations were higher downstream of the urban area and after precipitation. An increase in microplastic concentrations was discernible for the wastewater treatment plant located furthest upstream though not for the other two. (C) 2018 Elsevier Ltd. All rights reserved.
The potential link between erosion rates at the Earth’s surface and changes in global climate has intrigued geoscientists for decades1,2 because such a coupling has implications for the influence of silicate weathering3,4 and organic-carbon burial5 on climate and for the role of Quaternary glaciations in landscape evolution1,6. A global increase in late-Cenozoic erosion rates in response to a cooling, more variable climate has been proposed on the basis of worldwide sedimentation rates7. Other studies have indicated, however, that global erosion rates may have remained steady, suggesting that the reported increases in sediment-accumulation rates are due to preservation biases, depositional hiatuses and varying measurement intervals8,9,10. More recently, a global compilation of thermochronology data has been used to infer a nearly twofold increase in the erosion rate in mountainous landscapes over late-Cenozoic times6. It has been contended that this result is free of the biases that affect sedimentary records11, although others have argued that it contains biases related to how thermochronological data are averaged12 and to erosion hiatuses in glaciated landscapes13. Here we investigate the 30 locations with reported accelerated erosion during the late Cenozoic6. Our analysis shows that in 23 of these locations, the reported increases are a result of a spatial correlation bias—that is, combining data with disparate exhumation histories, thereby converting spatial erosion-rate variations into temporal increases. In four locations, the increases can be explained by changes in tectonic boundary conditions. In three cases, climatically induced accelerations are recorded, driven by localized glacial valley incision. Our findings suggest that thermochronology data currently have insufficient resolution to assess whether late-Cenozoic climate change affected erosion rates on a global scale. We suggest that a synthesis of local findings that include location-specific information may help to further investigate drivers of global erosion rates.
Changes in topography on Earth, particularly the growth of major mountain belts like the Central Andes, have a fundamental impact on regional and global atmospheric circulation patterns. These patterns, in turn, affect processes such as precipitation, erosion, and sedimentation. Over the last two decades, various geochemical, geomorphologic, and geologic approaches have helped identify when, where, and how quickly topography has risen in the past. The current spatio-temporal picture of Central Andean growth is now providing insight into which deep-Earth processes have left their imprint on the shape of the Earth's surface.
In areas characterized by many volcanoes, identifying the source of a deposit may not be trivial. This becomes much more complicated when looking for the source of a debris avalanche deposit (DAD), the common products of catastrophic volcanic edifice collapses. To overcome this problem, in this work a methodology is proposed based on the integration of texture features and areal distribution of the deposit, comparison between the petrography of the coarser clasts within the DAD and of the proximal products, grain-size analysis, and the volumetric estimations of the deposit and the volume missing from the volcanic edifice. This methodology has been tested to a DAD occurred near the city of Ibarra (Imbabura Province; Northern Ecuador), having a controversial source. Two main volcanic edifice are located in proximity of the DAD, the Cubilche volcano (3826 m.a.s.l.), located immediately south of and east of the colossal dormant Imbabura volcano. The former displays a sharp horseshoe shaped scar towards the north and inside this post-collapse edifice, that we name old Cubilche volcano (OCV), is located the young Cubilche volcano (YCV) that refilled a portion of the collapse scar and partially covered the southern flank of the OCV. Detailed knowledge of Cubilche volcano is critical because of its close proximity and interspersed activity with Imbabura volcano. In fact, Imbabura most recent edifice was built over the northwestern slope of the OCV and partially covered it. Recent studies linked the studied DAD to both Imbabura volcano as a product of its northern sector collapse, as well as neighboring Cubilche volcano. Our data points to Cubilche as the most likely source for this DAD. A perspective view of the shaded relief image of the present day OCV shows that the morphology of the volcano is well-preserved on its southern, eastern, and western flanks. This allows us to reconstruct the morphology of the OCV previous to the collapse through interpolation of elevation and altitude data of preserved flanks. A DEM of the present day topography was used for extrapolating the morphology. Using similar methodology, the post collapse base of the amphitheater was reconstructed by removing the relief of the present day YCV. The reconstructed topography of the OCV shows that it could have been a symmetric cone, reaching a maximum elevation of similar to 4100 m.a.s.l. with a lack volume of similar to 3.5 km(3). Based on this scenario, the deposit originated from the OCV main collapse should have a volume >3-3.5 km(3) in accordance to the volume calculated for the studied DAD. (C) 2018 Elsevier B.V. All rights reserved.
During the last years Urban Green Infrastructure (UGI) has evolved as a research focus across Europe. UGI can be understood as a multifunctional network of different urban green spaces and elements contributing to urban benefits. Urban agriculture has gained increasing research interest in this context. While a strong focus has been made on functions and benefits of small scale activities, the question is still open, whether these findings can be up-scaled and transferred to the farmland scale. Furthermore, multifunctionality of urban and peri-urban agriculture is rarely being considered in the landscape context. This research aims to address these gaps and harnesses the question if agricultural landscapes – which in many European metropolitan regions provide significant spatial potential – can contribute to UGI as multifunctional green spaces. This work considers multifunctionality qualitatively based on stakeholder opinion, using a participatory research approach. This study provides new insights in peri-urban farmland potentials for UGI development, resulting into a strategy framework. Furthermore, it reflects on the role of the stakeholder involvement for `multifunctionality planning´. It suggests that it helps to define meaningful bundles of intertwined functions that interact on different scales, helping to deal with non-linearity of multiple functions and to better manage them simultaneously.
After a comprehensive geophysical prospecting the Quaternary MA 1/2 tina Maar, located on a line between the two Quaternary scoria cones Komorni could be revealed by a scientific drilling at the German-Czech border in 2007. Further geophysical field investigations led to the discovery of another geological structure about 2.5 km ESE of the small town Neualbenreuth (NE-Bavaria, Germany), inferred to be also a maar structure, being the fourth volcanic feature aligned along the NW-SE trending Tachov fault zone. It is only faintly indicated as a partial circular rim in the digital elevation model. Though not expressed by a clear magnetic anomaly, geoelectric and refraction seismic tomography strongly indicates a bowl-shaped depression filled with low-resistivity and low-velocity material, correlating well with the well-defined negative gravity anomaly of - 2.5 mGal. Below ca. 15 m-thick debris layer, successions of mostly laminated sediments were recovered in a 100 m-long sediment core in 2015. Sections of finely laminated layers, likely varves, rich in organic matter and tree pollen, were recognized in the upper (22-30 m) and lower (70-86 m) part of the core, respectively, interpreted as interglacials, whereas mostly minerogenic laminated deposits, poor in organic matter, and (almost) barren of tree pollen are interpreted as clastic glacial deposits. According to a preliminary age model based on magnetostratigraphy, palynology, radiocarbon dating, and cyclostratigraphy, the recovered sediments span the time window from about 85 ka back to about 270 ka, covering marine isotope stages 5-8. Sedimentation rates are in the range of 10 cm ka(-1) in interglacials and up to 100 cm ka(-1) in glacial phases. The stratigraphic record resembles the one from MA 1/2 tina Maar, with its eruption date being derived from a nearby tephra deposit at 288 +/- 17 ka, thus supporting the age model of the inferred Neualbenreuth Maar.
The volumes, rates and grain size distributions of sediment supplied from hillslopes represent the initial input of sediment delivered from upland areas and propagated through sediment routing systems. Moreover, hillslope sediment supply has a significant impact on landscape response time to tectonic and climatic perturbations. However, there are very few detailed field studies characterizing hillslope sediment supply as a function of lithology and delivery process. Here, we present new empirical data from tectonically-active areas in southern Italy that quantifies how lithology and rock strength control the landslide fluxes and grain size distributions supplied from hillslopes. Landslides are the major source of hillslope sediment supply in this area, and our inventory of similar to 2800 landslides reveals that landslide sediment flux is dominated by small, shallow landslides. We find that lithology and rock strength modulate the abundance of steep slopes and landslides, and the distribution of landslide sizes. Outcrop-scale rock strength also controls the grain sizes supplied by bedrock weathering, and influences the degree of coarsening of landslide supply with respect to weathering supply. Finally, we show that hillslope sediment supply largely determines the grain sizes of fluvial export, from catchments and that catchments with greater long-term landslide rates deliver coarser material. Therefore, our results demonstrate a dual control of lithology on hillslope sediment supply, by modulating both the sediment fluxes from landslides and the grain sizes supplied by hillslopes to the fluvial system.
The Timpa delle Murge ophiolite in the North Calabrian Unit is part of the Liguride Complex (southern Apennines). The study is concentrated on the gabbroic part of the ophiolite of the Pollino area. They preserve the high-grade ocean floor metamorphic and locally developed flaser textures under ocean floor conditions. The primary magmatic assemblages are clinopyroxene, plagioclase, and opaques. Brown hornblende is a late magmatic phase. Green hornblende, actinolite, albite, chlorite and epidote display metamorphic recrystallization under lower amphibolite facies conditions, followed by greenschist facies. The gabbros show subalkaline near to alkaline character with a tendency to a more calkalkaline trend. The normalization to primitive mantle and mid-ocean ridge basalt (N-MORB) compositions indicates a considerable depletion in Nb, P, Zr and Ti and an enrichment in Ba, Rb, K, Sr and Eu. This points to a mantle source, which is not compatible with a "normal" mid-ocean ridge situation. Rather, the gabbros are generated from a N-MORB-like melt with a strong crustal component, which was influenced by subduction related fluids and episodic melting during mid-ocean-ridge processes. Plausible geodynamic settings of the Timpa delle Murge gabbros are oceanic back-arc positions with embryonic MORB-activities. Similar slab contaminated magmatism is also known from the early stage of island arc formation in supra-subduction zone environments like the Izu-Bonin-Mariana island arc.
Global patterns of declining temperature variability from the Last Glacial Maximum to the Holocene
(2018)
Changes in climate variability are as important for society to address as are changes in mean climate(1). Contrasting temperature variability during the Last Glacial Maximum and the Holocene can provide insights into the relationship between the mean state of the climate and its variability(2,3). However, although glacial-interglacial changes in variability have been quantified for Greenland(2), a global view remains elusive. Here we use a network of marine and terrestrial temperature proxies to show that temperature variability decreased globally by a factor of four as the climate warmed by 3-8 degrees Celsius from the Last Glacial Maximum (around 21,000 years ago) to the Holocene epoch (the past 11,500 years). This decrease had a clear zonal pattern, with little change in the tropics (by a factor of only 1.6-2.8) and greater change in the mid-latitudes of both hemispheres (by a factor of 3.3-14). By contrast, Greenland ice-core records show a reduction in temperature variability by a factor of 73, suggesting influences beyond local temperature or a decoupling of atmospheric and global surface temperature variability for Greenland. The overall pattern of reduced variability can be explained by changes in the meridional temperature gradient, a mechanism that points to further decreases in temperature variability in a warmer future.
Reading the sediment record in terms of past climates is challenging since linking climate change to the associated responses of sedimentary systems is not always straightforward. Here we analyze the erosional response of landscapes on the Tibetan Plateau to interglacial climate forcing. Using the theory of dynamical systems on Holocene time series of geochemical proxies, we derive a sedimentary response model that accurately simulates observed proxy variation in three lake records. The model suggests that millennial variations in sediment composition reflect a self-organization of landscapes in response to abrupt climate change between 11.6 and 11.9 ka BP. The self-organization is characterized by oscillations in sediment supply emerging from a feedback between physical and chemical erosion processes, with estimated response times between 3,000 to 18,000 years depending on catchment topography. The implications of our findings emphasize the need for landscape response models to decipher the paleoclimatic code in continental sediment records. Plain Language Summary Lake sediments are an important source of information on past climates. Reading the information is not always straightforward. Complex interactions in landscapes can affect the transmission of climatic signals to the sediment record. However, the exact nature of such complex interactions remains unknown. This study compares sediment deposits of three lakes on the Tibetan Plateau. The deposits are continuous records of landscape responses to climate change during the last 12,000 years. We identified a mathematical model that accurately simulates changes in sediment composition at all sites. The model simulations suggest that an abrupt warming at the end of the last glacial period destabilized the landscapes. This caused fluctuations in the transport of sediments, which persisted for several thousand years. Our findings present evidence for a long-lasting impact of abrupt climate change on fundamental Earth surface processes.
Reading the sediment record in terms of past climates is challenging since linking climate change to the associated responses of sedimentary systems is not always straightforward. Here we analyze the erosional response of landscapes on the Tibetan Plateau to interglacial climate forcing. Using the theory of dynamical systems on Holocene time series of geochemical proxies, we derive a sedimentary response model that accurately simulates observed proxy variation in three lake records. The model suggests that millennial variations in sediment composition reflect a self-organization of landscapes in response to abrupt climate change between 11.6 and 11.9 ka BP. The self-organization is characterized by oscillations in sediment supply emerging from a feedback between physical and chemical erosion processes, with estimated response times between 3,000 to 18,000 years depending on catchment topography. The implications of our findings emphasize the need for landscape response models to decipher the paleoclimatic code in continental sediment records. Plain Language Summary Lake sediments are an important source of information on past climates. Reading the information is not always straightforward. Complex interactions in landscapes can affect the transmission of climatic signals to the sediment record. However, the exact nature of such complex interactions remains unknown. This study compares sediment deposits of three lakes on the Tibetan Plateau. The deposits are continuous records of landscape responses to climate change during the last 12,000 years. We identified a mathematical model that accurately simulates changes in sediment composition at all sites. The model simulations suggest that an abrupt warming at the end of the last glacial period destabilized the landscapes. This caused fluctuations in the transport of sediments, which persisted for several thousand years. Our findings present evidence for a long-lasting impact of abrupt climate change on fundamental Earth surface processes.
Retrogressive thaw slumps (RTSs) are among the most active thermokarst landforms in the Arctic and deliver a large amount of material to the Arctic Ocean. However, their contribution to the organic carbon (OC) budget is unknown. We provide the first estimate of the contribution of RTSs to the nearshore OC budget of the Yukon Coast, Canada, and describe the evolution of coastal RTSs between 1952 and 2011 in this area. We (1) describe the evolution of RTSs between 1952 and 2011; (2) calculate the volume of eroded material and stocks of OC mobilized through slumping, including soil organic carbon (SOC) and dissolved organic carbon (DOC); and (3) estimate the OC fluxes mobilized through slumping between 1972 and 2011. We identified RTSs using high- resolution satellite imagery from 2011 and geocoded aerial photographs from 1952 and 1972. To estimate the volume of eroded material, we applied spline interpolation on an airborne lidar dataset acquired in July 2013. We inferred the stocks of mobilized SOC and DOC from existing related literature. Our results show a 73% increase in the number of RTSs and 14% areal expansion between 1952 and 2011. In the study area, RTSs displaced at least 16.6 x 10(6) m(3) of material, 53% of which was ice, and mobilized 145.9 x 10(6) kg of OC. Between 1972 and 2011, 49 RTSs displaced 8.6 x 10(3) m(3) yr(-1) of material, adding 0.6% to the OC flux released by coastal retreat along the Yukon Coast. Our results show that the contribution of RTSs to the nearshore OC budget is non-negligible and should be included when estimating the quantity of OC released from the Arctic coast to the ocean.
Chemostat experiments are employed to study predator-prey and other trophic interactions, frequently using phytoplankton-zooplankton systems. These experiments often use population dynamics as fingerprints of ecological and evolutionary processes, assuming that the contributions of all major actors to these dynamics are known. However, bacteria are often neglected although they are frequently present. We argue that even without external carbon input bacteria may affect the experimental outcomes depending on experimental conditions and the physiological traits of bacteria, phytoplankton, and zooplankton. Using a static carbon flux model and a dynamic simulation model, we predict the minimum and maximum impact of bacteria on phytoplankton-zooplankton population dynamics. Under bacteria-suppressing conditions, we find that the effect of bacteria is indeed negligible and their omission justified. Under bacteria-favoring conditions, however, bacteria may strongly affect average biomasses of phytoplankton and zooplankton. The population dynamics may become highly complex, which may result in wrong interpretations when inferring processes (e.g., trait changes) from population dynamic patterns without considering bacteria. We provide suggestions to reduce the bacterial impact experimentally. Besides optimizing experimental conditions (e.g., the dilution rate) the appropriate choice of the zooplankton predator is decisive. Counterintuitively, bacteria have a larger impact if the predator is not bacterivorous as high bacterial biomasses and complex population dynamics arise via competition for nutrients with the phytoplankton. Only at least partial bacterivory minimizes the impact of bacteria. Our results help to improve the design of chemostat experiments and their interpretation, and advance the study of ecological and evolutionary processes in aquatic food webs.
The Sierra Nevada de Santa Marta in NW Colombia is an isolated massif at the northernmost end of the Andes chain near the boundary with the Caribbean plate. Previous geologic mapping and K-Ar dating have shown that Jurassic plutonic and volcanic units make up a large part of the Santa Marta Massif (SMM). These rocks have been considered to be part of a Jurassic magmatic arc extending from NW Colombia to northern Chile, but without any geochemical basis for comparison. This paper reports on a geochemical and Sr-Nd-Pb isotope study of the Jurassic rocks in the SMM and provides 12 new U-Pb zircon ages from in-situ laser ICP-MS dating. The plutonic and volcanic units span a range from 45 to 78 wt.% SiO2, with a dominance of intermediate to felsic compositions with SiO2 > 57 wt.%. They classify as calc-alkaline, medium to high-K, metaluminous rocks with trace-element features typical for arc-derived magma series. In terms of their major and trace-element compositions, the SMM Jurassic units overlap with contemporary plutonic and volcanic rocks from other regions of the Central and Eastern Cordilleras of Colombia, and confirm an arc affinity. The new U-Pb ages range from 176 +/- 1 Ma to 192 +/- 2 Ma (n = 12), with most between 180 and 188 Ma (n = 7). The initial Sr isotope ratios (at 180 Ma) are between 0.7012 and 0.7071 (n = 29), with 3 outliers attributed to mobilization of Rb and/or Sr, Nearly all samples have negative( )epsilon Nd-(180) values between - 10.3 and 0.0 (n = 30), the two exceptions being only slightly positive (1.1 and 1.9). Measured Pb isotope ratios fall in a narrow range, with Pb-206/Pb-204 from 18.02 to 19.95, (207) Pb/(204) Pb from 15.56 to 15.67 and Pb-208/Pb-204 from 37.76 to 39.04 (n = 28). In the regional context of previous studies, these results confirm early Jurassic ages and an arc affinity for the widespread magmatism exposed in the eastern and northeastern Colombian Andes. We also note patterns in the distribution and composition of magmas. The magmatic activity in the Central Cordillera tends to be younger than in the Eastern Cordillera and is spatially more restricted to the vicinity of regional fault systems. In terms of composition, Jurassic igneous rocks in the Eastern Cordillera have systematically lower epsilon Nd-(180) values than those from the Central Cordillera, whereas the Pb isotope ratios overlap. We ascribe the Nd isotope variations to heterogeneity in the mantle source and/or degree of crustal contamination, whereas the Pb isotope ratios are crust-dominated and similar throughout the region. The spatio-temporal and compositional evolution of Jurassic magmatic rocks in the Northern Andes reflect the major plate kinematic readjustment between the Triassic and the Early Jurassic in the proto-Andean margin.
In the arctic and high mountains it is common to measure vertical changes of ice sheets and glaciers via digital elevation model (DEM) differencing. This requires the signal of change to outweigh the noise associated with the datasets. Excluding large landslides, on the ice-free earth the land-level change is smaller in vertical magnitude and thus requires more accurate DEMs for differencing and identification of change. Previously, this has required meter to submeter data at small spatial scales. Following careful corrections, we are able to measure land-level changes in gravel-bed channels and steep hillslopes in the south-central Andes using the SRTM-C (collected in 2000) and the TanDEM-X (collected from 2010 to 2015) near-global 12–30m DEMs. Long-standing errors in the SRTM-C are corrected using the TanDEM-X as a control surface and applying cosine-fit co-registration to remove ∼ 1∕10 pixel (∼ 3m) shifts, fast Fourier transform (FFT) and filtering to remove SRTM-C short- and long-wavelength stripes, and blocked shifting to remove remaining complex biases. The datasets are then differenced and outlier pixels are identified as a potential signal for the case of gravel-bed channels and hillslopes. We are able to identify signals of incision and aggradation (with magnitudes down to ∼ 3m in the best case) in two > 100km river reaches, with increased geomorphic activity downstream of knickpoints. Anthropogenic gravel excavation and piling is prominently measured, with magnitudes exceeding ±5m (up to > 10m for large piles). These values correspond to conservative average rates of 0.2 to > 0.5myr−1 for vertical changes in gravel-bed rivers. For hillslopes, since we require stricter cutoffs for noise, we are only able to identify one major landslide in the study area with a deposit volume of 16±0.15×106m3. Additional signals of change can be garnered from TanDEM-X auxiliary layers; however, these are more difficult to quantify. The methods presented can be extended to any region of the world with SRTM-C and TanDEM-X coverage where vertical land-level changes are of interest, with the caveat that remaining vertical uncertainties in primarily the SRTM-C limit detection in steep and complex topography.
The additional magnetic field produced by the ionospheric current system is a part of the Earth’s magnetic field. This current system is a highly variable part of a global electric circuit. The solar wind and interplanetary magnetic field (IMF) interaction with the Earth’s magnetosphere is the external driver for the global electric circuit in the ionosphere. The energy is transferred via the field-aligned currents (FACs) to the Earth’s ionosphere. The interactions between the neutral and charged particles in the ionosphere lead to the so-called thermospheric neutral wind dynamo which represents the second important driver for the global current system. Both processes are components of the magnetosphere–ionosphere–thermosphere (MIT) system, which depends on solar and geomagnetic conditions, and have significant seasonal and UT variations.
The modeling of the global dynamic Earth’s ionospheric current system is the first aim of this investigation. For our study, we use the Potsdam version of the Upper Atmosphere Model (UAM-P). The UAM is a first-principle, time-dependent, and fully self-consistent numerical global model. The model includes the thermosphere, ionosphere, plasmasphere, and inner magnetosphere as well as the electrodynamics of the coupled MIT system for the altitudinal range from 80 (60) km up to the 15 Earth radii. The UAM-P differs from the UAM by a new electric field block. For this study, the lower latitudinal and equatorial electrodynamics of the UAM-P model was improved.
The calculation of the ionospheric current system’s contribution to the Earth’s magnetic field is the second aim of this study. We present the method, which allows computing the additional magnetic field inside and outside the current layer as generated by the space current density distribution using the Biot-Savart law. Additionally, we perform a comparison of the additional magnetic field calculation using 2D (equivalent currents) and 3D current distribution.
Glaciers and rivers are the main agents of mountain erosion. While in the fluvial realm empirical relationships and their mathematical description, such as the stream power law, improved the understanding of fundamental controls on landscape evolution, simple constraints on glacial topography and governing scaling relations are widely lacking. We present a steady state solution for longitudinal profiles along eroding glaciers in a coupled system that includes tectonics and climate. We combined the shallow ice approximation and a glacial erosion rule to calculate ice surface and bed topography from prescribed glacier mass balance gradient and rock uplift rate. Our approach is inspired by the classic application of the stream power law for describing a fluvial steady state but with the striking difference that, in the glacial realm, glacier mass balance is added as an altitude-dependent variable. From our analyses we find that ice surface slope and glacial relief scale with uplift rate with scaling exponents indicating that glacial relief is less sensitive to uplift rate than relief in most fluvial landscapes. Basic scaling relations controlled by either basal sliding or internal deformation follow a power law with the exponent depending on the exponents for the glacial erosion rule and Glen's flow law. In a mixed scenario of sliding and deformation, complicated scaling relations with variable exponents emerge. Furthermore, a cutoff in glacier mass balance or cold ice in high elevations can lead to substantially larger scaling exponents which may provide an explanation for high relief in high latitudes.
Near-surface geophysical imaging of alluvial fan settings is a challenging task but crucial for understating geological processes in such settings. The alluvial fan of Ghor Al-Haditha at the southeast shore of the Dead Sea is strongly affected by localized subsidence and destructive sinkhole collapses, with a significantly increasing sinkhole formation rate since ca. 1983. A similar increase is observed also on the western shore of the Dead Sea, in correlation with an ongoing decline in the Dead Sea level. Since different structural models of the upper 50 m of the alluvial fan and varying hypothetical sinkhole processes have been suggested for the Ghor Al-Haditha area in the past, this study aimed to clarify the subsurface characteristics responsible for sinkhole development. For this purpose, high-frequency shear wave reflection vibratory seismic surveys were carried out in the Ghor Al-Haditha area along several crossing and parallel profiles with a total length of 1.8 and 2.1 km in 2013 and 2014, respectively. The sedimentary architecture of the alluvial fan at Ghor Al-Haditha is resolved down to a depth of nearly 200 m at a high resolution and is calibrated with the stratigraphic profiles of two boreholes located inside the survey area. The most surprising result of the survey is the absence of evidence of a thick (>2-10 m) compacted salt layer formerly suggested to lie at ca. 35-40 m depth. Instead, seismic reflection amplitudes and velocities image with good continuity a complex interlocking of alluvial fan deposits and lacustrine sediments of the Dead Sea between 0 and 200 m depth. Furthermore, the underground section of areas affected by sinkholes is characterized by highly scattering wave fields and reduced seismic interval velocities. We propose that the Dead Sea mud layers, which comprise distributed inclusions or lenses of evaporitic chloride, sulfate, and carbonate minerals as well as clay silicates, become increasingly exposed to unsaturated water as the sea level declines and are consequently destabilized and mobilized by both dissolution and physical erosion in the subsurface. This new interpretation of the underlying cause of sinkhole development is supported by surface observations in nearby channel systems. Overall, this study shows that shear wave seismic reflection technique is a promising method for enhanced near-surface imaging in such challenging alluvial fan settings.
Carbonate-rich silicate and carbonate melts play a crucial role in deep Earth magmatic processes and their melt structure is a key parameter, as it controls physical and transport properties. Carbon-rich melts can be strongly enriched in trace elements, but the structural incorporation mechanisms of these elements are difficult to study because such melts generally cannot be quenched to glasses. In this contribution we investigate the influence of CO2 on the local environments of trace elements contained in silicate glasses with variable CO2 concentrations and in silicate and carbonate melts. The melts were studied in-situ at high pressure and temperature conditions using the Paris-Edinburgh press (2.2 to 2.6 GPa and 1200 to 1500 degrees C). The compositions studied include sodium-rich peralkaline silicate melts and glasses and carbonate melts similar to those occurring naturally at Oldoinyo Lengai volcano. The local environments of yttrium (Y), lanthanum (La) and strontium (Sr) were investigated using extended X-ray absorption fine structure (EXAFS) spectroscopy. Main findings of the study suggest: (1) In peralkaline silicate glasses the local structure of Y is unaffected by the CO2 content. Contrary, a slight increase of oxygen bond lengths of Sr and La is inferred with increasing CO2 content in peralkaline glasses, while they remain constant in glasses of even higher peralkalinity independent of the CO2 content. (2) In silicate melts of different CO2 contents Y-O bond lengths are constant, while a slight increase within carbonate melt compositions is deduced. On the other hand, a steady bond lengths increase over the whole compositional range is inferred for La-O and Sr-O. This may well be explained by distinct preferences of these elements for specific local environments. Based on these new data, we suggest potential mechanisms for the structural incorporation of these elements, a key step towards understanding their partitioning behavior in natural magmatic systems.
To evaluate the spatiotemporal variations of ground motions in northern Chile, we built a high-quality rock seismic acceleration database and an interface earthquakes catalog. Two ground-motion prediction equation (GMPE) models for subduction zones have been tested and validated for the area. They were then used as backbone models to describe the time-space variations of earthquake frequency content (Fourier and response spectra). Consistent with previous studies of large subduction earthquakes, moderate interface earthquakes in northern Chile show an increase of the high-frequency energy released with depth. A regional variability of earthquake frequency content is also observed, which may be related to a lateral segmentation of the mechanical properties of the subduction interface. Finally, interface earthquakes show a temporal evolution of their frequency content in the earthquake sequence associated with the 2014 Iquique M-w 8.1 megathrust earthquake. Surprisingly, the change does not occur with the mainshock but is associated with an 8 month slow slip preceding the megathrust. Electronic Supplement: Strong-motion database.
Two ground motion prediction equation models for subduction zones have been tested using a public ground motion database of the KiK-net records obtained by automated processing protocols (Dawood et al., 2016, https://doi.org/10.1193/071214EQS106). The database contains records of more than 700 interface earthquakes that occurred on the Japan subduction between 1998 and 2012. The Zhao et al. (2006, https://doi.org/10.1785/0120050122) ground motion prediction equation was shown to be the best suited model for the region. It was then used as backbone to analyze the variability of ground motion records. The residuals between observed and predicted ground motions have been analyzed to study the spatial variation of the earthquakes' ground motion frequency content on the Japan megathrust. This analysis revealed a depth dependency of generated ground motions consistent with the downdip segmentation proposed for subduction interfaces (Lay et al., 2012, https://doi.org/10.1029/2011JB009133), a regional ground motion dependency that may be related with lateral variations of the mechanical properties of the subduction interface and a high-frequency radiations drop in the earthquake sequence that preceded the Tohoku-Oki earthquake Mw 9.0. The regional ground motion dependency suggests the existence of different domains along trench of the Japan subduction megathrust that control the ground motions and the wave radiation patterns of interface earthquakes. The location of their boundaries is consistent with the extension of the rupture of the 2011 Tohoku-Oki earthquake, with pre-Tohoku interseismic coupling, and with the free air gravity anomalies.
Quantifying the pace of ice-sheet growth is critical to understanding ice-age climate and dynamics. Here, we show that the diversion of the Hudson River (northeastern North America) late in the last glaciation phase (ca. 30 ka), which some previous studies have speculated was due to glacial isostatic adjustment (GIA), can be used to infer the timing of the Laurentide Ice Sheet’s growth to its maximum extent. Landscapes in the vicinity of glaciated regions have likely responded to crustal deformation produced by ice-sheet growth and decay through river drainage reorganization, given that rates of uplift and subsidence are on the order of tens of meters per thousand years. We perform global, gravitationally self-consistent simulations of GIA and input the predicted crustal deformation field into a landscape evolution model. Our calculations indicate that the eastward diversion of the Hudson River at 30 ka is consistent with exceptionally rapid growth of the Laurentide Ice Sheet late in the glaciation phase, beginning at 50–35 ka.
Relative abundances of 157 diatom taxa from Yakutian lake surface-sediments were investigated for their potential to indicate certain environmental conditions. Data from 206 sites from Arctic, sub-Arctic and boreal environments were included. Redundancy analyses were performed to assess the explanatory power of mean July temperature (T-July), conductivity, pH, dissolved silica concentration, phosphate concentration, lake depth and vegetation type on diatom species composition. Boosted regression tree analyses were performed to infer the most relevant environmental variables for abundances of individual taxa and weighted average regression was applied to infer their respective optimum and tolerance. Electrical conductivity was best indicated by diatom taxa. In contrast, only few taxa were indicative of Si and water depth. Few taxa were related to specific pH values. Although T-July, explained the highest proportion of variance in the diatom spectra and was, after conductivity, the second-most selected splitting variable, we a priori decided not to present indicator taxa because of the poorly understood relationship between diatom occurrences and T-July. In total, 92 diatom taxa were reliable indicators of a certain vegetation type or a combination of several types. The high numbers of indicative species for open vegetation sites and for forested sites suggest that the principal turnover is the transition from forest-tundra to northern taiga. Overall, our results reveal that preference ranges of diatom taxa for environmental variables are mostly broad, and the use of indicator taxa for the purposes of environmental reconstruction or environmental monitoring is therefore restricted to marked rather than subtle environmental transitions.
Nearly 13,000 years ago, the warming trend into the Holocene was sharply interrupted by a reversal to near glacial conditions. Climatic causes and ecological consequences of the Younger Dryas (YD) have been extensively studied, however proxy archives from the Mediterranean basin capturing this period are scarce and do not provide annual resolution. Here, we report a hydroclimatic reconstruction from stable isotopes (delta O-18, delta C-13) in subfossil pines from southern France. Growing before and during the transition period into the YD (12 900-12 600 cal BP), the trees provide an annually resolved, continuous sequence of atmospheric change. Isotopic signature of tree sourcewater (delta O-18(sw)) and estimates of relative air humidity were reconstructed as a proxy for variations in air mass origin and precipitation regime. We find a distinct increase in inter-annual variability of sourcewater isotopes (delta O-18(sw)), with three major downturn phases of increasing magnitude beginning at 12 740 cal BP. The observed variation most likely results from an amplified intensity of North Atlantic (low delta O-18(sw)) versus Mediterranean (high delta O-18(sw)) precipitation. This marked pattern of climate variability is not seen in records from higher latitudes and is likely a consequence of atmospheric circulation oscillations at the margin of the southward moving polar front.
Quaternary glaciations have repeatedly shaped large tracts of the Andean foreland. Its spectacular large glacial lakes, staircases of moraine ridges, and extensive outwash plains have inspired generations of scientists to reconstruct the processes, magnitude, and timing of ice build-up and decay at the mountain front. Surprisingly few of these studies noticed many dozens of giant (≥108 m3) mass-wasting deposits in the foreland. We report some of the world's largest terrestrial landslides in the eastern piedmont of the Patagonian Ice Sheet (PIS) along the traces of the former Lago Buenos Aires and Lago Puyerredón glacier lobes and lakes. More than 283 large rotational slides and lateral spreads followed by debris slides, earthflows, rotational and translational rockslides, complex slides and few large rock avalanches detached some 164 ± 56 km3 of material from the slopes of volcanic mesetas, lake-bounding moraines, and river-gorge walls. Many of these landslide deposits intersect with well-dated moraine ridges or former glacial-lake shorelines, and offer opportunities for relative dating of slope failure. We estimate that >60% of the landslide volume (∼96 km3) detached after the Last Glacial Maximum (LGM). Giant slope failures cross-cutting shorelines of a large Late Glacial to Early Holocene lake (“glacial lake PIS”) likely occurred during successive lake-level drop between ∼11.5 and 8 ka, and some of them are the largest hitherto documented landslides in moraines. We conclude that 1) large portions of terminal moraines can fail catastrophically several thousand years after emplacement; 2) slopes formed by weak bedrock or unconsolidated glacial deposits bordering glacial lakes can release extremely large landslides; and 3) landslides still occur in the piedmont, particularly along postglacial gorges cut in response to falling lake levels.
Flash floods and debris flows are iconic hazards inmountainous regions with steep relief, high rainfall intensities, rapid snowmelt events, and abundant sediments. The cuesta landscapes of southern Germany hardly come to mind when dealing with such hazards. A series of heavy rainstorms dumping up to 140mm in 2 h caused destructive flash floods and debris flows in May 2016. The most severe damage occurred in the Braunsbach municipality, which was partly buried by 42,000 m(3) of boulders, gravel, mud, and anthropogenic debris from the small catchment of Orlacher Bach (similar to 6 km(2)). We analysed this event by combining rainfall patterns, geological conditions, and geomorphic impacts to estimate an average sediment yield of 14,000 t/km(2) that mostly (similar to 95%) came from some 50 riparian landslides and channel-bed incision of similar to 2 m. This specific sediment yield ranks among the top 20% globally, while the intensity-duration curve of the rainstormis similarly in the upper percentile range of storms that had triggered landslides. Compared to similar-sized catchments in the greater region hit by the rainstorms, we find that the Orlacher Bach is above the 95th percentile in terms of steepness, storm-rainfall intensity, and topographic curvatures. The flash flood transported a sediment volume equal to as much as 20-40% of the Pleistocene sediment volume stored in the Orlacher Bach fan, andmay have had several predecessors in the Holocene. River control structures from 1903 and records of a debris flow in the 1920s in a nearby catchment indicate that the local inhabitants may have been aware of the debris-flow hazards earlier. Such recurring and destructive events elude flood-hazard appraisals in humid landscapes of gentle relief, and broaden mechanistic views of how landslides and debris flows contribute to shaping small and deeply cut tributaries in the southern Germany cuesta landscape.
Evidence for Quaternary climate change in East Africa has been derived from outcrops on land and lake cores and from marine dust, leaf wax, and pollen records. These data have previously been used to evaluate the impact of climate change on hominin evolution, but correlations have proved to be difficult, given poor data continuity and the great distances between marine cores and terrestrial basins where fossil evidence is located. Here, we present continental coring evidence for progressive aridification since about 575 thousand years before present (ka), based on Lake Magadi (Kenya) sediments. This long-term drying trend was interrupted by many wet-dry cycles, with the greatest variability developing during times of high eccentricity-modulated precession. Intense aridification apparent in the Magadi record took place between 525 and 400 ka, with relatively persistent arid conditions after 350 ka and through to the present. Arid conditions in the Magadi Basin coincide with the Mid-Brunhes Event and overlap with mammalian extinctions in the South Kenya Rift between 500 and 400 ka. The 525 to 400 ka arid phase developed in the South Kenya Rift between the period when the last Acheulean tools are reported (at about 500 ka) and before the appearance of Middle Stone Age artifacts (by about 320 ka). Our data suggest that increasing Middle- to Late-Pleistocene aridification and environmental variability may have been drivers in the physical and cultural evolution of Homo sapiens in East Africa.
The Northern Central Iranian Micro-continent (CIM) represents Neotethys-related oceanic crust remnants, emplaced due to convergence between CIM and Eurasia plates during Eocene. Mafic and ultramafic units are exposed along the northern part of the CIM in the Sabzevar area. The geology and field relation of Sabzevar ophiolite indicate northward subduction of the Sabzevar basin. The average whole rock chemistry of mafic (gabbros) and ultramafic samples (lherzolite, harzburgite and dunite) is characterized by a range of MgO of 11.16-31.88, CaO 5.22-11.53 and Al2O3 2.77-14.57, respectively. Low LREE/HREE ratio of ultramafic samples is accompanied by enrichment of large ion lithophile elements (LILE) such as Sr, Pb and K. Mafic samples show two distinct groups with low and high LREE/HREE ratios. The spider diagram of mafic samples indicates enrichment in Sr, Pb and K and depletion in REE. Petrological and geochemical evidence and field relations show that the mafic rocks formed in a supra-subduction zone setting. Petrological studies reveal the role of fractional crystallization and assimilation effect by released fluids during subduction related generation of the Sabzevar mafic rocks. We suggest that the studied mafic rocks likely represent the basement of an initial island arc, which was generated in a supra-subduction zone setting within the Neotethys branch of the Sabzevar Ocean at the north of CIM. Copper, gold and chromite mineralizations are studied in relation to island arc setting and supra-subduction environment. Similarities in lithology, ophiolite age and mineralization between Sabzevar ophiolite and Bardaskan-Torbat Heydariyeh ophiolites testify for their separation due to rotation (or faulting) of the Central Iranian Micro-continent.
Mapping Damage-Affected Areas after Natural Hazard Events Using Sentinel-1 Coherence Time Series
(2018)
The emergence of the Sentinel-1A and 1B satellites now offers freely available and widely accessible Synthetic Aperture Radar (SAR) data. Near-global coverage and rapid repeat time (6–12 days) gives Sentinel-1 data the potential to be widely used for monitoring the Earth’s surface. Subtle land-cover and land surface changes can affect the phase and amplitude of the C-band SAR signal, and thus the coherence between two images collected before and after such changes. Analysis of SAR coherence therefore serves as a rapidly deployable and powerful tool to track both seasonal changes and rapid surface disturbances following natural disasters. An advantage of using Sentinel-1 C-band radar data is the ability to easily construct time series of coherence for a region of interest at low cost. In this paper, we propose a new method for Potentially Affected Area (PAA) detection following a natural hazard event. Based on the coherence time series, the proposed method (1) determines the natural variability of coherence within each pixel in the region of interest, accounting for factors such as seasonality and the inherent noise of variable surfaces; and (2) compares pixel-by-pixel syn-event coherence to temporal coherence distributions to determine where statistically significant coherence loss has occurred. The user can determine to what degree the syn-event coherence value (e.g., 1st, 5th percentile of pre-event distribution) constitutes a PAA, and integrate pertinent regional data, such as population density, to rank and prioritise PAAs. We apply the method to two case studies, Sarpol-e, Iran following the 2017 Iran-Iraq earthquake, and a landslide-prone region of NW Argentina, to demonstrate how rapid identification and interpretation of potentially affected areas can be performed shortly following a natural hazard event.
West Africa has been afflicted by droughts since the declining rains of the 1970s. Therefore, this study examines the characteristics of drought over the Niger River Basin (NRB), investigates the influence of the drought on the river flow, and projects the impacts of future climate change on drought. A combination of observation data and regional climate simulations of past (1986-2005) and future climates (2046-2065 and 2081-2100) were analyzed. The standardized precipitation index (SPI) and standardized precipitation and evapotranspiration index (SPEI) were used to characterize drought while the standardized runoff index (SRI) was used to quantify river flow. Results of the study show that the historical pattern of drought is consistent with previous studies over the Basin and most part of West Africa. RCA4 ensemble gives realistic simulations of the climatology of the Basin in the past climate. Generally, an increase in drought intensity and frequency are projected over NRB. The coupling between SRI and drought indices was very strong (P < 0.05). The dominant peaks can be classified into three distinct drought cycles with periods 1-2, 2-4, 4-8 years. These cycles may be associated with Quasi-Biennial Oscillation (QBO) and El-Nino Southern Oscillation (ENSO). River flow was highly sensitive to precipitation in the NRB and a 1-3 month lead time was found between drought indices and SRI. Under RCP4.5, changes in the SPEI drought frequency range from 1.8 (2046-2065) to 2.4 (2081-2100) month year(-1) while under RCP8.5, the change ranges from 2.2 (2046-2065) to 3.0 month year(-1) (2081-2100). Niger Middle sub-basin is likely to be mostly impacted in the future while the Upper Niger was projected to be least impacted. Results of this study may guide policymakers to evolve strategies to facilitate vulnerability assessment and adaptive capacity of the basin in order to minimize the negative impacts of climate change.
From north to south, denudation rates from cosmogenic nuclides are similar to 10 t km(-2) yr(-1) at the arid Pan de Aziicar site, similar to 20 t km(2) yr(-1) at the semi-arid site of Santa Gracia, -60 t km(-2) yr(-1) at the Mediterranean climate site of La Campana, and similar to 30 t km(-2) yr(-1) at the humid site of Nahuelbuta. A and B horizons increase in thickness and elemental depletion or enrichment increases from north (similar to 26 degrees S) to south (similar to 38 degrees S) in these horizons. Differences in the degree of chemical weathering, quantified by the chemical depletion fraction (CDF), are significant only between the arid and sparsely vegetated site and the other three sites. Differences in the CDF between the sites, and elemental depletion within the sites are sometimes smaller than the variations induced by the bedrock heterogeneity. Microbial abundances (bacteria and archaea) in saprolite substantially increase from the arid to the semi-arid sites. With this study, we provide a comprehensive dataset characterizing the Critical Zone geochemistry in the Chilean Coastal Cordillera. This dataset confirms climatic controls on weathering and denudation rates and provides prerequisites to quantify the role of biota in future studies.
The location and magnitude of Himalayan tectonic activity has been debated for decades, and several aspects remain unknown. For instance, the spatial distribution of crustal shortening that ultimately sustains Himalayan topography and the activity of major fault zones remain unknown at Ma timescales. In this study, we address the spatial deformation pattern in the data-scarce western Himalaya. We calculated catchment averaged, normalized river-steepness indices of non-glaciated drainage basins with tributary catchment areas between 5 and 200 km(2) (n = 2138). We analyzed the spatial distribution of the relative change of river steepness both along and across strike to gain information about the regional distribution of differential uplift pattern and relate this to the activity of distinctive fault segments. For our study area, we observe a positive correlation of averaged k(sn) values with long-term exhumation rates derived from previously published thermochronologic datasets combined with thermal modeling as well as with millennial timescale denudation rates based on cosmogenic nuclide dating. Our results indicate three tectono-geomorphic segments with distinctive landscape morphology, structural architecture, and fault geometry along the western Himalaya: Garhwal-Sutlej, Chamba, and Kashmir Himalaya (from east to west). Moreover, our data recognize distinctive fault segments showing varying thrust activity along strike of the Main Frontal Thrust, the Main Boundary Thrust, and in the vicinity of the steep topographic transition between the Lesser and Greater Himalaya. In this region, we relate out-of-sequence deformation along major basement thrust ramps, such as the Munsiari Thrust with deformation along a mid-crustal ramp along the basal decollement. We suggest that during the Quaternary, all major fault zones in the Western Himalaya experienced out-of-sequence faulting and have accommodated some portion of crustal shortening.
The use of chironomids (Diptera: Chironomidae) as indicators of changes in natural and climatic settings allows one to reconstruct paleoclimatic and paleoenvironmental conditions. Based on the results of lithological and chironomid analyses, the loss on ignition (LOI) values, as well as the results of radiocarbon dating (C-14 AMS), new data on the paleoclimate in the Karelian Isthmus, in particular the quantitative reconstruction of the mean July temperature (T (July), A degrees C) in the Late Pleistocene and Holocene and the evolution of the ecosystem of Medvedevskoe Lake have been obtained.
As one of the 195 member countries of the United Nations, Germany signed the Sendai Framework for Disaster Risk Reduction 2015-2030 (SFDRR). Among other targets, the SFDRR aims at reducing direct economic losses caused by natural hazards by 2030. The United Nations Office for Disaster Risk Reduction (UNISDR) has hence proposed a methodology for estimating direct economic losses per event and country, based on experiences from developing countries. Since its usability in industrialized countries is unknown, this study presents the first implementation and validation of this approach in Germany. The methodology was tested for the three costliest natural hazard types in Germany, i.e. floods, wind and hail storms, considering 12 case studies between 1984 and 2016. Although the event-specific input data requirements are restricted to the number of damaged or destroyed units per sector, incomplete event documentations did not allow a full validation of all sectors necessary to describe the total direct economic loss. New modules (cars, forestry, paved roads, housing contents and overall costs of urban infrastructure) were developed to better adapt this methodology to German conditions. Whereas the original UNISDR methodology both over-and underestimates the losses of the tested events by a wide margin, the adapted methodology is able to calculate losses accounting well for all event types except for flash floods. Hence, this approach serves as a good starting point for macro-scale loss estimations. By implementing this approach into damage and event documentation and reporting standards, a consistent monitoring of the SFDRR could be achieved.
Rock deformation at depths in the Earth’s crust is often localized in high temperature shear zones occurring at different scales in a variety of lithologies. The presence of material heterogeneities is known to trigger shear zone development, but the mechanisms controlling initiation and evolution of localization are not fully understood. To investigate the effect of loading conditions on shear zone nucleation along heterogeneities, we performed torsion experiments under constant twist rate (CTR) and constant torque (CT) conditions in a Paterson-type deformation apparatus. The sample assemblage consisted of cylindrical Carrara marble specimens containing a thin plate of Solnhofen limestone perpendicular to the cylinder’s longitudinal axis. Under experimental conditions (900 °C, 400 MPa confining pressure), samples were plastically deformed and limestone is about 9 times weaker than marble, acting as a weak inclusion in a strong matrix. CTR experiments were performed at maximum bulk shear strain rates of ~ 2*10-4s-1, yielding peak shear stresses of ~ 20 MPa. CT tests were conducted at shear stresses of ~ 20 MPa resulting in bulk shear strain rates of 1-4*10-4s-1. Experiments were terminated at maximum bulk shear strains of ~ 0.3 and 1.0.Strain was localized within the Carrara marble in front of the inclusion in an area of strongly deformed grains and intense grain size reduction. Locally, evidences for coexisting brittle deformation are also observed regardless of the imposed loading conditions. The local shear strain at the inclusion tipis up to 30 times higher than the strain in the adjacent host rock, rapidly dropping to 5times higher at larger distance from the inclusion. At both investigated bulk strains, the evolution of microstructural and textural parameters is independent of loading conditions. Ourresults suggest that loading conditions do not significantly affect material heterogeneity-induced strain localization during its nucleation and transient stages.
Volcanic activities triggered or inhibited by resonance of volcanic edifices to large earthquakes
(2018)
The existence of a causal link between large earthquakes and volcanic unrest is widely accepted. Recent observations have also revealed counterintuitive negative responses of volcanoes to large earthquakes, including decreased gas emissions and subsidence in volcanic areas. In order to explore the mechanisms that could simultaneously explain both the positive and negative responses of volcanic activity to earthquakes, we here focus on the role played by topography. In the laboratory, we shook a volcanic edifice analogue, made of gel, previously injected with a buoyant fluid. We find that shaking triggers rapid migration of the buoyant fluid upward, downward, or laterally, depending on the fluid’s buoyancy and storage depth; bubbly fluids stored at shallow depth ascend, while low-buoyancy fluids descend or migrate laterally. The migration of fluids induced by shaking is two orders of magnitude faster than without shaking. Downward or lateral fluid migration may decrease volcanic gas emissions and cause subsidence as a negative response, while upward migration is consistent both with an increase in volcanic activity and immediate unrest (deformation and seismicity) after large earthquakes. The fluid migration is more efficient when the oscillation frequency is close to the resonance frequency of the edifice. The resonance frequency for a 30-km-wide volcanic mountain range, such as those where subsidence was observed, is ∼0.07 Hz. Only large earthquakes are able to cause oscillation at such low frequencies.
Many regions central to our understanding of tectonics and landscape evolution are active or ancient magmatic terranes, and robust interpretation of low-temperature thermochronologic ages in these settings requires careful attention to the drivers of rock heating and cooling, including magmatism. However, we currently lack a quantitative framework for evaluating the potential role of magmatic coolingthat is, post-magmatic thermal relaxationin shaping cooling age patterns in regions with a history of intrusive magmatism. Here we use analytical approximations and numerical models to characterize how low-temperature thermochronometers document cooling inside and around plutons in steadily exhuming environments. Our models predict that the thermal field a pluton intrudes into, specifically the ambient temperatures relative to the closure temperature of a given thermochronometer, is as important as the pluton size and temperature in controlling the pattern and extent of thermochronometer resetting in the country rocks around a pluton. We identify one advective and several conductive timescales that govern the relationship between the crystallization and cooling ages inside a pluton. In synthetic vertical age-elevation relationships (AERs), resetting next to plutons results in changes in AER slope that could be misinterpreted as past changes in exhumation rate if the history of magmatism is not accounted for. Finally, we find that large midcrustal plutons, such as those emplaced at similar to 10-15-km depth, can reset the low-temperature thermochronometers far above them in the upper crusta result with considerable consequences for thermochronology in arcs and regions with a history of magmatic activity that may not have a surface expression.
For attributing hydrological changes to anthropogenic climate change, catchment models are driven by climate model output. A widespread approach to bridge the spatial gap between global climate and hydrological catchment models is to use a weather generator conditioned on weather patterns (WPs). This approach assumes that changes in local climate are characterized by between-type changes of patterns. In this study we test this assumption by analyzing a previously developed WP classification for the Rhine basin, which is based on dynamic and thermodynamic variables. We quantify changes in pattern characteristics and associated climatic properties. The amount of between- and within-type changes is investigated by comparing observed trends to trends resulting solely from WP occurrence. To overcome uncertainties in trend detection resulting from the selected time period, all possible periods in 1901-2010 with a minimum length of 31 years are analyzed. Increasing frequency is found for some patterns associated with high precipitation, although the trend sign highly depends on the considered period. Trends and interannual variations of WP frequencies are related to the long-term variability of large-scale circulation modes. Long-term WP internal warming is evident for summer patterns and enhanced warming for spring/autumn patterns since the 1970s. Observed trends in temperature and partly in precipitation are mainly associated with frequency changes of specific WPs, but some amount of within-type changes remains. The classification can be used for downscaling of past changes considering this limitation, but the inclusion of thermodynamic variables into the classification impedes the downscaling of future climate projections.
For over a century, geomorphologists have attempted to unravel information about landscape evolution, and processes that drive it, using river profiles. Many studies have combined new topographic datasets with theoretical models of channel incision to infer erosion rates, identify rock types with different resistance to erosion, and detect potential regions of tectonic activity. The most common metric used to analyse river profile geometry is channel steepness, or k(s). However, the calculation of channel steepness requires the normalisation of channel gradient by drainage area. This normalisation requires a power law exponent that is referred to as the channel concavity index. Despite the concavity index being crucial in determining channel steepness, it is challenging to constrain. In this contribution, we compare both slope-area methods for calculating the concavity index and methods based on integrating drainage area along the length of the channel, using so-called "chi" (chi) analysis. We present a new chi-based method which directly compares chi values of tributary nodes to those on the main stem; this method allows us to constrain the concavity index in transient landscapes without assuming a linear relationship between chi and elevation. Patterns of the concavity index have been linked to the ratio of the area and slope exponents of the stream power incision model (m/n); we therefore construct simple numerical models obeying detachment-limited stream power and test the different methods against simulations with imposed m and n. We find that chi-based methods are better than slope-area methods at reproducing imposed m/n ratios when our numerical landscapes are subject to either transient uplift or spatially varying uplift and fluvial erodibility. We also test our methods on several real landscapes, including sites with both lithological and structural heterogeneity, to provide examples of the methods' performance and limitations. These methods are made available in a new software package so that other workers can explore how the concavity index varies across diverse landscapes, with the aim to improve our understanding of the physics behind bedrock channel incision.
The impressive number of stream gauges in Chile, combined with a suite of past and recent large earthquakes, makes Chile a unique natural laboratory to study several streams that recorded responses to multiple seismic events. We document changes in discharge in eight streams in Chile following two or more large earthquakes. In all cases, discharge increases. Changes in discharge occur for peak ground velocities greater than about 7-11cm/s. Above that threshold, the magnitude of both the increase in discharge and the total excess water do not increase with increasing peak ground velocities. While these observations are consistent with previous work in California, they conflict with lab experiments that show that the magnitude of permeability changes increases with increasing amplitude of ground motion. Instead, our study suggests that streamflow responses are binary. Plain Language Summary Earthquakes deform and shake the surface and the ground below. These changes may affect groundwater flows by increasing the permeability along newly formed cracks and/or clearing clogged pores. As a result, groundwater flow may substantially increase after earthquakes and remain elevated for several months. Here we document streamflow anomalies following multiple high magnitude earthquakes in multiple streams in one of the most earthquake prone regions worldwide, Chile. We take advantage of the dense monitoring network in Chile that recorded streamflow since the 1940s. We show that once a critical ground motion is exceeded, streamflow responses to earthquakes can be expected.