Refine
Year of publication
- 2018 (314) (remove)
Document Type
- Article (237)
- Doctoral Thesis (25)
- Postprint (20)
- Other (17)
- Review (10)
- Habilitation Thesis (2)
- Conference Proceeding (1)
- Master's Thesis (1)
- Report (1)
Keywords
- climate change (8)
- Himalaya (7)
- permafrost (7)
- Germany (5)
- landslides (4)
- Climate change (3)
- Kalahari (3)
- Naturgefahren (3)
- Palynology (3)
- Permafrost (3)
Institute
- Institut für Geowissenschaften (314) (remove)
Glacial deposits on the high-altitude, arid southern Central Andean Plateau (CAP), the Puna in northwestern Argentina, document past changes in climate, but the associated geomorphic features have rarely been directly dated. This study provides direct age control of glacial moraine deposits from the central Puna (24 degrees S) at elevations of 3900-5000 m through surface exposure dating with cosmogenic nuclides. Our results show that the most extensive glaciations occurred before 95 ka and an additional major advance occurred between 46 and 39 ka. The latter period is synchronous with the highest lake levels in the nearby Pozuelos basin and the Minchin (Inca Huasi) wet phase on the Altiplano in the northern CAP. None of the dated moraines produced boulder ages corresponding to the Tauca wet phase (24-15 ka). Additionally, the volcanic lithologies of the deposits allow us to establish production ratios at low latitude and high elevation for five different nuclide and mineral systems: Be-10, Ne-21, and Al-26 from quartz (11 or 12 samples) and He-3 and Ne-21 from pyroxene (10 samples). We present production ratios for all combinations of the measured nuclides and cross-calibrated production rates for 21Ne in pyroxene and quartz for the high, (sub-)tropical Andes. The production rates are based on our Be-10-normalized production ratios and a weighted mean of reference 10Be production rates calibrated in the high, tropical Andes (4.02 +/- 0.12 at g(-1) yr(-1)). These are, Ne-21(qtz): 18.1 +/- 1.2 at g(-1) yr(-1) and Ne-21(px): 36.6 +/- 1.8 at g(-1) yr(-1) (En(88-94)) scaled to sea level and high latitude using the Lal/Stone scheme, with 1 sigma uncertainties. As He-3 and Al-26 have been directly calibrated in the tropical Andes, we recommend using those rates. Finally, we compare exposure ages calculated using all measured cosmogenic nuclides from each sample, including 11 feldspar samples measured for Cl-36, and a suite of previously published production rates. (C) 2018 Published by Elsevier B.V.
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.
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.
Mechanical and/or chemical removal of material from the subsurface may generate large subsurface cavities, the destabilisation of which can lead to ground collapse and the formation of sinkholes. Numerical simulation of the interaction of cavity growth, host material deformation and overburden collapse is desirable to better understand the sinkhole hazard but is a challenging task due to the involved high strains and material discontinuities. Here, we present 2-D distinct element method numerical simulations of cavity growth and sinkhole development. Firstly, we simulate cavity formation by quasi-static, stepwise removal of material in a single growing zone of an arbitrary geometry and depth. We benchmark this approach against analytical and boundary element method models of a deep void space in a linear elastic material. Secondly, we explore the effects of properties of different uniform materials on cavity stability and sinkhole development. We perform simulated biaxial tests to calibrate macroscopic geotechnical parameters of three model materials representative of those in which sinkholes develop at the Dead Sea shoreline: mud, alluvium and salt. We show that weak materials do not support large cavities, leading to gradual sagging or suffusion-style subsidence. Strong materials support quasi-stable to stable cavities, the overburdens of which may fail suddenly in a caprock or bedrock collapse style. Thirdly, we examine the consequences of layered arrangements of weak and strong materials. We find that these are more susceptible to sinkhole collapse than uniform materials not only due to a lower integrated strength of the overburden but also due to an inhibition of stabilising stress arching. Finally, we compare our model sinkhole geometries to observations at the Ghor Al-Haditha sinkhole site in Jordan. Sinkhole depth ∕ diameter ratios of 0.15 in mud, 0.37 in alluvium and 0.33 in salt are reproduced successfully in the calibrated model materials. The model results suggest that the observed distribution of sinkhole depth ∕ diameter values in each material type may partly reflect sinkhole growth trends.
Mechanical and/or chemical removal of material from the subsurface may generate large subsurface cavities, the destabilisation of which can lead to ground collapse and the formation of sinkholes. Numerical simulation of the interaction of cavity growth, host material deformation and overburden collapse is desirable to better understand the sinkhole hazard but is a challenging task due to the involved high strains and material discontinuities. Here, we present 2-D distinct element method numerical simulations of cavity growth and sinkhole development. Firstly, we simulate cavity formation by quasi-static, stepwise removal of material in a single growing zone of an arbitrary geometry and depth. We benchmark this approach against analytical and boundary element method models of a deep void space in a linear elastic material. Secondly, we explore the effects of properties of different uniform materials on cavity stability and sinkhole development. We perform simulated biaxial tests to calibrate macroscopic geotechnical parameters of three model materials representative of those in which sinkholes develop at the Dead Sea shoreline: mud, alluvium and salt. We show that weak materials do not support large cavities, leading to gradual sagging or suffusion-style subsidence. Strong materials support quasi-stable to stable cavities, the overburdens of which may fail suddenly in a caprock or bedrock collapse style. Thirdly, we examine the consequences of layered arrangements of weak and strong materials. We find that these are more susceptible to sinkhole collapse than uniform materials not only due to a lower integrated strength of the overburden but also due to an inhibition of stabilising stress arching. Finally, we compare our model sinkhole geometries to observations at the Ghor Al-Haditha sinkhole site in Jordan. Sinkhole depth / diameter ratios of 0.15 in mud, 0.37 in alluvium and 0.33 in salt are reproduced successfully in the calibrated model materials. The model results suggest that the observed distribution of sinkhole depth / diameter values in each material type may partly reflect sinkhole growth trends.
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.
From precipitation to damage
(2018)
Flood risk assessments for large river basins often involve piecing together smaller-scale assessments leading to erroneous risk statements. We describe a coupled model chain for quantifying flood risk at the scale of 100,000 km(2). It consists of a catchment model, a 1D-2D river network model, and a loss model. We introduce the model chain and present two applications. The first application for the Elbe River basin with an area of 66,000 km(2) demonstrates that it is feasible to simulate the complete risk chain for large river basins in a continuous simulation mode with high temporal and spatial resolution. In the second application, RFM is coupled to a multisite weather generator and applied to the Mulde catchment with an area of 6,000 km(2). This approach is able to provide a very long time series of spatially heterogeneous patterns of precipitation, discharge, inundation, and damage. These patterns respect the spatial correlation of the different processes and are suitable to derive large-scale risk estimates. We discuss how the RFM approach can be transferred to the continental scale.
The correct estimation of site-specific attenuation is crucial for the assessment of seismic hazard. Downhole instruments provide in this context valuable information to constrain attenuation directly from data. In this study, we apply an interferometric approach to this problem by deconvolving seismic motions recorded at depth with those recorded at the surface. In doing so, incident and surface-reflected waves can be separated. We apply this technique not only to earthquake data but also to recordings of ambient vibrations. We compute the transfer function between incident and surface-reflected waves in order to infer frequency-dependent quality factors for S waves. The method is applied to a 87m deep borehole sensor and a colocated surface instrument situated at a hard-rock site in West Bohemia/Vogtland, Germany. We show that the described method provides comparable attenuation estimates using either earthquake data or ambient noise for frequencies between 5 and 15 Hz. Moreover, a single hour of noise recordings seems to be sufficient to yield stable deconvolution traces and quality factors, thus, offering a fast and easy way to derive attenuation estimates from borehole recordings even in low- to mid-seismicity regions.
Foreword
(2018)
The flash-flood in Braunsbach in the north-eastern part of Baden-Wuerttemberg/Germany was a particularly strong and concise event which took place during the floods in southern Germany at the end of May/early June 2016. This article presents a detailed analysis of the hydro-meteorological forcing and the hydrological consequences of this event. A specific approach, the "forensic hydrological analysis" was followed in order to include and combine retrospectively a variety of data from different disciplines. Such an approach investigates the origins, mechanisms and course of such natural events if possible in a "near real time" mode, in order to follow the most recent traces of the event. The results show that it was a very rare rainfall event with extreme intensities which, in combination with catchment properties, led to extreme runoff plus severe geomorphological hazards, i.e. great debris flows, which together resulted in immense damage in this small rural town Braunsbach. It was definitely a record-breaking event and greatly exceeded existing design guidelines for extreme flood discharge for this region, i.e. by a factor of about 10. Being such a rare or even unique event, it is not reliably feasible to put it into a crisp probabilistic context. However, one can conclude that a return period clearly above 100 years can be assigned for all event components: rainfall, peak discharge and sediment transport. Due to the complex and interacting processes, no single flood cause or reason for the very high damage can be identified, since only the interplay and the cascading characteristics of those led to such an event. The roles of different human activities on the origin and/or intensification of such an extreme event are finally discussed. (C) 2018 Elsevier B.V. All rights reserved.
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.
One of the main purposes of detrital thermochronology is to provide constraints on the regional-scale exhumation rate and its spatial variability in actively eroding mountain ranges. Procedures that use cooling age distributions coupled with hypsometry and thermal models have been developed in order to extract quantitative estimates of erosion rate and its spatial distribution, assuming steady state between tectonic uplift and erosion. This hypothesis precludes the use of these procedures to assess the likely transient response of mountain belts to changes in tectonic or climatic forcing. Other methods are based on an a priori knowledge of the in situ distribution of ages to interpret the detrital age distributions. In this paper, we describe a simple method that, using the observed detrital mineral age distributions collected along a river, allows us to extract information about the relative distribution of erosion rates in an eroding catchment without relying on a steady-state assumption, the value of thermal parameters or an a priori knowledge of in situ age distributions. The model is based on a relatively low number of parameters describing lithological variability among the various sub-catchments and their sizes and only uses the raw ages. The method we propose is tested against synthetic age distributions to demonstrate its accuracy and the optimum conditions for it use. In order to illustrate the method, we invert age distributions collected along the main trunk of the Tsangpo-Siang-Brahmaputra river system in the eastern Himalaya. From the inversion of the cooling age distributions we predict present-day erosion rates of the catchments along the Tsangpo-Siang-Brahmaputra river system, as well as some of its tributaries. We show that detrital age distributions contain dual information about present-day erosion rate, i. e., from the predicted distribution of surface ages within each catchment and from the relative contribution of any given catchment to the river distribution. The method additionally allows comparing modern erosion rates to long-term exhumation rates. We provide a simple implementation of the method in Python code within a Jupyter Notebook that includes the data used in this paper for illustration purposes.
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.
Arctic warming has implications for the functioning of terrestrial Arctic ecosystems, global climate and socioeconomic systems of northern communities. A research gap exists in high spatial resolution monitoring and understanding of the seasonality of permafrost degradation, spring snowmelt and vegetation phenology. This thesis explores the diversity and utility of dense TerraSAR-X (TSX) X-Band time series for monitoring ice-rich riverbank erosion, snowmelt, and phenology of Arctic vegetation at long-term study sites in the central Lena Delta, Russia and on Qikiqtaruk (Herschel Island), Canada. In the thesis the following three research questions are addressed:
• Is TSX time series capable of monitoring the dynamics of rapid permafrost degradation in ice-rich permafrost on an intra-seasonal scale and can these datasets in combination with climate data identify the climatic drivers of permafrost degradation?
• Can multi-pass and multi-polarized TSX time series adequately monitor seasonal snow cover and snowmelt in small Arctic catchments and how does it perform compared to optical satellite data and field-based measurements?
• Do TSX time series reflect the phenology of Arctic vegetation and how does the recorded signal compare to in-situ greenness data from RGB time-lapse camera data and vegetation height from field surveys?
To answer the research questions three years of TSX backscatter data from 2013 to 2015 for the Lena Delta study site and from 2015 to 2017 for the Qikiqtaruk study site were used in quantitative and qualitative analysis complimentary with optical satellite data and in-situ time-lapse imagery.
The dynamics of intra-seasonal ice-rich riverbank erosion in the central Lena Delta, Russia were quantified using TSX backscatter data at 2.4 m spatial resolution in HH polarization and validated with 0.5 m spatial resolution optical satellite data and field-based time-lapse camera data. Cliff top lines were automatically extracted from TSX intensity images using threshold-based segmentation and vectorization and combined in a geoinformation system with manually digitized cliff top lines from the optical satellite data and rates of erosion extracted from time-lapse cameras. The results suggest that the cliff top eroded at a constant rate throughout the entire erosional season. Linear mixed models confirmed that erosion was coupled with air temperature and precipitation at an annual scale, seasonal fluctuations did not influence 22-day erosion rates. The results highlight the potential of HH polarized X-Band backscatter data for high temporal resolution monitoring of rapid permafrost degradation.
The distinct signature of wet snow in backscatter intensity images of TSX data was exploited to generate wet snow cover extent (SCE) maps on Qikiqtaruk at high temporal resolution. TSX SCE showed high similarity to Landsat 8-derived SCE when using cross-polarized VH data. Fractional snow cover (FSC) time series were extracted from TSX and optical SCE and compared to FSC estimations from in-situ time-lapse imagery. The TSX products showed strong agreement with the in-situ data and significantly improved the temporal resolution compared to the Landsat 8 time series. The final combined FSC time series revealed two topography-dependent snowmelt patterns that corresponded to in-situ measurements. Additionally TSX was able to detect snow patches longer in the season than Landsat 8, underlining the advantage of TSX for detection of old snow. The TSX-derived snow information provided valuable insights into snowmelt dynamics on Qikiqtaruk previously not available.
The sensitivity of TSX to vegetation structure associated with phenological changes was explored on Qikiqtaruk. Backscatter and coherence time series were compared to greenness data extracted from in-situ digital time-lapse cameras and detailed vegetation parameters on 30 areas of interest. Supporting previous results, vegetation height corresponded to backscatter intensity in co-polarized HH/VV at an incidence angle of 31°. The dry, tall shrub dominated ecological class showed increasing backscatter with increasing greenness when using the cross polarized VH/HH channel at 32° incidence angle. This is likely driven by volume scattering of emerging and expanding leaves. Ecological classes with more prostrate vegetation and higher bare ground contributions showed decreasing backscatter trends over the growing season in the co-polarized VV/HH channels likely a result of surface drying instead of a vegetation structure signal. The results from shrub dominated areas are promising and provide a complementary data source for high temporal monitoring of vegetation phenology.
Overall this thesis demonstrates that dense time series of TSX with optical remote sensing and in-situ time-lapse data are complementary and can be used to monitor rapid and seasonal processes in Arctic landscapes at high spatial and temporal resolution.
Silicon (Si) is considered as a quasiessential element for higher plants as its uptake increases plant growth and resistance against abiotic as well as biotic stresses. Foliar application of fertilizers generally is assumed to be a comparably environment-friendly form of fertilization because only small quantities are needed. The interest in foliar fertilization and the use of Si as a fertilizer in general increased significantly within the last decades, but there are only few publications dealing with the foliar application of Si at all. In the present review, the effects of Si foliar fertilization, including nano-Si fertilizers, on the three most important crops on a global scale, that is, maize, rice, and wheat, are summarized. Additionally, different pathways (i.e., cuticular pathways, stomata, and trichomes) of foliar uptake and functioning of Si foliar fertilizers against biotic (i.e., fungal diseases and harmful insects), as well as abiotic (i.e., water stress, macronutrient imbalance, and heavy metal toxicity) stressors are discussed. Future research should especially focus on (1) the gathering of empirical data from field and greenhouse experiments, (2) the intensification of co-operations between practitioners and scientists, (3) interdisciplinary research, and (4) the analysis of results from multiple studies (meta-analysis, big data) to fully understand effects, uptake, and functioning of Si foliar fertilizers and to evaluate their potential in modern sustainable agriculture concepts.
Background. Metal recycling factories (MRFs) have developed rapidly in Nigeria as recycling policies have been increasingly embraced. These MRFs are point sources for introducing potentially toxic elements (PTEs) into environmental media. Objectives. The aim of this study was to determine the constituents (elemental and mineralogy) of the wastes (slag and particulate matter, (PM)) and soils around the MRFs and to determine the level of pollution within the area. Methods. Sixty samples (30 slag samples, 15 soil samples and 15 PM samples) were collected for this study. The soils, slag and PM samples were analyzed for elemental constituents using inductively coupled plasma optical emission spectrometry. Mineralogy of the PM was determined using scanning electron microscope-energy dispersive spectroscopy (SEM-EDS), and soil mineralogy was determined by an X-ray diffractometer (XRD). Results. The results of the soil analyses revealed the following concentrations for the selected metals in mg/kg include lead (Pb) (21.0-2399.0), zinc (Zn) (56.0-4188.0), copper (Cu) (10.0-1470.0), nickel (Ni) (6.0-215.0), chromium (Cr) (921.0-1737.0) and cadmium (Cd) (below detectable limit (Bdl)-18.1). For the slags the results were Pb (68.0-.333.0), Zn (1364.0-3062), Cu (119.0-1470.0), Ni (12.0-675.0), Cr (297-1737) and Cd (Bdl-15.8). The results in µg/g for the metal analysis in PM were Pb (4.6-160.0), Zn (18.0-471.0), Cu (2.5-11.0), Ni (0.8-4.2), and Cr (2.5-11.0), while Cd was undetected. The slags are currently utilized for filling the foundations of buildings and roads, providing additional pathways for the introduction of PTEs into the environment from the suspended materials generated from mechanical breakdown of the slags. Conclusions. The MRFs were found to have impacted the quality of environmental media through the introduction of PTEs, impairing soil quality, in addition to PM, which can have detrimental health consequences. Further studies on the health implications of these pollutants and their impacts on human health are needed. Competing Interests. The authors declare no competing financial interests
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.
Evaluating climate geoengineering proposals in the context of the Paris Agreement temperature goals
(2018)
Current mitigation efforts and existing future commitments are inadequate to accomplish the Paris Agreement temperature goals. In light of this, research and debate are intensifying on the possibilities of additionally employing proposed climate geoengineering technologies, either through atmospheric carbon dioxide removal or farther-reaching interventions altering the Earth’s radiative energy budget. Although research indicates that several techniques may eventually have the physical potential to contribute to limiting climate change, all are in early stages of development, involve substantial uncertainties and risks, and raise ethical and governance dilemmas. Based on present knowledge, climate geoengineering techniques cannot be relied on to significantly contribute to meeting the Paris Agreement temperature goals.
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.
We explore the potential of spaceborne radar (SR) observations from the Ku-band precipitation radars onboard the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) satellites as a reference to quantify the ground radar (GR) reflectivity bias. To this end, the 3-D volume-matching algorithm proposed by Schwaller and Morris (2011) is implemented and applied to 5 years (2012–2016) of observations. We further extend the procedure by a framework to take into account the data quality of each ground radar bin. Through these methods, we are able to assign a quality index to each matching SR–GR volume, and thus compute the GR calibration bias as a quality-weighted average of reflectivity differences in any sample of matching GR–SR volumes. We exemplify the idea of quality-weighted averaging by using the beam blockage fraction as the basis of a quality index. As a result, we can increase the consistency of SR and GR observations, and thus the precision of calibration bias estimates. The remaining scatter between GR and SR reflectivity as well as the variability of bias estimates between overpass events indicate, however, that other error sources are not yet fully addressed. Still, our study provides a framework to introduce any other quality variables that are considered relevant in a specific context. The code that implements our analysis is based on the wradlib open-source software library, and is, together with the data, publicly available to monitor radar calibration or to scrutinize long series of archived radar data back to December 1997, when TRMM became operational.
We explore the potential of spaceborne radar (SR) observations from the Ku-band precipitation radars onboard the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) satellites as a reference to quantify the ground radar (GR) reflectivity bias. To this end, the 3-D volume-matching algorithm proposed by Schwaller and Morris (2011) is implemented and applied to 5 years (2012–2016) of observations. We further extend the procedure by a framework to take into account the data quality of each ground radar bin. Through these methods, we are able to assign a quality index to each matching SR–GR volume, and thus compute the GR calibration bias as a quality-weighted average of reflectivity differences in any sample of matching GR–SR volumes. We exemplify the idea of quality-weighted averaging by using the beam blockage fraction as the basis of a quality index. As a result, we can increase the consistency of SR and GR observations, and thus the precision of calibration bias estimates. The remaining scatter between GR and SR reflectivity as well as the variability of bias estimates between overpass events indicate, however, that other error sources are not yet fully addressed. Still, our study provides a framework to introduce any other quality variables that are considered relevant in a specific context. The code that implements our analysis is based on the wradlib open-source software library, and is, together with the data, publicly available to monitor radar calibration or to scrutinize long series of archived radar data back to December 1997, when TRMM became operational.
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.
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.
In catastrophe risk modeling, a defensible estimation of impact severity and its likelihood of occurrence to a portfolio of assets can only be made through a rigorous treatment of uncertainty and the consideration of multiple alternative models. This approach, however, requires repeating lengthy calculations multiple times. To limit the demand on computational time and resources, a frequent practice in the industry is to estimate the distribution of earthquake-induced portfolio losses using a simulated catalog of events from a single representative mean ground motion hazard model for the region. This simplified approach is faster but may provide biased estimates of the likelihood of occurrence of the large and infrequent losses that drive many risk mitigation decisions. Investigation through case studies of different portfolios of assets located in the San Francisco Bay Region shows the potential for both a bias in the mean loss estimates and an underestimation of their central 70% interpercentile. We propose a simplified and computationally practical approach that reduces the bias in the mean portfolio loss estimates. This approach does not improve the estimate of the interpercentile range, however, a quantity of no direct practical use.
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.
To safeguard the sustainable use of ecosystems and their services, early detection of potentially damaging changes in functional capabilities is needed. To support a proper ecosystem management, the analysis of an ecosystem’s vulnerability provide information on its weaknesses as well as on its capacity to recover after suffering an impact. However, the application of the vulnerability concept to ecosystems is still an emerging topic. After providing background on the vulnerability concept, we summarize existing ecosystem vulnerability research on the basis of a systematic literature review with a special focus on ecosystem type, disciplinary background, and more detailed definition of the ecosystem vulnerability components. Using the Web of ScienceTM Core Collection, we overviewed the literature from 1991 onwards but used the 5 years from 2011 to 2015 for an in-depth analysis, including 129 articles. We found that ecosystem vulnerability analysis has been applied most notably in conservation biology, climate change research, and ecological risk assessments, pinpointing a limited spreading across the environmental sciences. It occurred primarily within marine and freshwater ecosystems. To avoid confusion, we recommend using the unambiguous term ecosystem vulnerability rather than ecological, environmental, population, or community vulnerability. Further, common ground has been identified, on which to define the ecosystem vulnerability components exposure, sensitivity, and adaptive capacity. We propose a framework for ecosystem assessments that coherently connects the concepts of vulnerability, resilience, and adaptability as different ecosystem responses. A short outlook on the possible operationalization of the concept by ecosystem vulnerabilty indices, and a conclusion section complete the review.
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.
A potential human footprint on Western Central African rainforests before the Common Era has become the focus of an ongoing controversy. Between 3,000 y ago and 2,000 y ago, regional pollen sequences indicate a replacement of mature rainforests by a forest-savannah mosaic including pioneer trees. Although some studies suggested an anthropogenic influence on this forest fragmentation, current interpretations based on pollen data attribute the "rainforest crisis" to climate change toward a drier, more seasonal climate. A rigorous test of this hypothesis, however, requires climate proxies independent of vegetation changes. Here we resolve this controversy through a continuous 10,500-y record of both vegetation and hydrological changes from Lake Barombi in Southwest Cameroon based on changes in carbon and hydrogen isotope compositions of plant waxes. delta C-13-inferred vegetation changes confirm a prominent and abrupt appearance of C-4 plants in the Lake Barombi catchment, at 2,600 calendar years before AD 1950 (cal y BP), followed by an equally sudden return to rainforest vegetation at 2,020 cal y BP. delta D values from the same plant wax compounds, however, show no simultaneous hydrological change. Based on the combination of these data with a comprehensive regional archaeological database we provide evidence that humans triggered the rainforest fragmentation 2,600 y ago. Our findings suggest that technological developments, including agricultural practices and iron metallurgy, possibly related to the large-scale Bantu expansion, significantly impacted the ecosystems before the Common Era.
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.
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.
Participants of the 2017 European Space Weather Week in Ostend, Belgium, discussed the stakeholder requirements for space weather-related models. It was emphasized that stakeholders show an increased interest in space weather-related models. Participants of the meeting discussed particular prediction indicators that can provide first-order estimates of the impact of space weather on engineering systems.
Crop yield variations are strongly influenced by the spatial and temporal availabilities of water and nitrogen in the soil during the crop growth season. To estimate the quantities and distributions of water and nitrogen within a given soil, process-oriented soil models have often been used. These models require detailed information about the soil characteristics and profile architecture (e.g., soil depth, clay content, bulk density, field capacity and wilting point), but high resolution information about these soil properties, both vertically and laterally, is difficult to obtain through conventional approaches. However, on-the-go electrical resistivity tomography (ERT) measurements of the soil and data inversion tools have recently improved the lateral resolutions of the vertically distributed measurable information. Using these techniques, nearly 19,000 virtual soil profiles with defined layer depths were successfully created for a 30 ha silty cropped soil over loamy and sandy substrates in Central Germany, which were used to initialise the CArbon and Nitrogen DYnamics (CANDY) model. The soil clay content was derived from the electrical resistivity (ER) and the collected soil samples using a simple linear regression approach (the mean R-2 of clay = 0.39). The additional required structural and hydrological properties were derived from pedotransfer functions. The modelling results, derived soil texture distributions and original ER data were compared with the spatial winter wheat yield distribution in a relatively dry year using regression and boundary line analysis. The yield variation was best explained by the simulated soil water content (R-2 = 0.18) during the grain filling and was additionally validated by the measured soil water content with a root mean square error (RMSE) of 7.5 Vol%.
The Gutenberg-Richter relation for earthquake magnitudes is the most famous empirical law in seismology. It states that the frequency of earthquake magnitudes follows an exponential distribution; this has been found to be a robust feature of seismicity above the completeness magnitude, and it is independent of whether global, regional, or local seismicity is analyzed. However, the exponent b of the distribution varies significantly in space and time, which is important for process understanding and seismic hazard assessment; this is particularly true because of the fact that the Gutenberg-Richter b-value acts as a proxy for the stress state and quantifies the ratio of large-to-small earthquakes. In our work, we focus on the automatic detection of statistically significant temporal changes of the b-value in seismicity data. In our approach, we use Bayes factors for model selection and estimate multiple change-points of the frequency-magnitude distribution in time. The method is first applied to synthetic data, showing its capability to detect change-points as function of the size of the sample and the b-value contrast. Finally, we apply this approach to examples of observational data sets for which b-value changes have previously been stated. Our analysis of foreshock and after-shock sequences related to mainshocks, as well as earthquake swarms, shows that only a portion of the b-value changes is statistically significant.
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.
Time series of groundwater and stream water quality often exhibit substantial temporal and spatial variability, whereas typical existing monitoring data sets, e.g. from environmental agencies, are usually characterized by relatively low sampling frequency and irregular sampling in space and/or time. This complicates the differentiation between anthropogenic influence and natural variability as well as the detection of changes in water quality which indicate changes in single drivers. We suggest the new term "dominant changes" for changes in multivariate water quality data which concern (1) multiple variables, (2) multiple sites and (3) long-term patterns and present an exploratory framework for the detection of such dominant changes in data sets with irregular sampling in space and time. Firstly, a non-linear dimension-reduction technique was used to summarize the dominant spatiotemporal dynamics in the multivariate water quality data set in a few components. Those were used to derive hypotheses on the dominant drivers influencing water quality. Secondly, different sampling sites were compared with respect to median component values. Thirdly, time series of the components at single sites were analysed for long-term patterns. We tested the approach with a joint stream water and groundwater data set quality consisting of 1572 samples, each comprising sixteen variables, sampled with a spatially and temporally irregular sampling scheme at 29 sites in northeast Germany from 1998 to 2009. The first four components were interpreted as (1) an agriculturally induced enhancement of the natural background level of solute concentration, (2) a redox sequence from reducing conditions in deep groundwater to post-oxic conditions in shallow groundwater and oxic conditions in stream water, (3) a mixing ratio of deep and shallow groundwater to the streamflow and (4) sporadic events of slurry application in the agricultural practice. Dominant changes were observed for the first two components. The changing intensity of the first component was interpreted as response to the temporal variability of the thickness of the unsaturated zone. A steady increase in the second component at most stream water sites pointed towards progressing depletion of the denitrification capacity of the deep aquifer.
More than 41% of the Earth’s land area is covered by permanent or seasonally arid dryland ecosystems. Global development and human activity have led to an increase in aridity, resulting in ecosystem degradation and desertification around the world. The objective of the present work was to investigate and compare the microbial community structure and geochemical characteristics of two geographically distinct saline pan sediments in the Kalahari Desert of southern Africa. Our data suggest that these microbial communities have been shaped by geochemical drivers, including water content, salinity, and the supply of organic matter. Using Illumina 16S rRNA gene sequencing, this study provides new insights into the diversity of bacteria and archaea in semi-arid, saline, and low-carbon environments. Many of the observed taxa are halophilic and adapted to water-limiting conditions. The analysis reveals a high relative abundance of halophilic archaea (primarily Halobacteria), and the bacterial diversity is marked by an abundance of Gemmatimonadetes and spore-forming Firmicutes. In the deeper, anoxic layers, candidate division MSBL1, and acetogenic bacteria (Acetothermia) are abundant. Together, the taxonomic information and geochemical data suggest that acetogenesis could be a prevalent form of metabolism in the deep layers of a saline pan.
More than 41% of the Earth’s land area is covered by permanent or seasonally arid dryland ecosystems. Global development and human activity have led to an increase in aridity, resulting in ecosystem degradation and desertification around the world. The objective of the present work was to investigate and compare the microbial community structure and geochemical characteristics of two geographically distinct saline pan sediments in the Kalahari Desert of southern Africa. Our data suggest that these microbial communities have been shaped by geochemical drivers, including water content, salinity, and the supply of organic matter. Using Illumina 16S rRNA gene sequencing, this study provides new insights into the diversity of bacteria and archaea in semi-arid, saline, and low-carbon environments. Many of the observed taxa are halophilic and adapted to water-limiting conditions. The analysis reveals a high relative abundance of halophilic archaea (primarily Halobacteria), and the bacterial diversity is marked by an abundance of Gemmatimonadetes and spore-forming Firmicutes. In the deeper, anoxic layers, candidate division MSBL1, and acetogenic bacteria (Acetothermia) are abundant. Together, the taxonomic information and geochemical data suggest that acetogenesis could be a prevalent form of metabolism in the deep layers of a saline pan.
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.
More than a billion people rely on water from rivers sourced in High Mountain Asia (HMA), a significant portion of which is derived from snow and glacier melt. Rural communities are heavily dependent on the consistency of runoff, and are highly vulnerable to shifts in their local environment brought on by climate change. Despite this dependence, the impacts of climate change in HMA remain poorly constrained due to poor process understanding, complex terrain, and insufficiently dense in-situ measurements.
HMA's glaciers contain more frozen water than any region outside of the poles. Their extensive retreat is a highly visible and much studied marker of regional and global climate change. However, in many catchments, snow and snowmelt represent a much larger fraction of the yearly water budget than glacial meltwaters. Despite their importance, climate-related changes in HMA's snow resources have not been well studied.
Changes in the volume and distribution of snowpack have complex and extensive impacts on both local and global climates. Eurasian snow cover has been shown to impact the strength and direction of the Indian Summer Monsoon -- which is responsible for much of the precipitation over the Indian Subcontinent -- by modulating earth-surface heating. Shifts in the timing of snowmelt have been shown to limit the productivity of major rangelands, reduce streamflow, modify sediment transport, and impact the spread of vector-borne diseases. However, a large-scale regional study of climate impacts on snow resources had yet to be undertaken.
Passive Microwave (PM) remote sensing is a well-established empirical method of studying snow resources over large areas. Since 1987, there have been consistent daily global PM measurements which can be used to derive an estimate of snow depth, and hence snow-water equivalent (SWE) -- the amount of water stored in snowpack. The SWE estimation algorithms were originally developed for flat and even terrain -- such as the Russian and Canadian Arctic -- and have rarely been used in complex terrain such as HMA.
This dissertation first examines factors present in HMA that could impact the reliability of SWE estimates. Forest cover, absolute snow depth, long-term average wind speeds, and hillslope angle were found to be the strongest controls on SWE measurement reliability. While forest density and snow depth are factors accounted for in modern SWE retrieval algorithms, wind speed and hillslope angle are not. Despite uncertainty in absolute SWE measurements and differences in the magnitude of SWE retrievals between sensors, single-instrument SWE time series were found to be internally consistent and suitable for trend analysis.
Building on this finding, this dissertation tracks changes in SWE across HMA using a statistical decomposition technique. An aggregate decrease in SWE was found (10.6 mm/yr), despite large spatial and seasonal heterogeneities. Winter SWE increased in almost half of HMA, despite general negative trends throughout the rest of the year. The elevation distribution of these negative trends indicates that while changes in SWE have likely impacted glaciers in the region, climate change impacts on these two pieces of the cryosphere are somewhat distinct.
Following the discussion of relative changes in SWE, this dissertation explores changes in the timing of the snowmelt season in HMA using a newly developed algorithm. The algorithm is shown to accurately track the onset and end of the snowmelt season (70% within 5 days of a control dataset, 89% within 10). Using a 29-year time series, changes in the onset, end, and duration of snowmelt are examined. While nearly the entirety of HMA has experienced an earlier end to the snowmelt season, large regions of HMA have seen a later start to the snowmelt season. Snowmelt periods have also decreased in almost all of HMA, indicating that the snowmelt season is generally shortening and ending earlier across HMA.
By examining shifts in both the spatio-temporal distribution of SWE and the timing of the snowmelt season across HMA, we provide a detailed accounting of changes in HMA's snow resources. The overall trend in HMA is towards less SWE storage and a shorter snowmelt season. However, long-term and regional trends conceal distinct seasonal, temporal, and spatial heterogeneity, indicating that changes in snow resources are strongly controlled by local climate and topography, and that inter-annual variability plays a significant role in HMA's snow regime.
Alluvial fans are important geomorphic markers and sedimentary archives of tectonic and climatic changes. Hence, basins providing perfect studying conditions can often be found in arid regions due to the low weathering impact and thus well preservation of sedimentary features. Twelve samples for optically/infrared stimulated luminescence (OSL/IRSL) dating and one depth profile for cosmogenic radionuclide dating (10Be) were collected in the Santa Maria Valley in NW Argentina, where the exceptional preservation of several generations of alluvial fans allow exploring the external forcing conditions that led to repeated cycles of incision and aggradation. The results of the OSL/IRSL dating yielded ages ranging between 0.4 ± 0.1 ka and 271.8 ± 24.5 ka. Previous studies next to the study area indicate a depositional age of 1.5-2 Mio years for the oldest generation of alluvial fans, which might still be supported by our ongoing 10Be dating. Due to field observations, sediment provenance, stratigraphic characteristics and the geomorphic pattern of erosion, seven (/eight) generations of alluvial fan deposits were recognized. Comparing my ages with global glaciation cycles as well as linking them to temperature proxies retrieved from a lake on the Altiplano Plateau, a good fit between alluvial fan accumulation phases and global glacial periods (corresponding to cold/wet phases within the central Andes) is observed. This suggests that aggradation occurs during the early stages of glacial periods, while incision is expected at the end of glacial phases. This pattern might be linked to variations in the vegetational cover (controlled by water availability), which will decrease/increase during hot and dry/cold and wet interglacial/glacial phases favoring/limiting sediment production and will increase/decrease during cold and wet/hot and dry glacial/interglacial phases. Even though the eastern Andean margin is showing neotectonic activities and is assumed to be active up to recent times, deformation and seismicity might most probably have played only a minor role in relation to the rather short timescale reflected by the data.
Dating growth strata and basin fill by combining 26Al/10Be burial dating and magnetostratigraphy
(2018)
Cosmogenic burial dating enables dating of coarse-grained, Pliocene-Pleistocene sedimentary units that are typically difficult to date with traditional methods, such as magnetostratigraphy. In the actively deforming western Tarim Basin in NW China, Pliocene-Pleistocene conglomerates were dated at eight sites, integrating Al-26/Be-10 burial dating with previously published magnetostratigraphic sections. These samples were collected from growth strata on the flanks of growing folds and from sedimentary units beneath active faults to place timing constraints on the initiation of deformation of structures within the basin and on shortening rates on active faults. These new basin-fill and growthstrata ages document the late Neogene and Quaternary growth of the Pamir and Tian Shan orogens between >5 and 1 Ma and delineate the eastward propagation of deformation at rates up to 115 km/m.y. and basinward growth of both mountain belts at rates up to 12 km/m.y.
Previous thermomechanical modeling studies indicated that variations in the temperature and strength of the crystalline crust might be responsible for the juxtaposition of domains with thin-skinned and thick-skinned crustal deformation along strike the foreland of the central Andes. However, there is no evidence supporting this hypothesis from data-integrative models. We aim to derive the density structure of the lithosphere by means of integrated 3-D density modeling, in order to provide a new basis for discussions of compositional variations within the crust and for future thermal and rheological modeling studies. Therefore, we utilize available geological and geophysical data to obtain a structural and density model of the uppermost 200km of the Earth. The derived model is consistent with the observed Bouguer gravity field. Our results indicate that the crystalline crust in northern Argentina can be represented by a lighter upper crust (2,800kg/m(3)) and a denser lower crust (3,100kg/m(3)). We find new evidence for high bulk crustal densities >3,000kg/m(3) in the northern Pampia terrane. These could originate from subducted Puncoviscana wackes or pelites that ponded to the base of the crystalline crust in the late Proterozoic or indicate increasing bulk content of mafic material. The precise composition of the northern foreland crust, whether mafic or felsic, has significant implications for further thermomechanical models and the rheological behavior of the lithosphere. A detailed sensitivity analysis of the input parameters indicates that the model results are robust with respect to the given uncertainties of the input data.
Linking deep seismic profiles with regional-scale gravity inversion is a powerful tool to deduce the architecture of rifted margins and their structural evolution. Here we map upper and lower crustal thicknesses of the northern South China Sea (SCS) margin in order to investigate the occurrence of depth-dependent crustal extension from the proximal to the distal margin. By comparing upper and lower crustal stretching factors, we find that the northern margin of the SCS is segmented in three parts: (1) sedimentary basins where upper crust is stretched more than lower crust, (2) distal margin where lower crust is stretched more than upper crust, (3) mostly proximal margin regions where the two layers have similar stretching factors. Our results suggest that sedimentary basins and distal margin prominently feature depth-dependent extension, however accommodated by different processes. While differential thinning within sedimentary basins appears to be governed by lateral pressure variations inducing lower crustal flow, we suggest the distal margin to be affected by a combination of mantle flow-induced lower crustal shearing and sequential fault activity during crustal hyper-extension.
Terrestrial environmental systems are characterised by numerous feedback links between their different compartments. However, scientific research is organized into disciplines that focus on processes within the respective compartments rather than on interdisciplinary links. Major feedback mechanisms between compartments might therefore have been systematically overlooked so far. Without identifying these gaps, initiatives on future comprehensive environmental monitoring schemes and experimental platforms might fail. We performed a comprehensive overview of feedbacks between compartments currently represented in environmental sciences and explores to what degree missing links have already been acknowledged in the literature. We focused on process models as they can be regarded as repositories of scientific knowledge that compile findings of numerous single studies. In total, 118 simulation models from 23 model types were analysed. Missing processes linking different environmental compartments were identified based on a meta-review of 346 published reviews, model inter-comparison studies, and model descriptions. Eight disciplines of environmental sciences were considered and 396 linking processes were identified and ascribed to the physical, chemical or biological domain. There were significant differences between model types and scientific disciplines regarding implemented interdisciplinary links. The most wide-spread interdisciplinary links were between physical processes in meteorology, hydrology and soil science that drive or set the boundary conditions for other processes (e.g., ecological processes). In contrast, most chemical and biological processes were restricted to links within the same compartment. Integration of multiple environmental compartments and interdisciplinary knowledge was scarce in most model types. There was a strong bias of suggested future research foci and model extensions towards reinforcing existing interdisciplinary knowledge rather than to open up new interdisciplinary pathways. No clear pattern across disciplines exists with respect to suggested future research efforts. There is no evidence that environmental research would clearly converge towards more integrated approaches or towards an overarching environmental systems theory. (c) 2017 Elsevier B.V. All rights reserved.
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.
Small repeating earthquakes are thought to represent rupture of isolated asperities loaded by surrounding creep. The observed scaling between recurrence interval and seismic moment, T-r approximate to M-1/6, contrasts with expectation assuming constant stress drop and no aseismic slip (T-r approximate to M-1/3). Here we demonstrate that simple crack models of velocity-weakening asperities in a velocity-strengthening fault predict the M-1/6 scaling; however, the mechanism depends on asperity radius, R. For small asperities ( , where R is the nucleation radius) numerical simulations with rate-state friction show interseismic creep penetrating inward from the edge, and earthquakes nucleate in the center and rupture the entire asperity. Creep penetration accounts for approximate to 25% of the slip budget, the nucleation phase takes up a larger fraction of slip. Stress drop increases with increasing R; the lack of self-similarity being due to the finite nucleation dimension. For 2R<R less than or similar to 6Rsimulations exhibit simple cycles with ruptures nucleating from the edge. Asperities with R6R exhibit complex cycles of partial and full ruptures. Here T-r is explained by an energy criterion: full rupture requires that the energy release rate everywhere on the asperity at least equals the fracture energy, leading to the scaling T-r approximate to M-1/6. Remarkably, in spite of the variability in behavior with source dimension, the scaling of T-r with stress drop , nucleation length and creep rate v(pl) is the same across all regimes: Tr approximate to vpl. This supports the use of repeating earthquakes as creepmeters and provides a physical interpretation for the scaling observed in nature. Plain Language Summary While most earthquake sequences have complex temporal patterns, some small earthquakes are quite predictable: they repeat periodically. The time between consecutive events (recurrence interval) grows with earthquake size: as intuitive, it takes longer to accumulate the mechanical energy for large earthquakes. However, the scaling between the recurrence interval and earthquake energy (seismic moment) is not what simple physical considerations predict. It is often assumed that faults are locked between events and seismic slip must therefore keep up with long-term plate motion. This leads to the scaling: Tr approximate to M01/3, but the observed scaling is . In fact, faults are not fully locked between earthquakes: they can slip slowly, or release part of the energy in smaller quakes between the larger ones. Here we use numerical simulations, and ideas from fracture mechanics, to understand what controls the time between repeating quakes. The main results are (1) analytical expressions of the recurrence interval as a function of earthquake size, predicting the observed scaling; (2) explanation of the differences between the cycle of small and large earthquakes (fraction of slow slip, direction of rupture propagation, and the occurrence of smaller quakes between large ones) and the quantities determining these transitions.
The Aral Sea desiccation and related changes in hydroclimatic conditions on a regional level is a hot topic for past decades. The key problem of scientific research projects devoted to an investigation of modern Aral Sea basin hydrological regime is its discontinuous nature - the only limited amount of papers takes into account the complex runoff formation system entirely. Addressing this challenge we have developed a continuous prediction system for assessing freshwater inflow into the Small Aral Sea based on coupling stack of hydrological and data-driven models. Results show a good prediction skill and approve the possibility to develop a valuable water assessment tool which utilizes the power of classical physically based and modern machine learning models both for territories with complex water management system and strong water-related data scarcity. The source code and data of the proposed system is available on a Github page (https://github.com/SMASHIproject/IWRM2018).
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.
We study the rupture processes of Iquique earthquake 8.1 (2014/04/01) and its largest aftershock 7.7 (2014/04/03) that ruptured the North Chile subduction zone. High-rate Global Positioning System (GPS) recordings and strong motion data are used to reconstruct the evolution of the slip amplitude, rise time and rupture time of both earthquakes. A two-step inversion scheme is assumed, by first building prior models for both earthquakes from the inversion of the estimated static displacements and then, kinematic inversions in the frequency domain are carried out taken into account this prior information. The preferred model for the mainshock exhibits a seismic moment of 1.73 × 1021 Nm ( 8.1) and maximum slip of ∼9 m, while the aftershock model has a seismic moment of 3.88 × 1020 ( 7.7) and a maximum slip of ∼3 m. For both earthquakes, the final slip distributions show two asperities (a shallow one and a deep one) separated by an area with significant slip deficit. This suggests a segmentation along-dip which might be related to a change of the dipping angle of the subducting slab inferred from gravimetric data. Along-strike, the areas where the seismic ruptures stopped seem to be well correlated with geological features observed from geophysical information (high-resolution bathymetry, gravimetry and coupling maps) that are representative of the long-term segmentation of the subduction margin. Considering the spatially limited portions that were broken by these two earthquakes, our results support the idea that the seismic gap is not filled yet.
One paragraph of the manuscript of the paper has been inadvertently omitted in the very final stage of its compilation due to a technical mistake. Since this paragraph discusses the declustering of the used earthquake catalogue and is therefore necessary for the understanding of the seismicity data preprocessing, the authors decided to provide this paragraph in form of a correction. The respective paragraph belongs to chapter 2 of the paper, where it was placed originally, and should be inserted into the published paper before the second to the last paragraph. The omitted text reads as follows:
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.
Kim et al. recently measured the structure factor of deeply supercooled water droplets (Reports, 22 December 2017, p. 1589). We raise several concerns about their data analysis and interpretation. In our opinion, the reported data do not lead to clear conclusions about the origins of water’s anomalies.
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.
Reducing uncertainties about carbon cycling is important in the Arctic where rapid environmental changes contribute to enhanced mobilization of carbon. Here we quantify soil organic carbon (SOC) contents of permafrost soils along the Yukon Coastal Plain and determine the annual fluxes from coastal erosion. Different terrain units were assessed based on surficial geology, morphology, and ground ice conditions. To account for the volume of wedge ice and massive ice in a unit, SOC contents were reduced by 19% and sediment contents by 16%. The SOC content in a 1m(2) column of soil varied according to the height of the bluff, ranging from 30 to 662kg, with a mean value of 183kg. Forty-four per cent of the SOC was within the top 1m of soil and values varied based on surficial materials, ranging from 30 to 53kg C/m(3), with a mean of 41kg. Eighty per cent of the shoreline was erosive with a mean annual rate of change of -0.7m/yr. This resulted in a SOC flux per meter of shoreline of 132kg C/m/yr, and a total flux for the entire 282km of the Yukon coast of 35.5 x 10(6) kg C/yr (0.036 Tg C/yr). The mean flux of sediment per meter of shoreline was 5.3 x 10(3) kg/m/yr, with a total flux of 1,832 x 10(6)kg/yr (1.832 Tg/yr). Sedimentation rates indicate that approximately 13% of the eroded carbon was sequestered in nearshore sediments, where the overwhelming majority of organic carbon was of terrestrial origin. Plain Language Summary The oceans help slow the buildup of carbon dioxide (CO2) because they absorb much of this greenhouse gas. However, if carbon from other sources is added to the oceans, it can affect their ability to absorb atmospheric CO2. Our study examines the organic carbon added to the Canadian Beaufort Sea from eroding permafrost along the Yukon coast, a region quite vulnerable to erosion. Understanding carbon cycling in this area is important because environmental changes in the Arctic such as longer open water seasons, rising sea levels, and warmer air, water and soil temperatures are likely to increase coastal erosion and, thus, carbon fluxes to the sea. We measured the carbon in different types of permafrost soils and applied corrections to account for the volume taken up by various types of ground ice. By determining how quickly the shoreline is eroding, we assessed how much organic carbon is being transferred to the ocean each year. Our results show that 36 x 10(6) kg of carbon is added annually from this section of the coast. If we extrapolate these results to other coastal areas along the Canadian Beaufort Sea, the flux of organic carbon is nearly 3 times what was previously thought.
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.
The article describes the surface modification of 3D printed poly(lactic acid) (PLA) scaffolds with calcium phosphate (CP)/gelatin and CP/chitosan hybrid coating layers. The presence of gelatin or chitosan significantly enhances CP co-deposition and adhesion of the mineral layer on the PLA scaffolds. The hydrogel/CP coating layers are fairly thick and the mineral is a mixture of brushite, octacalcium phosphate, and hydroxyapatite. Mineral formation is uniform throughout the printed architectures and all steps (printing, hydrogel deposition, and mineralization) are in principle amenable to automatization. Overall, the process reported here therefore has a high application potential for the controlled synthesis of biomimetic coatings on polymeric biomaterials.
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.
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.
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.
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.
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.
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.
Causes for slow weathering and erosion in the steep, warm, monsoon-subjected Highlands of Sri Lanka
(2018)
In the Highlands of Sri Lanka, erosion and chemical weathering rates are among the lowest for global mountain denudation. In this tropical humid setting, highly weathered deep saprolite profiles have developed from high-grade metamorphic charnockite during spheroidal weathering of the bedrock. The spheroidal weathering produces rounded corestones and spalled rindlets at the rock-saprolite interface. I used detailed textural, mineralogical, chemical, and electron-microscopic (SEM, FIB, TEM) analyses to identify the factors limiting the rate of weathering front advance in the profile, the sequence of weathering reactions, and the underlying mechanisms. The first mineral attacked by weathering was found to be pyroxene initiated by in situ Fe oxidation, followed by in situ biotite oxidation. Bulk dissolution of the primary minerals is best described with a dissolution – re-precipitation process, as no chemical gradients towards the mineral surface and sharp structural boundaries are observed at the nm scale. Only the local oxidation in pyroxene and biotite is better described with an ion by ion process. The first secondary phases are oxides and amorphous precipitates from which secondary minerals (mainly smectite and kaolinite) form. Only for biotite direct solid state transformation to kaolinite is likely. The initial oxidation of pyroxene and biotite takes place in locally restricted areas and is relatively fast: log J = -11 molmin/(m2 s). However, calculated corestone-scale mineral oxidation rates are comparable to corestone-scale mineral dissolution rates: log R = -13 molpx/(m2 s) and log R = -15 molbt/(m2 s). The oxidation reaction results in a volume increase. Volumetric calculations suggest that this observed oxidation leads to the generation of porosity due to the formation of micro-fractures in the minerals and the bedrock allowing for fluid transport and subsequent dissolution of plagioclase. At the scale of the corestone, this fracture reaction is responsible for the larger fractures that lead to spheroidal weathering and to the formation of rindlets. Since these fractures have their origin from the initial oxidational induced volume increase, oxidation is the rate limiting parameter for weathering to take place. The ensuing plagioclase weathering leads to formation of high secondary porosity in the corestone over a distance of only a few cm and eventually to the final disaggregation of bedrock to saprolite. As oxidation is the first weathering reaction, the supply of O2 is a rate-limiting factor for chemical weathering. Hence, the supply of O2 and its consumption at depth connects processes at the weathering front with erosion at the surface in a feedback mechanism. The strength of the feedback depends on the relative weight of advective versus diffusive transport of O2 through the weathering profile. The feedback will be stronger with dominating diffusive transport. The low weathering rate ultimately depends on the transport of O2 through the whole regolith, and on lithological factors such as low bedrock porosity and the amount of Fe-bearing primary minerals. In this regard the low-porosity charnockite with its low content of Fe(II) bearing minerals impedes fast weathering reactions. Fresh weatherable surfaces are a pre-requisite for chemical weathering. However, in the case of the charnockite found in the Sri Lankan Highlands, the only process that generates these surfaces is the fracturing induced by oxidation. Tectonic quiescence in this region and low pre-anthropogenic erosion rate (attributed to a dense vegetation cover) minimize the rejuvenation of the thick and cohesive regolith column, and lowers weathering through the feedback with erosion.
Organic or inorganic (A) metal (M) halide (X) perovskites (AMX(3)) are semiconductor materials setting the basis for the development of highly efficient, low-cost and multijunction solar energy conversion devices. The best efficiencies nowadays are obtained with mixed compositions containing methylammonium, formamidinium, Cs and Rb as well as iodine, bromine and chlorine as anions. The understanding of fundamental properties such as crystal structure and its effect on the band gap, as well as their phase stability is essential. In this systematic study X-ray diffraction and photoluminescense spectroscopy were applied to evaluate structural and optoelectronic properties of hybrid perovskites with mixed compositions.
Fluvial terraces, floodplains, and alluvial fans are the main landforms to store sediments and to decouple hillslopes from eroding mountain rivers. Such low-relief landforms are also preferred locations for humans to settle in otherwise steep and poorly accessible terrain. Abundant water and sediment as essential sources for buildings and infrastructure make these areas amenable places to live at. Yet valley floors are also prone to rare and catastrophic sedimentation that can overload river systems by abruptly increasing the volume of sediment supply, thus causing massive floodplain aggradation, lateral channel instability, and increased flooding. Some valley-fill sediments should thus record these catastrophic sediment pulses, allowing insights into their timing, magnitude, and consequences.
This thesis pursues this theme and focuses on a prominent ~150 km2 valley fill in the Pokhara Valley just south of the Annapurna Massif in central Nepal. The Pokhara Valley is conspicuously broad and gentle compared to the surrounding dissected mountain terrain,
and is filled with locally more than 70 m of clastic debris. The area’s main river, Seti Khola, descends from the Annapurna Sabche Cirque at 3500-4500 m asl down to 900 m asl where it incises into this valley fill. Humans began to settle on this extensive
fan surface in the 1750’s when the Trans-Himalayan trade route connected the Higher Himalayas, passing Pokhara city, with the subtropical lowlands of the Terai. High and unstable river terraces and steep gorges undermined by fast flowing rivers with highly seasonal (monsoon-driven) discharge, a high earthquake risk, and a growing population make the Pokhara Valley an ideal place to study the recent geological and geomorphic history of its sediments and the implication for natural hazard appraisals.
The objective of this thesis is to quantify the timing, the sedimentologic and geomorphic processes as well as the fluvial response to a series of strong sediment pulses. I report
diagnostic sedimentary archives, lithofacies of the fan terraces, their geochemical provenance, radiocarbon-age dating and the stratigraphic relationship between them. All these various and independent lines of evidence show consistently that multiple sediment pulses filled the Pokhara Valley in medieval times, most likely in connection with, if not triggered by, strong seismic ground shaking. The geomorphic and sedimentary evidence is
consistent with catastrophic fluvial aggradation tied to the timing of three medieval Himalayan earthquakes in ~1100, 1255, and 1344 AD. Sediment provenance and calibrated radiocarbon-age data are the key to distinguish three individual sediment pulses, as these are not evident from their sedimentology alone. I explore various measures of adjustment and fluvial response of the river system following these massive aggradation pulses. By using proxies such as net volumetric erosion, incision and erosion rates, clast provenance on active river banks, geomorphic markers such as re-exhumed tree trunks in growth position, and knickpoint locations in tributary valleys, I estimate the response of the river network in the Pokhara Valley to earthquake disturbance over several centuries. Estimates of the removed volumes since catastrophic valley filling began, require average net sediment
yields of up to 4200 t km−2 yr−1 since, rates that are consistent with those reported for Himalayan rivers. The lithological composition of active channel-bed load differs from that of local bedrock material, confirming that rivers have adjusted 30-50% depending on data of different tributary catchments, locally incising with rates of 160-220 mm yr−1. In many tributaries to the Seti Khola, most of the contemporary river loads come from a Higher Himalayan source, thus excluding local hillslopes as sources. This imbalance in sediment provenance emphasizes how the medieval sediment pulses must have rapidly traversed up to 70 km downstream to invade the downstream reaches of the tributaries
up to 8 km upstream, thereby blocking the local drainage and thus reinforcing, or locally creating new, floodplain lakes still visible in the landscape today.
Understanding the formation, origin, mechanism and geomorphic processes of this valley fill is crucial to understand the landscape evolution and response to catastrophic sediment pulses. Several earthquake-triggered long-runout rock-ice avalanches or catastrophic dam burst in the Higher Himalayas are the only plausible mechanisms to explain both the geomorphic and sedimentary legacy that I document here. In any case, the Pokhara Valley was most likely hit by a cascade of extremely rare processes over some two centuries starting in the early 11th century. Nowhere in the Himalayas do we find valley fills of
comparable size and equally well documented depositional history, making the Pokhara Valley one of the most extensively dated valley fill in the Himalayas to date. Judging from the growing record of historic Himalayan earthquakes in Nepal that were traced and
dated in fault trenches, this thesis shows that sedimentary archives can be used to directly aid reconstructions and predictions of both earthquake triggers and impacts from a sedimentary-response perspective. The knowledge about the timing, evolution, and response of the Pokhara Valley and its river system to earthquake triggered sediment pulses is important to address the seismic and geomorphic risk for the city of Pokhara. This
thesis demonstrates how geomorphic evidence on catastrophic valley infill can help to independently verify paleoseismological fault-trench records and may initiate re-thinking on post-seismic hazard assessments in active mountain regions.
The precise determination of the stable C and O isotope fractionation between water and calcite (CC) and water and aragonite (AR) is of special interest for climate reconstructions, e.g. paleotemperatures. Previous studies reported results from both laboratory and field experiments, but their results are only partly consistent. Here we present C and O isotope data of a stalagmite from the Swiss Alps, which shows CC-AR transitions along individual growth layers. Using detailed analyses both laterally and perpendicular to such layers we examined the difference in the C and O isotope fractionation factor of the HCO3- - CC and the HCO3- - AR system. For O this difference is similar to the water-CC and water-AR offset provided in experimental studies. The O isotope fractionation difference in the water-CC and water-AR system is comparable to those determined in laboratory studies but shows a statistically significant correlation with the CaCO3 precipitation rate. For C we found a fractionation difference, which is independent of CaCO3 precipitation rate and with slightly smaller values for the fractionation offset between HCO3- - CC and HCO3- - AR compared to literature values. However, we also found an unexpected decrease in delta C-13 along growth layers, which contradicts the widely used concept of Rayleigh fractionation during CO2 degassing and CaCO3 precipitation. The results of this study can be used e.g., to correct stable isotope time series of stalagmites showing CC-AR transitions along their growth axes. (C) 2018 Elsevier Ltd. All rights reserved.
Ice-rich yedoma-dominated landscapes store considerable amounts of organic carbon (C) and nitrogen (N) and are vulnerable to degradation under climate warming. We investigate the C and N pools in two thermokarst-affected yedoma landscapes - on Sobo-Sise Island and on Bykovsky Peninsula in the north of eastern Siberia. Soil cores up to 3m depth were collected along geomorphic gradients and analysed for organic C and N contents. A high vertical sampling density in the profiles allowed the calculation of C and N stocks for short soil column intervals and enhanced understanding of within-core parameter variability. Profile-level C and N stocks were scaled to the landscape level based on landform classifications from 5 m resolution, multispectral RapidEye satellite imagery. Mean landscape C and N storage in the first metre of soil for Sobo-Sise Island is estimated to be 20.2 kg C m(-2) and 1.8 kg N m(-2) and for Bykovsky Peninsula 25.9 kg C m(-2) and 2.2 kg N m(-2). Radiocarbon dating demonstrates the Holocene age of thermokarst basin deposits but also suggests the presence of thick Holoceneage cover layers which can reach up to 2 m on top of intact yedoma landforms. Reconstructed sedimentation rates of 0.10-0.57 mm yr(-1) suggest sustained mineral soil accumulation across all investigated landforms. Both yedoma and thermokarst landforms are characterized by limited accumulation of organic soil layers (peat). We further estimate that an active layer deepening of about 100 cm will increase organic C availability in a seasonally thawed state in the two study areas by similar to 5.8 Tg (13.2 kg C m(-2)). Our study demonstrates the importance of increasing the number of C and N storage inventories in ice-rich yedoma and thermokarst environments in order to account for high variability of permafrost and thermokarst environments in pan-permafrost soil C and N pool estimates.
For the first time, a rain signature in Global Navigation Satellite System Reflectometry (GNSS-R) observations is demonstrated. Based on the argument that the forward quasi-specular scattering relies upon surface gravity waves with lengths larger than several wavelengths of the reflected signal, a commonly made conclusion is that the scatterometric GNSS-R measurements are not sensitive to the surface small-scale roughness generated by raindrops impinging on the ocean surface. On the contrary, this study presents an evidence that the bistatic radar cross section sigma(0) derived from TechDemoSat-1 data is reduced due to rain at weak winds, lower than approximate to 6 m/s. The decrease is as large as approximate to 0.7 dB at the wind speed of 3 m/s due to a precipitation of 0-2 mm/hr. The simulations based on the recently published scattering theory provide a plausible explanation for this phenomenon which potentially enables the GNSS-R technique to detect precipitation over oceans at low winds.
Basaltic fissure eruptions, such as on Hawai'i or on Iceland, are thought to be driven by the lateral propagation of feeder dikes and graben subsidence. Associated solid earth processes, such as deformation and structural development, are well studied by means of geophysical and geodetic technologies. The eruptions themselves, lava fountaining and venting dynamics, in turn, have been much less investigated due to hazardous access, local dimension, fast processes, and resulting poor data availability.
This thesis provides a detailed quantitative understanding of the shape and dynamics of lava fountains and the morphological changes at their respective eruption sites. For this purpose, I apply image processing techniques, including drones and fixed installed cameras, to the sequence of frames of video records from two well-known fissure eruptions in Hawai'i and Iceland. This way I extract the dimensions of multiple lava fountains, visible in all frames. By putting these results together and considering the acquisition times of the frames I quantify the variations in height, width and eruption velocity of the lava fountains. Then I analyse these time-series in both time and frequency domains and investigate the similarities and correlations between adjacent lava fountains. Following this procedure, I am able to link the dynamics of the individual lava fountains to physical parameters of the magma transport in the feeder dyke of the fountains.
The first case study in this thesis focuses on the March 2011 Pu'u'O'o eruption, Hawai'i, where a continuous pulsating behaviour at all eight lava fountains has been observed. The lava fountains, even those from different parts of the fissure that are closely connected, show a similar frequency content and eruption behaviour. The regular pattern in the heights of lava fountain suggests a controlling process within the magma feeder system like a hydraulic connection in the underlying dyke, affecting or even controlling the pulsating behaviour.
The second case study addresses the 2014-2015 Holuhraun fissure eruption, Iceland. In this case, the feeder dyke is highlighted by the surface expressions of graben-like structures and fault systems. At the eruption site, the activity decreases from a continuous line of fire of ~60 vents to a limited number of lava fountains. This can be explained by preferred upwards magma movements through vertical structures of the pre-eruptive morphology. Seismic tremors during the eruption reveal vent opening at the surface and/or pressure changes in the feeder dyke. The evolving topography of the cinder cones during the eruption interacts with the lava fountain behaviour. Local variations in the lava fountain height and width are controlled by the conduit diameter, the depth of the lava pond and the shape of the crater. Modelling of the fountain heights shows that long-term eruption behaviour is controlled mainly by pressure changes in the feeder dyke.
This research consists of six chapters with four papers, including two first author and two co-author papers. It establishes a new method to analyse lava fountain dynamics by video monitoring. The comparison with the seismicity, geomorphologic and structural expressions of fissure eruptions shows a complex relationship between focussed flow through dykes, the morphology of the cinder cones, and the lava fountain dynamics at the vents of a fissure eruption.
Thermal conductivity (lambda) is an essential physical property of minerals and rocks and fundamental in constraining the thermal field of the lithosphere. In case that adequate samples to measure lambda are not available, it could be indirectly inferred from calculation. One of the most widely applied indirect methods for rocks involve modal mineralogy and porosity as parameters that are incorporated into mathematical mean or mixing models. Robust inferences from these approaches for crystalline rocks were impeded by a small number of studied samples or restriction to certain rock types. We employ this method and examine its applicability to low-porosity plutonic rocks by calculating bulk thermal conductivity lambda(b) for 45 samples covering the entire range from gabbro/diorite to granite. We show that the use of the harmonic-mean model for both rock matrix and porosity provided a good match between lambda(b.meas) and lambda(b.calc) of <10% deviation (2 sigma), with relative and absolute errors amounting to 1.49.7% and 4.44.9%, respectively. The results of our study constitute a big step forward to a robust conclusion on the overall applicability of the harmonic-mean model for inferring lambda(b) of isotropic, low-porosity, mafic to silicic plutonic and metamorphic rocks with an acceptable magnitude of error. Drill cuttings and enclaves form particularly interesting objects for application of this method, as they are poorly suited for direct measurement. Well-derived lambda values for those rocks would permit to calculate heat flow and to model more profoundly the thermal state of the deeper lithosphere.
Global climate change is one of the greatest challenges of the 21st century, with influence on the environment, societies, politics and economies. The (semi-)arid areas of Southern Africa already suffer from water scarcity. There is a great variety of ongoing research related to global climate history but important questions on regional differences still exist.
In southern African regions terrestrial climate archives are rare, which makes paleoclimate studies challenging. Based on the assumption that continental pans (sabkhas) represent a suitable geo-archive for the climate history, two different pans were studied in the southern and western Kalahari Desert. A combined approach of molecular biological and biogeochemical analyses is utilized to investigate the diversity and abundance of microorganisms and to trace temporal and spatial changes in paleoprecipitation in arid environments. The present PhD thesis demonstrates the applicability of pan sediments as a late Quaternary geo-archive based on microbial signature lipid biomarkers, such as archaeol, branched and isoprenoid glycerol dialkyl glycerol tetraethers (GDGTs) as well as phospholipid fatty acids (PLFA). The microbial signatures contained in the sediment provide information on the current or past microbial community from the Last Glacial Maximum to the recent epoch, the Holocene. The results are discussed in the context of regional climate evolution in southwestern Africa. The seasonal shift of the Innertropical Convergence Zone (ITCZ) along the equator influences the distribution of precipitation- and climate zones. The different expansion of the winter- and summer rainfall zones in southern Africa was confirmed by the frequency of certain microbial biomarkers. A period of increased precipitation in the south-western Kalahari could be described as a result of the extension of the winter rainfall zone during the last glacial maximum (21 ± 2 ka). Instead a period of increased paleoprecipitation in the western Kalahari was indicated during the Late Glacial to Holocene transition. This was possibly caused by a southwestern shift in the position of the summer rainfall zone associated to the southward movement of the ITCZ.
Furthermore, for the first time this study characterizes the bacterial and archaeal life based on 16S rRNA gene high-throughput sequencing in continental pan sediments and provides an insight into the recent microbial community structure. Near-surface processes play an important role for the modern microbial ecosystem in the pans. Water availability as well as salinity might determine the abundance and composition of the microbial communities. The microbial community of pan sediments is dominated by halophilic and dry-adapted archaea and bacteria. Frequently occurring microorganisms such as, Halobacteriaceae, Bacillus and Gemmatimonadetes are described in more detail in this study.
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.
Great megathrust earthquakes arise from the sudden release of energy accumulated during centuries of interseismic plate convergence. The moment deficit (energy available for future earthquakes) is commonly inferred by integrating the rate of interseismic plate locking over the time since the previous great earthquake. But accurate integration requires knowledge of how interseismic plate locking changes decades after earthquakes, measurements not available for most great earthquakes. Here we reconstruct the post-earthquake history of plate locking at Guafo Island, above the seismogenic zone of the giant 1960 (M-w = 9.5) Chile earthquake, through forward modeling of land-level changes inferred from aerial imagery (since 1974) and measured by GPS (since 1994). We find that interseismic locking increased to similar to 70% in the decade following the 1960 earthquake and then gradually to 100% by 2005. Our findings illustrate the transient evolution of plate locking in Chile, and suggest a similarly complex evolution elsewhere, with implications for the time- and magnitude-dependent probability of future events.
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.
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.
For several decades the Sant Llorenc de Montgai reservoir has experienced different problems that could affect the safety of the engineering structure. For this reason, several corrective actions have been taken over the years. Here, we present a study involving complementary geophysical methods including electrical resistivity tomography, seismic refraction tomography and frequency-domain electromagnetic surveys. The analysis of the inverted electrical resistivity tomography cross-sections combined with the seismic refraction results and land subsidence monitoring data show the likely mechanism of abnormal seepage. The areas where mortar injections were applied as a corrective measure are also clearly delineated. In addition, the evolution of the state of the embankment has been established from two successive electrical resistivity tomography surveys in the last two decades. The results show areas where corrective mortar injections have been effective, while in other areas new abnormal seepage has been detected. The lithological heterogeneity of the bedrock, especially the dissolution of gypsum-rich rocks, induced subsidence effects and caused abnormal seepage in different areas along the embankment. Our results indicate how corrective solutions can be optimized to reduce the cost of corrective engineering interventions.
The oldest ice core records are obtained from the East Antarctic Plateau. Water isotopes are key proxies to reconstructing past climatic conditions over the ice sheet and at the evaporation source. The accuracy of climate reconstructions depends on knowledge of all processes affecting water vapour, precipitation and snow isotopic compositions. Fractionation processes are well understood and can be integrated in trajectory-based Rayleigh distillation and isotope-enabled climate models. However, a quantitative understanding of processes potentially altering snow isotopic composition after deposition is still missing. In low-accumulation sites, such as those found in East Antarctica, these poorly constrained processes are likely to play a significant role and limit the interpretability of an ice core's isotopic composition.
By combining observations of isotopic composition in vapour, precipitation, surface snow and buried snow from Dome C, a deep ice core site on the East Antarctic Plateau, we found indications of a seasonal impact of metamorphism on the surface snow isotopic signal when compared to the initial precipitation. Particularly in summer, exchanges of water molecules between vapour and snow are driven by the diurnal sublimation-condensation cycles. Overall, we observe in between precipitation events modification of the surface snow isotopic composition. Using high-resolution water isotopic composition profiles from snow pits at five Antarctic sites with different accumulation rates, we identified common patterns which cannot be attributed to the seasonal variability of precipitation. These differences in the precipitation, surface snow and buried snow isotopic composition provide evidence of post-deposition processes affecting ice core records in low-accumulation areas.
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.
Fibroblast growth factor 23 (FGF23) is a proteohormone regulating renal phosphate transport and vitamin D metabolism as well as inducing left heart hypertrophy. FGF23-deficient mice suffer from severe tissue calcification, accelerated aging and a myriad of aging-associated diseases. Bone cells produce FGF23 upon store-operated calcium ion entry (SOCE) through the calcium selective ion channel Orai1. AMP-activated kinase (AMPK) is a powerful energy sensor helping cells survive states of energy deficiency, and AMPK down-regulates Orai1. Here we investigated the role of AMPK in FGF23 production. Fgf23 gene transcription was analyzed by qRT-PCR and SOCE by fluorescence optics in UMR106 osteoblast-like cells while the serum FGF23 concentration and phosphate metabolism were assessed in AMPKa1-knockout and wild-type mice. The AMPK activator, 5-aminoimidazole-4-carboxamide ribonucleotide (AICAR) down-regulated, whereas the AMPK inhibitor, dorsomorphin dihydrochloride (compound C) and AMPK gene silencing induced Fgf23 transcription. AICAR decreased membrane abundance of Orai1 and SOCE. SOCE inhibitors lowered Fgf23 gene expression induced by AMPK inhibition. AMPKa1-knockout mice had a higher serum FGF23 concentration compared to wild-type mice. Thus, AMPK participates in the regulation of FGF23 production in vitro and in vivo. The inhibitory effect of AMPK on FGF23 production is at least in part mediated by Orai1-involving SOCE.
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.
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.
Identifying abrupt transitions is a key question in various disciplines. Existing transition detection methods, however, do not rigorously account for time series uncertainties, often neglecting them altogether or assuming them to be independent and qualitatively similar. Here, we introduce a novel approach suited to handle uncertainties by representing the time series as a time-ordered sequence of probability density functions. We show how to detect abrupt transitions in such a sequence using the community structure of networks representing probabilities of recurrence. Using our approach, we detect transitions in global stock indices related to well-known periods of politico-economic volatility. We further uncover transitions in the El Niño-Southern Oscillation which coincide with periods of phase locking with the Pacific Decadal Oscillation. Finally, we provide for the first time an ‘uncertainty-aware’ framework which validates the hypothesis that ice-rafting events in the North Atlantic during the Holocene were synchronous with a weakened Asian summer monsoon.
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.
The central European Bohemian Massif has undergone over two centuries of scientific investigation which has made it a pivotal area for the development and testing of modern geological theories. The discovery of melt inclusions in high-grade rocks, either crystallized as nanogranitoids or as glassy inclusions, prompted the re-evaluation of the area with an ‘inclusionist’ eye. Melt inclusions have been identified in a wide range of rocks, including felsic/perpotassic granulites, migmatites, eclogites and garnet clinopyroxenites, all the result of melting events albeit over a wide range of pressure/temperature conditions (800–1000°C/0.5–5 GPa). This contribution provides an overview of such inclusions and discusses the qualitative and quantitative constraints they provide for melting processes, and the nature of melts and fluids involved in these processes. In particular, data on trace-element signatures of melt inclusions trapped at mantle depths are presented and discussed. Moreover, experimental re-homogenization of nanogranitoids provided microstructural criteria allowing assessment of the conditions at which melt and host are mutually stable during melting. Overall this work aims to provide guidelines and suggestions for petrologists wishing to explore the fascinating field of melt inclusions in metamorphic terranes worldwide, based on the newest discoveries from the still-enigmatic Bohemian Massif.
A transparent and data-driven global tectonic regionalization model for seismic hazard assessment
(2018)
A key concept that is common to many assumptions inherent within seismic hazard assessment is that of tectonic similarity. This recognizes that certain regions of the globe may display similar geophysical characteristics, such as in the attenuation of seismic waves, the magnitude scaling properties of seismogenic sources or the seismic coupling of the lithosphere. Previous attempts at tectonic regionalization, particularly within a seismic hazard assessment context, have often been based on expert judgements; in most of these cases, the process for delineating tectonic regions is neither reproducible nor consistent from location to location. In this work, the regionalization process is implemented in a scheme that is reproducible, comprehensible from a geophysical rationale, and revisable when new relevant data are published. A spatial classification-scheme is developed based on fuzzy logic, enabling the quantification of concepts that are approximate rather than precise. Using the proposed methodology, we obtain a transparent and data-driven global tectonic regionalization model for seismic hazard applications as well as the subjective probabilities (e.g. degree of being active/degree of being cratonic) that indicate the degree to which a site belongs in a tectonic category.
A submerged pine forest from the early Holocene in the Mecklenburg Lake District, northern Germany
(2018)
For the first time, evidence of a submerged pine forest from the early Holocene can be documented in a central European lake. Subaquatic tree stumps were discovered in Lake Giesenschlagsee at a depth of between 2 and 5m using scuba divers, side-scan sonar and a remotely operated vehicle. Several erect stumps, anchored to the ground by roots, represent an insitu record of this former forest. Botanical determination revealed the stumps to be Scots pine (Pinus sylvestris) with an individual tree age of about 80years. The trees could not be dated by means of dendrochronology, as they are older than the regional reference chronology for pine. Radiocarbon ages from the wood range from 10880 +/- 210 to 10370 +/- 130cal. a BP, which is equivalent to the mid-Preboreal to early Boreal biozones. The trees are rooted in sedge peat, which can be dated to this period as well, using pollen stratigraphical analysis. Tilting of the peat bed by 4m indicates subsidence of the ground due to local dead ice melting, causing the trees to become submerged and preserved for millennia. Together with recently detected Lateglacial insitu tree occurrences in nearby lakes, the submerged pine forest at Giesenschlagsee represents a new and highly promising type of geo-bio-archive for the wider region. Comparable insitu pine remnants occur at some terrestrial (buried setting) and marine (submerged setting) sites in northern central Europe and beyond, but they partly differ in age. In general, the insitu pine finds document shifts of the zonal boreal forest ecosystem during the late Quaternary.
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.
Like almost all fields of science, hydrology has benefited to a large extent from the tremendous improvements in scientific instruments that are able to collect long-time data series and an increase in available computational power and storage capabilities over the last decades. Many model applications and statistical analyses (e.g., extreme value analysis) are based on these time series. Consequently, the quality and the completeness of these time series are essential. Preprocessing of raw data sets by filling data gaps is thus a necessary procedure. Several interpolation techniques with different complexity are available ranging from rather simple to extremely challenging approaches. In this paper, various imputation methods available to the hydrological researchers are reviewed with regard to their suitability for filling gaps in the context of solving hydrological questions. The methodological approaches include arithmetic mean imputation, principal component analysis, regression-based methods and multiple imputation methods. In particular, autoregressive conditional heteroscedasticity (ARCH) models which originate from finance and econometrics will be discussed regarding their applicability to data series characterized by non-constant volatility and heteroscedasticity in hydrological contexts. The review shows that methodological advances driven by other fields of research bear relevance for a more intensive use of these methods in hydrology. Up to now, the hydrological community has paid little attention to the imputation ability of time series models in general and ARCH models in particular.
Python is at the forefront of scientific computation for seismologists and therefore should be introduced to students interested in becoming seismologists. On its own, Python is open source and well designed with extensive libraries. However, Python code can also be executed, visualized, and communicated to others with "Jupyter Notebooks". Thus, Jupyter Notebooks are ideal for teaching students Python and scientific computation. In this article, we designed an openly available Python library and collection of Jupyter Notebooks based on defined scientific computation learning goals for seismology students. The Notebooks cover topics from an introduction to Python to organizing data, earthquake catalog statistics, linear regression, and making maps. Our Python library and collection of Jupyter Notebooks are meant to be used as course materials for an upper-division data analysis course in an Earth Science Department, and the materials were tested in a Probabilistic Seismic Hazard course. However, seismologists or anyone else who is interested in Python for data analysis and map making can use these materials.
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.
Sulfur is an important component in volcanic gases at the Earth surface but also present in the deep Earth in hydrothermal or magmatic fluids. Little is known about the evolution of such fluids during ascent in the crust. A new optical cell was developed for in situ Raman spectroscopic investigations on fluids allowing abrupt or continuous changes of pressure up to 200 MPa at temperatures up to 750 degrees C. The concept is based on a flexible gold bellow, which separates the sample fluid from the pressure medium water. To avoid reactions between aggressive fluids and the pressure cell, steel components in contact with the fluid are shielded by gold foil. The cell was tested to study redox reactions in fluids using aqueous ammonium sulfate solutions as a model system. During heating at constant pressure of 130 MPa, sulfate ions transform first to HSO4- ions and then to molecular units such as H2SO4. Variation of pressure shows that the stability of sulfate species relies on fluid density, i.e., highly charged species are stable only in high-density fluids. Partial decomposition of ammonium was evident above 550 degrees C by the occurrence of a nitrogen peak in the Raman spectra. Reduced sulfur species were observed above 700 degrees C by Raman signals near 2590 cm(-1) assigned to HS- and H2S. No clear evidence for the formation of sulfur dioxide was found in contrary to previous studies on aqueous H2SO4, suggesting very reducing conditions in our experiments. Fluid-mineral interaction was studied by inserting into the cell a small, semi-open capsule filled with a mixture of pyrite and pyrrhotite. Oxidation of the sample assembly was evident by transformation of pyrite to pyrrhotite. As a consequence, sulfide species were observed in the fluid already at temperatures of similar to 600 degrees C.
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.
With increasing amount of strong motion data, Ground Motion Prediction Equation (GMPE) developers are able to quantify empirical site amplification functions (delta S2S(s)) from GMPE residuals, for use in site-specific Probabilistic Seismic Hazard Assessment. In this study, we first derive a GMPE for 5% damped Pseudo Spectral Acceleration (g) of Active Shallow Crustal earthquakes in Japan with 3.4 <= M-w <= 7.3 and 0 <= R-JB <= 600km. Using k-mean spectral clustering technique, we then classify our estimated delta S2S(s)(T = 0.01 - 2s) of 588 wellcharacterized sites, into 8 site clusters with distinct mean site amplification functions, and within-cluster site-tosite variability similar to 50% smaller than the overall dataset variability (phi(S2S)). Following an evaluation of existing schemes, we propose a revised data-driven site classification characterized by kernel density distributions of V-s30, V-s10, H-800, and predominant period (T-G) of the site clusters.
A fast and sensitive method for the continuous determination of methane (CH4) and its stable carbon isotopic values (delta C-13-CH4) in surface waters was developed by applying a vacuum to a gas/liquid exchange membrane and measuring the extracted gases by a portable cavity ring-down spectroscopy analyser (M-CRDS). The M-CRDS was calibrated and characterized for CH4 concentration and delta C-13-CH4 with synthetic water standards. The detection limit of the M-CRDS for the simultaneous determination of CH4 and delta C-13-CH4 is 3.6 nmol L-1 CH4. A measurement precision of CH4 concentrations and delta C-13-CH4 in the range of 1.1%, respectively, 1.7 parts per thousand (1 sigma) and accuracy (1.3%, respectively, 0.8 parts per thousand [1 sigma]) was achieved for single measurements and averaging times of 10 min. The response time tau of 57 +/- 5 s allow determination of delta C-13-CH4 values more than twice as fast than other methods. The demonstrated M-CRDS method was applied and tested for Lake Stechlin (Germany) and compared with the headspace-gas chromatography and fast membrane CH4 concentration methods. Maximum CH4 concentrations (577 nmol L-1) and lightest delta C-13-CH4 (-35.2 parts per thousand) were found around the thermocline in depth profile measurements. The M-CRDS-method was in good agreement with other methods. Temporal variations in CH4 concentration and delta C-13-CH4 obtained in 24 h measurements indicate either local methane production/oxidation or physical variations in the thermocline. Therefore, these results illustrate the need of fast and sensitive analyses to achieve a better understanding of different mechanisms and pathways of CH4 formation in aquatic environments.
Earthquake localization is both a necessity within the field of seismology, and a prerequisite for further analysis such as source studies and hazard assessment. Traditional localization methods often rely on manually picked phases. We present an alternative approach using deep learning that once trained can predict hypocenter locations efficiently. In seismology, neural networks have typically been trained with either single-station records or based on features that have been extracted previously from the waveforms. We use three-component full-waveform records of multiple stations directly. This means no information is lost during preprocessing and preparation of the data does not require expert knowledge. The first convolutional layer of our deep convolutional neural network (CNN) becomes sensitive to features that characterize the waveforms it is trained on. We show that this layer can therefore additionally be used as an event detector. As a test case, we trained our CNN using more than 2000 earthquake swarm events from West Bohemia, recorded by nine local three-component stations. The CNN successfully located 908 validation events with standard deviations of 56.4 m in east-west, 123.8 m in north-south, and 136.3 m in vertical direction compared to a double-difference relocated reference catalog. The detector is sensitive to events with magnitudes down to M-L = -0.8 with 3.5% false positive detections.