Refine
Year of publication
Document Type
- Article (2599) (remove)
Language
- English (2599) (remove)
Keywords
- Holocene (35)
- climate change (31)
- Earthquake source observations (21)
- Climate change (18)
- permafrost (18)
- Himalaya (17)
- Pollen (17)
- Tibetan Plateau (17)
- erosion (17)
- Seismicity and tectonics (15)
Institute
- Institut für Geowissenschaften (2599) (remove)
Stable isotope paleoaltimetry makes use of systematic trends in the distribution and isotopic composition of modern precipitation with climate and topography, and of the potential to estimate the isotopic composition of paleoprecipitation from authigenic (in-situ formed) minerals. To illustrate the usefulness as well as potential limitations of this method, we review (1) processes controlling the isotopic composition of modem precipitation, (2) stable isotope data from modern precipitation across regions of high topography, and (3) stable isotope data from authigenic minerals that have been used to infer paleotopography. From this we conclude that stable isotope studies of authigenic minerals can permit useful inferences on paleotopography, with uncertainties that critically depend on a detailed understanding of local- to global-scale paleoclimate during the time interval of interest
The up to similar to4 km high southern Patagonian Andes form a pronounced topographic barrier to atmospheric circulation in the southern hemisphere westerlies, and cause one of the most drastic orographic rain shadows on earth. Geologic data imply that this climatic pattern has been established or significantly enhanced during Miocene surface uplift of this Andean segment. We report evidence for important climatic and ecologic changes in the eastern foreland of the Patagonian Andes that appear to be the result of this uplift. To provide constraints on Miocene plant ecosystems and precipitation in the eastern (leeward) foreland of the Patagonian Andes, we determined carbon and oxygen isotope values of pedogenic carbonate nodules from a similar to500 m thick section of the continental Santa Cruz Formation. The age of these deposits was constrained by Ar/Ar dating of intercalated tuffs, which range from similar to22 to 14 Ma. At similar to16.5 Ma, the delta(13)C values increase by similar to3parts per thousand, the delta(18)O values decrease by >2parts per thousand, and the scatter in the oxygen isotope data increases significantly. We interpret these changes as the consequence of >1 km surface uplift in this Andean segment (from the delta(18)O values), and increased aridity to its east (from the delta(13)C values and the increased scatter in the delta(18)O values). Sediments overlying the Santa Cruz Formation are very limited in extent and volume, and dominated by coarse conglomerates related to Pleistocene and older glaciations. It thus seems that, by similar to14 Ma, deposition in the eastern foreland of the Southern Patagonian Andes had essentially ceased as the result of rain shadow formation. (C) 2004 Elsevier B.V. All rights reserved
Three dimensional modelling of fractured and faulted reservoirs : framework and implementation
(2010)
Modelling of coupled physical processes in fractured and faulted media is a major challenge for the geoscience community. Due to the complexity related to the geometry of real fracture networks and fault systems, modelling studies have been mainly restricted either to two dimensional cases or to simplified orthogonal fracture systems consisting of vertical and horizontal fractures. An approach to generate three dimensional meshes for realistic fault geometries is presented. The method enables representation of faults in an arbitrary incline as two dimensional planes within a three dimensional, stratified porous matrix of a generic geometry. Based on a structural geological model, the method creates three dimensional unstructured tetrahedral meshes. These meshes can be used for finite element and finite volume numerical simulations. A simulation of a coupled fluid flow and heat transport problem for a two layered porous medium cut by two crossing faults is presented to test the reliability of the method.
Late Quaternary valley infill and dissection in the Indus River, western Tibetan Plateau margin
(2014)
The Indus, one of Earth's major rivers, drains large parts of the NW Himalaya and the Transhimalayan ranges that form part of the western Tibetan Plateau margin. In the western Himalayan syntaxis, where local topographic relief exceeds 7 km, the Indus has incised a steep bedrock gorge at rates of several mm yr(-1). Upstream, however, the upper Indus and its tributaries alternate between bedrock gorges and broad alluvial flats flanked by the Ladakh and Zanskar ranges. We review the late Quaternary valley history in this region with a focus on the confluence of the Indus and Zanskar Rivers, where vast alluvial terrace staircases and lake sediments record major episodes of aggradation and incision. New absolute dating of high-level fluvial terrace remnants using cosmogenic Be-10, optically and infrared stimulated luminescence (OSL, IRSL) indicates at least two phases of late Quaternary valley infilling. These phases commenced before similar to 200 ka and similar to 50-20 ka, judging from terrace treads stranded >150 m and similar to 30-40 m above modern river levels, respectively. Numerous stacks of lacustrine sediments that straddle the Indus River >200 km between the city of Leh and the confluence with the Shyok River share a distinct horizontal alignment. Constraints from IRSL samples of lacustrine sequences from the Leh-Spituk area reveal a protracted lake phase from >177 ka to 72 ka, locally accumulating >50-m thick deposits. In the absence of tectonic faulting, major lithological differences, and stream capture, we attribute the formation of this and other large lakes in the region to natural damming by large landslides, glaciers, and alluvial fans. The overall patchy landform age constraints from earlier studies can be reconciled by postulating a major deglacial control on sediment flux, valley infilling, and subsequent incision that has been modulated locally by backwater effects of natural damming. While comparison with Pleistocene monsoon proxies reveals no obvious correlation, a lateor post-glacial sediment pulse seems a more likely source of this widespread sedimentation that has partly buried the dissected bedrock topography. Overall, the long residence times of fluvial, alluvial and lacustrine deposits in the region (>500 ka) support previous studies, but remain striking given the dominantly steep slopes and deeply carved valleys that characterise this high-altitude mountain desert. Recalculated late Quaternary rates of fluvial bedrock incision in the Indus and Zanskar of 1.5 +/- 0.2 mm yr(-1) are at odds with the longevity of juxtaposed valley-fill deposits, unless a lack of decisive lateral fluvial erosion helps to preserve these late Pleistocene sedimentary archives. We conclude that alternating, similar to 10(4)-yr long, phases of massive infilling and incision have dominated the late Quaternary history of the Indus valley below the western Tibetan Plateau margin. (C) 2014 Elsevier Ltd. All rights reserved.
Rock glaciers in semiarid mountains contain large amounts of ice and might be important water stores aside from glaciers, lakes, and rivers. Yet whether and how rock glaciers interact with river channels in mountain valleys remains largely unresolved. We examine the potential for rock glaciers to block or disrupt river channels, using a new inventory of more than 2000 intact rock glaciers that we mapped from remotely sensed imagery in the Karakoram (KR), Tien Shan (TS), and Altai (ALT) mountains. We find that between 5% and 14% of the rock glaciers partly buried, blocked, diverted or constricted at least 95 km of mountain rivers in the entire study area. We use a Bayesian robust logistic regression with multiple topographic and climatic inputs to discern those rock glaciers disrupting mountain rivers from those with no obvious impacts. We identify elevation and potential incoming solar radiation (PISR), together with the size of feeder basins, as dominant predictors, so that lower-lying and larger rock glaciers from larger basins are more likely to disrupt river channels. Given that elevation and PISR are key inputs for modelling the regional distribution of mountain permafrost from the positions of rock-glacier toes, we infer that river-blocking rock glaciers may be diagnostic of non-equilibrated permafrost. Principal component analysis adds temperature evenness and wet-season precipitation to the controls that characterise rock glaciers impacting on rivers. Depending on the choice of predictors, the accuracy of our classification is moderate to good with median posterior area-under-the-curve values of 0.71-0.89. Clarifying whether rapidly advancing rock glaciers can physically impound rivers, or fortify existing dams instead, deserves future field investigation. We suspect that rock-glacier dams are conspicuous features that have a polygenetic history and encourage more research on the geomorphic coupling between permafrost lobes, river channels, and the sediment cascades of semiarid mountain belts. (c) 2018 John Wiley & Sons, Ltd.
Any understanding of sediment routing from mountain belts to their forelands and offshore sinks remains incomplete without estimates of intermediate storage that decisively buffers sediment yields from erosion rates, attenuates water and sediment fluxes, and protects underlying bedrock from incision. We quantify for the first time the sediment stored in > 38000 mainly postglacial Himalayan valley fills, based on an empirical volume-area scaling of valley-fill outlines automatically extracted from digital topographic data. The estimated total volume of 690(+452/-242) km(3) is mostly contained in few large valley fills > 1 km(3), while catastrophic mass wasting adds another 177(31) km(3). Sediment storage volumes are highly disparate along the strike of the orogen. Much of the Himalaya's stock of sediment is sequestered in glacially scoured valleys that provide accommodation space for similar to 44% of the total volume upstream of the rapidly exhuming and incising syntaxes. Conversely, the step-like long-wave topography of the central Himalayas limits glacier extent, and thus any significant glacier-derived storage of sediment away from tectonic basins. We show that exclusive removal of Himalayan valley fills could nourish contemporary sediment flux from the Indus and Brahmaputra basins for > 1 kyr, though individual fills may attain residence times of > 100 kyr. These millennial lag times in the Himalayan sediment routing system may sufficiently buffer signals of short-term seismic as well as climatic disturbances, thus complicating simple correlation and interpretation of sedimentary archives from the Himalayan orogen, its foreland, and its submarine fan systems. (C) 2013 Elsevier B.V. All rights reserved.
Mass wasting is an important process for denuding hillslopes and lowering ridge crests in active mountain belts such as the Himalaya-Karakoram ranges (HKR). Such a high-relief landscape is likely to be at its mechanical threshold, maintained by competing rapid rock uplift, river incision, and pervasive slope failure. We introduce excess topography, Z(E), for quantifying potentially unstable rock-mass volumes inclined at angles greater than a specified threshold angle. We find that Z(E) peaks along major fluvial and glacial inner gorges, which is also where the majority of 492 large (>0.1 km(2)) rock-slope failures occur in the Himalaya's largest cluster of documented Pleistocene to Holocene bedrock landslides. Our data reveal that bedrock landslides in the HKR chiefly detached from near or below the median elevation, whereas glaciers and rock glaciers occupy higher-elevation bands almost exclusively. Less than 10% of the area of the HKR is upslope of glaciers, such that possible censoring of evidence of large bedrock landslides above the permanent snow line barely affects this finding. Bedrock landslides appear to preferentially undermine topographic relief in response to fluvial and glacial incision along inner gorges, unless more frequent and smaller undetected failures, or rigorous (peri-)glacial erosion, compensate for this role at higher elevation. Either way, the distinct patterns of excess topography and large bedrock landsliding in the HKR juxtapose two stacked domains of landslide and (peri-)glacial erosion that may respond to different time scales of perturbation. Our findings call for more detailed analysis of vertical erosional domains and their geomorphic coupling in active mountain belts.
In this paper, we present an efficient kinematic ground-penetrating radar (GPR) surveying setup using a self- tracking total station (TTS). This setup combines the ability of modern GPR systems to interface with Global Positioning System (GPS) and the capability of the employed TTS system to immediately make the positioning information available in a standardized GPS data format. Wireless communication between the GPR and the TTS system is established by using gain variable radio modems. Such a kinematic surveying setup faces two major potential limitations. First, possible crosstalk effects between the GPR and the positioning system have to be evaluated. Based on multiple walkaway experiments, we show that, for reasonable field setups, instrumental crosstalk has no significant impact on GPR data quality. Second, we investigate systematic latency (i.e., the time delay between the actual position measurement by TTS and its fusion with the GPR data) and its impact on the positional precision of kinematically acquired 2-D and 3-D GPR data. To quantify latency for our kinematic survey setup, we acquired forward-reverse profile pairs across a well-known subsurface target. Comparing the forward and reverse GPR images using three fidelity measures allows determining the optimum latency value and correcting for it. Accounting for both of these potential limitations allows us to kinematically acquire high- quality and high-precision GPR data using off-the-shelf instrumentation without further hardware modifications. Until now, these issues have not been investigated in detail, and thus, we believe that our findings have significant implications also for other geophysical surveying approaches.
We have collected magnetic, 3D ground-penetrating-radar (GPR), and topographic data at an archaeological site within the Palace Garden of Paretz, Germany. The survey site covers an area of approximately 35 x 40 m across a hill structure (dips of up to 15 degrees) that is partly covered by trees. The primary goal of this study was to detect and locate the remains of ancient architectural elements, which, from historical records, were expected to be buried in the subsurface at this site. To acquire our geophysical data, we used a recently developed surveying approach that combines the magnetic and GPR instrument with a tracking total station (TTS). Besides efficient data acquisition, this approach provides positional information at an accuracy within the centimeter range. At the Paretz field site, this information was critical for processing and analyzing our geophysical data (in particular, GPR data) and enabled us to generate a high-resolution digital terrain model (DTM) of the surveyed area. Integrated analysis and interpretation based on composite images of the magnetic, 3D GPR, and high-resolution DTM data as well as selected attributes derived from these data sets allowed us to outline the remains of an artificial grotto and temple. Our work illustrates the benefit of using multiple surveying technologies, analyzing and interpreting the resulting data in an integrated fashion. It further demonstrates how modern surveying solutions allow for efficient, accurate data acquisition even in difficult terrain.
Three-dimensional (3D) ground-penetrating radar (GPR) represents an efficient high-resolution geophysical surveying method allowing to explore archaeological sites in a non-destructive manner. To effectively analyze large 3D GPR data sets, their combination with modern visualization techniques (e.g., 3D isoamplitude displays) has been acknowledged to facilitate interpretation beyond classical time-slice analysis. In this study, we focus on the application of data attributes (namely energy, coherency, and similarity), originally developed for petroleum reservoir related problems addressed by reflection seismology, to emphasize temporal and spatial variations within GPR data cubes. Based on two case studies, we illustrate the potential of such attribute based analyses towards a more comprehensive 3D GPR data interpretation. The main goal of both case studies was to localize and potentially characterize tombs inside medieval chapels situated in the state of Brandenburg, Germany. By comparing the calculated data attributes to the conventionally processed data cubes, we demonstrate the superior interpretability of the coherency and the similarity attribute for target identification and characterization.
Decomposition of geophysical signals (e.g., seismic and ground-penetrating radar data) into the time-frequency domain can provide valuable information for advanced interpretation (e.g., tuning effects) and processing (e.g., inverse Q-filtering). The quality of these subsequent processing steps is strongly related to the resolution of the selected time-frequency representation (TFR). In this study, we introduce a high-resolution spectral decomposition approach representing an extension of the recently proposed Tree-Based Pursuit (TBP) method. TBP significantly reduces the computational cost compared to the well known Matching Pursuit (MP) technique by introducing a tree structure prior to the actual matching procedure. Following the original implementation of TBP, we additionally incorporate waveforms commonly used in geophysical data processing and present an alternative approach to take phase shifts into account. Application of the proposed method to synthetic data and comparison of the results with other typically used decomposition approaches, illustrate the ability of our approach to provide decomposition results highly localized in both time and frequency. Applying our procedure to field GPR data illustrates its applicability to real data and provides examples for potential applications such as analyzing thin-bed responses and modulating the data frequency content.
The vertical radar profiling (VRP) technique uses surface-to-borehole acquisition geometries comparable to vertical seismic profiling (VSP). Major differences between the two methods do arise due to the fundamentally differing nature of the velocity-depth gradients and transmitter/receiver directivities. Largely for this reason, VRP studies have so far essentially been limited to the reconstruction of velocity-depth profiles by inverting direct arrival times from single-offset VRP surveys. In this study, we investigate the potential to produce high-resolution subsurface reflection images from multi-offset VRP data. Two synthetic data sets are used to evaluate a processing strategy suitably adapted from VSP processing. Despite the fundamental differences between VRP and VSP data, we found that our processing approach is capable of reconstructing subsurface structures of comparable complexity to those routinely imaged by VSP data. Finally, we apply our processing flow to two multi-offset VRP data sets recorded at a well constrained hydrogeophysical test site in SW-Germany. The inferred VRP images are compared with high-quality surface georadar reflection images and lithological logs available at the borehole locations. We find that the VRP images are in good agreement with the surface georadar data and reliably detect the major lithological boundaries. Due to the significantly shorter ray-paths, the depth penetration of the VRP data is, however, considerably higher than that of the surface georadar data. VRP reflection images thus provide an effective means for the depth-calibration and extension of conventional surface georadar data in the vicinity of boreholes.
Based on high-spatiotemporal-resolution data, the authors perform a climatological study of strong rainfall events propagating from southeastern South America to the eastern slopes of the central Andes during the monsoon season. These events account for up to 70% of total seasonal rainfall in these areas. They are of societal relevance because of associated natural hazards in the form of floods and landslides, and they form an intriguing climatic phenomenon, because they propagate against the direction of the low-level moisture flow from the tropics. The responsible synoptic mechanism is analyzed using suitable composites of the relevant atmospheric variables with high temporal resolution. The results suggest that the low-level inflow from the tropics, while important for maintaining sufficient moisture in the area of rainfall, does not initiate the formation of rainfall clusters. Instead, alternating low and high pressure anomalies in midlatitudes, which are associated with an eastward-moving Rossby wave train, in combination with the northwestern Argentinean low, create favorable pressure and wind conditions for frontogenesis and subsequent precipitation events propagating from southeastern South America toward the Bolivian Andes.
Extreme Rainfall of the South American Monsoon System: A Dataset Comparison Using Complex Networks
(2015)
In this study, the authors compare six different rainfall datasets for South America with a focus on their representation of extreme rainfall during the monsoon season (December February): the gauge-calibrated TRMM 3B42 V7 satellite product; the near-real-time TRMM 3B42 V7 RT, the GPCP 1 degrees daily (1DD) V1.2 satellite gauge combination product, the Interim ECMWF Re-Analysis (ERA-Interim) product; output of a high-spatial-resolution run of the ECHAM6 global circulation model; and output of the regional climate model Eta. For the latter three, this study can be understood as a model evaluation. In addition to statistical values of local rainfall distributions, the authors focus on the spatial characteristics of extreme rainfall covariability. Since traditional approaches based on principal component analysis are not applicable in the context of extreme events, they apply and further develop methods based on complex network theory. This way, the authors uncover substantial differences in extreme rainfall patterns between the different datasets: (i) The three model-derived datasets yield very different results than the satellite gauge combinations regarding the main climatological propagation pathways of extreme events as well as the main convergence zones of the monsoon system. (ii) Large discrepancies are found for the development of mesoscale convective systems in southeastern South America. (iii) Both TRMM datasets and ECHAM6 indicate a linkage of extreme rainfall events between the central Amazon basin and the eastern slopes of the central Andes, but this pattern is not reproduced by the remaining datasets. The authors' study suggests that none of the three model-derived datasets adequately captures extreme rainfall patterns in South America.
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.
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.
Climatic observables are often correlated across long spatial distances, and extreme events, such as heatwaves or floods, are typically assumed to be related to such teleconnections(1,2). Revealing atmospheric teleconnection patterns and understanding their underlying mechanisms is of great importance for weather forecasting in general and extreme-event prediction in particular(3,4), especially considering that the characteristics of extreme events have been suggested to change under ongoing anthropogenic climate change(5-8). Here we reveal the global coupling pattern of extreme-rainfall events by applying complex-network methodology to high-resolution satellite data and introducing a technique that corrects for multiple-comparison bias in functional networks. We find that the distance distribution of significant connections (P < 0.005) around the globe decays according to a power law up to distances of about 2,500 kilometres. For longer distances, the probability of significant connections is much higher than expected from the scaling of the power law. We attribute the shorter, power-law-distributed connections to regional weather systems. The longer, super-power-law-distributed connections form a global rainfall teleconnection pattern that is probably controlled by upper-level Rossby waves. We show that extreme-rainfall events in the monsoon systems of south-central Asia, east Asia and Africa are significantly synchronized. Moreover, we uncover concise links between south-central Asia and the European and North American extratropics, as well as the Southern Hemisphere extratropics. Analysis of the atmospheric conditions that lead to these teleconnections confirms Rossby waves as the physical mechanism underlying these global teleconnection patterns and emphasizes their crucial role in dynamical tropical-extratropical couplings. Our results provide insights into the function of Rossby waves in creating stable, global-scale dependencies of extreme-rainfall events, and into the potential predictability of associated natural hazards.
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.
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.
Measuring (dis)similarity between ecosystem states is a key theme in ecology. Much of community and ecosystem ecology is devoted to searching for patterns in ecosystem similarity from an external observer's viewpoint, using variables such as species abundances, measures of diversity and complexity. However, from the point of view of organisms in the ecosystem, proportional population growth rates are the only relevant aspect of ecosystem state, because natural selection acts on groups of organisms with different proportional population growth rates. We therefore argue that two ecosystem states are equivalent if and only if, for each species they contain, the proportional population growth rate does not differ between the states. Based on this result, we develop species-level and aggregated summary measures of ecosystem state and discuss their ecological meaning. We illustrate our approach using a long-term dataset on the plankton community from the Central European Lake Constance. We show that the first three principal components of proportional population growth rates describe most of the variation in ecosystem state in Lake Constance. We strongly recommend using proportional population growth rates and the derived equivalence classes for comparative ecosystem studies. This opens up new perspectives on important existing topics such as alternative stable ecosystem states, community assembly, and the processes generating regularities in ecosystems.
The Upper Cretaceous La Cova limestones (southern Pyrenees, Spain) host a rich and diverse larger foraminiferal fauna, which represents the first diversification of K-strategists after the mass extinction at the Cenomanian-Turonian boundary.
The stratigraphic distribution of the main taxa of larger foraminifera defines two assemblages. The first assemblage is characterised by the first appearance of lacazinids (Pseudolacazina loeblichi) and mean-dropsinids (Eofallotia simplex), by the large agglutinated Montsechiana montsechiensis, and by several species of complex rotalids (Rotorbinella campaniola, Iberorotalia reicheli, Orbitokhatina wondersmitti and Calcarinella schaubi). The second assemblage is defined by the appearance of Lacazina pyrenaica, Palandrosina taxyae and Martiguesia cyclamminiformis.
A late Coniacian-early Santonian age was so far accepted for the La Cova limestones, based on indirect correlation with deep-water fades bearing planktic foraminifers of the Dicarinella concavata zone. Strontium isotope stratigraphy, based on many samples of pristine biotic calcite of rudists and ostreids, indicates that the La Cova limestones span from the early Coniacian to the early-middle Santonian boundary. The first assemblage of larger foraminifera appears very close to the early-middle Coniacian boundary and reaches its full diversity by the middle Coniacian. The originations defining the second assemblage are dated as earliest Santonian: they represent important bioevents to define the Coniacian-Santonian boundary in the shallow-water facies of the South Pyrenean province.
By means of the calibration of strontium isotope stratigraphy to the Geological Time Scale, the larger foraminiferal assemblages of the La Cova limestones can be correlated to the standard biozonal scheme of ammonites, planktonic foraminifers and calcareous nannoplankton. This correlation is a first step toward a larger foraminifera standard biozonation for Upper Cretaceous carbonate platform facies.
Soil pollutants such as hydrocarbons can induce toxic effects in plants and associated arbuscular mycorrhizal fungi (AMF). This study was conducted to evaluate if the legume Lotus corniculatus and the grass Elymus trachycaulus and arbuscular mycorrhizal fungi could grow in two oil sands processing by-products after bitumen extraction from the oil sands in northern Alberta, Canada. Substrate treatments were coarse tailings sand (CTS), a mix of dry mature fine tailings (MFT) with CTS (1: 1) and Pleistocene sandy soil (hydrocarbon free); microbial treatments were without AMF, with AMF and AMF plus soil bacteria isolated from oil sands reclamation sites. Plant biomass, root morphology, leaf water content, shoot tissue phosphorus content and mycorrhizal colonization were evaluated. Both plant species had reduced growth in CTS and tailings mix relative to sandy soil. AMF frequency and intensity in roots of E. trachycaulus was not influenced by soil hydrocarbons; however, it decreased significantly over time in roots of L. corniculatus without bacteria in CTS. Mycorrhizal inoculation alone did not significantly improve plant growth in CTS and tailings mix; however, inoculation with mycorrhizae plus bacteria led to a significantly positive response of both plant species in CTS. Thus, combined inoculation with selected mycorrhizae and bacteria led to synergistic effects. Such combinations may be used in future to improve plant growth in reclamation of CTS and tailings mix.
A facile one-pot synthesis for preparing thermosensitive, luminescent gold nanodots with diameters of 1-2 nm is presented. The influence of the alkyl chain length of the surface ligands (alkyl thiols) on the optical properties of the gold nanodots was investigated. The synthesized gold nanodots show strong thermosensitive photoluminescence. A photoluminescence quantum yield of 16.6% was observed at room temperature, which could be improved to a value of 28.6% when cooling the gold nanodot solutions to -7 degrees C. The synthesized thermosensitive, luminescent gold nanodots are interesting candidates for optoelectronic devices, medical imaging, sensing, or security labels.
A SSHAC Level 3 Probabilistic Seismic Hazard Analysis for a New-Build Nuclear Site in South Africa
(2015)
A probabilistic seismic hazard analysis has been conducted for a potential nuclear power plant site on the coast of South Africa, a country of low-to-moderate seismicity. The hazard study was conducted as a SSHAC Level 3 process, the first application of this approach outside North America. Extensive geological investigations identified five fault sources with a non-zero probability of being seismogenic. Five area sources were defined for distributed seismicity, the least active being the host zone for which the low recurrence rates for earthquakes were substantiated through investigations of historical seismicity. Empirical ground-motion prediction equations were adjusted to a horizon within the bedrock at the site using kappa values inferred from weak-motion analyses. These adjusted models were then scaled to create new equations capturing the range of epistemic uncertainty in this region with no strong motion recordings. Surface motions were obtained by convolving the bedrock motions with site amplification functions calculated using measured shear-wave velocity profiles.
Logic trees are widely used in probabilistic seismic hazard analysis as a tool to capture the epistemic uncertainty associated with the seismogenic sources and the ground-motion prediction models used in estimating the hazard. Combining two or more ground-motion relations within a logic tree will generally require several conversions to be made, because there are several definitions available for both the predicted ground-motion parameters and the explanatory parameters within the predictive ground-motion relations. Procedures for making conversions for each of these factors are presented, using a suite of predictive equations in current use for illustration. The sensitivity of the resulting ground-motion models to these conversions is shown to be pronounced for some of the parameters, especially the measure of source-to-site distance, highlighting the need to take into account any incompatibilities among the selected equations. Procedures are also presented for assigning weights to the branches in the ground-motion section of the logic tree in a transparent fashion, considering both intrinsic merits of the individual equations and their degree of applicability to the particular application
Major earthquakes ( M > 8) have repeatedly ruptured the Nazca-South America plate interface of south-central Chile involving meter scale land-level changes. Earthquake recurrence intervals, however, extending beyond limited historical records are virtually unknown, but would provide crucial data on the tectonic behavior of forearcs. We analyzed the spatiotemporal pattern of Holocene earthquakes on Santa Maria Island (SMI; 37 degrees S), located 20 km off the Chilean coast and approximately 70 km east of the trench. SMI hosts a minimum of 21 uplifted beach berms, of which a subset were dated to calculate a mean uplift rate of 2.3 +/- 0.2 m/ky and a tilting rate of 0.022 +/- 0.002 degrees/ky. The inferred recurrence interval of strandline-forming earthquakes is similar to 180 years. Combining coseismic uplift and aseismic subsidence during an earthquake cycle, the net gain in strandline elevation in this environment is similar to 0.4 m per event
Aluminum-26 and beryllium-10 surface exposure dating on cut-and-fill river-terrace surfaces from the lower Sutlej Valley (northwest Himalaya) documents the close link between Indian Summer Monsoon (ISM) oscillations and intervals of enhanced fluvial incision. During the early Holocene ISM optimum, precipitation was enhanced and reached far into the internal parts of the orogen. The amplified sediment flux from these usually dry but glaciated areas caused alluviation of downstream valleys up to 120 m above present grade at ca. 9.9 k.y. B.P. Terrace formation (i.e., incision) in the coarse deposits occurred during century-long weak ISM phases that resulted in reduced moisture availability and most likely in lower sediment flux. Here, we suggest that the lower sediment flux during weak ISM phases allowed rivers to incise episodically into the alluvial fill
Hydrological modelling of a Pleistocene landslide-dammed lake in the Santa Maria Basin, NW Argentina
(2001)
The tectonic and climatic boundary conditions of the broken foreland and the orogen interior of the southern Central Andes of northwestern Argentina cause strong contrasts in elevation, rainfall, and surface-process regimes. The climatic gradient in this region ranges from the wet, windward eastern flanks (similar to 2 m/yr rainfall) to progressively drier western basins and ranges (similar to 0.1 m/yr) bordering the arid Altiplano-Puna Plateau. In this study, we analyze the impact of spatiotemporal climatic gradients on surface erosion: First, we present 41 new catchment-mean erosion rates derived from cosmogenic nuclide inventories to document spatial erosion patterns. Second, we re-evaluate paleoclimatic records from the Calchaquies basin (66 W, 26 S), a large intermontane basin bordered by high (> 4.5 km) mountain ranges, to demonstrate temporal variations in erosion rates associated with changing climatic boundary conditions during the late Pleistocene and Holocene. Three key observations in this region emphasize the importance of climatic parameters on the efficiency of surface processes in space and time: (1) First-order spatial patterns of erosion rates can be explained by a simple specific stream power (SSP) approach. We explicitly account for discharge by routing high-resolution, satellite derived rainfall. This is important as the steep climatic gradient results in a highly non-linear relation between drainage area and discharge. This relation indicates that erosion rates (ER) scale with ER similar to SSP1.4 on cosmogenic-nuclide time scales. (2) We identify an intrinsic channel-slope behavior in different climatic compartments. Channel slopes in dry areas (< 0.25 m/yr rainfall) are slightly steeper than in wet areas (> 0.75 m/yr) with equal drainage areas, thus compensating lower amounts of discharge with steeper slopes. (3) Erosion rates can vary by an order of magnitude between presently dry (similar to 0.05 mm/yr) and well-defined late Pleistocene humid (similar to 0.5 mm/yr) conditions within an intemontane basin. Overall, we document a strong climatic impact on erosion rates and channel slopes. We suggest that rainfall reaching areas with steeper channel slopes in the orogen interior during wetter climate periods results in intensified sediment mass transport, which is primarily responsible for maintaining the balance between surface uplift, erosion, sediment routing and transient storage in the orogen.
The interplay between topography and Indian summer monsoon circulation profoundly controls precipitation distribution, sediment transport, and river discharge along the Southern Himalayan Mountain Front (SHF). The Higher Himalayas form a major orographic barrier that separates humid sectors to the south and and regions to the north. During the Indian summer monsoon, vortices transport moisture from the Bay of Bengal, swirl along the SHF to the northwest, and cause heavy rainfall when colliding with the mountain front. In the eastern and central parts of the Himalaya, precipitation measurements derived from passive microwave analysis (SSM/I) show a strong gradient, with high values at medium elevations and extensive penetration of moisture along major river valleys into the orogen. The end of the monsoonal conveyer belt is near the Sutlej Valley in the NW Himalaya, where precipitation is lower and rainfall maxima move to lower elevations. This region thus comprises a climatic transition zone that is very sensitive to changes in Indian summer monsoon strength. To constrain magnitude, temporal, and spatial distribution of precipitation, we analyzed high-resolution passive microwave data from the last decade and identified an abnormal monsoon year (AMY) in 2002. During the 2002 AMY, violent rainstorms conquered orographic barriers and penetrated far into otherwise and regions in the northwest Himalaya at elevations in excess of 3 km asl. While precipitation in these regions was significantly increased and triggered extensive erosional processes (i.e., debris flows) on sparsely vegetated, steep hillslopes, mean rainfall along the low to medium elevations was not significantly greater in magnitude. This shift may thus play an important role in the overall sediment flux toward the Himalayan foreland. Using extended precipitation and sediment flux records for the last century, we show that these events have a decadal recurrence interval during the present-day monsoon circulation. Hence, episodically occurring AMYs control geomorphic processes primarily in the high-elevation and sectors of the orogen, while annual recurring monsoonal rainfall distribution dominates erosion in the low- to medium- elevation parts along the SHF. (C) 2004 Elsevier B.V. All rights reserved
Late Quaternary intensified monsoon phases control landscape evolution in the northwest Himalaya
(2005)
The intensity of the Asian summer-monsoon circulation varies over decadal to millennial time scales and is reflected in changes in surface processes, terrestrial environments, and marine sediment records. However, the mechanisms of long-lived (2-5 k.y.) intensified monsoon phases, the related changes in precipitation distribution, and their effect on landscape evolution and sedimentation rates are not yet well understood. The and high-elevation sectors of the orogen correspond to a climatically sensitive zone that currently receives rain only during abnormal (i.e., strengthened) monsoon seasons. Analogous to present-day rainfall anomalies, enhanced precipitation during an intensified monsoon phase is expected to have penetrated far into these geomorphic threshold regions where hillslopes are close to the angle of failure. We associate landslide triggering during intensified monsoon phases with enhanced precipitation, discharge, and sediment flux leading to an increase in pore-water pressure, lateral scouring of rivers, and over- steepening of hillslopes, eventually resulting in failure of slopes and exceptionally large mass movements. Here we use lacustrine deposits related to spatially and temporally clustered large landslides (>0.5 km(3)) in the Sutlej Valley region of the northwest Himalaya to calculate sedimentation rates and to infer rainfall patterns during late Pleistocene (29-24 ka) and Holocene (10-4 ka) intensified monsoon phases. Compared to present-day sediment-flux measurements, a fivefold increase in sediment-transport rates recorded by sediments in landslide-dammed lakes characterized these episodes of high climatic variability. These changes thus emphasize the pronounced imprint of millennial-scale climate change on surface processes and landscape evolution
Adjustment of median ground motion prediction equations (GMPEs) from one region to another region is one of the major challenges within the current practice of seismic hazard analysis. In our approach of generating response spectra, we derive two separate empirical models for a) Fourier amplitude spectrum (FAS) and b) duration of ground motion. To calculate response spectra, the two models are combined within the random vibration theory (RVT) framework. The models are calibrated on recordings obtained from shallow crustal earthquakes in active tectonic regions. We use a subset of NGA-West2 database with M3.2-7.9 earthquakes at distances 0-300 km. The NGA-West2 database expanded over a wide magnitude and distance range facilitates a better constraint over derived models. A frequency-dependent duration model is derived to obtain adjustable response spectral ordinates. Excellent comparison of our approach with other NGA-West2 models implies that it can also be used as a stand-alone model.
We have analyzed the recently developed pan-European strong motion database, RESORCE-2012: spectral parameters, such as stress drop (stress parameter, Delta sigma), anelastic attenuation (Q), near surface attenuation (kappa(0)) and site amplification have been estimated from observed strong motion recordings. The selected dataset exhibits a bilinear distance-dependent Q model with average kappa(0) value 0.0308 s. Strong regional variations in inelastic attenuation were also observed: frequency-independent Q(0) of 1462 and 601 were estimated for Turkish and Italian data respectively. Due to the strong coupling between Q and kappa(0), the regional variations in Q have strong impact on the estimation of near surface attenuation kappa(0). kappa(0) was estimated as 0.0457 and 0.0261 s for Turkey and Italy respectively. Furthermore, a detailed analysis of the variability in estimated kappa(0) revealed significant within-station variability. The linear site amplification factors were constrained from residual analysis at each station and site-class type. Using the regional Q(0) model and a site-class specific kappa(0), seismic moments (M-0) and source corner frequencies f (c) were estimated from the site corrected empirical Fourier spectra. Delta sigma did not exhibit magnitude dependence. The median Delta sigma value was obtained as 5.75 and 5.65 MPa from inverted and database magnitudes respectively. A comparison of response spectra from the stochastic model (derived herein) with that from (regional) ground motion prediction equations (GMPEs) suggests that the presented seismological parameters can be used to represent the corresponding seismological attributes of the regional GMPEs in a host-to-target adjustment framework. The analysis presented herein can be considered as an update of that undertaken for the previous Euro-Mediterranean strong motion database presented by Edwards and Fah (Geophys J Int 194(2):1190-1202, 2013a).
The functional form of empirical response spectral ground-motion prediction equations (GMPEs) is often derived using concepts borrowed from Fourier spectral modeling of ground motion. As these GMPEs are subsequently calibrated with empirical observations, this may not appear to pose any major problems in the prediction of ground motion for a particular earthquake scenario. However, the assumption that Fourier spectral concepts persist for response spectra can lead to undesirable consequences when it comes to the adjustment of response spectral GMPEs to represent conditions not covered in the original empirical data set. In this context, a couple of important questions arise, for example, what are the distinctions and/or similarities between Fourier and response spectra of ground motions? And, if they are different, then what is the mechanism responsible for such differences and how do adjustments that are made to Fourier amplitude spectrum (FAS) manifest in response spectra? The present article explores the relationship between the Fourier and response spectrum of ground motion by using random vibration theory (RVT). With a simple Brune (1970, 1971) source model, RVT-generated acceleration spectra for a fixed magnitude and distance scenario are used. The RVT analyses reveal that the scaling of low oscillator-frequency response spectral ordinates can be treated as being equivalent to the scaling of the corresponding Fourier spectral ordinates. However, the high oscillator-frequency response spectral ordinates are controlled by a rather wide band of Fourier spectral ordinates. In fact, the peak ground acceleration, counter to the popular perception that it is a reflection of the high-frequency characteristics of ground motion, is controlled by the entire Fourier spectrum of ground motion. Additionally, this article demonstrates how an adjustment made to FAS is similar or different to the same adjustment made to response spectral ordinates. For this purpose, two cases: adjustments to the stress parameter (Delta sigma) (source term), and adjustments to the attributes reflecting site response (V-S - kappa(0)) are considered.
Empirical ground-motion prediction equations (GMPEs) require adjustment to make them appropriate for site-specific scenarios. However, the process of making such adjustments remains a challenge. This article presents a holistic framework for the development of a response spectral GMPE that is easily adjustable to different seismological conditions and does not suffer from the practical problems associated with adjustments in the response spectral domain. The approach for developing a response spectral GMPE is unique, because it combines the predictions of empirical models for the two model components that characterize the spectral and temporal behavior of the ground motion. Essentially, as described in its initial form by Bora et al. (2014), the approach consists of an empirical model for the Fourier amplitude spectrum (FAS) and a model for the ground-motion duration. These two components are combined within the random vibration theory framework to obtain predictions of response spectral ordinates. In addition, FAS corresponding to individual acceleration records are extrapolated beyond the useable frequencies using the stochastic FAS model, obtained by inversion as described in Edwards and Fah (2013a). To that end, a (oscillator) frequency-dependent duration model, consistent with the empirical FAS model, is also derived. This makes it possible to generate a response spectral model that is easily adjustable to different sets of seismological parameters, such as the stress parameter Delta sigma, quality factor Q, and kappa kappa(0). The dataset used in Bora et al. (2014), a subset of the RESORCE-2012 database, is considered for the present analysis. Based upon the range of the predictor variables in the selected dataset, the present response spectral GMPE should be considered applicable over the magnitude range of 4 <= M-w <= 7.6 at distances <= 200 km.
The actual evapotranspiration is an important, but difficult to determine, element in the water balance of lakes and their catchment areas. Reliable data on evapotranspiration are not available for most lake basins for which paleoclimate reconstructions and modeling have been performed, particularly those in remote parts of Africa. We have used thermal infrared multispectral data for 14 ASTER scenes from the TERRA satellite to estimate the actual evapotranspiration in the 12,800 km(2) catchment of the Suguta Valley, northern Kenya Rift Evidence from sediments and paleo-shorelines indicates that, during the African Humid Period (AHP, 14.8 to 5.5 kyrs BP), this valley contained a large lake, 280 m deep and covering similar to 2200 km(2), which has now virtually disappeared. Evapotranspiration estimates for the Suguta Basin were generated using the Surface Energy Balance Algorithm for Land (SEBAL). Climate data required for the model were extracted from a high-resolution gridded dataset obtained from the Climatic Research Unit (East Anglia, UK). Results suggest significant spatial variations in evapotranspiration within the catchment area (ranging from 450 mm/yr in the basin to the north to 2000 mm/yr in more elevated areas) and precipitation that was similar to 20% higher during the AHP than in recent times. These results are in agreement with other estimates of paleo-precipitation in East Africa. The extreme response of the lake system (similar to 280 m greater water depth than today, and a lake surface area of 2200 km(2)) to only moderately higher precipitation illustrates the possible sensitivity of this area to future climate change.
Barium, lanthanum, ytterbium, and yttrium partitioning experiments between fluid-saturated haplogranitic melts and aqueous solutions were conducted at 750 to 950 degrees C and 0.2 to 1 GPa to investigate the effects of melt and fluid composition, pressure, and temperature. Partition coefficients were determined using different experimental methods. On one hand quenched experiments were performed, and on the other hand, trace element contents in the aqueous fluid were determined directly using a hydrothermal diamond-anvil cell and synchrotron radiation X-ray fluorescence microanalysis of K-lines. The latter required a high excitation energy of 50 key due to the high energies necessary to excite the K-lines of the studied elements. The data from these two techniques showed good agreement for chloridic solutions, whereas quenching had a significant effect on results of the experiments with only water in the case of Ba. In Cl-free experiments, lanthanum and yttrium, trace element contents were even below detection limit in the quenched fluids, whereas small concentrations were detected in comparable in-situ experiments. This distinct difference is likely due to back reactions between fluid and melt upon cooling. The partitioning data of all elements show no dependence on the temperature and only small dependence on pressure. In contrast, the partitioning is strongly influenced by the composition of the starting fluid and melt. For chloridic fluids, there was a sharp increase in the Ba, La, Y and Yb partition coefficients with the alumina saturation index (ASI). The Ba partition coefficient increased from 0.002 at an ASI of 0.8 to 0.55 at an ASI of 1.07. At higher ASI, it decreased slightly to 0.2 at an ASI of similar to 1.3. Likewise, it was one to two orders of magnitude higher in chloridic fluids compared to those found in H2O experiments. Fluid-melt partition coefficients of La and Y increased from 0.002 at an ASI of similar to 0.8 to similar to 0.1 at an ASI of 1.2. In the same ASI range, the Yb partition coefficient increased to a maximum value of 0.02. Even at high salinities all elements fractionate into the melt. The compositional dependence of the partitioning data imply that both melt composition and fluid composition have a strong influence on trace element behavior and that complexation of Ba. REE and Y tin the fluid is not only controlled by the presence of Cl- in the fluid. Instead, interaction of these elements with major melt components dissolved in the fluid is very likely.
Trace element concentrations in aqueous fluids in equilibrium with haplogranitic melt were determined in situ at elevated P-T conditions using hydrothermal diamond-anvil cells and synchrotron-radiation XRF microanalyses. Time- resolved analyses showed that the Rb and Sr concentrations in the fluids became constant in less than 2000 s at all temperatures (500 to 780 degrees C). Although fluid-melt equilibration was very rapid, the change in the concentration of both elements in the fluid with temperature was fairly small (a slight increase for Rb and a slight decrease for Sr). This permitted partitioning data for Rb and Sr between haplogranitic melt and H2O or NaCl+KCl+HCl aqueous solutions at 750 degrees C and 200 to 700 MPa to be obtained from EMP analyses of the quenched melt and the in situ SR-XRF analyses of the equilibrated fluid. The resulting D-Rb(f/m) and D-Sr(f/m) were 0.01 +/- 0.002 and 0.006 +/- 0.001 for water as starting fluid, and increased to 0.47 +/- 0.08 and 0.23 +/- 0.03 for 3.56 m (NaCl+KCl)+0.04 in HCl at pressures of 224 to 360 MPa. In the experiments with H2O as starting fluid, the partition coefficients increased with pressure, i.e. D- Rb(f/m) from 0.01 +/- 0.002 to 0.22 +/- 0.02 and D-Sr(f/m) from 0.006 0.001 to 0.02 +/- 0.005 with a change in pressure from 360 to 700 MPa. At pressures to 360 MPa, the Rb/Sr ratio in the fluid was found to be independent of the initial salt concentration (Rb/Sr = 1.45 +/- 0.6). This ratio increased to 7.89 +/- 1.95 at 700 MPa in experiments with chloride free fluids, which indicates different changes in the Rb and Sr speciation with pressure.
We study optimal and equilibrium sizes of cities in a city system model with pollution. Pollution is a function of population size. If pollution is local or per-capita pollution increases with population, equilibrium cities are too large under symmetry; with asymmetric cities, the largest cities are too large and the smallest too small. When pollution is global and per-capita pollution declines with city size, cities may be too small under symmetry; with asymmetric cities, the largest cities are too small and the smallest too large if the marginal damage of pollution is large enough. We calibrate the model to US cities and find that the largest cities may be undersized by 3-4%.
Garnet of eclogite (formerly termed garnet clinopyroxenite) hosted in lenses of orogenic garnet peridotite from the Granulitgebirge, NW Bohemian Massif, contains unique inclusions of granitic melt, now either glassy or crystallized. Analysed glasses and re‐homogenized inclusions are hydrous, peraluminous, and enriched in highly incompatible elements characteristic of the continental crust such as Cs, Li, B, Pb, Rb, Th, and U. The original melt thus represents a pristine, chemically evolved metasomatic agent, which infiltrated the mantle via deep continental subduction during the Variscan orogeny. The bulk chemical composition of the studied eclogites is similar to that of Fe‐rich basalt and the enrichment in LILE and U suggest a subduction‐related component. All these geochemical features confirm metasomatism. In comparison with many other garnet+clinopyroxene‐bearing lenses in peridotites of the Bohemian Massif, the studied samples from Rubinberg and Klatschmühle are more akin to eclogite than pyroxenites, as reflected in high jadeite content in clinopyroxene, relatively low Mg, Cr, and Ni but relatively high Ti. However, trace elements of both bulk rock and individual mineral phases show also important differences making these samples rather unique. Metasomatism involving a melt requiring a trace element pattern very similar to the composition reported here has been suggested for the source region of rocks of the so‐called durbachite suite, that is, ultrapotassic melanosyenites, which are found throughout the high‐grade Variscan basement. Moreover, the Th, U, Pb, Nb, Ta, and Ti patterns of these newly studied melt inclusions (MI) strongly resemble those observed for peridotite and its enclosed pyroxenite from the T‐7 borehole (Staré, České Středhoři Mountains) in N Bohemia. This suggests that a similar kind of crustal‐derived melt also occurred here. This study of granitic MI in eclogites from peridotites has provided the first direct characterization of a preserved metasomatic melt, possibly responsible for the metasomatism of several parts of the mantle in the Variscides.
Granitic melt inclusions were found in layers of garnet clinopyroxenites from orogenic peridotites hosted in high-pressure felsic granulites of the Granulitgebirge, central Europe. The inclusions are both glassy and crystallized, and occur as clusters in the garnet. Microstructural features suggest that the inclusions formed while garnet was growing as a peritectic phase, likely alongside clinopyroxene. The chemistry of the melt, in particular its trace element signature, shows a crustal contribution, probably due to the involvement of phengite in the melt-producing reaction, most likely in the presence of a fluid. The presence of a granitoid melt in mantle rocks may be the result of localized melting of a phengite-bearing protolith either already present in the peridotites or, more likely, within the local deeply subducted crustal units. In the latter case, the melt would have infiltrated the peridotites and generated pyroxenite via metasomatism. In either case, the presence of granitoid inclusions in orogenic peridotite provides direct evidence for a genetic connection between a high-pressure crustal melt and garnet pyroxenites. The in situ characterization of these remnants of natural melt provides direct quantitative constraints on (one of) the agents responsible for the interaction between crust and mantle.
This paper analyses the effect of spatial resolution and distribution of model input data on the results of regional-scale land use scenarios using three different hydrological catchment models. A 25 m resolution data set of a mesoscale catchment and three land use scenarios are used. Data are systematically aggregated to resolutions up to 2 kill. Land use scenarios are spatially redistributed, both randomly and topography based. Using these data, water fluxes are calculated on a daily time step for a 16 year time period without further calibration. Simulation results are used to identify grid size, distribution and model dependent scenario effects. In the case of data aggregation, all applied models react sensitively to grid size. WASIM and TOPLATS simulate constant water balances for grid sizes from 50 m to 300-500 m, SWAT is more sensitive to input data aggregation, simulating constant water balances between 50 m and 200 m grid size. The calculation of scenario effects is less robust to data aggregation. The maximum acceptable grid size reduces to 200-300 m for TOPLATS and WASIM. In case of spatial distribution, SWAT and TOPLATS are slightly sensitive to a redistribution of land use (below 1.5% for water balance terms), whereas WASIM shows almost no reaction. Because the aggregation effects were stronger than the redistribution effects, it is concluded that spatial discretisation is more important than spatial distribution. As the aggregation effect was mainly associated with a change in land use fraction, it is concluded that accuracy of data sets is much more important than a high spatial resolution.
The effect of water activity on the oxidation and structural state of Fe in a ferro-basaltic melt
(2005)
Experimental investigations have been performed at T = 1200 degrees C, P = 200 MPa and fH(2) corresponding to H2O-MnO-Mn3O4 and H2O-QFM redox buffers to study the effect of H2O activity on the oxidation and structural state of Fe in an iron-rich basaltic melt. The analysis of Mossbauer and Fe K-edge X-ray absorption nearedge structure (XANES) spectra of the quenched hydrous ferrobasaltic glasses shows that the Fe3+/Sigma Fe ratio of the glass is directly related to aH(2)O in a H-2-buffered system and, consequently, to the prevailing oxygen fugacity (through the reaction of water dissociation H2O <-> H-2 + 1/2 O-2). However, water as a chemical component of the silicate melt has an indistinguishable effect on the redox state of iron at studied conditions. The experimentally obtained relationship between fO(2) and Fe3+/Fe2+ in the hydrous ferrobasaltic melt can be adequately predicted in the investigated range by the existing empiric and thermodynamic models. The ratio of ferric and ferrous Fe is proportional to the oxygen fugacity to the power of similar to 0.25 which agrees with the theoretical value from the stoichiometry of the Fe redox reaction (FeO + 1/4 O-2 = FeO1.5). The mean centre shifts for Fe2+ and Fe3+ absorption doublets in Mossbauer spectra show little change with increasing Fe3+/Sigma Fe, suggesting no significant change in the type of iron coordination. Similarly, XANES preedge spectra indicate a mixed (C3h, Td, and Oh, i.e., 5-, 4-, and sixfold) coordination of Fe in hydrous basaltic glasses. Copyright (c) 2005 Elsevier Ltd
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.
Subduction factory : 1. Theoretical mineralogy, densities, seismic wave speeds, and H2O contents
(2005)
The literature on the costing of mitigation measures for reducing impacts of natural hazards is rather fragmented. This paper provides a concise overview of the current state of knowledge in Europe on the costing of mitigation measures for the reduction of natural hazard risks (droughts, floods, storms and induced coastal hazards as well as alpine hazards) and identifies knowledge gaps and related research recommendations. Furthermore, it provides a taxonomy of related mitigation options, classifying them into nine categories: (1) management plans, land-use planning, and climate adaptation; (2) hazard modification; (3) infrastructure; (4) mitigation measures (stricto sensu); (5) communication in advance of events; (6) monitoring and early warning systems; (7) emergency response and evacuation; (8) financial incentives; and (9) risk transfer (including insurance). It is found that the costing of mitigation measures in European and in other countries has almost exclusively focused on estimating direct costs. A cost assessment framework that addresses a range of costs, possibly informed by multiple stakeholders, would provide more accurate estimates and could provide better guidance to decision makers. (C) 2014 American Society of Civil Engineers.
The aim of the study is to record the occurrence of sediment deformation structures in one of the tectonically most active areas on the globe, the Tien Shan range in Central Asia and to examine the significance of the deformations as indicators of palaeoseismicity. Soft-sediment deformation structures in form of balls and pseudo-nodules are exposed in the Issyk-Kul basin, within interfingering beds of shallow lacustrine, beach and fluviatile origin. Additional deformation structures that were encountered are: a complex and chaotic folded structure, giant balls and a "pillar" structure which has not been previously reported, where marl intrudes down into coarse pebbley sand and forms pillar morphology. Liquefaction features and bedforms related to storm and breaking waves were not encountered. Neither was there evidence of turbidites. Seven field criteria for relating soft-sediment deformation to palaeoseismic triggering provide strong evidence for a seismic origin of the deformation structures. Empirical relationships between magnitude and the maximum distance from an epicenter to liquefaction sites make the active epicentral zone north of Lake Issyk- Kul, with its frequent high magnitude events, the most favorable source for the deformation structures. Luminescence dating of the sediments gives a time window of 26 +/- 2.1 to 10.5 +/- 0.7 ka BP, indicating latest Pleistocene seismic activity.(C) 2003 Elsevier B.V. All rights reserved
The Tien Shan is a most active intracontinental mountain-building range with abundant Quaternary fault-related folding. In order to improve our understanding of Quaternary intermontane basin deformation, we investigated the intermontane Issyk-Kul Lake area, an anticline that was up-warped through the piedmont cover, causing partitioning of the alluvial fan veneer. To follow the morphological scenario during the warping process, we relied on surface-exposed and trenched structures and on alluvial fans and bajadas as reference surfaces. We used air photos and satellite images to analyze the spatial -temporal morphological record and determined the age of near surface sediments by luminescence dating. We demonstrate that the up-warped Ak-Teke hills are a thrust-generated subdued anticline with strong morphological asymmetry which results from the coupling of the competing processes of up-warp and erosional feedback. The active creeks across the up-warped anticline indicate that the antecedent drainage system kept pace with the rate of uplift. The rivers which once sourced the piedmont, like the Toru-Aygyr, Kultor and the Dyuresu, became deeply entrenched and gradually transformed the study area into an abandoned morphological surface. The up-warp caused local lateral drainage diversion in front of the northern backlimb and triggered the formation of a dendritic drainage pattern upfan. Luminescence dating suggest that the period of up-warp and antecedent entrenchment started after 157 ka. The morphologically mature study area demonstrates the response of fluvial systems to growing folds on piedmont areas, induced by a propagating frontal fold at a thrust belt edge, following shortening. (C) 2004 Elsevier B.V All rights reserved
While the Dead Sea basin has been studied for a long time, the available knowledge about the detailed seismicity distribution in the area, as well as the deeper structure of the basin, is limited. Therefore, within the framework of the international project DESIRE (DEad Sea Integrated REsearch project), a dense temporary local seismological network was operated in the southern Dead Sea area. We use 530 local earthquakes, having all together 26 730 P- and S-arrival times for a simultaneous inversion of 1-D velocity models, station corrections and precise earthquake locations. Jackknife tests suggest an accuracy of the derived hypocentre locations of about 1 km. Thus, the result is the first clear image of the absolute distribution of the microseismicity of the area, especially in depth. The seismicity is concentrated in the upper crust down to 20 km depth while the lower limit of the seismicity is reached at 31 km depth. The seismic events at the eastern boundary fault (EBF) in the southern part of the study area represent the northward transform motion of the Arabian Plate along the Dead Sea Transform. North of the Boqeq fault the seismic activity represents the transfer of the motion in the pull-apart basin from the eastern to the western boundary. We find that from the surface downward the seismic events are tracing the boundary faults of the basin. The western boundary is mapped down to 12 km depth while the EBF reaches about 17 km depth, forming an asymmetric basin. One fifth of the data set is related to a specific cluster in time and space, which occurred in 2007 February at the western border fault. This cluster is aligned vertically, that is, it is perpendicular to the direction of the dominating left-lateral strike-slip movement at the main transform fault.
Within the framework of the international DESIRE (DEad Sea Integrated REsearch) project, a dense temporary local seismological network was operated in the southern Dead Sea area. During 18 recording months, 648 events were detected. Based on an already published tomography study clustering, focal mechanisms, statistics and the distribution of the microseismicity in relation to the velocity models from the tomography are analysed. The determined b value of 0.74 leads to a relatively high risk of large earthquakes compared to the moderate microseismic activity. The distribution of the seismicity indicates an asymmetric basin with a vertical strike-slip fault forming the eastern boundary of the basin, and an inclined western boundary, made up of strike-slip and normal faults. Furthermore, significant differences between the area north and south of the Bokek fault were observed. South of the Bokek fault, the western boundary is inactive while the entire seismicity occurs on the eastern boundary and below the basin-fill sediments. The largest events occurred here, and their focal mechanisms represent the northwards transform motion of the Arabian plate along the Dead Sea Transform. The vertical extension of the spatial and temporal cluster from February 2007 is interpreted as being related to the locking of the region around the Bokek fault. North of the Bokek fault similar seismic activity occurs on both boundaries most notably within the basin-fill sediments, displaying mainly small events with strike-slip mechanism and normal faulting in EW direction. Therefore, we suggest that the Bokek fault forms the border between the single transform fault and the pull-apart basin with two active border faults.
Cities often have higher species diversity than the surrounding landscape. This diversity is important for both nature conservation and urban planning. The recreation of residents and the protection of species and habitats are simultaneous targets of maintaining urban green spaces. Data about the distribution and richness of species and their habitats have been compiled frequently; however, it is difficult and costly to measure the complete biodiversity of a region, necessitating useful surrogates. We tested species and habitat data in 27 protected areas in a Central German city and asked (1) whether the diversity of selected taxa acts as a surrogate for the diversity of other taxa and total investigated diversity, and (2) whether landscape structure and human impact explain species richness. Landscape structure metrics were based on soil and habitat types; human influence was measured as the degree of hemeroby. We tested and accounted for sample bias prior to analyses. (1) Vascular plant species richness explained total richness and single taxon richness best. (2) The size of a protected area was the most important predictor of species richness. After correcting for the effect of size, shape complexity, isolation, and matrix properties remained significant. Accordingly, the type of data frequently used for urban planning - collected over several years, by various persons for various purposes - is suitable regarding systematic conservation planning for species richness. The surrogate taxa concept applies in urban areas but with restrictions. Additionally, species richness should be examined in the context of both the city and its surrounding countryside.
Pollination is a key ecosystem service. Pollinators, however, are in decline and their service is increasingly threatened. The decline is driven by several factors, most of which are related to agricultural management. However, the complexity of the landscape system, consisting of both cropped and non cropped areas, makes it difficult to address or even quantify the role of farming practices in pollinator abundance. Therefore, the aim of this paper is to improve our understanding of the relationships between pollinator abundance and their habitat use. We intend to identify and quantify the driving environmental factors that determine pollinator abundance in agricultural landscapes on the crop and landscape scale. These information helps us to design algorithms that can be used as a basis for predicting pollinator abundance on agricultural fields. To integrate varying environmental conditions data sampling was performed on farms in three different regions in Germany. Pollinators were classified into different groups with three aggregation levels. We observed crop parameters as well as landscape parameters in the areas surrounding fields in addition to temporal aspects. Generalized linear models (GLMs) were then calculated. Our results showed that both crop and landscape parameters affect pollinator abundance on agricultural fields. However, the explanatory power of the included parameters varied strongly among the particular pollinator groups and between aggregation levels. Furthermore, differentiation between species groups improves the explanatory power compared to models that are more aggregated. We also found that the temporal match between main activity periods of the particular pollinator groups and resource supply by the crop species is a key factor when analysing pollinator abundance. In conclusion, we demonstrated that the assessment and support of pollination services should be carried out with regard to individual pollinator groups. When studying pollinator abundance, the crop as well as the landscape scale should be addressed. A range of different habitat requirements and different activity periods of the pollinator groups must be covered to maintain pollination services, and therefore both diverse landscapes and diverse crop rotations are of crucial importance. (C) 2017 Elsevier B.V. All rights reserved.
Here I present a comparison between two of the most widely used reduced-complexity models for the representation of sediment transport and deposition processes, namely the transport-limited (or TL) model and the under-capacity (or xi-q) model more recently developed by Davy and Lague (2009). Using both models, I investigate the behavior of a sedimentary continental system of length L fed by a fixed sedimentary flux from a catchment of size A(0) in a nearby active orogen through which sediments transit to a fixed base level representing a large river, a lake or an ocean. This comparison shows that the two models share the same steady-state solution, for which I derive a simple 1D analytical expression that reproduces the major features of such sedimentary systems: a steep fan that connects to a shallower alluvial plain. The resulting fan geometry obeys basic observational constraints on fan size and slope with respect to the upstream drainage area, A(0). The solution is strongly dependent on the size of the system, L, in comparison to a distance L-0, which is determined by the size of A(0), and gives rise to two fundamentally different types of sedimentary systems: a constrained system where L < L-0 and open systems where L > L-0. I derive simple expressions that show the dependence of the system response time on the system characteristics, such as its length, the size of the upstream catchment area, the amplitude of the incoming sedimentary flux and the respective rate parameters (diffusivity or erodibility) for each of the two models. I show that the xi-q model predicts longer response times. I demonstrate that although the manner in which signals propagates through the sedimentary system differs greatly between the two models, they both predict that perturbations that last longer than the response time of the system can be recorded in the stratigraphy of the sedimentary system and in particular of the fan. Interestingly, the xi-q model predicts that all perturbations in the incoming sedimentary flux will be transmitted through the system, whereas the TL model predicts that rapid perturbations cannot. I finally discuss why and under which conditions these differences are important and propose observational ways to determine which of the two models is most appropriate to represent natural systems.
One of the main purposes of detrital thermochronology is to provide constraints on the regional-scale exhumation rate and its spatial variability in actively eroding mountain ranges. Procedures that use cooling age distributions coupled with hypsometry and thermal models have been developed in order to extract quantitative estimates of erosion rate and its spatial distribution, assuming steady state between tectonic uplift and erosion. This hypothesis precludes the use of these procedures to assess the likely transient response of mountain belts to changes in tectonic or climatic forcing. Other methods are based on an a priori knowledge of the in situ distribution of ages to interpret the detrital age distributions. In this paper, we describe a simple method that, using the observed detrital mineral age distributions collected along a river, allows us to extract information about the relative distribution of erosion rates in an eroding catchment without relying on a steady-state assumption, the value of thermal parameters or an a priori knowledge of in situ age distributions. The model is based on a relatively low number of parameters describing lithological variability among the various sub-catchments and their sizes and only uses the raw ages. The method we propose is tested against synthetic age distributions to demonstrate its accuracy and the optimum conditions for it use. In order to illustrate the method, we invert age distributions collected along the main trunk of the Tsangpo-Siang-Brahmaputra river system in the eastern Himalaya. From the inversion of the cooling age distributions we predict present-day erosion rates of the catchments along the Tsangpo-Siang-Brahmaputra river system, as well as some of its tributaries. We show that detrital age distributions contain dual information about present-day erosion rate, i. e., from the predicted distribution of surface ages within each catchment and from the relative contribution of any given catchment to the river distribution. The method additionally allows comparing modern erosion rates to long-term exhumation rates. We provide a simple implementation of the method in Python code within a Jupyter Notebook that includes the data used in this paper for illustration purposes.
To calibrate delta O-18 time-series from speleothems in the eastern Indian summer monsoon (ISM) region of India, and to understand the moisture regime over the northern Bay of Bengal (BoB) we analyze the delta O-18 and delta D of rainwater, collected in 2007 and 2008 near Cherrapunji, India. delta D values range from + 18.5 parts per thousand to 144.4 parts per thousand, while delta O-18 varies between +0.8 parts per thousand and 18.8 parts per thousand. The Local Meteoric Water Line (LMWL) is found to be indistinguishable from the Global Meteoric Water Line (GMWL). Late ISM (September-October) rainfall exhibits lowest delta O-18 and delta D values, with little relationship to the local precipitation amount. There is a trend to lighter isotope values over the course of the ISM, but it does not correlate with the patterns of temperature and rainfall amount delta O-18 and delta D time-series have to be interpreted with caution in terms of the 'amount effect' in this subtropical region. We find that the temporal trend in delta O-18 reflects increasing transport distance during the ISM, isotopic changes in the northern BoB surface waters during late ISM, and vapor re-equilibration with rain droplets. Using an isotope box model for surface ocean waters, we quantify the potential influence of river runoff on the isotopic composition of the seasonal freshwater plume in the northern BoB. Temporal variations in this source can contribute up to 25% of the observed changes in stable isotopes of precipitation in NE India. To delineate other moisture sources, we use backward trajectory computations and find a strong correlation between source region and isotopic composition. Palaeoclimatic stable isotope time-series from northeast Indian speleothems likely reflect changes in moisture source and transport pathway, as well as the isotopic composition of the BoB surface water, all of which in turn reflect ISM strength. Stalagmite records from the region can therefore be interpreted as integrated measures of the ISM strength.
Reliable age models are fundamental for any palaeoclimate reconstruction. Available interpolation procedures between age control points are often inadequately reported, and very few translate age uncertainties to proxy uncertainties. Most available modeling algorithms do not allow incorporation of layer counted intervals to improve the confidence limits of the age model in question.
We present a framework that allows detection and interactive handling of age reversals and hiatuses, depth-age modeling, and proxy-record reconstruction. Monte Carlo simulation and a translation procedure are used to assign a precise time scale to climate proxies and to translate dating uncertainties to uncertainties in the proxy values. The presented framework allows integration of incremental relative dating information to improve the final age model. The free software package COPRA1.0 facilitates easy interactive usage.
The strongly peraluminous and P-rich, protolithionite and zinnwaldite leucogranites from Podlesi, western Krusne Hory Mts., Czech Republic, contain accessory zircon with extraordinary enrichment of several elements, which constitute trace elements in common zircon. Elements showing a not yet reported anomalous enrichment include P (up to 20.2 wt.% P2O5; equivalent to 0.60 apfu, formula calculated on the basis of 4 oxygen atoms), Bi (up to 9.0 wt.% Bi2O3; 0.086 apfu), Nb (up to 6.7 wt.% Nb2O5, 0.12 apfu), Sc (up to 3.45 wt.% Sc2O3; 0.10 apfu), U (up to 14.8 wt.% UO2; 0.12 apfu) and F (up to 3.81 wt.% F; 0.42 apfu). Strong enrichment of P preferentially involved the berlinite-type substitution (2 Si4+ double left right arrow P5+ + Al3+) implying that significant Al may enter the Si position in zircon. Incorporation of other exotic elements is primarily governed by the xenotime (Si4++Zr4+ double left right arrow P5++Y3+), pretulite (Sc3++P5+ double left right arrow Zr4++Si4+), brabantite-type (Ca2++(U, Th)(4+)+2P(5+) double left right arrow 2Zr(4+)+2Si(4+)), and ximengite-type (Bi3++P5+double left right arrow Zr4++Si4+) substitution reactions. One part of the anomalous zircons formed late-magmatically, from a strongly peraluminous, P-F-U-rich hydrous residual melt that gave rise to the zinnwaldite granite. Interaction with aggressive residual fluids and metamictization have further aided in element enrichment or depletion, particularly in altered parts of zircon contained in the protolithionite granite. Most of the zircon from F-rich greisens have a composition close to endmember ZrSiO4 and are chemically distinct from zircon in its granite parent. This discrepancy implies that at Podlesi, granitic zircon became unstable and completely dissolved during greisenization. Part of the mobilized elements was reprecipitated in newly grown, hydrothermal zircon.
Synergistic applications based on integrated hyperspectral and lidar data are receiving a growing interest from the remote-sensing community. A prerequisite for the optimum sensor fusion of hyperspectral and lidar data is an accurate geometric coalignment. The simple unadjusted integration of lidar elevation and hyperspectral reflectance causes a substantial loss of information and does not exploit the full potential of both sensors. This paper presents a novel approach for the geometric coalignment of hyperspectral and lidar airborne data, based on their respective adopted return intensity information. The complete approach incorporates ray tracing and subpixel procedures in order to overcome grid inherent discretization. It aims at the correction of extrinsic and intrinsic (camera resectioning) parameters of the hyperspectral sensor. In additional to a tie-point-based coregistration, we introduce a ray-tracing-based back projection of the lidar intensities for area-based cost aggregation. The approach consists of three processing steps. First is a coarse automatic tie-point-based boresight alignment. The second step coregisters the hyperspectral data to the lidar intensities. Third is a parametric coalignment refinement with an area-based cost aggregation. This hybrid approach of combining tie-point features and area-based cost aggregation methods for the parametric coregistration of hyperspectral intensity values to their corresponding lidar intensities results in a root-mean-square error of 1/3 pixel. It indicates that a highly integrated and stringent combination of different coalignment methods leads to an improvement of the multisensor coregistration.
The fusion of hyperspectral imaging (HSI) sensor and airborne lidar scanner (ALS) data provides promising potential for applications in environmental sciences. Standard fusion approaches use reflectance information from the HSI and distance measurements from the ALS to increase data dimen-sionality and geometric accuracy. However, the potential for data fusion based on the respective intensity information of the complementary active and passive sensor systems is high and not yet fully exploited. Here, an approach for the rigorous illumination correction of HSI data, based on the radiometric cross-calibrated return intensity information of ALS data, is presented. The cross calibration utilizes a ray tracing-based fusion of both sensor measurements by intersecting their particular beam shapes. The developed method is capable of compensating for the drawbacks of passive HSI systems, such as cast and cloud shadowing effects, illumination changes over time, across track illumination, and partly anisotropy effects. During processing, spatial and temporal differences in illumination patterns are detected and corrected over the entire HSI wavelength domain. The improvement in the classification accuracy of urban and vegetation surfaces demonstrates the benefit and potential of the proposed HSI illumination correction. The presented approach is the first step toward the rigorous in-flight fusion of passive and active system characteristics, enabling new capabilities for a variety of applications.
This paper introduces the project on 'Assessing the impact of land use change on hydrology by ensemble modeling (LUCHEM)' that aims at investigating the envelope of predictions on changes in hydrological fluxes due to land use change. As part of a series of four papers, this paper outlines the motivation and setup of LUCHEM, and presents a model intercomparison for the present-day simulation results. Such an intercomparison provides a valuable basis to investigate the effects of different model structures on model predictions and paves the ground for the analysis of the performance of multi-model ensembles and the reliability of the scenario predictions in companion papers. in this study, we applied a set of 10 lumped, semi-lumped and fully distributed hydrological models that have been previously used in land use change studies to the low mountainous Dill catchment. Germany. Substantial differences in model performance were observed with Nash-Sutcliffe efficiencies ranging from 0.53 to 0.92. Differences in model performance were attributed to (1) model input data, (2) model calibration and (3) the physical basis of the models. The models were applied with two sets of input data: an original and a homogenized data set. This homogenization of precipitation, temperature and leaf area index was performed to reduce the variation between the models. Homogenization improved the comparability of model simulations and resulted in a reduced average bias, although some variation in model data input remained. The effect of the physical differences between models on the long-term water balance was mainly attributed to differences in how models represent evapotranspiration. Semi-lumped and lumped conceptual models slightly outperformed the fully distributed and physically based models. This was attributed to the automatic model calibration typically used for this type of models. Overall, however, we conclude that there was no superior model if several measures of model performance are considered and that all models are suitable to participate in further multi-model ensemble set-ups and land use change scenario investigations.
Forest structure is a crucial component in the assessment of whether a forest is likely to act as a carbon sink under changing climate. Detailed 3D structural information about the tundra–taiga ecotone of Siberia is mostly missing and still underrepresented in current research due to the remoteness and restricted accessibility. Field based, high-resolution remote sensing can provide important knowledge for the understanding of vegetation properties and dynamics. In this study, we test the applicability of consumer-grade Unmanned Aerial Vehicles (UAVs) for rapid calculation of stand metrics in treeline forests. We reconstructed high-resolution photogrammetric point clouds and derived canopy height models for 10 study sites from NE Chukotka and SW Yakutia. Subsequently, we detected individual tree tops using a variable-window size local maximum filter and applied a marker-controlled watershed segmentation for the delineation of tree crowns. With this, we successfully detected 67.1% of the validation individuals. Simple linear regressions of observed and detected metrics show a better correlation (R2) and lower relative root mean square percentage error (RMSE%) for tree heights (mean R2 = 0.77, mean RMSE% = 18.46%) than for crown diameters (mean R2 = 0.46, mean RMSE% = 24.9%). The comparison between detected and observed tree height distributions revealed that our tree detection method was unable to representatively identify trees <2 m. Our results show that plot sizes for vegetation surveys in the tundra–taiga ecotone should be adapted to the forest structure and have a radius of >15–20 m to capture homogeneous and representative forest stands. Additionally, we identify sources of omission and commission errors and give recommendations for their mitigation. In summary, the efficiency of the used method depends on the complexity of the forest’s stand structure.
Compound natural hazards likeEl Ninoevents cause high damage to society, which to manage requires reliable risk assessments. Damage modelling is a prerequisite for quantitative risk estimations, yet many procedures still rely on expert knowledge, and empirical studies investigating damage from compound natural hazards hardly exist. A nationwide building survey in Peru after theEl Ninoevent 2017 - which caused intense rainfall, ponding water, flash floods and landslides - enables us to apply data-mining methods for statistical groundwork, using explanatory features generated from remote sensing products and open data. We separate regions of different dominant characteristics through unsupervised clustering, and investigate feature importance rankings for classifying damage via supervised machine learning. Besides the expected effect of precipitation, the classification algorithms select the topographic wetness index as most important feature, especially in low elevation areas. The slope length and steepness factor ranks high for mountains and canyons. Partial dependence plots further hint at amplified vulnerability in rural areas. An example of an empirical damage probability map, developed with a random forest model, is provided to demonstrate the technical feasibility.
The flash-flood in Braunsbach in the north-eastern part of Baden-Wuerttemberg/Germany was a particularly strong and concise event which took place during the floods in southern Germany at the end of May/early June 2016. This article presents a detailed analysis of the hydro-meteorological forcing and the hydrological consequences of this event. A specific approach, the "forensic hydrological analysis" was followed in order to include and combine retrospectively a variety of data from different disciplines. Such an approach investigates the origins, mechanisms and course of such natural events if possible in a "near real time" mode, in order to follow the most recent traces of the event. The results show that it was a very rare rainfall event with extreme intensities which, in combination with catchment properties, led to extreme runoff plus severe geomorphological hazards, i.e. great debris flows, which together resulted in immense damage in this small rural town Braunsbach. It was definitely a record-breaking event and greatly exceeded existing design guidelines for extreme flood discharge for this region, i.e. by a factor of about 10. Being such a rare or even unique event, it is not reliably feasible to put it into a crisp probabilistic context. However, one can conclude that a return period clearly above 100 years can be assigned for all event components: rainfall, peak discharge and sediment transport. Due to the complex and interacting processes, no single flood cause or reason for the very high damage can be identified, since only the interplay and the cascading characteristics of those led to such an event. The roles of different human activities on the origin and/or intensification of such an extreme event are finally discussed. (C) 2018 Elsevier B.V. All rights reserved.
Flood generation in mountainous headwater catchments is governed by rainfall intensities, by the spatial distribution of rainfall and by the state of the catchment prior to the rainfall, e. g. by the spatial pattern of the soil moisture, groundwater conditions and possibly snow. The work presented here explores the limits and potentials of measuring soil moisture with different methods and in different scales and their potential use for flood simulation. These measurements were obtained in 2007 and 2008 within a comprehensive multi-scale experiment in the Weisseritz headwater catchment in the Ore-Mountains, Germany. The following technologies have been applied jointly thermogravimetric method, frequency domain reflectometry (FDR) sensors, spatial time domain reflectometry (STDR) cluster, ground-penetrating radar (GPR), airborne polarimetric synthetic aperture radar (polarimetric SAR) and advanced synthetic aperture radar (ASAR) based on the satellite Envisat. We present exemplary soil measurement results, with spatial scales ranging from point scale, via hillslope and field scale, to the catchment scale. Only the spatial TDR cluster was able to record continuous data. The other methods are limited to the date of over-flights (airplane and satellite) or measurement campaigns on the ground. For possible use in flood simulation, the observation of soil moisture at multiple scales has to be combined with suitable hydrological modelling, using the hydrological model WaSiM-ETH. Therefore, several simulation experiments have been conducted in order to test both the usability of the recorded soil moisture data and the suitability of a distributed hydrological model to make use of this information. The measurement results show that airborne-based and satellite-based systems in particular provide information on the near-surface spatial distribution. However, there are still a variety of limitations, such as the need for parallel ground measurements (Envisat ASAR), uncertainties in polarimetric decomposition techniques (polarimetric SAR), very limited information from remote sensing methods about vegetated surfaces and the non-availability of continuous measurements. The model experiments showed the importance of soil moisture as an initial condition for physically based flood modelling. However, the observed moisture data reflect the surface or near-surface soil moisture only. Hence, only saturated overland flow might be related to these data. Other flood generation processes influenced by catchment wetness in the subsurface such as subsurface storm flow or quick groundwater drainage cannot be assessed by these data. One has to acknowledge that, in spite of innovative measuring techniques on all spatial scales, soil moisture data for entire vegetated catchments are still today not operationally available. Therefore, observations of soil moisture should primarily be used to improve the quality of continuous, distributed hydrological catchment models that simulate the spatial distribution of moisture internally. Thus, when and where soil moisture data are available, they should be compared with their simulated equivalents in order to improve the parameter estimates and possibly the structure of the hydrological model.
To support scientifically sound water management in dryland environments a modelling system has been developed for the quantitative assessment of water and sediment fluxes in catchments, transport in the river system, and retention in reservoirs. The spatial scale of interest is the mesoscale because this is the scale most relevant for management of water and land resources.
This modelling system comprises process-oriented hydrological components tailored for dryland characteristics coupled with components comprising hillslope erosion, sediment transport and reservoir deposition processes. The spatial discretization is hierarchically designed according to a multi-scale concept to account for particular relevant process scales. The non-linear and partly intermittent run-off generation and sediment dynamics are dealt with by accounting for connectivity phenomena at the intersections of landscape compartments. The modelling system has been developed by means of data from nested research catchments in NE-Spain and in NE-Brazil.
In the semi-arid NE of Brazil sediment retention along the topography is the main process for sediment retention at all scales, i.e. the sediment delivery is transport limited. This kind of deposition retains roughly 50 to 60 % of eroded sediment, maintaining a similar deposition proportion in all spatial scales investigated. On the other hand, the sediment retained in reservoirs is clearly related to the scale, increasing with catchment area. With increasing area, there are more reservoirs, increasing the possibility of deposition. Furthermore, the area increase also promotes an increase in flow volume, favouring the construction of larger reservoirs, which generally overflow less frequently and retain higher sediment fractions. The second example comprises a highly dynamic Mediterranean catchment in NE-Spain with nested sub-catchments and reveals the full dynamics of hydrological, erosion and deposition features. The run-off modelling performed well with only some overestimation during low-flow periods due to the neglect of water losses along the river. The simulated peaks in sediment flux are reproduced well, while low-flow sediment transport is less well captured, due to the disregard of sediment remobilization in the riverbed during low flow.
This combined observation and modelling study deepened the understanding of hydro-sedimentological systems characterized by flashy run-off generation and by erosion and sediment transport pulses through the different landscape compartments. The connectivity between the different landscape compartments plays a very relevant role, regarding both the total mass of water and sediment transport and the transport time through the catchment.
The hydrological cycle is a dynamic system by its nature, but sometimes accelerated through anthropogenic activity. A "hydrological change" (i.e. a water cycle that is significantly changing over a longer period of time) can be very different in character, depending on the specific natural conditions and the underlying spatial and temporal scales. Such changes may affect the availability and quality of water as essential pre-requisites for human development and ecosystem stability. Hydrological extremes, such as floods and droughts, may also be affected, what is also vitally important, because of their profound economic and societal impacts. Anthropogenically induced hydrological change can be attributed to three main external causes: first, the Earth's climate is changing significantly and thus directly affecting the terrestrial hydro-systems via the exchange of energy and heat. The second major issue is the land cover and its management that has been modified fundamentally by conversion of land for agriculture, forestry, and other purposes such as industrialisation and urbanisation. Finally, water resources are being used more than ever for human development, especially for agriculture, industrial activities, and navigation. If the regional terrestrial hydrological cycle is changing and counter-measures are desirable, it is from a scientific perspective mandatory to understand the extent and nature of such changes, and, especially, to identify their possible anthropogenic origin. There are, however, fundamental gaps in our knowledge, in particular about the role of feedbacks between individual processes and compartments of the hydrological cycle or the relevance of the interactions with other sub-systems of our planet, such as the atmosphere or the vegetation. This paper mentions several examples of hydrological change and discusses their identification, interaction processes, and feedback mechanisms, along with modelling issues. The possibilities and limitations of modelling are demonstrated by means of two studies: one from the river-lake system on the Middle-Havel River and one from the catchment of the Wahnbach Reservoir. The applied model systems comprise a series of consecutively coupled individual models (so-called one-way-coupling). Model systems that are able reflect feedback effects (two-way- coupling) are still in the development stage. It became clear that the applied model systems were able to reproduce the observed dynamics of the hydrological cycle and of selected matter fluxes. However, one has to be aware that the simulated time periods and scenarios represent rather moderately transient conditions, what is the justification why the one-way-coupling seems to be applicable. Furthermore, it was shown that the modelling uncertainty is considerably large. Nevertheless, this uncertainty can be distinguished from effects of changed internal systems dynamics or from changed boundary conditions, what is a basis for the usability of such model systems for prognostic purposes.
How information about sediment transport processes is transmitted to the sedimentary record remains a complex problem for the interpretation of fluvial stratigraphy. Alluvial fan deposits represent the condensed archive of sediment transport, which is at least partly controlled by tectonics and climate. For three coupled catchment-fan systems in northern Death Valley, California, we measure grain size across 12 well-preserved Holocene and late-Pleistocene surfaces, mapped in detail from field observations and remote sensing. Our results show that fan surfaces correlated to the late Pleistocene are, on average, 30-50% coarser than active or Holocene fan surfaces. We adopt a self-similar form of grain size distribution based on the observed stability of the ratio between mean grain size and standard deviation downstream. Using statistical analysis, we show that fan surface grain size distributions are self-similar. We derive a relative mobility function using our self-similar grain size distributions, which describes the relative probability of a given grain size being transported. We show that the largest mobile grain sizes are between 20 and 35mm, a value that varies over time and is clearly lower in the Holocene than in the Pleistocene; a change we suggest is due to a drier climate in the Holocene. These results support recent findings that alluvial fan sedimentology can record past environmental change and that these landscapes are potentially sensitive to climatic change over a glacial-interglacial cycle. We demonstrate that the self-similarity methodology offers a means to explore changes in relative mobility of grain sizes from preserved fluvial deposits. Plain Language Summary A key challenge in Earth Science is understanding how landscapes respond to climate. It may be possible to observe measurable differences in certain landscapes settings such as alluvial fans in desert regions. Alluvial fans are believed to be effective recorders of climate, representing a cumulative store of material transported downstream by rainfall-sensitive river systems. In northern Death Valley, California, we measure at high resolution grain size on three alluvial fans with surfaces that date from the Holocene and the arid climate of today to the 20-40% wetter late-Pleistocene epoch. We find that older late-Pleistocene surfaces are coarser on average than surfaces deposited during the modern and Holocene dry period, suggesting a changing sediment transport regime potentially in response to precipitation. We also show that measured grain size distributions within and between surfaces can be successfully normalized based on the decay in mean grain size and variance downstream, exhibiting a self-similar pattern. Finally, we employ a grain size relative mobility model using our field data to establish which grain sizes are likely to be in transport or locked in the substrate. This model predicts that during the wetter late-Pleistocene mobile grain sizes are up to 40% larger than during the Holocene.
Knowledge of the origin of suspended sediment is important for improving our understanding of sediment dynamics and thereupon support of sustainable watershed management. An direct approach to trace the origin of sediments is the fingerprinting technique. It is based on the assumption that potential sediment sources can be discriminated and that the contribution of these sources to the sediment can be determined on the basis of distinctive characteristics (fingerprints). Recent studies indicate that visible-near-infrared (VNIR) and shortwave-infrared (SWIR) reflectance characteristics of soil may be a rapid, inexpensive alternative to traditional fingerprint properties (e.g. geochemistry or mineral magnetism).
To further explore the applicability of VNIR-SWIR spectral data for sediment tracing purposes, source samples were collected in the Isabena watershed, a 445 km(2) dryland catchment in the central Spanish Pyrenees. Grab samples of the upper soil layer were collected from the main potential sediment source types along with in situ reflectance spectra. Samples were dried and sieved, and artificial mixtures of known proportions were produced for algorithm validation. Then, spectral readings of potential source and artificial mixture samples were taken in the laboratory. Colour coefficients and physically based parameters were calculated from in situ and laboratory-measured spectra. All parameters passing a number of prerequisite tests were subsequently applied in discriminant function analysis for source discrimination and mixing model analyses for source contribution assessment.
The three source types (i.e. badlands, forest/grassland and an aggregation of other sources, including agricultural land, shrubland, unpaved roads and open slopes) could be reliably identified based on spectral parameters. Laboratory-measured spectral fingerprints permitted the quantification of source contribution to artificial mixtures, and introduction of source heterogeneity into the mixing model decreased accuracies for some source types. Aggregation of source types that could not be discriminated did not improve mixing model results. Despite providing similar discrimination accuracies as laboratory source parameters, in situ derived source information was found to be insufficient for contribution modelling.
The laboratory mixture experiment provides valuable insights into the capabilities and limitations of spectral fingerprint properties. From this study, we conclude that combinations of spectral properties can be used for mixing model analyses of a restricted number of source groups, whereas more straightforward in situ measured source parameters do not seem suitable. However, modelling results based on laboratory parameters also need to be interpreted with care and should not rely on the estimates of mean values only but should consider uncertainty intervals as well.
Knowledge of sediment sources is a prerequisite for sustainable management practices and may furthermore improve our understanding of water and sediment fluxes. Investigations have shown that a number of characteristic soil properties can be used as "fingerprints" to trace back the sources of river sediments. Spectral properties have recently been successfully used as such characteristics in fingerprinting studies. Despite being less labour-intensive than geochemical analyses, for example, spectroscopy allows measurements of small amounts of sediment material (> 60 mg), thus enabling inexpensive analyses even of intra-event variability. The focus of this study is on the examination of spectral properties of fluvial sediment samples to detect changes in source contributions, both between and within individual flood events.
Sediment samples from the following three different origins were collected in the Isabena catchment (445 km(2)) in the central Spanish Pyrenees: (1) soil samples from the main potential source areas, (2) stored fine sediment from the channel bed once each season in 2011 and (3) suspended sediment samples during four flood events in autumn 2011 and spring 2012 at the catchment outlet as well as at several subcatchment outlets. All samples were dried and measured for spectral properties in the laboratory using an ASD spectroradiometer. Colour parameters and physically based features (e.g. organic carbon, iron oxide and clay content) were calculated from the spectra. Principal component analyses (PCA) were applied to all three types of samples to determine natural clustering of samples, and a mixing model was applied to determine source contributions.
We found that fine sediment stored in the river bed seems to be mainly influenced by grain size and seasonal variability, while sampling location-and thus the effect of individual tributaries or subcatchments-seem to be of minor importance. Suspended sediment sources were found to vary between, as well as within, flood events; although badlands were always the major source. Forests and grasslands contributed little (< 10 %), and other sources (not further determinable) contributed up to 40 %. The analyses further suggested that sediment sources differ among the subcatchments and that subcatchments comprising relatively large proportions of badlands contributed most to the four flood events analyzed.
Spectral fingerprints provide a rapid and cost-efficient alternative to conventional fingerprint properties. However, a combination of spectral and conventional fingerprint properties could potentially permit discrimination of a larger number of source types.
Background: GEOMAGIA50.v3 for sediments is a comprehensive online database providing access to published paleomagnetic, rock magnetic, and chronological data obtained from lake and marine sediments deposited over the past 50 ka. Its objective is to catalogue data that will improve our understanding of changes in the geomagnetic field, physical environments, and climate.
Findings: GEOMAGIA50.v3 for sediments builds upon the structure of the pre-existing GEOMAGIA50 database for magnetic data from archeological and volcanic materials. A strong emphasis has been placed on the storage of geochronological data, and it is the first magnetic archive that includes comprehensive radiocarbon age data from sediments. The database will be updated as new sediment data become available.
Conclusions: The web-based interface for the sediment database is located at http://geomagia.gfz-potsdam.de/geomagiav3/SDquery.php. This paper is a companion to Brown et al. (Earth Planets Space doi:10.1186/s40623-015-0232-0,2015) and describes the data types, structure, and functionality of the sediment database.
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.
Movements of tectonic plates often induce oblique deformation at divergent plate boundaries. This is in striking contrast with traditional conceptual models of rifting and rifted margin formation, which often assume 2-D deformation where the rift velocity is oriented perpendicular to the plate boundary. Here we quantify the validity of this assumption by analysing the kinematics of major continent-scale rift systems in a global plate tectonic reconstruction from the onset of Pangea breakup until the present day. We evaluate rift obliquity by joint examination of relative extension velocity and local rift trend using the script-based plate reconstruction software pyGPlates. Our results show that the global mean rift obliquity since 230 Ma amounts to 34 degrees with a standard deviation of 24 degrees, using the convention that the angle of obliquity is spanned by extension direction and rift trend normal. We find that more than similar to 70 % of all rift segments exceeded an obliquity of 20 degrees demonstrating that oblique rifting should be considered the rule, not the exception. In many cases, rift obliquity and extension velocity increase during rift evolution (e.g. Australia-Antarctica, Gulf of California, South Atlantic, India-Antarctica), which suggests an underlying geodynamic correlation via obliquity-dependent rift strength. Oblique rifting produces 3-D stress and strain fields that cannot be accounted for in simplified 2-D plane strain analysis. We therefore highlight the importance of 3-D approaches in modelling, surveying, and interpretation of most rift segments on Earth where oblique rifting is the dominant mode of deformation.
In piston cores from the open subarctic Pacific and the Okhotsk Sea, diatom-bound delta N-15 (delta N-15(db)), biogenic opal, calcium carbonate, and barium were measured from coretop to the previous glacial maximum (MIS 6). Glacial intervals are generally characterized by high delta N-15(db) (similar to 8 parts per thousand) and low productivity, whereas interglacial intervals have a lower delta N-15(db) (5.7-6.3 parts per thousand) and indicate high biogenic productivity. These data extend the regional swath of evidence for nearly complete surface nutrient utilization during glacial maxima, consistent with stronger upper water column stratification throughout the subarctic region during colder intervals. An early deglacial decline in delta N-15(db) of 2 parts per thousand at similar to 17.5 ka, previously observed in the Bering Sea, is found here in the open subarctic Pacific record and arguably also in the Okhotsk, and a case can be made that a similar decrease in delta N-15(db) occurred in both regions at the previous deglaciation as well. The early deglacial delta N-15(db) decrease, best explained by a decrease in surface nutrient utilization, appears synchronous with southern hemisphere-associated deglacial changes and with the Heinrich 1 event in the North Atlantic. This delta N-15(db) decrease may signal the initial deglacial weakening in subarctic North Pacific stratification and/or a deglacial increase in shallow subsurface nitrate concentration. If the former, it would be the North Pacific analogue to the increase in vertical exchange inferred for the Southern Ocean at the time of Heinrich Event 1. In either case, the lack of any clear change in paleoproductivity proxies during this interval would seem to require an early deglacial decrease in the iron-to-nitrate ratio of subsurface nutrient supply or the predominance of light limitation of phytoplankton growth during the deglaciation prior to Bolling-Allerod warming.
Water stable isotope ratios of tropical precipitation predominantly reflect moisture source and precipitation intensity. Trees can incorporate the isotopic signals into annual tree-ring cellulose records, permitting reconstruction of the temporal changes of hydroclimate over decades to millennia. This is especially valuable in the Himalayas where the understanding of monsoon dynamics is limited by the lack of a dense and representative observational network. We have analyzed tree ring delta O-18 records from two distinct physiographic sites along the upper Kali Gandaki valley in the central Nepal Himalayas, representing the wet High-Himalayas and the Trans-Himalayan dryland to the north. Empirical correlations and regression analyses were compared to an in-situ calibrated oxygen isotope fractionation model, exploring the relationships between tree ring delta O-18 and seasonal-mean variability of hydroclimatic forcing at the different locations. For this purpose, gridded precipitation data from the Asian rain gauge dataset APHRODITE, as well as high resolution onsite observations (relative humidity, air temperature, delta O-18 of precipitation and radial tree growth) were used. We found that two distinct sets of meteorological values, reflecting pre-monsoon and monsoon conditions, are needed to reproduce the measured tree ring delta O-18 values from the High-Himalayan site, but that a single set of monsoonal values performs best for the Trans-Himalayan site. We conclude that Trans-Himalayan trees capture long-term changes in strength of the Indian summer monsoon. In contrast, High-Himalayan tree ring delta(18)Orecords a more complex hydro-climatic signal reflecting both pre-monsoon and monsoon seasons with very contrasting isotopic signatures of precipitation. This difference in the two hydroclimatic proxy records offers an opportunity to reconstruct first-order hydroclimate conditions, such as local precipitation rates, and to gain new insights into monsoon timing and seasonal water source determination across the Himalayan orographic region. (C) 2019 Elsevier B.V. All rights reserved.
A local and flexible definition of the monsoon season based on hydrological evidence is important for the understanding and management of Himalayan water resources. Here, we present an objective statistical method to retrieve seasonal hydrometeorological transitions. Applied to daily rainfall data (1951-2015), this method shows an average longitudinal delay of similar to 15 days, with later monsoon onset and earlier withdrawal in the western Himalaya, consistent with the continental progression of wet air masses. This delay leads to seasons of different length along the Himalaya and biased precipitation amounts when using uniform calendric monsoon boundaries. In the Central Himalaya annual precipitation has increased, due primarily to an increase of premonsoon precipitation. These findings highlight issues associated with a static definition of monsoon boundaries and call for a deeper understanding of nonmonsoonal precipitation over the Himalayan water tower. <br /> Plain Language Summary Precipitation in the Himalayas determines water availability for the Indian foreland with large socioeconomic implications. Despite its importance, spatial and temporal patterns of precipitation are poorly understood. Here, we estimate the long-term average and trends of seasonal precipitation at the scale of individual catchments draining the Himalayas. We apply a statistical method to detect the timing of hydrometeorological seasons from local precipitation measurements, focusing on monsoon onset and withdrawal. We identify longitudinal and latitudinal delays, resulting in seasons of different length along and across the Himalayas. These spatial patterns and the annual variability of the monsoon boundaries mean that oft-used, fixed calendric dates, for example, 1 June to 30 September, may be inadequate for retrieving monsoon rainfall totals. Moreover, we find that, despite its prominent contribution to annual rainfall totals, the Indian summer monsoon cannot explain the increase of the annual precipitation over the Central Himalayas. Instead, this appears to be mostly driven by changes in premonsoon and winter rainfall. So far, little attention has been paid to premonsoon precipitation, but governed by evaporative processes and surface water availability, it may be enhanced by irrigation and changed land use in the Gangetic foreland.
Can short-term incentives induce long-lasting cooperation? Results from a public-goods experiment
(2014)
This paper investigates whether providing strong cooperation incentives only at the outset of a group interaction spills over to later periods to ensure cooperation in the long run. We study a repeated linear public-good game with punishment opportunities and a parameter change after the first ten (of twenty) rounds. Our data shows that cooperation among subjects who had experienced a higher marginal return on public-good contributions or low punishment costs in rounds 1-10 rapidly deteriorated in rounds 11-20 once these incentives were removed, eventually trending below the level of cooperation in the control group. This suggests the possibility of temporary incentives backfiring in the long run. This paper ties in with the literature highlighting the potentially adverse effects of the use of incentives. (C) 2014 Elsevier Inc. All rights reserved.
Local earthquake data from a dense temporary seismological network in the southern Dead Sea area have been analysed within the project DESIRE (Dead Sea Integrated Research Project). Local earthquakes are used for the first precise image of the distribution of the P-wave velocity and the vP/vS ratios. 65 stations registered 655 local events within 18 months of observation time. A subset of 530 well-locatable events with 26 730 P- and S-arrival times was used to calculate a tomographic model for the vP and vP/vS distribution. Since the study area is at first-order 2-D, a gradual approach was chosen, which compromised a 2-D inversion followed by a 3-D inversion. The sedimentary basin fill is clearly imaged through high vP/vS ratios and low vP. The basin fill shows an asymmetric structure with average depth of 7 km at the western boundary and depth between 10 and 14 km at the eastern boundary. This asymmetry is reflected by the vertical strike-slip eastern border fault, and the normal faulting at the western boundary, caused by the transtensional deformation within the last 5 Myr. Within the basin fill the Lisan salt diapir is imaged through low vP/vS ratios, reflecting its low fluid content. The extensions were determined to 12 km in EW and 17 km in NS direction while its depth is 56 km. The thickness of the pre-basin sediments below the basin fill cannot be derived from the tomography datait is estimated to less than 3 km from former investigations. Below the basin, down to 18 km depth very low P-wave velocities and low vP/vS ratios are observedmost likely caused by fluids from the surrounding crust or the upper mantle.
Protection motivation theory (PMT) has become a popular theory to explain the risk-reducing behavior of residents against natural hazards. PMT captures the two main cognitive processes that individuals undergo when faced with a threat, namely, threat appraisal and coping appraisal. The latter describes the evaluation of possible response measures that may reduce or avert the perceived threat. Although the coping appraisal component of PMT was found to be a better predictor of protective intentions and behavior, little is known about the factors that influence individuals’ coping appraisals of natural hazards. More insight into flood-coping appraisals of PMT, therefore, are needed to better understand the decision-making process of individuals and to develop effective risk communication strategies. This study presents the results of two surveys among more than 1,600 flood-prone households in Germany and France. Five hypotheses were tested using multivariate statistics regarding factors related to flood-coping appraisals, which were derived from the PMT framework, related literature, and the literature on social vulnerability. We found that socioeconomic characteristics alone are not sufficient to explain flood-coping appraisals. Particularly, observational learning from the social environment, such as friends and neighbors, is positively related to flood-coping appraisals. This suggests that social norms and networks play an important role in flood-preparedness decisions. Providing risk and coping information can also have a positive effect. Given the strong positive influence of the social environment on flood-coping appraisals, future research should investigate how risk communication can be enhanced by making use of the observed social norms and network effects.
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.
Receiver functions (RF) are used to investigate the upper mantle structure beneath the Eifel, the youngest volcanic area of Central Europe. Data from 96 teleseismic events recorded by 242 seismological stations from permanent and a temporary network has been analysed. The temporary network operated from 1997 November to 1998 June and covered an area of approximately 400 x 250 km(2) centred on the Eifel volcanic fields. The average Moho depth in the Eifel is approximately 30 km, thinning to ca. 28 km under the Eifel volcanic fields. RF images suggest the existence of a low velocity zone at about 60-90 km depth under the West Eifel. This observation is supported by P- and S-wave tomographic results and absorption (but the array aperture limits the resolution of the tomographic methods to the upper 400 km). There are also indications for a zone of elevated velocities at around 200 km depth, again in agreement with S-wave and absorption tomographic results. This anomaly is not visible in P-wave tomography and could be due to S-wave anisotropy. The RF anomalies at the Moho, at 60-90 km, and near 200 km depth have a lateral extent of about 100 km. The 410 km discontinuity under the Eifel is depressed by 15-25 km, which could be explained by a maximum temperature increase of +200 degrees C to +300 degrees C. In the 3-D RF image of the Eifel Plume we also notice two additional currently unexplained conversions between 410 and 550 km depth. They could represent remnants of previous subduction or anomalies due to delayed phase changes. The lateral extent of these conversions and the depression of the 410 km discontinuity is about 200 km. The 660 km discontinuity does not show any depth deviation from its expected value. Our observations are consistent with interpretation in terms of an upper mantle plume but they do not rule out connections to processes at larger depth
Understanding the evolution of continental deformation zones relies on quantifying spatial and temporal changes in deformation rates of tectonic structures. Along the eastern boundary of the Pamir-Tian Shan collision zone, we constrain secular variations of rock uplift rates for a series of five Quaternary detachment- and fault-related folds from their initiation to the modern day. When combined with GPS data, decomposition of interferometric synthetic aperture radar time series constrains the spatial pattern of surface and rock uplift on the folds deforming at decadal rates of 1-5mm/yr. These data confirm the previously proposed basinward propagation of structures during the Quaternary. By fitting our geodetic rates and previously published geologic uplift rates with piecewise linear functions, we find that gradual rate changes over >100kyr can explain the interferometric synthetic aperture radar observations where changes in average uplift rates are greater than similar to 1 mm/yr among different time intervals (similar to 10(1), 10(4-5), and 10(5-6) years).
Fluvial planation surfaces, such as straths, commonly serve as recorders of climatic and tectonic changes and are formed by the lateral erosion of rivers, a process that remains poorly understood. Here we present a study of kilometer-wide, fluvially eroded, low-relief surfaces on rapidly uplifting folds in the foreland of the southwestern Tian Shan. A combination of field work, digital elevation model analysis, and dating of fluvial deposits reveals that despite an arid climate and rapid average rock-uplift rates of 1-3mm/yr, rivers cut extensive (>1-2km wide) surfaces with typical height variations of <6m over periods of >2-6kyr. The extent of this beveling varies in space and time, such that different beveling episodes affect individual structures. Between times of planation, beveled surfaces are abandoned, incised, and deformed across the folds. In a challenge to models that link strath cutting and abandonment primarily to changes in river incision rates, we demonstrate that lateral erosion rates of antecedent streams crossing the folds have to vary by more than 1 order of magnitude to explain the creation of beveled platforms in the past and their incision at the present day. These variations do not appear to covary with climate variability and might be caused by relatively small (much less than an order of magnitude) changes in sediment or water fluxes. It remains uncertain in which settings variations in lateral bedrock erosion rates predominate over changes in vertical erosion rates. Therefore, when studying fluvial planation and strath terraces, variability of both lateral and vertical erosion rates should be considered.
The effect of lithology on the relationship between denudation rate and chemical weathering pathways
(2022)
The denudation of rocks in mountain belts exposes a range of fresh minerals to the surface of the Earth that are chemically weathered by acidic and oxygenated fluids. The impact of the resulting coupling between denudation and weathering rates fundamentally depends on the types of minerals that are weathering. Whereas silicate weathering sequesters CO2, the combination of sulfide oxidation and carbonate dissolution emits CO2 to the atmosphere. Here, we combine the concentrations of dissolved major elements in stream waters with Be-10 basin-wide denudation rates from 35 small catchments in eastern Tibet to elucidate the importance of lithology in modulating the relationships between denudation rate, chemical weathering pathways, and CO2 consumption or release. Our catchments span 3 orders of magnitude in denudation rate in low-grade flysch, high-grade metapelites, and granitoid rocks. For each stream, we estimate the concentrations of solutes sourced from silicate weathering, carbonate dissolution, and sulfide oxidation using a mixing model. We find that for all lithologies, cation concentrations from silicate weathering are largely independent of denudation rate, but solute concentrations from carbonates and, where present, sulfides increase with increasing denudation rate. With increasing denudation rates, weathering may therefore shift from consuming to releasing CO2 in both (meta)sedimentary and granitoid lithologies. For a given denudation rate, we report dissolved solid concentrations and inferred weathering fluxes in catchments underlain by (meta)sedimentary rock that are 2-10 times higher compared to catchments containing granitoid lithologies, even though climatic and topographic parameters do not vary systematically between these catchments. Thus, varying proportions of exposed (meta)sedimentary and igneous rocks during orogenesis could lead to changes in the sequestration and release of CO2 that are independent of denudation rate.
Lateral movements of alluvial river channels control the extent and reworking rates of alluvial fans, floodplains, deltas, and alluvial sections of bedrock rivers. These lateral movements can occur by gradual channel migration or by sudden changes in channel position (avulsions). Whereas models exist for rates of river avulsion, we lack a detailed understanding of the rates of lateral channel migration on the scale of a channel belt. In a two-step process, we develop here an expression for the lateral migration rate of braided channel systems in coarse, non-cohesive sediment. On the basis of photographic and topographic data from laboratory experiments of braided channels performed under constant external boundary conditions, we first explore the impact of autogenic variations of the channel-system geometry (i.e. channel-bank heights, water depths, channel-system width, and channel slope) on channel-migration rates. In agreement with theoretical expectations, we find that, under such constant boundary conditions, the laterally reworked volume of sediment is constant and lateral channel-migration rates scale inversely with the channel-bank height. Furthermore, when channel-bank heights are accounted for, lateral migration rates are independent of the remaining channel geometry parameters. These constraints allow us, in a second step, to derive two alternative expressions for lateral channel-migration rates under different boundary conditions using dimensional analysis. Fits of a compilation of laboratory experiments to these expressions suggest that, for a given channel bank-height, migration rates are strongly sensitive to water discharges and more weakly sensitive to sediment discharges. In addition, external perturbations, such as changes in sediment and water discharges or base level fall, can indirectly affect lateral channel-migration rates by modulating channel-bank heights.