Refine
Has Fulltext
- no (254) (remove)
Year of publication
- 2017 (254) (remove)
Document Type
- Article (229)
- Other (11)
- Doctoral Thesis (7)
- Review (5)
- Part of a Book (2)
Is part of the Bibliography
- yes (254)
Keywords
- Holocene (8)
- Arctic (3)
- Chile (3)
- Climate change (3)
- Palaeoclimate (3)
- Tibetan Plateau (3)
- Alkylpyridinium salts (2)
- Bolboschoenus maritimus (2)
- Central Andes (2)
- Chironomids (2)
Institute
- Institut für Geowissenschaften (254) (remove)
Fluvial fill terraces in intermontane basins are valuable geomorphic archives that can record tectonically and/or climatically driven changes of the Earth-surface process system. However, often the preservation of fill terrace sequences is incomplete and/or they may form far away from their source areas, complicating the identification of causal links between forcing mechanisms and landscape response, especially over multi-millennial timescales. The intermontane Toro Basin in the southern Central Andes exhibits at least five generations of fluvial terraces that have been sculpted into several-hundred-meter-thick Quaternary valley-fill conglomerates. New surface-exposure dating using nine cosmogenic Be-10 depth profiles reveals the successive abandonment of these terraces with a 100 kyr cyclicity between 75 +/- 7 and 487 +/- 34 ka. Depositional ages of the conglomerates, determined by four Al-26/Be-10 burial samples and U-Pb zircon ages of three intercalated volcanic ash beds, range from 18 +/- 141 to 936 +/- 170 ka, indicating that there were multiple cut-and-fill episodes. Although the initial onset of aggradation at similar to 1 Ma and the overall net incision since ca. 500 ka can be linked to tectonic processes at the narrow basin outlet, the superimposed 100 kyr cycles of aggradation and incision are best explained by eccentricity-driven climate change. Within these cycles, the onset of river incision can be correlated with global cold periods and enhanced humid phases recorded in paleoclimate archives on the adjacent Bolivian Altiplano, whereas deposition occurred mainly during more arid phases on the Altiplano and global interglacial periods. We suggest that enhanced runoff during global cold phases - due to increased regional precipitation rates, reduced evapotranspiration, or both - resulted in an increased sediment-transport capacity in the Toro Basin, which outweighed any possible increases in upstream sediment supply and thus triggered incision. Compared with two nearby basins that record precessional (21-kyr) and long-eccentricity (400-kyr) forcing within sedimentary and geomorphic archives, the recorded cyclicity scales with the square of the drainage basin length. (C) 2017 Elsevier B.V. All rights reserved.
Direct current systems employing a kinematic surveying strategy allow to analyze the electrical resistivity of the subsurface for large areas (i.e., several hectares). Typical applications are found in precision agriculture, archaeological prospecting and soil sciences. With the typical survey setting, the collected data sets are often characterized by a rather high level of noise and a rather coarse lateral sampling compared to data acquired with fixed electrodes. We therefore present an efficient one-dimensional inversion approach in which we put special attention on modeling the effects of noise. We apply this method to data recorded with a five-offset equatorial dipole-dipole system employing rolling electrodes. By performing several synthetic tests with realistic noise levels, we found that the considered five-configuration soundings allow for a reliable imaging of two-layer cases in the uppermost two meters of the subsurface, where the subsurface can be assumed to follow a horizontally layered geometry within 3 m around the system. By analyzing the corresponding sensitivity functions, we also show that the equatorial dipole-dipole array is relatively well suited for a 1D inversion approach compared to standard in-line electrode arrays. To illustrate this aspect, we show that our method can provide results similar to those obtained with a 2D Wenner imaging procedure for data recorded across a well-constrained 2D target. We finally apply our method to a large five-offset data set acquired in an agricultural study. The final pseudo-3D model of electrical resistivity is in accordance with borehole data available for the surveyed area. Our results demonstrate the applicability and the versatility of the presented inversion approach for large-scale data sets as they are typically collected with such rolling electrode systems. (C) 2017 Elsevier B.V. All rights reserved.
Accurate time series representation of paleoclimatic proxy records is challenging because such records involve dating errors in addition to proxy measurement errors. Rigorous attention is rarely given to age uncertainties in paleoclimatic research, although the latter can severely bias the results of proxy record analysis. Here, we introduce a Bayesian approach to represent layer-counted proxy records - such as ice cores, sediments, corals, or tree rings - as sequences of probability distributions on absolute, error-free time axes. The method accounts for both proxy measurement errors and uncertainties arising from layer-counting-based dating of the records. An application to oxygen isotope ratios from the North Greenland Ice Core Project (NGRIP) record reveals that the counting errors, although seemingly small, lead to substantial uncertainties in the final representation of the oxygen isotope ratios. In particular, for the older parts of the NGRIP record, our results show that the total uncertainty originating from dating errors has been seriously underestimated. Our method is next applied to deriving the overall uncertainties of the Suigetsu radiocarbon comparison curve, which was recently obtained from varved sediment cores at Lake Suigetsu, Japan. This curve provides the only terrestrial radiocarbon comparison for the time interval 12.5-52.8 kyr BP. The uncertainties derived here can be readily employed to obtain complete error estimates for arbitrary radiometrically dated proxy records of this recent part of the last glacial interval.
Protection of natural or semi-natural ecosystems is an important part of societal strategies for maintaining biodiversity, ecosystem services, and achieving overall sustainable development. The assessment of multiple emerging land use trade-offs is complicated by the fact that land use changes occur and have consequences at local, regional, and even global scale. Outcomes also depend on the underlying socio-economic trends. We apply a coupled, multi-scale modelling system to assess an increase in nature protection areas as a key policy option in the European Union (EU). The main goal of the analysis is to understand the interactions between policy-induced land use changes across different scales and sectors under two contrasting future socio-economic pathways. We demonstrate how complementary insights into land system change can be gained by coupling land use models for agriculture, forestry, and urban areas for Europe, in connection with other world regions. The simulated policy case of nature protection shows how the allocation of a certain share of total available land to newly protected areas, with specific management restrictions imposed, may have a range of impacts on different land-based sectors until the year 2040. Agricultural land in Europe is slightly reduced, which is partly compensated for by higher management intensity. As a consequence of higher costs, total calorie supply per capita is reduced within the EU. While wood harvest is projected to decrease, carbon sequestration rates increase in European forests. At the same time, imports of industrial roundwood from other world regions are expected to increase. Some of the aggregate effects of nature protection have very different implications at the local to regional scale in different parts of Europe. Due to nature protection measures, agricultural production is shifted from more productive land in Europe to on average less productive land in other parts of the world. This increases, at the global level, the allocation of land resources for agriculture, leading to a decrease in tropical forest areas, reduced carbon stocks, and higher greenhouse gas emissions outside of Europe. The integrated modelling framework provides a method to assess the land use effects of a single policy option while accounting for the trade-offs between locations, and between regional, European, and global scales.
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.
A comprehensive hydrometeorological dataset is presented spanning the period 1 January 201131 December 2014 to improve the understanding of the hydrological processes leading to flash floods and the relation between rainfall, runoff, erosion and sediment transport in a mesoscale catchment (Auzon, 116 km(2)) of the Mediterranean region. Badlands are present in the Auzon catchment and well connected to high-gradient channels of bedrock rivers which promotes the transfer of suspended solids downstream. The number of observed variables, the various sensors involved (both in situ and remote) and the space-time resolution (similar to km(2), similar to min) of this comprehensive dataset make it a unique contribution to research communities focused on hydrometeorology, surface hydrology and erosion. Given that rainfall is highly variable in space and time in this region, the observation system enables assessment of the hydrological response to rainfall fields. Indeed, (i) rainfall data are provided by rain gauges (both a research network of 21 rain gauges with a 5 min time step and an operational network of 10 rain gauges with a 5 min or 1 h time step), S-band Doppler dual-polarization radars (1 km(2), 5 min resolution), disdrometers (16 sensors working at 30 s or 1 min time step) and Micro Rain Radars (5 sensors, 100m height resolution). Additionally, during the special observation period (SOP-1) of the HyMeX (Hydrological Cycle in the Mediterranean Experiment) project, two X-band radars provided precipitation measurements at very fine spatial and temporal scales (1 ha, 5 min). (ii) Other meteorological data are taken from the operational surface weather observation stations of Meteo-France (including 2m air temperature, atmospheric pressure, 2 m relative humidity, 10m wind speed and direction, global radiation) at the hourly time resolution (six stations in the region of interest). (iii) The monitoring of surface hydrology and suspended sediment is multi-scale and based on nested catchments. Three hydrometric stations estimate water discharge at a 2-10 min time resolution. Two of these stations also measure additional physico-chemical variables (turbidity, temperature, conductivity) and water samples are collected automatically during floods, allowing further geochemical characterization of water and suspended solids. Two experimental plots monitor overland flow and erosion at 1 min time resolution on a hillslope with vineyard. A network of 11 sensors installed in the intermittent hydrographic network continuously measures water level and water temperature in headwater subcatchments (from 0.17 to 116 km(2)) at a time resolution of 2-5 min. A network of soil moisture sensors enables the continuous measurement of soil volumetric water content at 20 min time resolution at 9 sites. Additionally, concomitant observations (soil moisture measurements and stream gauging) were performed during floods between 2012 and 2014. Finally, this dataset is considered appropriate for understanding the rainfall variability in time and space at fine scales, improving areal rainfall estimations and progressing in distributed hydrological and erosion modelling.
A hydrochemical approach to quantify the role of return flow in a surface flow-dominated catchment
(2017)
Stormflow generation in headwater catchments dominated by subsurface flow has been studied extensively, yet catchments dominated by surface flow have received less attention. We addressed this by testing whether stormflow chemistry is controlled by either (a) the event-water signature of overland flow, or (b) the pre-event water signature of return flow. We used a high-resolution hydrochemical data set of stormflow and end-members of multiple storms in an end-member mixing analysis to determine the number of end-members needed to explain stormflow, characterize and identify potential end-members, calculate their contributions to stormflow, and develop a conceptual model of stormflow. The arrangement and relative positioning of end-members in stormflow mixing space suggest that saturation excess overland flow (26-48%) and return flow from two different subsurface storage pools (17-53%) are both similarly important for stormflow. These results suggest that pipes and fractures are important flow paths to rapidly release stored water and highlight the value of within-event resolution hydrochemical data to assess the full range and dynamics of flow paths.
Lakes are dominant and diverse landscape features in the Arctic, but conventional land cover classification schemes typically map them as a single uniform class. Here, we present a detailed lake-centric geospatial database for an Arctic watershed in northern Alaska. We developed a GIS dataset consisting of 4362 lakes that provides information on lake morphometry, hydrologic connectivity, surface area dynamics, surrounding terrestrial ecotypes, and other important conditions describing Arctic lakes. Analyzing the geospatial database relative to fish and bird survey data shows relations to lake depth and hydrologic connectivity, which are being used to guide research and aid in the management of aquatic resources in the National Petroleum Reserve in Alaska. Further development of similar geospatial databases is needed to better understand and plan for the impacts of ongoing climate and land-use changes occurring across lake-rich landscapes in the Arctic.
The features of Middle Miocene deposits in the Puna-Eastern Cordillera transition (Valles Calchaquies) indicate that Cenozoic deformation, sedimentation and volcanism follow a complex spatiotemporal relationship. The intense volcanic activity recorded in the eastern Puna border between 14 and 11.5 Ma coincides with the occurrence of one of the most important deformation events of the Neogene tectonic evolution in the region. Studies performed across the Puna-Eastern Cordillera transition show different relationships between volcanic deposits of ca. 13.5-12.1 Ma and the Oligocene-Miocene Angastaco Formation. In this paper we describe the ash-flow tuff deposits which are the first of this type found concordant in the sedimentary fill of Valles Calchaquies. Several analyses performed on these pyroclastic deposits allow a correlation to be made with the Alto de Las Lagunas Ignimbrite (ca. 13.5 Ma) of the Pucarilla-Cerro Tipillas Volcanic Complex located in the Puna. Outcrops of the ca. 13.5 Ma pyroclastic deposits are recognised within the Puna and the Valle Calchaqui. However, in the southern prolongation of the Valle de Hualfin (Tiopampa-Pucarilla depression) that separates the Puna from the Valle Calchaqui at these latitudes, these deposits are partially eroded and buried, and thus their occurrence is recorded only by abundant volcanic clasts included in conglomerates of the Angastaco Formation. The sedimentation of the Angastaco Formation was aborted at ca. 12 Ma in the Tiopampa-Pucarilla depression by the Pucarilla Ignimbrite, which unconformably covers the synorogenic units. On the contrary, in the Valle Calchaqui the sedimentation of the Angastaco Formation continued until the Late Miocene. The different relationships between the Miocene Angastaco Formation and the ignimbrites with ages of ca. 13.5 and ca. 12 Ma reveal that in this short period (-1.5 m.y.) a significant deformation event took place and resulted in marked palaeogeographic changes, as evidenced by stratigraphic-sedimentological and chronological records in the Angastaco Formation. (C) 2017 Elsevier B.V. All rights reserved.
The polar and subtropical jet streams are strong upper-level winds with a crucial influence on weather throughout the Northern Hemisphere midlatitudes. In particular, the polar jet is located between cold arctic air to the north and warmer subtropical air to the south. Strongly meandering states therefore often lead to extreme surface weather. The parameter values of the detection scheme are optimized using simulated annealing and a skill function that accounts for the zonal-mean jet stream position (Rikus, 2015). After the successful optimization process, we apply our scheme to reanalysis data covering 1979-2015 and calculate seasonal-mean probabilistic maps and trends in wind strength and position of jet streams. We present longitudinally defined probability distributions of the positions for both jets for all on the Northern Hemisphere seasons. This shows that winter is characterized by two well-separated jets over Europe and Asia (ca. 20 degrees W to 140 degrees E). In contrast, summer normally has a single merged jet over the western hemisphere but can have both merged and separated jet states in the eastern hemisphere.
Processes driving the production, transformation and transport of methane (CH4 / in wetland ecosystems are highly complex. We present a simple calculation algorithm to separate open-water CH4 fluxes measured with automatic chambers into diffusion-and ebullition-derived components. This helps to reveal underlying dynamics, to identify potential environmental drivers and, thus, to calculate reliable CH4 emission estimates. The flux separation is based on identification of ebullition-related sudden concentration changes during single measurements. Therefore, a variable ebullition filter is applied, using the lower and upper quartile and the interquartile range (IQR). Automation of data processing is achieved by using an established R script, adjusted for the purpose of CH4 flux calculation. The algorithm was validated by performing a laboratory experiment and tested using flux measurement data (July to September 2013) from a former fen grassland site, which converted into a shallow lake as a result of rewetting. Ebullition and diffusion contributed equally (46 and 55 %) to total CH4 emissions, which is comparable to ratios given in the literature. Moreover, the separation algorithm revealed a concealed shift in the diurnal trend of diffusive fluxes throughout the measurement period. The water temperature gradient was identified as one of the major drivers of diffusive CH4 emissions, whereas no significant driver was found in the case of erratic CH4 ebullition events.
Sustainable development goals (SDGs) have set the 2030 agenda to transform our world by tackling multiple challenges humankind is facing to ensure well-being, economic prosperity, and environmental protection. In contrast to conventional development agendas focusing on a restricted set of dimensions, the SDGs provide a holistic and multidimensional view on development. Hence, interactions among the SDGs may cause diverging results. To analyze the SDG interactions we systematize the identification of synergies and trade-offs using official SDG indicator data for 227 countries. A significant positive correlation between a pair of SDG indicators is classified as a synergy while a significant negative correlation is classified as a trade-off. We rank synergies and trade-offs between SDGs pairs on global and country scales in order to identify the most frequent SDG interactions. For a given SDG, positive correlations between indicator pairs were found to outweigh the negative ones in most countries. Among SDGs the positive and negative correlations between indicator pairs allowed for the identification of particular global patterns. SDG 1 (No poverty) has synergetic relationship with most of the other goals, whereas SDG 12 (Responsible consumption and production) is the goal most commonly associated with trade-offs. The attainment of the SDG agenda will greatly depend on whether the identified synergies among the goals can be leveraged. In addition, the highlighted trade-offs, which constitute obstacles in achieving the SDGs, need to be negotiated and made structurally nonobstructive by deeper changes in the current strategies.
Landscape evolution models (LEMs) allow the study of earth surface responses to changing climatic and tectonic forcings. While much effort has been devoted to the development of LEMs that simulate a wide range of processes, the numerical accuracy of these models has received less attention. Most LEMs use first-order accurate numerical methods that suffer from substantial numerical diffusion. Numerical diffusion particularly affects the solution of the advection equation and thus the simulation of retreating landforms such as cliffs and river knickpoints. This has potential consequences for the integrated response of the simulated landscape. Here we test a higher-order flux-limiting finite volume method that is total variation diminishing (TVD-FVM) to solve the partial differential equations of river incision and tectonic displacement. We show that using the TVD-FVM to simulate river incision significantly influences the evolution of simulated landscapes and the spatial and temporal variability of catchment-wide erosion rates. Furthermore, a two-dimensional TVD-FVM accurately simulates the evolution of landscapes affected by lateral tectonic displacement, a process whose simulation was hitherto largely limited to LEMs with flexible spatial discretization. We implement the scheme in TTLEM (TopoToolbox Landscape Evolution Model), a spatially explicit, raster-based LEM for the study of fluvially eroding landscapes in TopoToolbox 2.
The Lower Rhine Graben (Central Europe) is a prime example of a seismically active low-strain rift zone characterized by pronounced anthropogenic and climatic overprint of structures, and long recurrence intervals of large earthquakes. These factors render the identification of active faults and surface ruptures difficult. We investigated two fault scarps in the Lower Rhine Graben, to decipher their structural character, offset and potential seismogenic origin. Both scarps were modified by anthropogenic activity. The Hemmerich site lies c. 20 km SW of Cologne, along the Erft Fault. The Untermaubach site lies SW of Duren, where the Schafberg Fault projects into the Rur River valley. At the Hemmerich site, geomorphic and geophysical data, as well as exploratory coring reveal evidence of repeated normal faulting. Geophysical analysis and palaeoseismological excavation at the Untermaubach site reveal a complex fault zone in Holocene gravels characterized by subtle gravel deformation. Differentiation of tectonic and fluvial features was only possible with trenching, because fault structures and grain sizes of the sediments were below the resolution of the geophysical data. Despite these issues, our investigation demonstrates that valuable insight into past earthquakes and seismogenic deformation in a low-strain environment can be revealed using a multidisciplinary approach.
Adaptation to flood risk
(2017)
As flood impacts are increasing in large parts of the world, understanding the primary drivers of changes in risk is essential for effective adaptation. To gain more knowledge on the basis of empirical case studies, we analyze eight paired floods, that is, consecutive flood events that occurred in the same region, with the second flood causing significantly lower damage. These success stories of risk reduction were selected across different socioeconomic and hydro-climatic contexts. The potential of societies to adapt is uncovered by describing triggered societal changes, as well as formal measures and spontaneous processes that reduced flood risk. This novel approach has the potential to build the basis for an international data collection and analysis effort to better understand and attribute changes in risk due to hydrological extremes in the framework of the IAHSs Panta Rhei initiative. Across all case studies, we find that lower damage caused by the second event was mainly due to significant reductions in vulnerability, for example, via raised risk awareness, preparedness, and improvements of organizational emergency management. Thus, vulnerability reduction plays an essential role for successful adaptation. Our work shows that there is a high potential to adapt, but there remains the challenge to stimulate measures that reduce vulnerability and risk in periods in which extreme events do not occur.
The early Cretaceous Paraná–Etendeka Large Igneous Province is attributed to the impact of the Tristan mantle plume on the base of the continental lithosphere and the associated opening of the South Atlantic Ocean during the breakup of West Gondwana. Although the geochemistry of the Paraná and Etendeka volcanic rocks has been extensively studied, there is still disagreement on the role of the mantle plume in the production of the magma types observed, because some of their primary compositions are obscured by continental crustal contamination. However, there are related plutonic rocks that preserve mantle signatures. The Doros Complex is a shallow-level mafic intrusion within the Etendeka Province of Namibia. New 39Ar/40Ar step-heating ages for Doros gabbros from this study (weighted mean of 130 ± 1 Ma; 2σ error) confirm contemporaneity with the Paraná–Etendeka magmatic event. The Doros suite yields mean ɛNd values of +5·3 ± 1·0 (1σ; n = 11), initial 87Sr/86Sr = 0·70418 ± 0·00017 (n = 11) and 206Pb/204Pb = 18·11 ± 0·06 (n = 13) at 132 Ma. The clustering of isotopic data and trends in incompatible trace element ratios indicate that all the magmas in the complex were derived from the same mantle source components, during the same melting episode. By quantitative isotopic modelling of mixing processes, we constrain the Doros parental magma to comprise 60–80% melt of a depleted asthenospheric mantle component and 20–40% melt of a more enriched, Tristan plume-derived, asthenospheric component. No lithospheric mantle component is required to explain the Doros magma compositions. The chilled margin to the complex is the only rock type that shows evidence of significant continental crustal contamination, by assimilation of the metasedimentary host-rock upon emplacement. The identification of a substantial Tristan plume component in the Doros source confirms the integral role of the deep-seated thermal anomaly in Paraná–Etendeka magmatism. We show, in addition, that the Doros suite has consistent, strong geochemical affinities with the Tafelkop group ‘ferropicrite’ lavas of the Etendeka Province. This provides crucial evidence in support of Doros as the eruptive site for the Tafelkop lavas, thereby linking the Doros magmatism to the earliest eruptive phase in the Etendeka event. The distinctive chemistry of this magma group has been attributed to relatively deep decompression melting of pyroxenite-bearing material in the heterogeneous Tristan plume head, related to the initial impact of the plume on the base of the lithosphere.
Six N-alkylpyridinium salts [CnPy](2)[MCl4] (n = 4 or 12 and M = Co, Cu, Zn) were synthesized, and their structure and thermal properties were studied. The [C4Py](2)[MCl4] compounds are monoclinic and crystallize in the space group P2(1)/n. The crystals of the longer chain analogues [C12Py](2)[MCl4] are triclinic and crystallize in the space group P (1) over bar. Above the melting temperature, all compounds are ionic liquids (ILs). The derivatives with the longer C12 chain exhibit liquid crystallinity and the shorter chain compounds only show a melting transition. Consistent with single-crystal analysis, electron paramagnetic resonance spectroscopy suggests that the [CuCl4](2-) ions in the Cu-based ILs have a distorted tetrahedral geometry.
Six N-alkylpyridinium salts [CnPy](2)[MCl4] (n = 4 or 12 and M = Co, Cu, Zn) were synthesized, and their structure and thermal properties were studied. The [C4Py](2)[MCl4] compounds are monoclinic and crystallize in the space group P2(1)/n. The crystals of the longer chain analogues [C12Py](2)[MCl4] are triclinic and crystallize in the space group P (1) over bar. Above the melting temperature, all compounds are ionic liquids (ILs). The derivatives with the longer C12 chain exhibit liquid crystallinity and the shorter chain compounds only show a melting transition. Consistent with single-crystal analysis, electron paramagnetic resonance spectroscopy suggests that the [CuCl4](2-) ions in the Cu-based ILs have a distorted tetrahedral geometry.
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.
We introduce the Annual Magnetospheric Currents index as long-term proxy for the geomagnetic signal of magnetospheric currents on Earth valid within the time span 1900–2010. Similar to the widely used disturbance storm time and ‘Ring Current’ indices, it is based on geomagnetic observatory data, but provides a realistic absolute level and uncertainty estimates. Crucial aspects to this end are the revision of observatory crustal biases as well as the implementation of a Bayesian inversion accounting for uncertainties in the main field estimate, both required for the index derivation. The observatory choice is based on a minimization of index variance during a reference period spanning 1960–2010. The new index is capable of correcting observatory time series from large-scale external signals in a user-friendly manner. At present the index is only available as annual mean values. An extension to hourly values for the same time span is in progress.
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.
In Porisdalur valley, a small relict of a sedimentary body was identified in southeastern Iceland. It probably represents a remnant of the deep, tectonically arranged paleolake (Late Miocene, 8-9 Ma), and filled by volcaniclastic material from nearby, active volcanic centers. In the profile of tuffitic sandstone, siltstone and claystone, the ripple-bedding layers, molds and flute casts indicate periodic mass flow episodes. In the sedimentary profile, the characteristic arrangement of sediments is evident, showing features of the Bouma sequences. In the claystone layers, deposited during episodes of lowest kinetic energy, a specific ichnoassemblage was found, represented by Thorichnus ramosus igen. et isp. nov., T. corniculatus igen. et isp. nov., Mammillichnis jakubi isp. nov., Helminthoidichnites multilaqueatus comb. nov., Vamaspor jachymi igen. et isp. nov. and five preliminarily identified trace fossils. The assemblage belongs to Mermia ichnofacies, the nonmarine representative of an ichnofacies, developed in a turbiditic environment; most of identified trace fossils are so far endemic. (C) 2017 Elsevier B.V. All rights reserved.
This study combines spaceborne multitemporal and hyperspectral data to analyze the spatial distribution of surface evaporite minerals and changes in a semi-arid depositional environment associated with episodic flooding events, the Omongwa salt pan (Kalahari, Namibia). The dynamic of the surface crust is evaluated by a change-detection approach using the Iterative-reweighted Multivariate Alteration Detection (IR-MAD) based on the Landsat archive imagery from 1984 to 2015. The results show that the salt pan is a highly dynamic and heterogeneous landform. A change gradient is observed from very stable pan border to a highly dynamic central pan. On the basis of hyperspectral EO-1 Hyperion images, the current distribution of surface evaporite minerals is characterized using Spectral Mixture Analysis (SMA). Assessment of field and image endmembers revealed that the pan surface can be categorized into three major crust types based on diagnostic absorption features and mineralogical ground truth data. The mineralogical crust types are related to different zones of surface change as well as pan morphology that influences brine flow during the pan inundation and desiccation cycles. These combined information are used to spatially map depositional environments where the more dynamic halite crust concentrates in lower areas although stable gypsum and calcite/sepiolite crusts appear in higher elevated areas.
The Misho complex in Northwest Iran is a prominent topographic massif bounded by well known active faults. Our new structural analysis of this area indicates that faulting has important role in the exhumation of this complex. The conjugate orientation of the North and South Misho Faults caused uplift in the Misho and exhumation of the Precambrian crystalline basement. Our structural and stratigraphic data shows that rapid uplift could have been initiation since the 21-22 Ma and exhumation rate was about 0.16 to 0.24 km/Ma. To refine this age, we performed U/Pb analysis of detrital zircon from the Upper Red Formation using LA-ICP-MS. We conducted AFT analysis on 6 basement samples from the hanging wall and 1 sample from the Upper Red Formation in the footwall NMF. Uplift in the hanging wall of NMF led to resting of sample 916 marl. This geochronologic and thermochronologic data shows that exhumation in the MC is diachronously along strike and affected by faults. The phase of exhumation is documented in the study area and entire Iranian plateau is related to the final closure of the Neo-Tethys and northward motion of the Arabian Plate.
The shape of trait distributions may inform about the selective forces that structure ecological communities. Here, we present a new moment-based approach to classify the shape of observed biomass-weighted trait distributions into normal, peaked, skewed, or bimodal that facilitates spatio-temporal and cross-system comparisons. Our observed phytoplankton trait distributions exhibited substantial variance and were mostly skewed or bimodal rather than normal. Additionally, mean, variance, skewness und kurtosis were strongly correlated. This is in conflict with trait-based aggregate models that often assume normally distributed trait values and small variances. Given these discrepancies between our data and general model assumptions we used the observed trait distributions to test how well different aggregate models with first- or second-order approximations and different types of moment closure predict the biomass, mean trait, and trait variance dynamics using weakly or moderately nonlinear fitness functions. For weakly non-linear fitness functions aggregate models with a second-order approximation and a data-based moment closure that relied on the observed correlations between skewness and mean, and kurtosis and variance predicted biomass and often also mean trait changes fairly well and better than models with first-order approximations or a normal-based moment closure. In contrast, none of the models reflected the changes of the trait variances reliably. Aggregate model performance was often also poor for moderately nonlinear fitness functions. This questions a general applicability of the normal-based approach, in particular for predicting variance dynamics determining the speed of trait changes and maintenance of biodiversity. We evaluate in detail how and why better approximations can be obtained.
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.
The information about climate change impact on river discharge is vitally important for planning adaptation measures. The future changes can affect different water-related sectors. The main goal of this study was to investigate the potential water resource changes in Ukraine, focusing on three mesoscale river catchments (Teteriv, UpperWestern Bug, and Samara) characteristic for different geographical zones. The catchment scale watershed model—Soil and Water Integrated Model (SWIM)—was setup, calibrated, and validated for the three catchments under consideration. A set of seven GCM-RCM (General Circulation Model-Regional Climate Model) coupled climate scenarios corresponding to RCPs (Representative Concentration Pathways) 4.5 and 8.5 were used to drive the hydrological catchment model. The climate projections, used in the study, were considered as three combinations of low, intermediate, and high end scenarios. Our results indicate the shifts in the seasonal distribution of runoff in all three catchments. The spring high flow occurs earlier as a result of temperature increases and earlier snowmelt. The fairly robust trend is an increase in river discharge in the winter season, and most of the scenarios show a potential decrease in river discharge in the spring.
Avalanche activity is an important factor when estimating the regional avalanche danger. Moreover, a complete and detailed picture of avalanche activity is needed to understand the processes that lead to natural avalanche release. Currently, information on avalanche activity is mainly obtained through visual observations. However, this involves large uncertainties in the number and release times, influencing the subsequent analysis. Therefore, alternative methods for the remote detection of snow avalanches in particular in non-observed areas are highly desirable. In this study, we use the excited ground vibration to identify avalanches automatically. The specific seismic signature of avalanches facilitates the objective detection by a recently developed classification procedure. A probabilistic description of the signals, called hidden Markov models, allows the robust identification of corresponding signals in the continuous data stream. The procedure is based upon learning a general background model from continuous seismic data. Then, a single reference waveform is used to update an event-specific classifier. Thus, a minimum amount of training data is required by constructing such a classifier on the fly. In this study, we processed five days of continuous data recorded in the Swiss Alps during the avalanche winter 1999. With the restriction of testing large wet-snow avalanches only, the presented approach achieved very convincing results. We successfully detect avalanches over a large volume and distance range. Ninety-two percentage of all detections (43 out of 47) could be confirmed as avalanche events; only four false alarms are reported. We see a clear dependence of recognition capability on run-out distance and source-receiver distance of the observed events: Avalanches are detectable up to a source-receiver distance of eight times the avalanche length. Implications for analyzing a more comprehensive data set (smaller events and different flow regimes) are discussed in detail.
River-bed disturbance and associated sedimentary processes such as particle mobility are central elements to assess river ecosystem functioning. Dams change river dynamics affecting the magnitude and frequency of biophysical elements that depend on them. This paper examines the effects of two dams different in size, management, and location, on the flow regime, flood competence, and bed disturbance in two contrasting Mediterranean rivers, the Esera and the Siurana. For this purpose, two reaches on each river were monitored upstream and downstream from reservoirs. Several monitoring and modeling techniques were used to characterize flow competence, particle entrainment, and the volumes of sediments eroded and deposited after floods. The flow regime of the Esera has been modified from nivo-pluvial regime, typical of humid mountainous environments, to that observed in dry semiarid regions, in which high magnitude but low frequency floods are the dominant processes. Conversely, the flow regime of the Siurana has changed from a typical Mediterranean stream to a regime observed in more temperate environments, with more permanent and stable flows. Both rivers show notably physical changes, with channels clearly less dynamic below the dams. The lack of competent flows together with the sediment deficit associated with the dams has led to less active fluvial environments (reduced particle mobility and bed scour dynamics), a fact that affects instream habitat structure (more uniform grain size distribution, less physical heterogeneity, more stable flows), overall contributing to the degradation of the stream corridor and the subsequent environmental deterioration of the whole fluvial landscape. Copyright (c) 2017 John Wiley & Sons, Ltd.
During the second phase of the Alpine Fault, Deep Fault Drilling Project (DFDP) in the Whataroa River, South Westland, New Zealand, bedrock was encountered in the DFDP-2B borehole from 238.5–893.2 m Measured Depth (MD). Continuous sampling and meso- to microscale characterisation of whole rock cuttings established that, in sequence, the borehole sampled amphibolite facies, Torlesse Composite Terrane-derived schists, protomylonites and mylonites, terminating 200–400 m above an Alpine Fault Principal Slip Zone (PSZ) with a maximum dip of 62°. The most diagnostic structural features of increasing PSZ proximity were the occurrence of shear bands and reduction in mean quartz grain sizes. A change in composition to greater mica:quartz + feldspar, most markedly below c. 700 m MD, is inferred to result from either heterogeneous sampling or a change in lithology related to alteration. Major oxide variations suggest the fault-proximal Alpine Fault alteration zone, as previously defined in DFDP-1 core, was not sampled.
The large variety of atmospheric circulation systems affecting the eastern Asian climate is reflected by the complex Asian vegetation distribution. Particularly in the transition zones of these circulation systems, vegetation is supposed to be very sensitive to climate change. Since proxy records are scarce, hitherto a mechanistic understanding of the past spatio-temporal climate-vegetation relationship is lacking. To assess the Holocene vegetation change and to obtain an ensemble of potential mid-Holocene biome distributions for eastern Asia, we forced the diagnostic biome model BIOME4 with climate anomalies of different transient Holocene climate simulations performed in coupled atmosphere-ocean(-vegetation) models. The simulated biome changes are compared with pollen-based biome records for different key regions. In all simulations, substantial biome shifts during the last 6000 years are confined to the high northern latitudes and the monsoon-westerly wind transition zone, but the temporal evolution and amplitude of change strongly depend on the climate forcing. Large parts of the southern tundra are replaced by taiga during the mid-Holocene due to a warmer growing season and the boreal treeline in northern Asia is shifted northward by approx. 4 degrees in the ensemble mean, ranging from 1.5 to 6 degrees in the individual simulations, respectively. This simulated treeline shift is in agreement with pollen-based reconstructions from northern Siberia. The desert fraction in the transition zone is reduced by 21% during the mid-Holocene compared to pre-industrial due to enhanced precipitation. The desert-steppe margin is shifted westward by 5 degrees (1-9 degrees in the individual simulations). The forest biomes are expanded north-westward by 2 degrees, ranging from 0 to 4 degrees in the single simulations. These results corroborate pollen-based reconstructions indicating an extended forest area in north-central China during the mid-Holocene. According to the model, the forest-to-non-forest and steppe-to-desert changes in the climate transition zones are spatially not uniform and not linear since the mid-Holocene.
Sedimentary lipid biomarkers have become widely used tools for reconstructing past climatic and ecological changes due to their ubiquitous occurrence in lake sediments. In particular, the hydrogen isotopic composition (expressed as delta D values) of leaf wax lipids derived from terrestrial plants has been a focus of research during the last two decades and the understanding of competing environmental and plant physiological factors influencing the delta D values has greatly improved. Comparatively less attention has been paid to lipid biomarkers derived from aquatic plants, although these compounds are abundant in many lacustrine sediments. We therefore conducted a field and laboratory experiment to study the effect of salinity and groundwater discharge on the isotopic composition of aquatic plant biomarkers. We analyzed samples of the common submerged plant species, Potamogeton pectinatus (sago pondweed), which has a wide geographic distribution and can tolerate high salinity. We tested the effect of groundwater discharge (characterized by more negative delta D values relative to lake water) and salinity on the delta D values of n-alkanes from P. pectinatus by comparing plants (i) collected from the oligotrophic freshwater Lake Stechlin (Germany) at shallow littoral depth from locations with and without groundwater discharge, and (ii) plants grown from tubers collected from the eutrophic Lake Muggelsee in nutrient solution at four salinity levels. Isotopically depleted groundwater did not have a significant influence on the delta D values of n-alkanes in Lake Stechlin P. pectinatus and calculated isotopic fractionation factors epsilon(l/w) between lake water and n-alkanes averaged -137 +/- 9%(n-C-23), -136 +/- 7%(n-C-25) and -131 +/- 6%(n-C-27), respectively. Similar epsilon values were calculated for plants from Lake Muggelsee grown in freshwater nutrient solution (-134 +/- 11% for n-C-23), while greater fractionation was observed at increased salinity values of 10 (163 +/- 12%) and 15(-172 +/- 15%). We therefore suggest an average e value of -136 +/- 9% between source water and the major n-alkanes in P. pectinatus grown under freshwater conditions. Our results demonstrate that isotopic fractionation can increase by 30-40% at salinity values 10 and 15. These results could be explained either by inhibited plant growth at higher salinity, or by metabolic adaptation to salt stress that remain to be elucidated. A potential salinity effect on dD values of aquatic lipids requires further examination, since this would impact on the interpretation of downcore isotopic data in paleohydrologic studies. (C) 2017 Elsevier Ltd. All rights reserved.
Bumps in river profiles: uncertainty assessment and smoothing using quantile regression techniques
(2017)
The analysis of longitudinal river profiles is an important tool for studying landscape evolution. However, characterizing river profiles based on digital elevation models (DEMs) suffers from errors and artifacts that particularly prevail along valley bottoms. The aim of this study is to characterize uncertainties that arise from the analysis of river profiles derived from different, near-globally available DEMs. We devised new algorithms quantile carving and the CRS algorithm - that rely on quantile regression to enable hydrological correction and the uncertainty quantification of river profiles. We find that globally available DEMs commonly overestimate river elevations in steep topography. The distributions of elevation errors become increasingly wider and right skewed if adjacent hillslope gradients are steep. Our analysis indicates that the AW3D DEM has the highest precision and lowest bias for the analysis of river profiles in mountainous topography. The new 12m resolution TanDEM-X DEM has a very low precision, most likely due to the combined effect of steep valley walls and the presence of water surfaces in valley bottoms. Compared to the conventional approaches of carving and filling, we find that our new approach is able to reduce the elevation bias and errors in longitudinal river profiles.
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.
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.
Nitrogen lipid regulator (NlpR) is a pleiotropic regulator that positively controls genes associated with both nitrogen and lipid metabolism in the oleaginous bacterium Rhodococcus jostii RHA1. In this study, we investigated the effect of nlpR disruption and overexpression on the assimilation of C-13-labeled glucose as carbon source, during cultivation of cells under nitrogen-limiting and nitrogen-rich conditions, respectively. Label incorporation into the total lipid extract (TLE) fraction was about 30% lower in the mutant strain in comparison with the wild type strain under low-nitrogen conditions. Moreover, a higher C-13 abundance (similar to 60%) into the extracellular polymeric substance fraction was observed in the mutant strain, nlpR disruption also promoted a decrease in the label incorporation into several TLE-derivative fractions including neutral lipids (NL), glycolipids (GL), phospholipids (PL), triacylglycerols (TAG), diacylglycerols (DAG), and free fatty acids (FFA), with the DAG being the most affected. In contrast, the nlpR overexpression in RHA1 cells under nitrogen-rich conditions produced an increase of the label incorporation into the TLE and its derivative NL and PL fractions, the last one being the highest C-13 enriched. In addition, a higher C-13 enrichment occurred in the TAG, DAG, and FFA fractions after nlpR induction, with the FFA fraction being the most affected within the TLE. Isotopic-labeling experiments demonstrated that NlpR regulator is contributing in oleaginous phenotype of R. jostii RHA1 to the allocation of carbon into the different lipid fractions in response to nitrogen levels, increasing the rate of carbon flux into lipid metabolism.
We present a study of the local atomic environment of the oxygen atoms in the aqueous solutions of NaOH and HCl under simultaneous high-temperature and high-pressure conditions. Experimental nonresonant X-ray Raman scattering core-level spectra at the oxygen K-edge show systematic changes as a function of temperature and pressure. These systematic changes are distinct for the two different solutes and are described well by calculations within the Bethe- Salpeter formalism for snapshots from ab initio molecular dynamics simulations. The agreement between experimental and simulation results allows us to use the computations for a detailed fingerprinting analysis in an effort to elucidate the local atomic structure and hydrogen-bonding topology in these relevant solutions. We observe that both electrolytes, especially NaOH, enhance hydrogen bonding and tetrahedrality in the water structure at supercritical conditions, in particular in the vicinity of the hydration shells. This effect is accompanied with the association of the HCl and NaOH molecules at elevated temperatures.
Erosion and tectonic uplift are widely thought to be coupled through feedbacks involving orographic precipitation, relief development, and crustal weakening. In many orogenic systems, it can be difficult to distinguish whether true feedbacks exist, or whether observed features are a consequence of tectonic forcing. To help elucidate these interactions, we examine Gongga Shan, a 7556 m peak on the eastern margin of the Tibetan Plateau where cosmogenic Be-10 basin-wide erosion rates reach >5 mm/yr, defining a region of localized rapid erosion associated with a restraining bend in the left-lateral Xianshuihe Fault. Erosion rates are consistent with topography, thermochronometry, and geodetic data, suggesting a stable pattern of uplift and exhumation over at least the past 2-3 My. Transpression along the Xianshuihe Fault, orographically enhanced precipitation, thermally weakened crust, and substantial local relief all developed independently in the Gongga region and existed there prior to the uplift of Gongga Shan. However, only where all of these conditions are present do the observed topographic and erosional extremes exist, and their relative timing indicates that these conditions are not a consequence of rapid uplift. We conclude that their collocation at 3-4 Ma set into motion a series of feedbacks between erosion and uplift that has resulted in the exceptionally high topography and rapid erosion rates observed today.
Glacial periods are recognized to be dustier than interglacials, but the conditions leading to greater dust mobilization are poorly defined. Here we present a new high-resolution dust record based on Th-230-normalized He-4 flux from Ocean Drilling Program site 882 in the Subarctic North Pacific covering the last 170,000years. By analogy with modern relationships, we infer the mechanisms controlling orbital-scale dust storm variability in East Asia. We propose that orbital-scale dust flux variability is the result of an expansion of the dust season into summer, in addition to more intense dust storms during spring and fall. The primary drivers influencing dust flux include summer insolation at subarctic latitudes and variable Siberian alpine glaciation, which together control the cold air reservoir in Siberia. Changes in the extent of the Northern Hemisphere ice sheets may be a secondary control.
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.
Northeastern (NE) India experiences extraordinarily pronounced seasonal climate, governed by the Indian summer monsoon (ISM). The vulnerability of this region to floods and droughts calls for detailed and highly resolved paleoclimate reconstructions to assess the recurrence rate and driving factors of ISM changes. We use stable oxygen and carbon isotope ratios (delta O-18 and delta C-13) from stalagmite MAW-6 from Mawmluh Cave to infer climate and environmental conditions in NE India over the last deglaciation (16-6ka). We interpret stalagmite delta O-18 as reflecting ISM strength, whereas delta C-13 appears to be driven by local hydroclimate conditions. Pronounced shifts in ISM strength over the deglaciation are apparent from the delta O-18 record, similarly to other records from monsoonal Asia. The ISM is weaker during the late glacial (LG) period and the Younger Dryas, and stronger during the BOlling-Allerod and Holocene. Local conditions inferred from the delta C-13 record appear to have changed less substantially over time, possibly related to the masking effect of changing precipitation seasonality. Time series analysis of the delta O-18 record reveals more chaotic conditions during the late glacial and higher predictability during the Holocene, likely related to the strengthening of the seasonal recurrence of the ISM with the onset of the Holocene.
Climatic and limnological changes at Lake Karakul (Tajikistan) during the last similar to 29 cal ka
(2017)
We present results of analyses on a sediment core from Lake Karakul, located in the eastern Pamir Mountains, Tajikistan. The core spans the last similar to 29 cal ka. We investigated and assessed processes internal and external to the lake to infer changes in past moisture availability. Among the variables used to infer lake-external processes, high values of grain-size end-member (EM) 3 (wide grain-size distribution that reflects fluvial input) and high Sr/Rb and Zr/Rb ratios (coinciding with coarse grain sizes), are indicative of moister conditions. High values in EM1, EM2 (peaks of small grain sizes that reflect long-distance dust transport or fine, glacially derived clastic input) and TiO2 (terrigenous input) are thought to reflect greater influence of dry air masses, most likely of Westerly origin. High input of dust from distant sources, beginning before the Last Glacial Maximum (LGM) and continuing to the late glacial, reflects the influence of dry Westerlies, whereas peaks in fluvial input suggest increased moisture availability. The early to early-middle Holocene is characterised by coarse mean grain sizes, indicating constant, high fluvial input and moister conditions in the region. A steady increase in terrigenous dust and a decrease in fluvial input from 6.6 cal ka BP onwards points to the Westerlies as the predominant atmospheric circulation through to present, and marks a return to drier and even arid conditions in the area. Proxies for productivity (TOC, TOC/TN, TOCBr), redox potential (Fe/Mn) and changes in the endogenic carbonate precipitation (TIC, delta(18) OCarb) indicate changes within the lake. Low productivity characterised the lake from the late Pleistocene until 6.6 cal ka BP, and increased rapidly afterwards. Lake level remained low until the LGM, but water depth increased to a maximum during the late glacial and remained high into the early Holocene. Subsequently, the water level decreased to its present stage. Today the lake system is mainly climatically controlled, but the depositional regime is also driven by internal limnogeological processes.
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.
Sensitivity and identifiability analyses are common diagnostic tools to address over-parametrization in complex environmental models, but a combined application of the two analyses is rarely conducted. In this study, we performed a temporal global sensitivity analysis using the variance-based method of Sobol’ and a temporal identifiability analysis of model parameters using the dynamic identifiability method (DYNIA). We discuss the relationship between the two analyses with a focus on parameter identification and output uncertainty reduction. The hydrological model HydroGeoSphere was used to simulate daily evapotranspiration, water content, and seepage at the lysimeter scale. We found that identifiability of a parameter does not necessarily reduce output uncertainty. It was also found that the information from the main and total effects (main Sobol' sensitivity indices) is required to allow uncertainty reduction in the model output. Overall, the study highlights the role of combined temporal diagnostic tools for improving our understanding of model behavior.
Long-term estimates of the dissolution of silicate rock are generally derived from a range of isotopic proxies, such as the radiogenic strontium isotope ratio (Sr-87/Sr-86), which are preserved in sediment archives. For these systems to fairly represent silicate weathering, the changes in isotopic ratios in terrestrial surface waters should correspond to changes in the overall silicate dissolution. This assumes that the silicate mineral phases that act as sources of a given isotope dissolve at a rate that is proportional to the overall silicate weathering. Bedrock landsliding exhumes large quantities of fresh rock for weathering in transient storage, and rapid weathering in these deposits is controlled primarily by dissolution of the most reactive phases. In this study, we test the hypothesis that preferential weathering of these labile minerals can decouple the dissolution of strontium sources from the actual silicate weathering rates in the rapidly eroding Western Southern Alps (WSA) of New Zealand. We find that rapid dissolution of relatively radiogenic calcite and biotite in landslides leads to high local fluxes in strontium with isotopic ratios that offer no clear discrimination between sources. These higher fluxes of radiogenic strontium are in contrast to silicate weathering rates in landslides that are not systematically elevated. On a mountain belt scale, radiogenic strontium fluxes are not coupled to volumes of recent landslides in large (>100 km(2)) catchments, but silicate weathering fluxes are. Such decoupling is likely due first to the broad variability in the strontium content of carbonate minerals, and second to the combination of radiogenic strontium released from both biotite and carbonate in recent landslides. This study supports previous work suggesting the limited utility of strontium isotopes as a system to study silicate weathering in the WSA. Crucially however, in settings where bedrock landsliding is a dominant erosive process there is potential for both random and systematic bias in isotope proxies if the most reactive phases exposed for dissolution by landslides disproportionately contribute to the proxy of choice. This clearly suggests that the isotopic composition of marine Sr is a proxy for periods of rapid mountain uplift and erosion rather than for the associated enhanced silicate weathering. (C) 2017 Elsevier Ltd. All rights reserved.
To better assess ecosystem C budgets of croplands and understand their potential response to climate and management changes, detailed information on the mechanisms and environmental controls driving the individual C flux components are needed. This accounts in particular for the ecosystem respiration (R-eco) and its components, the autotrophic (R-a) and heterotrophic respiration (R-h) which vary tremendously in time and space. This study presents a method to separate R-eco into R-a [as the sum of R-a (shoot) and R-a (root)] and R-h in order to detect temporal and small-scale spatial dynamics within their relative contribution to overall R-eco. Thus, predominant environmental drivers and underlying mechanisms can be revealed. R-eco was derived during nighttime by automatic chamber CO2 flux measurements on plant covered plots. R-h was derived from CO2 efflux measurements, which were performed in parallel to R-eco measurements on a fallow plot using CO2 sampling tubes in 10 cm soil depth. R-a (root) was calculated as the difference between sampling tube CO2 efflux measurements on a plant covered plot and R-h. R-a (shoot) was calculated as R-eco - R-a (root) - R-h. Measurements were carried out for winter wheat (Triticum aestivum L.) during the crop season 2015 at an experimental plot located in the hummocky ground moraine landscape of NE Germany. R-eco varied seasonally from < 1 to 9.5 g C m(-2) d(-1), and was higher in adult (a) and reproductive (r) than juvenile (j) stands (gC m(-2) d(-1): j = 1.2, a = 4.6, r = 5.3). Observed R-a and R-h were in general smaller compared to the independently measured R-eco, contributing in average 58% and 42% to R-eco. However, both varied strongly regarding their environmental drivers and particular contribution throughout the study period, following the seasonal development of soil temperature and moisture (R-h) as well as crop development (R-a). Thus, our results consistently revealed temporal dynamics regarding the relative contribution of R-a (root) and R-a (shoot) to R-a, as well as of R-a and R-h to R-eco. Based on the observed results, implications for partitioning of R-eco in croplands are given.
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.
Water keeps puzzling scientists because of its numerous properties which behave oppositely to those of usual liquids: for instance, water expands upon cooling, and liquid water is denser than ice. To explain this anomalous behavior, several theories have been proposed, with different predictions for the properties of supercooled water (liquid at conditions where ice is stable). However, discriminating between those theories with experiments has remained elusive because of spontaneous ice nucleation. Here we measure the sound velocity in liquid water stretched to negative pressure and derive an experimental equation of state, which reveals compressibility anomalies. We show by rigorous thermodynamic relations how these anomalies are intricately linked with the density anomaly. Some features we observe are necessary conditions for the validity of two theories of water.
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.
In 2002, a severe flood caused financial losses of EUR 11.6 billion in Germany and triggered many changes in flood risk management. This chapter focuses on flood insurance, which is a voluntary supplementary insurance in Germany: it is explored how flood insurance has contributed to enhance resilience of flood-prone residents. The analyses are based on empirical data collected by post-event surveys in the federal states of Saxony and Bavaria and refer to the three pillars upon which the concept of flood resilience usually builds in the natural hazards context: recovery, adaptive capacity and resistance. Overall, the penetration of flood insurance has increased since 2002 and there is strong empirical evidence that losses of insured residents are more often and better compensated than those of uninsured despite the provision of governmental financial disaster assistance after big floods. This facilitation of recovery is, however, not the only contribution to flood resilience. Insured residents tend to invest more in further flood mitigation measures at their properties than uninsured. Obviously, flood insurance is embedded in a complex safety strategy of property owners that needs more investigation in order to be addressed more effectively in risk communication and integrated risk management strategies.
Mountain channels can be strongly coupled with adjacent hillslopes, exchanging both mass and energy. However, hypotheses of the underlying cause and effect relations are based on indirect observations that do not resolve the mechanics of channel-hillslope coupling at the process scale. Here we present direct observational data of a coupled channel-hillslope system in the catchment area of the Erlenbach, a mountain stream in Switzerland. A slow-moving landslide flanking this alpine stream failed after a flood had eroded an alluvial step in the channel at its base, representing evidence for an upsystem link in channel-hillslope coupling. Progressive accumulation of landslide debris in the channel resulted in a renewed step, stabilizing the hillslope and restoring the channel step in a downsystem link. Thus, upsystem and downsystem coupling mechanisms are joined in a negative feedback cycle. In this cycle, debuttressing and rebuttressing due to channel bed erosion and alluviation are the dominant controls on hillslope stability. Based on an order of magnitude estimate it is plausible that the observed feedback mechanism is a relevant process in the production of coarse (>2 mm) sediment in the Erlenbach.
Inherited rheological structures in the lithosphere are expected to have large impact on the architecture of continental rifts. The Turkana depression in the East African Rift connects the Main Ethiopian Rift to the north with the Kenya rift in the south. This region is characterized by a NW-SE trending band of thinned crust inherited from a Mesozoic rifting event, which is cutting the present-day N-S rift trend at high angle. In striking contrast to the narrow rifts in Ethiopia and Kenya, extension in the Turkana region is accommodated in subparallel deformation domains that are laterally distributed over several hundred kilometers. We present both analog experiments and numerical models that reproduce the along-axis transition from narrow rifting in Ethiopia and Kenya to a distributed deformation within the Turkana depression. Similarly to natural observations, our models show that the Ethiopian and Kenyan rifts bend away from each other within the Turkana region, thus forming a right-lateral step over and avoiding a direct link to form a continuous N-S depression. The models reveal five potential types of rift linkage across the preexisting basin: three types where rifts bend away from the inherited structure connecting via a (1) wide or (2) narrow rift or by (3) forming a rotating microplate, (4) a type where rifts bend towards it, and (5) straight rift linkage. The fact that linkage type 1 is realized in the Turkana region provides new insights on the rheological configuration of the Mesozoic rift system at the onset of the recent rift episode. Plain Language Summary The Turkana depression in the Kenya/Ethiopia borderland is most famous for its several million years old human fossils, but it also holds a rich geological history of continental separation. The Turkana region is a lowland located between the East African and Ethiopian domes because its crust and mantle have been stretched in a continent-wide rift event during Cretaceous times about 140-120 Ma ago. This thin lithosphere exerted paramount control on the dynamics of East African rifting in this area, which commenced around 15 Ma ago and persists until today. Combining analog "sandbox" experiments with numerical geodynamic modeling, we find that inherited rift structures explain the dramatic change in rift style from deep, narrow rift basins north and south of the Turkana area to wide, distributed deformation within the Turkana depression. The failed Cretaceous rift is also responsible for the eastward bend of the Ethiopian rift and the westward bend of the Kenyan rift when entering the Turkana depression, which generated the characteristic hook-shaped form of present-day Lake Turkana. Combing two independent modeling techniques-analog and numerical experiments-is a very promising approach allowing to draw robust conclusions about the processes that shape the surface of our planet.
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.
Flow patterns in conjunction with seasonal and diurnal temperature variations control ecological and biogeochemical conditions in hyporheic sediments. In particular, hyporheic temperatures have a great impact on many temperature-sensitive microbial processes. In this study, we used 3-D coupled water flow and heat transport simulations applying the HydroGeoSphere code in combination with high-resolution observations of hydraulic heads and temperatures to quantify reach-scale water and heat flux across the river-groundwater interface and hyporheic temperature dynamics of a lowland gravel bed river. The model was calibrated in order to constrain estimates of the most sensitive model parameters. The magnitude and variations of the simulated temperatures matched the observed ones, with an average mean absolute error of 0.7 degrees C and an average Nash Sutcliffe efficiency of 0.87. Our results indicate that nonsubmerged streambed structures such as gravel bars cause substantial thermal heterogeneity within the saturated sediment at the reach scale. Individual hyporheic flow path temperatures strongly depend on the flow path residence time, flow path depth, river, and groundwater temperature. Variations in individual hyporheic flow path temperatures were up to 7.9 degrees C, significantly higher than the daily average (2.8 degrees C), but still lower than the average seasonal hyporheic temperature difference (19.2 degrees C). The distribution between flow path temperatures and residence times follows a power law relationship with exponent of about 0.37. Based on this empirical relation, we further estimated the influence of hyporheic flow path residence time and temperature on oxygen consumption which was found to partly increase by up to 29% in simulations.
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.
Permafrost is a distinct feature of the terrestrial Arctic and is vulnerable to climate warming. Permafrost degrades in different ways, including deepening of a seasonally unfrozen surface and localized but rapid development of deep thaw features. Pleistocene ice-rich permafrost with syngenetic ice-wedges, termed Yedoma deposits, are widespread in Siberia, Alaska, and Yukon, Canada and may be especially prone to rapid-thaw processes. Freeze-locked organic matter in such deposits can be re-mobilized on short time-scales and contribute to a carbon-cycle climate feedback. Here we synthesize the characteristics and vulnerability of Yedoma deposits by synthesizing studies on the Yedoma origin and the associated organic carbon pool. We suggest that Yedoma deposits accumulated under periglacial weathering, transport, and deposition dynamics in non-glaciated regions during the late Pleistocene until the beginning of late glacial warming. The deposits formed due to a combination of aeolian, colluvial, nival, and alluvial deposition and simultaneous ground ice accumulation. We found up to 130 gigatons organic carbon in Yedoma, parts of which are well-preserved and available for fast decomposition after thaw. Based on incubation experiments, up to 10% of the Yedoma carbon is considered especially decomposable and may be released upon thaw. The substantial amount of ground ice in Yedoma makes it highly vulnerable to disturbances such as thermokarst and thermo-erosion processes. Mobilization of permafrost carbon is expected to increase under future climate warming. Our synthesis results underline the need of accounting for Yedoma carbon stocks in next generation Earth-System-Models for a more complete representation of the permafrost-carbon feedback.
Seismic-hazard assessment is of great importance within the field of engineering seismology. Nowadays, it is common practice to define future seismic demands using probabilistic seismic-hazard analysis (PSHA). Often it is neither obvious nor transparent how PSHA responds to changes in its inputs. In addition, PSHA relies on many uncertain inputs. Sensitivity analysis (SA) is concerned with the assessment and quantification of how changes in the model inputs affect the model response and how input uncertainties influence the distribution of the model response. Sensitivity studies are challenging primarily for computational reasons; hence, the development of efficient methods is of major importance. Powerful local (deterministic) methods widely used in other fields can make SA feasible, even for complex models with a large number of inputs; for example, automatic/algorithmic differentiation (AD)-based adjoint methods. Recently developed derivative-based global sensitivity measures can combine the advantages of such local SA methods with efficient sampling strategies facilitating quantitative global sensitivity analysis (GSA) for complex models. In our study, we propose and implement exactly this combination. It allows an upper bounding of the sensitivities involved in PSHA globally and, therefore, an identification of the noninfluential and the most important uncertain inputs. To the best of our knowledge, it is the first time that derivative-based GSA measures are combined with AD in practice. In addition, we show that first-order uncertainty propagation using the delta method can give satisfactory approximations of global sensitivity measures and allow a rough characterization of the model output distribution in the case of PSHA. An illustrative example is shown for the suggested derivative-based GSA of a PSHA that uses stochastic ground-motion simulations.
Carbon (C) sequestration in soils plays a key role in the global C cycle. It is therefore crucial to adequately monitor dynamics in soil organic carbon (Delta SOC) stocks when aiming to reveal underlying processes and potential drivers. However, small-scale spatial (10-30 m) and temporal changes in SOC stocks, particularly pronounced in arable lands, are hard to assess. The main reasons for this are limitations of the well-established methods. On the one hand, repeated soil inventories, often used in long-term field trials, reveal spatial patterns and trends in Delta SOC but require a longer observation period and a sufficient number of repetitions. On the other hand, eddy covariance measurements of C fluxes towards a complete C budget of the soil-plant-atmosphere system may help to obtain temporal Delta SOC patterns but lack small-scale spatial resolution. To overcome these limitations, this study presents a reliable method to detect both short-term temporal dynamics as well as small-scale spatial differences of Delta SOC using measurements of the net ecosystem carbon balance (NECB) as a proxy. To estimate the NECB, a combination of automatic chamber (AC) measurements of CO2 exchange and empirically modeled aboveground biomass development (NPPshoot / were used. To verify our method, results were compared with Delta SOC observed by soil resampling. Soil resampling and AC measurements were performed from 2010 to 2014 at a colluvial depression located in the hummocky ground moraine landscape of northeastern Germany. The measurement site is characterized by a variable groundwater level (GWL) and pronounced small-scale spatial heterogeneity regarding SOC and nitrogen (Nt) stocks. Tendencies and magnitude of Delta SOC values derived by AC measurements and repeated soil inventories corresponded well. The period of maximum plant growth was identified as being most important for the development of spatial differences in annual Delta SOC. Hence, we were able to confirm that AC-based C budgets are able to reveal small-scale spatial differences and short-term temporal dynamics of Delta SOC.
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.
Am Abend des 29. Mai 2016 wurde der Ort Braunsbach im Landkreis Schwäbisch-Hall (Baden-Württemberg) von einer Sturzflut getroffen, bei der mehrere Häuser stark beschädigt oder zerstört wurden. Die Sturzflut war eine der Unwetterfolgen, die im Frühsommer 2016 vom Tiefdruckgebiet Elvira ausgelöst wurden. Der vorliegende Bericht ist der zweite Teil einer Doppelveröffentlichung, welche die Ergebnisse zur Untersuchung des Sturzflutereignisses im Rahmen des DFG-Graduiertenkollegs “Naturgefahren und Risiken in einer sich verändernden Welt” (NatRiskChange, GRK 2043/1) der Universität Potsdam präsentiert. Während Teil 1 die meteorologischen und hydrologischen Ereignisse analysiert, fokussiert Teil 2 auf die geomorphologischen Prozesse und die verursachten Gebäudeschäden. Dazu wurden Ursprung und Ausmaß des während des Sturzflutereignisses mobilisierten und in den Ort getragenen Materials untersucht. Des Weiteren wurden zu 96 betroffenen Gebäuden Daten zum Schadensgrad sowie Prozess- und Gebäudecharakteristika aufgenommen und ausgewertet. Die Untersuchungen zeigen, dass bei der Betrachtung von Hochwassergefährdung die Berücksichtigung von Sturzfluten und ihrer speziellen Charakteristika, wie hoher Feststofftransport und sprunghaftes Verhalten insbesondere in bebautem Gelände, wesentlich ist, um effektive Schutzmaßnahmen ergreifen zu können.
Disentangling shallow‐water bulk carbonate carbon isotope archives into primary and diagenetic components is a notoriously difficult task and even diagenetically screened records often provide chemostratigraphic patterns that significantly differ from global signals. This is mainly caused by the polygenetic nature of shallow‐water carbonate substrates, local carbon cycle processes causing considerable neritic–pelagic isotope gradients and the presence of hiatal surfaces resulting in extremely low carbonate preservation rates. Provided here is an in‐depth petrographic and geochemical evaluation of different carbonate phases of a mid‐Cretaceous (Barremian–Aptian) shallow‐water limestone succession (Jabal Madar section) deposited on the tropical Arabian carbonate platform in Oman. The superposition of stable isotope signatures of identified carbonate phases causes a complex and often noisy bulk carbon isotope pattern. Blocky sparite cements filling intergranular pores and bioclastic voids evidence intermediate to (arguably) deep burial diagenetic conditions during their formation, owing to different timing or differential faulting promoting the circulation of fluids from variable sources. In contrast, sparite cements filling sub‐vertical veins reveal a rock‐buffered diagenetic fluid composition with an intriguing moderate enrichment in 13C, probably due to fractionation during pressure release in the context of the Miocene exhumation of the carbonate platform under study. The presence of abundant, replacive dedolomite in mud‐supported limestone samples forced negative carbon and oxygen isotope changes that are either associated with the thermal breakdown of organic matter in the deep burial realm or the expulsion of buried meteoric water in the intermediate burial realm. Notwithstanding the documented stratigraphically variable and often facies‐related impact of different diagenetic fluids on the bulk‐rock stable isotope signature, the identification of diagenetic end‐members defined δ13C and δ18O threshold values that allowed the most reliable ‘primary’ bulk carbon isotope signatures to be extracted. Most importantly, this approach exemplifies how to place regional shallow‐water stable isotope patterns with evidence for a complex multi‐stage diagenetic history into a supraregional or even global context.
Arctic and alpine treelines worldwide differ in their reactions to climate change. A northward advance of or densification within the treeline ecotone will likely influence climate-vegetation feedback mechanisms. In our study, which was conducted in the Taimyr Depression in the North Siberian Lowlands, w present a combined field-and model-based approach helping us to better understand the population processes involved in the responses of the whole treeline ecotone, spanning from closed forest to single-tree tundra, to climate warming. Using information on stand structure, tree age, and seed quality and quantity from seven sites, we investigate effects of intra-specific competition and seed availability on the specific impact of recent climate warming on larch stands. Field data show that tree density is highest in the forest-tundra, and average tree size decreases from closed forest to single-tree tundra. Age-structure analyses indicate that the trees in the closed forest and forest-tundra have been present for at least similar to 240 yr. At all sites except the most southerly ones, past establishment is positively correlated with regional temperature increase. In the single-tree tundra, however, a change in growth form from krummholz to erect trees, beginning similar to 130 yr ago, rather than establishment date has been recorded. Seed mass decreases from south to north, while seed quantity increases. Simulations with LAVESI (Larix Vegetation Simulator) further suggest that relative density changes strongly in response to a warming signal in the forest-tundra while intra-specific competition limits densification in the closed forest and seed limitation hinders densification in the single-tree tundra. We find striking differences in strength and timing of responses to recent climate warming. While forest-tundra stands recently densified, recruitment is almost non-existent at the southern and northern end of the ecotone due to autecological processes. Palaeo-treelines may therefore be inappropriate to infer past temperature changes at a fine scale. Moreover, a lagged treeline response to past warming will, via feedback mechanisms, influence climate change in the future.
Dissolved CO2 storage in geological formations with low pressure, low risk and large capacities
(2017)
Geological CO2 storage is a mitigation technology to reduce CO2 emissions from fossil fuel combustion. However, major concerns are the pressure increase and saltwater displacement in the mainly targeted deep groundwater aquifers due to injection of supercritical CO2. The suggested solution is storage of CO2 exclusively in the dissolved state. In our exemplary regional case study of the North East German Basin based on a highly resolved temperature and pressure distribution model and a newly developed reactive transport coupling, we have quantified that 4.7 Gt of CO2 can be stored in solution compared to 1.5 Gt in the supercritical state.
This special issue is the result of several fruitful conference sessions on disturbance hydrology, which started at the 2013 AGU Fall Meeting in San Francisco and have continued every year since. The stimulating presentations and discussions surrounding those sessions have focused on understanding both the disruption of hydrologic functioning following discrete disturbances, as well as the subsequent recovery or change within the affected watershed system. Whereas some hydrologic disturbances are directly linked to anthropogenic activities, such as resource extraction, the contributions to this special issue focus primarily on those with indirect or less pronounced human involvement, such as bark-beetle infestation, wildfire, and other natural hazards. However, human activities are enhancing the severity and frequency of these seemingly natural disturbances, thereby contributing to acute hydrologic problems and hazards. Major research challenges for our increasingly disturbed planet include the lack of continuous pre and postdisturbance monitoring, hydrologic impacts that vary spatially and temporally based on environmental and hydroclimatic conditions, and the preponderance of overlapping or compounding disturbance sequences. In addition, a conceptual framework for characterizing commonalities and differences among hydrologic disturbances is still in its infancy. In this introduction to the special issue, we advance the fusion of concepts and terminology from ecology and hydrology to begin filling this gap. We briefly explore some preliminary approaches for comparing different disturbances and their hydrologic impacts, which provides a starting point for further dialogue and research progress.
Tree stands in the boreal treeline ecotone are, in addition to climate change, impacted by disturbances such as fire, water-related disturbances and logging. We aim to understand how these disturbances affect growth, age structure, and spatial patterns of larch stands in the north-eastern Siberian treeline ecotone (lower Kolyma River region), an insufficiently researched region. Stand structure of Larix cajanderi Mayr was studied at seven sites impacted by disturbances. Maximum tree age ranged from 44 to 300 years. Young to medium-aged stands had, independent of disturbance type, the highest stand densities with over 4000 larch trees per ha. These sites also had the highest growth rates for tree height and stem diameter. Overall lowest stand densities were found in a polygonal field at the northern end of the study area, with larches growing in distinct " tree islands". At all sites, saplings are significantly clustered. Differences in fire severity led to contrasting stand structures with respect to tree, recruit, and overall stand densities. While a low severity fire resulted in low-density stands with high proportions of small and young larches, high severity fires resulted in high-density stands with high proportions of big trees. At waterdisturbed sites, stand structure varied between waterlogged and drained sites and latitude. These mixed effects of climate and disturbance make it difficult to predict future stand characteristics and the treeline position.
River restoration is a main emphasis of river management in European countries. Cross-national comparisons of its implementation are still rare in scientific literature. Based on French and German national censuses, this study compares river restoration practices and monitoring by analysing 102 French and 270 German projects. This comparison aims to draw a spatial and temporal framework of restoration practices in both countries to identify potential drivers of cross-national similarities and differences. The results underline four major trends: (1) a lag of almost 15 years in river restoration implementation between France and Germany, with a consequently higher share of projects in Germany than in France, (2) substantial similarities in restored reach characteristics, short reach length, small rivers, and in "agricultural" areas, (3) good correspondences between stressors identified and restoration measures implemented. Morphological alterations were the most important highlighted stressors. River morphology enhancement, especially instream enhancements, were the most frequently implemented restoration measures. Some differences exist in specific restoration practices, as river continuity restoration were most frequently implemented in French projects, while large wood introduction or channel re-braiding were most frequently implemented in German projects, and (4) some quantitative and qualitative differences in monitoring practices and a significant lack of project monitoring, especially in Germany compared to France. These similarities and differences between Germany and France in restoration application and monitoring possibly result from a complex set of drivers that might be difficult to untangle (e.g., environmental, technical, political, cultural).
High precipitation quantiles tend to rise with temperature, following the so-called Clausius-Clapeyron (CC) scaling. It is often reported that the CC-scaling relation breaks down and even reverts for very high temperatures. In our study, we investigate this reversal using observational climate data from 142 stations across Germany. One of the suggested meteorological explanations for the breakdown is limited moisture supply. Here we argue that, instead, it could simply originate from undersampling. As rainfall frequency generally decreases with higher temperatures, rainfall intensities as dictated by CC scaling are less likely to be recorded than for moderate temperatures. Empirical quantiles are conventionally estimated from order statistics via various forms of plotting position formulas. They have in common that their largest representable return period is given by the sample size. In small samples, high quantiles are underestimated accordingly. The small-sample effect is weaker, or disappears completely, when using parametric quantile estimates from a generalized Pareto distribution (GPD) fitted with L moments. For those, we obtain quantiles of rainfall intensities that continue to rise with temperature.
High precipitation quantiles tend to rise with temperature, following the so-called Clausius–Clapeyron (CC) scaling. It is often reported that the CC-scaling relation breaks down and even reverts for very high temperatures. In our study, we investigate this reversal using observational climate data from 142 stations across Germany. One of the suggested meteorological explanations for the breakdown is limited moisture supply. Here we argue that, instead, it could simply originate from undersampling. As rainfall frequency generally decreases with higher temperatures, rainfall intensities as dictated by CC scaling are less likely to be recorded than for moderate temperatures. Empirical quantiles are conventionally estimated from order statistics via various forms of plotting position formulas. They have in common that their largest representable return period is given by the sample size. In small samples, high quantiles are underestimated accordingly. The small-sample effect is weaker, or disappears completely, when using parametric quantile estimates from a generalized Pareto distribution (GPD) fitted with L moments. For those, we obtain quantiles of rainfall intensities that continue to rise with temperature.
Effects of topographic smoothing on the simulation of winter precipitation in High Mountain Asia
(2017)
Numerous studies have projected future changes in High Mountain Asia water resources based on temperature and precipitation from global circulation models (GCMs) under future climate scenarios. Although the potential benefit of such studies is immense, coarse grid-scale GCMs are unable to resolve High Mountain Asia's complex topography and thus have a biased representation of regional weather and climate. This study investigates biases in the simulation of physical mechanisms that generate snowfall and contribute to snowpack in High Mountain Asia in coarse topography experiments using the Weather Research and Forecasting model. Regional snowpack is event driven, thus 33 extreme winter orographic precipitation events are simulated at fine atmospheric resolution with 6.67 km resolution topography and smoothed 1.85° × 1.25° GCM topography. As with many modified topography experiments performed in other regions, the distribution of precipitation is highly dependent on first-order orographic effects, which dominate regional meteorology. However, we demonstrate that topographic smoothing enhances circulation in simulated extratropical cyclones, with significant impacts on orographic precipitation. Despite precipitation reductions of 28% over the highest ranges, due to reduced ascent on windward slopes, total precipitation over the study domain increased by an average of 9% in smoothed topography experiments on account of intensified extratropical cyclone dynamics and cross-barrier moisture flux. These findings identify an important source of bias in coarse-resolution simulated precipitation in High Mountain Asia, with important implications for the application of GCMs toward projecting future hydroclimate in the region.
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.
Understanding evolution and ecology of endemic plants is of great importance for conservation of those rare and endangered species. Pairwise comparisons of plant functional traits could be an adequate method to get insights in evolutionary and ecological processes. We examined whether morphological traits representing competitive ability and habitat specificity differ between endemics and common plants. Therefore, we performed pairwise comparison analyses of 9 plant functional traits in 36 congeneric pairs of endemics and their common congeners on the East Mediterranean island of Cyprus, i.e., the first such study conducted on a Mediterranean island. We found that endemic species prefer higher elevations and more extreme habitats. Endemics were smaller and they had smaller flowers than their common congeners. Common species had higher chromosome numbers than endemic ones. Endemic and common species showed no significant differences in canopy height, inflorescence height, leaf length and width, and flowering period. Our study showed that the situation on a large oceanic island does not differ from results in mainland research areas.
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.
An understanding of the relationship between canopy structure and the water balance is needed for predicting how forest structure changes affect rainfall partitioning and, consequently, water resources. The objective of this study was to predict rainfall interception (I) and canopy storage capacity (S) using canopy structure variables and to investigate how seasonal changes influence their relationship. The study was conducted in twelve 50 m x 50 m plots in the Zagros forest in the western Iranian state of Ilam, protected forests of Dalab region. Average cumulative I was 84.2mm, accounting for 10.2% of cumulative gross precipitation (GP) over a 1-year period. Using a regression based method, S averaged similar to 1 mm and 0.1 mm in the leafed and leafless periods, respectively. There were no relationships between tree density and I: GP or S, but I: GP and S increased with leaf area index, canopy cover fraction, basal area, tree height, and diameter at breast height in the leafed period. In addition, wood area index and canopy cover fraction were related to I: GP or S in the leafless period. (C) 2017 Elsevier B.V. All rights reserved.
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 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.
Background. Metal recycling factories (MRFs) have developed rapidly in Nigeria as recycling policies have been increasingly embraced. These MRFs are point sources for introducing potentially toxic elements (PTEs) into environmental media. Objectives. The aim of this study was to determine the constituents (elemental and mineralogy) of the wastes (slag and particulate matter, (PM)) and soils around the MRFs and to determine the level of pollution within the area. Methods. Sixty samples (30 slag samples, 15 soil samples and 15 PM samples) were collected for this study. The soils, slag and PM samples were analyzed for elemental constituents using inductively coupled plasma optical emission spectrometry. Mineralogy of the PM was determined using scanning electron microscope-energy dispersive spectroscopy (SEM-EDS), and soil mineralogy was determined by an X-ray diffractometer (XRD). Results. The results of the soil analyses revealed the following concentrations for the selected metals in mg/kg include lead (Pb) (21.0-2399.0), zinc (Zn) (56.0-4188.0), copper (Cu) (10.0-1470.0), nickel (Ni) (6.0-215.0), chromium (Cr) (921.0-1737.0) and cadmium (Cd) (below detectable limit (Bdl)-18.1). For the slags the results were Pb (68.0-.333.0), Zn (1364.0-3062), Cu (119.0-1470.0), Ni (12.0-675.0), Cr (297-1737) and Cd (Bdl-15.8). The results in mu g/g for the metal analysis in PM were Pb (4.6-160.0), Zn (18.0-471.0), Cu (2.5-11.0), Ni (0.8-4.2), and Cr (2.5-11.0), while Cd was undetected. The slags are currently utilized for filling the foundations of buildings and roads, providing additional pathways for the introduction of PTEs into the environment from the suspended materials generated from mechanical breakdown of the slags. Conclusions. The MRFs were found to have impacted the quality of environmental media through the introduction of PTEs, impairing soil quality, in addition to PM, which can have detrimental health consequences. Further studies on the health implications of these pollutants and their impacts on human health are needed. Competing Interests. The authors declare no competing financial interests
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.
The Acheulian site of Gesher Benot Ya‘aqov (GBY) in the Upper Jordan Valley revealed important data on environment and material culture, as well as evidence for hominin behavioural and cognitive patterns documented at the margins of the Hula Palaeo-lake. A 50 m long core (GBY#2) drilled at the archaeological site has provided a long Pleistocene geological, environmental and climatological record, which expands the existing knowledge of hominin-habitat relationships. Bracketed by two basalt flows dated by 40Ar/39Ar and based on the identification of the Matuyama-Brunhes Boundary (MBB) and correlation with the GBY excavation site, the sedimentary sequence provides the climatic history around the MBB. Multi-proxy data including pollen and non-pollen palynomorphs, macro-botanical remains, molluscs and ostracods provide evidence for lake and lake-margin environments during Marine Isotope Stages (MIS) 20 and 19. Semi-moist conditions were followed by a pronounced dry phase during MIS 20, and warm and moist conditions with Quercus-Pistacia woodlands prevailed during MIS 19. In contrast to the reconstructed climate change from relatively dry to moister conditions, the depositional environment developed from an open-water lake during MIS 20 to a lake margin environment in MIS 19. Generally shallower conditions at the core site in MIS 19 resulted from the progradation of the lake shore due to the filling of the basin. Micro-charcoal analysis suggests a likelihood of human-induced fire in some parts of the core, which can be correlated with artefact-containing layers of the GBY excavation site. The Hula Palaeo-lake region provided an ideal niche for hominins and other vertebrates during global glacial-interglacial climate fluctuations at the end of the Early Pleistocene.
The crystallization kinetics of silicate liquids were studied experimentally in the system haplogranite-B-Li-H2O, at variable degrees of undercooling and variable water concentration. We investigated the kinetics of nucleation and crystallization of unseeded synthetic hydrous haplogranite with 1 wt % Li2O + 2 center dot 3 wt % B2O3 added (composition C1) and 2 wt % Li2O + 4 center dot 6 wt % B2O3 added (composition C2). Compositions C1 and C2 are simplified representative bulk compositions of Li-rich pegmatites and their highly differentiated cores, respectively. Starting water contents varied between 3 and 9 wt %. With few exceptions, the system remained water-undersaturated. About 86 isothermal runs of 1-60 days duration, grouped in 25 time series of constant temperature and initial H2O content, were carried out at temperatures from 400 to 700A degrees C at 300 MPa, corresponding to variable degrees of undercooling between the liquidus and glass transition. Viscosity measurements indicate that the glass transition for both compositions is below 400A degrees C for 3 wt % water and below 300A degrees C for 6 center dot 5 wt % water. The melts remained virtually crystal free at 400A degrees C, about 100A degrees C and 120A degrees C above the glass transition for compositions C1 and C2, respectively, in experiments up to 30 days long. This result is consistent with the existence of low-temperature, undercooled melts in the crust. At lower values of undercooling the runs crystallized partially, up to about 70% volume fraction. Undercooling and the amount of water are the main factors controlling nucleation and growth rates, and therefore textures. Minerals nucleate and grow sequentially according to mineral-specific nucleation delays. The mineral assemblage started with Li-Al stuffed quartz (in C1) and virgilite (in C2), solid-solutions between quartz and gamma-spodumene. The quartz-like phases were typically followed by spherulitic alkali feldspar-quartz intergrowths, euhedral petalite, and fine-grained muscovite. Nearly pure quartz formed as rims and replacement of metastable virgilite and stuffed quartz, in particular at the boron- and water-rich crystallization front of large feldspar or petalite. With the exception of muscovite, all minerals nucleated heterogeneously, on the capsule wall or on pre-existing minerals, and grew inwards, towards the capsule center. Experimental textures resembled the textures of zoned pegmatites, including skeletal, graphic, unidirectional, radiating, spherulitic, massive, and replacement textures. In some cases, when fluid saturation was reached, miarolitic cavities developed containing euhedral crystals. Although unidirectional growth rates appeared to slow down in time, volumetric rates for stable graphic alkali-feldspar quartz intergrowths and petalite remained constant for up to 60 days and similar to 70% crystallization. Metastable stuffed quartz and virgilite diminished in their growth rates in runs of 30 days or longer, were resorbed in the melt, and were partially replaced by second-generation quartz. Unobstructed, self-sustained crystal growth in conditions of very low nucleation density appears to be the dominant mechanism to form giant pegmatitic crystals, although experimental growth rates are much slower than predicted in nature based on conductive-cooling models.
Flood risk management in Europe and worldwide is not static but constantly in a state of flux. There has been a trend towards more integrated flood risk management in many countries. However, the initial situation and the pace and direction of change is very different in the various countries. In this paper, we will present a conceptual framework that seeks to explain why countries opt for different flood risk management portfolios. The developed framework utilises insights from a range of policy science concepts in an integrated way and considers, among others, factors such as geographical characteristics, the experience with flood disasters, as well as human behavioural aspects.
Electromagnetic induction (EMI) sensors using sufficiently low-frequency harmonic sources and sufficiently small loop separations operate in the low-induction-number (LIN) domain for a relatively wide range of background conductivity. These systems are used in diverse near-surface investigations including applications from soil sciences, hydrology, and archaeology. The special case of portable multiconfiguration EMI sensors operating at frequencies <= 20 kHz offers the possibility of using a fast linear deconvolution method to interpret multichannel data sets in three dimensions. Here, we have developed a fast 3D inversion/deconvolution method regularized with 3D smoothness constraints and formulated in the hybrid spectral-spatial domain. Compared with other linear approaches, the spectral-spatial domain formulation significantly reduces the computational cost of the processing and opens the door for real-time 3D interpretation of large data sets consisting of more than 100,000 data points. First, we test our proposed algorithm on synthetic data sets computed with the full Maxwell theory. Then, we apply our method to a real four-configuration EMI data set acquired to map the thickness of peat layers embedded in a sandy environment. For the synthetic and the field example, we compared our result with the result obtained using a standard point-by-point 1D nonlinear inversion approach. This comparison demonstrates that the proposed methodology provides superior lateral resolution compared with the 1D nonlinear inversion, at the same time significantly reducing the computational cost of the processing.
Eleven inverse isochron ages from total fusion and three from stepwise heating analyses fit the age model. Four experiments resulted in older inverse isochron ages that do not concur with the model within 2 sigma uncertainties and that deviate from 1 ka to 17 ka minimum. C-and R-type zoning are interpreted as representing growth in magma chamber cupolas, as wall mushes, or in narrow conduits. Persistent compositions of PO-type crystals and abundant surfaces recording dissolution features correspond to formation within a magma chamber. C-type zoning and R-type zoning have revealed an irregular incorporation of melt and fluid inclusions. These two types of zoning in feldspar are interpreted as preferentially contributing either heterogeneously distributed excess Ar-40 or inherited Ar-40 to the deviating Ar-40/Ar-39 ages that are discussed in this study. (C) 2017 Elsevier Ltd. All rights reserved.
We analyse current and possible future reforms of the Indian food policies for the most important staple grains, wheat and rice, within a two-commodity dynamic partial equilibrium model with stochastic shocks. The model is empirically grounded and reproduces past values well. It uses a new reduced-form approach to capture private storage dynamics. We evaluate the implementation of the National Food Security Act (NFSA) under several policy measures with the current regime as well as two scenarios with a regime change - implementation of cash transfers and deficiency payments. Implications for market fundamentals and fiscal costs are simulated in the medium term - until 2020/21. The NFSA puts a high pressure on fiscal costs and public stocks. Relying on imports with low support prices results in low fiscal costs and stable, but higher domestic and international prices, and a high risk of zero stocks. A policy strategy to manipulate procurement prices in order to maintain public stocks close to the norms leads to slightly higher fiscal costs with lower, but more volatile prices. The highest domestic price volatility occurs under a strategy which uses export bans in order to maintain sufficient public stocks. A cash-based regime can bring considerable savings and curb fiscal costs, particularly if targeted to the poor, and would leave sufficient stocks due to higher private stocks.
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).
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).
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.
The increasing numbers of recordings at individual sites allows quantification of empirical linear site-response adjustment factors (delta S2S(s)) from the ground motion prediction equation (GMPE) residuals. The delta S2S(s) are then used to linearly scale the ergodic GMPE predictions to obtain site-specific ground motion predictions in a partially non-ergodic Probabilistic Seismic Hazard Assessment (PSHA). To address key statistical and conceptual issues in the current practice, we introduce a novel empirical region-and site-specific PSHA methodology wherein, (1) site-to-site variability (phi(S2S)) is first estimated as a random-variance in a mixed-effects GMPE regression, (2) delta S2S(s) at new sites with strong motion are estimated using the a priori phi(S2S), and (3) the GMPE site-specific single-site aleatory variability sigma(ss,s) is replaced with a generic site-corrected aleatory variability sigma(0). Comparison of region- and site-specific hazard curves from our method against the traditional ergodic estimates at 225 sites in Europe and Middle East shows an approximate 50% difference in predicted ground motions over a range of hazard levels-a strong motivation to increase seismological monitoring of critical facilities and enrich regional ground motion data sets.
Through the zmir-Ankara-Erzincan and the Vardar oceans suture zones, convergence between the Eurasian and African plates played a key role in controlling Palaeogene magmatism in north-western Anatolia, northern Aegean, and eastern Balkans. LA-ICP-MS dating of U and Pb isotopes on zircon separates from the tuffs of the Harmankaya Volcanic Rocks, which are inter-fingered with the lower-middle Eocene deposits of the Gazikoy Formation to the north of the Ganos Fault and the Karaaac Formation in the Gelibolu Peninsula, yielded a late Ypresian (51Ma) age. The chemical characteristics suggest that the lavas and tuffs of the Harmankaya Volcanic Rocks are products of syn- or post-collision magmas. These volcanic rocks show also close affinities to the subduction-related magmas. In addition to the already known andesitic volcanic rocks, our field observations in Gokceada Island indicate also the existence of granitic and rhyolitic rocks (Marmaros Magmatic Assemblage). Our U-Pb zircon age data has shown that the newly discovered Marmaros granitic plutons intruded during late Oligocene (26Ma) into the deposits of the Karaaac Formation in Gokceada Island. LA-ICP-MS dating of U and Pb isotopes on zircon separates from the Marmaros rhyolitic rocks yielded a late Oligocene (26Ma) crystallization age. Geochemical characteristics indicate that the more-evolved Oligocene granitic and rhyolitic rock of the Marmaros Magmatic Assemblage possibly assimilated a greater amount of crustal material than the lower Eocene Harmankaya Volcanic Rocks. Geochemical features and age relationships suggest increasing amounts of crustal contamination and a decreasing subduction signature during the evolution of magmas in NW Turkey from the early Eocene to the Oligocene. The magmatic activity developed following the northward subduction of the zmir-Ankara-Erzincan oceanic lithosphere and the earliest Palaeocene final continental collision between the Sakarya and Anatolide-Tauride zones.