Refine
Year of publication
- 2014 (225) (remove)
Document Type
- Article (195)
- Doctoral Thesis (17)
- Preprint (8)
- Review (4)
- Postprint (1)
Language
- English (225) (remove)
Keywords
- Earthquake source observations (7)
- Holocene (6)
- Tibetan Plateau (4)
- Erosion (3)
- Himalaya (3)
- Seismicity and tectonics (3)
- Aleatory variability (2)
- Andes (2)
- Baseline shift (2)
- Body waves (2)
Institute
- Institut für Geowissenschaften (225) (remove)
For a detailed characterization of near-surface environments, geophysical techniques are increasingly used to support more conventional point-based techniques such as borehole and direct-push logging. Because the underlying parameter relations are often complex, site-specific, or even poorly understood, a remaining challenging task is to link the geophysical parameter models to the actual geotechnical target parameters measured only at selected points. We propose a workflow based on nonparametric regression to establish functional relationships between jointly inverted geophysical parameters and selected geotechnical parameters as measured, for example, by different borehole and direct-push tools. To illustrate our workflow, we present field data collected to characterize a near-surface sedimentary environment Our field data base includes crosshole ground penetrating radar (GPR), seismic P-, and S-wave data sets collected between 25 m deep boreholes penetrating sand- and gravel dominated sediments. Furthermore, different typical borehole and direct-push logs are available. We perform a global joint inversion of traveltimes extracted from the crosshole geophysical data using a recently proposed approach based on particle swarm optimization. Our inversion strategy allows for generating consistent models of GPR, P-wave, and S-wave velocities including an appraisal of uncertainties. We analyze the observed complex relationships between geophysical velocities and target parameter logs using the alternating conditional expectation (ACE) algorithm. This nonparametric statistical tool allows us to perform multivariate regression analysis without assuming a specific functional relation between the variables. We are able to explain selected target parameters such as characteristic grain size values or natural gamma activity by our inverted geophysical data and to extrapolate these parameters to the inter-borehole plane covered by our crosshole experiments. We conclude that the ACE algorithm is a powerful tool to analyze a multivariate petrophysical data base and to develop an understanding of how a multi-parameter geophysical model can be linked and translated to selected geotechnical parameters.
This paper presents a Bayesian non-parametric method based on Gaussian Process (GP) regression to derive ground-motion models for peak-ground parameters and response spectral ordinates. Due to its non-parametric nature there is no need to specify any fixed functional form as in parametric regression models. A GP defines a distribution over functions, which implicitly expresses the uncertainty over the underlying data generating process. An advantage of GP regression is that it is possible to capture the whole uncertainty involved in ground-motion modeling, both in terms of aleatory variability as well as epistemic uncertainty associated with the underlying functional form and data coverage. The distribution over functions is updated in a Bayesian way by computing the posterior distribution of the GP after observing ground-motion data, which in turn can be used to make predictions. The proposed GP regression models is evaluated on a subset of the RESORCE data base for the SIGMA project. The experiments show that GP models have a better generalization error than a simple parametric regression model. A visual assessment of different scenarios demonstrates that the inferred GP models are physically plausible.
This article presents comparisons among the five ground-motion models described in other articles within this special issue, in terms of data selection criteria, characteristics of the models and predicted peak ground and response spectral accelerations. Comparisons are also made with predictions from the Next Generation Attenuation (NGA) models to which the models presented here have similarities (e.g. a common master database has been used) but also differences (e.g. some models in this issue are nonparametric). As a result of the differing data selection criteria and derivation techniques the predicted median ground motions show considerable differences (up to a factor of two for certain scenarios), particularly for magnitudes and distances close to or beyond the range of the available observations. The predicted influence of style-of-faulting shows much variation among models whereas site amplification factors are more similar, with peak amplification at around 1s. These differences are greater than those among predictions from the NGA models. The models for aleatory variability (sigma), however, are similar and suggest that ground-motion variability from this region is slightly higher than that predicted by the NGA models, based primarily on data from California and Taiwan.
The Dansgaard-Oeschger oscillations and Heinrich events described in North Atlantic sediments and Greenland ice are expressed in the climate of the tropics, for example, as documented in Arabian Sea sediments. Given the strength of this teleconnection, we seek to reconstruct its range of environmental impacts. We present geochemical and sedimentological data from core SO130-289KL from the Indus submarine slope spanning the last similar to 80 kyr. Elemental and grain size analyses consistently indicate that interstadials are characterized by an increased contribution of fluvial suspension from the Indus River. In contrast, stadials are characterized by an increased contribution of aeolian dust from the Arabian Peninsula. Decadal-scale shifts at climate transitions, such as onsets of interstadials, were coeval with changes in productivity-related proxies. Heinrich events stand out as especially dry and dusty events, indicating a dramatically weakened Indian summer monsoon, potentially increased winter monsoon circulation, and increased aridity on the Arabian Peninsula. This finding is consistent with other paleoclimate evidence for continental aridity in the northern tropics during these events. Our results strengthen the evidence that circum-North Atlantic temperature variations translate to hydrological shifts in the tropics, with major impacts on regional environmental conditions such as rainfall, river discharge, aeolian dust transport, and ocean margin anoxia.
New low-temperature thermochronological data from 80 samples in eastern Kyrgyzstan are combined with previously published data from 61 samples to constrain exhumation in a number of mountain ranges in the Central Kyrgyz Tien Shan. All sampled ranges are found to have a broadly consistent Cenozoic exhumation history, characterized by initially low cooling rates (<1 degrees C/Myr) followed by a series of increases in exhumation that occurred diachronously across the region in the late Cenozoic that are interpreted to record the onset of deformation in different mountain ranges. Combined with geological estimates for the onset of proximal deformation, our data suggest that the Central Kyrgyz Tien Shan started deforming in the late Oligocene-early Miocene, leading to the development of several, widely spaced mountain ranges separated by large intermontane basins. Subsequently, more ranges have been constructed in response to significant shortening increases across the Central Kyrgyz Tien Shan, notably in the late Miocene. The order of range construction is interpreted to reflect variations in the susceptibility of inherited structures to reactivation. Reactivated structures are also shown to have significance along strike variations in fault vergence and displacement, which have influenced the development and growth of individual mountain ranges. Moreover, the timing of deformation allows the former extent of many intermontane basins that have since been partitioned to be inferred; this can be linked to the highly time-transgressive onset of late Cenozoic coarse clastic sedimentation.
Soils in various places of the Panama Canal Watershed feature a low saturated hydraulic conductivity (K-s) at shallow depth, which promotes overland-flow generation and associated flashy catchment responses. In undisturbed forests of these areas, overland flow is concentrated in flow lines that extend the channel network and provide hydrological connectivity between hillslopes and streams. To understand the dynamics of overland-flow connectivity, as well as the impact of connectivity on catchment response, we studied an undisturbed headwater catchment by monitoring overland-flow occurrence in all flow lines and discharge, suspended sediment, and total phosphorus at the catchment outlet. We find that connectivity is strongly influenced by seasonal variation in antecedent wetness and can develop even under light rainfall conditions. Connectivity increased rapidly as rainfall frequency increased, eventually leading to full connectivity and surficial drainage of entire hillslopes. Connectivity was nonlinearly related to catchment response. However, additional information on factors such as overland-flow volume would be required to constrain relationships between connectivity, stormflow, and the export of suspended sediment and phosphorus. The effort to monitor those factors would be substantial, so we advocate applying the established links between rain event characteristics, drainage network expansion by flow lines, and catchment response for predictive modeling and catchment classification in forests of the Panama Canal Watershed and in similar regions elsewhere.
Garnet brought to the surface by late Miocene granitoids at La Galite Archipelago (Central Mediterranean, Tunisia) contains abundant primary melt and fluid inclusions. Microstructural observations and mineral chemistry define the host garnet as a peritectic phase produced by biotite incongruent melting at ~800 degrees C and 0.5GPa, under fluid-present conditions. The trapped melt is leucogranitic with an unexpected metaluminous and almost peralkaline character. Fluid inclusions are one phase at room temperature, and contain a CO2-dominated fluid, with minor H2O, N-2 and CH4. Siderite and an OH-bearing phase were identified by Raman and IR spectroscopy within every analysed inclusion, and are interpreted as products of a post-entrapment carbonation/hydration reaction between the fluid and the host during cooling. The fluid present during anatexis is therefore inferred to have been originally richer in both H2O and CO2. The production of anatectic melt with a metaluminous signature can be explained as the result of partial melting of relatively Al-poor protoliths assisted by CO2-rich fluids.
Erosion in the Himalaya is responsible for one of the greatest mass redistributions on Earth and has fueled models of feedback loops between climate and tectonics. Although the general trends of erosion across the Himalaya are reasonably well known, the relative importance of factors controlling erosion is less well constrained. Here we present 25 Be-10-derived catchment-averaged erosion rates from the Yamuna catchment in the Garhwal Himalaya, northern India. Tributary erosion rates range between similar to 0.1 and 0.5mmyr(-1) in the Lesser Himalaya and similar to 1 and 2mmyr(-1) in the High Himalaya, despite uniform hillslope angles. The erosion-rate data correlate with catchment-averaged values of 5 km radius relief, channel steepness indices, and specific stream power but to varying degrees of nonlinearity. Similar nonlinear relationships and coefficients of determination suggest that topographic steepness is the major control on the spatial variability of erosion and that twofold to threefold differences in annual runoff are of minor importance in this area. Instead, the spatial distribution of erosion in the study area is consistent with a tectonic model in which the rock uplift pattern is largely controlled by the shortening rate and the geometry of the Main Himalayan Thrust fault (MHT). Our data support a shallow dip of the MHT underneath the Lesser Himalaya, followed by a midcrustal ramp underneath the High Himalaya, as indicated by geophysical data. Finally, analysis of sample results from larger main stem rivers indicates significant variability of Be-10-derived erosion rates, possibly related to nonproportional sediment supply from different tributaries and incomplete mixing in main stem channels.
The Norwegian traffic network is impacted by about 2000 landslides, avalanches, and debris flows each year that incur high economic losses. Despite the urgent need to mitigate future losses, efforts to locate potential debris flow source areas have been rare at the regional scale. We tackle this research gap by exploring a minimal set of possible topographic predictors of debris flow initiation that we input to a Weights-of-Evidence (WofE) model for mapping the regional susceptibility to debris flows in western Norway. We use an inventory of 429 debris flows that were recorded between 1979 and 2008, and use the terrain variables of slope, total curvature, and contributing area (flow accumulation) to compute the posterior probabilities of local debris flow occurrence. The novelty of our approach is that we quantify the uncertainties in the WofE approach arising from different predictor classification schemes and data input, while estimating model accuracy and predictive performance from independent test data. Our results show that a percentile-based classification scheme excels over a manual classification of the predictor variables because differing abundances in manually defined bins reduce the reliability of the conditional independence tests, a key, and often neglected, prerequisite for the WofE method. The conditional dependence between total curvature and flow accumulation precludes their joint use in the model. Slope gradient has the highest true positive rate (88%), although the fraction of area classified as susceptible is very large (37%). The predictive performance, i.e. the reduction of false positives, is improved when combined with either total curvature or flow accumulation. Bootstrapping shows that the combination of slope and flow accumulation provides more reliable predictions than the combination of slope and total curvature, and helps refining the use of slope-area plots for identifying morphometric fingerprints of debris flow source areas, an approach used outside the field of landslide susceptibility assessments.
During the last glacial period lake Kinneret (the Sea of Galilee) fluctuated between high and low water levels reflecting the hydrological conditions of the lake watershed. Here, we focus on the hydrology of the lake after its retreat from the last glacial MIS2 (similar to 27-25 ka BP) highest stand of similar to 170 m below mean sea level (m bsl) to the low stand of similar to 214 m bsl at similar to 24-21 ka BP. The limnological-hydrological history of this time interval is recovered from trench and borehole that were dug and drilled in the southwestern shore of the lake at Ohalo-II archeological site. Cyprideis torosa (Ostracoda) recovered from the trench yielded elemental, Sr-87/Sr-86 and delta O-18 isotope data that provide information on the shore environment during the low stand period. The Sr-87/Sr-86 and Sr/Ca ratios in the ostracods, varying between similar to 0.70789 and similar to 0.70815 and 0.0017 and 0.0030, respectively indicate contributions of waters from the last glacial lake and regional runoff. The increase in the Sr-87/Sr-86 ratios reflects the decreasing effect of the last glacial Lake Kinneret waters and enhanced contribution of local runoff that washed down dried mountain soils that were previously developed during the wet and vegetated glacial. The lake retreat at similar to 24 ka BP coincided with the Heinrich event H2 at the northern Atlantic. H2 was expressed by severe aridity in Lake Kinneret-Dead Sea watershed. The limnological-hydrological change at post H2 was accompanied by 2% decrease in the delta O-18 value from -2% to -4% reflecting the change in the composition of the east Mediterranean rain sources.
The last glacial lake Sr-87/Sr-86 ratio is similar to the Tiberias Spa saline waters and distinctly different from the modern Lake Kinneret fresh waters: Sr-87/Sr-86 similar to 0.70785 compared to similar to 0.70760, respectively. This difference is explained by enhanced contribution of Ca-chloride brines with high Sr-87/Sr-86 values to the last glacial lake and reduced Jordan River contribution due to cold freezing conditions at its headwaters, while the modern Lake Kinneret is more affected by low Sr-87/Sr-86 freshwater from the Jordan watershed.
The water-level record from the 300 m deep paleo-lake Suguta (Northern Kenya Rift) during the African Humid Period (AHP, 15-5 ka BP) helps to explain decadal to centennial intensity variations in the West African Monsoon (WAM) and the Indian Summer Monsoon (ISM). This water-level record was derived from three different sources: (1) grain size variations in radiocarbon dated and reservoir corrected lacustrine sediments, (2) the altitudes and ages of paleo-shorelines within the basin, and (3) the results of hydro-balance modeling, providing important insights into the character of water level variations (abrupt or gradual) in the amplifier paleo-Lake Suguta. The results of these comprehensive analyses suggest that the AHP highstand in the Suguta Valley was the direct consequence of a northeastwards shift in the Congo Air Boundary (CAB), which was in turn caused by an enhanced atmospheric pressure gradient between East Africa and India during a northern hemisphere insolation maximum. Rapidly decreasing water levels of up to 90 m over less than a hundred years are best explained by changes in solar irradiation either reducing the East African-Indian atmospheric pressure gradient and preventing the CAB from reaching the study area, or reducing the overall humidity in the atmosphere, or a combination of both these effects. In contrast, although not well documented in our record we hypothesize a gradual end of the AHP despite an abrupt change in the source of precipitation when a decreasing pressure gradient between Asia and Africa prevented the CAB from reaching the Suguta Valley. The abruptness was probably buffered by a contemporaneous change in precession producing an insolation maximum at the equator during October. Whether or not this is the case, the water-level record from the Suguta Valley demonstrates the importance of both orbitally-controlled insolation variations and short-term changes in solar irradiation as factors affecting the significant water level variations in East African rift lakes.
We suggest a new clustering approach to classify focal mechanisms from large moment tensor catalogues, with the purpose of automatically identify families of earthquakes with similar source geometry, recognize the orientation of most active faults, and detect temporal variations of the rupture processes. The approach differs in comparison to waveform similarity methods since clusters are detected even if they occur in large spatial distances. This approach is particularly helpful to analyse large moment tensor catalogues, as in microseismicity applications, where a manual analysis and classification is not feasible. A flexible algorithm is here proposed: it can handle different metrics, norms, and focal mechanism representations. In particular, the method can handle full moment tensor or constrained source model catalogues, for which different metrics are suggested. The method can account for variable uncertainties of different moment tensor components. We verify the method with synthetic catalogues. An application to real data from mining induced seismicity illustrates possible applications of the method and demonstrate the cluster detection and event classification performance with different moment tensor catalogues. Results proof that main earthquake source types occur on spatially separated faults, and that temporal changes in the number and characterization of focal mechanism clusters are detected. We suggest that moment tensor clustering can help assessing time dependent hazard in mines.
Automated location of seismic events is a very important task in microseismic monitoring operations as well for local and regional seismic monitoring. Since microseismic records are generally characterized by low signal-to-noise ratio, automated location methods are requested to be noise robust and sufficiently accurate. Most of the standard automated location routines are based on the automated picking, identification and association of the first arrivals of P and S waves and on the minimization of the residuals between theoretical and observed arrival times of the considered seismic phases. Although current methods can accurately pick P onsets, the automatic picking of the S onset is still problematic, especially when the P coda overlaps the S wave onset. In this paper, we propose a picking free earthquake location method based on the use of the short-term-average/long-term-average (STA/LTA) traces at different stations as observed data. For the P phases, we use the STA/LTA traces of the vertical energy function, whereas for the S phases, we use the STA/LTA traces of a second characteristic function, which is obtained using the principal component analysis technique. In order to locate the seismic event, we scan the space of possible hypocentral locations and origin times, and stack the STA/LTA traces along the theoretical arrival time surface for both P and S phases. Iterating this procedure on a 3-D grid, we retrieve a multidimensional matrix whose absolute maximum corresponds to the spatial coordinates of the seismic event. A pilot application was performed in the Campania-Lucania region (southern Italy) using a seismic network (Irpinia Seismic Network) with an aperture of about 150 km. We located 196 crustal earthquakes (depth < 20 km) with magnitude range 1.1 < M-L < 2.7. A subset of these locations were compared with accurate manual locations refined by using a double-difference technique. Our results indicate a good agreement with manual locations. Moreover, our method is noise robust and performs better than classical location methods based on the automatic picking of the P and S waves first arrivals.
Global Positioning System (GPS) has been proven to be an effective tool to retrieve high-precision displacement for the natural hazard monitoring. The network positioning and Precise Point Positioning (PPP) are the two basic approaches for its data solution, but the former one can only get a relative displacement within the local reference frame and requires a complex and continuously linked infrastructure, and the latter one with a long convergence time to obtain the absolute displacements within the global reference frame. To overcome these drawbacks, this paper proposed a method of fast determining the displacement by PPP velocity estimation (PPPVE). The key of the approach is that the velocity vector parameters are not correlated with other unknown parameters, such as ambiguities and atmosphere, so they can be fast and accurately estimated and integrated into displacements. The validation shows that the displacement can be provided with a precision of 1-2 cm in 1 min by PPPVE. In additional, the Kalman smoothing estimation can be used to improve the PPP solution.
The growing body of research on large-scale mass wasting events so far has only scarcely investigated the sedimentology of chaotic deposits from non-volcanic terrestrial landslides such that any overarching and systematic terminological framework remains elusive. Yet recent work has emphasized the need for better understanding the internal structure and composition of rockslide deposits as a means to characterise the mechanics during the final stages of runout and emplacement. We offer a comprehensive overview on the occurrence of rock fragmentation and frictional melt both at different geographic locations, and different sections within large (>10(6) m(3)) rockslide masses. We argue that exposures of pervasively fragmented and interlocked jigsaw-cracked rock masses; basal melange containing rip-up clasts and phantom blocks; micro-breccia; and thin bands of basal frictionite are indispensable clues for identifying deposits from giant rockslides that may remain morphologically inconspicuous otherwise. These sedimentary assemblages are diagnostic tools for distinguishing large rockslide debris from macro and microscopically similar glacial deposits, tectonic fault-zone breccias, and impact breccias, and thus help avoid palaeoclimatic and tectonic misinterpretations, let alone misestimates of the hazard from giant rockslides. Moreover, experimental results from Mossbauer spectroscopy of frictionite samples support visual interpretations of thin sections, and demonstrate that short-lived (<10 s) friction-induced partial melting at temperatures >1500 degrees C in the absence of water occurred at the base of several giant moving rockslides. This finding supports previous theories of dry excess runout accompanied by comminution of rock masses down to gm-scale, and indicates that catastrophic motion of large fragmenting rock masses does not require water as a potential lubricant.
We apply and evaluate a recent machine learning method for the automatic classification of seismic waveforms. The method relies on Dynamic Bayesian Networks (DBN) and supervised learning to improve the detection capabilities at 3C seismic stations. A time-frequency decomposition provides the basis for the required signal characteristics we need in order to derive the features defining typical "signal" and "noise" patterns. Each pattern class is modeled by a DBN, specifying the interrelationships of the derived features in the time-frequency plane. Subsequently, the models are trained using previously labeled segments of seismic data. The DBN models can now be compared against in order to determine the likelihood of new incoming seismic waveform segments to be either signal or noise. As the noise characteristics of seismic stations varies smoothly in time (seasonal variation as well as anthropogenic influence), we accommodate in our approach for a continuous adaptation of the DBN model that is associated with the noise class. Given the difficulty for obtaining a golden standard for real data (ground truth) the proof of concept and evaluation is shown by conducting experiments based on 3C seismic data from the International Monitoring Stations, BOSA and LPAZ.
Coastal flood damage and adaptation costs under 21st century sea-level rise are assessed on a global scale taking into account a wide range of uncertainties in continental topography data, population data, protection strategies, socioeconomic development and sea-level rise. Uncertainty in global mean and regional sea level was derived from four different climate models from the Coupled Model Intercomparison Project Phase 5, each combined with three land-ice scenarios based on the published range of contributions from ice sheets and glaciers. Without adaptation, 0.2-4.6% of global population is expected to be flooded annually in 2100 under 25-123 cm of global mean sea-level rise, with expected annual losses of 0.3-9.3% of global gross domestic product. Damages of this magnitude are very unlikely to be tolerated by society and adaptation will be widespread. The global costs of protecting the coast with dikes are significant with annual investment and maintenance costs of US$ 12-71 billion in 2100, but much smaller than the global cost of avoided damages even without accounting for indirect costs of damage to regional production supply. Flood damages by the end of this century are much more sensitive to the applied protection strategy than to variations in climate and socioeconomic scenarios as well as in physical data sources (topography and climate model). Our results emphasize the central role of long-term coastal adaptation strategies. These should also take into account that protecting large parts of the developed coast increases the risk of catastrophic consequences in the case of defense failure.
Sedimentation in the floodplains of the Mekong Delta, Vietnam Part II: deposition and erosion
(2014)
Deposition and erosion play a key role in the determination of the sediment budget of a river basin, as well as for floodplain sedimentation. Floodplain sedimentation, in turn, is a relevant factor for the design of flood protection measures, productivity of agro-ecosystems, and for ecological rehabilitation plans. In the Mekong Delta, erosion and deposition are important factors for geomorphological processes like the compensation of deltaic subsidence as well as for agricultural productivity. Floodplain deposition is also counteracting the increasing climate change induced hazard by sea level rise in the delta. Despite this importance, a sediment database of the Mekong Delta is lacking, and the knowledge about erosion and deposition processes is limited. In the Vietnamese part of the Delta, the annually flooded natural floodplains have been replaced by a dense system of channels, dikes, paddy fields, and aquaculture ponds, resulting in floodplain compartments protected by ring dikes. The agricultural productivity depends on the sediment and associated nutrient input to the floodplains by the annual floods. However, no quantitative information regarding their sediment trapping efficiency has been reported yet. The present study investigates deposition and erosion based on intensive field measurements in three consecutive years (2008, 2009, and 2010). Optical backscatter sensors are used in combination with sediment traps for interpreting deposition and erosion processes in different locations. In our study area, the mean calculated deposition rate is 6.86kg/m(2) (approximate to 6mm/year). The key parameters for calculating erosion and deposition are estimated, i.e. the critical bed shear stress for deposition and erosion and the surface constant erosion rate. The bulk of the floodplain sediment deposition is found to occur during the initial stage of floodplain inundation. This finding has direct implications on the operation of sluice gates in order to optimize sediment input and distribution in the floodplains.
Sedimentation in the floodplains of the Mekong Delta, Vietnam. Part I: suspended sediment dynamics
(2014)
Suspended sediment is the primary source for a sustainable agro-ecosystem in the Mekong Delta by providing nutrient input for the subsequent cropping season. In addition, the suspended sediment concentration (SSC) plays an important role in the erosion and deposition processes in the Delta; that is, it influences the morphologic development and may counteract the deltaic subsidence and sea level rise. Despite this importance, little is known about the dynamics of suspended sediment in the floodplains of the Mekong Delta. In particular, quantitative analyses are lacking mainly because of data scarcity with respect to the inundation processes in the floodplains. In 2008, therefore, a comprehensive in situ system to monitor the dynamics of suspended sediment in a study area located in the Plain of Reeds was established, aiming at the characterization and quantification of suspended sediment dynamics in the deeply inundated parts of the Vietnamese part of the Mekong Delta. The monitoring system was equipped with seven water quality-monitoring stations. They have a robust design and autonomous power supply suitable for operation on inundated floodplains, enabling the collection of reliable data over a long period of time with a high temporal resolution. The data analysis shows that the general seasonal dynamics of suspended sediment transport in the Delta is controlled by two main mechanisms: the flood wave of the Mekong River and the tidal backwater influences from the coast. In the channel network, SSC decreases exponentially with distance from the Mekong River. The anthropogenic influence on SSC could also be identified for two periods: at the start of the floodplain inundation and at the end of the flood period, when subsequent paddy rice crops are prepared. Based on the results, we recommend an operation scheme for the sluice gates, which intends to distribute the sediment and thus the nutrients equally over the floodplain.
The soils on four lithologies (basaltic conglomerates, Bohio; Andesite; volcanoclastic sediments with basaltic agglomerates, Caimito volcanic; foraminiferal limestone, Caimito marine) on Barro Colorado Island (BCI have high exchangeable Ca concentrations and cation-exchange capacities (CEC) compared to other tropical soils on similar parent material. In the 0-10 cm layer of 24 mineral soils, pH values ranged from 5.7 (Caimito volcanic and Andesite) to 6.5 (Caimito marine), concentrations of exchangeable Ca from 134 mmol(c) kg(-1) (Caimito volcanic) to 585 mmolc kg-1 (Caimito marine), and cation exchange capacities from 317 mmol(c) kg(-1) (Caimito volcanic) to 933 mmol(c) kg(-1) (Caimito marine). X-ray diffractometry of the fraction <2 mu m revealed that smectites dominated the clay mineral assemblage in soil except on Caimito volcanic, where kaolinite was the dominant clay mineral. Exchangeable Ca concentrations decreased with increasing soil depth except on Caimito marine. The weathering indices Chemical Index of Alteration (CIA), Plagioclase Index of Alteration (PIA) and Weathering Index of Parker (WIP) determined for five soils on all geological formations, suggested that in contrast to expectation the topsoil (0-10 cm) appeared to be the least and the subsoil (50-70 cm) and saprolite (isomorphically weathered rock in the soil matrix) the most weathered. Additionally, the weathering indices indicated depletion of base cations and enrichment of Al-(hydr)oxides throughout the soil profile. Tree species did not have an effect on soil properties. Impeded leaching and the related occurrence of overland flow seem to be important in determining clay mineralogy. Our results suggest that (i) edaphic conditions favor the formation of smectites on most lithologies resulting in high CEC and thus high retention capacity for Ca and (ii) that there is an external source such as dust or sea spray deposition supplying Ca to the soils.
Our study aims at gaining insights into the processes determining the current treeline dynamics in Finnish Lapland. Using forest surveys conducted in 1978 and 2003 we modelled the occurrence and abundance of three dominant tree species in Finnish Lapland, i.e. Pinus sylvestris, Picea abies and Betula pubescens, with boosted regression trees. We assessed the importance of climatic, biotic and topographic variables in predicting tree occurrence and abundance based on their relative importance and response curves. We compared temporal and spatial transferability by using an extended transferability index.
Site fertility, the abundance of co-occurring species and growing degree days were generally the most important predictors for both occurrence and abundance across all species and datasets. Climatic predictors were more important for modelling occurrences than for modelling abundances. Occurrence models were able to reproduce the observed treeline pattern within one time period or region. Abundance models underestimated basal area but captured the general pattern of low and high values. Model performance as well as transferability differed considerably between species and datasets. Pinus sylvestris was modelled more successfully than P. abies and B. pubescens. Generally, spatial transferability was greater than temporal transferability. Comparing the environmental space between datasets revealed that transferring models means extrapolating to novel environments, providing a plausible explanation for limited transferability.
Our study illustrates how climate change can shift the environmental space and lead to limited model transferability. We identified non-climatic factors to be important in predicting the distribution of dominant tree species, contesting the widespread assumption of climatically induced range expansion.
The ratio of unsupported protactinium-231 to thorium-230 in marine sediments, (Pa/Th)(xs), is potentially sensitive to several processes of oceanographic and climatological interest: deep ocean circulation, marine biological productivity (as it relates to total particle flux) and particle composition (specifically, biogenic opal and authigenic Mn). In order to attribute variations in (Pa/Th)(xs) observed in sediment records to changes in specific processes through time, a better understanding of the chemical cycling of these elements in the modern ocean is necessary. To this end, a survey was undertaken of (Pa/Th)(xs) in surface sediments from the subarctic Pacific (SO202-INOPEX expedition) in combination with a Pacific-wide compilation of published data. Throughout the Pacific, (Pa/Th)(xs) is robustly correlated with the opal content of sediments. In the North and equatorial Pacific, simultaneous positive correlations with productivity indicators suggest that boundary scavenging and opal scavenging combine to enhance the removal of Pa in the eastern equatorial Pacific and subarctic Pacific. Deep ocean water mass ageing (>3.5 km) associated with the Pacific overturning appears to play a secondary role in determining the basin scale distribution of (Pa/Th)(xs). A basin-wide extrapolation of Pa removal is performed which suggests that the Pacific Pa budget is nearly in balance. We hypothesize that through time (Pa/Th)(xs) distributions in the Pacific could define the evolving boundaries of contrasting biogeographic provinces in the North Pacific, while the influence of hydrothermal scavenging of Pa potentially confounds this approach in the South Pacific.
Concerns about the regional impact of global climate change in a warming scenario have highlighted the gaps in our understanding of the Indian Summer Monsoon (ISM, also referred to as the Indian Ocean summer monsoon) and the absence of long term palaeoclimate data from the central Indian core monsoon zone (CMZ). Here we present the first high resolution, well-dated, multiproxy reconstruction of Holocene palaeoclimate from a 10 m long sediment core raised from the Lonar Lake in central India. We show that while the early Holocene onset of-intensified monsoon in the CMZ is similar to that reported from other ISM records, the Lonar data shows two prolonged droughts (PD, multidecadal to centennial periods of weaker monsoon) between 4.6-3.9 and 2-0.6 cal ka. A comparison of our record with available data from other ISM influenced sites shows that the impact of these PD was observed in varying degrees throughout the ISM realm and coincides with intervals of higher solar irradiance. We demonstrate that (i) the regional warming in the Indo-Pacific Warm Pool (IPWP) plays an important role in causing ISM PD through changes in meridional overturning circulation and position of the anomalous Walker cell; (ii) the long term influence of conditions like El Nino-Southern Oscillation (ENSO) on the ISM began only ca. 2 cal ka BP and is coincident with the warming of the southern IPWP; (iii) the first settlements in central India coincided with the onset of the first PD and agricultural populations flourished between the two PD, highlighting the significance of natural climate variability and PD as major environmental factors affecting human settlements.
Due to its high spatial and temporal variability, preferential flow is difficult to measure and quantify. Earthworms create macropores that provide common pathways for preferential flow. Therefore in this article, we link earthworm abundance to macropore numbers and hydrological effectiveness, with the future aim to use species distribution models of earthworms for the spatial parameterization of preferential flow.
Earthworms are generally categorized into three ecological types with varying burrowing behaviour, resulting in a different impact on soil hydrological processes. Therefore, we studied the relationships between the abundance of the earthworm ecological types and macropores of different size classes and in different soil depths. The abundance and biomass of earthworms were well correlated to different sizes of macropores in different soil depths. This is mainly the case for the larger, vertically oriented macropores (>6mm diameter), which are generally connected to the soil surface and hydrologically most effective. The correlation of total earthworm biomass and macropores ranges from 072 to 089 for different soil depths.
Although there is quite some variation in infiltration patterns, infiltration from macropores into the matrix is profile-specific, as it varies strongly between profiles, but not within one profile. Macropore coating seems to have a larger effect on this macropore matrix interaction than the soil physical properties of the matrix. Although the amount of macropores and their effectiveness are clearly related to the earthworm distribution, the variation in infiltration from macropores to soil matrix should be further studied.
We investigate the usefulness of complex flood damage models for predicting relative damage to residential buildings in a spatial and temporal transfer context. We apply eight different flood damage models to predict relative building damage for five historic flood events in two different regions of Germany. Model complexity is measured in terms of the number of explanatory variables which varies from 1 variable up to 10 variables which are singled out from 28 candidate variables. Model validation is based on empirical damage data, whereas observation uncertainty is taken into consideration. The comparison of model predictive performance shows that additional explanatory variables besides the water depth improve the predictive capability in a spatial and temporal transfer context, i.e., when the models are transferred to different regions and different flood events. Concerning the trade-off between predictive capability and reliability the model structure seem more important than the number of explanatory variables. Among the models considered, the reliability of Bayesian network-based predictions in space-time transfer is larger than for the remaining models, and the uncertainties associated with damage predictions are reflected more completely.
We present pollen-based reconstructions of the spatio-temporal dynamics of northern European regional vegetation abundance through the Holocene. We apply the Regional Estimates of VEgetation Abundance from Large Sites (REVEALS) model using fossil pollen records from eighteen sites within five modern biomes in the region. The eighteen sites are classified into four time-trajectory types on the basis of principal components analysis of both the REVEALS-based vegetation estimates (RVs) and the pollen percentage (PPs). The four trajectory types are more clearly separated for RVs than PPs. Further, the timing of major Holocene shifts, rates of compositional change, and diversity indices (turnover and evenness) differ between RVs and PPs. The differences are due to the reduction by REVEALS of biases in fossil pollen assemblages caused by different basin size, and inter-taxonomic differences in pollen productivity and dispersal properties. For example, in comparison to the PPs, the RVs show an earlier increase in Corylus and Ulmus in the early-Holocene and a more pronounced increase in grassland and deforested areas since the mid-Holocene. The results suggest that the influence of deforestation and agricultural activities on plant composition and abundance from Neolithic times was stronger than previously inferred from PPs. Relative to PPs, RVs show a more rapid compositional change, a largest decrease in turnover, and less variable evenness in most of northern Europe since 5200 cal yr BP. All these changes are primarily related to the strong impact of human activities on the vegetation. This study demonstrates that RV-based estimates of diversity indices, timing of shifts, and rates of change in reconstructed vegetation provide new insights into the timing and magnitude of major human distribution on Holocene regional, vegetation, feature that are critical in the assessment of human impact on vegetation, land-cover, biodiversity, and climate in the past.
Deep seepage estimation is important for water balance investigations of groundwater and the vadose zone. A simplified Buckingham-Darcy method to assess time series of deep seepage fluxes was proposed by Schindler and Muller (1998). In the method dynamics of water fluxes are calculated by a soil hydraulic conductivity function. Measured soil moistures and matric heads are used as input data. Resulting time series of flux dynamics are scaled to realistic absolute levels by calibrating the method with the areal water balance. An assumption of the method is that water fluxes at different positions exhibit identical dynamics although their absolute values can differ. The aim of this study was to investigate uncertainties of that method depending on the particle size distribution and textural heterogeneity in non-layered soils. We performed a numerical experiment using the two-dimensional Richards Equation. A basic model of transient water fluxes beneath the root and capillary zone was setup and used to simulate time series of soil moisture, matric head, and seepage fluxes for 4221 different cases of particle size distribution and intensities of textural heterogeneity. Soil hydraulic parameters were predicted by the pedotransfer function Rosetta. Textural heterogeneity was modeled with Miller and Miller scaling factors arranged in spatial random fields. Seepage fluxes were calculated with the Buckingham-Darcy method from simulated soil moisture and matric head time series and compared with simulated reference fluxes. The median of Root Mean Square Error was about 0.026 cm d(-1) and the median of maximum cross correlation was 0.96 when the method was calibrated adequately. The method's performance was mainly influenced by (i) the soil textural class and (ii) the time period used for flux calibration. It performed best in sandy loam while hotspots of errors occurred in sand and silty texture. Calibrating the method with time periods that exhibit high variance of seepage fluxes yielded the best performance. The geostatistical properties of the Miller and Miller scaling field influenced the performance only slightly. However, the Miller and Miller scaling procedure generated heterogeneous flow fields that were addressed as main reason for mismatches of simulated reference fluxes and fluxes obtained with the Buckingham-Darcy method.
Many coasts feature sequences of Quaternary and Neogene shorelines that are shaped by a combination of sea-level oscillations and tectonics. We compiled a global synthesis of sea-level changes for the following highstands: MIS 1, MIS 3, MIS 5e and MIS 11. Also, we date the apparent onset of sequences of paleoshorelines either from published data or tentatively extrapolating an age for the uppermost, purported oldest shoreline in each sequence. Including the most documented MIS 5e benchmark, we identify 926 sequences out of which 185 also feature Holocene shorelines. Six areas are identified where elevations of the MIS 3 shorelines are known, and 31 feature elevation data for MIS 11 shorelines. Genetic relationships to regional geodynamics are further explored based on the elevations of the MIS 5e benchmark. Mean apparent uplift rates range from 0.01 0.01 mm/yr (hotspots) to 1.47 0.08 mm/yr (continental collision). Passive margins appear as ubiquitously uplifting, while tectonic segmentation is more important on active margins. From the literature and our extrapolations, we infer ages for the onset of formation for -180 coastal sequences. Sea level fingerprinting on coastal sequences started at least during mid Miocene and locally as early as Eocene. Whether due to the changes in the bulk volume of seawater or to the temporal variations in the shape of ocean basins, estimates of eustasy fail to explain the magnitude of the apparent sea level drop. Thus, vertical ground motion is invoked, and we interpret the longlasting development of those paleoshore sequences as the imprint of glacial cycles on globally uplifted margins in response to continental compression. The geomorphological expression of the sequences matches the amplitude and frequency of glacial cyclicity. From middle Pleistocene to present-day, moderately fast (100,000 yrs) oscillating sea levels favor the development of well identified strandlines that are distinct from one another. Pliocene and Lower Pleistocene strandlines associated with faster cyclicity (40,000 yrs) are more compact and easily merge into rasas, whereas older Cenozoic low-frequency eustatic changes generally led to widespread flat-lying coastal plains.
The Mohorovicic discontinuity, Moho for short, which marks the boundary between crust and mantle, is the main first-order structure within the lithosphere. Geodynamics and tectonic evolution determine its depth level and properties. Here, we present a map of the Moho in central Europe across the Teisseyre-Tornquist Zone, a region for which a number of previous studies are available. Our results are based on homogeneous and consistent processing of P-and S-receiver functions for the largest passive seismological data set in this region yet, consisting of more than 40 000 receiver functions from almost 500 station. Besides, we also provide new results for the crustal vP/vS ratio for the whole area.
Our results are in good agreement with previous, more localized receiver function studies, as well as with the interpretation of seismic profiles, while at the same time resolving a higher level of detail than previous maps covering the area, for example regarding the Eifel Plume region, Rhine Graben and northern Alps. The close correspondence with the seismic data regarding crustal structure also increases confidence in use of the data in crustal corrections and the imaging of deeper structure, for which no independent seismic information is available.
In addition to the pronounced, stepwise transition from crustal thicknesses of 30 km in Phanerozoic Europe to more than 45 beneath the East European Craton, we can distinguish other terrane boundaries based on Moho depth as well as average crustal v(P)/v(S) ratio and Moho phase amplitudes. The terranes with distinct crustal properties span a wide range of ages, from Palaeoproterozoic in Lithuania to Cenozoic in the Alps, reflecting the complex tectonic history of Europe. Crustal thickness and properties in the study area are also markedly influenced by tectonic overprinting, for example the formation of the Central European Basin System, and the European Cenozoic Rift System. In the areas affected by Cenozoic rifting and volcanism, thinning of the crust corresponds to lithospheric updoming reported in recent surface wave and S-receiver function studies, as expected for thermally induced deformation. The same correlation applies for crustal thickening, not only across the Trans-European Suture Zone, but also within the southern part of the Bohemian Massif.
Costing natural hazards
(2014)
Sedimentary proxies used to reconstruct marine productivity suffer from variable preservation and are sensitive to factors other than productivity. Therefore, proxy calibration is warranted. Here we map the spatial patterns of two paleoproductivity proxies, biogenic opal and barium fluxes, from a set of core-top sediments recovered in the Subarctic North Pacific. Comparisons of the proxy data with independent estimates of primary and export production, surface water macronutrient concentrations, and biological pCO(2) drawdown indicate that neither proxy shows a significant correlation with primary or export productivity for the entire region. Biogenic opal fluxes, when corrected for preservation using Th-230-normalized accumulation rates, show a good correlation with primary productivity along the volcanic arcs (tau = 0.71, p = 0.0024) and with export productivity throughout the western Subarctic North Pacific (tau = 0.71, p = 0.0107). Moderate and good correlations of biogenic barium flux with export production (tau = 0.57, p = 0.0022) and with surface water silicate concentrations (tau = 0.70, p = 0.0002) are observed for the central and eastern Subarctic North Pacific. For reasons unknown, however, no correlation is found in the western Subarctic North Pacific between biogenic barium flux and the reference data. Nonetheless, we show that barite saturation, uncertainty in the lithogenic barium corrections, and problems with the reference data sets are not responsible for the lack of a significant correlation between biogenic barium flux and the reference data. Further studies evaluating the factors controlling the variability of the biogenic constituents in the sediments are desirable in this region.
The century-long debate over the origins of inner gorges that were repeatedly covered by Quaternary glaciers hinges upon whether the gorges are fluvial forms eroded by subaerial rivers, or subglacial forms cut beneath ice. Here we apply cosmogenic nuclide exposure dating to seven inner gorges along similar to 500 km of the former Fennoscandian ice sheet margin in combination with a new deglaciation map. We show that the timing of exposure matches the advent of ice-free conditions, strongly suggesting that gorges were cut by channelized subglacial meltwater while simultaneously being shielded from cosmic rays by overlying ice. Given the exceptional hydraulic efficiency required for meltwater channels to erode bedrock and evacuate debris, we deduce that inner gorges are the product of ice sheets undergoing intense surface melting. The lack of postglacial river erosion in our seven gorges implicates subglacial meltwater as a key driver of valley deepening on the Baltic Shield over multiple glacial cycles.
In this paper we quantify the sediment dynamics in the formerly glaciated Zielbach catchment in the Italian Alps from the end of the Last Glacial Maximum (LGM) until today. As a basis for our quantification, we use the stratigraphic record offered by a 3.5 km(2) large fan that we explore with a seismic survey, stratigraphic analyses of drillhole material, and C-14 ages measured on organic matter encountered in these drillings. In addition, we calculate past denudation rate variability in the fan deposits using concentrations of cosmogenic Be-10. We merge this information into a scenario of how the sediment flux has changed through time and how this variability can be related to climatic variations, framed within well-known paraglacial models. The results document a highly complex natural system. From the LGM to the very early Holocene, ice-melted discharge and climate variability promoted a high sediment flux (sedimentation rate up to 40 mm/yr). This flux then dramatically decreased toward interglacial values (0.8 mm/yr at 5-4 calibrated kyr B.P.). However, in contrast to the trend of classic paraglacial models, the flux recorded at Zielbach shows secondary peaks at 6.5 ka and 2.5 ka, with values of 13 mm/yr and 1.5 mm/yr, respectively. Paleo-denudation rates also decrease from similar to 33 mm/yr at the beginning of the Holocene to 0.42 mm/yr at 5 ka, with peaks of similar to 6 mm/yr and 1.1 mm/yr at 6.5 ka and 2.5 ka. High-amplitude climate change is the most likely cause of the secondary peaks, but anthropogenic activities may have contributed as well. The good correlation between paleo-sedimentation and paleo-denudation rates suggests that the majority of the deglaciated material destocked from the Zielbach catchment is stored in the alluvial fan.
Some of the largest and most erosive floods on Earth result from the failure of glacial dams. While potentially cataclysmic ice dams are recognized to have repeatedly formed along ice-sheet margins, much less is known about the frequency and longevity of ice dams caused by mountain glaciers, and their impact on landscape evolution. Here we present field observations and results from cosmogenic nuclide dating that allow reconstructing a > 100-k.y.-long history of glacial damming in the Shyok Valley, eastern Karakoram (South Asia). Our field observations provide evidence that Asia's second-longest glacier, the Siachen, once extended for over 180 km and blocked the Shyok River during the penultimate glacial period, leading to upstream deposition of a more than 400-m-thick fluvio-lacustrine valley fill. Be-10-depth profile modeling indicates that glacial damming ended with the onset of the Eemian interglacial and that the Shyok River subsequently incised the valley fill at an average rate of similar to 4-7 m k.y.(-1). Comparison with contemporary ice-dammed lakes in the Karakoram and elsewhere suggests recurring outburst floods during the aggradation period, while over 25 cycles of fining-upward lake deposits within the valley fill indicate impounding of floods from farther upstream. Despite prolonged damming, the net effect of this and probably earlier damming episodes by the Siachen Glacier is dominated by glacial erosion in excess of fluvial incision, as evidenced by a pronounced overdeepening that follows the glaciated valley reach. Strikingly similar overdeepened valleys at all major confluences of the Shyok and Indus Rivers with Karakoram tributaries indicate that glacial dams and subsequent outburst floods have been widespread and frequent in this region during the Quaternary. Our study suggests that the interaction of Karakoram glaciers with the Shyok and Indus Rivers promoted valley incision and headward erosion into the western margin of the Tibetan Plateau.
Tile drains strongly influence the water cycle in agricultural catchment in terms of water quantity and quality. The connectivity of preferential flow to tile drains can create shortcuts for rapid transport of solutes into surface waters. The leaching of pesticides can be linked to a set of main factors including, rainfall characteristics, soil moisture, chemical properties of the pesticides, soil properties, and preferential flow paths. The connectivity of the macropore system to the tile drain is crucial for pesticide leaching. Concurring influences of the main factors, threshold responses and the role of flow paths are still poorly understood. The objective of this study is to investigate these influences by a replica series of three irrigation experiments on a tile drain field site using natural and artificial tracers together with applied pesticides. We found a clear threshold behavior in the initialization of pesticide transport that was different between the replica experiments. Pre-event soil water contributed significantly to the tile drain flow, and creates a flow path for stored pesticides from the soil matrix to the tile drain. This threshold is controlled by antecedent soil moisture and precipitation characteristics, and the interaction between the soil matrix and preferential flow system. Fast transport of pesticides without retardation and the remobilization could be attributed to this threshold and the interaction between the soil matrix and the preferential flow system. Thus, understanding of the detailed preferential flow processes clearly enhances the understanding of pesticide leaching on event and long term scale, and can further improve risk assessment and modeling approaches. (C) 2014 Elsevier B.V. All rights reserved.
In this paper we present the sedimentary features and the luminescence chronology for Nahal Mahanayeem Outlet (NMO), an archaeological open air site at the southern margin of the Hula Basin (Northern Jordan Rift Valley, Israel). The site is characterized by a lithic assemblage ascribed primarily to the Middle Paleolithic Mousterian tradition, and by an excellent preservation of floral and faunal remains.
Six geological units forming the stratigraphic sequence of the site were distinguished: (Unit 6) archaeologically sterile, light-colored limnic carbonates; (Unit 5) conglomerates of rounded basalt boulders and cobbles forming a hill-like topography; (Unit 4 and Unit 3) a sequence of similar dark silty sediments, attached to and overlaying the conglomerates, containing the archaeological horizons of the site; (Unit 2) a number of channels cutting into the top of Unit 3, filled with coarse sand and rounded basalt and limestone gravels of fluvial origin; and (Unit 1) a thin sand layer laid down by the present-day Jordan River covering another unconformity as a result of heavy machinery drainage operations in 1999.
The OSL age for Unit 6 yielded a minimum age older than 460 ka. Sedimentary features and the embedded fossils suggest that Unit 1 can be linked to the Early Pleistocene Gadot Chalk. Unit 5 represents a local geological feature and could be an indicator for a period of increased erosion with formation of coarse grained sediments. The archaeological horizons form the lower parts of Unit 4 and yielded OSL-ages between 55 and 65 ka, indicating an affiliation to the sediments called "Ashmura Formation" with an Upper Pleistocene age for the site. The channel fills of Unit 2 can be dated by the recovered artifacts. These range in age from the Upper Paleolithic (Aurignacian) to historic times. Unit 1 is recent.
The study of the complex NMO stratigraphy, combined with coherent OSL chronology, has enabled us to reconstruct parts of the geological history of the Hula Basin during the Late Pleistocene. It is this history that forms the background for the human migration and utilization of natural resources in the Upper Jordan Rift Valley. (C) 2013 Elsevier Ltd and INQUA. All rights reserved.
(Paleo-)climatologists are challenged to identify mechanisms that cause the observed abrupt Holocene monsoon events despite the fact that monsoonal circulation is assumed to be driven by gradual insolation changes. Here we provide proxy and model evidence to show that moisture-advection feedback can lead to a non-linear relationship between sea-surface and continental temperatures and monsoonal precipitation. A pollen record from Lake Ximencuo (Nianbaoyeze Mountains) indicates that vegetation from the eastern margin of the Tibetan Plateau was characterized by alpine deserts and glacial flora after the Last Glacial Maximum (LGM) (21-15.5 cal kyr BP), by alpine meadows during the Late Glacial (15.5-10.4 cal kyr BP) and second half of the Holocene (5.0 cal kyr BP to present) and by mixed forests during the first half of the Holocene (10.4-5.0 cal kyr BP). The application of pollen-based transfer functions yields an abrupt temperature increase at 10.4 cal kyr BP and a decrease at 5.0 cal kyr BP of about 3 degrees C. By applying endmember modeling to grain-size data from the same sediment core we infer that frequent fluvial events (probably originating from high-magnitude precipitation events) were more common in the early and mid Holocene. We assign the inferred exceptional strong monsoonal circulation to the initiation of moisture-advection feedback, a result supported by a simple model that reproduces this feedback pattern over the same time period. (C) 2014 Published by Elsevier B.V.
Kinematic interaction of faults is an important issue for detailed seismic hazard assessments in seismically active regions. The Eastern Mosha Fault (EMF) and the North Tehran Fault (NTF) are two major active faults of the southern central Alborz mountains, located in proximity of Tehran (population similar to 9 million). We used field, geomorphological and paleoseismological data to explore the kinematic transition between the faults, and compare their short-term and long-term history of activity. We introduce the Niknamdeh segment of the NTF along which the strike-slip kinematics of EMF is transferred onto the NTF, and which is also responsible for the development of a pull-apart basin between the eastern segments of the NTF. The Ira trench site at the linkage zone between the two faults reveals the history of interaction between rock avalanches, active faulting and sag-pond development. The kinematic continuity between the EMF and NTF requires updating of seismic hazard models for the NTF, the most active fault adjacent to the Tehran Metropolis. Study of offsets of large-scale morphological features along the EMF, and comparison with estimated slip rates along the fault indicates that the EMF has started its left-lateral kinematics between 3.2 and 4.7 Ma. According to our paleoseismological data and the morphology of the nearby EMF and NTF, we suggest minimum and maximum values of about 1.8 and 3.0 mm/year for the left-lateral kinematics on the two faults in their linkage zone, averaged over Holocene time scales. Our study provides a partial interpretation, based on available data, for the fault activity in northeastern Tehran region, which may be completed with studies of other active faults of the region to evaluate a more realistic seismic hazard analysis for this heavily populated major city. (C) 2014 Elsevier B.V. All rights reserved.
A wide range of basic and applied problems in water resources research requires high-quality estimates of the spatial mean of throughfall. Many throughfall sampling schemes, however, are not optimally adapted to the system under study. The application of inappropriate sampling schemes may partly reflect the lack of generally applicable guidelines on throughfall sampling strategies. In this study we conducted virtual sampling experiments using simulated fields which are based on empirical throughfall data from three structurally distinct forests (a 12-year old teak plantation, a 5-year old young secondary forest, and a 130-year old secondary forest). In the virtual sampling experiments we assessed the relative error of mean throughfall estimates for 38 different throughfall sampling schemes comprising a variety of funnel- and trough-type collectors and a large range of sample sizes. Moreover, we tested the performance of each scheme for both event-based and accumulated throughfall data. The key findings of our study are threefold. First, as errors of mean throughfall estimates vary as a function of throughfall depth, the decision on which temporal scale (i.e. event-based versus accumulated data) to sample strongly influences the required sampling effort. Second, given a chosen temporal scale throughfall estimates can vary considerably as a function of canopy complexity. Accordingly, throughfall sampling in simply structured forests requires a comparatively modest effort, whereas heterogeneous forests can be extreme in terms of sampling requirements, particularly if the focus is on reliable data of small events. Third, the efficiency of trough-type collectors depends on the spatial structure of throughfall. Strong, long-ranging throughfall patterns decrease the efficiency of troughs substantially. Based on the results of our virtual sampling experiments, which we evaluated by applying two contrasting sampling approaches simultaneously, we derive readily applicable guidelines for throughfall monitoring. (C) 2014 Elsevier B.V. All rights reserved.
The Relative Pollen Productivities (RPPs) of common steppe species are estimated using Extended R-value (ERV) model based on pollen analysis and vegetation survey of 30 surface soil samples from typical steppe area of northern China. Artemisia, Chenopodiaceae, Poaceae, Cyperaceae, and Asteraceae are the dominant pollen types in pollen assemblages, reflecting the typical steppe communities well. The five dominant pollen types and six common types (Thalictrum, Iridaceae, Potentilla, Ephedra, Brassicaceae, and Ulmus) have strong wind transport abilities; the estimated Relevant Source Area of Pollen (RSAP) is ca. 1000 m when the sediment basin radius is set at 0.5 m. Ulmus, Artemisia, Brassicaceae, Chenopodiaceae, and Thalictrum have relative high RPPs; Poaceae, Cyperaceae, Potentilla, and Ephedra pollen have moderate RPPs; Asteraceae and Iridaceae have low RPPs. The reliability test of RPPs revealed that most of the RPPs are reliable in past vegetation reconstruction. However, the RPPs of Asteraceae and Iridaceae are obviously underestimated, and those of Poaceae, Chenopodiaceae, and Ephedra are either slightly underestimated or slightly overestimated, suggesting that those RPPs should be considered with caution. These RPPs were applied to estimating plant abundances for two fossil pollen spectra (from the Lake Bayanchagan and Lake Haoluku) covering the Holocene in typical steppe area, using the "Regional Estimates of Vegetation Abundance from Large Sites" (REVEALS) model. The RPPs-based vegetation reconstruction revealed that meadow-steppe dominated by Poaceae, Cyperaceae, and Artemisia plants flourished in this area before 6500-5600 cal yr BP, and then was replaced by present typical steppe.
Two lines of research are combined in this study: first, the development of tools for the temporal disaggregation of precipitation, and second, some newer results on the exponential scaling of heavy short-term precipitation with temperature, roughly following the Clausius-Clapeyron (CC) relation. Having no extra temperature dependence, the traditional disaggregation schemes are shown to lack the crucial CC-type temperature dependence. The authors introduce a proof-of-concept adjustment of an existing disaggregation tool, the multiplicative cascade model of Olsson, and show that, in principal, it is possible to include temperature dependence in the disaggregation step, resulting in a fairly realistic temperature dependence of the CC type. They conclude by outlining the main calibration steps necessary to develop a full-fledged CC disaggregation scheme and discuss possible applications.
Due to large uncertainties and non-uniqueness in fault slip inversion, the investigation of stress coupling based on the direct comparison of independent slip inversions, for example, between the coseismic slip distribution and the interseismic slip deficit, may lead to ambiguous conclusions. In this study, we therefore adopt the stress-constrained joint inversion in the Bayesian approach of Wang et al., and implement the physical hypothesis of stress coupling as a prior. We test the hypothesis that interseismic locking is coupled with the coseismic rupture, and the early post-seismic deformation is a stress relaxation process in response to the coseismic stress perturbation. We characterize the role of stress coupling in the seismic cycle by evaluating the efficiency of the model to explain the available data. Taking the 2004 M6 Parkfield earthquake as a study case, we find that the stress coupling hypothesis is in agreement with the data. The coseismic rupture zone is found to be strongly locked during the interseismic phase and the post-seismic slip zone is indicated to be weakly creeping. The post-seismic deformation plays an important role to rebuild stress in the coseismic rupture zone. Based on our results for the stress accumulation during both inter- and post-seismic phase in the coseismic rupture zone, together with the coseismic stress drop, we estimate a recurrence time of M6 earthquake in Parkfield around 23-41 yr, suggesting that the duration of 38 yr between the two recent M6 events in Parkfield is not a surprise.
An expanded ostracod-based conductivity transfer function for climate reconstruction in the Levant
(2014)
We present the first modern calibration dataset linking ostracod assemblage composition to water chemistry, and other site-specific variables, in the hydrologically and geopolitically sensitive southern Levant region. A total of 42 ostracod taxa were recorded from the 178 sampled sites in Israel and Jordan. Ilyocypris spp., Heterocypris salina and Cypridopsis vidua are the most abundant taxa. Species strictly confined to freshwater conditions are Prionocypris zenkeri, Gomphocythere ortali and Prionocypris olivaceus. In contrast, H. sauna, Bradleytriebella lineata and Cyprideis torosa show high frequencies in brackish waters (waters with higher conductivity). Humphcypris subterranea, G. ortali, P. olivaceus and Cypridopsis elongata apparently prefer flowing waters. Specific conductivity optima and tolerance ranges were calculated for the recorded ostracod species and may be used for the palaeoenvironmental assessment of fossil ostracod assemblages. In addition, a transfer-function for quantitative specific conductivity estimation based on 141 samples was established with weighted averaging partial least squares regression (WA-PLS). The resulting coefficient of determination r(2) between observed and predicted conductivity values (0.72) and the root-mean-square error of prediction (RMSEP) in % gradient length (13.1) indicate that conductivity may be reliably estimated from ostracod assemblage data. The transfer function was first applied to last glacial ostracod assemblage data from an archaeological trench in the Sea of Galilee (northern Israel). Relatively large conductivity fluctuations between ca 1 and 7 mS cm(-1) were inferred for the period 24-20 cal ka BP. In addition, four episodes of freshwater influx near the site of the trench were identified from the presence of shells of freshwater and stream-dwelling species intermingled with very abundant shells of Cyprideis torosa. The results of our study allow a better use of Quaternary ostracods from the Levant as palaeoenvironmental indicators of water-body types and past conductivity levels and will contribute to a better understanding of Quaternary environmental and climate change in the Levant. (C) 2014 Elsevier Ltd. All rights reserved.
In this study, we suggest a novel approach for the retrieval of regional moment tensors for earthquakes with small to moderate magnitudes. The first modification is the combined inversion of broadband and short-period waveform data. The broadband waveforms are inverted in a frequency range suitable for surface waves, whereas for the short-period data a frequency range suitable for body waves is applied. The second modification is the use of first-motion body-wave polarities to select the most probable solution out of all solutions from inversion. To combine three different criteria for selecting the most probable solution (i.e., residual from inversion, double-couple content of solution, number of nonmatching first-motion body-wave polarities), the L2 norm is applied to the normalized parameters. We chose five earthquakes within the Alborz mountains, Iran, as a case study (3.1 <= M-w <= 4.1). In this area, several factors exacerbate the difficulty of performing inversion for moment tensors, for example, a heterogeneous station network and large azimuthal gaps. We have demonstrated that our approach supplies reliable moment tensors when inversion from broadband data alone fails. In one case, we successfully retrieved a stable solution from short-period waveform data alone. Thus, our approach enables successful determination of seismic moment tensors wherever a sparse network of broadband stations has thus far prevented it.
We examine the use of ambient noise cross-correlation tomography for shallow site characterization using a modified two-step approach. Initially, we extract Rayleigh wave traveltimes from correlation traces of vertical component seismic recordings from a local network installed in Mygdonia basin, northern Greece. The obtained Rayleigh wave traveltimes show significant spatial variability, as well as distance and frequency dependence due to the 3-D structure of the area, dispersion, and anelastic attenuation effects. The traveltime data sets are inverted through a surface wave tomography approach to determine group velocity maps for each frequency. The proposed tomographic inversion involves the use of approximate Fresnel volumes and interfrequency smoothing constraints to stabilize the results. In the last step, we determine a final 3-D velocity model using a node-based Monte Carlo 1-D dispersion curve inversion. The reliability of the final 3-D velocity model is examined by spatial and depth resolution analysis, as well as by inversion for different model parameterizations. The obtained results are in very good agreement with previous findings from seismic and other geophysical methods. The new 3-D VS model provides additional structural constraints for the shallow sediments and bedrock structure of the northern Mygdonia basin up to the depth of similar to 200-250 m. Present work results suggest that the migration of ambient tomography techniques from large scales (tens or hundreds of km) to local scales (few hundred meters) is possible but cannot be used as a black box technique for 3-D modeling and detailed geotechnical site characterization.
In a study from 2008, Lariviere and colleagues showed, for the field of natural sciences and engineering, that the median age of cited references is increasing over time. This result was considered counterintuitive: with the advent of electronic search engines, online journal issues and open access publications, one could have expected that cited literature is becoming younger. That study has motivated us to take a closer look at the changes in the age distribution of references that have been cited in water resources journals since 1965. Not only could we confirm the findings of Lariviere and colleagues. We were also able to show that the aging is mainly happening in the oldest 10-25% of an average reference list. This is consistent with our analysis of top-cited papers in the field of water resources. Rankings based on total citations since 1965 consistently show the dominance of old literature, including text books and research papers in equal shares. For most top-cited old-timers, citations are still growing exponentially. There is strong evidence that most citations are attracted by publications that introduced methods which meanwhile belong to the standard toolset of researchers and practitioners in the field of water resources. Although we think that this trend should not be overinterpreted as a sign of stagnancy, there might be cause for concern regarding how authors select their references. We question the increasing citation of textbook knowledge as it holds the risk that reference lists become overcrowded, and that the readability of papers deteriorates.
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.
In Lutzito catchment on Barro Colorado Island, Panama, extraordinarily high suspended-sediment yields of 1-2Mgha-1year-1 were generated despite the dense forest cover coinciding with erosion-resistant soils. We hypothesized that ant mounding activity is an important zoogeomorphological mechanism in this area, providing relevant quantities of easily transportable material at the soil surface. To test this hypothesis, all ant mound material was collected collected for dry mass determination from thirty 4m2 plots installed in the study area every 1-3days during the 39-day sampling period. Additionally, three ground-nesting ant species responsible for mounds in the study area, Ectatomma ruidum, Trachymyrmex cornetzi and Strumigenys marginiventris, were identified. On the basis of the total of 1.38kg of material collected in the wet season of 2011, the estimate for the whole 8months wet season amounts to 725kgha-1. As this value is in the same order of magnitude as sediment output, it shows that ants may act as important ecosystem engineers and contribute to sediment production here by providing large quantities of fine-grained, readily erodible material at the soil surface for subsequent transport to the streambed. Copyright (c) 2014 John Wiley & Sons, Ltd.