Refine
Has Fulltext
- no (140)
Year of publication
- 2013 (140) (remove)
Document Type
- Article (140) (remove)
Language
- English (140)
Is part of the Bibliography
- yes (140)
Keywords
- Climate change (4)
- Land use change (4)
- Iran (3)
- exhumation (3)
- Alborz Mountains (2)
- Asian monsoon (2)
- Earthquake source observations (2)
- Geostatistics (2)
- Holocene (2)
- Mountain basins (2)
Institute
- Institut für Geowissenschaften (140) (remove)
Astronomically tuned cyclic sedimentary successions provide unprecedented insight into the temporal evolution of depositional systems and major geologic events. However, placing astronomically calibrated records into an absolute time frame with confidence requires independent and precise geochronologic constraints. Astronomical tuning of the precessionally modulated sedimentary cycles of the Mediterranean Basin deposited during the Messinian Salinity Crisis (5.96-5.33 Ma) has indicated an similar to 90 k.y. "Messinian gap", corresponding to the evaporative drawdown of the Mediterranean following the closure of the Mediterranean-Atlantic gateway. In the Messinian deposits, a volcanic ash dated by Ar-40/Ar-39 geochronology was used to anchor the sedimentary cycles to the insolation curve. However, the uncertainty of the Ar-40/Ar-39 date introduces a potential two-cycle (similar to 40 k.y.) uncertainty in the tuning. Using high-precision chemical abrasion-thermal ionization mass spectrometry (CA-TIMS) U-Pb geochronology on single zircon grains from two Messinian ash layers in Italy, we obtained dates of 5.5320 +/- 0.0046 Ma and 5.5320 +/- 0.0074 Ma with sub-precessional resolution. Combined with our astronomical tuning of the Messinian Lower Evaporites, the results refine the duration of the "Messinian gap" to at most 28 or 58 +/- 9.6 k.y., which correlates with either the TG12 glacial interval alone, or both TG12 and TG14 glacial intervals, supporting the hypothesis of a glacio-eustatic contribution in fully isolating the Mediterranean from the Atlantic Ocean. Our new U-Pb dates also allow us to infer a precessionally modulated cyclicity for the post-evaporitic deposits, and hence enable us to tune those successions to the insolation curve.
Improving our understanding of biodiversity and ecosystem functioning and our capacity to inform ecosystem management requires an integrated framework for functional biodiversity research (FBR). However, adequate integration among empirical approaches (monitoring and experimental) and modelling has rarely been achieved in FBR. We offer an appraisal of the issues involved and chart a course towards enhanced integration. A major element of this path is the joint orientation towards the continuous refinement of a theoretical framework for FBR that links theory testing and generalization with applied research oriented towards the conservation of biodiversity and ecosystem functioning. We further emphasize existing decision-making frameworks as suitable instruments to practically merge these different aims of FBR and bring them into application. This integrated framework requires joint research planning, and should improve communication and stimulate collaboration between modellers and empiricists, thereby overcoming existing reservations and prejudices. The implementation of this integrative research agenda for FBR requires an adaptation in most national and international funding schemes in order to accommodate such joint teams and their more complex structures and data needs.
While sea level rise is one of the most likely consequences of climate change, the provoked costs remain highly uncertain. Based on a block-maxima approach, we provide a stochastic framework to estimate the increase of expected damages with sea level rise as well as with meteorological changes and demonstrate the application to two case studies. In addition, the uncertainty of the damage estimations due to the stochastic nature of extreme events is studied. Starting with the probability distribution of extreme flood levels, we calculate the distribution of implied damages in a specific region employing stage-damage functions. Universal relations of the expected damages and their standard deviation, which demonstrate the importance of the shape of the damage function, are provided. We also calculate how flood protection reduces the damages leading to a more complex picture, where the extreme value behavior plays a fundamental role. Citation: Boettle, M., D. Rybski, and J. P. Kropp (2013), How changing sea level extremes and protection measures alter coastal flood damages, Water Resour. Res., 49, 1199-1210, doi: 10.1002/wrcr.20108.
The Tibetan Plateau (TP), including its surrounding mountain ranges, represents the largest store of ice outside the polar regions. It hosts numerous lakes as well as the head waters of major Asian rivers, on which billions of people depend, and it is particularly sensitive to climate change. The moisture transport to the TP is controlled by the Indian and Pacific monsoon and the Westerlies. Understanding the evolution of the interaction of these circulation systems requires studies on climate archives in different spatial and temporal contexts. The objective of this study is to learn more about the interannual variability of precipitation patterns across the TP and how different hydrologic systems react to different climatic factors.
Aragonite shells of the aquatic gastropod Radix, which is widely distributed in the region, may represent suitable archives for inferring hydrologic and climatic signals in particularly high resolution. Therefore, sclerochronological studies of delta O-18 and delta C-13 ratios in Radix shells from seven lakes were conducted, each representing a different hydrologic and climatic setting, on a transect from the Pamirs across the TP.
The shell patterns exhibit an increasing influence of precipitation and a decreasing influence of evaporation on the isotope compositions from west to east. delta O-18 values of shells from lakes on the eastern and central TP (Donggi Cona, Yamdrok Yumco, Tarab Co) mirror monsoon signals, indicated by more negative values and higher variabilities compared to the more western lakes (Karakul, Bangong/Nyak, Manasarovar). In Yadang Co, located on the central southern TP, the monsoon rains did not reach the lake in the sampling year, although it is located in a region which is usually affected by monsoon circulation. The delta O-18 values are used to differentiate the annual hydrological cycle into ice cover period, melt water period, precipitation period and evaporation period. delta C-13 compositions in the shells particularly depend on specific habitats, which vary in biological productivity and in carbon sources. delta O-18 and delta C-13 patterns show a positive covariance in shells originating from large closed basins. The results show that Radix shells mirror general climatic differences between the seven lake regions. These differences reflect both regional and local climate signals in sub-seasonal resolution, without noticeable dependence on the particular lake system.
Scarcity of hydrological data, especially streamflow discharge and groundwater level series, restricts the understanding of channel transmission losses (TL) in drylands. Furthermore, the lack of information on spatial river dynamics encompasses high uncertainty on TL analysis in large rivers. The objective of this study was to combine the information from streamflow and groundwater level series with multi-temporal satellite data to derive a hydrological concept of TL for a reach of the Middle Jaguaribe River (MJR) in semi-arid north-eastern Brazil. Based on this analysis, we proposed strategies for its modelling and simulation. TL take place in an alluvium, where river and groundwater can be considered to be hydraulically connected. Most losses certainly infiltrated only through streambed and levees and not through the flood plains, as could be shown by satellite image analysis. TL events whose input river flows were smaller than a threshold did not reach the outlet of the MJR. TL events whose input flows were higher than this threshold reached the outlet losing on average 30% of their input. During the dry seasons (DS) and at the beginning of rainy seasons (DS/BRS), no river flow is expected for pre-events, and events have vertical infiltration into the alluvium. At the middle and the end of the rainy seasons (MRS/ERS), river flow sustained by base flow occurs before/after events, and lateral infiltration into the alluvium plays a major role. Thus, the MJR shifts from being a losing river at DS/BRS to become a losing/gaining (mostly losing) river at MRS/ERS. A model of this system has to include the coupling of river and groundwater flow processes linked by a leakage approach.
In this study, we evaluated the potential of the hydrogen isotopic composition of algal lipid biomarkers as a proxy for past hydroclimatic variability in hypersaline Isabel Lake, Mexico (Eastern Pacific). We compared rainfall variability recorded in the region over the last 65 years with changes in delta D values of the most abundant compounds preserved in the uppermost 16 cm of lake sediment. Changes in delta D values of the 1,15-C-32 diol (delta D-diol), a specific biomarker of algal populations, were related to rainfall variability; specifically, n-alkyl diols were more deuterium-enriched (depleted) during wetter (drier) periods. Strikingly, neither the magnitude of lipid biomarker isotopic changes over interannual timescales (of up to 70-80 parts per thousand) nor the direction of that variability can be explained by changes in delta D values of the water source or salinity fluctuations (approximately 30 on the practical salinity scale) controlled by seasonal rainfall. However, changes in sedimentary biomarker composition, higher total organic carbon content and less negative delta C-13 values of the 1,15-C-32 diol indicate enhanced algal growth during wetter periods. We find that these conditions result in less negative delta D values of n-alkyl diols. We hypothesize that due to higher lipid demand during enhanced algal growth, an increasing proportion of hydrogen for lipid synthesis is derived from the cytosol via oxidation of polysaccharides, which may cause a deuterium enrichment of the acetogenic compounds. This study has significant implications for paleohydrological reconstructions using algal lipid delta D values, particularly in highly seasonal environments such as Isabel Lake. In such environments, delta D values of specific algal lipid biomarkers may not record the full seasonal cycle in source water delta D but appear to be mainly controlled by the physiological state of algal populations. Our data provide the first evidence that changes in D/H fractionation due to algal growth conditions can be recorded in sediments. For paleoclimate reconstructions in such environments, algal growth conditions should be constrained with additional proxy data (delta C-13 values of the same biomarkers), as the net D/H fractionation between water and lipids may not have been constant over time.
In many hydrological applications, ground-wave velocity measurements are increasingly used to map and monitor shallow soil water content. In this study, we propose an automated spectral velocity analysis method to determine the direct ground-wave (DGW) velocity from common midpoint (CMP) or multi-offset ground-penetrating radar (GPR) data. The method introduced in this paper is a variation of the well-known spectral velocity analysis for seismic and GPR reflection events where velocity spectra are computed using different coherency measures along hyperbolas following the normal moveout model. Here, the unnormalized cross-correlation is computed between waveforms across data gathers that are corrected with a linear moveout equation using a predefined range of velocities. Peaks in the resulting velocity spectra identify linear events in the GPR data gathers like DGW events and allow for estimating the corresponding velocities. In addition to obtaining a DGW velocity measurement, we propose a robust method to estimate the associated velocity uncertainties based on the width of the peak in the calculated velocity spectrum. Our proposed method is tested on synthetic data examples to evaluate the influence of subsurface velocity, surveying geometry and signal frequency on the accuracy of estimated ground-wave velocities. In addition, we investigate the influence of such velocity uncertainties on subsequent soil water content estimates using an established petrophysical relationship. Furthermore, we apply our approach to analyse field data, which were collected across a test site in Canada to monitor a wide range of seasonal soil moisture variations. A comparison between our spectral velocity estimates and results derived from manually picked ground-wave arrivals shows good agreement, which illustrates that our spectral velocity analysis is a feasible tool to analyse DGW arrivals in multi-offset GPR data gathers in an objective and more automated manner.
Beta diversity is a conceptual link between diversity at local and regional scales. Various additional methodologies of quantifying this and related phenomena have been applied. Among them, measures of pairwise (dis)similarity of sites are particularly popular. Undersampling, i.e. not recording all taxa present at a site, is a common situation in ecological data. Bias in many metrics related to beta diversity must be expected, but only few studies have explicitly investigated the properties of various measures under undersampling conditions. On the basis of an empirical data set, representing near-complete local inventories of the Lepidoptera from an isolated Pacific island, as well as simulated communities with varying properties, we mimicked different levels of undersampling. We used 14 different approaches to quantify beta diversity, among them dataset-wide multiplicative partitioning (i.e. true beta diversity') and pairwise site x site dissimilarities. We compared their values from incomplete samples to true results from the full data. We used these comparisons to quantify undersampling bias and we calculated correlations of the dissimilarity measures of undersampled data with complete data of sites. Almost all tested metrics showed bias and low correlations under moderate to severe undersampling conditions (as well as deteriorating precision, i.e. large chance effects on results). Measures that used only species incidence were very sensitive to undersampling, while abundance-based metrics with high dependency on the distribution of the most common taxa were particularly robust. Simulated data showed sensitivity of results to the abundance distribution, confirming that data sets of high evenness and/or the application of metrics that are strongly affected by rare species are particularly sensitive to undersampling. The class of beta measure to be used should depend on the research question being asked as different metrics can lead to quite different conclusions even without undersampling effects. For each class of metric, there is a trade-off between robustness to undersampling and sensitivity to rare species. In consequence, using incidence-based metrics carries a particular risk of false conclusions when undersampled data are involved. Developing bias corrections for such metrics would be desirable.
Soil moisture estimation under low vegetation cover using a multi-angular polarimetric decomposition
(2013)
The estimation of volumetric soil moisture under low agricultural vegetation from fully polarimetric synthetic aperture radar (SAR) data at L-band using a multi-angular polarimetric decomposition is investigated. Radar polarimetry provides the framework to decompose the backscattered signal into different canonical scattering mechanisms referring to scattering contributions from the underlying soil and the vegetation cover. Multiangular observation diversity further increases the information space for soil moisture inversion enabling higher inversion rates and a stable inversion performance. The developed approach was applied on the multi-angular L-band data set acquired by German Aerospace Center's ESAR sensor as part of the OPAQUE campaign in 2008. The obtained results are compared against ground measurements collected by the OPAQUE team over a variety of vegetated agricultural fields. The validation of the estimated against ground measured soil moisture results in an root mean square error level of 6-8 vol.% including all test fields with a variety of crop types.
Facies bodies in geostatistical models of deep-water depositional environments generally represent channel-levee-overbank-lobe morphologies. Such models adequately capture one set of the erosional and depositional processes resulting from turbidity currents traveling downslope to the ocean basin floor. However, depositional morphologies diverge from the straight forward channel-levee-overbank-lobe paradigm when the topography of the slope or the shape of the basin impacts the timing and magnitude of turbidity current deposition. Subaqueous mass-transport-deposits (MTDs) present the need for an exception to the channel-levee-overbank-lobe archetype. Irregular surface topography of subaqueous MTDs can play a primary role in controlling sand deposition from turbidity currents. MTD topography creates mini-basins in which sand accumulates in irregularly-shaped deposits. These accumulations are difficult to laterally correlate using well-log data due to their variable and unpredictable shape and size. Prediction is further complicated because sandstone bodies typical of this setting are difficult to resolve in seismic-reflection data. An event-based model is presented, called DFTopoSim, which simulates debris flows and turbidity currents. The accommodation space on top of and between debris flow lobes is filled in by sand from turbidity currents. When applied to a subsurface case in the Molasse Basin of Upper Austria, DFTopoSim predicts sand packages consistent with observations from core, well, and seismic data and the interpretation of the sedimentologic processes. DFTopoSim expands the set of available geostatistical deep-water depositional models beyond the standard channel-levee-overbank-lobe model.
In the Western Alps, the Piemont-Ligurian oceanic domain records blueschist to eclogite metamorphic conditions during the Alpine orogeny. This domain is classically divided into two "zones" (Combin and Zermatt-Saas), with contrasting metamorphic evolution, and separated tectonically by the Combin fault. This study presents new metamorphic and temperature (RSCM thermometry) data obtained in Piemont-Ligurian metasediments and proposes a reevaluation of the P-T evolution of this domain. In the upper unit (or "Combin zone") temperatures are in the range of 420-530 A degrees C, with an increase of temperature from upper to lower structural levels. Petrological evidences show that these temperatures are related to the retrograde path and to deformation at greenschist metamorphic conditions. This highlights heating during exhumation of HP metamorphic rocks. In the lower unit (or "Zermatt-Saas zone"), temperatures are very homogeneous in the range of 500-540 A degrees C. This shows almost continuous downward temperature increase in the Piemont-Ligurian domain. The observed thermal structure is interpreted as the result of the upper and lower unit juxtaposition along shear zones at a temperature of similar to 500 A degrees C during the Middle Eocene. This juxtaposition probably occurred at shallow crustal levels (similar to 15-20 km) within a subduction channel. We finally propose that the Piemont-Ligurian Domain should not be viewed as two distinct "zones", but rather as a stack of several tectonic slices.
This study examines patterns of climate variability by mapping shoreline changes between 1973 and 2008 for two shallow lakes in Malawi, East Africa. Multi-temporal LANDSAT data covering eight different years within the period of investigation were utilised for lake-area mapping, using image classification techniques. The approach was verified using simple comparisons with local rainfall data and satellite altimetry data. Results indicated that the lake areas varied between the different years investigated and that, although the lakes were affected by the same climate-forcing mechanisms, the individual basin characteristics had a critical effect on their responses. The most likely drivers for the climate-related shoreline variations are the El Nino/Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD), causing synchronous changes to the water levels in both lakes, albeit showing differences in climate signal amplification.
This study presents a reconstruction of the Late Holocene climate in Kamchatka based on chironomid remains from a 332 cm long composite sediment core recovered from Dvuyurtochnoe Lake (Two-Yurts Lake, TYL) in central Kamchatka. The oldest recovered sediments date to about 4500 cal years BP. Chironomid head capsules from TYL reflect a rich and diverse fauna. An unknown morphotype of Tanytarsini, Tanytarsus type klein, was found in the lake sediments. Our analysis reveals four chironomid assemblage zones reflecting four different climatic periods in the Late Holocene. Between 4500 and 4000 cal years BP, the chironomid composition indicates a high lake level, well-oxygenated lake water conditions and close to modern temperatures (similar to 13 degrees C). From 4000 to 1000 cal years BP, two consecutive warm intervals were recorded, with the highest reconstructed temperature reaching 16.8 degrees C between 3700 and 2800 cal years BP. Cooling trend, started around 1100 cal years BP led to low temperatures during the last stage of the Holocene. Comparison with other regional studies has shown that termination of cooling at the beginning of late Holocene is relatively synchronous in central Kamchatka, South Kurile, Bering and Japanese Islands and take place around 3700 cal years BP. From ca 3700 cal years BP to the last millennium, a newly strengthened climate continentality accompanied by general warming trend with minor cool excursions led to apparent spatial heterogeneity of climatic patterns in the region. Some timing differences in climatic changes reconstructed from chironomid record of TYL sediments and late Holocene events reconstructed from other sites and other proxies might be linked to differences in local forcing mechanisms or caused by the different degree of dating precision, the different temporal resolution, and the different sensitive responses of climate proxies to the climate variations. Further high-resolution stratigraphic studies in this region are needed to understand the spatially complex pattern of climate change in Holocene in Kamchatka and the surrounding region.
GEOPHILUS ELECTRICUS (nickname GEOPHILUS) is a novel system for mapping the complex electrical bulk resistivity of soils. Rolling electrodes simultaneously measure amplitude and phase data at frequencies ranging from 1 mHz to 1 kHz. The sensor's design and technical specifications allow for measuring these parameters at five depths of up to ca. 1.5 m. Data inversion techniques can be employed to determine resistivity models instead of apparent values and to image soil layers and their geometry with depth. When used in combination with a global positioning system (GPS) and a suitable cross-country vehicle, it is possible to map about 100 ha/day (assuming 1 data point is recorded per second and the line spacing is 18 m). The applicability of the GEOPHILUS system has been demonstrated on several sites, where soils show variations in texture, stratification, and thus electrical characteristics. The data quality has been studied by comparison with 'static' electrodes, by repeated measurements, and by comparison with other mobile conductivity mapping devices (VERIS3100 and EM38). The high quality of the conductivity data produced by the GEOPHILUS system is evident and demonstrated by the overall consistency of the individual maps, and in the clear stratification also confirmed by independent data.
The GEOPHILUS system measures complex values of electrical resistivity in terms of amplitude and phase. Whereas electrical conductivity data (amplitude) are well established in soil science, the interpretation of phase data is a topic of current research. Whether phase data are able to provide additional information depends on the site-specific settings. Here, we present examples, where phase data provide complementary information on man-made structures such as metal pipes and soil compaction.
Temperate forest soils of central Europe are regarded as important pools for soil organic carbon (SOC) and thought to have a high potential for carbon (C) sequestration. However, comprehensive data on total SOC storage, particularly under different forest types, and its drivers is limited. In this study, we analyzed a forest data set of 596 completely sampled soil profiles down to the parent material or to a depth of 1 m within Bavaria in southeast Germany in order to determine representative SOC stocks under different forest types in central Europe and the impact of different environmental parameters. We calculated a total median SOC stock of 9.8 kg m(-2) which is considerably lower compared with many other inventories within central Europe that used modelled instead of measured soil properties. Statistical analyses revealed climate as controlling parameter for the storage of SOC with increasing stocks in cool, humid mountainous regions and a strong decrease in areas with higher temperatures. No significant differences of total SOC storage were found between broadleaf, coniferous and mixed forests. However, coniferous forests stored around 35% of total SOC in the labile organic layer that is prone to human disturbance, forest fires and rising temperatures. In contrast, mixed and broadleaf forests stored the major part of SOC in the mineral soil. Moreover, these two forest types showed unchanged or even slightly increased mineral SOC stocks with higher temperatures, whereas SOC stocks in mineral soils under coniferous forest were distinctly lower. We conclude that mixed and broadleaf forests are more advantageous for C sequestration than coniferous forests. An intensified incorporation of broadleaf species in extent coniferous forests of Bavaria would prevent substantial SOC losses as a result of rising temperatures in the course of climate change.
Considerable effort has been devoted to the development of simulation algorithms for facies modeling, whereas a discussion of how to combine those techniques has not existed. The integration of multiple geologic data into a three-dimensional model, which requires the combination of simulation techniques, is yet a current challenge for reservoir modeling. This article presents a thought process that guides the acquisition and modeling of geologic data at various scales. Our work is based on outcrop data collected from a Jurassic carbonate ramp located in the High Atlas mountain range of Morocco. The study window is 1 km (0.6 mi) wide and 100 m (328.1 ft) thick. We describe and model the spatial and hierarchical arrangement of carbonate bodies spanning from largest to smallest: (1) stacking pattern of high-frequency depositional sequences, (2) facies association, and (3) lithofacies. Five sequence boundaries were modeled using differential global position system mapping and light detection and ranging data. The surface-based model shows a low-angle profile with modest paleotopographic relief at the inner-to-middle ramp transition. Facies associations were populated using truncated Gaussian simulation to preserve ordered trends between the inner, middle, and outer ramps. At the lithofacies scale, field observations and statistical analysis show a mosaiclike distribution that was simulated using a fully stochastic approach with sequential indicator simulation.
This study observes that the use of one single simulation technique is unlikely to correctly model the natural patterns and variability of carbonate rocks. The selection and implementation of different techniques customized for each level of the stratigraphic hierarchy will provide the essential computing flexibility to model carbonate settings. This study demonstrates that a scale-dependent modeling approach should be a common procedure when building subsurface and outcrop models.
Palaeoenvironmental implications of evaporative gaylussite crystals from Lonar Lake, central India
(2013)
We have undertaken petrographic, mineralogical, geochemical and isotopic investigations on carbonate minerals found within a 10-m-long core from Lonar Lake, central India, with the aim of evaluating their potential as palaeoenvironmental proxies. The core encompasses the entire Holocene and is the first well-dated high-resolution record from central India. While calcite and/or aragonite were found throughout the core, the mineral gaylussite was found only in two specific intervals (46303890 and 2040560 cal a BP). Hydrochemical and isotope data from inflowing streams and lake waters indicate that evaporitic processes play a dominant role in the precipitation of carbonates within this lake. Isotopic (18O and 13C) studies on the evaporative gaylussite crystals and residual bulk carbonates (calcite) from the long core show that evaporation is the major control on 18O enrichment in both the minerals. However, in case of 13C additional mechanisms, for example methanogenesis (gaylussite) and phytoplankton productivity (calcium carbonate), play an additional important role in some intervals. We also discuss the relevance of our investigation for palaeoclimate reconstruction and late Holocene monsoon variability.
The persistence of topography within ancient orogens remains one of the outstanding questions in landscape evolution. In the eastern North American Appalachians, this question is manifest in the outstanding problem of whether topographic relief is in a quasi-equilibrium state, decaying slowly over many millennia, or whether relief has increased during the late Cenozoic. Here we present quantitative geomorphic data from the nonglaciated portion of the Susquehanna River drainage basin that provide insight into these end-member models. Analysis of channel profiles draining upland catchments in the northern Valley and Ridge, Appalachian Plateau, Blue Ridge, and Piedmont provinces reveals that a large number of streams have well defined knickpoints clustered at 300-600 m elevation but not systematically associated with transitions from weak to resistant substrate. Cosmogenic Be-10 inventories of modern stream sediment indicate that erosion rates are spatially variable, ranging from similar to 5-30 m/Myr above knickpoints to similar to 50-100 m/Myr below knickpoints. Overall, channel gradients, normalized for drainage area, scale linearly with catchment-averaged erosion rates. Collectively, regionally consistent spatial relationships among erosion rate, channel steepness, and knickpoints reveal an ongoing wave of transient channel adjustment to a change in relative base level. Reconstructions of relict channel profiles above knickpoints suggest that higher rates of incision are associated with similar to 100-150 m of relative base level fall that accompanied epierogenic rock uplift rather than a change to a more erosive climate or drainage reorganization. Channel response timescales imply that the onset of relative base level change predates similar to 3.5 Ma and may have begun as early as similar to 15 Ma. We suggest that adjustment of the channel network was likely driven by changes in mantle dynamics along the eastern seaboard of North America during the Neogene.
Soil moisture dynamics are affected by complex interactions among several factors. Understanding the relative importance of these factors is still an important challenge in the study of water fluxes and solute transport in unsaturated media. In this study, the spatio-temporal variability of surface soil moisture was investigated in a 10 ha flat cropped field located in northern Italy. Soil moisture was measured on a regular 50 x 50 m grid on seven dates during the growing season. For each measurement campaign, the spatial variability of the soil moisture was compared with the spatial variability of the soil texture and crop properties. In particular, to better understand the role of the vegetation, the spatio-temporal variability of two different parameters - leaf area index and crop height - was monitored on eight dates at different crop development stages. Statistical and geostatistical analysis was then applied to explore the interactions between these variables. In agreement with other studies, the results show that the soil moisture variability changes according to the average value within the field, with the standard deviation reaching a maximum value under intermediate mean soil moisture conditions and the coefficient of variation decreasing exponentially with increasing mean soil moisture. The controls of soil moisture variability change according to the average soil moisture within the field. Under wet conditions, the spatial distribution of the soil moisture reflects the variability of the soil texture. Under dry conditions, the spatial distribution of the soil moisture is affected mostly by the spatial variability of the vegetation. The interaction between these two factors is more important under intermediate soil moisture conditions. These results confirm the importance of considering the average soil moisture conditions within a field when investigating the controls affecting the spatial variability of soil moisture. This study highlights the importance of considering the spatio-temporal variability of the vegetation in investigating soil moisture dynamics, especially under intermediate and dry soil moisture conditions. The results of this study have important implications in different hydrological applications, such as for sampling design, ranking stability application, indirect measurements of soil properties and model parameterisation.
The northeastern margin of the Tibetan Plateau is a tectonically active region consisting of a series of faults with bounded intermountain basins and is located in the transition zone between the Tibetan Plateau and the Loess Plateau. Active deformation that may affect the topography in this region can be quantified using geomorphic indices. Therefore, we applied geomorphic indices such as the hypsometric integral and the stream length gradient index to infer neo-tectonics in the northeastern margin of the Tibetan Plateau. Different time-scaled geodetic leveling data and river incision rates were also integrated into the investigation. The results show that the hypsometric integrals are not significantly affected by lithology but spatially correspond to the hanging walls of thrust faults. The hypsometric integrals are also positively correlated with the leveling data. Although the stream length gradient index is influenced by lithology, its most pronounced anomalies of the stream length gradient are associated with the thrust faults. Consequently, the uplift in the northeast margin of the Tibetan Plateau appeared to be concentrated along the hanging walls of the thrust faults.
Changing food consumption patterns and associated greenhouse gas (GHG) emissions have been a matter of scientific debate for decades. The agricultural sector is one of the major GHG emitters and thus holds a large potential for climate change mitigation through optimal management and dietary changes. We assess this potential, project emissions, and investigate dietary patterns and their changes globally on a per country basis between 1961 and 2007. Sixteen representative and spatially differentiated patterns with a per capita calorie intake ranging from 1,870 to >3,400 kcal/day were derived. Detailed analyses show that low calorie diets are decreasing worldwide, while in parallel diet composition is changing as well: a discernable shift towards more balanced diets in developing countries can be observed and steps towards more meat rich diets as a typical characteristics in developed countries. Low calorie diets which are mainly observable in developing countries show a similar emission burden than moderate and high calorie diets. This can be explained by a less efficient calorie production per unit of GHG emissions in developing countries. Very high calorie diets are common in the developed world and exhibit high total per capita emissions of 3.7-6.1 kg CO2eq./day due to high carbon intensity and high intake of animal products. In case of an unbridled demographic growth and changing dietary patterns the projected emissions from agriculture will approach 20 Gt CO2eq./yr by 2050.
This study provides a detailed analysis of the mid-Holocene to present-day precipitation change in the Asian monsoon region. We compare for the first time results of high resolution climate model simulations with a standardised set of mid-Holocene moisture reconstructions. Changes in the simulated summer monsoon characteristics (onset, withdrawal, length and associated rainfall) and the mechanisms causing the Holocene precipitation changes are investigated. According to the model, most parts of the Indian subcontinent received more precipitation (up to 5 mm/day) at mid-Holocene than at present-day. This is related to a stronger Indian summer monsoon accompanied by an intensified vertically integrated moisture flux convergence. The East Asian monsoon region exhibits local inhomogeneities in the simulated annual precipitation signal. The sign of this signal depends on the balance of decreased pre-monsoon and increased monsoon precipitation at mid-Holocene compared to present-day. Hence, rainfall changes in the East Asian monsoon domain are not solely associated with modifications in the summer monsoon circulation but also depend on changes in the mid-latitudinal westerly wind system that dominates the circulation during the pre-monsoon season. The proxy-based climate reconstructions confirm the regional dissimilarities in the annual precipitation signal and agree well with the model results. Our results highlight the importance of including the pre-monsoon season in climate studies of the Asian monsoon system and point out the complex response of this system to the Holocene insolation forcing. The comparison with a coarse climate model simulation reveals that this complex response can only be resolved in high resolution simulations.
On 2012 May 20 and 29, two damaging earthquakes with magnitudes M-w 6.1 and 5.9, respectively, struck the Emilia-Romagna region in the sedimentary Po Plain, Northern Italy, causing 26 fatalities, significant damage to historical buildings and substantial impact to the economy of the region. The earthquake sequence included four more aftershocks with M-w, >= 5.0, all at shallow depths (about 7-9 km), with similar WNW-ESE striking reverse mechanism. The timeline of the sequence suggests significant static stress interaction between the largest events. We perform here a detailed source inversion, first adopting a point source approximation and considering pure double couple and full moment tensor source models. We compare different extended source inversion approaches for the two largest events, and find that the rupture occurred in both cases along a subhorizontal plane, dipping towards SSW Directivity is well detected for the May 20 main shock, indicating that the rupture propagated unilaterally towards SE. Based on the focal mechanism solution, we further estimate the co-seismic static stress change induced by the May 20 event. By using the rate-and-state model and a Poissonian earthquake occurrence, we infer that the second largest event of May 29 was induced with a probability in the range 0.2-0.4. This suggests that the segment of fault was already prone to rupture. Finally, we estimate peak ground accelerations for the two main events as occurred separately or simultaneously. For the scenario involving hypothetical rupture areas of both main events, we estimate M-w = 6.3 and an increase of ground acceleration by 50 per cent. The approach we propose may help to quantify rapidly which regions are invested by a significant increase of the hazard, bearing the potential for large aftershocks or even a second main shock.
Reliable estimations of magnitude of completeness (M-c) are essential for a correct interpretation of seismic catalogues. The spatial distribution of M-c may be strongly variable and difficult to assess in mining environments, owing to the presence of galleries, cavities, fractured regions, porous media and different mineralogical bodies, as well as in consequence of inhomogeneous spatial distribution of the seismicity. We apply a 3-D modification of the probabilistic magnitude of completeness (PMC) method, which relies on the analysis of network detection capabilities. In our approach, the probability to detect an event depends on its magnitude, source receiver Euclidian distance and source receiver direction. The suggested method is proposed for study of the spatial distribution of the magnitude of completeness in a mining environment and here is applied to a 2-months acoustic emission (AE) data set recorded at the Morsleben salt mine, Germany. The dense seismic network and the large data set, which includes more than one million events, enable a detailed testing of the method. This method is proposed specifically for strongly heterogeneous media. Besides, it can also be used for specific network installations, with sensors with a sensitivity, dependent on the direction of the incoming wave (e.g. some piezoelectric sensors). In absence of strong heterogeneities, the standards PMC approach should be used. We show that the PMC estimations in mines strongly depend on the source receiver direction, and cannot be correctly accounted using a standard PMC approach. However, results can be improved, when adopting the proposed 3-D modification of the PMC method. Our analysis of one central horizontal and vertical section yields a magnitude of completeness of about M-c approximate to 1 (AE magnitude) at the centre of the network, which increases up to M-c approximate to 4 at further distances outside the network; the best detection performance is estimated for a NNE-SSE elongated region, which corresponds to the strike direction of the low-attenuating salt body. Our approach provides us with small-scale details about the capability of sensors to detect an earthquake, which can be linked to the presence of heterogeneities in specific directions. Reduced detection performance in presence of strong structural heterogeneities (cavities) is confirmed by synthetic waveform modelling in heterogeneous media.
Supra-glacial deposition and flux of catastrophic rock-slope failure debris, south-central Alaska
(2013)
The ongoing debate over the effects of global environmental change on Earth's cryosphere calls for detailed knowledge about process rates and their variability in cold environments. In this context, appraisals of the coupling between glacier dynamics and para-glacial erosion rates in tectonically active mountains remain rare. We contribute to filling this knowledge gap and present an unprecedented regional-scale inventory of supra-glacial sediment flux and hillslope erosion rates inferred from an analysis of 123 large (> 0 center dot 1km2) catastrophic bedrock landslides that fell onto glaciers in the Chugach Mountains, Alaska, as documented by satellite images obtained between 1972 to 2008. Assuming these supra-glacial landslide deposits to be passive strain markers we infer minimum decadal-scale sediment yields of 190 to 7400tkm-2yr-1 for a given glacier-surface cross-section impacted by episodic rock-slope failure. These rates compare to reported fluvial sediment yields in many mountain rivers, but are an order of magnitude below the extreme sediment yields measured at the snouts of Alaskan glaciers, indicating that the bulk of debris discharged derives from en-glacial, sub-glacial or ice-proximal sources. We estimate an average minimum para-glacial erosion rate by large, episodic rock-slope failures at 0 center dot 5-0 center dot 7mmyr-1 in the Chugach Mountains over a 50-yr period, with earthquakes likely being responsible for up to 73% of this rate. Though ranking amongst the highest decadal landslide erosion rates for this size of study area worldwide, our inferred rates of hillslope erosion in the Chugach Mountains remain an order of magnitude below the pace of extremely rapid glacial sediment export and glacio-isostatic surface uplift previously reported from the region.
This paper presents a simple approach for the temporal disaggregation from daily to 3-hourly observed gridded temperature and precipitation (1 x 1km(2)) on the national scale. The intended use of the disaggregated 3-hourly data is to recalibrate the hydrological model currently used by the Norwegian Water Resources and Energy Directorate (NVE) for daily flood forecasting. By adapting the hydrological model to a 3-hourly temporal scale, the flood forecasting can benefit from available meteorological forecasts with finer temporal resolution and can better represent critical events of short duration and at small spatial scales. By consulting the temporal patterns of a High-Resolution Limited-Area Model (HIRLAM) hindcast series for northern Europe with an hourly temporal and a 0.1 degrees spatial resolution, existing daily 1 x 1km(2) grids for temperature and precipitation covering all of Norway (the seNorge data) were disaggregated into 3-hourly values for the time period September 1957 to December 2010. For the period 2000-05, the disaggregated 3-hourly temperature and precipitation data are validated against observed values from five meteorological stations and against 3-hourly data from the HIRLAM hindcast and daily seNorge data simply split into eight fractions. The results show that the disaggregated data perform best with anomaly correlation coefficients between 0.89 and 0.92 for temperature. With regard to precipitation, the disaggregated data also provide the highest correlations and the lowest errors. In addition, the disaggregated data prove to be best in estimating intervals without precipitation and tend to be most appropriate in estimating extreme precipitation with low occurrence probability (<20%).
The Pleistocene archeological record in East Africa has revealed unusual accumulations of Acheulean handaxes at prehistoric sites. In particular, there has been intensive debate concerning whether the artifact accumulation at the Middle Pleistocene Olorgesailie (Southern Kenya Rift) and Kariandusi (Central Kenya Rift) sites were a result of fluvial reworking or of in situ deposition by hominids. We used a two-step approach to test the hypothesis of fluvial reworking. Firstly, the behavior of handaxes in water currents was investigated in a current flume and the flow threshold required to reorientate the handaxes was determined. The results of these experiments suggested that, in relatively high energy and non-steady flow conditions, handaxes will reorientate themselves perpendicular to the current direction. Secondly, an automated image analysis routine was developed and applied to archeological plans from three Acheulean sites, two at Olorgesailie and one at Kariandusi, in order to determine the orientations of the handaxes. A Rayleigh test was then applied to the orientation data to test for a preferred orientation. The results revealed that the handaxes at the Upper Kariandusi Site and the Olorgesailie Main Site Mid Trench had a preferential orientation, suggesting reworking by a paleocurrent. The handaxes from the Olorgesailie Main Site H/6A, however, appeared to be randomly oriented and in situ deposition by the producers therefore remains a possibility.
In the last decade growing interest has emerged in quantifying the spatial and temporal variations in mountain building. Until recently, insufficient data have been available to attempt such a task at the scale of large orogens such as the Himalaya. The Himalaya accommodates ongoing convergence between India and Eurasia and is a focal point for studying orogen evolution and hypothesized interactions between tectonics and climate. Here we integrate 1126 published bedrock mineral cooling ages with a transient 1D Monte-Carlo thermal-kinematic erosion model to quantify the denudation histories along similar to 2700 km of the Himalaya. The model free parameter is a temporally variable denudation rate from 50 Ma to present. Thermophysical material properties and boundary conditions were tuned to individual study areas. Monte-Carlo simulations were conducted to identify the range of denudation histories that can reproduce the observed cooling ages. Results indicate large temporal and spatial variations in denudation and these are resolvable across different tectonic units of the Himalaya. More specifically, across > 1000 km of the southern Greater Himalaya denudation rates were highest (similar to 1.5-3 mm/yr) between similar to 10 and 2 Ma and lower (0.5-2.6 mm/yr) over the last 2 My. These differences are best determined in the NW-Himalaya. In contrast to this, across the similar to 2500 km length of the northern Greater Himalaya denudation rates vary over length scales of similar to 300-1700 km. Slower denudation (<1 mm/yr) occurred between 10 and 4 Ma followed by a large increase (1.2-2.6 mm/yr) in the last similar to 4 Ma. We find that only the southern Greater Himalayan Sequence clearly supports a continuous co-evolution of tectonics, climate and denudation. Results from the higher elevation northern Greater Himalaya suggest either tectonic driven variations in denudation due to a ramp-flat geometry in the main decollement and/or recent glacially enhanced denudation.
The northward motion of the Pamir indenter with respect to Eurasia has resulted in coeval thrusting, strike-slip faulting, and normal faulting. The eastern Pamir is currently deformed by east-west oriented extension, accompanied by uplift and exhumation of the Kongur Shan (7719m) and Muztagh Ata (7546m) gneiss domes. Both domes are an integral part of the footwall of the Kongur Shan extensional fault system (KES), a 250 km long, north-south oriented graben. Why active normal faulting within the Pamir is primarily localized along the KES and not distributed more widely throughout the orogen has remained unclear. In addition, relatively little is known about how deformation has evolved throughout the Cenozoic, despite refined estimates on present-day crustal deformation rates and microseismicity, which indicate where crustal deformation is presently being accommodated. To better constrain the spatiotemporal evolution of faulting along the KES, we present 39 new apatite fission track, zircon U-Th-Sm/He, and Ar-40/Ar-39 cooling ages from a series of footwall transects along the KES graben shoulder. Combining these data with present-day topographic relief, 1-D thermokinematic and exhumational modeling documents successive stages, rather than synchronous deformation and gneiss dome exhumation. While the exhumation of the Kongur Shan commenced during the late Miocene, extensional processes in the Muztagh Ata massif began earlier and have slowed down since the late Miocene. We present a new model of synorogenic extension suggesting that thermal and density effects associated with a lithospheric tear fault along the eastern margin of the subducting Alai slab localize extensional upper plate deformation along the KES and decouple crustal motion between the central/western Pamir and eastern Pamir/Tarim basin.
Basement-cored ranges formed by reverse faulting within intracontinental mountain belts are often composed of poly-deformed lithologies. Geological data capable of constraining the timing, magnitude, and distribution of the most recent deformational phase are usually missing in such ranges. In this paper, we present new low temperature thermochronological and geological data from a transect through the basement-cored Terskey Range, located in the Kyrgyz Tien Shan. Using these data, we are able to investigate the range's late Cenozoic deformation for the first time. Displacements on reactivated faults are constrained and deformation of thermochronologically derived structural markers is assessed. These structural markers postdate the earlier deformational phases, providing the only record of Cenozoic deformation and of the reactivation of structures within the Terskey Range. Overall, these structural markers have a southern inclination, interpreted to reflect the decreasing inclination of the reverse fault bounding the Terskey Range. Our thermochronological data are also used to investigate spatial and temporal variations in the exhumation of the Terskey Range, identifying a three-stage Cenozoic exhumation history: (1) virtually no exhumation in the Paleogene, (2) increase to slightly higher exhumation rates at similar to 26-20Ma, and (3) significant increase in exhumation starting at similar to 10Ma.
The Sierra de Aconquija and Cumbres Calchaquies in the thick-skinned northern Sierras Pampeanas, NW Argentina present an ideal setting to investigate the tectonically and erosionally controlled exhumation and uplift history of mountain ranges using thermochronological methods. Although these ranges are located along strike of one another, their spatiotemporal evolution varies significantly. Integrating modeled cooling histories constrained by K-Ar ages of muscovite and biotite, apatite fission track data as well as (U-Th)/He measurement of zircon and apatite reveal the structural evolution of these ranges beginning in the late stage of the Paleozoic Famatinian Orogeny. Following localized rift-related exhumation in the central part of the study area and slow erosion elsewhere, growth of the modern topography commenced in the Cenozoic during Andean deformation. The main activity occurred during the late Miocene, with varying magnitudes of rock uplift, surface uplift, and exhumation in the two mountain ranges. The Cumbres Calchaquies is characterized by a total of 5-7km of vertical rock uplift, around 3km of crestal surface uplift, and a maximum exhumation of 2-4km since that time. The Sierra de Aconquija experienced 10-13km of vertical rock uplift, similar to 4-5km of peak surface uplift, and 6-8km of exhumation since around 9Ma. Much of this exhumation occurred along a previously poorly recognized fault. Miocene reactivation of Cretaceous rift structures may explain along-strike variations within these ranges. Dating of sedimentary samples from adjacent basins supports the evolutionary model developed for the mountain ranges.
North of Naran in the Kaghan Valley (NE Pakistan), the metabasites of the melange units lying below the Kohistan Arc, contain glaucophane. Typically they reflect blueschist-metamorphic conditions (0.7 GPa, 400 degrees C). Associated graphite-rich metapelites with quartz veins underwent upper greenschist to amphibolite conditions. Near Naran we observed in quartz grains of type 3 veins first minute relics of Fe-Mg carpholite indicating earlier blueschist metamorphic conditions. P-T estimates indicate 1.2-1.6 GPa at 380-410 degrees C, pressure and temperature values typical for blueschist metamorphic conditions. Changes in mineral assemblages and abundant sudoite component in associated chlorite point to a pressure drop after peak I conditions. We assign the observed changes to peak I conditions occurring during a Cretaceous subduction event. Temperatures estimated with Raman graphite-thermometry clearly indicate a significant subsequent rise of post-peak I temperatures up to 500 degrees C. This is compatible with the amphibolite peak II assigned to the Tertiary continental collision that leads to subduction of the Indian Plate and ultra-high-pressure metamorphism. During subduction the blueschist metamorphic metapelites underwent dehydration, which caused alteration in the overlying lithospheric mantle. In a hydrated lithospheric mantle density is significantly reduced which enhanced subduction of continental crust in the Higher Himalaya. This P-T evolution is typical for a collision orogen with a high plateau but remarkably contrasting findings from Eastern Anatolia, where plateau building is in "statu nascendi" (e.g., Oberhänsli et al., 2010).
Timber harvesting by clear cutting is known to impose environmental impacts, including severe disturbance of the soil hydraulic properties which intensify the frequency and magnitude of surface runoff and soil erosion. However, it remains unanswered if harvest areas act as sources or sinks for runoff and soil erosion and whether such behavior operates in a steady state or evolves through time. For this purpose, 92 small-scale rainfall simulations of different intensities were carried out under pine plantation conditions and on two clear-cut harvest areas of different age. Nonparametrical Random Forest statistical models were set up to quantify the impact of environmental variables on the hydrological and erosion response. Regardless of the applied rainfall intensity, runoff always initiated first and yielded most under plantation cover. Counter to expectations, infiltration rates increased after logging activities. Once a threshold rainfall intensity of 20mm/h was exceeded, the younger harvest area started to act as a source for both runoff and erosion after connectivity was established, whereas it remained a sink under lower applied rainfall intensities. The results suggest that the impact of microtopography on surface runoff connectivity and water-repellent properties of the topsoil act as first-order controls for the hydrological and erosion processes in such environments. Fast rainfall-runoff response, sediment-discharge-hystereses, and enhanced postlogging groundwater recharge at catchment scale support our interpretation. At the end, we show the need to account for nonstationary hydrological and erosional behavior of harvest areas, a fact previously unappreciated in predictive models.
The expansion and intensification of soya bean agriculture in southeastern Amazonia can alter watershed hydrology and biogeochemistry by changing the land cover, water balance and nutrient inputs. Several new insights on the responses of watershed hydrology and biogeochemistry to deforestation in Mato Grosso have emerged from recent intensive field campaigns in this region. Because of reduced evapotranspiration, total water export increases threefold to fourfold in soya bean watersheds compared with forest. However, the deep and highly permeable soils on the broad plateaus on which much of the soya bean cultivation has expanded buffer small soya bean watersheds against increased stormflows. Concentrations of nitrate and phosphate do not differ between forest or soya bean watersheds because fixation of phosphorus fertilizer by iron and aluminium oxides and anion exchange of nitrate in deep soils restrict nutrient movement. Despite resistance to biogeochemical change, streams in soya bean watersheds have higher temperatures caused by impoundments and reduction of bordering riparian forest. In larger rivers, increased water flow, current velocities and sediment flux following deforestation can reshape stream morphology, suggesting that cumulative impacts of deforestation in small watersheds will occur at larger scales.
Badlands have long been considered as model landscapes due to their perceived close relationship between form and process. The often intense features of erosion have also attracted many geomorphologists because the associated high rates of erosion appeared to offer the opportunity for studying surface processes and the resulting forms. Recently, the perceived simplicity of badlands has been questioned because the expected relationships between driving forces for erosion and the resulting sediment yield could not be observed. Further, a high variability in erosion and sediment yield has been observed across scales. Finally, denudation based on currently observed erosion rates would have lead to the destruction of most badlands a long time ago. While the perceived simplicity of badlands has sparked a disproportional (compared to the land surface they cover) amount of research, our increasing amount of information has not necessarily increased our understanding of badlands in equal terms. Overall, badlands appear to be more complex than initially assumed. In this paper, we review 40 years of research in the Zin Valley Badlands in Israel to reconcile some of the conflicting results observed there and develop a perspective on the function of badlands as model landscapes. While the data collected in the Zin Valley clearly confirm that spatial and temporal patterns of geomorphic processes and their interaction with topography and surface properties have to be understood, we still conclude that the process of realizing complexity in the "simple" badlands has a model function both for our understanding as well as perspective on all landscape systems.
The M (w) 6.2 Baladeh earthquake occurred on 28 May 2004 in the Alborz Mountains, northern Iran. This earthquake was the first strong shock in this intracontinental orogen for which digital regional broadband data are available. The Baladeh event provides a rare opportunity to study fault geometry and ongoing deformation processes using modern seismological methods. A joint inversion for hypocentres and a velocity model plus a surface-wave group dispersion curve analysis were used to obtain an adapted velocity model, customised for mid- and long-period waveform modelling. Based on the new velocity model, regional waveform data of the mainshock and larger aftershocks (M (w) a parts per thousand yen3.3) were inverted for moment tensors. For the Baladeh mainshock, this included inversion for kinematic parameters. All analysed earthquakes show dominant thrust mechanisms at depths between 14 and 26 km, with NW-SE striking fault planes. The mainshock ruptured a 28A degrees south-dipping area of 24 x 21 km along a north-easterly direction. The rupture plane of the mainshock does not coincide with the aftershock distribution, neither in map view nor with respect to depth. The considered aftershocks form two main clusters. The eastern cluster is associated with the mainshock. The western cluster does not appear to be connected with the rupture plane of the mainshock but, instead, indicates a second activated fault plane dipping at 85A degrees towards the north.
Extremely rare veinlets and reaction textures composed of symplectites of olivine (similar to Fo(43-55)) + plagioclase +/- spinel +/- ilmenite, associated with more common pyroxene + plagioclase and amphibole + plagioclase varieties, are preserved within eclogites and garnet pyroxenites in the Moldanubian Zone of the Bohemian Massif. Thermodynamic modelling integrated with conventional geothermometry conducted on an eclogite reveals that the symplectite-forming stage occurred at high T (similar to 850 degrees C) and low P (< 6 and > 2 center dot 5 kbar). The development of the different symplectite types reflects reactions that took place in micro-scale domains. The breakdown of high-P garnet controlled the formation of olivine-bearing and amphibole + plagioclase symplectites, whereas breakdown of high-P omphacite led to formation of pyroxene + plagioclase symplectites. In addition, post-eclogite facies but pre-symplectite stage porphyroblastic amphibole and phlogopite were also replaced by olivine-bearing symplectites. Material transfer calculations and thermodynamic modelling indicate that the formation of different symplectite types was linked despite their different bulk compositions. For example, the olivine-bearing symplectites gained Fe +/- Mg, whereas adjacent amphibole + plagioclase and pyroxene + plagioclase symplectites show losses in Fe and Mg; Al, Si and Ca were also variably exchanged. The olivine-bearing symplectites were particularly sensitive to Na despite the small concentration of this element. In eclogites where Na was readily available, the plagioclase composition in the olivine-bearing symplectites shifted from pure anorthite to bytownite, with the less calcic feldspar partitioning Si and inhibiting the formation of orthopyroxene. This regional high-T, low-P granulite-facies symplectite overprint may have been caused by advective heat loss from rapidly exhumed high-T, high-P granulitic bodies (Gfohl Unit) that were emplaced into and over the middle crust (Monotonous and Varied Series) during Carboniferous continent-continent collision.
In order to speed up Precise Point Positioning (PPP)'s convergence, a combined PPP method with GPS and GLONASS which is based on using raw observations is proposed, and the positioning results and convergence time have been compared with that of single system. The ionospheric delays and receiver's Differential Code Bias (DCB) corrections are estimated as unknown parameters in this method. The numerical results show that the combined PPP has not caused significant impacts on the final solutions, but it greatly improved Position Dilution of Precision (PDOP) and convergence speed and enhanced the reliability of the solution. Meanwhile, the convergence speed is greatly influenced by the receiver's DCB, positioning results in horizontal which are better than 10 cm can be realized within 10 min. In addition, the ionosphere and DCB products can be provided with high precision.
Automated seismic event location by travel-time stacking an application to mining induced seismicity
(2013)
The Sabzevar ophiolites mark the Neotethys suture in east-north-central Iran. The Sabzevar metamorphic rocks, as part of the Cretaceous Sabzevar ophiolitic complex, consist of blueschist, amphibolite and greenschist. The Sabzevar blueschists contain sodic amphibole, epidote, phengite, calcite +/- A omphacite +/- A quartz. The epidote amphibolite is composed of sodic-calcic amphibole, epidote, albite, phengite, quartz +/- A omphacite, ilmenite and titanite. The greenschist contains chlorite, plagioclase and pyrite, as main minerals. Thermobarometry of a blueschist yields a pressure of 13-15.5 kbar at temperatures of 420-500 A degrees C. Peak metamorphic temperature/depth ratios were low (similar to 12 A degrees C/km), consistent with metamorphism in a subduction zone. The presence of epidote in the blueschist shows that the rocks were metamorphosed entirely within the epidote stability field. Amphibole schist samples experienced pressures of 5-7 kbar and temperatures between 450 and 550 A degrees C. The presence of chlorite, actinolite, biotite and titanite indicate greenschist facies metamorphism. Chlorite, albite and biotite replacing garnet or glaucophane suggests temperatures of > 300 A degrees C for greenschist facies. The formation of high-pressure metamorphic rocks is related to north-east-dipping subduction of the Neotethys oceanic crust and subsequent closure during lower Eocene between the Central Iranian Micro-continent and Eurasia (North Iran).
Clear S-to-P converted waves from the crust-mantle boundary (Moho) and lithosphere-asthenosphere boundary (LAB) have been observed on the eastern part of the Dead Sea Basin (DSB), and are used for the determination of the depth of the Moho and the LAB. A temporary network consisting of 18 seismic broad-band stations was operated in the DSB region as part of the DEad Sea Integrated REsearch project for 1.5 years beginning in September 2006. The obtained Moho depth (similar to 35 km) from S-to-P receiver functions agrees well with the results from P-to-S receiver functions and other geophysical data. The thickness of the lithosphere on the eastern part of the DSB is about 75 km. The results obtained here support and confirm previous studies, based on xenolith data, geodynamic modeling, heat flow observations, and S-to-P receiver functions. Therefore, the lithosphere on the eastern part of the DSB and along Wadi Araba has been thinned in the Late Cenozoic, following rifting and spreading of the Red Sea. The thinning of the lithosphere occurred without a concomitant change in the crustal thickness and thus an upwelling of the asthenosphere in the study area is invoked as the cause of the lithosphere thinning.
We integrate the information of multiple tomographic models acquired from the earth's surface by modifying a statistical approach recently developed for the integration of cross-borehole tomographic models. In doing so, we introduce spectral cluster analysis as the new core of the model integration procedure to capture the spatial heterogeneity present in all considered tomographic models and describe this heterogeneity in a fuzzy sense. Because spectral cluster algorithms analyze model structure locally, they are considered relatively robust with regard to systematically and spatially varying imaging capabilities typical for geophysical tomographic surveys conducted on the earth's surface. Using a synthetic aquifer example, a fuzzy spectral cluster algorithm can be used to integrate the information provided by 2D tomographic refraction seismic and DC resistivity surveys. The integrated information in the fuzzy membership domain is then used to derive an integrated zonal geophysical model outlining the major structural units present in both input models. We also explain how the fuzzy membership information can be used to identify optimal locations for sparse logging of additional target parameters, i.e., porosity information in our synthetic example. We demonstrate how this sparse porosity information can be extrapolated based on all tomographic input models. The resultant 2D porosity model matches the original porosity distribution reasonably well within the spatial resolution limits of the underlying tomographic models. Consecutively, we apply this approach to a field data base acquired over a former river channel. Sparse information about natural gamma radiation is available and extrapolated on the basis of the fuzzy membership information obtained by spectral cluster analysis of 2D P-wave velocity and electrical resistivity models. This field data shows that the presented parameter extrapolation procedure is robust, even if the locations of target parameter acquisition have not been optimized with regard to the fuzzy membership information.
A total of 271 pollen records were selected from a large collection of both raw and digitized pollen spectra from eastern continental Asia (70 degrees-135 degrees E and 18 degrees-55 degrees N). Following pollen percentage recalculations, taxonomic homogenization, and age-depth model revision, the pollen spectra were interpolated at a 500-year resolution and a taxonomically harmonized and temporally standardized fossil pollen dataset established with 226 pollen taxa, covering the last 22 cal lea. Of the 271 pollen records, 85% were published since 1990, with reliable chronologies and high temporal resolutions; of these, 50% have raw data with complete pollen assemblages, ensuring the quality of this dataset The pollen records available for each 500-year time slice are well distributed over all main vegetation types and climatic zones of the study area, making their pollen spectra suitable for paleovegetation and paleoclimate research. Such a dataset can be used as an example for the development of similar datasets for other regions of the world.
Deserts are a major source of loess and may undergo substantial wind-erosion as evidenced by yardang fields, deflation pans, and wind-scoured bedrock landscapes. However, there are few quantitative estimates of bedrock removal by wind abrasion and deflation. Here, we report wind-erosion rates in the western Qaidam Basin in central China based on measurements of cosmogenic Be-10 in exhumed Miocene sedimentary bedrock. Sedimentary bedrock erosion rates range from 0.05 to 0.4 mm/yr, although the majority of measurements cluster at 0.125 +/- 0.05 mm/yr. These results, combined with previous work, indicate that strong winds, hyper-aridity, exposure of friable Neogene strata, and ongoing rock deformation and uplift in the western Qaidam Basin have created an environment where wind, instead of water, is the dominant agent of erosion and sediment transport. Its geographic location (upwind) combined with volumetric estimates suggest that the Qaidam Basin is a major source (up to 50%) of dust to the Chinese Loess Plateau to the east. The cosmogenically derived wind erosion rates are within the range of erosion rates determined from glacial and fluvial dominated landscapes worldwide, exemplifying the effectiveness of wind to erode and transport significant quantities of bedrock.
As part of the International Continental Drilling Program deep lake drilling project Paleo Van, we investigated sulfate reduction (SR) in deep sediment cores of the saline, alkaline (salinity 21.4%0, alkalinity 155 m mEq-1, pH 9.81) Lake Van, Turkey. The cores were retrieved in the Northern Basin (NB) and at Ahlat Ridge (AR) and reached a maximum depth of 220 m. Additionally, 65-75 cm long gravity cores were taken at both sites. SR rates (SRR) were low (<22 nmol cm-3 day-1) compared to lakes with higher salinity and alkalinity, indicating that salinity and alkalinity are not limiting SR in Lake Van. Both sites differ significantly in rates and depth distribution of SR. In NB, SRR are up to 10 times higher than at AR. SR could be detected down to 19 mblf (meters below lake floor) at NB and down to 13 mblf at AR. Although SRR were lower at AR than at NB, organic matter (OM) concentrations were higher. In contrast, dissolved OM in the pore water at AR contained more macromolecular OM and less low molecular weight OM.VVe thus suggest, that OM content alone cannot be used to infer microbial activity at Lake Van but that quality of OM has an important impact as well. These differences suggest that biogeochemical processes in lacustrine sediments are reacting very sensitively to small variations in geological, physical, or chemical parameters over relatively short distances.
Convexities in the longitudinal profiles of actively incising rivers are typically considered to represent the morphologic signal of a transient response to external perturbations in tectonic or climatic forcing. Distinguishing such knickzones from those that may be anchored to the channel network by spatial variations in rock uplift, however, can be challenging. Here, we combine stream profile analysis, Be-10 watershed-averaged erosion rates, and numerical modeling of stream profile evolution to evaluate whether knickzones in the Abukuma massif of northeast Japan represent a temporal or spatial change in rock uplift rate in relation to forearc shortening. Knickzones in channels that drain the eastern flank of the Abukuma massif are characterized by breaks in slope-area scaling and separate low-gradient, alluvial upper-channel segments from high-gradient, deeply-incised lower channel segments. Average erosion rates inferred from Be-10 concentrations in modern sediment below knickzones exceed erosion rates above knickzones by 20-50%. Although profile convexities could be interpreted as a transient response to an increase in rock uplift rate associated with slip on the range-bounding fault, geologic constraints on the initiation of fault slip and the magnitude of displacement cannot be reconciled with a recent, spatially uniform increase in slip rate. Rather, we find that knickzone position, stream profile gradients, and basin averaged erosion rates are best explained by a relatively abrupt spatial increase in uplift rate localized above a flat-ramp transition in the fault system. These analyses highlight the importance of considering spatially non-uniform uplift in the interpretation of stream profile evolution and demonstrate that the adjustment of river profiles to fault displacement can provide constraints on fault geometry in actively eroding landscapes. (C) 2013 Elsevier B.V. All rights reserved.
The rhizosphere is a dynamic system strongly influenced by root activity. Roots modify the pH of their surrounding soil causing the soil pH to vary as a function of distance from root surface, location along root axes, and root maturity. Non-invasive imaging techniques provide the possibility to capture pH patterns around the roots as they develop.
We developed a novel fluorescence imaging set up and applied to the root system of two lupin (Lupinus albus L., Lupinus angustifolius L.) and one soft-rush (Juncus effusus L.) species. We grew plants in glass containers filled with soil and equipped with fluorescence sensor foils on the container side walls. We gained highly-resolved data on the spatial distribution of H+ around the roots by taking time-lapse images of the samples over the course of several days.
We showed how the soil pH in the vicinity of roots developed over time to different values from that of the original bulk soil. The soil pH in the immediate vicinity of the root surface varied greatly along the root length, with the most acidic point being at 0.56-3.36 mm behind the root tip. Indications were also found for temporal soil pH changes due to root maturity.
In conclusion, this study shows that this novel optical fluorescence imaging set up is a powerful tool for studying pH developments around roots in situ.
A high resolution multi proxy (pollen, grain size, total organic carbon) record from a small mountain lake (Lake Khuisiin; 46.6 degrees N, 101.8 degrees E; 2270 m a.s.l.) in the south eastern Khangai Mountains of central Mongolia has been used to explore changes in vegetation and climate over the last 1200 years. The pollen data indicates that the vegetation changed from dry steppe dominated by Poaceae and Artemisia (ca AD 760-950), to Larix forest steppe (ca AD 950-1170), Larix Betula forest steppe (ca AD 1170-1380), meadow dominated by Cyperaceae and Poaceae (ca AD 1380-1830), and Larix Betula forest steppe (after similar to AD 1830). The cold-wet period between AD 1380 and 1830 may relate to the Little Ice Age. Environmental changes were generally subtle and climate change seems to have been the major driver of variations in vegetation until at least the early part of the 20th century, suggesting that either the level of human activity was generally low, or the relationship between human activity and vegetation did not alter substantially between AD 760 and 1830. A review of centennial scale moisture records from China and Mongolia revealed that most areas experienced major changes at ca AD 1500 and AD 1900. However, the moisture availability since AD 1500 varied between sites, with no clear regional pattern or relationship to present day conditions. Both the reconstructions and the moisture levels simulation on a millennium scale performed in the MPI Earth System Model indicate that the monsoon-westerlies transition area shows a greater climate variability than those areas influenced by the westerlies, or by the summer monsoon only.
Shear folding, which is also referred to as slip folding, involves shear along planes that are oriented approximately parallel to the axial plane of the fold structure. These planes, which are typically axial-planar cleavage planes, facilitate high-angle reverse slip leading to fold limb rotation and amplification. This study builds on recent advances in our understanding of the role of weak fault zones in facilitating slip on misoriented faults; i.e., faults at a high angle to the maximum principal tectonic stress (sigma(1)). Analysis of folded marine sedimentary rocks from the Variscan of southern Ireland provides unambiguous microstructural evidence for reverse shear on chemically weakened cleavage domains. Significant silica loss in these cleavage domains, and as a consequence marked mechanical weakening, is seen as the primary cause for the reverse slip associated with the shear folding of these sedimentary rocks.