Refine
Year of publication
- 2016 (291) (remove)
Document Type
- Article (221)
- Doctoral Thesis (37)
- Other (14)
- Postprint (12)
- Review (5)
- Habilitation Thesis (1)
- Report (1)
Is part of the Bibliography
- yes (291) (remove)
Keywords
- climate change (5)
- erosion (5)
- water balance (5)
- Pollen (4)
- August 2002 flood (3)
- Body waves (3)
- Central Europe (3)
- Floods Directive (3)
- Himalaya (3)
- Himalayas (3)
Institute
- Institut für Geowissenschaften (291) (remove)
The problem of estimating the maximum possible earthquake magnitude m(max) has attracted growing attention in recent years. Due to sparse data, the role of uncertainties becomes crucial. In this work, we determine the uncertainties related to the maximum magnitude in terms of confidence intervals. Using an earthquake catalog of Iran, m(max) is estimated for different predefined levels of confidence in six seismotectonic zones. Assuming the doubly truncated Gutenberg-Richter distribution as a statistical model for earthquake magnitudes, confidence intervals for the maximum possible magnitude of earthquakes are calculated in each zone. While the lower limit of the confidence interval is the magnitude of the maximum observed event, the upper limit is calculated from the catalog and the statistical model. For this aim, we use the original catalog which no declustering methods applied on as well as a declustered version of the catalog. Based on the study by Holschneider et al. (Bull Seismol Soc Am 101(4): 1649-1659, 2011), the confidence interval for m(max) is frequently unbounded, especially if high levels of confidence are required. In this case, no information is gained from the data. Therefore, we elaborate for which settings finite confidence levels are obtained. In this work, Iran is divided into six seismotectonic zones, namely Alborz, Azerbaijan, Zagros, Makran, Kopet Dagh, Central Iran. Although calculations of the confidence interval in Central Iran and Zagros seismotectonic zones are relatively acceptable for meaningful levels of confidence, results in Kopet Dagh, Alborz, Azerbaijan and Makran are not that much promising. The results indicate that estimating mmax from an earthquake catalog for reasonable levels of confidence alone is almost impossible.
Flood loss modeling is an important component for risk analyses and decision support in flood risk management. Commonly, flood loss models describe complex damaging processes by simple, deterministic approaches like depth-damage functions and are associated with large uncertainty. To improve flood loss estimation and to provide quantitative information about the uncertainty associated with loss modeling, a probabilistic, multivariable Bagging decision Tree Flood Loss Estimation MOdel (BT-FLEMO) for residential buildings was developed. The application of BT-FLEMO provides a probability distribution of estimated losses to residential buildings per municipality. BT-FLEMO was applied and validated at the mesoscale in 19 municipalities that were affected during the 2002 flood by the River Mulde in Saxony, Germany. Validation was undertaken on the one hand via a comparison with six deterministic loss models, including both depth-damage functions and multivariable models. On the other hand, the results were compared with official loss data. BT-FLEMO outperforms deterministic, univariable, and multivariable models with regard to model accuracy, although the prediction uncertainty remains high. An important advantage of BT-FLEMO is the quantification of prediction uncertainty. The probability distribution of loss estimates by BT-FLEMO well represents the variation range of loss estimates of the other models in the case study.
In this study, we investigated the scale sizes of equatorial plasma irregularities (EPIs) using measurements from the Swarm satellites during its early mission and final constellation phases. We found that with longitudinal separation between Swarm satellites larger than 0.4°, no significant correlation was found any more. This result suggests that EPI structures include plasma density scale sizes less than 44 km in the zonal direction. During the Swarm earlier mission phase, clearly better EPI correlations are obtained in the northern hemisphere, implying more fragmented irregularities in the southern hemisphere where the ambient magnetic field is low. The previously reported inverted-C shell structure of EPIs is generally confirmed by the Swarm observations in the northern hemisphere, but with various tilt angles. From the Swarm spacecrafts with zonal separations of about 150 km, we conclude that larger zonal scale sizes of irregularities exist during the early evening hours (around 1900 LT).
The reservoir sediments are important sinks for organic carbon (OC), the OC burial being dependent on two opposite processes, deposition and mineralization. Hence factors such as severe water level fluctuations are expected to influence the rate of OC accumulation as they may affect both deposition and mineralization. The Barasona Reservoir has been historically threatened by siltation, whilst the use of water for irrigation involves a drastic decrease of the water level. In this context, we have studied the physical and chemical characteristics (grain size, major and minor elemental compositions, organic and inorganic carbon, and nitrogen) of the recent sediments of the Barasona Reservoir and the relationships among them in order to: a) elucidate the main processes governing OC accumulation, b) evaluate the rate of OC mineralization and c) approach the effect of drought on the sediment characteristics in this system. Our results indicated that Barasona sediments were dominated by fine silts (>60%) and clays (>20%), the mean particle size decreasing from tail to dam. Desiccation increased particle sorting and size distribution became bimodal, but no effect on average size was observed. Attending to the composition, Barasona sediments were very homogeneous with low concentrations of nitrogen (TN) and phosphorus (<12 g kg(-1) dw and <0.6 g kg(-1) dw, respectively) and high concentration of OC (approximate to 36 kg(-1) dw). TN was negatively related to dry weight Sediment mixing due to drastic changes in water level may have favoured the observed homogeneity of Barasona sediments affecting carbon, major ions and grain size. The high amount of OC deposited in Barasona sediment suggested that the adsorption of OC onto fine particles was more important than in boreal lakes. The rate of oxygen consumption by wet sediment ranged from 2.26 to 3.15 mg O-2 m(-2) day(-1), values close to those compiled for Mediterranean running waters. (C) 2015 Elsevier B.V. All rights reserved.
Regulation alters the characteristics of riversty transforming parts of them into lakes, affecting their hydrology and also the physical, chemical, and biological characteristics and dynamics. Reservoirs have proven to be very effective retaining particulate materials, thereby avoiding the downstream transport of suspended sediment and the chemical substances associated with it (e.g. Carbon, C or Nitrogen, N). The study of fluvial transport of C and N is of great interest since river load represents a major link to the global C and N cycles. Moreover, reservoirs are the most important sinks for organic carbon among inland waters and have a potential significance as nitrogen sinks. In this respect, this paper investigates the effects of a Pyrenean reservoir on the runoff, suspended sediment, C and N derived from the highly active Esera and Isabena rivets. Key findings indicate that the reservoir causes a considerable impact on the Esera-Isabena river fluxes, reducing them dramatically as almost all the inputs are retained within the reservoir. Despite the very dry study year (2011-2012), it can be calculated that almost 300,000 t of suspended sediment were deposited into the Barasona Reservoir, from which more than 16,000 were C (i.e. 2200 t as organic C) and 222 t were N. These values may not be seen as remarkable in a wider global context but, assuming that around 30 hm(3) of sediment are currently stored in the reservoir, figures would increase up to ca. 2.6 x 10(6) t of C (i.e. 360,000 t of organic C) and 35,000 t of N. Nevertheless, these values are indicative and should be treated with caution as there is incomplete understanding of all the processes which affect C and N. Further investigation to establish a more complete picture of C and N yields and budgets by monitoring the different processes involved is essential. (C) 2015 Elsevier B.V. All rights reserved.
Mediterranean climate is characterized by highly irregular rainfall patterns with marked differences between wet and dry seasons which lead to highly variable hydrological fluvial regimes. As a result, and in order to ensure water availability and reduce its temporal variability, a high number of large dams were built during the 20th century (more than 3500 located in Mediterranean rivers). Dams modify the flow regime but also interrupt the continuity of sediment transfer along the river network, thereby changing its functioning as an ecosystem. Within this context, the present paper aims to assess the suspended sediment loads and dynamics of two climatically contrasting Mediterranean regulated rivers (i.e. the Esera and Siurana) during a 2-yr period. Key findings indicate that floods were responsible for 92% of the total suspended sediment load in the River Siurana, while this percentage falls to 70% for the Esera, indicating the importance of baseflows on sediment transport in this river. This fact is related to the high sediment availability, with the Esera acting as a non-supply-limited catchment due to the high productivity of the sources (i.e. badlands). In contrast, the Siurana can be considered a supply-limited system due to its low geomorphic activity and reduced sediment availability, with suspended sediment concentration remaining low even for high magnitude flood events. Reservoirs in both rivers reduce sediment load up to 90%, although total runoff is only reduced in the case of the River Esera. A remarkable fact is the change of the hydrological character of the River Lem downstream for the dam, shifting from a humid mountainous river regime to a quasi-invariable pattern, whereas the Siurana experiences the opposite effect, changing from a flashy Mediterranean river to a more constant flow regime below the dam. (C) 2015 Elsevier B.V. All rights reserved.
The uppermost Cretaceous (upper Campanian-Maastrichtian) marine deposits of the central south Pyrenees host a rich larger benthic foraminiferal fauna and several rudist-rich levels. These marine deposits are directly overlain by the continental facies of the Aren and Tremp Formations, which are famous for their fossil dinosaur remains. Larger benthic foraminiferal distribution documents an important faunal turnover in all the carbonate platform environments within the photic zone, from open marine to littoral areas. Biostratigraphy indicates that this turnover occurred close to the Campanian-Maastrichtian boundary. This is also confirmed by strontium isotope stratigraphy which indicates an earliest Maastrichtian age for the appearance of the larger benthic foraminiferal assemblage constituted by Lepidorbitoides socialis, Clypeorbis mammillata, Wannierina cataluniensis, Orbitoides gruenbachensis, Siderolites aff. calcitrapoides, Fascispira colomi, Omphalocyclus macroporus and Laffiteina mengaudi. In particular, a numerical age of 71 Ma is obtained for the Hippurites radiosus level, just a few meters below the first continental deposits of the Aren sensu stricto Formation. The youngest marine sediments of the central south Pyrenees are early Maastrichtian in age. This is also an important constraint for the age of the end-Cretaceous dinosaur fossil localities of the Tremp basin. (C) 2015 Elsevier Ltd. All rights reserved.
Sarmentofascis zamparelliae n. sp., a new demosponge from the lower Campanian of southern Italy
(2016)
A new coralline sponge, exhibiting typical "stromatoporoid" bodyplan, is described as Sarmentofascis zamparelliae n. sp. from the lower Campanian of the southern Apennines, Italy. It is differentiated from Sarmentofascis cretacea (Turnsek) (Hauterivian of Montenegro) and Sarmentofascis chabrieri Termier, Termier and Vachard (Santonian of France) above all by its slender arborescent skeleton, exhibiting longitudinally distributed astrorhizae-like canals. S. zamparelliae n. sp. is the youngest representative of the genus and is reported from a period exhibiting a distinct decline of "stromatoporoid" sponges. With its clinogonal microstructure and occurrence in inner platform stromatoporoid-foraminiferan floatstones it can be considered a Late Cretaceous environmental analog to the Late Jurassic Cladocoropsis. (C) 2015 Elsevier Ltd. All rights reserved.
The water cycle of sites with shallow groundwater tables is characterized by complex interactions of hydrological and ecological processes. The water balance components, which are subject to diurnal fluctuations, are best measured with groundwater lysimeters. However, the lower boundary condition of such lysimeters affects most of the hydrological variables, particularly when considering short time scales, and has to be defined in such a way as to facilitate realistic simulations. In this paper, different means of controlling the lower boundary condition of groundwater lysimeters were compared with respect to their ability to simulate the behavior of the water balance components properly. Measurements of rain-free periods from a lysimeter station installed in the Spreewald wetland in north-east Germany were evaluated. The most common groundwater lysimeter type is controlled using a Mariotte bottle and sets the groundwater level in the soil monolith to a constant level, which here caused an alteration of the inflow to the lysimeter, with respect to both its value and diurnal behavior. Still, daily evapotranspiration values were realistic and this simple and robust approach may be used for time intervals not shorter than one day. High-resolution measurements can be gained from lysimeters that automatically adjust the groundwater level by a system of pumps and valves on an hourly basis. Still, reliable results were only obtained when the conditions in the lysimeter and the surrounding field, where the target groundwater level was measured, were in accordance. Otherwise (e.g., when the groundwater level differed) an unrealistic inflow behavior evolved. Reasonable results, even for slightly diverging conditions, were gained with a new approach that defined the lower boundary conditions by controlling the inflows and outflows of the lysimeter. This approach further enabled the groundwater level itself to be the study subject, thereby enlarging the field of possible applications of groundwater lysimeters. (C) 2015 Elsevier B.V. All rights reserved.
Fossil oyster shells are well-suited to provide palaeotemperature proxies from geologic to seasonal timescales due to their ubiquitous occurrence from Triassic to Quaternary sediments, the seasonal nature of their shell growth and their relative strong resistance to post-mortem alteration. However, the common use to translate calcitic oxygen isotopes into palaeotemperatures is challenged by uncertainties in accounting for past seawater delta O-18, especially in shallow coastal environment where oysters calcify. In principle, the Mg/Ca ratio in oyster shells can provide an alternative palaeothermometer. Several studies provided temperature calibrations for this potential proxy based on modem species, nevertheless their application to palaeo-studies remains hitherto unexplored. Here, we show that past temperature variability in seawater can be obtained from Mg/Ca analyses from selected fossil oyster species and specimens. High-resolution Mg/Ca profiles, combined with delta O-18, were obtained along 41 fossil oyster shells of seven different species from the Palaeogene Proto-Paratethys sea (Central Asia) found in similar as well as different depositional age and environments providing comparison. Suitable Mg/Ca profiles, defined by continuous cyclicity and reproducibility within one shell, are found to be consistent for specimens of the same species but differ systematically between species, implying a dominant species-specific effect on the Mg/Ca signal. Two species studied here (Ostrea (Turkostrea) strictiplicata and Sokolowia buhsii) provide an excellent proxy for palaeoclimate reconstruction from China to Europe in Palaeogene marine sediments. More generally, the protocol developed here can be applied to identify other fossil oyster species suitable for palaeoclimate reconstructions. (C) 2015 Elsevier B.V. All rights reserved.
A link between chemical weathering and physical erosion exists at the catchment scale over a wide range of erosion rates(1,2). However, in mountain environments, where erosion rates are highest, weathering may be kinetically limited(3-5) and therefore decoupled from erosion. In active mountain belts, erosion is driven by bedrock landsliding(6) at rates that depend strongly on the occurrence of extreme rainfall or seismicity(7). Although landslides affect only a small proportion of the landscape, bedrock landsliding can promote the collection and slow percolation of surface runoff in highly fragmented rock debris and create favourable conditions for weathering. Here we show from analysis of surface water chemistry in the Southern Alps of New Zealand that weathering in bedrock landslides controls the variability in solute load of these mountain rivers. We find that systematic patterns in surface water chemistry are strongly associated with landslide occurrence at scales from a single hillslope to an entire mountain belt, and that landslides boost weathering rates and river solute loads over decades. We conclude that landslides couple erosion and weathering in fast-eroding uplands and, thus, mountain weathering is a stochastic process that is sensitive to climatic and tectonic controls on mass wasting processes.
The Puna Plateau, adjacent Eastern Cordillera and the Sierras Pampeanas of the central Andes are largely characterized by thick-skinned, basement-involved deformation. The Puna Plateau hosts similar to N-S trending bedrock ranges bounded by deep-seated reverse faults and sedimentary basins. We contribute to the understanding of thick-skinned dynamics in the Puna Plateau by constraining regional kinematics of the poorly understood southern Puna Plateau through a multidisciplinary approach. On the southeastern plateau, sandstone modal composition and detrital zircon U-Pb and apatite fission-track data from Cenozoic strata indicate basin accumulation during the late Eocene to early Oligocene (similar to 38-28 Ma). Provenance analysis reveals the existence of a regional-scale basin covering the southern Puna Plateau during late Eocene to early Oligocene time (similar to 38-28 Ma) that was sourced from both the western plateau and the eastern plateau margin and had a depocenter located to the west. Petrographic and detrital zircon U-Pb data reveal erosion of proximal western and eastern sources after 12 Ma, in mid-late Miocene time. This indicates that the regional basin was compartmentalized into small-scale depocenters by the growth of basement-cored ranges continuing into the late Miocene (similar to 12-8 Ma). We suggest that the Cenozoic history of the southern Puna Plateau records the formation of a regional basin that was possibly driven by lithospheric flexure during the late Eocene to early Oligocene, before the growth of distributed basement-cored ranges starting as early as the late Oligocene. (C) 2015 Elsevier B.V. All rights reserved.
New age data have been obtained to time constrain the recent Quaternary volcanism of El Hierro (Canary Islands) and to estimate its recurrence rate. We have carried out Ar-40/Ar-39 geochronology on samples spanning the entire volcanostratigraphic sequence of the island and C-14 geochronology on the most recent eruption on the northeast rift of the island: 2280 +/- 30 yr BP. We combine the new absolute data with a revision of published ages onshore, some of which were identified through geomorphological criteria (relative data). We present a revised and updated chronology of volcanism for the last 33 ka that we use to estimate the maximum eruptive recurrence of the island. The number of events per year determined is 9.7 x 10(-4) for the emerged part of the island, which means that, as a minimum, one eruption has occurred approximately every 1000 years. This highlights the need of more geochronological data to better constrain the eruptive recurrence of El Hierro. (C) 2015 Elsevier Ltd. All rights reserved.
The Alanya Massif, which is located to the south of central Taurides in Turkey, presents a typical nappe pile consisting of thrust sheets with contrasting metamorphic histories. In two thrust sheets, Sugozu and GundogmuAY nappes, HP metamorphism under eclogite (550-567 A degrees C/14-18 kbar) and blueschist facies (435-480 A degrees C/11-13 kbar) conditions have been recognized, respectively. Whereas the rest of the Massif underwent MP metamorphism under greenschist to amphibolite facies (525-555 A degrees C/6.5-7.5 kbar) conditions. Eclogite facies metamorphism in Sugozu nappe, which consists of homogeneous garnet-glaucophane-phengite schists with eclogite lenses is dated at 84.8 +/- A 0.8, 84.7 +/- A 1.5 and 82 +/- A 3 Ma (Santonian-Campanian) by Ar-40/Ar-39 phengite, U/Pb zircon and rutile dating methods, respectively. Similarly, phengites in GundogmuAY nappe representing an accretionary complex yield 82-80 Ma (Campanian) ages for blueschist facies metamorphism. During the exhumation, the retrograde overprint of the HP units under greenschist-amphibolite facies conditions and tectonic juxtaposition with the Barrovian units occurred during Campanian (75-78 Ma). Petrological and geochronological data clearly indicate a similar Late Cretaceous tectonometamorphic evolution for both Alanya (84-75 Ma) and Bitlis (84-72 Ma) Massifs. They form part of a single continental sliver (Alanya-Bitlis microcontinent), which was rifted from the southern part of the Anatolide-Tauride platform. The P-T-t coherence between two Massifs suggests that both Massifs have been derived from the closure of the same ocean (Alanya-Bitlis Ocean) located to the south of the Anatolide-Tauride block by a northward subduction. The boundary separating the autochthonous Tauride platform to the north from both the Alanya and Bitlis Massifs to the south represents a suture zone, the Pamphylian-Alanya-Bitlis suture.
Anatolia’s high-pressure metamorphic belts are characterized in part by a Neotethyan stratigraphic succession that includes a mid-Cretaceous hemi-pelagic marble sequence. This unit contains, towards its stratigraphic top, dm-to-m-long radiating calcitic rods forming rosette-like textures. Here, we refer to these features as “Rosetta Marble”. The remarkable textural similarity of non-metamorphic selenite crystals and radiating calcite rods in the Rosetta Marble strongly suggests that these textures represent pseudomorphs after selenites. Metamorphosed hemi-pelagic limestones, dominated by Rosetta selenite pseudomorphs, are alternating with siliceous meta-sediments containing relictic radiolaria tests. This stratigraphic pattern is indicative of transient phases characterized by evaporites precipitated from basinal brines alternating with non-evaporative hemi-pelagic deposition from normal-marine seawater. The regional distribution of Rosetta Marble exposures over 600 km is indicative of basin-scale evaporitic intervals. High-pressure, low-temperature metamorphism of these rocks is witnessed by Sr-rich (up to 3500 ppm), fibrous calcite pseudomorphs after aragonite and isolated aragonite inclusions in quartz. Peak metamorphic conditions of 1.2 GPa and 300–350 °C are attested by high-Si white mica thermobarometry. The Rosetta Marble case example examines the potential to unravel the complete history from deposition to diagenesis and metamorphism of meta-sedimentary rocks.
This paper addresses the lithosphere-scale subduction-collision history of the eastern termination of the Aegean retreating subduction system, i.e. western Anatolia. Although there is some general consensus on the protracted subduction evolution of the Aegean since the early Cenozoic at least, correlation with western Anatolia has been widely debated for more than several decades. In western Anatolia, three main tectonic configurations have been envisaged in the past years to reconstruct slab dynamics during the closure of the Neotethyan oceanic realm since the Late Cretaceous. Some authors have suggested an Aegean-type scenario, with the continuous subduction of a single lithospheric slab, punctuated by episodic slab roll-back and trench retreat, whereas others assumed a discontinuous subduction history marked by intermittent slab break-off during either the Campanian (ca. 75 Ma) or the Early Eocene (ca. 55-50 Ma). The third view implies three partly contemporaneous subduction zones. Our review of these models points to key debated aspects that can be re-evaluated in the light of multidisciplinary constraints from the literature. Our discussion leads us to address the timing of subduction initiation, the existence of hypothetical ocean basins, the number of intervening subduction zones between the Taurides and the Pontides, the palaeogeographic origin of tectonic units and the possibility for slab break-off during either the Campanian or the Early Eocene. Thence, we put forward a favoured tectonic scenario featuring two successive phases of subduction of a single lithospheric slab and episodic accretion of two continental domains separated by a continental trough, representing the eastern end of the Cycladic Ocean of the Aegean. The lack of univocal evidence for slab break-off in western Anatolia and southward-younging HP/LT metamorphism in continental tectonic units (from similar to 85, 70 to 50 Ma) in the Late Cretaceous-Palaeogene period suggests continuous subduction since similar to 110 Ma, marked by roll-back episodes in the Palaeocene and the Oligo-Miocene, and slab tearing below western Anatolia during the Miocene.
A unique assemblage including kumdykolite and kokchetavite, polymorphs of albite and K-feldspar, respectively, together with cristobalite, micas, and calcite has been identified in high-pressure granulites of the Orlica-Snieznik dome (Bohemian Massif) as the product of partial melt crystallization in preserved nanogranites. Previous reports of both kumdykolite and kokchetavite in natural rocks are mainly from samples that passed through the diamond stability field. However, because the maximum pressure recorded in these host rocks is <3 GPa, our observations indicate that high pressure is not required for the formation of kumdykolite and kokchetavite, and their presence is not therefore an indicator of ultrahigh-pressure conditions. Detailed microstructural and microchemical investigation of these inclusions indicates that such phases should instead be regarded as (1) a direct mineralogical criteria to identify former melt inclusions with preserved original compositions, including H2O and CO2 contents and (2) indicators of rapid cooling of the host rocks. Thus, the present study provides novel criteria for the interpretation of melt inclusions in natural rocks and allows a more rigorous characterization of partial melts during deep subduction to mantle depth as well as their behavior on exhumation.
The South American Andes are frequently exposed to intense rainfall events with varying moisture sources and precipitation-forming processes. In this study, we assess the spatiotemporal characteristics and geographical origins of rainfall over the South American continent. Using high-spatiotemporal resolution satellite data (TRMM 3B42 V7), we define four different types of rainfall events based on their (1) high magnitude, (2) long temporal extent, (3) large spatial extent, and (4) high magnitude, long temporal and large spatial extent combined. In a first step, we analyze the spatiotemporal characteristics of these events over the entire South American continent and integrate their impact for the main Andean hydrologic catchments. Our results indicate that events of type 1 make the overall highest contributions to total seasonal rainfall (up to 50%). However, each consecutive episode of the infrequent events of type 4 still accounts for up to 20% of total seasonal rainfall in the subtropical Argentinean plains. In a second step, we employ complex network theory to unravel possibly non-linear and long-ranged climatic linkages for these four event types on the high-elevation Altiplano-Puna Plateau as well as in the main river catchments along the foothills of the Andes. Our results suggest that one to two particularly large squall lines per season, originating from northern Brazil, indirectly trigger large, long-lasting thunderstorms on the Altiplano Plateau. In general, we observe that extreme rainfall in the catchments north of approximately 20 degrees S typically originates from the Amazon Basin, while extreme rainfall at the eastern Andean foothills south of 20 degrees S and the Puna Plateau originates from southeastern South America.
Freshwater ostracods (Crustacea, Ostracoda) are valuable biological indicators. In Arctic environments, their habitat conditions are barely known and the abundance and diversity of ostracods is documented only in scattered records with incomplete ecological characterization. To determine the taxonomic range of ostracod assemblages and their habitat conditions in polygon ponds in the Indigirka Lowland, north-east Siberia, we collected more than 100 living ostracod individuals per site with a plankton net (mesh size 65 mm) and an exhaustor system from 27 water bodies and studied them in the context of substrate and hydrochemical data. During the summer of 2011, a single pond site and its ostracod population was selected for special study. This first record of the ostracod fauna in the Indigirka Lowland comprises eight species and three additional taxa. Fabaeformiscandona krochini and F. groenlandica were documented for the first time in continental Siberia. Repeated sampling of a low-centre polygon pond yielded insights into the population dynamics of F. pedata. We identified air temperature and precipitation as the main external drivers of water temperatures, water levels, ion concentrations and water stable isotope composition on diurnal and seasonal scales.
Permafrost-related processes drive regional landscape dynamics in the Arctic terrestrial system. A better understanding of past periods indicative of permafrost degradation and aggradation is important for predicting the future response of Arctic landscapes to climate change. Here, we used a multi-proxy approach to analyse a4m long sediment core from a drained thermokarst lake basin on the northern Seward Peninsula in western Arctic Alaska (USA). Sedimentological, biogeochemical, geochronological, micropalaeontological (ostracoda, testate amoebae) and tephra analyses were used to determine the long-term environmental Early-Wisconsin to Holocene history preserved in our core for central Beringia. Yedoma accumulation dominated throughout the Early to Late-Wisconsin but was interrupted by wetland formation from 44.5 to 41.5ka BP. The latter was terminated by the deposition of 1m of volcanic tephra, most likely originating from the South Killeak Maar eruption at about 42ka BP. Yedoma deposition continued until 22.5ka BP and was followed by a depositional hiatus in the sediment core between 22.5 and 0.23ka BP. We interpret this hiatus as due to intense thermokarst activity in the areas surrounding the site, which served as a sediment source during the Late-Wisconsin to Holocene climate transition. The lake forming the modern basin on the upland initiated around 0.23ka BP and drained catastrophically in spring 2005. The present study emphasises that Arctic lake systems and periglacial landscapes are highly dynamic and that permafrost formation as well as degradation in central Beringia was controlled by regional to global climate patterns as well as by local disturbances. Copyright (c) 2015 John Wiley & Sons, Ltd.
The new sediment record from the deep Dead Sea basin (ICDP core 5017-1) provides a unique archive for hydroclimatic variability in the Levant. Here, we present high-resolution sediment facies analysis and elemental composition by micro-X-ray fluorescence (mu XRF) scanning of core 5017-1 to trace lake levels and responses of the regional hydroclimatology during the time interval from ca. 117 to 75 ka, i. e. the transition between the last interglacial and the onset of the last glaciation. We distinguished six major micro-facies types and interpreted these and their alterations in the core in terms of relative lake level changes. The two end-member facies for highest and lowest lake levels are (a) up to several metres thick, greenish sediments of alternating aragonite and detrital marl laminae (aad) and (b) thick halite facies, respectively. Intermediate lake levels are characterised by detrital marls with varying amounts of aragonite, gypsum or halite, reflecting lower-amplitude, shorter-term variability. Two intervals of pronounced lake level drops occurred at similar to 110-108 +/- 5 and similar to 93-87 +/- 7 ka. They likely coincide with stadial conditions in the central Mediterranean (Melisey I and II pollen zones in Monticchio) and low global sea levels during Marine Isotope Stage (MIS) 5d and 5b. However, our data do not support the current hypothesis of an almost complete desiccation of the Dead Sea during the earlier of these lake level low stands based on a recovered gravel layer. Based on new petrographic analyses, we propose that, although it was a low stand, this well-sorted gravel layer may be a vestige of a thick turbidite that has been washed out during drilling rather than an in situ beach deposit. Two intervals of higher lake stands at similar to 108-93 +/- 6 and similar to 87-75 +/- 7 ka correspond to interstadial conditions in the central Mediterranean, i. e. pollen zones St. Germain I and II in Monticchio, and Greenland interstadials (GI) 24+23 and 21 in Greenland, as well as to sapropels S4 and S3 in the Mediterranean Sea. These apparent correlations suggest a close link of the climate in the Levant to North Atlantic and Mediterranean climates during the time of the build-up of Northern Hemisphere ice shields in the early last glacial period.
Pollen influx (number of pollen grains cm−2 year−1) can objectively reflect the dispersal and deposition features of pollen within a certain time and space, and is often used as a basis for the quantitative reconstruction of palaeovegetation; however, little is known about the features and mechanisms of vertical dispersal of pollen. Here we present the results from a 5 year (2006–2010) monitoring program using pollen traps placed at different heights from ground level up to 60 m and surface soil samples in a mixed coniferous and deciduous broad-leaved woodland in the Changbai mountains, northeastern China. The pollen percentages and pollen influx from the traps have very similar characteristics to the highest values for Betula, Fraxinus, Quercus and Pinus, among the tree taxa and Artemisia, Chenopodiaceae and Asteraceae among the herb taxa. Pollen influx values vary significantly with height and show major differences between three distinct layers, above-canopy (≥32 m), within the trunk layer (8 ≤ 32 m) and on the ground (0 m). These differences in pollen influx are explained by differences in (i) the air flows in each of these layers and (ii) the fall speed of pollen of the various taxa. We found that the pollen recorded on the ground surface is a good representation of the major part of the pollen transported in the trunk space of the woodland. Comparison of the pollen influx values with the theoretical, calculated “characteristic pollen source area” (CPSA) of 12 selected taxa indicates that the pollen deposited on the ground surface of the woodland is a fair representation with 85–90 % of the total pollen deposited at a wind speed of 2.4 m s−1 coming from within ca. 1–5 km for Pinus and Quercus, ca. 5–10 km for Ulmus, Tilia, Oleaceae and Betula, ca. 20–40 km for Fraxinus, Poaceae, Chenopodiaceae, Populus and Salix, and ca. 30–60 km for Artemisia; it is also a good representation with 90–98 % of the total pollen deposited coming from within 60 km at a wind speed of 2.4 m s−1, or 100 km at a wind speed: 6 m s−1, for the 12 selected taxa used in the CPSA calculation. Furthermore, comparison with the vegetation map of the area around the sampling site shows that the pollen deposited on the ground represents all plant communities which grow in the study area within 70 km radius of the sampling site. In this study, the pollen percentages obtained from the soil surface samples are significantly biased towards pollen taxa with good preservation due to thick and robust pollen walls. Therefore, if mosses are available instead, soil samples should be avoided for pollen studies, in particular for the study of pollen-vegetation relationships, the estimation of pollen productivities and quantitative reconstruction of past vegetation. The results also indicate that the existing model of pollen dispersal and deposition, Prentice’s model, provides a fair description of the actual pollen dispersal and deposition in this kind of woodland, which suggests that the application of the landscape reconstruction algorithm would be relevant for reconstruction of this type of woodland in the past.
A great number of Central Asian wall paintings, archeological materials, architectural fragments, and textiles, as well as painting fragments on silk and paper, make up the so called Turfan Collection at the Asian Art Museum in Berlin. The largest part of the collection comes from the Kucha region, a very important cultural center in the third to ninth centuries. Between 1902 and 1914, four German expeditions traveled along the northern Silk Road. During these expeditions, wall paintings were detached from their original settings in Buddhist cave complexes. This paper reports a technical study of a wall painting, existing in eight fragments, from the Buddhist cave no. 40 (Ritterhohle). Its original painted surface is soot blackened and largely illegible. Gruwedel, leader of the first and third expeditions, described the almost complete destruction of the rediscovered temple complex and evidence of fire damage. The aim of this case study is to identify the materials used for the wall paintings. Furthermore, soot deposits as well as materials from conservation interventions were of interest. Non-invasive analyses were preferred but a limited number of samples were taken to provide more precise information on the painting technique. By employing optical and scanning electron microscopy, energy dispersive X-ray spectroscopy, micro X-ray fluorescence spectroscopy, X-ray diffraction analysis, and Raman spectroscopy, a layer sequence of earthen render, a ground layer made of gypsum, and a paint layer containing a variety of inorganic pigments were identified.
In contrast to recent advances in projecting sea levels, estimations about the economic impact of sea level rise are vague. Nonetheless, they are of great importance for policy making with regard to adaptation and greenhouse-gas mitigation. Since the damage is mainly caused by extreme events, we propose a stochastic framework to estimate the monetary losses from coastal floods in a confined region. For this purpose, we follow a Peak-over-Threshold approach employing a Poisson point process and the Generalised Pareto Distribution. By considering the effect of sea level rise as well as potential adaptation scenarios on the involved parameters, we are able to study the development of the annual damage. An application to the city of Copenhagen shows that a doubling of losses can be expected from a mean sea level increase of only 11 cm. In general, we find that for varying parameters the expected losses can be well approximated by one of three analytical expressions depending on the extreme value parameters. These findings reveal the complex interplay of the involved parameters and allow conclusions of fundamental relevance. For instance, we show that the damage typically increases faster than the sea level rise itself. This in turn can be of great importance for the assessment of sea level rise impacts on the global scale. Our results are accompanied by an assessment of uncertainty, which reflects the stochastic nature of extreme events. While the absolute value of uncertainty about the flood damage increases with rising mean sea levels, we find that it decreases in relation to the expected damage.
This case study evaluates the suitability of radar-based quantitative precipitation estimates (QPEs) for the simulation of streamflow in the Marikina River Basin (MRB), the Philippines. Hourly radar-based QPEs were produced from reflectivity that had been observed by an S-band radar located about 90 km from the MRB. Radar data processing and precipitation estimation were carried out using the open source library wradlib. To assess the added value of the radar-based QPE, we used spatially interpolated rain gauge observations (gauge-only (GO) product) as a benchmark. Rain gauge observations were also used to quantify rainfall estimation errors at the point scale. At the point scale, the radar-based QPE outperformed the GO product in 2012, while for 2013, the performance was similar. For both periods, estimation errors substantially increased from daily to the hourly accumulation intervals. Despite this fact, both rainfall estimation methods allowed for a good representation of observed streamflow when used to force a hydrological simulation model of the MRB. Furthermore, the results of the hydrological simulation were consistent with rainfall verification at the point scale: the radar-based QPE performed better than the GO product in 2012, and equivalently in 2013. Altogether, we could demonstrate that, in terms of streamflow simulation, the radar-based QPE can perform as good as or even better than the GO product - even for a basin such as the MRB which has a comparatively dense rain gauge network. This suggests good prospects for using radar-based QPE to simulate and forecast streamflow in other parts of the Philippines where rain gauge networks are not as dense.
Isostasy is one of the oldest and most widely applied concepts in the geosciences, but the geoscientific community lacks a coherent, easy-to-use tool to simulate flexure of a realistic (i.e., laterally heterogeneous) lithosphere under an arbitrary set of surface loads. Such a model is needed for studies of mountain building, sedimentary basin formation, glaciation, sea-level change, and other tectonic, geodynamic, and surface processes. Here I present gFlex (for GNU flexure), an open-source model that can produce analytical and finite difference solutions for lithospheric flexure in one (profile) and two (map view) dimensions. To simulate the flexural isostatic response to an imposed load, it can be used by itself or within GRASS GIS for better integration with field data. gFlex is also a component with the Community Surface Dynamics Modeling System (CSDMS) and Landlab modeling frameworks for coupling with a wide range of Earth-surface-related models, and can be coupled to additional models within Python scripts. As an example of this in-script coupling, I simulate the effects of spatially variable lithospheric thickness on a modeled Iceland ice cap. Finite difference solutions in gFlex can use any of five types of boundary conditions: 0-displacement, 0-slope (i.e., clamped); 0-slope, 0-shear; 0-moment, 0-shear (i.e., broken plate); mirror symmetry; and periodic. Typical calculations with gFlex require << 1 s to similar to 1 min on a personal laptop computer. These characteristics - multiple ways to run the model, multiple solution methods, multiple boundary conditions, and short compute time - make gFlex an effective tool for flexural isostatic modeling across the geosciences.
Injection of fluids into deep saline aquifers causes a pore pressure increase in the storage formation, and thus displacement of resident brine. Via hydraulically conductive faults, brine may migrate upwards into shallower aquifers and lead to unwanted salinisation of potable groundwater resources. In the present study, we investigated different scenarios for a potential storage site in the Northeast German Basin using a three-dimensional (3-D) regional-scale model that includes four major fault zones. The focus was on assessing the impact of fault length and the effect of a secondary reservoir above the storage formation, as well as model boundary conditions and initial salinity distribution on the potential salinisation of shallow groundwater resources. We employed numerical simulations of brine injection as a representative fluid. Our simulation results demonstrate that the lateral model boundary settings and the effective fault damage zone volume have the greatest influence on pressure build-up and development within the reservoir, and thus intensity and duration of fluid flow through the faults. Higher vertical pressure gradients for short fault segments or a small effective fault damage zone volume result in the highest salinisation potential due to a larger vertical fault height affected by fluid displacement. Consequently, it has a strong impact on the degree of shallow aquifer salinisation, whether a gradient in salinity exists or the saltwater-freshwater interface lies below the fluid displacement depth in the faults. A small effective fault damage zone volume or low fault permeability further extend the duration of fluid flow, which can persist for several tens to hundreds of years, if the reservoir is laterally confined. Laterally open reservoir boundaries, large effective fault damage zone volumes and intermediate reservoirs significantly reduce vertical brine migration and the potential of freshwater salinisation because the origin depth of displaced brine is located only a few decametres below the shallow aquifer in maximum. The present study demonstrates that the existence of hydraulically conductive faults is not necessarily an exclusion criterion for potential injection sites, because salinisation of shallower aquifers strongly depends on initial salinity distribution, location of hydraulically conductive faults and their effective damage zone volumes as well as geological boundary conditions.
Ice complex deposits are characteristic, ice-rich formations in northern East Siberia and represent an important part in the arctic carbon pool. Recently, these late Quaternary deposits are the objective of numerous investigations typically relying on outcrop and borehole data. Many of these studies can benefit from a 3D structural model of the subsurface for upscaling their observations or for constraining estimations of inventories, such as the local carbon stock. We have addressed this problem of structural imaging by 3D ground-penetrating radar (GPR), which, in permafrost studies, has been primarily used for 2D profiling. We have used a 3D kinematic GPR surveying strategy at a field site located in the New Siberian Archipelago on top of an ice complex. After applying a 3D GPR processing sequence, we were able to trace two horizons at depths below 20 m. Taking available borehole and outcrop data into account, we have interpreted these two features as interfaces of major lithologic units and derived a 3D cryostratigraphic model of the subsurface. Our data example demonstrated that a 3D surveying and processing strategy was crucial at our field site and showed the potential of 3D GPR to image geologic structures in complex ice-rich permafrost landscapes.
Analysis of time-lapse ground-penetrating radar (GPR) data can provide information regarding subsurface hydrological processes, such as preferential flow. However, the analysis of time-lapse data is often limited by data quality; for example, for noisy input data, the interpretation of difference images is often difficult. Motivated by modern image-processing tools, we have developed two robust GPR attributes, which allow us to distinguish amplitude (contrast similarity) and time-shift (structural similarity) variations related to differences between individual time-lapse GPR data sets. We tested and evaluated our attributes using synthetic data of different complexity. Afterward, we applied them to a field data example, in which subsurface flow was induced by an artificial rainfall event. For all examples, we identified our structural similarity attribute to be a robust measure for highlighting time-lapse changes also in data with low signal-to-noise ratios. We determined that our new attribute-based workflow is a promising tool to analyze time-lapse GPR data, especially for imaging subsurface hydrological processes.
Most climate change impacts manifest in the form of natural hazards. Damage assessment typically relies on damage functions that translate the magnitude of extreme events to a quantifiable damage. In practice, the availability of damage functions is limited due to a lack of data sources and a lack of understanding of damage processes. The study of the characteristics of damage functions for different hazards could strengthen the theoretical foundation of damage functions and support their development and validation. Accordingly, we investigate analogies of damage functions for coastal flooding and for wind storms and identify a unified approach. This approach has general applicability for granular portfolios and may also be applied, for example, to heat-related mortality. Moreover, the unification enables the transfer of methodology between hazards and a consistent treatment of uncertainty. This is demonstrated by a sensitivity analysis on the basis of two simple case studies (for coastal flood and storm damage). The analysis reveals the relevance of the various uncertainty sources at varying hazard magnitude and on both the microscale and the macroscale level. Main findings are the dominance of uncertainty from the hazard magnitude and the persistent behaviour of intrinsic uncertainties on both scale levels. Our results shed light on the general role of uncertainties and provide useful insight for the application of the unified approach.
Climate change increases riverine carbon outgassing, while export to the ocean remains uncertain
(2016)
Any regular interaction of land and river during flooding affects carbon pools within the terrestrial system, riverine carbon and carbon exported from the system. In the Amazon basin carbon fluxes are considerably influenced by annual flooding, during which terrigenous organic material is imported to the river. The Amazon basin therefore represents an excellent example of a tightly coupled terrestrial-riverine system. The processes of generation, conversion and transport of organic carbon in such a coupled terrigenous-riverine system strongly interact and are climate-sensitive, yet their functioning is rarely considered in Earth system models and their response to climate change is still largely unknown. To quantify regional and global carbon budgets and climate change effects on carbon pools and carbon fluxes, it is important to account for the coupling between the land, the river, the ocean and the atmosphere. We developed the RIVerine Carbon Model (RivCM), which is directly coupled to the well-established dynamic vegetation and hydrology model LPJmL, in order to account for this large-scale coupling. We evaluate RivCM with observational data and show that some of the values are reproduced quite well by the model, while we see large deviations for other variables. This is mainly caused by some simplifications we assumed. Our evaluation shows that it is possible to reproduce large-scale carbon transport across a river system but that this involves large uncertainties. Acknowledging these uncertainties, we estimate the potential changes in riverine carbon by applying RivCM for climate forcing from five climate models and three CO2 emission scenarios (Special Report on Emissions Scenarios, SRES). We find that climate change causes a doubling of riverine organic carbon in the southern and western basin while reducing it by 20% in the eastern and northern parts. In contrast, the amount of riverine inorganic carbon shows a 2- to 3-fold increase in the entire basin, independent of the SRES scenario. The export of carbon to the atmosphere increases as well, with an average of about 30 %. In contrast, changes in future export of organic carbon to the Atlantic Ocean depend on the SRES scenario and are projected to either decrease by about 8.9% (SRES A1B) or increase by about 9.1% (SRES A2). Such changes in the terrigenous-riverine system could have local and regional impacts on the carbon budget of the whole Amazon basin and parts of the Atlantic Ocean. Changes in riverine carbon could lead to a shift in the riverine nutrient supply and pH, while changes in the exported carbon to the ocean lead to changes in the supply of organic material that acts as a food source in the Atlantic. On larger scales the increased outgassing of CO2 could turn the Amazon basin from a sink of carbon to a considerable source. Therefore, we propose that the coupling of terrestrial and riverine carbon budgets should be included in subsequent analysis of the future regional carbon budget.
Robust appraisals of climate impacts at different levels of global-mean temperature increase are vital to guide assessments of dangerous anthropogenic interference with the climate system. The 2015 Paris Agreement includes a two-headed temperature goal: "holding the increase in the global average temperature to well below 2 degrees C above pre-industrial levels and pursuing efforts to limit the temperature increase to 1.5 degrees C". Despite the prominence of these two temperature limits, a comprehensive overview of the differences in climate impacts at these levels is still missing. Here we provide an assessment of key impacts of climate change at warming levels of 1.5 degrees C and 2 degrees C, including extreme weather events, water availability, agricultural yields, sea-level rise and risk of coral reef loss. Our results reveal substantial differences in impacts between a 1.5 degrees C and 2 degrees C warming that are highly relevant for the assessment of dangerous anthropogenic interference with the climate system. For heat-related extremes, the additional 0.5 degrees C increase in global-mean temperature marks the difference between events at the upper limit of present-day natural variability and a new climate regime, particularly in tropical regions. Similarly, this warming difference is likely to be decisive for the future of tropical coral reefs. In a scenario with an end-of-century warming of 2 degrees C, virtually all tropical coral reefs are projected to be at risk of severe degradation due to temperature-induced bleaching from 2050 onwards. This fraction is reduced to about 90% in 2050 and projected to decline to 70% by 2100 for a 1.5 degrees C scenario. Analyses of precipitation-related impacts reveal distinct regional differences and hot-spots of change emerge. Regional reduction in median water availability for the Mediterranean is found to nearly double from 9% to 17% between 1.5 degrees C and 2 degrees C, and the projected lengthening of regional dry spells increases from 7 to 11%. Projections for agricultural yields differ between crop types as well as world regions. While some (in particular high-latitude) regions may benefit, tropical regions like West Africa, South-East Asia, as well as Central and northern South America are projected to face substantial local yield reductions, particularly for wheat and maize. Best estimate sea-level rise projections based on two illustrative scenarios indicate a 50cm rise by 2100 relative to year 2000-levels for a 2 degrees C scenario, and about 10 cm lower levels for a 1.5 degrees C scenario. In a 1.5 degrees C scenario, the rate of sea-level rise in 2100 would be reduced by about 30% compared to a 2 degrees C scenario. Our findings highlight the importance of regional differentiation to assess both future climate risks and different vulnerabilities to incremental increases in global-mean temperature. The article provides a consistent and comprehensive assessment of existing projections and a good basis for future work on refining our understanding of the difference between impacts at 1.5 degrees C and 2 degrees C warming.
Widespread flooding in June 2013 caused damage costs of (sic)6 to 8 billion in Germany, and awoke many memories of the floods in August 2002, which resulted in total damage of (sic)11.6 billion and hence was the most expensive natural hazard event in Germany up to now. The event of 2002 does, however, also mark a reorientation toward an integrated flood risk management system in Germany. Therefore, the flood of 2013 offered the opportunity to review how the measures that politics, administration, and civil society have implemented since 2002 helped to cope with the flood and what still needs to be done to achieve effective and more integrated flood risk management. The review highlights considerable improvements on many levels, in particular (1) an increased consideration of flood hazards in spatial planning and urban development, (2) comprehensive property-level mitigation and preparedness measures, (3) more effective flood warnings and improved coordination of disaster response, and (4) a more targeted maintenance of flood defense systems. In 2013, this led to more effective flood management and to a reduction of damage. Nevertheless, important aspects remain unclear and need to be clarified. This particularly holds for balanced and coordinated strategies for reducing and overcoming the impacts of flooding in large catchments, cross-border and interdisciplinary cooperation, the role of the general public in the different phases of flood risk management, as well as a transparent risk transfer system. Recurring flood events reveal that flood risk management is a continuous task. Hence, risk drivers, such as climate change, land-use changes, economic developments, or demographic change and the resultant risks must be investigated at regular intervals, and risk reduction strategies and processes must be reassessed as well as adapted and implemented in a dialogue with all stakeholders.
Rainfall-induced attenuation is a major source of underestimation for radar-based precipitation estimation at C-band. Unconstrained gate-by-gate correction procedures are known to be inherently unstable and thus not suited for unsupervised attenuation correction. In this study, we evaluate three different procedures to constrain gate-by-gate attenuation correction using reflectivity as the only input. These procedures are benchmarked against rainfall estimates from uncorrected radar data, using six years of radar observations from the single-polarized C-band radar in South-West Germany. The precipitation estimation error is obtained by comparing the radar-based estimates to rain gauge observations. All attenuation correction procedures benchmarked in this study lead to an effective improvement of precipitation estimation. The first method caps the corrections if the rain intensity increase exceeds a factor of two. The second method decreases the parameters of the attenuation correction iteratively for every radar beam calculation until attaining a stability criterion. The second method outperforms the first method and leads to a consistent distribution of path-integrated attenuation along the radar beam. As a third method, we propose a slight modification of Kraemer's approach which allows users to exert better control over attenuation correction by introducing an additional constraint that prevents unplausible corrections in cases of dramatic signal losses.
Ice-rich permafrost coasts often undergo rapid erosion, which results in land loss and release of considerable amounts of sediment, organic carbon and nutrients, impacting the near-shore ecosystems. Because of the lack of volumetric erosion data, Arctic coastal erosion studies typically report on planimetric erosion. Our aim is to explore the relationship between planimetric and volumetric coastal erosion measurements and to update the coastal erosion rates on Herschel Island in the Canadian Arctic. We used high-resolution digital elevation models to compute sediment release and compare volumetric data to planimetric estimations of coastline movements digitized from satellite imagery. Our results show that volumetric erosion is locally less variable and likely corresponds better with environmental forcing than planimetric erosion. Average sediment release volumes are in the same range as sediment release volumes calculated from coastline movements combined with cliff height. However, the differences between these estimates are significant for small coastal sections. We attribute the differences between planimetric and volumetric coastal erosion measurements to mass wasting, which is abundant along the coasts of Herschel Island. The average recorded coastline retreat on Herschel Island was 0.68m a(-1) for the period 2000-2011. Erosion rates increased by more than 50% in comparison with the period 1970-2000, which is in accordance with a recently observed increase along the Alaskan Beaufort Sea. The estimated annual sediment release was 28.2 m(3) m(-1) with resulting fluxes of 590 kg C m(-1) and 104 kg N m(-1).
Peatlands are scarce and threatened ecosystems in the semiarid region of Kyrgyzstan. Knowledge about the Kyrgyz peatlands is still poor and, especially, their genesis has hardly been investigated so far. Typically, the peatland substrates are characterised by the admixture of silt-sized particles in various quantities. In this work we report the abiotic properties and genesis of three peatlands within different altitudinal zones in southern Kyrgyzstan. We surveyed the stratification of the peatlands and their water chemistry. In addition, we investigated whether the silt found in the peatland substrates was deposited by wind, rivers or springs. The mineral constituents of the peatland substrates were analysed for particle size distribution and their elemental composition was compared with that of nearby loess, river and spring sediments using the immobile trace element titanium. One peatland shows a high abundance of different peatland substrates, indicating a frequent change of ecological conditions in the past. All three peatlands are fed by groundwater. Overgrazing and trampling by cattle has led to recent degradation of the upper peat layer. The resulting compaction of the peats prevents water from seeping into the substrates of the peatlands and subsequently changes their hydrology. Our results indicate that both wind and rivers have deposited silt in the peatlands, depending on their positions in the relief. Silts may also have been relocated by springs within the peatlands.
Steppe vegetation represents a key marker of past Asian aridification and is associated with monsoonal intensification. Little is, however, known about the origin of this pre-Oligocene vegetation, its specific composition and how it changed over time and responded to climatic variations. Here, we describe the morphological characters of Ephedraceae pollen in Eocene strata of the Xining Basin and compare the pollen composition with the palynological composition of Late Cretaceous and Paleocene deposits of the Xining Basin and the Quaternary deposits of the Qaidam Basin. We find that the Late Cretaceous steppe was dominated by Gnetaceaepollenites; in the transition from the Cretaceous to the Paleocene, Gnetaceaepollenites became extinct and Ephedripites subgenus Ephedripites dominated the flora with rare occurrences of Ephedripites subgen. Distachyapites; the middle to late Eocene presents a strong increase of Ephedripites subgen. Distachyapites; and the Quaternary/Recent is marked by a significantly lower diversity of Ephedraceae (and Nitrariaceae) compared to the Eocene. In the modern landscape of China, only a fraction of the Paleogene species diversity of Ephedraceae remains and we propose that these alterations in Ephedreaceae composition occurred in response to the climatic changes at least since the Eocene. In particular, the strong Eocene monsoons that enhanced the continental aridification may have played an important role in the evolution of Ephedripites subgen. Distachyapites triggering an evolutionary shift to wind-pollination in this group. Conceivably, the Ephedraceae/Nitrariaceae dominated steppe ended during the Eocene/Oligocene climatic cooling and aridification, which favoured other plant taxa.
Geomorphic footprints of past large Himalayan earthquakes are elusive, although they are urgently needed for gauging and predicting recovery times of seismically perturbed mountain landscapes. We present evidence of catastrophic valley infill following at least three medieval earthquakes in the Nepal Himalaya. Radiocarbon dates from peat beds, plant macrofossils, and humic silts in fine-grained tributary sediments near Pokhara, Nepal’s second-largest city, match the timing of nearby M > 8 earthquakes in ~1100, 1255, and 1344 C.E. The upstream dip of tributary valley fills and x-ray fluorescence spectrometry of their provenance rule out local sources. Instead, geomorphic and sedimentary evidence is consistent with catastrophic fluvial aggradation and debris flows that had plugged several tributaries with tens of meters of calcareous sediment from a Higher Himalayan source >60 kilometers away.
Successively smaller glacial extents have been proposed for continental Eurasia during the stadials of the last glacial period leading up to the Last Glacial Maximum (LGM). At the same time the large mountainous region east of Lake Baikal, Transbaikalia, has remained unexplored in terms of glacial chronology despite clear geomorphological evidence of substantial past glaciations. We have applied cosmogenic Be-10 exposure dating and optically stimulated luminescence to establish the first quantitative glacial chronology for this region. Based on eighteen exposure ages from five moraine complexes, we propose that large mountain ice fields existed in the Kodar and Udokan mountains during Oxygen Isotope Stage 2, commensurate with the global LGM. These ice fields fed valley glaciers (>100 km in length) reaching down to the Chara Depression between the Kodar and Udokan mountains and to the valley of the Vitim River northwest of the Kodar Mountains. Two of the investigated moraines date to the Late Glacial, but indications of incomplete exposure among some of the sampled boulders obscure the specific details of the post-LGM glacial history. In addition to the LGM ice fields in the highest mountains of Transbaikalia, we report geomorphological evidence of a much more extensive, ice-cap type glaciation at a time that is yet to be firmly resolved. (C) 2015 Elsevier Ltd. All rights reserved.
Subsurface microbial communities undertake many terminal electron-accepting processes, often simultaneously. Using a tritium-based assay, we measured the potential hydrogen oxidation catalyzed by hydrogenase enzymes in several subsurface sedimentary environments (Lake Van, Barents Sea, Equatorial Pacific, and Gulf of Mexico) with different predominant electron-acceptors. Hydrogenases constitute a diverse family of enzymes expressed by microorganisms that utilize molecular hydrogen as a metabolic substrate, product, or intermediate. The assay reveals the potential for utilizing molecular hydrogen and allows qualitative detection of microbial activity irrespective of the predominant electron-accepting process. Because the method only requires samples frozen immediately after recovery, the assay can be used for identifying microbial activity in subsurface ecosystems without the need to preserve live material. We measured potential hydrogen oxidation rates in all samples from multiple depths at several sites that collectively span a wide range of environmental conditions and biogeochemical zones. Potential activity normalized to total cell abundance ranges over five orders of magnitude and varies, dependent upon the predominant terminal electron acceptor. Lowest per-cell potential rates characterize the zone of nitrate reduction and highest per-cell potential rates occur in the methanogenic zone. Possible reasons for this relationship to predominant electron acceptor include (i) increasing importance of fermentation in successively deeper biogeochemical zones and (ii) adaptation of H(2)ases to successively higher concentrations of H-2 in successively deeper zones.
On 16 September 2015, the M-W = 8.2 Illapel megathrust earthquake ruptured the Central Chilean margin. Combining inversions of displacement measurements and seismic waveforms with high frequency (HF) teleseismic backprojection, we derive a comprehensive description of the rupture, which also predicts deep ocean tsunami wave heights. We further determine moment tensors and obtain accurate depth estimates for the aftershock sequence. The earthquake nucleated near the coast but then propagated to the north and updip, attaining a peak slip of 5-6 m. In contrast, HF seismic radiation is mostly emitted downdip of the region of intense slip and arrests earlier than the long period rupture, indicating smooth slip along the shallow plate interface in the final phase. A superficially similar earthquake in 1943 with a similar aftershock zone had a much shorter source time function, which matches the duration of HF seismic radiation in the recent event, indicating that the 1943 event lacked the shallow slip.
Mineral topsoils possess large organic carbon (OC) contents but there is only limited knowledge on the mechanisms controlling the preservation of organic matter (OM) against microbial decay. Samples were taken from the uppermost mineral topsoil horizon (0 to 5 cm) of seven sites under mature deciduous forest showing OC contents between 69 and 164 g kg(-1) and a wide range in mineral characteristics. At first, organic particles and the water-extractable OM were removed from the soil samples. Thereafter, Na-pyrophosphate extractable organic matter (OM(PY)), assumed to be indicative for OM bound via cation mediated interactions, and the OM remaining in the extraction residue (OM(ER)), supposed to be indicative for OM occluded in mechanically highly stable micro-aggregates, were sequentially separated and quantified. The composition of OM(PY) and OM(ER) was analyzed by FTIR and their stability by C-14 measurements. The OC remaining in the extraction residues accounted for 38 to 59% of the bulk soil OC (SOC) suggesting a much larger relevance of OM(ER) for the OM dynamic in the analyzed soils as compared with OM(PY) that accounted for 1.6 to 7.5% of the SOC. The FUR analyses revealed a lower relative proportion of C=O groups in OM(ER) compared to OM(PY) indicating differences in the degree of microbial processing between these fractions. Correlation analyses suggest an increase in the stability of OM(PY) with the soil pH and contents of Na-pyrophosphate soluble Fe, Al, and Mg and an increase in the stability of OM(ER) with the soil pH and the contents of clay and oxalate-soluble Fe and Al. Despite the detected influence of soil mineral characteristics on the turnover of OM(PY) and OM(ER), the Delta C-14 signatures indicated mean residence times less than 100 years. The presence of less stabilized OM in these fractions can be derived from methodological uncertainties and/or the fast cycling compartment of mineral-associated OM. (C) 2015 Elsevier B.V. All rights reserved.
Storm runoff from the Marikina River Basin frequently causes flood events in the Philippine capital region Metro Manila. This paper presents and evaluates a system to predict short-term runoff from the upper part of that basin (380km(2)). It was designed as a possible component of an operational warning system yet to be installed. For the purpose of forecast verification, hindcasts of streamflow were generated for a period of 15 months with a time-continuous, conceptual hydrological model. The latter was fed with real-time observations of rainfall. Both ground observations and weather radar data were tested as rainfall forcings. The radar-based precipitation estimates clearly outperformed the raingauge-based estimates in the hydrological verification. Nevertheless, the quality of the deterministic short-term runoff forecasts was found to be limited. For the radar-based predictions, the reduction of variance for lead times of 1, 2 and 3hours was 0.61, 0.62 and 0.54, respectively, with reference to a no-forecast scenario, i.e. persistence. The probability of detection for major increases in streamflow was typically less than 0.5. Given the significance of flood events in the Marikina Basin, more effort needs to be put into the reduction of forecast errors and the quantification of remaining uncertainties.
As field data on in-stream nitrate retention is scarce at catchment scales, this study aimed at quantifying net retention of nitrate within the entire river network of a fourth-order stream. For this purpose, a practical mass balance approach combined with a Lagrangian sampling scheme was applied and seasonally repeated to estimate daily in-stream net retention of nitrate for a 17.4 km long, agriculturally influenced, segment of the Steinlach River in southwestern Germany. This river segment represents approximately 70 % of the length of the main stem and about 32 % of the streambed area of the entire river network. Sampling days in spring and summer were biogeochemically more active than in autumn and winter. Results obtained for the main stem of Steinlach River were subsequently extrapolated to the stream network in the catchment. It was demonstrated that, for baseflow conditions in spring and summer, in-stream nitrate retention could sum up to a relevant term of the catchment’s nitrogen balance if the entire stream network was considered.
High-resolution topographic data greatly facilitate the remote identification of geomorphic features, furnishing valuable information concerning surface processes and characterization of reference markers for quantifying tectonic deformation. Marine terraces have been used as long baseline geodetic markers of relative past sea-level positions, reflecting the interplay between vertical crustal movements and sea-level oscillations. Uplift rates may be determined from the terrace age and the elevation of its shoreline angle, a geomorphic feature that can be correlated with past sea-levels positions. A precise definition of the shoreline angle in time and space is essential to obtain reliable uplift rates with coherent spatial correlation. To improve our ability to rapidly assess and map shoreline angles at regional and local scales, we have developed TerraceM, a MATLAB (R) graphical user interface that allows the shoreline angle and its associated error to be estimated using high-resolution topography. TerraceM uses topographic swath profiles oriented orthogonally to the terrace riser. Four functions are included to analyze the swath profiles and extract the shoreline angle, from both staircase sequences of multiple terraces and rough coasts characterized by eroded remnants of emerged terrace surfaces. The former are measured by outlining the paleocliffs and paieo-platforms and finding their intersection by extrapolating linear regressions, whereas the latter are assessed by automatically detecting peaks of sea-stack tops and back-projecting them to the modern sea cliff. In the absence of rigorous absolute age determinations of marine terraces, their geomorphic age may be estimated using previously published diffusion models. Postprocessing functions are included to obtain first-order statistics of shoreline-angle elevations and their spatial distribution. TerraceM has the ability to process series of profiles from several sites in an efficient and structured workflow. Results may be exported in Google Earth and ESRI shapefile formats. The precision and accuracy of the method have been estimated from a case study at Santa Cruz, California, by comparing TerraceM results with published field measurements. The repeatability was evaluated using multiple measurements made by inexperienced users. TerraceM will improve the efficiency and precision of estimating shoreline-angle elevations in wave-cut terraces in both marine and lacustrine environments.