Refine
Year of publication
- 2014 (228) (remove)
Document Type
- Article (196)
- Doctoral Thesis (19)
- Preprint (8)
- Review (4)
- Postprint (1)
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 (228) (remove)
The hydrology of western Central Asia is highly sensitive to climatic perturbations. In order to understand its long-term variability and to infer linkages between precipitation and atmospheric and oceanic systems, we conducted a thorough sedimentary and geochemical study on a composite core retrieved in lake Son Kul (central Kyrgyzstan). A multi-proxy approach was conducted on lake sediments based on grain size analyses, magnetic susceptibility, total organic carbon (TOC), total nitrogen (TN) and carbon and oxygen isotope analyses on bulk and biogenic materials (ostracoda and molluscs shells) at a resolution equivalent to ca 40 years, aiming to characterise the sequence of palaeolimnological changes in Son Kul.
As indicated by delta O-18 record of bulk carbonates, mainly consisting of aragonite, the Holocene hydrological balance was negative during most of time, suggesting an excess of evaporation (E) over precipitation (P). Limnological conditions fluctuated rapidly before 5000 cal yr BP indicating significant changes in regional hydrology and climate. In particular, the long-term negative hydrological balance was impeded by several short stages with marked increase of precipitation, lasting several decades to a few centuries (e.g., 8300-8200, 6900-6700, 6300-6100, 5500-5400, 5300-5200 and 3100 -3000 cal yr BP). Precipitation changes as inferred from 8180 data are also documented by increased minerogenic detritus and higher TOC. We propose that the seasonal pattern of precipitation varied transiently in western Central Asia during the Holocene, although evaporation changes may also account for the rapid changes observed in delta O-18 data. When the annual water balance was less critical (P <= E), the excess of water might be ascribed to increased precipitation during cold seasons mainly because winter precipitation has more negative delta O-18 than its summer equivalent. Conversely, when the annual water balance is negative (P E), the moisture was mainly delivered during the warm season, as between 5000 and 2000 cal yr BP. Our results thus imply that moisture sources could have changed as well during the Holocene. Moisture was delivered as today mainly during summer from the extended Caspian-Aral Basin and eastern Mediterranean, although Arctic and even North Atlantic seas might be important moisture sources when seasonal precipitation was dominated by winter precipitation. We hypothesise that warming Arctic and North Atlantic seas were important for the North Hemisphere circulation during the cold season. (C) 2014 Elsevier Ltd. All rights reserved.
Submerged macrophytes can stabilise clear water conditions in shallow lakes. However, many existing models for deep lakes neglect their impact. Here, we tested the hypothesis that submerged macrophytes can affect the water clarity in deep lakes. A one-dimensional, vertically resolved macrophyte model was developed based on PCLake and coupled to SALMO-1D and GOTM hydrophysics and validated against field data. Validation showed good coherence in dynamic growth patterns and colonisation depths. In our simulations the presence of submerged macrophytes resulted in up to 50% less phytoplankton biomass in the shallowest simulated lake (11 m) and still 15% less phytoplankton was predicted in 100 m deep oligotrophic lakes. Nutrient loading, lake depth, and lake shape had a strong influence on macrophyte effects. Nutrient competition was found to be the strongest biological interaction. Despite a number of limitations, the derived dynamic lake model suggests significant effects of submerged macrophytes on deep lake water quality. (C) 2014 Elsevier Ltd. All rights reserved.
Sediment-discharge hysteresis loops are frequently analyzed to facilitate the understanding of sediment transport processes. Hysteresis patterns, however, are often complex and their interpretation can be complicated. Particularly, quantifying hysteresis patterns remains a problematic issue. Moreover, it is currently unknown how much data is required for analyzing sediment-discharge hysteresis loops in a given area. These open questions and challenges motivated us to develop a new method for quantifying suspended-sediment hysteresis. Subsequently, we applied the new hysteresis index to three suspended-sediment and discharge datasets from a small tropical rainforest catchment. The datasets comprised a different number of events and sampling sites. Our analyses show three main findings: (1) datasets restricted to only few events, which is typical for rapid assessment surveys, were always sufficient to identify the dominating hysteresis pattern in our research area. Furthermore, some of these small datasets contained multiple-peak events that allowed identifying intra-event exhaustion effects and hence, limitations in sediment supply. (2) Datasets comprising complete hydrological years were particularly useful for analyzing seasonal dynamics of hysteresis. These analyses revealed an exhaustion of hysteresis on the inter-event scale which also points to a limited sediment supply. (3) Datasets comprising measurements from two consecutive gauges installed at the catchment outlet and on a slope within that catchment allowed analyzing the change of hysteresis patterns along the flowpath. On the slope, multiple-peak events showed a stronger intra-event exhaustion of hysteresis than at the catchment outlet. Furthermore, exhaustion of hysteresis on the inter-event scale was not evident on the slope but occurred at the catchment outlet. Our results indicate that even small sediment datasets can provide valuable insights into sediment transport processes of small catchments. Furthermore, our results may serve as a first guideline on what to expect from an analysis of hysteresis patterns for datasets of varying quality and quantity. (c) 2014 Elsevier B.V. All rights reserved.
Shallow inland water bodies in Malawi continue to be threatened by various environmental challenges despite their importance to the fisheries industry. Due to the complex interaction between natural and anthropogenic disturbances, disentangling the effect of the two may be a complicated process. The littoral zone of most water bodies is important in environmental reconstructions including pollution and lake level monitoring. This study used a littoral zone, transect-based approach employing multi-proxy palaeolimnological techniques to reconstruct recent environmental change (ca. 100 yrs.) in Lake Malombe in the Malawi Rift, East Africa. The results of the study could inform fisheries management in Lake Malombe, which experienced a catastrophic decline in fish stocks. Results support documentary evidence for the complete desiccation of the lake less than 100 years ago. Subsequently, there is evidence for accelerated eutrophication in the recent past. In light of these results, it is concluded that transect sampling approaches rather than relying on single core measurements, and the need for careful consideration of the types of proxy, are significant considerations in palaeoenvironmental reconstructions.
We assessed tropical montane cloud forest (TMCF) sensitivity to natural disturbance by drought, fire, and dieback with a 7300-year-long paleorecord. We analyzed pollen assemblages, charcoal accumulation rates, and higher plant biomarker compounds (average chain length [ACL] of n-alkanes) in sediments from Wai 'anapanapa, a small lake near the upper forest limit and the mean trade wind inversion ('IWI) in Hawai`i. The paleorecord of ACL suggests increased drought frequency and a lower awl elevation from 2555-1323 cal yr B.P. and 606-334 cal yr B.P. Charcoal began to accumulate and a novel fire regime was initiated ca. 880 cal yr B.P., followed by a decreased fire return interval at ca. 550 cal yr B.P. Diebacks occurred at 2931, 2161, 1162, and 306 cal yr B.P., and two of these were independent of drought or fire. Pollen assemblages indicate that on average species composition changed only 2.8% per decade. These dynamics, though slight, were significantly associated with disturbance. The direction of species composition change varied with disturbance type. Drought was associated with significantly more vines and lianas; fire was associated with an increase in the tree fern Sadleria and indicators of open, disturbed landscapes at the expense of epiphytic ferns; whereas stand-scale dieback was associated with an increase in the tree fern Cibotium. Though this cloud forest was dynamic in response to past disturbance, it has recovered, suggesting a resilient TMCF with no evidence of state change in vegetation type (e.g., grassland or shrubland).
Soil in a changing world is subject to both anthropogenic and environmental stresses. Soil monitoring is essential to assess the magnitude of changes in soil variables and how they affect ecosystem processes and human livelihoods. However, we cannot always be sure which sampling design is best for a given monitoring task. We employed a rotational stratified simple random sampling (rotStRS) for the estimation of temporal changes in the spatial mean of saturated hydraulic conductivity (K-s) at three sites in central Panama in 2009, 2010 and 2011. To assess this design's efficiency we compared the resulting estimates of the spatial mean and variance for 2009 with those gained from stratified simple random sampling (StRS), which was effectively the data obtained on the first sampling time, and with an equivalent unexecuted simple random sampling (SRS). The poor performance of geometrical stratification and the weak predictive relationship between measurements of successive years yielded no advantage of sampling designs more complex than SRS. The failure of stratification may be attributed to the small large-scale variability of K-s. Revisiting previously sampled locations was not beneficial because of the large small-scale variability in combination with destructive sampling, resulting in poor consistency between revisited samples. We conclude that for our K-s monitoring scheme, repeated SRS is equally effective as rotStRS. Some problems of small-scale variability might be overcome by collecting several samples at close range to reduce the effect of small-scale variation. Finally, we give recommendations on the key factors to consider when deciding whether to use stratification and rotation in a soil monitoring scheme.
Biomass allometries and coarse root biomass distribution of mountain birch in southern Iceland
(2014)
Root systems are an important pool of biomass and carbon in forest ecosystems. However, most allometric studies on forest trees focus only on the aboveground components. When estimated, root biomass has most often been calculated by using a fixed conversion factor from aboveground biomass. In order to study the size-related development of the root system of native mountain birch (Betula pubescens Ehrh. ssp. czerepanovii), we collected the coarse root system of 25 different aged birch trees (stem diameter at 50 cm length between 0.2 and 14.1 cm) and characterized them by penetration depth (< 1 m) and root thickness. Based on this dataset, allometric functions for coarse roots (> 5 mm and > 2 mm), root stock, total belowground biomass and aboveground biomass components were calculated by a nonlinear and a linear fitting approach. The study showed that coarse root biomass of mountain birch was almost exclusively (> 95 weight-%) located in the top 30 cm, even in a natural old-growth woodland. By using a cross-validation approach, we found that the nonlinear fitting procedure performed better than the linear approach with respect to predictive power. In addition, our results underscore that general assumptions of fixed conversion factors lead to an underestimation of the belowground biomass. Thus, our results provide allometric functions for a more accurate root biomass estimation to be utilized in inventory reports and ecological studies.
We have investigated the influence that megathrust earthquake slip has on the activation of splay faults using a 2-D finite element method (FEM), taking into account the effects of gravity and variations in the frictional strength properties of splay faults. We simulated both landward-dipping and seaward-dipping splay fault geometries, and imposed depth-variable slip distributions of subduction events. Our results indicate that the two types of splay fault exhibit a similar behavior, with variations in frictional properties along the faults affecting only the seismic magnitude. The triggering process is controlled by a critical depth. Megathrust slip concentrated at depths shallower than the critical depth will favor normal displacement, while megathrust slip concentrated at depths deeper than the critical depth is likely to result in reverse motion. Our results thus provide a useful tool for predicting the activation of secondary faults and may have direct implications for tsunami hazard research.
TopoToolbox is a MATLAB program for the analysis of digital elevation models (DEMs). With the release of version 2, the software adopts an object-oriented programming (OOP) approach to work with gridded DEMs and derived data such as flow directions and stream networks. The introduction of a novel technique to store flow directions as topologically ordered vectors of indices enables calculation of flow-related attributes such as flow accumulation similar to 20 times faster than conventional algorithms while at the same time reducing memory overhead to 33% of that required by the previous version. Graphical user interfaces (GUIs) enable visual exploration and interaction with DEMs and derivatives and provide access to tools targeted at fluvial and tectonic geomorphologists. With its new release, TopoToolbox has become a more memory-efficient and faster tool for basic and advanced digital terrain analysis that can be used as a framework for building hydrological and geomorphological models in MATLAB.
In the context of examining the potential usage of safe and sustainable geothermal energy in the Alberta Basin, whether in deep sediments or crystalline rock, the understanding of the in situ stress state is crucial. It is a key challenge to estimate the 3-D stress state at an arbitrarily chosen point in the crust, based on sparsely distributed in situ stress data.
To address this challenge, we present a large-scale 3-D geomechanical-numerical model (700 km x 1200 km x 80 km) from a large portion of the Alberta Basin, to provide a 3-D continuous quantification of the contemporary stress orientations and stress magnitudes. To calibrate the model, we use a large database of in situ stress orientation (321 S-Hmax) as well as stress magnitude data (981 S-V, 1720 S-hmin and 2 (+11) S-Hmax) from the Alberta Basin. To find the best-fit model, we vary the material properties and primarily the displacement boundary conditions of the model. This study focusses in detail on the statistical calibration procedure, because of the large amount of available data, the diversity of data types, and the importance of the order of data tests.
The best-fit model provides the total 3-D stress tensor for nearly the whole Alberta Basin, and allows estimation of stress orientation and stress magnitudes in advance of any well. First-order implications for the well design and configuration of enhanced geothermal systems are revealed. Systematic deviations of the modelled stress from the in situ data are found for stress orientations in the Peace River and the Bow Island Arch as well as for leak-off test magnitudes.
The Canadian database on contemporary crustal stress has not been revised systematically in the past two decades. Here we present the results of our new compilation that contains 514 new data records for the orientation data of maximum compressive horizontal stress and 188 data records that were re-assessed. In total the Canadian stress database has now 1667 data records, which is an increase of about 45%. From these data, a new Canadian Stress map as well as one for the Province of Alberta is presented.
To analyse the stress pattern, we use the quasi median on the circle as a smoothing algorithm that generates a smoothed stress map of the maximum compressive horizontal stress orientation on a regular grid. The newly introduced quasi interquartile range on the circle estimates the spreading of the data and is used as a measure for the wave-length of the stress pattern. The result of the hybrid wavelength analysis confirms that long spatial wavelength stress patterns (>= 1000 km) exist in large areas in Canada. The observed stress pattern is transmitted through the intra-plate regions.
The results reveal that shorter spatial wave length variation of the maximum compressive horizontal stress orientation of less than 200 km, prevails particularly in south-eastern and western Canada. Regional stress sources such as density contrasts, active fault systems, crustal structures, etc. might have a significant impact in these regions. In contrast to these variations, the observed stress pattern in the Alberta Basin is very homogeneous and mainly controlled by plate boundary forces and body forces. The influence of curvature of the Rocky Mountains salient in southern Alberta is minimal. The present-day horizontal stress orientations determined herein have important implications for the production of hydrocarbons and geothermal energy in the Alberta Basin. (C) 2014 Elsevier B.V. All rights reserved.
Intraplate seismicity across the Cape Verde swell: A contribution from a temporary seismic network
(2014)
We present an analysis and characterization of the regional seismicity recorded by a temporary broadband seismic network deployed in the Cape Verde archipelago between November 2007 and September 2008.
The detection of earthquakes was based on spectrograms, allowing the discrimination from low-frequency volcanic signals, resulting in 358 events of which 265 were located, the magnitudes usually being smaller than 3. For the location, a new 1-D P-velocity model was derived for the region showing a crust consistent with an oceanic crustal structure. The seismicity is located mostly offshore the westernmost and geologically youngest areas of the archipelago, near the islands of Santo Antao and Sao Vicente in the NW and Brava and Fogo in the SW. The SW cluster has a lower occurrence rate and corresponds to seismicity concentrated mainly along an alignment between Brava and the Cadamosto seamount presenting normal faulting mechanisms. The existence of the NW cluster, located offshore SW of Santo Antao, was so far unknown and concentrates around a recently recognized submarine cone field; this cluster presents focal depths extending from the crust to the upper mantle and suggests volcanic unrest No evident temporal behaviour could be perceived, although the events tend to occur in bursts of activity lasting a few days. In this recording period, no significant activity was detected at Fogo volcano, the most active volcanic edifice in Cape Verde.
The seismicity characteristics point mainly to a volcanic origin. The correlation of the recorded seismicity with active volcanic structures agrees with the tendency for a westward migration of volcanic activity in the archipelago as indicated by the geologic record. (C) 2014 Elsevier B.V. All rights reserved.
We provide new insights into the prograde evolution of HP/LT metasedimentary rocks on the basis of detailed petrologic examination, element-partitioning analysis, and thermodynamic modelling of well-preserved Fe-Mg-carpholite- and Fe-Mg-chloritoid-bearing rocks from the Afyon Zone (Anatolia). We document continuous and discontinuous compositional (ferromagnesian substitution) zoning of carpholite (overall X-Mg = 0.27-0.73) and chloritoid (overall X-Mg = 0.07-0.30), as well as clear equilibrium and disequilibrium (i.e., reaction-related) textures involving carpholite and chloritoid, which consistently account for the consistent enrichment in Mg of both minerals through time, and the progressive replacement of carpholite by chloritoid. Mg/Fe distribution coefficients calculated between carpholite and chloritoid vary widely within samples (2.2-20.0). Among this range, only values of 7-11 correlate with equilibrium textures, in agreement with data from the literature. Equilibrium phase diagrams for metapelitic compositions are calculated using a newly modified thermodynamic dataset, including most recent data for carpholite, chloritoid, chlorite, and white mica, as well as further refinements for Fe-carpholite, and both chloritoid end-members, as required to reproduce accurately petrologic observations (phase relations, experimental constraints, Mg/Fe partitioning). Modelling reveals that Mg/Fe partitioning between carpholite and chloritoid is greatly sensitive to temperature and calls for a future evaluation of possible use as a thermometer. In addition, calculations show significant effective bulk composition changes during prograde metamorphism due to the fractionation of chloritoid formed at the expense of carpholite. We retrieve P-T conditions for several carpholite and chloritoid growth stages (1) during prograde stages using unfractionated, bulk-rock XRF analyses, and (2) at peak conditions using compositions fractionated for chloritoid. The P-T paths reconstructed for the Kutahya and Afyon areas shed light on contrasting temperature conditions for these areas during prograde and peak stages.
The Himalayan mountains are dissected by some of the deepest and most impressive gorges on Earth. Constraining the interplay between river incision and rock uplift is important for understanding tectonic deformation in this region. We report here the discovery of a deeply incised canyon of the Yarlung Tsangpo River, at the eastern end of the Himalaya, which is now buried under more than 500 meters of sediments. By reconstructing the former valley bottom and dating sediments at the base of the valley fill, we show that steepening of the Tsangpo Gorge started at about 2 million to 2.5 million years ago as a consequence of an increase in rock uplift rates. The high erosion rates within the gorge are therefore a direct consequence of rapid rock uplift.
Recultivation of disturbed oil sand mining areas is an issue of increasing importance. Nevertheless only little is known about the fate of organic matter, cell abundances and microbial community structures during oil sand processing, tailings management and initial soil development on reclamation sites. Thus the focus of this work is on biogeochemical changes of mined oil sands through the entire process chain until its use as substratum for newly developing soils on reclamation sites. Therefore, oil sand, mature fine tailings (MFTs) from tailings ponds and drying cells and tailings sand covered with peat-mineral mix (PMM) as part of land reclamation were analyzed. The sample set was selected to address the question whether changes in the above-mentioned biogeochemical parameters can be related to oil sand processing or biological processes and how these changes influence microbial activities and soil development.
GC-MS analyses of oil-derived biomarkers reveal that these compounds remain unaffected by oil sand processing and biological activity. In contrast, changes in polycyclic aromatic hydrocarbon (PAH) abundance and pattern can be observed along the process chain. Especially naphthalenes, phenanthrenes and chrysenes are altered or absent on reclamation sites, Furthermore, root-bearing horizons on reclamation sites exhibit cell abundances at least ten times higher (10(8) to 10(9) cells g(-1)) than in oil sand and MFF samples (10(7) cells g(-1)) and show a higher diversity in their microbial community structure. Nitrate in the pore water and roots derived from the PMM seem to be the most important stimulants for microbial growth. The combined data show that the observed compositional changes are mostly related to biological activity and the addition of exogenous organic components (PMM), whereas oil extraction, tailings dewatering and compaction do not have significant influences on the evaluated compounds. Microbial community composition remains relatively stable through the entire process chain. (C) 2014 Elsevier B.V. All rights reserved.
We report the results of our investigations on the catchment area, lake surface sediments, and hydrology of the high altitude alpine Tso Moriri Lake, NW Himalayas (India). Our results indicate that the lake is currently alkaline, and thermally stratified with an oxic bottom layer. Results from hydrochemistry and isotopic composition (delta O-18 and delta D) of inflowing streams and lake waters show that Tso Moriri Lake is an evaporative lake with contributions from both westerly source (snow melt) and Indian summer monsoon precipitation. Geochemical and mineralogical investigations on the catchment and lake surface sediments reveal the presence of authigenic aragonite in modern lake sediment. The lithogenic components reflect the inflow and sorting processes during transport into the lake, whereas the authigenic carbonate fraction can be linked to the changes in ([precipitation+meltwater]/evaporation) (I/E) balance within the lake. The spatial variability in grain size distribution within the lake surface sediments shows that the grain size data can be utilised as a proxy for transport energy and shoreline proximity in the lake basin. We have evaluated the applicability of commonly applied environmentally sensitive proxies (isotopes, mineralogy, weathering indices) for palaeoenvironmental reconstruction in the Tso Moriri Lake. Our results show that the commonly used weathering index (Rb/Sr) is not applicable due to Sr contribution from authigenic carbonates. The useful weathering indices in Tso Moriri Lake are the Si/Al and the Chemical Proxy of Alteration (CPA). Since the carbonates are formed by evaporative processes, their presence and isotopic values can be used as indicators of I/E changes in the lake.
Rates and style of Cenozoic deformation around the Gonghe Basin, northeastern Tibetan Plateau
(2014)
The northeastern Tibetan Plateau constitutes a transitional region between the low-relief physiographic plateau to the south and the high-relief ranges of the Qilian Shan to the north. Cenozoic deformation across this margin of the plateau is associated with localized growth of fault-cored mountain ranges and associated basins. Herein, we combine detailed structural analysis of the geometry of range-bounding faults and deformation of foreland basin strata with geomorphic and exhumational records of erosion in hangingwall ranges in order to investigate the magnitude, timing, and style of deformation along the two primary fault systems, the Qinghai Nan Shan and the Gonghe Nan Shan. Structural mapping shows that both ranges have developed above imbricate fans of listric thrust faults, which sole into decollements in the middle crust. Restoration of shortening along balanced cross sections suggests a minimum of 0.8-2.2 km and 5.1-6.9 km of shortening, respectively. Growth strata in the associated foreland basin record the onset of deformation on the two fault systems at ca. 6-10 Ma and ca. 7-10 Ma, respectively, and thus our analysis suggests late Cenozoic shortening rates of 0.2 +0.2/-0.1 km/m.y. and 0.7 +0.3/-0.2 km/m.y. along the north and south sides of Gonghe Basin. Along the Qinghai Nan Shan, these rates are similar to late Pleistocene slip rates of similar to 0.10 +/- 0.04 mm/yr, derived from restoration and dating of a deformed alluvial-fan surface. Collectively, our results imply that deformation along both flanks of the doubly vergent Qilian Shan-Nan Shan initiated by ca. 10 Ma and that subsequent shortening has been relatively steady since that time.
Due to increasing demands and competition for high quality groundwater resources in many parts of the world, there is an urgent need for efficient methods that shed light on the interplay between complex natural settings and anthropogenic impacts. Thus a new approach is introduced, that aims to identify and quantify the predominant processes or factors of influence that drive groundwater and lake water dynamics on a catchment scale. The approach involves a non-linear dimension reduction method called Isometric feature mapping (Isomap). This method is applied to time series of groundwater head and lake water level data from a complex geological setting in Northeastern Germany. Two factors explaining more than 95% of the observed spatial variations are identified: (1) the anthropogenic impact of a waterworks in the study area and (2) natural groundwater recharge with different degrees of dampening at the respective sites of observation. The approach enables a presumption-free assessment to be made of the existing geological conception in the catchment, leading to an extension of the conception. Previously unknown hydraulic connections between two aquifers are identified, and connections revealed between surface water bodies and groundwater. (C) 2014 Elsevier B.V. All rights reserved.
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.
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.
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.
Eolian dust is a significant source of iron and other nutrients that are essential for the health of marine ecosystems and potentially a controlling factor of the high nutrient-low chlorophyll status of the Subarctic North Pacific. We map the spatial distribution of dust input using three different geochemical tracers of eolian dust, He-4, Th-232 and rare earth elements, in combination with grain size distribution data, from a set of core-top sediments covering the entire Subarctic North Pacific. Using the suite of geochemical proxies to fingerprint different lithogenic components, we deconvolve eolian dust input from other lithogenic inputs such as volcanic ash, ice-rafted debris, riverine and hemipelagic input. While the open ocean sites far away from the volcanic arcs are dominantly composed of pure eolian dust, lithogenic components other than eolian dust play a more crucial role along the arcs. In sites dominated by dust, eolian dust input appears to be characterized by a nearly uniform grain size mode at similar to 4 mu m.
Applying the Th-230-normalization technique, our proxies yield a consistent pattern of uniform dust fluxes of 1-2 g/m(2)/yr across the Subarctic North Pacific. Elevated eolian dust fluxes of 2-4 g/m(2)/yr characterize the westernmost region off Japan and the southern Kurile Islands south of 45 degrees N and west of 165 degrees E along the main pathway of the westerly winds. The core-top based dust flux reconstruction is consistent with recent estimates based on dissolved thorium isotope concentrations in seawater from the Subarctic North Pacific. The dust flux pattern compares well with state-of-the-art dust model predictions in the western and central Subarctic North Pacific, but we find that dust fluxes are higher than modeled fluxes by 0.5-1 g/m(2)/yr in the northwest, northeast and eastern Subarctic North Pacific. Our results provide an important benchmark for biogeochemical models and a robust approach for downcore studies testing dust-induced iron fertilization of past changes in biological productivity in the Subarctic North Pacific.
The general warming trend of the last deglaciation was interrupted by the Younger Dryas, a period of abrupt cooling and widespread environmental change(1-10). Ice core records suggest the abrupt cooling began 12,846 years ago in Greenland(10), about 170 years before the significant environmental and vegetation change in western Europe(7) classically defined as the Younger Dryas. However, this difference in timing falls within age model uncertainties. Here we use the hydrogen isotope composition of lipid biomarkers from precisely dated varved sediments from Lake Meerfelder Maar to reconstruct hydroclimate over western Europe. We observe a decrease in the hydrogen isotope values of both aquatic and terrestrial lipids 12,850 years ago, indicating cooling climate in this region synchronous with the abrupt cooling in Greenland. A second drop occurs 170 years later, mainly in the hydrogen isotope record of aquatic lipids but to a lesser extent in the terrestrial lipids, which we attribute to aridification, as well as a change in moisture source and pathway. We thus confirm that there was indeed a lag between cooling and substantial hydrologic and environmental change in western Europe. We suggest the delay is related to the expansion of sea ice in the North Atlantic Ocean and the subsequent southward migration of the westerly wind system(9). We further suggest that these hydrological changes amplified environmental change in western Europe at the onset of the Younger Dryas.
The Shell Bar in the Qaidam Basin, China, is a prominent geological feature composed of millions of densely packed Corbicula shells. Since the mid 1980s, it has been regarded as evidence for existence of a large lake during Marine Isotope Stage (MIS) 3 in the presently hyper-arid Qaidam Basin. Early studies suggested the bivalve shells accumulated at the shore of a large lake, whereas more recent work led to the conclusion that the Shell Bar was formed within a deeper water body. Based on our re-assessment of sediments and fossils from the Shell Bar, investigation of exposed fluvio-lacustrine sections upstream of the Shell Bar and study of nearby modern streams, we infer that the Shell Bar represents a stream deposit. Corbicula is a typical stream-dweller around the world. Preservation of Corbicula shells of different sizes, as well as occurrence of many articulated shells, provide evidence against post-mortem transport and accumulation along a lake shore. Additionally, the SE-NW alignment of the Shell Bar is similar to modern intermittent stream beds in its vicinity and corresponds to the present-day slope towards the basin centre further NW, and furthermore, the predominantly sandy sediments also indicate that the Shell Bar was formed in a stream. Abundant ostracod shells in the Shell Bar sediments originated from stream-dwelling species that are abundant in modern streams in the vicinity of the Shell Bar, or in part from fluvio-lacustrine sediments exposed upstream of the Shell Bar, as a result of erosion and re-deposition. Deflation of alluvial fine-grained sediments in the Shell Bar region and protection of the stream deposits by the large and thick-walled Corbicula shells reversed the former channel relief and yielded the modern exposure, which is a prominent morphological feature. Occurrence of Corbicula shells in the Qaidam Basin indicates climate was apparently warmer than present during the formation of the Shell Bar because Corbicula does not live at similar or higher altitudes in the region today. Because the Shell Bar is no longer considered a deposit formed within a lake, its presence does not indicate paleoclimate conditions wetter than today.
A geological feature in the Qaidam Basin known as the "Shell Bar" contains millions of freshwater clam shells buried in situ. Since the 1980s, this feature in the now hyper-arid basin has been interpreted to be lake deposits that provide evidence for a warmer and more humid climate than present during late marine isotope stage 3 (MIS 3). Global climate during late MIS 3 and the last glacial maximum, however, was cold and dry, with much lower sea levels. We re-investigated the feature geomorphologically and sedimentologically, and employed optically stimulated luminescence (OSL) dating to verify the chronology of the sediments. We interpret the Shell Bar to be a remnant of a river channel formed by a stream that ran across an exposed lake bed during a regressive lake phase. Deflation of the surrounding older, fine-grained lacustrine deposits has left the fluvial channel sediments topographically inverted, indicating the erosive nature of the landscape. Luminescence ages place the formation of the Shell Bar in MIS 5 (similar to 113-99 ka), much older than previous radiocarbon ages of < 40 ka BP, but place the paleoclimatic inferences more in accord with other regional and global climate proxy records. We present a brief review of the age differences derived from C-14 and OSL dating of some critical sections that were thought to represent a warmer and more humid climate than present during late MIS 3. We attribute the differences to underestimation of C-14 ages. We suggest that C-14 ages older than similar to 25 ka BP may require re-investigation, especially dates on samples from arid regions.
Advance in geocomputation
(2014)
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.
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.
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.
Aim To understand the role and significance of the reindeer, Rangifer tarandus (Linnaeus, 1758), as a specific indicator in terms of late Quaternary biogeography and to determine the effects of global climate change on its range and local extinction dynamics at the end of the Ice Age.
Location Late Pleistocene/early Holocene range of reindeer over all of central and western Europe, including southern Scandinavia and northern Iberia, but excluding Russia, Belarus and the Ukraine.
Methods Radiocarbon-dated subfossil records of R. tarandus from both archaeological and natural deposits younger than 25,000 years were assembled in a database. The distribution area was divided into six representative regions. The C-14 dates were calibrated and plotted chronologically in maps in order to compare presence and absence and regional extinction patterns from one region to another.
Main conclusions The late Quaternary record for reindeer in Europe during the last 25 kyr shows a climate-driven dispersal and retreat in response to climate change, with regional variations. The collapse of the mammoth steppe biome did not lead to the local extinction in Europe, as in the case of other megafaunal species. Rangifer tarandus co-existed for about 3000 years during the Late Glacial and early Holocene with typical temperate species such as red deer and roe deer in non-analogue faunal communities. The regional extinction at the end of the Pleistocene coincides with the transition from light open birch/pine forests to pine/deciduous forests.
Velocity models are essential to process two-and three-dimensional ground-penetrating radar (GPR) data. Furthermore, velocity information aids the interpretation of such data sets because velocity variations reflect important material properties such as water content. In many GPR applications, common midpoint (CMP) surveys are routinely collected to determine one-dimensional velocity models at selected locations. To analyse CMP data gathers, spectral velocity analyses relying on the normal-moveout (NMO) model are commonly employed. Using Dix's formula, the derived NMO velocities can be further converted to interval velocities which are needed for processing and interpretation. Because of the inherent assumptions and limitations of such approaches, we investigate and propose an alternative procedure based on the global inversion of reflection travel-times. We use a finite-difference solver of the Eikonal equation to accurately solve the forward problem in combination with particle swarm optimization (PSO) to find one-dimensional GPR velocity models explaining our data. Because PSO is a robust and efficient global optimization tool, our inversion approach includes generating an ensemble of representative solutions that allows us to analyse uncertainties in the model space. Using synthetic data examples, we test and evaluate our inversion approach to analyse CMP data collected across typical near-surface environments. Application to a field data set recorded at a well-constrained test site including a comparison to independent borehole and direct-push data, further illustrates the potential of the proposed approach, which includes a straightforward and understandable appraisal of non-uniqueness and uncertainty issues, respectively. We conclude that our methodology is a feasible and powerful tool to analyse GPR CMP data and allows practitioners and researchers to evaluate the reliability of CMP derived velocity models.
Inferring a ground-motion prediction equation (GMPE) for a region in which only a small number of seismic events has been observed is a challenging task. A response to this data scarcity is to utilise data from other regions in the hope that there exist common patterns in the generation of ground motion that can contribute to the development of a GMPE for the region in question. This is not an unreasonable course of action since we expect regional GMPEs to be related to each other. In this work we model this relatedness by assuming that the regional GMPEs occupy a common low-dimensional manifold in the space of all possible GMPEs. As a consequence, the GMPEs are fitted in a joint manner and not independent of each other, borrowing predictive strength from each other's regional datasets. Experimentation on a real dataset shows that the manifold assumption displays better predictive performance over fitting regional GMPEs independent of each other. (C) 2014 Elsevier Ltd. All rights reserved.
We report delta O-18 and delta C-13 values of 21 fossil shells from the aquatic gastropod Radix from a sediment core taken in the eastern basin of Lake Karakul, Tajikistan (38.86-39.16A degrees N, 73.26-73.56A degrees E, 3,928 m above sea level) and covering the last 4,200 cal yr BP. The lake is surrounded by many palaeoshorelines evidencing former lake-level changes, most likely triggered by changes in meltwater flux. This hypothesis was tested by interpreting the isotope ratios of Radix shells together with delta O-18 values of Ostracoda and of authigenic aragonite. The mean delta O-18 values of Radix and Ostracoda fall along the same long-term trend indicating a change in the isotopic composition of precipitation, which contributed to the glaciers in the catchment as snow and finally as melt water to the lake. The sclerochronological delta O-18 and delta C-13 patterns in Radix shells provide seasonal weather information, which is discussed in context with previously proposed climatic changes during the last 4,200 cal yr BP. The period between similar to 4,200 and 3,000 cal yr BP was characterized by stepwise glacier advance in the catchment most likely due to a precipitation surplus. Subsequently the climate remained relatively cold but the lake level fluctuated, as indicated by ostracod shell isotope data. From similar to 1,800 cal yr BP the sclerochronological patterns provide evidence for increasing melt water flux and transport of allochthonous carbon into the lake, most likely due to an accelerated glacier retreat. The period around 1,500 cal yr BP was characterized by strong warming, increasing meltwater flux, glacier retreat and an increasing lake level. Warm conditions continued until similar to 500 cal yr I'P probably representing the end of the Medieval Warm Period. A short relatively cold (dry?) period and a lower lake level are assumed for similar to 350 cal yr BP, possibly an analogue to the Maunder Minimum cooling in the North Atlantic region. Our results show that the lake system is complex, and that changes were triggered by external forcing and feedbacks. The similarity of delta O-18 values in Radix and ostracod shells demonstrates that both archives provide complementary information.
Leaching of dissolved C in arable hummocky ground moraine soil landscapes is characterized by a spatial continuum of more or less erosion-affected Luvisols, Calcaric Regosols at exposed positions, and Colluvic Regosols in depressions. Our objective was to estimate the fluxes of dissolved C in four differently eroded soils as affected by erosion-induced pedological and soil structural alterations. In this model study, we considered landscape position effects by adapting the water table as the bottom boundary condition and erosion effects by using pedon-specific soil hydraulic properties. The one-dimensional vertical water movement was described with the Richards equation using HYDRUS-1D. Solute fluxes were obtained by combining calculated water fluxes with concentrations of dissolved organic and inorganic C (DOC and DIC, respectively) measured from soil solution extracted by suction cups at biweekly intervals. In the 3-yr period (2010-2012), DOC fluxes in the 2-m soil depth were similar at the three non-colluvic locations with -0.8 +/- 0.1 g m(-2) yr(-1) (i.e., outflow) but were 0.4 g m(-2) yr(-1) (i.e., input) in the depression. The DIC fluxes ranged from -10.2 g m(-2) yr(-1) for the eroded Luvisol, -9.2 g m(-2) yr(-1) for the Luvisol, and -6.1 g m(-2) yr(-1) for the Calcaric Regosol to 3.2 g m(-2) yr(-1) for the Colluvic Regosol. The temporal variations in DOC and DIC fluxes were controlled by water fluxes. The spatially distributed leaching results corroborate the hypothesis that the effects of soil erosion influence fluxes through modified hydraulic and transport properties and terrain-dependent boundary conditions.
Vertical radar profiling (VRP) is a single-borehole geophysical technique, in which the receiver antenna is located within a borehole and the transmitter antenna is placed at one or various offsets from the borehole. Today, VRP surveying is primarily used to derive 1D velocity models by inverting the arrival times of direct waves. Using field data collected at a well-constrained test site in Germany, we evaluated a VRP workflow relying on the analysis of direct-arrival traveltimes and amplitudes as well as on imaging reflection events. To invert our VRP traveltime data, we used a global inversion strategy resulting in an ensemble of acceptable velocity models, and thus, it allowed us to appraise uncertainty issues in the estimated velocities as well as in porosity models derived via petrophysical translations. In addition to traveltime inversion, the analysis of direct-wave amplitudes and reflection events provided further valuable information regarding subsurface properties and architecture. The used VRP amplitude preprocessing and inversion procedures were adapted from raybased crosshole ground-penetrating radar (GPR) attenuation tomography and resulted in an attenuation model, which can be used to estimate variations in electrical resistivity. Our VRP reflection imaging approach relied on corridor stacking, which is a well-established processing sequence in vertical seismic profiling. The resulting reflection image outlines bounding layers and can be directly compared to surface-based GPR reflection profiling. Our results of the combined analysis of VRP, traveltimes, amplitudes, and reflections were consistent with independent core and borehole logs as well as GPR reflection profiles, which enabled us to derive a detailed hydro-stratigraphic model as needed, for example, to understand and model groundwater flow and transport.
Geophysical techniques offer the potential to tomographically image physical parameter variations in the ground in two or three dimensions. Due to the limited number and accuracy of the recorded data, geophysical model generation by inversion suffers ambiguity. Linking the model generation process of disparate data by jointly inverting two or more data sets allows for improved model reconstruction. Fully nonlinear inversion using optimization techniques searching the solution space of the inverse problem globally enables quantitative assessment of the ambiguity inherent to the model reconstruction. We used two different multiobjective particle swarm optimization approaches to jointly invert synthetic crosshole tomographic data sets comprising radar and P-wave traveltimes, respectively. Beginning with a nonlinear joint inversion founded on the principle of Pareto optimality and game theoretic concepts, we obtained a set of Pareto-optimal solutions comprising commonly structured radar and P-wave velocity models for low computational costs. However, the efficiency of the approach goes along with some risk of achieving a final model ensemble not adequately illustrating the ambiguity inherent to the model reconstruction process. Taking advantage of the results of the first approach, we inverted the database using a different nonlinear joint-inversion approach reducing the multiobjective optimization problem to a single-objective one. Computational costs were significantly higher, but the final models were obtained mutually independently allowing for objective appraisal of model parameter determination. Despite the high computational effort, the approach was found to be an efficient nonlinear joint-inversion formulation compared to what could be extracted from individual nonlinear inversions of both data sets.
The goal of China's sloping land conversion programme (SLCP) is to combat soil erosion and to reduce rural poverty. An ex-ante assessment of possible SLCP impacts was conducted with a focus on rural sustainability, taking the drought-prone region of Guyuan in Western China as an example. The Framework for Participatory Impact Assessment (FoPIA) was used to conduct two complementary impact assessments, one assessing SLCP impacts at regional level and a second one assessing alternative forest management options, to explore possible trade-offs among the economic, social and environmental dimensions of sustainability. Regional stakeholders assessed the SLCP to be capable of reducing soil erosion but felt it negatively affected rural employment, and a further continuation of the Programme was advocated. Assessment of three forest management scenarios by scientists showed that an orientation towards energy forests is potentially beneficial to all three sustainability dimensions. Ecological forests had disproportionate positive impacts on environmental functions and adverse impact on the other two sustainability dimensions. Economic forests were assessed to serve primarily the economic and social sustainability dimensions, while environmental impacts were still tolerable. The FoPIA results were evaluated against the available literature on the SLCP. Overall, the assessment results appeared to be reasonable, but the results of the regional stakeholders appeared to be too optimistic compared with the more critical assessment of the scientists. The SLCP seems to have the potential to tackle soil erosion but requires integrated forest management to minimize the risk of water stress while contributing to economic and social benefits in Guyuan. Copyright (C) 2012 John Wiley & Sons, Ltd.
On 11 May 2011, a M (w) = 5.1 earthquake shook the town of Lorca (SE Spain) causing a disproportionately large damage for its magnitude. In order to contribute to knowledge of the behavior of the active faults present in the region and define the parameters which control their motion, we made a detailed study of the rupture process of this earthquake from inversion of body waves at regional and teleseismic distances. Ground motion displacements obtained in this way are in agreement with near-field strong motion data and GPS observations recorded in Lorca. We have obtained a partly bilateral rupture propagating to WSW (238A degrees, 54A degrees, 59A degrees) with 27 cm of maximum slip and shallow focus (4 km). The fault plane orientation corresponds to that of the Cejo de los Enamorados Fault located NE of the Lorca town and parallel to the Alhama de Murcia Fault. The distribution of slip on the fault plane can explain the lack of any observed surface rupture as we found that the rupture started at 4-km depth along a plane dipping at 54A degrees, with motion propagating upward to stop at 1.5 km below the surface. The strong motion and GPS data recorded near the epicenter are in agreement with the maximum slip on the fault. Directivity effects and the extreme shallowness of the rupture could explain the considerable damage that the earthquake caused in the town of Lorca.
Scenario analysis is a widely used approach to incorporate uncertainties in global change research. In the context of regional ecosystem service and landscape management where global IPCC climate simulations and their downscaled derivates are applied, it can be useful to work with regional sodo-economic scenarios that are coherent with the global IPCC scenarios. The consistency with the original source scenarios, transparency and reproducibility of the methods used as well as the internal consistency of the derived scenarios are important methodological prerequisites for coherently downscaling pre-existing source scenarios. In contrast to well-established systematic-qualitative scenario techniques, we employ here a formal technique of scenario construction which combines expert judgement with a quantitative, indicator-based selection algorithm in order to deduce a formally consistent set of focus scenario. In our case study, these focus scenarios reflect the potential development pathways of major national-level drivers for ecosystem service management in Swiss mountain regions. The integration of an extra impact factor ("Global Trends") directly referring to the four principle SRES scenario families, helped us to formally internalise base assumptions of IPCC SRES scenarios to regional scenarios that address a different thematic focus (ecosystem service management), spatial level (national) and time horizon (2050). Compared to the well-established systematic-qualitative approach, we find strong similarities between the two methods, including the susceptibility to personal judgement which is only partly reduced by the formal method. However, the formalised scenario approach conveys four clear advantages, (1) the better documentation of the process, (2) its reproducibility, (3) the openness in terms of the number and directions of the finally selected set of scenarios, and (4) its analytical power. (C) 2013 Elsevier Inc. All rights reserved.
The age models of fluvio-lacustrine sedimentary sequences are often subject of discussions in paleoclimate research. The techniques employed to build an age model are very diverse, ranging from visual or intuitive estimation of the age-depth relationship over linear or spline interpolations between age control points to sophisticated Bayesian techniques also taking into account the most likely deposition times of the type of sediment within the sequence. All these methods, however, fail in detecting abrupt variations in sedimentation rates, including the possibility of episodes of no deposition (hiatus), which is the strength of the method presented in this work. The new technique simply compares the deposition time of equally thick sediment slices from the differences of subsequent radiometric age dates and the unit deposition times of the various sediment types. The percentage overlap of the distributions of these two sources of information, together with the evidence from the sedimentary record, helps to build an age model of complex sequences including abrupt variations in the rate of deposition including one or many hiatuses. (C) 2014 Elsevier B.V. All rights reserved.
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.
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.
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.
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.
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.
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.
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.