Refine
Year of publication
- 2017 (298) (remove)
Document Type
- Article (229)
- Doctoral Thesis (28)
- Postprint (20)
- Other (11)
- Review (5)
- Part of a Book (2)
- Bachelor Thesis (1)
- Conference Proceeding (1)
- Master's Thesis (1)
Keywords
- Holocene (9)
- climate change (5)
- Arctic (3)
- Chile (3)
- Climate change (3)
- Erosion (3)
- Palaeoclimate (3)
- Tibetan Plateau (3)
- change detection (3)
- earthquake (3)
Institute
- Institut für Geowissenschaften (298) (remove)
The Makran accretionary prism in SE Iran and SW Pakistan is one of the most extensive subduction accretions on Earth. It is characterized by intense folding, thrust faulting and dislocation of the Cenozoic units that consist of sedimentary, igneous and metamorphic rocks. Rock units forming the northern Makran ophiolites are amalgamated as a melange. Metamorphic rocks, including greenschist, amphibolite and blueschist, resulted from metamorphism of mafic rocks and serpentinites. In spite of the geodynamic significance of blueschist in this area, it has been rarely studied. Peak metamorphic phases of the northern Makran mafic blueschist in the Iranshahr area are glaucophane, phengite, quartz +/- omphacite+epidote. Post peak minerals are chlorite, albite and calcic amphibole. Blueschist facies metasedimentary rocks contain garnet, phengite, albite and epidote in the matrix and as inclusions in glaucophane. The calculated P-T pseudosection for a representative metabasic glaucophane schist yields peak pressure and temperature of 11.5-15kbar at 400-510 degrees C. These rocks experienced retrograde metamorphism from blueschist to greenschist facies (350-450 degrees C and 7-8kbar) during exhumation. A back arc basin was formed due to northward subduction of Neotethys under Eurasia (Lut block). Exhumation of the high-pressure metamorphic rocks in northern Makran occurred contemporarily with subduction. Several reverse faults played an important role in exhumation of the ophiolitic and HP-LT rocks. The presence of serpentinite shows the possible role of a serpentinite diapir for exhumation of the blueschist. A tectonic model is proposed here for metamorphism and exhumation of oceanic crust and accretionary sedimentary rocks of the Makran area. Vast accretion of subducted materials caused southward migration of the shore.
The fusion of hyperspectral imaging (HSI) sensor and airborne lidar scanner (ALS) data provides promising potential for applications in environmental sciences. Standard fusion approaches use reflectance information from the HSI and distance measurements from the ALS to increase data dimen-sionality and geometric accuracy. However, the potential for data fusion based on the respective intensity information of the complementary active and passive sensor systems is high and not yet fully exploited. Here, an approach for the rigorous illumination correction of HSI data, based on the radiometric cross-calibrated return intensity information of ALS data, is presented. The cross calibration utilizes a ray tracing-based fusion of both sensor measurements by intersecting their particular beam shapes. The developed method is capable of compensating for the drawbacks of passive HSI systems, such as cast and cloud shadowing effects, illumination changes over time, across track illumination, and partly anisotropy effects. During processing, spatial and temporal differences in illumination patterns are detected and corrected over the entire HSI wavelength domain. The improvement in the classification accuracy of urban and vegetation surfaces demonstrates the benefit and potential of the proposed HSI illumination correction. The presented approach is the first step toward the rigorous in-flight fusion of passive and active system characteristics, enabling new capabilities for a variety of applications.
Scenario-neutral response surfaces illustrate the sensitivity of a simulated natural system, represented by a specific impact variable, to systematic perturbations of climatic parameters. This type of approach has recently been developed as an alternative to top-down approaches for the assessment of climate change impacts. A major limitation of this approach is the underrepresentation of changes in the temporal structure of the climate input data (i.e., the seasonal and day-to-day variability) since this is not altered by the perturbation. This paper presents a framework that aims to examine this limitation by perturbing both observed and projected climate data time series for a future period, which both serve as input into a hydrological model (the HBV model). The resulting multiple response surfaces are compared at a common domain, the standardized runoff response surface (SRRS). We apply this approach in a case study catchment in Norway to (i) analyze possible changes in mean and extreme runoff and (ii) quantify the influence of changes in the temporal structure represented by 17 different climate input sets using linear mixed-effect models. Results suggest that climate change induced increases in mean and peak flow runoff and only small changes in low flow. They further suggest that the effect of the different temporal structures of the climate input data considerably affects low flows and floods (at least 21% influence), while it is negligible for mean runoff.
This article offers a reconstruction of the vegetation and climate of the south-western Siberian Baraba forest-steppe area during the last ca. 8000 years. The analysis of palynological data from the sediment core of Lake Bolshie Toroki using quantitative methods has made it possible to reconstruct changes of the dominant types of vegetation and mean July air temperatures. Coniferous forests grew in the vicinity of the lake, and mean July air temperatures were similar to present-day ones between 7.9 and 7.0 kyr BP. The warmest and driest climate occurred at 7.0-5.0 kyr BP. At that time, the region had open steppe landscapes; birch groves began to spread. A cooling trend is seen after 5.5 kyr BP, when forest-steppe began to emerge. Steppe communities started to dominate again after 1.5 kyr BP. Mean July air temperatures lower than now are reconstructed for the period of 1.9-1 kyr BP, and then the temperatures became similar to present-day ones. Comparing the archaeological data on the types of economy of the population which inhabited the Baraba forest-steppe with the data on changes in the natural environment revealed a connection between the gradual transition from hunting and fishing to livestock breeding and the development of forest-steppe landscapes with a decrease in the area covered by forests. The development of the forest-steppe as an ecotonic landscape starting around 5 kyr BP might have contributed to the coexistence of several archaeological cultures with different types of economy on the same territory. (C) 2017 Elsevier Ltd. All rights reserved.
The Vicuna Pampa volcanic complex, at the SE edge of the arid Puna Plateau of the Central Andes, records the interplay between volcanic construction and degra-dational processes. The low-sloping Vicuna Pampa volcanic complex, with a 1200-m-deep, southeastward-opening depression, was previously interpreted as a collapse caldera based on morphological considerations. However, characteristic features associated with collapse calderas do not exist, and close inspection instead suggests that the Vicuna Pampa volcanic complex is a strongly eroded, broad, massif-type composite volcano of mainly basaltic to trachyandesitic composition. Construction of the Vicuna Pampa volcanic complex occurred during two distinct cycles separated by the development of the depression. The first and main cycle took place at ca. 12 Ma and was dominated by lava flows and subordinate scoria cones and domes. The second cycle, possibly late Miocene in age, affected the SW portion of the depression with the emplacement of domes. We interpret the central depression as the result of a possible sector collapse and subsequent intense fluvial erosion during middle to late Miocene time, facilitated by faulting, steepened topography, and wetter climate conditions compared to today. We estimate that similar to 65% of the initial edifice of similar to 240 km(3) was degraded. The efficiency of degradation processes for removing mass from the Vicuna Pampa volcanic complex is surprising, considering that today the region is arid, and the stream channels within the complex are predominantly transport limited, forming a series of coalesced, aggraded alluvial fans and eolian infill. Hence, the Vicuna Pampa volcanic complex records the effects of past degradation efficiency that differs substantially from that of today.
The characteristics of an aboveground cosmic-ray neutron sensor (CRNS) are evaluated for monitoring a mountain snowpack in the Austrian Alps from March 2014 to June 2016. Neutron counts were compared to continuous point-scale snow depth (SD) and snow-water-equivalent (SWE) measurements from an automatic weather station with a maximum SWE of 600 mm (April 2014). Several spatially distributed Terrestrial Laser Scanning (TLS)-based SD and SWE maps were additionally used. A strong nonlinear correlation is found for both SD and SWE. The representative footprint of the CRNS is in the range of 230-270 m. In contrast to previous studies suggesting signal saturation at around 100 mm of SWE, no complete signal saturation was observed. These results imply that CRNS could be transferred into an unprecedented method for continuous detection of spatially averaged SD and SWE for alpine snowpacks, though with sensitivity decreasing with increasing SWE. While initially different functions were found for accumulation and melting season conditions, this could be resolved by accounting for a limited measurement depth. This depth limit is in the range of 200 mm of SWE for dense snowpacks with high liquid water contents and associated snow density values around 450 kg m(-3) and above. In contrast to prior studies with shallow snowpacks, interannual transferability of the results is very high regardless of presnowfall soil moisture conditions. This underlines the unexpectedly high potential of CRNS to close the gap between point-scale measurements, hydrological models, and remote sensing of the cryosphere in alpine terrain.
Information about the hydrological behaviour of a river basin prior to setting up, calibrating and validating a distributed hydrological model requires extensive datasets that are hardly available for many parts of the world due to insufficient monitoring networks. In this study, the focus was on prevailing spatio-temporal patterns of remotely sensed evapotranspiration (ET) that enabled conclusions to be drawn about the hydrological behaviour and spatial peculiarities of a river basin at rather high spatial resolution. The prevailing spatio-temporal patterns of ET were identified using a principal component analysis of a time series of 644 images of MODIS ET covering the Wami River basin (Tanzania) between the years 2000 and 2013. The time series of the loadings on the principal components were analysed for seasonality and significant long-term trends. The spatial patterns of principal component scores were tested for significant correlation with elevations and slopes, and for differences between different soil texture and land use classes. The results inferred that the temporal and spatial patterns of ET were related to those of preceding rainfalls. At the end of the dry season, high ET was maintained only in areas of shallow groundwater and in cloud forest nature reserves. A region of clear reduction of ET in the long-term was related to massive land use change. The results also confirmed that most soil texture and land use classes differed significantly. Moreover, ET was exceptionally high in natural forests and loam soil, and very low in bushland and sandy-loam soil. Clearly, this approach has shown great potential of publicly available remote sensing data in providing a sound basis for water resources management as well as for distributed hydrological models in data-scarce river basins at lower latitudes.
This study is focused on the analysis of volcanic deposits that crop out at the middle portion of the Las Planchadas range, northern part of the Famatina System in Argentina. These volcanic rocks are records of an Ordovician effusive basaltic volcanism that took place under subaqueous marine conditions. Along the study area crop out an Arenigian volcanic and volcaniclastic rocks succession with massive and autoclastic lavas, hyaloclastites of basaltic composition and volcaniclastic sandstones and mudstones. Large volumes of the volcanic deposits were strongly affected by fragmentation processes during their subaqueous emplacement and in situ accumulated as basaltic breccias. The same volcanic-volcaniclastic association crops out to the south of the Las Planchadas range, forming a basaltic belt with similar characteristics. The geochemical features of the basalts are compatible with depleted mid-ocean ridge basalt (MORB)-like source for the magma, with contribution of subducted related components such as water rich marine hemipelagic sediments, compatible with a back arc geotectonic setting developed along the northern part of the Famatina System during the Arenigian.
Soil properties show high heterogeneity at different spatial scales and their correct characterization remains a crucial challenge over large areas. The aim of the study is to quantify the impact of different types of uncertainties that arise from the unresolved soil spatial variability on simulated hydrological states and fluxes. Three perturbation methods are presented for the characterization of uncertainties in soil properties. The methods are applied on the soil map of the upper Neckar catchment (Germany), as an example. The uncertainties are propagated through the distributed mesoscale hydrological model (mHM) to assess the impact on the simulated states and fluxes. The model outputs are analysed by aggregating the results at different spatial and temporal scales. These results show that the impact of the different uncertainties introduced in the original soil map is equivalent when the simulated model outputs are analysed at the model grid resolution (i.e. 500 m). However, several differences are identified by aggregating states and fluxes at different spatial scales (by subcatchments of different sizes or coarsening the grid resolution). Streamflow is only sensitive to the perturbation of long spatial structures while distributed states and fluxes (e.g. soil moisture and groundwater recharge) are only sensitive to the local noise introduced to the original soil properties. A clear identification of the temporal and spatial scale for which finer-resolution soil information is (or is not) relevant is unlikely to be universal. However, the comparison of the impacts on the different hydrological components can be used to prioritize the model improvements in specific applications, either by collecting new measurements or by calibration and data assimilation approaches. In conclusion, the study underlines the importance of a correct characterization of uncertainty in soil properties. With that, soil maps with additional information regarding the unresolved soil spatial variability would provide strong support to hydrological modelling applications.
The mechanisms leading to large earthquakes are poorly understood and documented. Here we characterize the long-term precursory phase of the 1 April 2014 M(w)8.1 North Chile megathrust. We show that a group of coastal GPS stations accelerated westward 8months before the main shock, corresponding to a M(w)6.5 slow slip event on the subduction interface, 80% of which was aseismic. Concurrent interface foreshocks underwent a diminution of their radiation at high frequency, as shown by the temporal evolution of Fourier spectra and residuals with respect to ground motions predicted by recent subduction models. Such ground motions change suggests that in response to the slow sliding of the subduction interface, seismic ruptures are progressively becoming smoother and/or slower. The gradual propagation of seismic ruptures beyond seismic asperities into surrounding metastable areas could explain these observations and might be the precursory mechanism eventually leading to the main shock.
Failure to consider the costs of adaptation strategies can be seen by decision makers as a barrier to implementing coastal protection measures. In order to validate adaptation strategies to sea-level rise in the form of coastal protection, a consistent and repeatable assessment of the costs is necessary. This paper significantly extends current knowledge on cost estimates by developing - and implementing using real coastal dike data - probabilistic functions of dike costs. Data from Canada and the Netherlands are analysed and related to published studies from the US, UK, and Vietnam in order to provide a reproducible estimate of typical sea dike costs and their uncertainty. We plot the costs divided by dike length as a function of height and test four different regression models. Our analysis shows that a linear function without intercept is sufficient to model the costs, i.e. fixed costs and higher-order contributions such as that due to the volume of core fill material are less significant. We also characterise the spread around the regression models which represents an uncertainty stemming from factors beyond dike length and height. Drawing an analogy with project cost overruns, we employ log-normal distributions and calculate that the range between 3x and x/3 contains 95% of the data, where x represents the corresponding regression value. We compare our estimates with previously published unit costs for other countries. We note that the unit costs depend not only on the country and land use (urban/non-urban) of the sites where the dikes are being constructed but also on characteristics included in the costs, e.g. property acquisition, utility relocation, and project management. This paper gives decision makers an order of magnitude on the protection costs, which can help to remove potential barriers to develop-ing adaptation strategies. Although the focus of this research is sea dikes, our approach is applicable and transferable to other adaptation measures.
Pollination is a key ecosystem service. Pollinators, however, are in decline and their service is increasingly threatened. The decline is driven by several factors, most of which are related to agricultural management. However, the complexity of the landscape system, consisting of both cropped and non cropped areas, makes it difficult to address or even quantify the role of farming practices in pollinator abundance. Therefore, the aim of this paper is to improve our understanding of the relationships between pollinator abundance and their habitat use. We intend to identify and quantify the driving environmental factors that determine pollinator abundance in agricultural landscapes on the crop and landscape scale. These information helps us to design algorithms that can be used as a basis for predicting pollinator abundance on agricultural fields. To integrate varying environmental conditions data sampling was performed on farms in three different regions in Germany. Pollinators were classified into different groups with three aggregation levels. We observed crop parameters as well as landscape parameters in the areas surrounding fields in addition to temporal aspects. Generalized linear models (GLMs) were then calculated. Our results showed that both crop and landscape parameters affect pollinator abundance on agricultural fields. However, the explanatory power of the included parameters varied strongly among the particular pollinator groups and between aggregation levels. Furthermore, differentiation between species groups improves the explanatory power compared to models that are more aggregated. We also found that the temporal match between main activity periods of the particular pollinator groups and resource supply by the crop species is a key factor when analysing pollinator abundance. In conclusion, we demonstrated that the assessment and support of pollination services should be carried out with regard to individual pollinator groups. When studying pollinator abundance, the crop as well as the landscape scale should be addressed. A range of different habitat requirements and different activity periods of the pollinator groups must be covered to maintain pollination services, and therefore both diverse landscapes and diverse crop rotations are of crucial importance. (C) 2017 Elsevier B.V. All rights reserved.
Rhizosphere processes are highly dynamic in time and space and strongly depend on each other. Key factors influencing pH changes in the rhizosphere are root exudation, respiration, and nutrient supply, which are influenced by soil water content levels. In this study, we measured the real-time distribution of soil water, pH changes, and oxygen distribution in the rhizosphere of young maize plants using a recently developed imaging approach. Neutron radiography was used to capture the root system and soil water distribution, while fluorescence imaging was employed to map soil pH and soil oxygen changes. Germinated seeds of maize (Zea mays L.) were planted in glass rhizotrons equipped with pH and oxygen-sensitive sensor foils. After 20 d, the rhizotrons were wetted from the bottom and time-lapsed images via fluorescence and neutron imaging were taken during the subsequent day and night cycles for 5 d. We found higher water content and stronger acidification in the first 0.5 mm from the root surface compared to the bulk soil, which could be a consequence of root exudation. While lateral roots only slightly acidified their rhizosphere, crown roots induced stronger acidification of up to 1 pH unit. We observed changing oxygen patterns at different soil moisture conditions and increasing towards lateral as well as crown roots while extending laterally with ongoing water logging. Our work indicates that plants alter the rhizosphere pH and oxygen also depending on root type, which may indirectly arise also from differences in age and water content changes. The results presented here were possible only by combining different imaging techniques to examine profiles at the root-soil interface in a comprehensive way during wetting and drying.
Streams commonly respond to base-level fall by localizing erosion within steepened, convex knickzone reaches. Localized incision causes knickzone reaches to migrate upstream. Such migrating knickzones dictate the pace of landscape response to changes in tectonics or erosional efficiency and can help quantify the timing and source of base-level fall. Identification of knickzones typically requires individual selection of steepened reaches: a process that is tedious and subjective and has no efficient means to measure knickzone size. We construct an algorithm to automate this procedure by selecting the bounds of knickzone reaches in a -space (drainage-area normalized) framework. An automated feature calibrates algorithm parameters to a subset of knickzones handpicked by the user. The algorithm uses these parameters as consistent criteria to identify knickzones objectively, and then the algorithm measures the height, length, and slope of each knickzone reach. We test the algorithm on 1, 10, and 30m resolution digital elevation models (DEMs) of six catchments (trunk-stream lengths: 2.1-5.4km) on Santa Cruz Island, southern California. On the 1m DEM, algorithm-selected knickzones confirm 93% of handpicked knickzone positions (n=178) to a spatial accuracy of 100m, 88% to an accuracy within 50m, and 46% to an accuracy within 10m. Using 10 and 30m DEMs, accuracy is similar: 88-86% to 100m and 82% to 50m (n=38 and 36, respectively). The algorithm enables efficient regional comparison of the size and location of knickzones with geologic structures, mapped landforms, and hillslope morphology, thereby facilitating approaches to characterize the dynamics of transient landscapes. Plain Language Summary The shape of rivers reflects the environments that they flow through and the environments that they link together: mountains and oceans. Anywhere along the length of a river, changes in environmental conditions are propagated upstream and downstream as the river changes its morphology to match the new environmental conditions. Commonly, rivers steepen as land uplifts faster in regions of high tectonic convergence. The steepening of river gradients is propagated upstream and can be mapped to trace zones of high tectonic activity across landscapes and estimate the source and timing of environmental change. Such insights may indicate regions where earthquakes have become more frequent in the recent past and how rivers respond to these changes. In this submission, we detail an algorithm that can use digital topographic data (similar to google earth), to automatically map and measure anomalously steep river reaches across continental scales. This technology can highlight areas that have experienced recent sustained changes in environmental conditions, evident by changes in the morphology of rivers. Such environmental conditions could be changes in tectonic uplift and earthquake activity, changes in sea level, changes in land-use, or changes in climate, all factors that can produce measurable differences in river morphology over time.
Most of the deformation associated with the seismic cycle in subduction zones occurs offshore and has been therefore difficult to quantify with direct observations at millennial timescales. Here we study millennial deformation associated with an active splay-fault system in the Arauco Bay area off south central Chile. We describe hitherto unrecognized drowned shorelines using high-resolution multibeam bathymetry, geomorphic, sedimentologic, and paleontologic observations and quantify uplift rates using a Landscape Evolution Model. Along a margin-normal profile, uplift rates are 1.3m/ka near the edge of the continental shelf, 1.5m/ka at the emerged Santa Maria Island, -0.1m/ka at the center of the Arauco Bay, and 0.3m/ka in the mainland. The bathymetry images a complex pattern of folds and faults representing the surface expression of the crustal-scale Santa Maria splay-fault system. We modeled surface deformation using two different structural scenarios: deep-reaching normal faults and deep-reaching reverse faults with shallow extensional structures. Our preferred model comprises a blind reverse fault extending from 3km depth down to the plate interface at 16km that slips at a rate between 3.0 and 3.7m/ka. If all the splay-fault slip occurs during every great megathrust earthquake, with a recurrence of similar to 150-200years, the fault would slip similar to 0.5m per event, equivalent to a magnitude similar to 6.4 earthquake. However, if the splay-fault slips only with a megathrust earthquake every similar to 1000years, the fault would slip similar to 3.7m per event, equivalent to a magnitude similar to 7.5 earthquake.
Central Anatolia is a low-relief, high-elevation region where decadal-scale deformation rates estimated from space geodesy suggest low strain rates within a stiff microplate. However, numerous Quaternary faults have been mapped within this low-strain region and estimating their slip rate and seismic potential is important for hazard assessments in an area of increasing infrastructural development. Here we focus on the Sultanhani Fault (SF), which constitutes an integral part of the Eskisehir-Cihanbeyli Fault System, and use deformed maximum highstand shorelines of palaeo-lake Konya to estimate tectonic slip rates at millennial scale. Some of these shorelines were previously interpreted as fault scarps, but we provide conclusive evidence for their erosional origin. We found that shoreline-angle elevations estimated from differential GPS profiles record vertical displacements of 10.2 m across the SF. New radiocarbon ages of lacustrine molluscs suggest 22.4 m of relative lake-level fall between 22.1 +/- 0.3 and 21.7 +/- 0.4 cal. kaBP, constraining the timing of abrupt abandonment of the highstand shoreline. Models of lithospheric rebound associated with regressions of the Tuz Golu and Konya palaeolakes predict only similar to 1 m of regional-scale uplift across the Konya Basin. Dislocation models of displaced shorelines suggest fault-slip rates of 1.5 and 1.8 mm yr(-1) for planar and listric fault geometries, respectively, providing reasonable results for the latter. We found fault scarps in the Nasuhpinar mudflat that likely represent the most recent ground-breaking rupture of the SF, with an average vertical displacement of 1.2 +/- 0.5 m estimated from 54 topographic profiles, equivalent to a M similar to 6.5-6.9 earthquake based on empirical scaling laws. If such events were characteristic during the ultimate 21 ka, a relatively short recurrence time of similar to 800-900 yr would be needed to account for the millennial slip rate. Alternatively, the fault scarp at Nasuhpinar might represent a larger earthquake requiring more frequent smaller events to account for the millennial rate. The relatively fast slip rate of the SF over the past 21 ka is unlikely to have persisted over longer timescales and might reflect spatiotemporal variations in deformation rates within kinematically-linked fault systems within Central Anatolia, or a transient perturbation to the local stress field or fault strength. Such perturbation might have been related to climatically controlled changes in surface and near-surface loads and by interactions among the different tectonic processes that have been proposed to drive the overall slow uplift and associated extension in the Central Anatolian Plateau.
Tropical peatlands now emit hundreds of megatons of carbon dioxide per year because of human disruption of the feedbacks that link peat accumulation and groundwater hydrology. However, no quantitative theory has existed for how patterns of carbon storage and release accompanying growth and subsidence of tropical peatlands are affected by climate and disturbance. Using comprehensive data from a pristine peatland in Brunei Darussalam, we show how rainfall and groundwater flow determine a shape parameter (the Laplacian of the peat surface elevation) that specifies, under a given rainfall regime, the ultimate, stable morphology, and hence carbon storage, of a tropical peatland within a network of rivers or canals. We find that peatlands reach their ultimate shape first at the edges of peat domes where they are bounded by rivers, so that the rate of carbon uptake accompanying their growth is proportional to the area of the still-growing dome interior. We use this model to study how tropical peatland carbon storage and fluxes are controlled by changes in climate, sea level, and drainage networks. We find that fluctuations in net precipitation on timescales from hours to years can reduce long-term peat accumulation. Our mathematical and numerical models can be used to predict long-term effects of changes in temporal rainfall patterns and drainage networks on tropical peatland geomorphology and carbon storage.
Magnetic and electric field observations from the European Space Agency Swarm mission are used to report the direction of electromagnetic energy flux associated with equatorial plasma depletions. Contrary to expectations, the observations suggest a general interhemispheric Poynting flux rather than concurrent flows at both hemispheres toward or away from the equator. Of high interest is a particular behavior noticed over the region with the largest variation in the magnetic declination. This is a Poynting flux flowing mainly into the southern magnetic hemisphere about between 60 degrees W and 30 degrees E and into the northern magnetic hemisphere between 110 degrees W and 60 degrees W. The abrupt change in the flow direction at 60 degrees W is suggested to be caused by an asymmetry between the hemispheres on the ionospheric conductivity, likely due to the influence of thermospheric winds and the presence of the South Atlantic Anomaly.
Lakes are a ubiquitous landscape feature in northern permafrost regions. They have a strong impact on carbon, energy and water fluxes and can be quite responsive to climate change. The monitoring of lake change in northern high latitudes, at a sufficiently accurate spatial and temporal resolution, is crucial for understanding the underlying processes driving lake change. To date, lake change studies in permafrost regions were based on a variety of different sources, image acquisition periods and single snapshots, and localized analysis, which hinders the comparison of different regions. Here, we present a methodology based on machine-learning based classification of robust trends of multi-spectral indices of Landsat data (TM, ETM+, OLI) and object-based lake detection, to analyze and compare the individual, local and regional lake dynamics of four different study sites (Alaska North Slope, Western Alaska, Central Yakutia, Kolyma Lowland) in the northern permafrost zone from 1999 to 2014. Regional patterns of lake area change on the Alaska North Slope (-0.69%), Western Alaska (-2.82%), and Kolyma Lowland (-0.51%) largely include increases due to thermokarst lake expansion, but more dominant lake area losses due to catastrophic lake drainage events. In contrast, Central Yakutia showed a remarkable increase in lake area of 48.48%, likely resulting from warmer and wetter climate conditions over the latter half of the study period. Within all study regions, variability in lake dynamics was associated with differences in permafrost characteristics, landscape position (i.e., upland vs. lowland), and surface geology. With the global availability of Landsat data and a consistent methodology for processing the input data derived from robust trends of multi-spectral indices, we demonstrate a transferability, scalability and consistency of lake change analysis within the northern permafrost region.
Reliable flood risk analyses, including the estimation of damage, are an important prerequisite for efficient risk management. However, not much is known about flood damage processes affecting companies. Thus, we conduct a flood damage assessment of companies in Germany with regard to two aspects. First, we identify relevant damage-influencing variables. Second, we assess the prediction performance of the developed damage models with respect to the gain by using an increasing amount of training data and a sector-specific evaluation of the data. Random forests are trained with data from two postevent surveys after flood events occurring in the years 2002 and 2013. For a sector-specific consideration, the data set is split into four subsets corresponding to the manufacturing, commercial, financial, and service sectors. Further, separate models are derived for three different company assets: buildings, equipment, and goods and stock. Calculated variable importance values reveal different variable sets relevant for the damage estimation, indicating significant differences in the damage process for various company sectors and assets. With an increasing number of data used to build the models, prediction errors decrease. Yet the effect is rather small and seems to saturate for a data set size of several hundred observations. In contrast, the prediction improvement achieved by a sector-specific consideration is more distinct, especially for damage to equipment and goods and stock. Consequently, sector-specific data acquisition and a consideration of sector-specific company characteristics in future flood damage assessments is expected to improve the model performance more than a mere increase in data.
Past climatic change can be reconstructed from sedimentary archives by a number of proxies. However, few methods exist to directly estimate hydrological changes and even fewer result in quantitative data, impeding our understanding of the timing, magnitude and mechanisms of hydrological changes. Here we present a novel approach based on delta H-2 values of sedimentary lipid biomarkers in combination with plant physiological modeling to extract quantitative information on past changes in relative humidity. Our initial application to an annually laminated lacustrine sediment sequence from western Europe deposited during the Younger Dryas cold period revealed relative humidity changes of up to 15% over sub-centennial timescales, leading to major ecosystem changes, in agreement with palynological data from the region. We show that by combining organic geochemical methods and mechanistic plant physiological models on well characterized lacustrine archives it is possible to extract quantitative ecohydrological parameters from sedimentary lipid biomarker delta H-2 data.
Urban climate is determined by a variety of factors, whose knowledge can help to attenuate heat stress in the context of ongoing urbanization and climate change. We study the influence of city size and urban form on the Urban Heat Island (UHI) phenomenon in Europe and find a complex interplay between UHI intensity and city size, fractality, and anisometry. Due to correlations among these urban factors, interactions in the multi-linear regression need to be taken into account. We find that among the largest 5,000 cities, the UHI intensity increases with the logarithm of the city size and with the fractal dimension, but decreases with the logarithm of the anisometry. Typically, the size has the strongest influence, followed by the compactness, and the smallest is the influence of the degree to which the cities stretch. Accordingly, from the point of view of UHI alleviation, small, disperse, and stretched cities are preferable. However, such recommendations need to be balanced against e.g. positive agglomeration effects of large cities. Therefore, trade-offs must be made regarding local and global aims.
Prediction of the area affected by earthquake-induced landsliding based on seismological parameters
(2017)
We present an analytical, seismologically consistent expression for the surface area of the region within which most landslides triggered by an earthquake are located (landslide distribution area). This expression is based on scaling laws relating seismic moment, source depth, and focal mechanism with ground shaking and fault rupture length and assumes a globally constant threshold of acceleration for onset of systematic mass wasting. The seismological assumptions are identical to those recently used to propose a seismologically consistent expression for the total volume and area of landslides triggered by an earthquake. To test the accuracy of the model we gathered geophysical information and estimates of the landslide distribution area for 83 earthquakes. To reduce uncertainties and inconsistencies in the estimation of the landslide distribution area, we propose an objective definition based on the shortest distance from the seismic wave emission line containing 95% of the total landslide area. Without any empirical calibration the model explains 56% of the variance in our dataset, and predicts 35 to 49 out of 83 cases within a factor of 2, depending on how we account for uncertainties on the seismic source depth. For most cases with comprehensive landslide inventories we show that our prediction compares well with the smallest region around the fault containing 95% of the total landslide area. Aspects ignored by the model that could explain the residuals include local variations of the threshold of acceleration and processes modulating the surface ground shaking, such as the distribution of seismic energy release on the fault plane, the dynamic stress drop, and rupture directivity. Nevertheless, its simplicity and first-order accuracy suggest that the model can yield plausible and useful estimates of the landslide distribution area in near-real time, with earthquake parameters issued by standard detection routines.
In 2009, a group of prominent Earth scientists introduced the "planetary boundaries" (PB) framework: they suggested nine global control variables, and defined corresponding "thresholds which, if crossed, could generate unacceptable environmental change". The concept builds on systems theory, and views Earth as a complex adaptive system in which anthropogenic disturbances may trigger nonlinear, abrupt, and irreversible changes at the global scale, and "push the Earth system outside the stable environmental state of the Holocene". While the idea has been remarkably successful in both science and policy circles, it has also raised fundamental concerns, as the majority of suggested processes and their corresponding planetary boundaries do not operate at the global scale, and thus apparently lack the potential to trigger abrupt planetary changes. This paper picks up the debate with specific regard to the planetary boundary on "global freshwater use". While the bio-physical impacts of excessive water consumption are typically confined to the river basin scale, the PB proponents argue that water-induced environmental disasters could build up to planetary-scale feedbacks and system failures. So far, however, no evidence has been presented to corroborate that hypothesis. Furthermore, no coherent approach has been presented to what extent a planetary threshold value could reflect the risk of regional environmental disaster. To be sure, the PB framework was revised in 2015, extending the planetary freshwater boundary with a set of basin-level boundaries inferred from environmental water flow assumptions. Yet, no new evidence was presented, either with respect to the ability of those basin-level boundaries to reflect the risk of regional regime shifts or with respect to a potential mechanism linking river basins to the planetary scale. So while the idea of a planetary boundary on freshwater use appears intriguing, the line of arguments presented so far remains speculative and implicatory. As long as Earth system science does not present compelling evidence, the exercise of assigning actual numbers to such a boundary is arbitrary, premature, and misleading. Taken as a basis for waterrelated policy and management decisions, though, the idea transforms from misleading to dangerous, as it implies that we can globally offset water-related environmental impacts. A planetary boundary on freshwater use should thus be disapproved and actively refuted by the hydrological and water resources community.
The application of heat as a hydrological tracer has become a standard method for quantifying water fluxes between groundwater and surface water. The typical application is to estimate vertical water fluxes in the shallow subsurface beneath streams or lakes. For this purpose, time series of temperatures in the surface water and in the sediment are measured and evaluated by a vertical 1D representation of heat transport by advection and conduction. Several analytical solutions exist to calculate the vertical water flux from the measured temperatures. Although analytical solutions can be easily implemented, they are restricted to specific boundary conditions such as a sinusoidal upper temperature boundary. Numerical solutions offer higher flexibility in the selection of the boundary conditions. This, in turn, reduces the effort of data preprocessing, such as the extraction of the diurnal temperature variation from the raw data. Here, we present software to estimate water fluxes based on temperaturesFLUX-BOT. FLUX-BOT is a numerical code written in MATLAB that calculates vertical water fluxes in saturated sediments based on the inversion of measured temperature time series observed at multiple depths. FLUX-BOT applies a centred Crank-Nicolson implicit finite difference scheme to solve the one-dimensional heat advection-conduction equation. FLUX-BOT includes functions for the inverse numerical routines, functions for visualizing the results, and a function for performing uncertainty analysis. We present applications of FLUX-BOT to synthetic and to real temperature data to demonstrate its performance.
Rotations of the principal stress axes are observed as a result of fluid injection into reservoirs. We use a generic, fully coupled 3-D thermo-hydro-mechanical model to investigate systematically the dependence of this stress rotation on different reservoir properties and injection scenarios. We find that permeability, injection rate, and initial differential stress are the key factors, while other reservoir properties only play a negligible role. In particular, we find that thermal effects do not significantly contribute to stress rotations. For reservoir types with usual differential stress and reservoir treatment the occurrence of significant stress rotations is limited to reservoirs with a permeability of less than approximately 10(-12)m(2). Higher permeability effectively prevents stress rotations to occur. Thus, according to these general findings, the observed principal stress axes rotation can be used as a proxy of the initial differential stress provided that rock permeability and fluid injection rate are known a priori.
The study deals with the identification and characterization of rapid subsurface flow structures through pedo- and geo-physical measurements and irrigation experiments at the point, plot and hillslope scale. Our investigation of flow-relevant structures and hydrological responses refers to the general interplay of form and function, respectively. To obtain a holistic picture of the subsurface, a large set of different laboratory, exploratory and experimental methods was used at the different scales. For exploration these methods included drilled soil core profiles, in situ measurements of infiltration capacity and saturated hydraulic conductivity, and laboratory analyses of soil water retention and saturated hydraulic conductivity. The irrigation experiments at the plot scale were monitored through a combination of dye tracer, salt tracer, soil moisture dynamics, and 3-D time-lapse ground penetrating radar (GPR) methods. At the hillslope scale the subsurface was explored by a 3-D GPR survey. A natural storm event and an irrigation experiment were monitored by a dense network of soil moisture observations and a cascade of 2-D time-lapse GPR "trenches". We show that the shift between activated and non-activated state of the flow paths is needed to distinguish structures from overall heterogeneity. Pedo-physical analyses of point-scale samples are the basis for sub-scale structure inference. At the plot and hillslope scale 3-D and 2-D time-lapse GPR applications are successfully employed as non-invasive means to image subsurface response patterns and to identify flow-relevant paths. Tracer recovery and soil water responses from irrigation experiments deliver a consistent estimate of response velocities. The combined observation of form and function under active conditions provides the means to localize and characterize the structures (this study) and the hydrological processes (companion study Angermann et al., 2017, this issue).
The phrase form and function was established in architecture and biology and refers to the idea that form and functionality are closely correlated, influence each other, and co-evolve. We suggest transferring this idea to hydrological systems to separate and analyze their two main characteristics: their form, which is equivalent to the spatial structure and static properties, and their function, equivalent to internal responses and hydrological behavior. While this approach is not particularly new to hydrological field research, we want to employ this concept to explicitly pursue the question of what information is most advantageous to understand a hydrological system. We applied this concept to subsurface flow within a hillslope, with a methodological focus on function: we conducted observations during a natural storm event and followed this with a hillslope-scale irrigation experiment. The results are used to infer hydrological processes of the monitored system. Based on these findings, the explanatory power and conclusiveness of the data are discussed. The measurements included basic hydrological monitoring methods, like piezometers, soil moisture, and discharge measurements. These were accompanied by isotope sampling and a novel application of 2-D time-lapse GPR (ground-penetrating radar). The main finding regarding the processes in the hillslope was that preferential flow paths were established quickly, despite unsaturated conditions. These flow paths also caused a detectable signal in the catchment response following a natural rainfall event, showing that these processes are relevant also at the catchment scale. Thus, we conclude that response observations (dynamics and patterns, i.e., indicators of function) were well suited to describing processes at the observational scale. Especially the use of 2-D time-lapse GPR measurements, providing detailed subsurface response patterns, as well as the combination of stream-centered and hillslope-centered approaches, allowed us to link processes and put them in a larger context. Transfer to other scales beyond observational scale and generalizations, however, rely on the knowledge of structures (form) and remain speculative. The complementary approach with a methodological focus on form (i.e., structure exploration) is presented and discussed in the companion paper by Jackisch et al. (2017).
Water infiltration in soil is not only affected by the inherent heterogeneities of soil, but even more by the interaction with plant roots and their water uptake. Neutron tomography is a unique non-invasive 3D tool to visualize plant root systems together with the soil water distribution in situ. So far, acquisition times in the range of hours have been the major limitation for imaging 3D water dynamics. Implementing an alternative acquisition procedure we boosted the speed of acquisition capturing an entire tomogram within 10 s. This allows, for the first time, tracking of a water front ascending in a rooted soil column upon infiltration of deuterated water time-resolved in 3D. Image quality and resolution could be sustained to a level allowing for capturing the root system in high detail. Good signal-to-noise ratio and contrast were the key to visualize dynamic changes in water content and to localize the root uptake. We demonstrated the ability of ultra-fast tomography to quantitatively image quick changes of water content in the rhizosphere and outlined the value of such imaging data for 3D water uptake modelling. The presented method paves the way for time-resolved studies of various 3D flow and transport phenomena in porous systems.
The paleoclimate during the Early Eocene in Maritime Antarctica is characterized by cool conditions without a pronounced dry season. Soils formed on volcanic material under such climate conditions in modern analogue environments are usually Andosols rich in nanocrystalline minerals without pedogenic smectite. The paleosols formed on volcanic material on King Georges Island are covered by basalts, dated by 6 new 40Ar/39Ar datings to 51-48 Ma, and are rich in smectite. A pedogenic origin of the smectites would suggest a semi-arid rather than a wet non-seasonal humid paleoclimate. To investigate the origin of the smectites in these paleosols we used X-ray diffraction and microscopic techniques. Minor mineralogical changes between the volcanic parent material and the paleosols and a homogenous distribution of smectites throughout the paleosol horizons indicate that these smectites were mainly inherited from the pyroclastic parent material, which was altered prior to surficial weathering. Nevertheless, the mineralogical properties, such as degree of crystallinity and octahedral site occupancy, of these smectites were modified during the ancient soil formation. Our findings highlight that trioctahedral smectites were a product of deuteric alteration of pyroclastic rocks and were progressively transformed to dioctahedral smectites during weathering in a soil environment on King George Island.
Vast areas of the terrestrial Subarctic and Arctic are underlain by permafrost. Landscape evolution is therefore largely controlled by climate-driven periglacial processes. The response of the frozen ground to late Quaternary warm and cold stages is preserved in permafrost sequences, and deducible by multi-proxy palaeoenvironmental approaches. Here, we analyse radiocarbon-dated mid-Wisconsin Interstadial and Holocene lacustrine deposits preserved in the Kit-1 pingo permafrost sequence combined with water and surface sediment samples from nine modern water bodies on Seward Peninsula (NW Alaska) to reconstruct thermokarst dynamics and determine major abiotic factors that controlled the aquatic ecosystem variability. Our methods comprise taxonomical diatom analyses as well as Detrended Correspondence Analysis (DCA) and Redundancy Analysis (RDA). Our results show, that the fossil diatom record reflects thermokarst lake succession since about 42 C-14 kyr BP. Different thermolcarst lake stages during the mid-Wisconsin Interstadial, the late Wisconsin and the early Holocene are mirrored by changes in diatom abundance, diversity, and ecology. We interpret the taxonomical changes in the fossil diatom assemblages in combination with both modern diatom data from surrounding ponds and existing micropalaeontological, sedimentological and mineralogical data from the pingo sequence. A diatom based quantitative reconstruction of lake water pH indicates changing lake environments during mid-Wisconsin to early Holocene stages. Mineralogical analyses indicate presence of tephra fallout and its impact on fossil diatom communities. Our comparison of modern and fossil diatom communities shows the highest floristic similarity of modern polygon ponds to the corresponding initial (shallow water) development stages of thermolcarst lakes. We conclude, that mid-Wisconsin thermokarst processes in the study area could establish during relatively warm interstadial climate conditions accompanied by increased precipitation due to approaching coasts, while still high continentality and hence high seasonal temperature gradients led to warm summers in the central part of Beringia. (C) 2017 Elsevier B.V. All rights reserved.
The Himalaya has a major influence on global and regional climate, in particular on the Asian monsoon system. The foreland basin of the Himalaya contains a record of tectonics and paleoclimate since the Miocene. Previous work on the evolution of vegetation and climate has focused on the central and western Himalaya, where a shift from C3 to C4 vegetation has been observed at similar to 7 Ma and linked to increased seasonality, but the climatic evolution of the eastern part of the orogen is less well understood. In order to track vegetation as a marker of monsoon intensity and seasonality, we analyzed delta C-13 and 8180 values of soil carbonate and associated delta C-13 values of bulk organic carbon from previously dated sedimentary sections exposing the syn-orogenic detrital Dharamsala and Siwalik Groups in the west, and, for the first time, the Siwalik Group in the east of the Himalayan foreland basin. Sedimentary records span from 20 to 1 Myr in the west (Joginder Nagar, Jawalamukhi, and Haripur Kolar sections) and from 13 to 1 Myr in the east (Kameng section), respectively. The presence of soil carbonate in the west and its absence in the east is a first indication of long-term lateral climatic variation, as soil carbonate requires seasonally arid conditions to develop. delta C-13 values in soil carbonate show a shift from around -10 parts per thousand to -2 parts per thousand at similar to 7 Ma in the west, which is confirmed by delta C-13 analyses on bulk organic carbon that show a shift from around -23 parts per thousand to -19 parts per thousand at the same time. Such a shift in isotopic values is likely to be associated with a change from C3 to C4 vegetation. In contrast, delta C-13 values of bulk organic carbon remain at 23 parts per thousand o in the east. Thus, our data show that the current east -west variation in climate was established at similar to 7 Ma. We propose that the regional change towards a more seasonal climate in the west is linked to a decrease of the influence of the Westerlies, delivering less winter precipitation to the western Himalaya, while the east remained annually humid due to its proximity to the monsoonal moisture source. (C) 2017 Elsevier B.V. All rights reserved.
Located along the Izmir-Ankara-Erzincan Suture (IAES), the Maastrichtian - Paleogene Orhaniye Basin has yielded a highly enigmatic-yet poorly dated- Paleogene mammal fauna, the endemic character of which has suggested high faunal provincialism associated with paleogeographic isolation of the Anatolian landmass during the early Cenozoic. Despite its biogeographic significance, the tectono-stratigraphic history of the Orhaniye Basin has been poorly documented; Here, we combine sedimentary, magnetostratigraphic, and geochronological data to infer the chronology and depositional history of the Orhaniye Basin. We then assess how our new data and interpretations for the Orhaniye Basin impact (1) the timing and mechanisms of seaway closure along the IAES and (2) the biogeographic evolution of Anatolia. Our results show that the Orhaniye Basin initially developed as a forearc basin during the Maastrichtian, before shifting to a retroarc foreland basin setting sometime between the early Paleocene and 44 Ma. This chronology supports a two-step scenario for the assemblage of the central Anatolian landmass, with incipient collision during the Paleocene - Early Eocene and final seaway retreat along the IAES during the earliest Late Eocene after the last marine incursion into the foreland basin. Our dating for the Orhaniye mammal fauna (44-43 Ma) indicates the persistence of faunal endemism in northern Anatolia until at least the late Lutetian despite the advanced stage of IAES closure. The tectonic evolution of dispersal corridors linking northern Anatolia with adjacent parts of Eurasia was not directly associated with IAES closure and consecutive uplifts, but rather with the build-up of continental bridges on the margins of Anatolia, in the Alpine and Tibetan-Himalayan orogens.
Sedimentary basins in the interior of orogenic plateaus can provide unique insights into the early history of plateau evolution and related geodynamic processes. The northern sectors of the Iranian Plateau of the Arabia-Eurasia collision zone offer the unique possibility to study middle-late Miocene terrestrial clastic and volcaniclastic sediments that allow assessing the nascent stages of collisional plateau formation. In particular, these sedimentary archives allow investigating several debated and poorly understood issues associated with the long-term evolution of the Iranian Plateau, including the regional spatio-temporal characteristics of sedimentation and deformation and the mechanisms of plateau growth. We document that middle-late Miocene crustal shortening and thickening processes led to the growth of a basement-cored range (Takab Range Complex) in the interior of the plateau. This triggered the development of a foreland-basin (Great Pari Basin) to the east between 16.5 and 10.7Ma. By 10.7Ma, a fast progradation of conglomerates over the foreland strata occurred, most likely during a decrease in flexural subsidence triggered by rock uplift along an intraforeland basement-cored range (Mahneshan Range Complex). This was in turn followed by the final incorporation of the foreland deposits into the orogenic system and ensuing compartmentalization of the formerly contiguous foreland into several intermontane basins. Overall, our data suggest that shortening and thickening processes led to the outward and vertical growth of the northern sectors of the Iranian Plateau starting from the middle Miocene. This implies that mantle-flow processes may have had a limited contribution toward building the Iranian Plateau in NW Iran.
The concept of a Global Monsoon (GM) has been proposed based on modern precipitation observations, but its application over a wide range of temporal scales is still under debate. Here, we present a synthesis of 268 continental paleo-moisture records collected from monsoonal systems in the Eastern Hemisphere, including the East Asian Monsoon (EAsM), the Indian Monsoon (IM), the East African Monsoon (EAfM), and the Australian Monsoon (AuM) covering the last 18,000 years. The overall pattern of late Glacial to Holocene moisture change is consistent with those inferred from ice cores and marine records. With respect to the last 10,000 years (10 ka), i.e. a period that has high spatial coverage, a Fuzzy c-Means clustering analysis of the moisture index records together with "Xie-Beni" index reveals four clusters of our data set. The paleoclimatic meaning of each cluster is interpreted considering the temporal evolution and spatial distribution patterns. The major trend in the tropical AuM, EAfM, and IM regions is a gradual decrease in moisture conditions since the early Holocene. Moisture changes in the EAsM regions show maximum index values between 8 and 6 ka. However, records located in nearby subtropical areas, i.e. in regions not influenced by the intertropical convergence zone, show an opposite trend compared to the tropical monsoon regions (AuM, EAfM and IM), i.e. a gradual increase. Analyses of modern meteorological data reveal the same spatial patterns as in the paleoclimate records such that, in times of overall monsoon strengthening, lower precipitation rates are observed in the nearby subtropical areas. We explain this pattern as the effect of a strong monsoon circulation suppressing air uplift in nearby subtropical areas, and hence hindering precipitation. By analogy to the modern system, this would mean that during the early Holocene strong monsoon period, the intensified ascending airflows within the monsoon domains led to relatively weaker ascending or even descending airflows in the adjacent subtropical regions, resulting in a precipitation deficit compared to the late Holocene. Our conceptual model therefore integrates regionally contrasting moisture changes into the Global Monsoon hypothesis. (C) 2017 Elsevier Ltd. All rights reserved.
Low-relief plateaus separated by deeply incised fjords are hallmarks of glaciated, passive continental margins. Spectacular examples fringe the once ice-covered North Atlantic coasts of Greenland, Norway and Canada, but low-relief plateau landscapes also underlie present-day ice sheets in Antarctica and Greenland. Dissected plateaus have long been viewed as the outcome of selective linear erosion by ice sheets that focus incision in glacial troughs, leaving the intervening landscapes essentially unaffected. According to this hypothesis, the plateaus are remnants of preglacial low-relief topography. However, here we use computational experiments to show that, like fjords, plateaus are emergent properties of long-term ice-sheet erosion. Ice sheets can either increase or decrease subglacial relief depending on the wavelength of the underlying topography, and plateau topography arises dynamically from evolving feedbacks between topography, ice dynamics and erosion over million-year timescales. This new mechanistic explanation for plateau formation opens the possibility of plateaus contributing significantly to accelerated sediment flux at the onset of the late Cenozoic glaciations, before becoming stable later in the Quaternary.
Radiocarbon and optically stimulated luminescence dating of sediments from Lake Karakul, Tajikistan
(2017)
Lake Karakul in the eastern Pamirs is a large and closed-basin lake in a partly glaciated catchment. Two parallel sediment cores were collected from 12 m water depth. The cores were correlated using XRF analysis and dated using radiocarbon and OSL techniques. The age results of the two dating methods are generally in agreement. The correlated composite core of 12.26 m length represents continuous accumulation of sediments in the lake basin since 31 ka. The lake reservoir effect (LRE) remained relatively constant over this period. High sediment accumulation rates (SedARs) were recorded before 23 ka and after 6.5 ka. The relatively close position of the coring location near the eastern shore of the lake implies that high SedARs resulted from low lake levels. Thus, high SedARs and lower lake levels before 23 ka probably reflect cold and dry climate conditions that inhibited the arrival of moist air at high elevation in the eastern Pamirs. Low lake levels after 6.5 ka were probably caused by declining temperatures after the warmer early Holocene, which had caused a reduction in water resources stored as snow, ice and frozen ground in the catchment. Low SedARs during 23-6.5 ka suggest increased lake levels in Lake Karakul. A short-lived increase of SedARs at 15 ka probably corresponds to the rapid melting of glaciers in the Karakul catchment during the Greenland Interstadial le, shortly after glaciers in the catchment had reached their maximum extents. The sediment cores from Lake Karakul represent an important climate archive with robust chronology for the last glacial interglacial cycle from Central Asia. (C) 2017 Elsevier B.V. All rights reserved.
This study demonstrates the potential of using single-pass TanDEM-X (TDX) radar imagery to analyse inter- and intra-annual glacier changes in mountainous terrain. Based on SAR images acquired in February 2012, March 2013 and November 2013 over the Inylchek Glacier, Kyrgyzstan, we discuss in detail the processing steps required to generate three reliable digital elevation models (DEMs) with a spatial resolution of 10 m that can be used for glacial mass balance studies. We describe the interferometric processing steps and the influence of a priori elevation information that is required to model long wavelength topographic effects. We also focus on DEM alignment to allow optimal DEM comparisons and on the effects of radar signal penetration on ice and snow surface elevations. We finally compare glacier elevation changes between the three TDX DEMs and the C-band shuttle radar topography mission (SRTM) DEM from February 2000. We introduce a new approach for glacier elevation change calculations that depends on the elevation and slope of the terrain. We highlight the superior quality of the TDX DEMs compared to the SRTM DEM, describe remaining DEM uncertainties and discuss the limitations that arise due to the side-looking nature of the radar sensor. (C) 2017 International Society for Photogrammetry and Remote Sensing, Inc. (ISPRS). Published by Elsevier B.V. All rights reserved.
This study examines the characteristics of drought in the Volta River Basin (VRB), investigates the influence of drought on the streamflow, and projects the impacts of future climate change on the drought. A combination of observation data and regional climate simulations of past and future climates (1970-2013, 2046-2065, and 2081-2100) were analyzed for the study. The Standardized Precipitation Index (SPI) and Standardized Precipitation and Evapotranspiration (SPEI) were used to characterize drought while the Standardized Runoff Index (SRI) were used to quantify runoff. Results of the study show that the historical pattern of drought is generally consistent with previous studies over the Basin and most part of West Africa. RCA ensemble medians (RMED) give realistic simulations of drought characteristics and area extent over the Basin and the sub-catchments in the past climate. Generally, an increase in drought intensity and spatial extent are projected over VRB for SPEI and SPI, but the magnitude of increase is higher with SPEI than with SPI. Drought frequency (events per decade) may be magnified by a factor of 1.2, (2046-2065) to 1.6 (2081-2100) compared to the present day episodes in the basin. The coupling between streamflow and drought episodes was very strong (P < 0.05) for the 1-16-year band before the 1970 but showed strong correlation all through the time series period for the 4-8 -years band. Runoff was highly sensitive to precipitation in the VRB and a 2-3 month time lag was found between drought indices and streamflow in the Volta River Basin. Results of this study may guide policymakers in planning how to minimize the negative impacts of future climate change that could have consequences on agriculture, water resources and energy supply.
Kijko et al. (2016) present various methods to estimate parameters that are relevant for probabilistic seismic-hazard assessment. One of these parameters, although not the most influential, is the maximum possible earthquake magnitude m(max). I show that the proposed estimation of m(max) is based on an erroneous equation related to a misuse of the estimator in Cooke (1979) and leads to unstable results. So far, reported finite estimations of m(max) arise from data selection, because the estimator in Kijko et al. (2016) diverges with finite probability. This finding is independent of the assumed distribution of earthquake magnitudes. For the specific choice of the doubly truncated Gutenberg-Richter distribution, I illustrate the problems by deriving explicit equations. Finally, I conclude that point estimators are generally not a suitable approach to constrain m(max).
We use recently deployed seismological arrays in Africa to sample a 2D cross section through the mantle down to the core-mantle boundary (CMB). By making use of travel-time residuals of S, ScS, and SKS phases, a new shear-velocity model of the African low-velocity zone (ALVZ) is derived. Our model suggests between 1.2% shear-velocity reduction at the top and 5% at the bottom with respect to 1D reference models. The average reduction over the whole low-velocity zone (LVZ) amounts to 2% in the presented model and is therefore about twice as strong as values found in global tomographic models. The top of the LVZ reaches up to 1200-km depth, and its lateral extent at the CMB is about 35 degrees. We propose the existence of a gap of 300 km, splitting the structure into two blocks. Our results are based on remarkable differences in SK(K) S travel-time residuals over a few degrees distance. The complexity of the structure could provide a key to an improved understanding of the deep-mantle LVZ dynamics and composition by comparison to geodynamic models. The gap in the model might suggest that the 2D cross section is cutting through a 3D indentation in the boundary of the ALVZ but may also be interpreted as a sign of two individual plumes, rather than one large homogeneous upwelling.
Testing socio-cultural valuation methods of ecosystem services to explain land use preferences
(2017)
Socio-cultural valuation still emerges as a methodological field in ecosystem service (ES) research and until now lacks consistent formalisation and balanced application in ES assessments. In this study, we examine the explanatory value of ES values for land use preferences. We use 563 responses to a survey about the Pentland Hills regional park in Scotland. Specifically, we aim to (1) identify clusters of land use preferences by using a novel visualisation tool, (2) test if socio-cultural values of ESs or (3) user characteristics are linked with land use preferences, and (4) determine whether both socio-cultural values of ESs and user characteristics can predict land use preferences. Our results suggest that there are five groups of people with different land use preferences, ranging from forest and nature enthusiasts to traditionalists, multi-functionalists and recreation seekers. Rating and weighting of ESs and user characteristics were associated with different clusters. Neither socio-cultural values nor user characteristics were suitable predictors for land use preferences. While several studies have explored land use preferences by identifying socio-cultural values in the past, our findings imply that in this case study ES values inform about general perceptions but do not replace the assessment of land use preferences. (C) 2017 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license.
The characteristics of a landscape pose essential factors for hydrological processes. Therefore, an adequate representation of the landscape of a catchment in hydrological models is vital. However, many of such models exist differing, amongst others, in spatial concept and discretisation. The latter constitutes an essential pre-processing step, for which many different algorithms along with numerous software implementations exist. In that context, existing solutions are often model specific, commercial, or depend on commercial back-end software, and allow only a limited or no workflow automation at all. Consequently, a new package for the scientific software and scripting environment R, called lumpR, was developed. lumpR employs an algorithm for hillslope-based landscape discretisation directed to large-scale application via a hierarchical multi-scale approach. The package addresses existing limitations as it is free and open source, easily extendible to other hydrological models, and the workflow can be fully automated. Moreover, it is user-friendly as the direct coupling to a GIS allows for immediate visual inspection and manual adjustment. Sufficient control is furthermore retained via parameter specification and the option to include expert knowledge. Conversely, completely automatic operation also allows for extensive analysis of aspects related to landscape discretisation. In a case study, the application of the package is presented. A sensitivity analysis of the most important discretisation parameters demonstrates its efficient workflow automation. Considering multiple streamflow metrics, the employed model proved reasonably robust to the discretisation parameters. However, parameters determining the sizes of subbasins and hillslopes proved to be more important than the others, including the number of representative hillslopes, the number of attributes employed for the lumping algorithm, and the number of sub-discretisations of the representative hillslopes.
We investigated chironomid assemblages of a well-dated sediment core from a small seepage lake situated at the eastern slope of the Central Kamchatka Mountain Chain, Far East Russia. The chironomid fauna of the investigated Sigrid Lake is dominated by littoral taxa that are sensitive to fluctuations of the water level. Two groups of taxa interchangeably dominate the record responding to the changes in the lake environment during the past 2800 years. The first group of littoral phytophilic taxa includes Psectrocladius sordidellus-type, Corynoneura arctica-type and Dicrotendipes nervosus-type. The abundances of the taxa from this group have the strongest influence on the variations of PCA 1, and these taxa mostly correspond to low water levels, moderate temperatures and slightly acidified conditions. The second group of taxa includes Microtendipes pedellus-type, Tanytarsus lugens-type, and Tanytarsus pallidicornistype. The variations in the abundances of these taxa, and especially of M. pedellus-type, are in accordance with PCA 2 and correspond to the higher water level in the lake, more oligotrophic and neutral pH conditions. Water depths (WD) were reconstructed, using a modern chironomid-based temperature and water depth calibration data set (training set) and inference model from East Siberia (Nazarova et al., 2011). Mean July air temperatures (T July) were inferred using a chironomid-based temperature inference model based on a modern calibration data set for the Far East (Nazarova et al., 2015). The application of transfer functions resulted in reconstructed T July fluctuations of approximately 3 degrees C over the last 2800 years. Low temperatures (11.0-12.0 degrees C) were reconstructed for the periods between ca 1700 and 1500 cal yr BP (corresponding to the Kofun cold stage) and between ca 1200 and 150 cal yr BP (partly corresponding to the Little Ice Age [LIA]). Warm periods (modern T July or higher) were reconstructed for the periods between ca 2700 and 1800 cal yr BP, 1500 and 1300 cal yr BP and after 150 cal yr BP. WD reconstruction revealed that the lake level was lower than its present level at the beginning of the record between ca 2600 and 2300 cal yr BP and ca 550 cal yr BP. Between ca 2300 and 700 cal yr BP as well as between 450 and 150 cal yr BP, the lake level was higher than it is today, most probably reflecting more humid conditions. (C) 2016 Elsevier Ltd and INQUA. All rights reserved.
In order to understand present day earthquake kinematics at the Indian plate boundary, we analyse seismic broadband data recorded between 2007 and 2015 by the regional network in the Garhwal-Kumaun region, northwest Himalaya. We first estimate a local 1-D velocity model for the computation of reliable Green's functions, based on 2837 P-wave and 2680 S-wave arrivals from 251 well located earthquakes. The resulting 1-D crustal structure yields a 4-layer velocity model down to the depths of 20 km. A fifth homogeneous layer extends down to 46 km, constraining the Moho using travel-time distance curve method. We then employ a multistep moment tensor (MT) inversion algorithm to infer seismic moment tensors of 11 moderate earthquakes with Mw magnitude in the range 4.0–5.0. The method provides a fast MT inversion for future monitoring of local seismicity, since Green's functions database has been prepared. To further support the moment tensor solutions, we additionally model P phase beams at seismic arrays at teleseismic distances. The MT inversion result reveals the presence of dominant thrust fault kinematics persisting along the Himalayan belt. Shallow low and high angle thrust faulting is the dominating mechanism in the Garhwal-Kumaun Himalaya. The centroid depths for these moderate earthquakes are shallow between 1 and 12 km. The beam modeling result confirm hypocentral depth estimates between 1 and 7 km. The updated seismicity, constrained source mechanism and depth results indicate typical setting of duplexes above the mid crustal ramp where slip is confirmed along out-of-sequence thrusting. The involvement of Tons thrust sheet in out-of-sequence thrusting indicate Tons thrust to be the principal active thrust at shallow depth in the Himalayan region. Our results thus support the critical taper wedge theory, where we infer the microseismicity cluster as a result of intense activity within the Lesser Himalayan Duplex (LHD) system.