Refine
Year of publication
Document Type
- Article (2826)
- Doctoral Thesis (508)
- Postprint (139)
- Other (76)
- Review (52)
- Monograph/Edited Volume (34)
- Preprint (17)
- Conference Proceeding (13)
- Habilitation Thesis (12)
- Master's Thesis (6)
Keywords
- climate change (52)
- Holocene (44)
- erosion (28)
- permafrost (27)
- Himalaya (26)
- remote sensing (24)
- Climate change (23)
- Tibetan Plateau (22)
- Earthquake source observations (21)
- Pollen (21)
Institute
- Institut für Geowissenschaften (3691) (remove)
The sedimentary record of the Dead Sea provides an exceptional high-resolution archive of past climate changes in the drought-sensitive eastern Mediterranean-Levant, a key region for the development of humankind at the boundary of global climate belts. Moreover, it is the only deep hypersaline lake known to have deposited long sequences of finely laminated, annually deposited sediments (i.e. varves) of varied compositions, including aragonite, gypsum, halite and clastic sediments. Vast efforts have been made over the years to decipher the environmental information stored in these evaporitic-clastic sequences spanning from the Pleistocene Lake Amora to the Holocene Dead Sea. A general characterisation of sediment facies has been derived from exposed sediment sections, as well as from shallow- and deep-water sediment cores. During high lake stands and episodes of positive water budget, mostly during glacial times, alternating aragonite and detritus laminae (‘aad’ facies) were accumulated, whereas during low lake stands and droughts, prevailing during interglacials, laminated detritus (‘ld’ facies) and laminated halite (‘lh’ facies) dominate the sequence. In this paper, we (i) review the three types of laminated sediments of the Dead Sea sedimentary record (‘aad’, ‘ld’ and ‘lh’ facies), (ii) discuss their modes of formation, deposition and accumulation, and their interpretation as varves, and (iii) illustrate how Dead Sea varves are utilized for palaeoclimate reconstructions and for establishing floating chronologies.
The main intention of the PhD project was to create a varve chronology for the Suigetsu Varves 2006' (SG06) composite profile from Lake Suigetsu (Japan) by thin section microscopy. The chronology was not only to provide an age-scale for the various palaeo-environmental proxies analysed within the SG06 project, but also and foremost to contribute, in combination with the SG06 14C chronology, to the international atmospheric radiocarbon calibration curve (IntCal). The SG06 14C data are based on terrestrial leaf fossils and therefore record atmospheric 14C values directly, avoiding the corrections necessary for the reservoir ages of the marine datasets, which are currently used beyond the tree-ring limit in the IntCal09 dataset (Reimer et al., 2009). The SG06 project is a follow up of the SG93 project (Kitagawa & van der Plicht, 2000), which aimed to produce an atmospheric calibration dataset, too, but suffered from incomplete core recovery and varve count uncertainties. For the SG06 project the complete Lake Suigetsu sediment sequence was recovered continuously, leaving the task to produce an improved varve count. Varve counting was carried out using a dual method approach utilizing thin section microscopy and micro X-Ray Fluorescence (µXRF). The latter was carried out by Dr. Michael Marshall in cooperation with the PhD candidate. The varve count covers 19 m of composite core, which corresponds to the time frame from ≈10 to ≈40 kyr BP. The count result showed that seasonal layers did not form in every year. Hence, the varve counts from either method were incomplete. This rather common problem in varve counting is usually solved by manual varve interpolation. But manual interpolation often suffers from subjectivity. Furthermore, sedimentation rate estimates (which are the basis for interpolation) are generally derived from neighbouring, well varved intervals. This assumes that the sedimentation rates in neighbouring intervals are identical to those in the incompletely varved section, which is not necessarily true. To overcome these problems a novel interpolation method was devised. It is computer based and automated (i.e. avoids subjectivity and ensures reproducibility) and derives the sedimentation rate estimate directly from the incompletely varved interval by statistically analysing distances between successive seasonal layers. Therefore, the interpolation approach is also suitable for sediments which do not contain well varved intervals. Another benefit of the novel method is that it provides objective interpolation error estimates. Interpolation results from the two counting methods were combined and the resulting chronology compared to the 14C chronology from Lake Suigetsu, calibrated with the tree-ring derived section of IntCal09 (which is considered accurate). The varve and 14C chronology showed a high degree of similarity, demonstrating that the novel interpolation method produces reliable results. In order to constrain the uncertainties of the varve chronology, especially the cumulative error estimates, U-Th dated speleothem data were used by linking the low frequency 14C signal of Lake Suigetsu and the speleothems, increasing the accuracy and precision of the Suigetsu calibration dataset. The resulting chronology also represents the age-scale for the various palaeo-environmental proxies analysed in the SG06 project. One proxy analysed within the PhD project was the distribution of event layers, which are often representatives of past floods or earthquakes. A detailed microfacies analysis revealed three different types of event layers, two of which are described here for the first time for the Suigetsu sediment. The types are: matrix supported layers produced as result of subaqueous slope failures, turbidites produced as result of landslides and turbidites produced as result of flood events. The former two are likely to have been triggered by earthquakes. The vast majority of event layers was related to floods (362 out of 369), which allowed the construction of a respective chronology for the last 40 kyr. Flood frequencies were highly variable, reaching their greatest values during the global sea level low-stand of the Glacial, their lowest values during Heinrich Event 1. Typhoons affecting the region represent the most likely control on the flood frequency, especially during the Glacial. However, also local, non-climatic controls are suggested by the data. In summary, the work presented here expands and revises knowledge on the Lake Suigetsu sediment and enabls the construction of a far more precise varve chronology. The 14C calibration dataset is the first such derived from lacustrine sediments to be included into the (next) IntCal dataset. References: Kitagawa & van der Plicht, 2000, Radiocarbon, Vol 42(3), 370-381 Reimer et al., 2009, Radiocarbon, Vol 51(4), 1111-1150
Rheology describes the flow of matter under the influence of stress, and - related to solids- it investigates how solids subjected to stresses deform. As the deformation of the Earth’s outer layers, the lithosphere and the crust, is a major focus of rheological studies, rheology in the geosciences describes how strain evolves in rocks of variable composition and temperature under tectonic stresses. It is here where deformation processes shape the form of ocean basins and mountain belts that ultimately result from the complex interplay between lithospheric plate motion and the susceptibility of rocks to the influence of plate-tectonic forces. A rigorous study of the strength of the lithosphere and deformation phenomena thus requires in-depth studies of the rheological characteristics of the involved materials and the temporal framework of deformation processes.
This dissertation aims at analyzing the influence of the physical configuration of the lithosphere on the present-day thermal field and the overall rheological characteristics of the lithosphere to better understand variable expressions in the formation of passive continental margins and the behavior of strike-slip fault zones. The main methodological approach chosen is to estimate the present-day thermal field and the strength of the lithosphere by 3-D numerical modeling. The distribution of rock properties is provided by 3-D structural models, which are used as the basis for the thermal and rheological modeling. The structural models are based on geophysical and geological data integration, additionally constrained by 3-D density modeling. More specifically, to decipher the thermal and rheological characteristics of the lithosphere in both oceanic and continental domains, sedimentary basins in the Sea of Marmara (continental transform setting), the SW African passive margin (old oceanic crust), and the Norwegian passive margin (young oceanic crust) were selected for this study.
The Sea of Marmara, in northwestern Turkey, is located where the dextral North Anatolian Fault zone (NAFZ) accommodates the westward escape of the Anatolian Plate toward the Aegean. Geophysical observations indicate that the crust is heterogeneous beneath the Marmara basin, but a detailed characterization of the lateral crustal heterogeneities is presented for the first time in this study. Here, I use different gravity datasets and the general non-uniqueness in potential field modeling, to propose three possible end-member scenarios of crustal configuration. The models suggest that pronounced gravitational anomalies in the basin originate from significant density heterogeneities within the crust. The rheological modeling reveals that associated variations in lithospheric strength control the mechanical segmentation of the NAFZ. Importantly, a strong crust that is mechanically coupled to the upper mantle spatially correlates with aseismic patches where the fault bends and changes its strike in response to the presence of high-density lower crustal bodies. Between the bends, mechanically weaker crustal domains that are decoupled from the mantle are characterized by creep.
For the passive margins of SW Africa and Norway, two previously published 3-D conductive and lithospheric-scale thermal models were analyzed. These 3-D models differentiate various sedimentary, crustal, and mantle units and integrate different geophysical data, such as seismic observations and the gravity field. Here, the rheological modeling suggests that the present-day lithospheric strength across the oceanic domain is ultimately affected by the age and past thermal and tectonic processes as well as the depth of the thermal lithosphere-asthenosphere boundary, while the configuration of the crystalline crust dominantly controls the rheological behavior of the lithosphere beneath the continental domains of both passive margins.
The thermal and rheological models show that the variations of lithospheric strength are fundamentally influenced by the temperature distribution within the lithosphere. Moreover, as the composition of the lithosphere significantly influences the present-day thermal field, it therefore also affects the rheological characteristics of the lithosphere. Overall my studies add to our understanding of regional tectonic deformation processes and the long-term behavior of sedimentary basins; they confirm other analyses that have pointed out that crustal heterogeneities in the continents result in diverse lithospheric thermal characteristics, which in turn results in higher complexity and variations of rheological behavior compared to oceanic domains with a thinner, more homogeneous crust.
Fluvial planation surfaces, such as straths, commonly serve as recorders of climatic and tectonic changes and are formed by the lateral erosion of rivers, a process that remains poorly understood. Here we present a study of kilometer-wide, fluvially eroded, low-relief surfaces on rapidly uplifting folds in the foreland of the southwestern Tian Shan. A combination of field work, digital elevation model analysis, and dating of fluvial deposits reveals that despite an arid climate and rapid average rock-uplift rates of 1-3mm/yr, rivers cut extensive (>1-2km wide) surfaces with typical height variations of <6m over periods of >2-6kyr. The extent of this beveling varies in space and time, such that different beveling episodes affect individual structures. Between times of planation, beveled surfaces are abandoned, incised, and deformed across the folds. In a challenge to models that link strath cutting and abandonment primarily to changes in river incision rates, we demonstrate that lateral erosion rates of antecedent streams crossing the folds have to vary by more than 1 order of magnitude to explain the creation of beveled platforms in the past and their incision at the present day. These variations do not appear to covary with climate variability and might be caused by relatively small (much less than an order of magnitude) changes in sediment or water fluxes. It remains uncertain in which settings variations in lateral bedrock erosion rates predominate over changes in vertical erosion rates. Therefore, when studying fluvial planation and strath terraces, variability of both lateral and vertical erosion rates should be considered.
Rivers transfer particulate organic carbon (POC) from eroding mountains into geological sinks. Organic carbon source composition and selective mobilization have been shown to affect the type and quantity of POC export, but their combined effects across complex mountain ranges remain underexplored. Here, we examine the variation in organic carbon sourcing and transport in the trans-Himalayan Kali Gandaki River catchment, along strong gradients in precipitation, rock type and vegetation. Combining bulk stable nitrogen, and stable and radioactive organic carbon isotopic composition of bedrock, litter, soil and river sediment samples with a Bayesian end-member mixing approach, we differentiate POC sources along the river and quantify their export. Our analysis shows that POC export from the Tibetan segment of the catchment, where carbon bearing shales are partially covered by aged and modern soils, is dominated by petrogenic POC. Based on our data we re-assess the presence of aged biospheric OC in this part of the catchment, and its contribution to the river load. In the High Himalayan segment, we observed low inputs of petrogenic and biospheric POC, likely due to very low organic carbon concentrations in the metamorphic bedrock, combined with erosion dominated by deep-seated landslides. Our findings show that along the Kali Gandaki River, the sourcing of sediment and organic carbon are decoupled, due to differences in rock organic carbon content, soil and above ground carbon stocks, and geomorphic process activity. While the fast eroding High Himalayas are the principal source of river sediment, the Tibetan headwaters, where erosion rates are lower, are the principal source of organic carbon. To robustly estimate organic carbon export from the Himalayas, the mountain range should be divided into tectono-physiographic zones with distinct organic carbon yields due to differences in substrate and erosion processes and rates.
Small-scale variations in mineral chemistry, textures, and platinum group element (PGE) mineralization were investigated in the Lower and Middle Group chromitite layers LG6, LG6a, MG1, MG2, and MG2 II from vertical drill core profiles at the Thaba mine in the northwestern limb of the Bushveld Complex. We present detailed geochemical profiles of chromite composition and chromite crystal size distribution curves to shed light on the processes of chromite accumulation and textural modification as well as mineralization. Multiple samples within each layer were assayed for PGE concentrations, and the respective platinum group mineral association was determined by mineral liberation analysis (MLA). There is strong evidence for postcumulus changes in the chromitites. The crystal size distribution curves suggest that the primary chromite texture was coarsened by a combination of adcumulus growth and textural equilibration, while compaction of the crystal mush played only a minor role. Mineral compositions were also modified by postcumulus processes, but because of the very high modal amount of chromite and its local preservation in orthopyroxene oikocrysts, that phase retained much primary information. Vertical variations of chromite composition within chromitite layers and from one layer to another do not support the idea of chromite accumulation from crystal-rich slurries or crystal settling from a large magma chamber. Instead, we favor a successive buildup of chromitite layers by repeated injections of relatively thin layers of chromite-saturated magmas, with in situ crystallization occurring at the crystal mush-magma interface. The adcumulus growth of chromite grains to form massive chromitite required addition of Cr to the layers, which we attribute to downward percolation from the overlying magma. The PGE concentrations are elevated in all chromitite layers compared to adjacent silicate rocks and show a systematic increase upward from LG6 (avg 807 ppb Ir + Ru + Rh + Pt + Pd + Au) to MG2 II (avg 2,062 ppb). There are also significant internal variations in all layers, with enrichments at hanging and/or footwalls. The enriched nature of chromitites in PGEs compared to host pyroxenites is a general feature, independent of the layer thickness. The MLA results distinguish two principal groups of PGE mineral associations: the LG6, LG6, and MG1 are dominated by the malanite series, laurite, and PGE sulfarsenides, while the MG2 and MG2 II layers are characterized by laurite and PGE sulfides as well as Pt-Fe-Sn and PGE-Sb-Bi-Pb alloys. Differences in the PGE associations are attributed to postcumulus alteration of the MG2 and MG2 II layer, while the chromitites below, particularly LG6 and LG6a, contain a more pristine association.
With the growing size and use of night light time series from the Visible Infrared Imaging Radiometer Suite Day/Night Band (DNB), it is important to understand the stability of the dataset. All satellites observe differences in pixel values during repeat observations. In the case of night light data, these changes can be due to both environmental effects and changes in light emission. Here we examine the stability of individual locations of particular large scale light sources (e.g., airports and prisons) in the monthly composites of DNB data from April 2012 to September 2017. The radiances for individual pixels of most large light emitters are approximately normally distributed, with a standard deviation of typically 15-20% of the mean. Greenhouses and flares, however, are not stable sources. We observe geospatial autocorrelation in the monthly variations for nearby sites, while the correlation for sites separated by large distances is small. This suggests that local factors contribute most to the variation in the pixel radiances and furthermore that averaging radiances over large areas will reduce the total variation. A better understanding of the causes of temporal variation would improve the sensitivity of DNB to lighting changes.
Variation of deuterium excess in surface waters across a 5000-m elevation gradient in eastern Nepal
(2020)
The strong elevation gradient of the Himalaya allows for investigation of altitude and orographic impacts on surface water delta O-18 and delta D stable isotope values. This study differentiates the time- and altitude-variable contributions of source waters to the Arun River in eastern Nepal. It provides isotope data along a 5000-m gradient collected from tributaries as well as groundwater, snow, and glacial-sourced surface waters and time-series data from April to October 2016. We find nonlinear trends in delta O-18 and delta D lapse rates with high-elevation lapse rates (4000-6000 masl) 5-7 times more negative than low-elevation lapse rates (1000-3000 masl). A distinct seasonal signal in delta O-18 and delta D lapse rates indicates time-variable source-water contributions from glacial and snow meltwater as well as precipitation transitions between the Indian Summer Monsoon and Winter Westerly Disturbances. Deuterium excess correlates with the extent of snowpack and tracks melt events during the Indian Summer Monsoon season. Our analysis identifies the influence of snow and glacial melt waters on river composition during low-flow conditions before the monsoon (April/May 2016) followed by a 5-week transition to the Indian Summer Monsoon-sourced rainfall around mid-June 2016. In the post-monsoon season, we find continued influence from glacial melt waters as well as ISM-sourced groundwater.
The southern margin of the Central Anatolian Plateau (CAP) records a strong uplift phase after the early Middle Pleistocene, which has been related to the slab break-off of the subducting Arabian plate beneath the Anatolian microplate. During the last 450 kyr the area underwent an uplift phase at a mean rate of similar to 3.2 m/kyr, as suggested by Middle Pleistocene marine sediments exposed at similar to 1,500 m above sea level. These values are significantly higher than the 1.0-1.5 m/kyr estimated since the Late Pleistocene, suggesting temporal variations in uplift rate. To estimate changes in uplift rate during the Pleistocene we studied the marine terraces along the CAP southern margin, mapping the remnants of the platforms and their associated deposits in the field, and used the TerraceM software to identify the position and elevation of associated shoreline angles. We used shoreline angles and the timing of Quaternary marine sedimentation as constrains for a Landscape Evolution Model that simulates wave erosion of an uplifting coast. We applied random optimization algorithms and minimization statistics to find the input parameters that better reproduce the morphology of CAP marine terraces. The best-fitting uplift rate history suggests a significative increase from 1.9 to 3.5 m/kyr between 500 and 200 kyr, followed by an abrupt decrease to 1.4 m/kyr until the present. Our results agree with slab break-off models, which suggest a strong uplift pulse during slab rupture followed by a smoother decrease.
Abstract. The aim of this study is to investigate the shallow thermal field differences for two differently aged passive continental margins by analyzing regional variations in geothermal gradient and exploring the controlling factors for these variations. Hence, we analyzed two previously published 3-D conductive and lithospheric-scale thermal models of the Southwest African and the Norwegian passive margins. These 3-D models differentiate various sedimentary, crustal, and mantle units and integrate different geophysical data such as seismic observations and the gravity field. We extracted the temperature–depth distributions in 1 km intervals down to 6 km below the upper thermal boundary condition. The geothermal gradient was then calculated for these intervals between the upper thermal boundary condition and the respective depth levels (1, 2, 3, 4, 5, and 6 km below the upper thermal boundary condition). According to our results, the geothermal gradient decreases with increasing depth and shows varying lateral trends and values for these two different margins. We compare the 3-D geological structural models and the geothermal gradient variations for both thermal models and show how radiogenic heat production, sediment insulating effect, and thermal lithosphere–asthenosphere boundary (LAB) depth influence the shallow thermal field pattern. The results indicate an ongoing process of oceanic mantle cooling at the young Norwegian margin compared with the old SW African passive margin that seems to be thermally equilibrated in the present day.
Abstract. The aim of this study is to investigate the shallow thermal field differences for two differently aged passive continental margins by analyzing regional variations in geothermal gradient and exploring the controlling factors for these variations. Hence, we analyzed two previously published 3-D conductive and lithospheric-scale thermal models of the Southwest African and the Norwegian passive margins. These 3-D models differentiate various sedimentary, crustal, and mantle units and integrate different geophysical data such as seismic observations and the gravity field. We extracted the temperature–depth distributions in 1 km intervals down to 6 km below the upper thermal boundary condition. The geothermal gradient was then calculated for these intervals between the upper thermal boundary condition and the respective depth levels (1, 2, 3, 4, 5, and 6 km below the upper thermal boundary condition). According to our results, the geothermal gradient decreases with increasing depth and shows varying lateral trends and values for these two different margins. We compare the 3-D geological structural models and the geothermal gradient variations for both thermal models and show how radiogenic heat production, sediment insulating effect, and thermal lithosphere–asthenosphere boundary (LAB) depth influence the shallow thermal field pattern. The results indicate an ongoing process of oceanic mantle cooling at the young Norwegian margin compared with the old SW African passive margin that seems to be thermally equilibrated in the present day.
Variability of the Cold Season Climate in Central Asia. Part II: Hydroclimatic Predictability
(2019)
Central Asia (CA) is subjected to a large variability of precipitation. This study presents a statistical model, relating precipitation anomalies in three subregions of CA in the cold season (November-March) with various predictors in the preceding October. Promising forecast skill is achieved for two subregions covering 1) Uzbekistan, Turkmenistan, Kyrgyzstan, Tajikistan, and southern Kazakhstan and 2) Iran, Afghanistan, and Pakistan. ENSO in October is identified as the major predictor. Eurasian snow cover and the quasi-biennial oscillation further improve the forecast performance. To understand the physical mechanisms, an analysis of teleconnections between these predictors and the wintertime circulation over CA is conducted. The correlation analysis of predictors and large-scale circulation indices suggests a seasonal persistence of tropical circulation modes and a dynamical forcing of the westerly circulation by snow cover variations over Eurasia. An EOF analysis of pressure and humidity patterns allows separating the circulation variability over CA into westerly and tropical modes and confirms that the identified predictors affect the respective circulation characteristics. Based on the previously established weather type classification for CA, the predictors are investigated with regard to their effect on the regional circulation. The results suggest a modification of the Hadley cell due to ENSO variations, with enhanced moisture supply from the Arabian Gulf during El Nino. They further indicate an influence of Eurasian snow cover on the wintertime Arctic Oscillation (AO) and Northern Hemispheric Rossby wave tracks. Positive anomalies favor weather types associated with dry conditions, while negative anomalies promote the formation of a quasi-stationary trough over CA, which typically occurs during positive AO conditions.
A reconstruction of Milankovitch to millennial-scale variability of sea surface temperature (SST) and sea surface productivity in the Pleistocene midlatitude North Atlantic Ocean (marine isotope stage (MIS) 16-9) and its relationship to ice sheet instability was carried out on sediments from Integrated Ocean Drilling Program (IODP) Site U1313. This reconstruction is based on alkenone and n-alkane concentrations, U-37(K)' index, total organic carbon (TOC) and carbonate contents, X-ray diffraction data, magnetic susceptibility, and accumulation rates. Increased input of ice-rafted debris occurred during MIS 16, 12, and 10, characterized by high concentrations of dolomite, quartz, and feldspars and elevated accumulation rates of terrigenous matter. Minimum input values of terrigenous matter, on the other hand, were determined for MIS 13 and 11. Peak values of dolomite, coinciding with quartz, plagioclase, and kalifeldspar peaks and maxima in long-chain n-alkanes indicative for land plants, are interpreted as Heinrich-like events related to sudden instability of the Laurentide Ice Sheet during early and late (deglacial) phases of the glacials. The coincidence of increased TOC values with elevated absolute concentrations of alkenones suggests increased glacial productivity, probably due to a more southern position of the Polar Front. Alkenone-based SST reached absolute maxima of about 19 degrees C during MIS 11.3 and absolute minima of <10 degrees C during MIS 12 and 10. Within MIS 11, prominent cooling events (MIS 11.22 and 11.24) occurred. The absolute SST minima recorded directly before and after the glacial maxima MIS 10.2 and 12.2 are related to Heinrich-like event meltwater pulses, as supported by the coincidence of SST minima and maxima in C-37:4 alkenones and dolomite. These sudden meltwater pulses, especially during terminations IV and V, probably caused a collapse of phytoplankton productivity as indicated by the distinct drop in alkenone concentrations. Ice sheet disintegration and subsequent surges and outbursts of icebergs and meltwater discharge may have been triggered by increased insolation in the northern high latitudes.
To better understand the reaction of Arctic coasts to increasing environmental pressure, coastal changes along a 210-km length of the Yukon Territory coast in north-west Canada were investigated. Shoreline positions were acquired from aerial and satellite images between 1951 and 2011. Shoreline change rates were calculated for multiple time periods along the entire coast and at six key sites. Additionally, Differential Global Positioning System (DGPS) measurements of shoreline positions from seven field sites were used to analyze coastal dynamics from 1991 to 2015 at higher spatial resolution. The whole coast has a consistent, spatially averaged mean rate of shoreline change of 0.7 +/- 0.2 m/a with a general trend of decreasing erosion from west to east. Additional data from six key sites shows that the mean shoreline change rate decreased from -1.3 +/- 0.8 (1950s-1970s) to -0.5 +/- 0.6 m/a (1970s-1990s). This was followed by a significant increase in shoreline change to -1.3 +/- 0.3 m/a in the 1990s to 2011. This increase is confirmed by DGPS measurements that indicate increased erosion rates at local rates up to -8.9 m/a since 2006. Ground surveys and observations with remote sensing data indicate that the current rate of shoreline retreat along some parts of the Yukon coast is higher than at any time before in the 64-year-long observation record. Enhanced availability of material in turn might favor the buildup of gravel features, which have been growing in extent throughout the last six decades. Plain Language Summary The Arctic is warming, but the impacts on its coasts are not well documented. To better understand the reaction of Arctic coasts to increasing environmental pressure, shoreline position changes along a 210-km length of the Yukon Territory coast in northwest Canada were investigated for the time period from 1951 to 2015. Shoreline positions were extracted from historical aerial images from the 1950s, 1970s, and 1990s and from satellite images from 2011. Additionally, measurements of shoreline positions from field sites were used to analyze coastal dynamics from 1991 to 2015. The mean shoreline change rate was -1.3 m/a between the 1950s and 1970s and followed by a decrease to -0.5 m/a between the 1970s to 1990s. This was followed by a significant increase in mean shoreline change rates again to -1.3 m/a in the 1990s to 2011 time period. This acceleration in erosion is confirmed by field measurements that indicate increased erosion rates at high local rates up to -8.9 m/a since 2006. Enhanced coastal erosion might, in turn, favor the buildup of gravel features, which have been growing in extent throughout the last six decades.
In this study, we validate and compare elevation accuracy and geomorphic metrics of satellite-derived digital elevation models (DEMs) on the southern Central Andean Plateau. The plateau has an average elevation of 3.7 km and is characterized by diverse topography and relief, lack of vegetation, and clear skies that create ideal conditions for remote sensing. At 30m resolution, SRTM-C, ASTER GDEM2, stacked ASTER L1A stereopair DEM, ALOS World 3D, and TanDEM-X have been analyzed. The higher-resolution datasets include 12m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X DEMs, and 5m ALOS World 3D. These DEMs are state of the art for optical (ASTER and ALOS) and radar (SRTM-C and TanDEM-X) spaceborne sensors. We assessed vertical accuracy by comparing standard deviations of the DEM elevation versus 307 509 differential GPS measurements across 4000m of elevation. For the 30m DEMs, the ASTER datasets had the highest vertical standard deviation at > 6.5 m, whereas the SRTM-C, ALOS World 3D, and TanDEM-X were all < 3.5 m. Higher-resolution DEMs generally had lower uncertainty, with both the 12m TanDEM-X and 5m ALOSWorld 3D having < 2m vertical standard deviation. Analysis of vertical uncertainty with respect to terrain elevation, slope, and aspect revealed the low uncertainty across these attributes for SRTM-C (30 m), TanDEM-X (12–30 m), and ALOS World 3D (5–30 m). Single-CoSSC TerraSAR-X/TanDEM-X 10m DEMs and the 30m ASTER GDEM2 displayed slight aspect biases, which were removed in their stacked counterparts (TanDEM-X and ASTER Stack). Based on low vertical standard deviations and visual inspection alongside optical satellite data, we selected the 30m SRTM-C, 12–30m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X, and 5m ALOS World 3D for geomorphic metric comparison in a 66 km2 catchment with a distinct river knickpoint. Consistent m=n values were found using chi plot channel profile analysis, regardless of DEM type and spatial resolution. Slope, curvature, and drainage area were calculated and plotting schemes were used to assess basin-wide differences in the hillslope-to-valley transition related to the knickpoint. While slope and hillslope length measurements vary little between datasets, curvature displays higher magnitude measurements with fining resolution. This is especially true for the optical 5m ALOS World 3D DEM, which demonstrated high-frequency noise in 2–8 pixel steps through a Fourier frequency analysis. The improvements in accurate space-radar DEMs (e.g., TanDEM-X) for geomorphometry are promising, but airborne or terrestrial data are still necessary for meter-scale analysis.
In this study, we validate and compare elevation accuracy and geomorphic metrics of satellite-derived digital elevation models (DEMs) on the southern Central Andean Plateau. The plateau has an average elevation of 3.7 km and is characterized by diverse topography and relief, lack of vegetation, and clear skies that create ideal conditions for remote sensing. At 30m resolution, SRTM-C, ASTER GDEM2, stacked ASTER L1A stereopair DEM, ALOS World 3D, and TanDEM-X have been analyzed. The higher-resolution datasets include 12m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X DEMs, and 5m ALOS World 3D. These DEMs are state of the art for optical (ASTER and ALOS) and radar (SRTM-C and TanDEM-X) spaceborne sensors. We assessed vertical accuracy by comparing standard deviations of the DEM elevation versus 307 509 differential GPS measurements across 4000m of elevation. For the 30m DEMs, the ASTER datasets had the highest vertical standard deviation at > 6.5 m, whereas the SRTM-C, ALOS World 3D, and TanDEM-X were all < 3.5 m. Higher-resolution DEMs generally had lower uncertainty, with both the 12m TanDEM-X and 5m ALOSWorld 3D having < 2m vertical standard deviation. Analysis of vertical uncertainty with respect to terrain elevation, slope, and aspect revealed the low uncertainty across these attributes for SRTM-C (30 m), TanDEM-X (12–30 m), and ALOS World 3D (5–30 m). Single-CoSSC TerraSAR-X/TanDEM-X 10m DEMs and the 30m ASTER GDEM2 displayed slight aspect biases, which were removed in their stacked counterparts (TanDEM-X and ASTER Stack). Based on low vertical standard deviations and visual inspection alongside optical satellite data, we selected the 30m SRTM-C, 12–30m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X, and 5m ALOS World 3D for geomorphic metric comparison in a 66 km2 catchment with a distinct river knickpoint. Consistent m=n values were found using chi plot channel profile analysis, regardless of DEM type and spatial resolution. Slope, curvature, and drainage area were calculated and plotting schemes were used to assess basin-wide differences in the hillslope-to-valley transition related to the knickpoint. While slope and hillslope length measurements vary little between datasets, curvature displays higher magnitude measurements with fining resolution. This is especially true for the optical 5m ALOS World 3D DEM, which demonstrated high-frequency noise in 2–8 pixel steps through a Fourier frequency analysis. The improvements in accurate space-radar DEMs (e.g., TanDEM-X) for geomorphometry are promising, but airborne or terrestrial data are still necessary for meter-scale analysis.
In this study, we validate and compare elevation accuracy and geomorphic metrics of satellite-derived digital elevation models (DEMs) on the southern Central Andean Plateau. The plateau has an average elevation of 3.7 km and is characterized by diverse topography and relief, lack of vegetation, and clear skies that create ideal conditions for remote sensing. At 30m resolution, SRTM-C, ASTER GDEM2, stacked ASTER L1A stereopair DEM, ALOS World 3D, and TanDEM-X have been analyzed. The higher-resolution datasets include 12m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X DEMs, and 5m ALOS World 3D. These DEMs are state of the art for optical (ASTER and ALOS) and radar (SRTM-C and TanDEM-X) spaceborne sensors. We assessed vertical accuracy by comparing standard deviations of the DEM elevation versus 307 509 differential GPS measurements across 4000m of elevation. For the 30m DEMs, the ASTER datasets had the highest vertical standard deviation at > 6.5 m, whereas the SRTM-C, ALOS World 3D, and TanDEM-X were all < 3.5 m. Higher-resolution DEMs generally had lower uncertainty, with both the 12m TanDEM-X and 5m ALOSWorld 3D having < 2m vertical standard deviation. Analysis of vertical uncertainty with respect to terrain elevation, slope, and aspect revealed the low uncertainty across these attributes for SRTM-C (30 m), TanDEM-X (12-30 m), and ALOS World 3D (5-30 m). Single-CoSSC TerraSAR-X/TanDEM-X 10m DEMs and the 30m ASTER GDEM2 displayed slight aspect biases, which were removed in their stacked counterparts (TanDEM-X and ASTER Stack). Based on low vertical standard deviations and visual inspection alongside optical satellite data, we selected the 30m SRTM-C, 12-30m TanDEM-X, 10m single-CoSSC TerraSAR-X/TanDEM-X, and 5m ALOS World 3D for geomorphic metric comparison in a 66 km2 catchment with a distinct river knickpoint. Consistent m = n values were found using chi plot channel profile analysis, regardless of DEM type and spatial resolution. Slope, curvature, and drainage area were calculated and plotting schemes were used to assess basin-wide differences in the hillslope-to-valley transition related to the knickpoint. While slope and hillslope length measurements vary little between datasets, curvature displays higher magnitude measurements with fining resolution. This is especially true for the optical 5m ALOS World 3D DEM, which demonstrated high-frequency noise in 2-8 pixel steps through a Fourier frequency analysis. The improvements in accurate space-radar DEMs (e. g., TanDEM-X) for geomorphometry are promising, but airborne or terrestrial data are still necessary for meter-scale analysis.
The aim of this paper is to investigate the ability of various site-condition proxies (SCPs) to reduce ground-motion aleatory variability and evaluate how SCPs capture nonlinearity site effects. The SCPs used here are time-averaged shear-wave velocity in the top 30 m (V-S30), the topographical slope (slope), the fundamental resonance frequency (f(0)) and the depth beyond which V-s exceeds 800 m/s (H800). We considered first the performance of each SCP taken alone and then the combined performance of the 6 SCP pairs [V-S30-f(0)], [V-S30-H-800], [f(0)-slope], [H-800-slope], [V-S30-slope] and [f(0)-H-800]. This analysis is performed using a neural network approach including a random effect applied on a KiK-net subset for derivation of ground-motion prediction equations setting the relationship between various ground-motion parameters such as peak ground acceleration, peak ground velocity and pseudo-spectral acceleration PSA (T), and Mw, RJB, focal depth and SCPs. While the choice of SCP is found to have almost no impact on the median groundmotion prediction, it does impact the level of aleatory uncertainty. VS30 is found to perform the best of single proxies at short periods (T < 0.6 s), while f(0) and H-800 perform better at longer periods; considering SCP pairs leads to significant improvements, with particular emphasis on [V-S30-H-800] and [f(0)-slope] pairs. The results also indicate significant nonlinearity on the site terms for soft sites and that the most relevant loading parameter for characterising nonlinear site response is the "stiff" spectral ordinate at the considered period.
The apparent isotope enrichment factor epsilon(macrophyte) of submerged plants (epsilon(macrophyte-DIC) = delta C-13(macrophyte) - delta C-13(DIC)) is indicative of dissolved inorganic carbon (DIC) supply in neutral to alkaline waters and is related to variations in aquatic productivity (Papadimitriou et al. in Limnol Oceanogr 50:1084-1095, 2005). This paper aims to evaluate the usage of epsilon(macrophyte) inferred from isotopic analyses of submerged plant fossils in addition to analyses of lake carbonate as a palaeolimnological proxy for former HCO3 (-) concentrations. Stable carbon isotopic analysis of modern Potamogeton pectinatus leaves and its host water DIC from the Tibetan Plateau and Central Yakutia (Russia) yielded values between -23.3 and +0.4aEuro degrees and between +14.0 and +6.5aEuro degrees, respectively. Values of epsilon (Potamogeton-DIC) (range -15.4 to +1.1aEuro degrees) from these lakes are significantly correlated with host water HCO3 (-) concentration (range 78-2,200 mg/l) (r = -0.86; P < 0.001), thus allowing for the development of a transfer function. Palaeo-epsilon (Potamogeton-ostracods) values from Luanhaizi Lake on the NE Tibetan Plateau, as inferred from the stable carbon isotope measurement of fossil Potamogeton pectinatus seeds (range -24 to +2.8aEuro degrees) and ostracods (range -7.8 to +7.5%) range between -14.8 and 1.6aEuro degrees. Phases of assumed disequilibrium between delta C-13(DIC) and delta C-13(ostracods) known to occur in charophyte swards (as indicated by the deposition of charophyte fossils) were excluded from the analysis of palaeo-epsilon. The application of the epsilon (Potamogeton-DIC)-HCO3 (-) transfer function yielded a median palaeo-HCO3 (-) -concentration of 290 mg/l. Variations in the dissolved organic carbon supply compare well with aquatic plant productivity changes and lake level variability as inferred from a multiproxy study of the same record including analyses of plant macrofossils, ostracods, carbonate and organic content.
Major earthquakes ( M > 8) have repeatedly ruptured the Nazca-South America plate interface of south-central Chile involving meter scale land-level changes. Earthquake recurrence intervals, however, extending beyond limited historical records are virtually unknown, but would provide crucial data on the tectonic behavior of forearcs. We analyzed the spatiotemporal pattern of Holocene earthquakes on Santa Maria Island (SMI; 37 degrees S), located 20 km off the Chilean coast and approximately 70 km east of the trench. SMI hosts a minimum of 21 uplifted beach berms, of which a subset were dated to calculate a mean uplift rate of 2.3 +/- 0.2 m/ky and a tilting rate of 0.022 +/- 0.002 degrees/ky. The inferred recurrence interval of strandline-forming earthquakes is similar to 180 years. Combining coseismic uplift and aseismic subsidence during an earthquake cycle, the net gain in strandline elevation in this environment is similar to 0.4 m per event
Sedimentary proxies used to reconstruct marine productivity suffer from variable preservation and are sensitive to factors other than productivity. Therefore, proxy calibration is warranted. Here we map the spatial patterns of two paleoproductivity proxies, biogenic opal and barium fluxes, from a set of core-top sediments recovered in the Subarctic North Pacific. Comparisons of the proxy data with independent estimates of primary and export production, surface water macronutrient concentrations, and biological pCO(2) drawdown indicate that neither proxy shows a significant correlation with primary or export productivity for the entire region. Biogenic opal fluxes, when corrected for preservation using Th-230-normalized accumulation rates, show a good correlation with primary productivity along the volcanic arcs (tau = 0.71, p = 0.0024) and with export productivity throughout the western Subarctic North Pacific (tau = 0.71, p = 0.0107). Moderate and good correlations of biogenic barium flux with export production (tau = 0.57, p = 0.0022) and with surface water silicate concentrations (tau = 0.70, p = 0.0002) are observed for the central and eastern Subarctic North Pacific. For reasons unknown, however, no correlation is found in the western Subarctic North Pacific between biogenic barium flux and the reference data. Nonetheless, we show that barite saturation, uncertainty in the lithogenic barium corrections, and problems with the reference data sets are not responsible for the lack of a significant correlation between biogenic barium flux and the reference data. Further studies evaluating the factors controlling the variability of the biogenic constituents in the sediments are desirable in this region.
High Mountain Asia provides water for more than a billion downstream users. Many catchments receive the majority of their yearly water budget in the form of snow - the vast majority of which is not monitored by sparse weather networks. We leverage passive microwave data from the SSMI series of satellites (SSMI, SSMI/S, 1987-2016), reprocessed to 3.125 km resolution, to examine trends in the volume and spatial distribution of snow-water equivalent (SWE) in the Indus Basin. We find that the majority of the Indus has seen an increase in snow-water storage. There exists a strong elevation-trend relationship, where high-elevation zones have more positive SWE trends. Negative trends are confined to the Himalayan foreland and deeply-incised valleys which run into the Upper Indus. This implies a temperature-dependent cutoff below which precipitation increases are not translated into increased SWE. Earlier snowmelt or a higher percentage of liquid precipitation could both explain this cutoff.(1) Earlier work 2 found a negative snow-water storage trend for the entire Indus catchment over the time period 1987-2009 (-4 x 10(-3) mm/yr). In this study based on an additional seven years of data, the average trend reverses to 1.4 x 10(-3). This implies that the decade since the mid-2000s was likely wetter, and positively impacted long-term SWE trends. This conclusion is supported by an analysis of snowmelt onset and end dates which found that while long-term trends are negative, more recent (since 2005) trends are positive (moving later in the year).(3)
Due to increasing demands and competition for high quality groundwater resources in many parts of the world, there is an urgent need for efficient methods that shed light on the interplay between complex natural settings and anthropogenic impacts. Thus a new approach is introduced, that aims to identify and quantify the predominant processes or factors of influence that drive groundwater and lake water dynamics on a catchment scale. The approach involves a non-linear dimension reduction method called Isometric feature mapping (Isomap). This method is applied to time series of groundwater head and lake water level data from a complex geological setting in Northeastern Germany. Two factors explaining more than 95% of the observed spatial variations are identified: (1) the anthropogenic impact of a waterworks in the study area and (2) natural groundwater recharge with different degrees of dampening at the respective sites of observation. The approach enables a presumption-free assessment to be made of the existing geological conception in the catchment, leading to an extension of the conception. Previously unknown hydraulic connections between two aquifers are identified, and connections revealed between surface water bodies and groundwater. (C) 2014 Elsevier B.V. All rights reserved.
Many institutions struggle to tap into the potential of their large archives of radar reflectivity: these data are often affected by miscalibration, yet the bias is typically unknown and temporally volatile. Still, relative calibration techniques can be used to correct the measurements a posteriori. For that purpose, the usage of spaceborne reflectivity observations from the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) platforms has become increasingly popular: the calibration bias of a ground radar (GR) is estimated from its average reflectivity difference to the spaceborne radar (SR). Recently, Crisologo et al. (2018) introduced a formal procedure to enhance the reliability of such estimates: each match between SR and GR observations is assigned a quality index, and the calibration bias is inferred as a quality-weighted average of the differences between SR and GR. The relevance of quality was exemplified for the Subic S-band radar in the Philippines, which is greatly affected by partial beam blockage.
The present study extends the concept of quality-weighted averaging by accounting for path-integrated attenuation (PIA) in addition to beam blockage. This extension becomes vital for radars that operate at the C or X band. Correspondingly, the study setup includes a C-band radar that substantially overlaps with the S-band radar. Based on the extended quality-weighting approach, we retrieve, for each of the two ground radars, a time series of calibration bias estimates from suitable SR overpasses. As a result of applying these estimates to correct the ground radar observations, the consistency between the ground radars in the region of overlap increased substantially. Furthermore, we investigated if the bias estimates can be interpolated in time, so that ground radar observations can be corrected even in the absence of prompt SR overpasses. We found that a moving average approach was most suitable for that purpose, although limited by the absence of explicit records of radar maintenance operations.
Many institutions struggle to tap into the potential of their large archives of radar reflectivity: these data are often affected by miscalibration, yet the bias is typically unknown and temporally volatile. Still, relative calibration techniques can be used to correct the measurements a posteriori. For that purpose, the usage of spaceborne reflectivity observations from the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) platforms has become increasingly popular: the calibration bias of a ground radar (GR) is estimated from its average reflectivity difference to the spaceborne radar (SR). Recently, Crisologo et al. (2018) introduced a formal procedure to enhance the reliability of such estimates: each match between SR and GR observations is assigned a quality index, and the calibration bias is inferred as a quality-weighted average of the differences between SR and GR. The relevance of quality was exemplified for the Subic S-band radar in the Philippines, which is greatly affected by partial beam blockage. The present study extends the concept of quality-weighted averaging by accounting for path-integrated attenuation (PIA) in addition to beam blockage. This extension becomes vital for radars that operate at the C or X band. Correspondingly, the study setup includes a C-band radar that substantially overlaps with the S-band radar. Based on the extended quality-weighting approach, we retrieve, for each of the two ground radars, a time series of calibration bias estimates from suitable SR overpasses. As a result of applying these estimates to correct the ground radar observations, the consistency between the ground radars in the region of overlap increased substantially. Furthermore, we investigated if the bias estimates can be interpolated in time, so that ground radar observations can be corrected even in the absence of prompt SR overpasses. We found that a moving average approach was most suitable for that purpose, although limited by the absence of explicit records of radar maintenance operations.
Many institutions struggle to tap into the potential of their large archives of radar reflectivity: these data are often affected by miscalibration, yet the bias is typically unknown and temporally volatile. Still, relative calibration techniques can be used to correct the measurements a posteriori. For that purpose, the usage of spaceborne reflectivity observations from the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) platforms has become increasingly popular: the calibration bias of a ground radar (GR) is estimated from its average reflectivity difference to the spaceborne radar (SR). Recently, Crisologo et al. (2018) introduced a formal procedure to enhance the reliability of such estimates: each match between SR and GR observations is assigned a quality index, and the calibration bias is inferred as a quality-weighted average of the differences between SR and GR. The relevance of quality was exemplified for the Subic S-band radar in the Philippines, which is greatly affected by partial beam blockage.
The present study extends the concept of quality-weighted averaging by accounting for path-integrated attenuation (PIA) in addition to beam blockage. This extension becomes vital for radars that operate at the C or X band. Correspondingly, the study setup includes a C-band radar that substantially overlaps with the S-band radar. Based on the extended quality-weighting approach, we retrieve, for each of the two ground radars, a time series of calibration bias estimates from suitable SR overpasses. As a result of applying these estimates to correct the ground radar observations, the consistency between the ground radars in the region of overlap increased substantially. Furthermore, we investigated if the bias estimates can be interpolated in time, so that ground radar observations can be corrected even in the absence of prompt SR overpasses. We found that a moving average approach was most suitable for that purpose, although limited by the absence of explicit records of radar maintenance operations.
The aims of this study were to identify areas of wind erosion and dust deposition and to quantify the effects of different grazing intensities on soil redistribution rates in grasslands based on the Cs-137 technique. Because the method uses a reference inventory as threshold for erosion or deposition, the classification of any other site as source or sink for dust depends on the accurate selection of this reference site.
Measurements of Cs-137 inventories and depth distributions were carried out at pasture sites with predominant species of Stipa grandis and Leymus chinensis which are grazed with different intensities. Additional measurements were made at arable land, plant-covered sand dunes and alluvial plains. Wind-induced soil erosion and dust deposition rates were calculated from Cs-137 inventories by means of the "Profile-Distribution" and the "Mass Balance II" models.
The selection of the reference site was based on fluid dynamical and process-determining parameters. The chosen site should meet the following four conditions: (i) located at a summit position with obviously low deposition rates, (ii) sufficient vegetation cover to prevent wind erosion, (iii) plane to exclude water erosion and (iv) in the wind/dust shadow of a higher elevation. The measured reference inventory of Cs-137 was 1967(+/- 102) Bqm(-2) located at a summit position of moderately grazed Leymus chinensis steppe. The Cs-137 inventories at other sites ranged from 1330 Bqm(-2) at heavily grazed sites to 5119 Bqm(-2) at river deposits, representing annual average soil losses of up to 130 tkm(-2) and deposits of up to 540 tkm(-2), respectively. The calculated annual averages of dust depositions at ungrazed Leymus chinensis sites were related to the dust storm frequencies of the last 50 years resulting in a description of the temporal variability of annual dust depositions from about 154 tkm(-2) in the 1960s to 26 tkm(-2) at recent times. Based on this quantification already 80% of the total dust depositions can be related to the 20 years between the 1960s and the end of the 1970s and only 20% to the time between 1980 and 2001.
Cs-137 technique is a promising method to assess the effect of grazing intensity and land use types on the spatial variability of wind-induced soil and dust redistribution processes in semi-arid grasslands. However, considerable efforts are needed to identify a reliable reference site, because erosion and deposition induced by wind may occur at the same places. The combination of the dust deposition rates derived from Cs-137 profile data with the dust storm frequencies is helpful for a better reconstruction of the temporal variability of dust deposition and wind erosion in this region. The calculated recent deposition rates of about 20 tkm(-2) are in good agreement with data of other authors.
The interactions between atmosphere and steep topography in the eastern south–central Andes result in complex relations with inhomogenous rainfall distributions. The atmospheric conditions leading to deep convection and extreme rainfall and their spatial patterns—both at the valley and mountain-belt scales—are not well understood. In this study, we aim to identify the dominant atmospheric conditions and their spatial variability by analyzing the convective available potential energy (CAPE) and dew-point temperature (Td). We explain the crucial effect of temperature on extreme rainfall generation along the steep climatic and topographic gradients in the NW Argentine Andes stretching from the low-elevation eastern foreland to the high-elevation central Andean Plateau in the west. Our analysis relies on version 2.0 of the ECMWF’s (European Centre for Medium-RangeWeather Forecasts) Re-Analysis (ERA-interim) data and TRMM (Tropical Rainfall Measuring Mission) data. We make the following key observations: First, we observe distinctive gradients along and across strike of the Andes in dew-point temperature and CAPE that both control rainfall distributions. Second, we identify a nonlinear correlation between rainfall and a combination of dew-point temperature and CAPE through a multivariable regression analysis. The correlation changes in space along the climatic and topographic gradients and helps to explain controlling factors for extreme-rainfall generation. Third, we observe more contribution (or higher importance) of Td in the tropical low-elevation foreland and intermediate-elevation areas as compared to the high-elevation central Andean Plateau for 90th percentile rainfall. In contrast, we observe a higher contribution of CAPE in the intermediate-elevation area between low and high elevation, especially in the transition zone between the tropical and subtropical areas for the 90th percentile rainfall. Fourth, we find that the parameters of the multivariable regression using CAPE and Td can explain rainfall with higher statistical significance for the 90th percentile compared to lower rainfall percentiles. Based on our results, the spatial pattern of rainfall-extreme events during the past ∼16 years can be described by a combination of dew-point temperature and CAPE in the south–central Andes.
The interactions between atmosphere and steep topography in the eastern south–central Andes result in complex relations with inhomogenous rainfall distributions. The atmospheric conditions leading to deep convection and extreme rainfall and their spatial patterns—both at the valley and mountain-belt scales—are not well understood. In this study, we aim to identify the dominant atmospheric conditions and their spatial variability by analyzing the convective available potential energy (CAPE) and dew-point temperature (Td). We explain the crucial effect of temperature on extreme rainfall generation along the steep climatic and topographic gradients in the NW Argentine Andes stretching from the low-elevation eastern foreland to the high-elevation central Andean Plateau in the west. Our analysis relies on version 2.0 of the ECMWF’s (European Centre for Medium-RangeWeather Forecasts) Re-Analysis (ERA-interim) data and TRMM (Tropical Rainfall Measuring Mission) data. We make the following key observations: First, we observe distinctive gradients along and across strike of the Andes in dew-point temperature and CAPE that both control rainfall distributions. Second, we identify a nonlinear correlation between rainfall and a combination of dew-point temperature and CAPE through a multivariable regression analysis. The correlation changes in space along the climatic and topographic gradients and helps to explain controlling factors for extreme-rainfall generation. Third, we observe more contribution (or higher importance) of Td in the tropical low-elevation foreland and intermediate-elevation areas as compared to the high-elevation central Andean Plateau for 90th percentile rainfall. In contrast, we observe a higher contribution of CAPE in the intermediate-elevation area between low and high elevation, especially in the transition zone between the tropical and subtropical areas for the 90th percentile rainfall. Fourth, we find that the parameters of the multivariable regression using CAPE and Td can explain rainfall with higher statistical significance for the 90th percentile compared to lower rainfall percentiles. Based on our results, the spatial pattern of rainfall-extreme events during the past ∼16 years can be described by a combination of dew-point temperature and CAPE in the south–central Andes.
Weltweit verursachen Städte etwa 70 % der Treibhausgasemissionen und sind daher wichtige Akteure im Klimaschutz bzw. eine wichtige Zielgruppe von Klimapolitiken. Gleichzeitig sind Städte besonders stark von möglichen Auswirkungen des Klimawandels betroffen: Insbesondere extreme Wetterereignisse wie Hitzewellen oder Starkregenereignisse mit Überflutungen verursachen in Städten hohe Sachschäden und wirken sich negativ auf die Gesundheit der städtischen Bevölkerung aus. Daher verfolgt das Projekt ExTrass das Ziel, die städtische Resilienz gegenüber extremen Wetterereignissen in enger Zusammenarbeit mit Stadtverwaltungen, Strukturen des Bevölkerungsschutzes und der Zivilgesellschaft zu stärken. Im Fokus stehen dabei (kreisfreie) Groß- und Mittelstädte mit 50.000 bis 500.000 Einwohnern, insbesondere die Fallstudienstädte Potsdam, Remscheid und Würzburg.
Der vorliegende Bericht beinhaltet die Ergebnisse der 14-monatigen Definitionsphase von ExTrass, in der vor allem die Abstimmung eines Arbeitsprogramms im Mittelpunkt stand, das in einem nachfolgenden dreijährigen Forschungsprojekt (F+E-Phase) gemeinsam von Wissenschaft und Praxispartnern umgesetzt werden soll. Begleitend wurde eine Bestandsaufnahme von Klimaanpassungs- und Klimaschutzstrategien/-plänen in 99 deutschen Groß- und Mittelstädten vorgenommen. Zudem wurden für Potsdam und Würzburg Pfadanalysen für die Klimapolitik durchgeführt. Darin wird insbesondere die Bedeutung von Schlüsselakteuren deutlich. Weiterhin wurden im Rahmen von Stakeholder-Workshops Anpassungsherausforderungen und aktuelle Handlungsbedarfe in den Fallstudienstädten identifiziert und Lösungsansätze erarbeitet, die in der F+E-Phase entwickelt und getestet werden sollen. Neben Maßnahmen auf gesamtstädtischer Ebene und auf Stadtteilebene wurden Maßnahmen angestrebt, die die Risikowahrnehmung, Vorsorge und Selbsthilfefähigkeit von Unternehmen und Bevölkerung stärken können. Daher wurde der Stand der Risikokommunikation in Deutschland für das Projekt aufgearbeitet und eine erste Evaluation von Risikokommunikationswerkzeugen durchgeführt. Der Bericht endet mit einer Kurzfassung des Arbeitsprogramms 2018-2021.
Cities will play a key role in the grand challenge of nourishing a growing global population, because, due to their population density, they set the demand. To ensure that food systems are sustainable, as well as nourishing, one solution often suggested is to shorten their supply chains toward a regional rather than a global basis. While such regional systems may have a range of costs and benefits, we investigate the mitigation potential of regionalized urban food systems by examining the greenhouse gas emissions associated with food transport. Using data on food consumption for 7108 urban administrative units (UAUs), we simulate total transport emissions for both regionalized and globalized supply chains. In regionalized systems, the UAUs' demands are fulfilled by peripheral food production, whereas to simulate global supply chains, food demand is met from an international pool (where the origin can be any location globally). We estimate that regionalized systems could reduce current emissions from food transport. However, because longer supply chains benefit from maximizing comparative advantage, this emission reduction would require closing yield gaps, reducing food waste, shifting toward diversified farming, and consuming seasonal produce. Regionalization of food systems will be an essential component to limit global warming to well below 2 degrees C in the future.
The chapter illustrates how cities can be understood from a system–theory perspective as complex social systems. It argues that the classical and linear intervention methods are often no (longer) suitable for the complex structure, temporal dynamics and multifaceted processuality of urban development. It offers a systemic and systems theory-inspired method as an alternative approach, which allows for extended possibilities that are more appropriate for dealing with urban development processes. The method was developed on the basis of practical experience and theoretical insights. The approach should demonstrate for local decision-makers potential areas of activity for organising urban changes through co-design.
Humankind and their environment need to be protected from the harmful effects of spent nuclear fuel, and therefore disposal in deep geological formations is favoured worldwide. Suitability of potential host rocks is evaluated, among others, by the retention capacity with respect to radionuclides. Safety assessments are based on the quantification of radionuclide migration lengths with numerical simulations as experiments cannot cover the required temporal (1 Ma) and spatial scales (>100 m).
Aim of the present thesis is to assess the migration of uranium, a geochemically complex radionuclide, in the potential host rock Opalinus Clay. Radionuclide migration in clay formations is governed by diffusion due to their low permeability and retarded by sorption. Both processes highly depend on pore water geochemistry and mineralogy that vary between different facies. Diffusion is quantified with the single-component (SC) approach using one diffusion coefficient for all species and the process-based multi-component (MC) option. With this, each species is assigned its own diffusion coefficient and the interaction with the diffuse double layer is taken into account. Sorption is integrated via a bottom-up approach using mechanistic surface complexation models and cation exchange. Therefore, reactive transport simulations are conducted with the geochemical code PHREEQC to quantify uranium migration, i.e. diffusion and sorption, as a function of mineralogical and geochemical heterogeneities on the host rock scale.
Sorption processes are facies dependent. Migration lengths vary between the Opalinus Clay facies by up to 10 m. Thereby, the geochemistry of the pore water, in particular the partial pressure of carbon dioxide (pCO2), is more decisive for the sorption capacity than the amount of clay minerals. Nevertheless, higher clay mineral quantities compensate geochemical variations. Consequently, sorption processes must be quantified as a function of pore water geochemistry in contact with the mineral assemblage.
Uranium diffusion in the Opalinus Clay is facies independent. Speciation is dominated by aqueous ternary complexes of U(VI) with calcium and carbonate. Differences in the migration lengths between SC and MC diffusion are with +/-5 m negligible. Further, the application of the MC approach highly depends on the quality and availability of the underlying data. Therefore, diffusion processes can be adequately quantified with the SC approach using experimentally determined diffusion coefficients.
The hydrogeological system governs pore water geochemistry within the formation rather than the mineralogy. Diffusive exchange with the adjacent aquifers established geochemical gradients over geological time scales that can enhance migration by up to 25 m. Consequently, uranium sorption processes must be quantified following the identified priority: pCO2 > hydrogeology > mineralogy.
The presented research provides a workflow and orientation for other potential disposal sites with similar pore water geochemistry due to the identified mechanisms and dependencies. With a maximum migration length of 70 m, the retention capacity of the Opalinus Clay with respect to uranium is sufficient to fulfill the German legal minimum requirement of a thickness of at least 100 m.
Marked along-strike changes in stratigraphy, mountain belt morphology, basement exhumation, and deformation styles characterize the Andean retroarc; these changes have previously been related to spatiotemporal variations in the subduction angle. We modeled new apatite fission track and apatite (U-Th-Sm)/He data from nine ranges located between 26 degrees S and 28 degrees S. Using new and previously published data, we constructed a Cretaceous to Pliocene paleogeographic model that delineates a four-stage tectonic evolution: extensional tectonics during the Cretaceous (120-75 Ma), the formation of a broken foreland basin between 55 and 30 Ma, reheating due to burial beneath sedimentary rocks (18-13 Ma), and deformation, exhumation, and surface uplift during the Late Miocene and the Pliocene (13-3 Ma). Our model highlights how preexisting upper plate structures control the deformation patterns of broken foreland basins. Because retroarc deformation predates flat-slab subduction, we propose that slab anchoring may have been the precursor of Eocene-Oligocene compression in the Andean retroarc. Our model challenges models which consider broken foreland basins and retroarc deformation in the NW Argentinian Andes to be directly related to Miocene flat subduction.
Marked along-strike changes in stratigraphy, mountain belt morphology, basement exhumation, and deformation styles characterize the Andean retroarc; these changes have previously been related to spatiotemporal variations in the subduction angle. We modeled new apatite fission track and apatite (U-Th-Sm)/He data from nine ranges located between 26 degrees S and 28 degrees S. Using new and previously published data, we constructed a Cretaceous to Pliocene paleogeographic model that delineates a four-stage tectonic evolution: extensional tectonics during the Cretaceous (120-75 Ma), the formation of a broken foreland basin between 55 and 30 Ma, reheating due to burial beneath sedimentary rocks (18-13 Ma), and deformation, exhumation, and surface uplift during the Late Miocene and the Pliocene (13-3 Ma). Our model highlights how preexisting upper plate structures control the deformation patterns of broken foreland basins. Because retroarc deformation predates flat-slab subduction, we propose that slab anchoring may have been the precursor of Eocene-Oligocene compression in the Andean retroarc. Our model challenges models which consider broken foreland basins and retroarc deformation in the NW Argentinian Andes to be directly related to Miocene flat subduction.
Passive continental margins offer the unique opportunity to study the processes involved in continental extension and break-up. Within the LISPWAL (Llthospheric Structure of the Namibian continental Passive margin at the intersection with the Walvis Ridge from amphibious seismic investigations) project, combined on- and offshore seismic experiments were designed to characterize the Southern African passive margin at the Walvis Ridge in northern Namibia. In addition to extensive analysis of the crustal structures, we carried out seismic investigations targeting the velocity structure of the upper mantle in the landfall region of the Walvis Ridge with the Namibian coast. Upper mantle P-n travel time tomography from controlled source, amphibious seismic data was used to investigate the sub-Moho upper mantle seismic velocity. We succeeded in imaging upper mantle structures potentially associated with continental break-up and/or the Tristan da Cunha hotspot track. We found mostly coast-parallel sub-Moho velocity anomalies, interpreted as structures which were created during Gondwana break-up. (C)2016 Published by Elsevier B.V.
We present a high-resolution study of the upper mantle structure of Central Europe, including the western part of the East European Platform, based on S-receiver functions of 345 stations. A distinct contrast is found between Phanerozoic Europe and the East European Craton across the Trans-European Suture Zone. To the west, a pronounced velocity reduction with depth interpreted as lithosphere-asthenosphere boundary (LAB) is found at an average depth of 90 km. Beneath the craton, no strong and continuous LAB conversion is observed. Instead we find a distinct velocity reduction within the lithosphere, at 80-120 km depth. This mid-lithospheric discontinuity (MLD) is attributed to a compositional boundary between depleted and more fertile lithosphere created by late Proterozoic metasomatism. A potential LAB phase beneath the craton is very weak and varies in depth between 180 and 250 km, consistent with a reduced velocity contrast between the lower lithosphere and the asthenosphere. Within the Trans-European Suture Zone, lithospheric structure is characterized by strong heterogeneity. A dipping or step-wise increase to LAB depth of 150 km is imaged from Phanerozoic Europe to 20-22 degrees E, whereas no direct connection to the cratonic LAB or MLD to the east is apparent. At larger depths, a positive conversion associated with the lower boundary of the asthenosphere is imaged at 210-250 km depth beneath Phanerozoic Europe, continuing down to 300 km depth beneath the craton. Conversions from both 410 km and 660 km discontinuities are found at their nominal depth beneath Phanerozoic Europe, and the discontinuity at 410 km depth can also be traced into the craton. A potential negative conversion on top of the 410 km discontinuity found in migrated images is analyzed by modeling and attributed to interference with other converted phases.
We present the first image of the Madeira upper crustal structure, using ambient seismic noise tomography. 16 months of ambient noise, recorded in a dense network of 26 seismometers deployed across Madeira, allowed reconstructing Rayleigh wave Green's functions between receivers. Dispersion analysis was performed in the short period band from 1.0 to 4.0 s. Group velocity measurements were regionalized to obtain 20 tomographic images, with a lateral resolution of 2.0 km in central Madeira. Afterwards, the dispersion curves, extracted from each cell of the 2D group velocity maps, were inverted as a function of depth to obtain a 3D shear wave velocity model of the upper crust, from the surface to a depth of 2.0 km. The obtained 3D velocity model reveals features throughout the island that correlates well with surface geology and island evolution. (C) 2015 Elsevier B.V. All rights reserved.
Das Menderes Massiv im Westen der Türkei stellt eine große Kulmination metamorpher Gesteine dar. Das Untersuchungsgebiet ist im Zentralen Menderes Massiv (Ödemis Submassiv) gelegen, das von den beiden aktiven Gräben, dem Gediz Graben im Norden und dem Büyük Menderes Graben im Süden begrenzt wird. Die Untersuchungen der Eklogit Relikte im zentralen Menderes Massiv haben ergeben, dass sich im Menderes Massiv Hochdruckrelikte in unterschiedlichen tektonischen Positionen befinden. Zum einen existieren Eklogit-Blöcke in der obersten Einheit (Selcuk Einheit) des zentralen Menderes Massivs und zum anderen Hochdruck-Relikte in der strukturell mittleren Birgi - Tire Decke. Die Granate der quarzfreien Eklogit-Blöcke weisen große Ähnlichkeiten mit denen der HP/LT Gesteine von Sifnos und Syros auf. Die Entwicklung der Eklogit-Blöcke in der Olistostrom-Einheit lässt sich jedoch nicht mit den Eklogit Relikten in der strukturell mittleren Birgi Tire Decke vergleichen. Für die Eklogit-Relikte in der Birgi Tire Decke wurde eine polymetamorphe Entwicklung mithilfe petrologischer Untersuchungen und chemischen und Pb-Pb Datierungen herausgearbeitet. Die Eklogit Relikte gehören zu einem metamorphen Teilpfad, der durch eine Amphibolitfazies 1 - Hochdruck - Amphibolitfazies 2/Granulitfazies charakterisiert ist. Der Endpunkt dieses Teilpfades ist mit Temperaturen zwischen 700 und 750 °C und Drücken von 1.2 - 1.4 GPa belegt. Für diese Bedingungen konnte ein minimales Alter von 520 Ma durch chemische Datierungen an Monaziten einer Augengneisprobe und Pb-Pb Datierungen an Zirkonen einer Augengneis- und Metagabbroprobe bestimmt werden. Dieser amphibolit/granulitfazieller Endpunkt wird mit den Granitintrusionen des zentralen und südlichen Menderes Massiv korreliert, die in einem Zeitraum zwischen 520 Ma bis 550 Ma stattfanden. Sowohl die Amphibolitfazies 1 als auch das Hochdruckereignis werden der Panafrikanischen Orogenese zugeordnet. Für die Hochdruckbedingungen wurden maximale Temperaturen zwischen 680°C und 720°C und bei einem Druck von 2.2 GPa bestimmt. In den untersuchten Metasedimenten konnte eine prograde metamorphe Entwicklung abgeleitet werden, die amphibolitfazielle Bedingungen von 660°C bei 0.6 GPa erreichte. Das Metamorphosealter dieser Metasedimente konnte mit < 100 Ma mittels chemischer Mikrosondendatierung bestimmt werden. Die in den Metasedimenten herausgearbeiteten Druck- und Temperaturbedingungen wurden ebenfalls in den metabasischen Gesteinen bestimmt. Diese Ergebnisse werden als Krustenstapelung der metabasischen Gesteine, Augengneise und Metasedimente interpretiert, die mit der alpinen Orogenese im Zusammenhang stehen. Durch die Ergebnisse dieser Arbeit lässt sich die Birgi-Tire Decke im zentralen Menderes Massiv genauer charakterisieren. Sie besteht aus Metasedimenten, pelitischen Gneisen, Augengneisen und metabasichen Gesteinen. Die Gneise (pelitische und Augengneise) und die metabasischen Gesteine stellen panafrikanische Relikte dar, die einen amphibolit- eklogit- amphibolit/granulitfaziellen Metamorphosepfad gespeichert haben. Die amphibolit- bis granulitfazielle Metamorphose hängt mit den Granitintrusionen zusammen und fand in einem Zeitraum zwischen 520 - 550 Ma statt. Große Teile der Metasedimente der Birgi Tire Decke haben jedoch nur eine alpine metamorphe Entwicklung durchlaufen, wo sie unter amphibolitfazielle Bedingungen Krustentiefen erreichten, bei denen sie mit den panafrikanischen Relikten zusammen gestapelt wurden und eine gemeinsame Exhumierung erfahren haben.
This study presents an unsupervised feature selection and learning approach for the discovery and intuitive imaging of significant temporal patterns in seismic single-station or network recordings. For this purpose, the data are parametrized by real-valued feature vectors for short time windows using standard analysis tools for seismic data, such as frequency-wavenumber, polarization, and spectral analysis. We use Self-Organizing Maps (SOMs) for a data-driven feature selection, visualization and clustering procedure, which is in particular suitable for high-dimensional data sets. Our feature selection method is based on significance testing using the Wald-Wolfowitz runs test for-individual features and on correlation hunting with SOMs in feature subsets. Using synthetics composed of Rayleigh and Love waves and real-world data, we show the robustness and the improved discriminative power of that approach compared to feature subsets manually selected from individual wavefield parametrization methods. Furthermore, the capability of the clustering and visualization techniques to investigate the discrimination of wave phases is shown by means of synthetic waveforms and regional earthquake recordings.
Salt marshes filter pollutants, protect coastlines against storm surges, and sequester carbon, yet are under threat from sea level rise and anthropogenic modification. The sustained existence of the salt marsh ecosystem depends on the topographic evolution of marsh platforms. Quantifying marsh platform topography is vital for improving the management of these valuable landscapes. The determination of platform boundaries currently relies on supervised classification methods requiring near-infrared data to detect vegetation, or demands labour-intensive field surveys and digitisation. We propose a novel, unsupervised method to reproducibly isolate salt marsh scarps and platforms from a digital elevation model (DEM), referred to as Topographic Identification of Platforms (TIP). Field observations and numerical models show that salt marshes mature into subhorizontal platforms delineated by subvertical scarps. Based on this premise, we identify scarps as lines of local maxima on a slope raster, then fill landmasses from the scarps upward, thus isolating mature marsh platforms. We test the TIP method using lidar-derived DEMs from six salt marshes in England with varying tidal ranges and geometries, for which topographic platforms were manually isolated from tidal flats. Agreement between manual and unsupervised classification exceeds 94% for DEM resolutions of 1 m, with all but one site maintaining an accuracy superior to 90% for resolutions up to 3 m. For resolutions of 1 m, platforms detected with the TIP method are comparable in surface area to digitised platforms and have similar elevation distributions. We also find that our method allows for the accurate detection of local block failures as small as 3 times the DEM resolution. Detailed inspection reveals that although tidal creeks were digitised as part of the marsh platform, unsupervised classification categorises them as part of the tidal flat, causing an increase in false negatives and overall platform perimeter. This suggests our method may benefit from combination with existing creek detection algorithms. Fallen blocks and high tidal flat portions, associated with potential pioneer zones, can also lead to differences between our method and supervised mapping. Although pioneer zones prove difficult to classify using a topographic method, we suggest that these transition areas should be considered when analysing erosion and accretion processes, particularly in the case of incipient marsh platforms. Ultimately, we have shown that unsupervised classification of marsh platforms from high-resolution topography is possible and sufficient to monitor and analyse topographic evolution.
The Subandean fold and thrust belt of Bolivia constitutes the easternmost part of the Andean orogen that reflects thin-skinned shortening and eastward propagation of the Andean deformation front. The exact interplay of tectonics, climate, and erosion in the deposition of up to 7.5 km of late Cenozoic strata exposed in the Subandes remains unclear. To better constrain these relationships, we use four W-E industry seismic reflection profiles, eight new zircon U-Pb ages from Mio-Pliocene sedimentary strata, and cross-section balancing to evaluate the rates of thrust propagation, shortening, and deposition pinch-out migration. Eastward thrusting arrived in the Subandean belt at similar to 12.4 +/- 0.5 Ma and propagated rapidly toward the foreland unit approximately 6 Ma. This was followed by out-of- sequence deformation from ca. 4 to 2.1 Ma and by renewed eastward propagation thereafter. Our results show that the thrust-front propagation- and deposition pinch-out migration rates mimic the sediment accumulation rate. The rates of deposition pinchout migration and thrust propagation increased three- and two fold, respectively (8 mm/a; 3.3 mm/a) at 86 Ma. The three-fold increase in deposition pinch-out migration rate at this time is an indication of enhanced erosional efficiency in the hinterland, probably coupled with flexural rebound of the basin. Following the pulse of pinch-out migration, the Subandean belt witnessed rapid similar to 80 km eastward propagation of thrusting to the La Vertiente structure at 6 Ma. As there is no evidence for this event of thrust front migration being linked to an increase in shortening rate, the enhanced frontal accretion suggests a shift to supercritical wedge taper conditions. We propose that the supercritical state was due to a drop in basal strength, caused by sediment loading and pore fluid overpressure. This scenario implies that climate-controlled variation in erosional efficiency was the driver of late Miocene mass redistribution, which induced flexural rebound of the Subandean thrust belt, spreading of a large clastic wedge across the basin, and subsequent thrust-front propagation.
The Indus Molasse records orogenic sedimentation associated with uplift and erosion of the southern margin of Asia in the course of ongoing India-Eurasia collision. Detailed field investigation clarifies the nature and extent of the depositional contact between this molasse and the underlying basement units. We report the first dataset on detrital zircon U-Pb ages, Hf isotopes and apatite U-Pb ages for the autochthonous molasse in the Indus Suture Zone. A latest Oligocene depositional age is proposed on the basis of the youngest detrital zircon U-Pb age peak and is consistent with published biostratigraphic data. Multiple provenance indicators suggest exclusively northerly derivation with no input from India in the lowermost parts of the section. The results provide constraints on the uplift and erosion history of the Ladakh Range following the initial India-Asia collision.
The climate is a complex dynamical system involving interactions and feedbacks among different processes at multiple temporal and spatial scales. Although numerous studies have attempted to understand the climate system, nonetheless, the studies investigating the multiscale characteristics of the climate are scarce. Further, the present set of techniques are limited in their ability to unravel the multi-scale variability of the climate system. It is completely plausible that extreme events and abrupt transitions, which are of great interest to climate community, are resultant of interactions among processes operating at multi-scale. For instance, storms, weather patterns, seasonal irregularities such as El Niño, floods and droughts, and decades-long climate variations can be better understood and even predicted by quantifying their multi-scale dynamics. This makes a strong argument to unravel the interaction and patterns of climatic processes at different scales. With this background, the thesis aims at developing measures to understand and quantify multi-scale interactions within the climate system.
In the first part of the thesis, I proposed two new methods, viz, multi-scale event synchronization (MSES) and wavelet multi-scale correlation (WMC) to capture the scale-specific features present in the climatic processes. The proposed methods were tested on various synthetic and real-world time series in order to check their applicability and replicability. The results indicate that both methods (WMC and MSES) are able to capture scale-specific associations that exist between processes at different time scales in a more detailed manner as compared to the traditional single scale counterparts.
In the second part of the thesis, the proposed multi-scale similarity measures were used in constructing climate networks to investigate the evolution of spatial connections within climatic processes at multiple timescales. The proposed methods WMC and MSES, together with complex network were applied to two different datasets.
In the first application, climate networks based on WMC were constructed for the univariate global sea surface temperature (SST) data to identify and visualize the SSTs patterns that develop very similarly over time and distinguish them from those that have long-range teleconnections to other ocean regions. Further investigations of climate networks on different timescales revealed (i) various high variability and co-variability regions, and (ii) short and long-range teleconnection regions with varying spatial distance. The outcomes of the study not only re-confirmed the existing knowledge on the link between SST patterns like El Niño Southern Oscillation and the Pacific Decadal Oscillation, but also suggested new insights into the characteristics and origins of long-range teleconnections.
In the second application, I used the developed non-linear MSES similarity measure to quantify the multivariate teleconnections between extreme Indian precipitation and climatic patterns with the highest relevance for Indian sub-continent. The results confirmed significant non-linear influences that were not well captured by the traditional methods. Further, there was a substantial variation in the strength and nature of teleconnection across India, and across time scales.
Overall, the results from investigations conducted in the thesis strongly highlight the need for considering the multi-scale aspects in climatic processes, and the proposed methods provide robust framework for quantifying the multi-scale characteristics.
The architecture of coastal sequences in tectonically-active regions results mostly from a combination of sea-level and land-level changes. The objective of this study is to unravel these signals by combining sequence stratigraphy and sedimentology of near-shore sedimentary sequences in wave-built terraces. We focus on Santa Maria Island at the south-central Chile margin, which hosts excellent exposures of coastal sediments from Marine Isotope Stage 3. A novel method based on statistical analysis of grain-size distributions coupled with fades descriptions provided a detailed account of transgressive-regressive cycles. Radiocarbon ages from paleosols constrain the chronology between >53 and similar to 31 cal ka BP. Because the influence of glaciations can be neglected, we calculated relative sea-level curves by tying the onset of deposition on a bedrock abrasion platform to a global sea-level curve. The observed depositional cycles match those predicted for uplift rates between 1.2 and 1.8 m/ka. The studied sedimentary units represent depositional cycles that resulted in reoccupation events of an existing marine terrace. Our study demonstrates wave-built marine terrace deposits along clastic shorelines in temperate regions can be used to distinguish between tectonic uplift and climate-induced sea-level changes.
Reductions in emissions have successfully led to a regional decline in atmospheric nitrogen depositions over the past 20 years. By analyzing long-term data from 110 mountainous streams draining into German drinking water reservoirs, nitrate concentrations indeed declined in the majority of catchments. Furthermore, our meta-analysis indicates that the declining nitrate levels are linked to the release of dissolved iron to streams likely due to a reductive dissolution of iron(III) minerals in riparian wetland soils. This dissolution process mobilized adsorbed compounds, such as phosphate, dissolved organic carbon and arsenic, resulting in concentration increases in the streams and higher inputs to receiving drinking water reservoirs. Reductive mobilization was most significant in catchments with stream nitrate concentrations < 6 mg L-1. Here, nitrate, as a competing electron acceptor, was too low in concentration to inhibit microbial iron(III) reduction. Consequently, observed trends were strongest in forested catchments, where nitrate concentrations were unaffected by agricultural and urban sources and which were therefore sensitive to reductions of atmospheric nitrogen depositions. We conclude that there is strong evidence that the decline in nitrogen deposition toward pre-industrial conditions lowers the redox buffer in riparian soils, destabilizing formerly fixed problematic compounds, and results in serious implications for water quality.
Due to the challenges in upscaling daily climatic forcing to geological time, physically realistic models describing how rainfall drives fluvial erosion are lacking. To bridge this gap between short-term hydrology and long-term geomorphology, we derive a theoretical framework for long-term fluvial erosion rates driven by realistic climate by integrating an established stochastic-mechanistic model of hydrology into a threshold-stochastic formulation of stream power. The hydrological theory provides equations for the daily streamflow probability distribution as a function of climatic boundary conditions. The new parameters introduced are rooted firmly in established climatic and hydrological theory. This allows us to account for how fluvial erosion rates respond to changes in rainfall intensity, frequency, evapotranspiration rates, and soil moisture dynamics in a way that is consistent with existing theories. We use this framework to demonstrate how hydroclimatic conditions and erosion threshold magnitude control the degree of nonlinearity between steepness index and erosion rate. We find that hydrological processes can have a significant influence on how erosive a particular climatic forcing will be. By accounting for the influence of hydrology on fluvial erosion, we conclude that climate is an important control on erosion rates and long-term landscape evolution.
In recent years, nature-based solutions are receiving increasing attention in the field of disaster risk reduction and climate change adaptation as inclusive, no regret approaches. Ecosystem-based adaptation (EbA) can mitigate the impacts of climate change, build resilience and tackle environmental degradation thereby supporting the targets set by the 2030 Agenda, the Paris Agreement and the Sendai Framework. Despite these benefits, EbA is still rarely implemented in practice. To better understand the barriers to implementation, this research examines policy-makers' perceptions of EbA, using an extended version of Protection Motivation Theory as an analytical framework. Through semi-structured interviews with policy-makers at regional and provincial level in Central Vietnam, it was found that EbA is generally considered a promising response option, mainly due to its multiple ecosystem-service benefits. The demand for EbA measures was largely driven by the perceived consequences of natural hazards and climate change. Insufficient perceived response efficacy and time-lags in effectiveness for disaster risk reduction were identified as key impediments for implementation. Pilot projects and capacity building on EbA are important means to overcome these perceptual barriers. This paper contributes to bridging the knowledge-gap on political decision-making regarding EbA and can, thereby, promote its mainstreaming into policy plans.
Understanding the geologic evolution of Northern Tibetan Plateau with multiple thermochronometers
(2018)
The early onset of deformation following the India-Asia collision, Neogene expanse of uplift, and complex systems that comprise strike-slip faults, thrust faults, and intermontane basins characterize the Cenozoic tectonism of Northern Tibetan Plateau and raise two prominent questions in orogenic geodynamics: 1) What mechanism(s) control(s) the transfer of stress related to the India-Asia collision across the distance of >2000 km; and 2) Why the development of high topography was delayed in the Northern Tibetan Plateau and what does it reveal about how the internal forces and external boundary conditions evolved. To address these two questions, we reconstruct a holistic spatial-temporal deformation history of the Northern Tibetan Plateau by using a range of thermochronometers, with closure temperature spanning from 350 degrees C to-60-70 degrees C. This multi-thermochronometer study reveals three stages of faulting related cooling, in the early Cretaceous, in Paleocene-Eocene and in middle-late Miocene. We observe that Paleocene-Eocene deformation was spatially restricted and mostly occurred on reactivated Cretaceous structures, indicating a control of pre-existing weakness on early Cenozoic deformation. Extensive Neogene deformation contrasts with restricted Paleocene-Eocene deformation and relatively quiescent shortening during the Oligocene-early Miocene, which implies a change in the regional tectonics regime. Global plate reconstructions show that this tectonic reorganization is coeval with an increase in Pacific-Asia plate convergence rates. We argue that this change in regional tectonics is a result of increasing constrictive environment of the eastern plate boundary, which changed the behavior of the Altyn Tagh fault the boundary fault of Northern Tibetan Plateau, causing it to change from feeding slip into structures out of the plateau to feeding slip into structures at plateau margins.
Accelerometric data from the well-studied valley EUROSEISTEST are used to investigate ground motion uncertainty and variability. We define a simple local ground motion prediction equation (GMPE) and investigate changes in standard deviation (σ) and its components, the between-event variability (τ) and within-event variability (φ). Improving seismological metadata significantly reduces τ (30–50%), which in turn reduces the total σ. Improving site information reduces the systematic site-to-site variability, φ S2S (20–30%), in turn reducing φ, and ultimately, σ. Our values of standard deviations are lower than global values from literature, and closer to path-specific than site-specific values. However, our data have insufficient azimuthal coverage for single-path analysis. Certain stations have higher ground-motion variability, possibly due to topography, basin edge or downgoing wave effects. Sensitivity checks show that 3 recordings per event is a sufficient data selection criterion, however, one of the dataset’s advantages is the large number of recordings per station (9–90) that yields good site term estimates. We examine uncertainty components binning our data with magnitude from 0.01 to 2 s; at smaller magnitudes, τ decreases and φ SS increases, possibly due to κ and source-site trade-offs Finally, we investigate the alternative approach of computing φ SS using existing GMPEs instead of creating an ad hoc local GMPE. This is important where data are insufficient to create one, or when site-specific PSHA is performed. We show that global GMPEs may still capture φ SS , provided that: (1) the magnitude scaling errors are accommodated by the event terms; (2) there are no distance scaling errors (use of a regionally applicable model). Site terms (φ S2S ) computed by different global GMPEs (using different site-proxies) vary significantly, especially for hard-rock sites. This indicates that GMPEs may be poorly constrained where they are sometimes most needed, i.e., for hard rock.
Groundwater is critical in supporting current and future reliable water supply throughout Africa. Although continental maps of groundwater storage and recharge have been developed, we currently lack a clear understanding on how the controls on groundwater recharge vary across the entire continent. Reviewing the existing literature, we synthesize information on reported groundwater recharge controls in Africa. We find that 15 out of 22 of these controls can be characterised using global datasets. We develop 11 descriptors of climatic, topographic, vegetation, soil and geologic properties using global datasets, to characterise groundwater recharge controls in Africa. These descriptors cluster Africa into 15 Recharge Landscape Units for which we expect recharge controls to be similar. Over 80% of the continents land area is organized by just nine of these units. We also find that aggregating the Units by similarity into four broader Recharge Landscapes (Desert, Dryland, Wet tropical and Wet tropical forest) provides a suitable level of landscape organisation to explain differences in ground-based long-term mean annual recharge and recharge ratio (annual recharge / annual precipitation) estimates. Furthermore, wetter Recharge Landscapes are more efficient in converting rainfall to recharge than drier Recharge Landscapes as well as having higher annual recharge rates. In Dryland Recharge Landscapes, we found that annual recharge rates largely varied according to mean annual precipitation, whereas recharge ratio estimates increase with increasing monthly variability in P-PET. However, we were unable to explain why ground based estimates of recharge signatures vary across other Recharge Landscapes, in which there are fewer ground based recharge estimates, using global datasets alone. Even in dryland regions, there is still considerable unexplained variability in the estimates of annual recharge and recharge ratio, stressing the limitations of global datasets for investigating ground-based information.
Understanding the rates and processes of denudation is key to unraveling the dynamic processes that shape active orogens. This includes decoding the roles of tectonic and climate-driven processes in the long-term evolution of high- mountain landscapes in regions with pronounced tectonic activity and steep climatic and surface-process gradients. Well-constrained denudation rates can be used to address a wide range of geologic problems. In steady-state landscapes, denudation rates are argued to be proportional to tectonic or isostatic uplift rates and provide valuable insight into the tectonic regimes underlying surface denudation. The use of denudation rates based on terrestrial cosmogenic nuclide (TCN) such as 10Beryllium has become a widely-used method to quantify catchment-mean denudation rates. Because such measurements are averaged over timescales of 102 to 105 years, they are not as susceptible to stochastic changes as shorter-term denudation rate estimates (e.g., from suspended sediment measurements) and are therefore considered more reliable for a comparison to long-term processes that operate on geologic timescales. However, the impact of various climatic, biotic, and surface processes on 10Be concentrations and the resultant denudation rates remains unclear and is subject to ongoing discussion. In this thesis, I explore the interaction of climate, the biosphere, topography, and geology in forcing and modulating denudation rates on catchment to orogen scales.
There are many processes in highly dynamic active orogens that may effect 10Be concentrations in modern river sands and therefore impact 10Be-derived denudation rates. The calculation of denudation rates from 10Be concentrations, however, requires a suite of simplifying assumptions that may not be valid or applicable in many orogens. I investigate how these processes affect 10Be concentrations in the Arun Valley of Eastern Nepal using 34 new 10Be measurements from the main stem Arun River and its tributaries. The Arun Valley is characterized by steep gradients in climate and topography, with elevations ranging from <100 m asl in the foreland basin to >8,000 asl in the high sectors to the north. This is coupled with a five-fold increase in mean annual rainfall across strike of the orogen. Denudation rates from tributary samples increase toward the core of the orogen, from <0.2 to >5 mm/yr from the Lesser to Higher Himalaya. Very high denudation rates (>2 mm/yr), however, are likely the result of 10Be TCN dilution by surface and climatic processes, such as large landsliding and glaciation, and thus may not be representative of long-term denudation rates. Mainstem Arun denudation rates increase downstream from ~0.2 mm/yr at the border with Tibet to 0.91 mm/yr at its outlet into the Sapt Kosi. However, the downstream 10Be concentrations may not be representative of the entire upstream catchment. Instead, I document evidence for downstream fining of grains from the Tibetan Plateau, resulting in an order-of-magnitude apparent decrease in the measured 10Be concentration.
In the Arun Valley and across the Himalaya, topography, climate, and vegetation are strongly interrelated. The observed increase in denudation rates at the transition from the Lesser to Higher Himalaya corresponds to abrupt increases in elevation, hillslope gradient, and mean annual rainfall. Thus, across strike (N-S), it is difficult to decipher the potential impacts of climate and vegetation cover on denudation rates. To further evaluate these relationships I instead took advantage of an along-strike west-to-east increase of mean annual rainfall and vegetation density in the Himalaya. An analysis of 136 published 10Be denudation rates from along strike of the revealed that median denudation rates do not vary considerably along strike of the Himalaya, ~1500 km E-W. However, the range of denudation rates generally decreases from west to east, with more variable denudation rates in the northwestern regions of the orogen than in the eastern regions. This denudation rate variability decreases as vegetation density increases (R=- 0.90), and increases proportionately to the annual seasonality of vegetation (R=0.99). Moreover, rainfall and vegetation modulate the relationship between topographic steepness and denudation rates such that in the wet, densely vegetated regions of the Himalaya, topography responds more linearly to changes in denudation rates than in dry, sparsely vegetated regions, where the response of topographic steepness to denudation rates is highly nonlinear. Understanding the relationships between denudation rates, topography, and climate is also critical for interpreting sedimentary archives. However, there is a lack of understanding of how terrestrial organic matter is transported out of orogens and into sedimentary archives. Plant wax lipid biomarkers derived from terrestrial and marine sedimentary records are commonly used as paleo- hydrologic proxy to help elucidate these problems. I address the issue of how to interpret the biomarker record by using the plant wax isotopic composition of modern suspended and riverbank organic matter to identify and quantify organic matter source regions in the Arun Valley. Topographic and geomorphic analysis, provided by the 10Be catchment-mean denudation rates, reveals that a combination of topographic steepness (as a proxy for denudation) and vegetation density is required to capture organic matter sourcing in the Arun River.
My studies highlight the importance of a rigorous and careful interpretation of denudation rates in tectonically active orogens that are furthermore characterized by strong climatic and biotic gradients. Unambiguous information about these issues is critical for correctly decoding and interpreting the possible tectonic and climatic forces that drive erosion and denudation, and the manifestation of the erosion products in sedimentary archives.
Statistical distributions of flood peak discharge often show heavy tail behavior, that is, extreme floods are more likely to occur than would be predicted by commonly used distributions that have exponential asymptotic behavior.
This heavy tail behavior may surprise flood managers and citizens, as human intuition tends to expect light tail behavior, and the heaviness of the tails is very difficult to predict, which may lead to unnecessarily high flood damage.
Despite its high importance, the literature on the heavy tail behavior of flood distributions is rather fragmented.
In this review, we provide a coherent overview of the processes causing heavy flood tails and the implications for science and practice.
Specifically, we propose nine hypotheses on the mechanisms causing heavy tails in flood peak distributions related to processes in the atmosphere, the catchment, and the river system.
We then discuss to which extent the current knowledge supports or contradicts these hypotheses.
We also discuss the statistical conditions for the emergence of heavy tail behavior based on derived distribution theory and relate them to the hypotheses and flood generation mechanisms.
We review the degree to which the heaviness of the tails can be predicted from process knowledge and data. Finally, we recommend further research toward testing the hypotheses and improving the prediction of heavy tails.
Understanding the rates and pattern of erosion is a key aspect of deciphering the impacts of climate and tectonics on landscape evolution. Denudation rates derived from terrestrial cosmogenic nuclides (TCNs) are commonly used to quantify erosion and bridge tectonic (Myr) and climatic (up to several kiloyears) time scales. However, how the processes of erosion in active orogens are ultimately reflected in Be-10 TCN samples remains a topic of discussion. We investigate this problem in the Arun Valley of eastern Nepal with 34 new Be-10-derived catchment-mean denudation rates. The Arun Valley is characterized by steep north-south gradients in topography and climate. Locally, denudation rates increase northward, from <0.2mmyr(-1) to similar to 1.5mmyr(-1) in tributary samples, while main stem samples appear to increase downstream from similar to 0.2mmyr(-1) at the border with Tibet to 0.91mmyr(-1) in the foreland. Denudation rates most strongly correlate with normalized channel steepness (R-2=0.67), which has been commonly interpreted to indicate tectonic activity. Significant downstream decrease of Be-10 concentration in the main stem Arun suggests that upstream sediment grains are fining to the point that they are operationally excluded from the processed sample. This results in Be-10 concentrations and denudation rates that do not uniformly represent the upstream catchment area. We observe strong impacts on Be-10 concentrations from local, nonfluvial geomorphic processes, such as glaciation and landsliding coinciding with areas of peak rainfall rates, pointing toward climatic modulation of predominantly tectonically driven denudation rates.
Beta diversity is a conceptual link between diversity at local and regional scales. Various additional methodologies of quantifying this and related phenomena have been applied. Among them, measures of pairwise (dis)similarity of sites are particularly popular. Undersampling, i.e. not recording all taxa present at a site, is a common situation in ecological data. Bias in many metrics related to beta diversity must be expected, but only few studies have explicitly investigated the properties of various measures under undersampling conditions. On the basis of an empirical data set, representing near-complete local inventories of the Lepidoptera from an isolated Pacific island, as well as simulated communities with varying properties, we mimicked different levels of undersampling. We used 14 different approaches to quantify beta diversity, among them dataset-wide multiplicative partitioning (i.e. true beta diversity') and pairwise site x site dissimilarities. We compared their values from incomplete samples to true results from the full data. We used these comparisons to quantify undersampling bias and we calculated correlations of the dissimilarity measures of undersampled data with complete data of sites. Almost all tested metrics showed bias and low correlations under moderate to severe undersampling conditions (as well as deteriorating precision, i.e. large chance effects on results). Measures that used only species incidence were very sensitive to undersampling, while abundance-based metrics with high dependency on the distribution of the most common taxa were particularly robust. Simulated data showed sensitivity of results to the abundance distribution, confirming that data sets of high evenness and/or the application of metrics that are strongly affected by rare species are particularly sensitive to undersampling. The class of beta measure to be used should depend on the research question being asked as different metrics can lead to quite different conclusions even without undersampling effects. For each class of metric, there is a trade-off between robustness to undersampling and sensitivity to rare species. In consequence, using incidence-based metrics carries a particular risk of false conclusions when undersampled data are involved. Developing bias corrections for such metrics would be desirable.
Saturated hydraulic conductivity (K-s) is an important soil characteristic affecting soil water storage, runoff generation and erosion processes. In some areas where high-intensity rainfall coincides with low K-s values at shallow soil depths, frequent overland flow entails dense drainage networks. Consequently, linear structures such as flowlines alternate with inter-flowline areas. So far, investigations of the spatial variability of K-s mainly relied on isotropic covariance models which are unsuitable to reveal patterns resulting from linear structures. In the present study, we applied two sampling approaches so as to adequately characterize K-s spatial variability in a tropical forest catchment that features a high density of flowlines: A classical nested sampling survey and a purposive sampling strategy adapted to the presence of flowlines. The nested sampling approach revealed the dominance of small-scale variability, which is in line with previous findings. Our purposive sampling, however, detected a strong spatial gradient: surface K-s increased substantially as a function of distance to flowline; 10 m off flowlines, values were similar to the spatial mean of K-s. This deterministic trend can be included as a fixed effect in a linear mixed modeling framework to obtain realistic spatial fields of K-s. In a next step we used probability maps based on those fields and prevailing rainfall intensities to assess the hydrological relevance of the detected pattern. This approach suggests a particularly good agreement between the probability statements of K-s exceedance and observed overland flow occurrence during wet stages of the rainy season.
The Late Permian Zechstein Group in northeastern Germany is characterized by shelf and slope carbonates that rimmed a basin extending from eastern England through the Netherlands and Germany to Poland. Conventional reservoirs are found in grainstones rimming islands created by pre-existing paleohighs and platform-rimming shoals that compose steep margins in the north and ramp deposits in the southern part. The slope and basin deposits are characterized by debris flows and organic-rich mudstones. Lagoonal and basinal evaporites formed the seal for these carbonate and underlying sandstone reservoirs. The objective of this investigation is to evaluate potential unconventional reservoirs in organic-rich, fine-grained and/or tight mudrocks in slope and basin as well as platform carbonates occurring in this stratigraphic interval. Therefore, a comprehensive study was conducted that included sedimentology, sequence stratigraphy, petrography, and geochemistry. Sequence stratigraphic correlations from shelf to basin are crucial in establishing a framework that allows correlation of potential productive facies in fine-grained, organic-rich basinal siliceous and calcareous mudstones or interfingering tight carbonates and siltstones, ranging from the lagoon, to slope to basin, which might be candidates for forming an unconventional reservoir. Most organic-rich shales worldwide are associated with eustatic transgressions. The basal Zechstein cycles, Z1 and Z2, contain organic-rich siliceous and calcareous mudstones and carbonates that form major transgressive deposits in the basin. Maturities range from over-mature (gas) in the basin to oil-generation on the slope with variable TOC contents. This sequence stratigraphic and sedimentologic evaluation of the transgressive facies in the Z1 and Z2 assesses the potential for shale-gas/oil and hybrid unconventional plays. Potential unconventional reservoirs might be explored in laminated organic-rich mudstones within the oil window along the northern and southern slopes of the basin. Although the Zechstein Z1 and Z2 cycles might have limited shale-gas potential because of low thickness and deep burial depth to be economic at this point, unconventional reservoir opportunities that include hybrid and shale-oil potential are possible in the study area.
Shrinking glaciers in the Hindu Kush-Karakoram-Himalaya-Nyainqentanglha (HKKHN) region have formed several thousand moraine-dammed glacial lakes(1-3), some of these having grown rapidly in past decades(3,4). This growth may promote more frequent and potentially destructive glacial lake outburst floods (GLOFs)(5-7). Testing this hypothesis, however, is confounded by incomplete databases of the few reliable, though selective, case studies. Here we present a consistent Himalayan GLOF inventory derived automatically from all available Landsat imagery since the late 1980s. We more than double the known GLOF count and identify the southern Himalayas as a hotspot region, compared to the more rarely affected Hindu Kush-Karakoram ranges. Nevertheless, the average annual frequency of 1.3 GLOFs has no credible posterior trend despite reported increases in glacial lake areas in most of the HKKHN3,8, so that GLOF activity per unit lake area has decreased since the late 1980s. We conclude that learning more about the frequency and magnitude of outburst triggers, rather than focusing solely on rapidly growing glacial lakes, might improve the appraisal of GLOF hazards.
Cosmic-ray neutron sensing (CRNS) is a promising proximal soil sensing technique to estimate soil moisture at intermediate scale and high temporal resolution. However, the signal shows complex and non-unique response to all hydrogen pools near the land surface, providing some challenges for soil moisture estimation in practical applications. Aims of the study were 1) to assess the uncertainty of CRNS as a stand-alone approach to estimate volumetric soil moisture in cropped field 2) to identify the causes of this uncertainty 3) and possible improvements. Two experimental sites in Germany were equipped with a CRNS probe and point-scale soil moisture network. Additional monitoring activities were conducted during the crop growing season to characterize the soil-plant systems. This data is used to identify and quantify the different sources of uncertainty (factors). An uncertainty analysis, based on Monte Carlo approach, is applied to propagate these uncertainties to CRNS soil moisture estimations. In addition, a sensitivity analysis based on the Sobol’ method is performed to identify the most important factors explaining this uncertainty. Results show that CRNS soil moisture compares well to the soil moisture network when these point-scale values are weighted to account for the spatial sensitivity of the signal and other sources of hydrogen (lattice water and organic carbon) are added to the water content. However, the performance decreases when CRNS is considered as a stand-alone method to retrieve the actual (non-weighted) volumetric soil moisture. The support volume (penetration depth and radius) shows also a considerable uncertainty, especially in relatively dry soil moisture conditions. Four of the seven factors analyzed (the vertical soil moisture profile, bulk density, incoming neutron correction and the calibrated parameter N0) were found to play an important role. Among the possible improvements identified, a simple correction factor based on vertical point-scale soil moisture profiles shows to be a promising approach to account for the sensitivity of the CRNS signal to the upper soil layers.
The selection of earthquake focal mechanisms (FMs) for stress tensor inversion (STI) is commonly done on a spatial basis, that is, hypocentres. However, this selection approach may include data that are undesired, for example, by mixing events that are caused by different stress tensors when for the STI a single stress tensor is assumed. Due to the significant increase of FM data in the past decades, objective data-driven data selection is feasible, allowing more refined FM catalogues that avoid these issues and provide data weights for the STI routines. We present the application of angular classification with expectation-maximization (ACE) as a tool for data selection. ACE identifies clusters of FM without a priori information. The identified clusters can be used for the classification of the style-of-faulting and as weights of the FM data. We demonstrate that ACE effectively selects data that can be associated with a single stress tensor. Two application examples are given for weighted STI from South America. We use the resulting clusters and weights as a priori information for an STI for these regions and show that uncertainties of the stress tensor estimates are reduced significantly.
Himalayan water resources attract a rapidly growing number of hydroelectric power projects (HPP) to satisfy Asia's soaring energy demands. Yet HPP operating or planned in steep, glacier-fed mountain rivers face hazards of glacial lake outburst floods (GLOFs) that can damage hydropower infrastructure, alter water and sediment yields, and compromise livelihoods downstream. Detailed appraisals of such GLOF hazards are limited to case studies, however, and a more comprehensive, systematic analysis remains elusive. To this end we estimate the regional exposure of 257 Himalayan HPP to GLOFs, using a flood-wave propagation model fed by Monte Carlo-derived outburst volumes of >2300 glacial lakes. We interpret the spread of thus modeled peak discharges as a predictive uncertainty that arises mainly from outburst volumes and dam-breach rates that are difficult to assess before dams fail. With 66% of sampled HPP are on potential GLOF tracks, up to one third of these HPP could experience GLOF discharges well above local design floods, as hydropower development continues to seek higher sites closer to glacial lakes. We compute that this systematic push of HPP into headwaters effectively doubles the uncertainty about GLOF peak discharge in these locations. Peak discharges farther downstream, in contrast, are easier to predict because GLOF waves attenuate rapidly. Considering this systematic pattern of regional GLOF exposure might aid the site selection of future Himalayan HPP. Our method can augment, and help to regularly update, current hazard assessments, given that global warming is likely changing the number and size of Himalayan meltwater lakes.
Himalayan water resources attract a rapidly growing number of hydroelectric power projects (HPP) to satisfy Asia's soaring energy demands. Yet HPP operating or planned in steep, glacier-fed mountain rivers face hazards of glacial lake outburst floods (GLOFs) that can damage hydropower infrastructure, alter water and sediment yields, and compromise livelihoods downstream. Detailed appraisals of such GLOF hazards are limited to case studies, however, and a more comprehensive, systematic analysis remains elusive. To this end we estimate the regional exposure of 257 Himalayan HPP to GLOFs, using a flood-wave propagation model fed by Monte Carlo-derived outburst volumes of >2300 glacial lakes. We interpret the spread of thus modeled peak discharges as a predictive uncertainty that arises mainly from outburst volumes and dam-breach rates that are difficult to assess before dams fail. With 66% of sampled HPP are on potential GLOF tracks, up to one third of these HPP could experience GLOF discharges well above local design floods, as hydropower development continues to seek higher sites closer to glacial lakes. We compute that this systematic push of HPP into headwaters effectively doubles the uncertainty about GLOF peak discharge in these locations. Peak discharges farther downstream, in contrast, are easier to predict because GLOF waves attenuate rapidly. Considering this systematic pattern of regional GLOF exposure might aid the site selection of future Himalayan HPP. Our method can augment, and help to regularly update, current hazard assessments, given that global warming is likely changing the number and size of Himalayan meltwater lakes.
Himalayan water resources attract a rapidly growing number of hydroelectric power projects (HPP) to satisfy Asia's soaring energy demands. Yet HPP operating or planned in steep, glacier-fed mountain rivers face hazards of glacial lake outburst floods (GLOFs) that can damage hydropower infrastructure, alter water and sediment yields, and compromise livelihoods downstream. Detailed appraisals of such GLOF hazards are limited to case studies, however, and a more comprehensive, systematic analysis remains elusive. To this end we estimate the regional exposure of 257 Himalayan HPP to GLOFs, using a flood-wave propagation model fed by Monte Carlo-derived outburst volumes of >2300 glacial lakes. We interpret the spread of thus modeled peak discharges as a predictive uncertainty that arises mainly from outburst volumes and dam-breach rates that are difficult to assess before dams fail. With 66% of sampled HPP are on potential GLOF tracks, up to one third of these HPP could experience GLOF discharges well above local design floods, as hydropower development continues to seek higher sites closer to glacial lakes. We compute that this systematic push of HPP into headwaters effectively doubles the uncertainty about GLOF peak discharge in these locations. Peak discharges farther downstream, in contrast, are easier to predict because GLOF waves attenuate rapidly. Considering this systematic pattern of regional GLOF exposure might aid the site selection of future Himalayan HPP. Our method can augment, and help to regularly update, current hazard assessments, given that global warming is likely changing the number and size of Himalayan meltwater lakes.
Empirical species distribution models (SDMs) constitute often the tool of choice for the assessment of rapid climate change effects on species vulnerability. Conclusions regarding extinction risks might be misleading, however, because SDMs do not explicitly incorporate dispersal or other demographic processes. Here, we supplement SDMs with a dynamic population model 1) to predict climate-induced range dynamics for black grouse in Switzerland, 2) to compare direct and indirect measures of extinction risks, and 3) to quantify uncertainty in predictions as well as the sources of that uncertainty. To this end, we linked models of habitat suitability to a spatially explicit, individual-based model. In an extensive sensitivity analysis, we quantified uncertainty in various model outputs introduced by different SDM algorithms, by different climate scenarios and by demographic model parameters. Potentially suitable habitats were predicted to shift uphill and eastwards. By the end of the 21st century, abrupt habitat losses were predicted in the western Prealps for some climate scenarios. In contrast, population size and occupied area were primarily controlled by currently negative population growth and gradually declined from the beginning of the century across all climate scenarios and SDM algorithms. However, predictions of population dynamic features were highly variable across simulations. Results indicate that inferring extinction probabilities simply from the quantity of suitable habitat may underestimate extinction risks because this may ignore important interactions between life history traits and available habitat. Also, in dynamic range predictions uncertainty in SDM algorithms and climate scenarios can become secondary to uncertainty in dynamic model components. Our study emphasises the need for principal evaluation tools like sensitivity analysis in order to assess uncertainty and robustness in dynamic range predictions. A more direct benefit of such robustness analysis is an improved mechanistic understanding of dynamic species responses to climate change.
Tillage erosion on arable land is a very important process leading to a net downslope movement of soil and soil constitutes. Tillage erosion rates are commonly in the same order of magnitude as water erosion rates and can be even higher, especially under highly mechanized agricultural soil management. Despite its prevalence and magnitude, tillage erosion is still understudied compared to water erosion. The goal of this study was to bring together experts using different techniques to determine tillage erosion and use the different results to discuss and quantify uncertainties associated with tillage erosion measurements. The study was performed in northeastern Germany on a 10 m by 50 m plot with a mean slope of 8%. Tillage erosion was determined after two sequences of seven tillage operations. Two different micro-tracers (magnetic iron oxide mixed with soil and fluorescent sand) and one macro-tracer (passive radio-frequency identification transponders (RFIDs), size: 4 x 22 mm) were used to directly determine soil fluxes. Moreover, tillage induced changes in topography were measured for the entire plot with two different terrestrial laser scanners and an unmanned aerial system for structure from motion topography analysis. Based on these elevation differences, corresponding soil fluxes were calculated. The mean translocation distance of all techniques was 0.57 m per tillage pass, with a relatively wide range of mean soil translocation distances ranging from 039 to 0.72 m per pass. A benchmark technique could not be identified as all used techniques have individual error sources, which could not be quantified. However, the translocation distances of the macro-tracers used were consistently smaller than the translocation distances of the micro-tracers (mean difference = -26 +/- 12%), which questions the widely used assumption of non-selective soil transport via tillage operations. This study points out that tillage erosion measurements, carried out under almost optimal conditions, are subject to major uncertainties that are far from negligible. (C) 2018 Elsevier B.V. All rights reserved.
The advantages of remote sensing using Unmanned Aerial Vehicles (UAVs) are a high spatial resolution of images, temporal flexibility and narrow-band spectral data from different wavelengths domains. This enables the detection of spatio-temporal dynamics of environmental variables, like plant-related carbon dynamics in agricultural landscapes. In this paper, we quantify spatial patterns of fresh phytomass and related carbon (C) export using imagery captured by a 12-band multispectral camera mounted on the fixed wing UAV Carolo P360. The study was performed in 2014 at the experimental area CarboZALF-D in NE Germany. From radiometrically corrected and calibrated images of lucerne (Medicago sativa), the performance of four commonly used vegetation indices (VIs) was tested using band combinations of six near-infrared bands. The highest correlation between ground-based measurements of fresh phytomass of lucerne and VIs was obtained for the Enhanced Vegetation Index (EVI) using near-infrared band b(899). The resulting map was transformed into dry phytomass and finally upscaled to total C export by harvest. The observed spatial variability at field- and plot-scale could be attributed to small-scale soil heterogeneity in part.
The advantages of remote sensing using Unmanned Aerial Vehicles (UAVs) are a high spatial resolution of images, temporal flexibility and narrow-band spectral data from different wavelengths domains. This enables the detection of spatio-temporal dynamics of environmental variables, like plant-related carbon dynamics in agricultural landscapes. In this paper, we quantify spatial patterns of fresh phytomass and related carbon (C) export using imagery captured by a 12-band multispectral camera mounted on the fixed wing UAV Carolo P360. The study was performed in 2014 at the experimental area CarboZALF-D in NE Germany. From radiometrically corrected and calibrated images of lucerne (Medicago sativa), the performance of four commonly used vegetation indices (VIs) was tested using band combinations of six near-infrared bands. The highest correlation between ground-based measurements of fresh phytomass of lucerne and VIs was obtained for the Enhanced Vegetation Index (EVI) using near-infrared band b(899). The resulting map was transformed into dry phytomass and finally upscaled to total C export by harvest. The observed spatial variability at field- and plot-scale could be attributed to small-scale soil heterogeneity in part.
The age of Proterozoic granulite facies metamorphism and deformation in the Strangways Metamorphic Complex (SMC) of central Australia is determined on zircon grown in syn-metamorphic and syn-deformational orthopyroxene-bearing, enderbitic, veins. SHRIMP zircon studies suggest that M1-M2 and the correlated periods of intense deformation (D1-D2) are part of a single tectonothermal event between 1,717-2 and 1,732-7 Ma. It is considered unlikely that the two metamorphic phases (M1, M2) suggested by earlier work represent separate events occurring within 10-25 Ma of each other. Previous higher estimates for the age of M1 granulite metamorphism in the SMC (Early Strangways event at ca. 1,770 Ma) based on U-Pb zircon dating of granitic, intrusive rocks, are not believed to relate to the metamorphism, but to represent pre-metamorphic intrusion ages. Conventional multi-grain U-Pb monazite analyses on high-grade metasediments from three widely spaced localities in the western SMC yield 207Pb/235U ages between 1,728-11 and 1,712-2 Ma. The age range of the monazites corresponds to the SHRIMP zircon ages in the granulitic veins and is interpreted to record monazite growth (prograde in the metasedimentary rocks). The data imply a maximum time-span of 30 Ma for high-grade metamorphism and deformation in the SMC. There is, thus, no evidence for an extremely long period of continuous high- temperature conditions from 1,770 to ca. 1,720 Ma as previously proposed. The results firmly establish that the SMC has a very different high-grade metamorphic history than the neighbouring Harts Range, where upper amphibolite facies metamorphism in the Palaeozoic caused widespread growth or recrystallization of monazite.
The type-locality granulites from the Granulitgebirge of Saxony, Germany, are rocks of broadly granitic composition containing minor garnet and kyanite within a commonly mylonitised matrix of feldspars and quartz. Petrographic evidence indicates a primary assemblage of ternary feldspar + quartz + garnet + kyanite + rutile, most likely resulting from partial melting of a granitic protolith, for which equilibrium temperature and pressure conditions of > 1000 degrees C and > 1.5 GPa have been deduced. These extreme (for crustal rocks) conditions, and the inferred peak assemblage, are supported by the newly-developed Zr-in-rutile geothermometer and experimental studies on the same bulk composition, respectively. As these conditions lie above those required for plagioclase stability in quartz tholeiites, they are thus in the eclogite facies. Widespread modification of the peak assemblage, for example mesoperthite formation after ternary feldspar, deformation-induced recrystallisation of perthites to two-feldspar + quartz aggregates, biotite replacing garnet, Ca-loss at garnet rims, sillimanite replacing kyanite or secondary garnet growth, makes reliable interpretation of equilibrium assemblages and compositions very difficult and explains the spread of published pressure- temperature values and consequent confusion about formation depths and the validity of tectonometamorphic models. Such extreme metamorphic conditions in rock compositions typical for the upper continental crust, reflecting a hot subduction environment, has important consequences for understanding some collisional orogens
Tectonic and geological processes on Earth often result in structural anisotropy of the subsurface, which can be imaged by various geophysical methods. In order to achieve appropriate and realistic Earth models for interpretation, inversion algorithms have to allow for an anisotropic subsurface. Within the framework of this thesis, I analyzed a magnetotelluric (MT) data set taken from the Cape Fold Belt in South Africa. This data set exhibited strong indications for crustal anisotropy, e.g. MT phases out of the expected quadrant, which are beyond of fitting and interpreting with standard isotropic inversion algorithms. To overcome this obstacle, I have developed a two-dimensional inversion method for reconstructing anisotropic electrical conductivity distributions. The MT inverse problem represents in general a non-linear and ill-posed minimization problem with many degrees of freedom: In isotropic case, we have to assign an electrical conductivity value to each cell of a large grid to assimilate the Earth's subsurface, e.g. a grid with 100 x 50 cells results in 5000 unknown model parameters in an isotropic case; in contrast, we have the sixfold in an anisotropic scenario where the single value of electrical conductivity becomes a symmetric, real-valued tensor while the number of the data remains unchanged. In order to successfully invert for anisotropic conductivities and to overcome the non-uniqueness of the solution of the inverse problem it is necessary to use appropriate constraints on the class of allowed models. This becomes even more important as MT data is not equally sensitive to all anisotropic parameters. In this thesis, I have developed an algorithm through which the solution of the anisotropic inversion problem is calculated by minimization of a global penalty functional consisting of three entries: the data misfit, the model roughness constraint and the anisotropy constraint. For comparison, in an isotropic approach only the first two entries are minimized. The newly defined anisotropy term is measured by the sum of the square difference of the principal conductivity values of the model. The basic idea of this constraint is straightforward. If an isotropic model is already adequate to explain the data, there is no need to introduce electrical anisotropy at all. In order to ensure successful inversion, appropriate trade-off parameters, also known as regularization parameters, have to be chosen for the different model constraints. Synthetic tests show that using fixed trade-off parameters usually causes the inversion to end up by either a smooth model with large RMS error or a rough model with small RMS error. Using of a relaxation approach on the regularization parameters after each successful inversion iteration will result in smoother inversion model and a better convergence. This approach seems to be a sophisticated way for the selection of trade-off parameters. In general, the proposed inversion method is adequate for resolving the principal conductivities defined in horizontal plane. Once none of the principal directions of the anisotropic structure is coincided with the predefined strike direction, only the corresponding effective conductivities, which is the projection of the principal conductivities onto the model coordinate axes direction, can be resolved and the information about the rotation angles is lost. In the end the MT data from the Cape Fold Belt in South Africa has been analyzed. The MT data exhibits an area (> 10 km) where MT phases over 90 degrees occur. This part of data cannot be modeled by standard isotropic modeling procedures and hence can not be properly interpreted. The proposed inversion method, however, could not reproduce the anomalous large phases as desired because of losing the information about rotation angles. MT phases outside the first quadrant are usually obtained by different anisotropic anomalies with oblique anisotropy strike. In order to achieve this challenge, the algorithm needs further developments. However, forward modeling studies with the MT data have shown that surface highly conductive heterogeneity in combination with a mid-crustal electrically anisotropic zone are required to fit the data. According to known geological and tectonic information the mid-crustal zone is interpreted as a deep aquifer related to the fractured Table Mountain Group rocks in the Cape Fold Belt.
Skarn deposits are found on every continents and were formed at different times from Precambrian to Tertiary. Typically, the formation of a skarn is induced by a granitic intrusion in carbonates-rich sedimentary rocks. During contact metamorphism, fluids derived from the granite interact with the sedimentary host rocks, which results in the formation of calc-silicate minerals at the expense of carbonates. Those newly formed minerals generally develop in a metamorphic zoned aureole with garnet in the proximal and pyroxene in the distal zone. Ore elements contained in magmatic fluids are precipitated due to the change in fluid composition. The temperature decrease of the entire system, due to the cooling of magmatic fluids and the entering of meteoric water, allows retrogression of some prograde minerals.
The Hämmerlein skarn deposit has a multi-stage history with a skarn formation during regional metamorphism and a retrogression of primary skarn minerals during the granitic intrusion. Tin was mobilized during both events. The 340 Ma old tin-bearing skarn minerals show that tin was present in sediments before the granite intrusion, and that the first Sn enrichment occurred during the skarn formation by regional metamorphism fluids. In a second step at ca. 320 Ma, tin-bearing fluids were produced with the intrusion of the Eibenstock granite. Tin, which has been added by the granite and remobilized from skarn calc-silicates, precipitated as cassiterite.
Compared to clay or marl, the skarn is enriched in Sn, W, In, Zn, and Cu. These metals have been supplied during both regional metamorphism and granite emplacement. In addition, the several isotopic and chemical data of skarn samples show that the granite selectively added elements such as Sn, and that there was no visible granitic contribution to the sedimentary signature of the skarn
The example of Hämmerlein shows that it is possible to form a tin-rich skarn without associated granite when tin has already been transported from tin-bearing sediments during regional metamorphism by aqueous metamorphic fluids. These skarns are economically not interesting if tin is only contained in the skarn minerals. Later alteration of the skarn (the heat and fluid source is not necessarily a granite), however, can lead to the formation of secondary cassiterite (SnO2), with which the skarn can become economically highly interesting.
It is commonly recognized that soil moisture exhibits spatial heterogeneities occurring in a wide range of scales. These heterogeneities are caused by different factors ranging from soil structure at the plot scale to land use at the landscape scale. There is an urgent need for effi-cient approaches to deal with soil moisture heterogeneity at large scales, where manage-ment decisions are usually made. The aim of this dissertation was to test innovative ap-proaches for making efficient use of standard soil hydrological data in order to assess seep-age rates and main controls on observed hydrological behavior, including the role of soil het-erogeneities.
As a first step, the applicability of a simplified Buckingham-Darcy method to estimate deep seepage fluxes from point information of soil moisture dynamics was assessed. This was done in a numerical experiment considering a broad range of soil textures and textural het-erogeneities. The method performed well for most soil texture classes. However, in pure sand where seepage fluxes were dominated by heterogeneous flow fields it turned out to be not applicable, because it simply neglects the effect of water flow heterogeneity. In this study a need for new efficient approaches to handle heterogeneities in one-dimensional water flux models was identified.
As a further step, an approach to turn the problem of soil moisture heterogeneity into a solu-tion was presented: Principal component analysis was applied to make use of the variability among soil moisture time series for analyzing apparently complex soil hydrological systems. It can be used for identifying the main controls on the hydrological behavior, quantifying their relevance, and describing their particular effects by functional averaged time series. The ap-proach was firstly tested with soil moisture time series simulated for different texture classes in homogeneous and heterogeneous model domains. Afterwards, it was applied to 57 mois-ture time series measured in a multifactorial long term field experiment in Northeast Germa-ny.
The dimensionality of both data sets was rather low, because more than 85 % of the total moisture variance could already be explained by the hydrological input signal and by signal transformation with soil depth. The perspective of signal transformation, i.e. analyzing how hydrological input signals (e.g., rainfall, snow melt) propagate through the vadose zone, turned out to be a valuable supplement to the common mass flux considerations. Neither different textures nor spatial heterogeneities affected the general kind of signal transfor-mation showing that complex spatial structures do not necessarily evoke a complex hydro-logical behavior. In case of the field measured data another 3.6% of the total variance was unambiguously explained by different cropping systems. Additionally, it was shown that dif-ferent soil tillage practices did not affect the soil moisture dynamics at all.
The presented approach does not require a priori assumptions about the nature of physical processes, and it is not restricted to specific scales. Thus, it opens various possibilities to in-corporate the key information from monitoring data sets into the modeling exercise and thereby reduce model uncertainties.
TURBO2 - a MATLAB simulation to study the effects of bioturbation on paleoceanographic time series
(2013)
Bioturbation (or benthic mixing) causes significant distortions in marine stable isotope signals and other palaeoceanographic records. Although the influence of bioturbation on these records is well known it has rarely been dealt systematically. The MATLAB program called TURBO2 can be used to simulate the effect of bioturbation on individual sediment particles. It can therefore be used to model the distortion of all physical, chemical, and biological signals in deep-sea sediments, such as Mg/Ca ratios and UK37-based sea-surface temperature (SST) variations. In particular, it can be used to study the distortions in paleoceanographic records that are based on individual sediment particles, such as SST records based on foraminifera assemblages. Furthermore. TURBO2 provides a tool to study the effect of benthic mixing of isotope signals such as C-14, delta O-18, and delta C-13, measured in a stratigraphic carrier such as foraminifera shells.
Much progress has been made in estimating recurrence intervals of great and giant subduction earthquakes using terrestrial, lacustrine, and marine paleoseismic archives. Recent detailed records suggest these earthquakes may have variable recurrence periods and magnitudes forming supercycles. Understanding seismic supercycles requires long paleoseismic archives that record timing and magnitude of such events. Turbidite paleoseismic archives may potentially extend past earthquake records to the Pleistocene and can thus complement commonly shorter-term terrestrial archives. However, in order to unambiguously establish recurring seismicity as a trigger mechanism for turbidity currents, synchronous deposition of turbidites in widely spaced, isolated depocenters has to be ascertained. Furthermore, characteristics that predispose a seismically active continental margin to turbidite paleoseismology and the correct sample site selection have to be taken into account.
Here we analyze 8 marine sediment cores along 950 km of the Chile margin to test for the feasibility of compiling detailed and continuous paleoseismic records based on turbidites. Our results suggest that the deposition of areally widespread, synchronous turbidites triggered by seismicity is largely controlled by sediment supply and, hence, the climatic and geomorphic conditions of the adjacent subaerial setting. The feasibility of compiling a turbidite paleoseismic record depends on the delicate balance between sufficient sediment supply providing material to fail frequently during seismic shaking and sufficiently low sedimentation rates to allow for coeval accumulation of planktonic foraminifera for high-resolution radiocarbon dating.
We conclude that offshore northern central Chile (29-32.5 degrees S) Holocene turbidite paleoseismology is not feasible, because sediment supply from the semi-arid mainland is low and almost no Holocene turbidity-current deposits are found in the cores. In contrast, in the humid region between 36 and 38 degrees S frequent Holocene turbidite deposition may generally correspond to paleoseismic events. However, high terrigenous sedimentation rates prevent high-resolution radiocarbon dating. The climatic transition region between 32.5 and 36 degrees S appears to be best suited for turbidite paleoseismology. (C) 2015 Elsevier Ltd. All rights reserved.
Tuning systems of traditional Georgian singing determined from a new corpus of field recordings
(2022)
In this study we examine the tonal organization of the 2016 GVM dataset, a newly-created corpus of high-quality multimedia field recordings of traditional Georgian singing with a focus on Svaneti. For this purpose, we developed a new processing pipeline for the computational analysis of non-western polyphonic music which was subsequently applied to the complete 2016 GVM dataset.
To evaluate under what conditions a single tuning system is representative of current Svan performance practice, we examined the stability of the obtained tuning systems from an ensemble-, a song-, and a corpus-related perspective.
Furthermore, we compared the resulting Svan tuning systems with the tuning systems obtained for the Erkomaishvili dataset (Rosenzweig et al., 2020) in the study by Scherbaum et al. (2020). In comparison to a 12-TET (12-tone-equal-temperament) system, the Erkomaishvili and the Svan tuning systems are surprisingly similar.
Both systems show a strong presence of pure fourths (500 cents) and fifths (700 cents), and 'neutral' thirds (peaking around 350 cents) as well as 'neutral' sixths.
In addition, the sizes of the melodic and the harmonic seconds in both tuning systems differ systematically from each other, with the size of the harmonic second being systematically larger than the melodic one.
Arctic tundra landscapes are composed of a complex mosaic of patterned ground features, varying in soil moisture, vegetation composition, and surface hydrology over small spatial scales (10-100 m). The importance of microtopography and associated geomorphic landforms in influencing ecosystem structure and function is well founded, however, spatial data products describing local to regional scale distribution of patterned ground or polygonal tundra geomorphology are largely unavailable. Thus, our understanding of local impacts on regional scale processes (e.g., carbon dynamics) may be limited. We produced two key spatiotemporal datasets spanning the Arctic Coastal Plain of northern Alaska (similar to 60,000 km(2)) to evaluate climate-geomorphological controls on arctic tundra productivity change, using (1) a novel 30m classification of polygonal tundra geomorphology and (2) decadal-trends in surface greenness using the Landsat archive (1999-2014). These datasets can be easily integrated and adapted in an array of local to regional applications such as (1) upscaling plot-level measurements (e.g., carbon/energy fluxes), (2) mapping of soils, vegetation, or permafrost, and/or (3) initializing ecosystem biogeochemistry, hydrology, and/or habitat modeling.
Tundra be dammed
(2018)
Increasing air temperatures are changing the arctic tundra biome. Permafrost is thawing, snow duration is decreasing, shrub vegetation is proliferating, and boreal wildlife is encroaching. Here we present evidence of the recent range expansion of North American beaver (Castor canadensis) into the Arctic, and consider how this ecosystem engineer might reshape the landscape, biodiversity, and ecosystem processes. We developed a remote sensing approach that maps formation and disappearance of ponds associated with beaver activity. Since 1999, 56 new beaver pond complexes were identified, indicating that beavers are colonizing a predominantly tundra region (18,293km(2)) of northwest Alaska. It is unclear how improved tundra stream habitat, population rebound following overtrapping for furs, or other factors are contributing to beaver range expansion. We discuss rates and likely routes of tundra beaver colonization, as well as effects on permafrost, stream ice regimes, and freshwater and riparian habitat. Beaver ponds and associated hydrologic changes are thawing permafrost. Pond formation increases winter water temperatures in the pond and downstream, likely creating new and more varied aquatic habitat, but specific biological implications are unknown. Beavers create dynamic wetlands and are agents of disturbance that may enhance ecosystem responses to warming in the Arctic.
Process life cycle assessment (PLCA) is widely used to quantify environmental flows associated with the manufacturing of products and other processes. As PLCA always depends on defining a system boundary, its application involves truncation errors. Different methods of estimating truncation errors are proposed in the literature; most of these are based on artificially constructed system complete counterfactuals. In this article, we review the literature on truncation errors and their estimates and systematically explore factors that influence truncation error estimates. We classify estimation approaches, together with underlying factors influencing estimation results according to where in the estimation procedure they occur. By contrasting different PLCA truncation/error modeling frameworks using the same underlying input-output (I-O) data set and varying cut-off criteria, we show that modeling choices can significantly influence estimates for PLCA truncation errors. In addition, we find that differences in I-O and process inventory databases, such as missing service sector activities, can significantly affect estimates of PLCA truncation errors. Our results expose the challenges related to explicit statements on the magnitude of PLCA truncation errors. They also indicate that increasing the strictness of cut-off criteria in PLCA has only limited influence on the resulting truncation errors. We conclude that applying an additional I-O life cycle assessment or a path exchange hybrid life cycle assessment to identify where significant contributions are located in upstream layers could significantly reduce PLCA truncation errors.