Refine
Year of publication
- 2014 (225) (remove)
Document Type
- Article (195)
- Doctoral Thesis (17)
- Preprint (8)
- Review (4)
- Postprint (1)
Language
- English (225) (remove)
Keywords
- Earthquake source observations (7)
- Holocene (6)
- Tibetan Plateau (4)
- Erosion (3)
- Himalaya (3)
- Seismicity and tectonics (3)
- Aleatory variability (2)
- Andes (2)
- Baseline shift (2)
- Body waves (2)
Institute
- Institut für Geowissenschaften (225) (remove)
We report the results of our investigations on the catchment area, surface sediments, and hydrology of the monsoonal Lonar Lake, central India. Our results indicate that the lake is currently stratified with an anoxic bottom layer, and there is a spatial heterogeneity in the sensitivity of sediment parameters to different environmental processes. In the shallow (0-5 m) near shore oxic-suboxic environments the lithogenic and terrestrial organic content is high and spatially variable, and the organics show degradation in the oxic part. Due to aerial exposure resulting from lake level changes of at least 3m, the evaporitic carbonates are not completely preserved. In the deep water (>5 m) anoxic environment the lithogenics are uniformly distributed and the delta C-13 is an indicator not only for aquatic vs. terrestrial plants but also of lake pH and salinity. The isotopic composition of the evaporites is dependent not only on the isotopic composition of source water (monsoon rainfall and stream inflow) and evaporation, but is also influenced by proximity to the isotopically depleted stream inflow. We conclude that in the deep water environment lithogenic content, and isotopic composition of organic matter can be used for palaeoenvironmental reconstruction.
In wetlands or riparian areas, water withdrawal by plants with access to groundwater or the capillary fringe often causes diurnal groundwater fluctuations. Various approaches use the characteristics of these fluctuations for estimation of daily groundwater evapotranspiration rates. The objective of this paper was to review the available methods, compare them with measured evapotranspiration and assess their recharge assumptions. For this purpose, we employed data of 85 rain-free days of a weighable groundwater lysimeter situated at a grassland site in the Spreewald wetland in north-east Germany. Measurements of hourly recharge and daily evapotranspiration rates were used to assess the different approaches. Our results showed that a maximum of 50% of the day to day variance of the daily evapotranspiration rates could be explained by the approaches based on groundwater fluctuations. Simple and more complex methods performed similarly. For some of the approaches, there were indications that erroneous assumptions compensated each other (e.g., when overestimated recharge counteracted underestimated storage change). We found that the usage of longer time spans resulted in improved estimates of the daily recharge rates and that the estimates were further enhanced by including two night averages. When derived from fitting estimates of recharge or evapotranspiration with according measurements the specific yield, needed to convert changes in water level to water volumes, differed considerably among the methods (from 0.022 to 0.064). Thus, the specific yield can be seen as correction factor that compensates for inadequate process descriptions.
Grain-size distributions offer powerful proxies of past environmental conditions that are related to sediment sorting processes. However, they are often of multimodal character because sediments can get mixed during deposition. To facilitate the use of grain size as palaeoenvironmental proxy, this study aims to distinguish the main detrital processes that contribute to lacustrine sedimentation across the Tibetan Plateau using grain-size end-member modelling analysis. Between three and five robust grain-size end-member subpopulations were distinguished at different sites from similarly-likely end-member model runs. Their main modes were grouped and linked to common sediment transport and depositional processes that can be associated with contemporary Tibetan climate (precipitation patterns and lake ice phenology, gridded wind and shear stress data from the High Asia Reanalysis) and local catchment configurations. The coarse sands and clays with grain-size modes > 250 mu m and < 2 mu m were probably transported by fluvial processes. Aeolian sands (similar to 200 mu m) and coarse local dust (similar to 60 mu m), transported by saltation and in near-surface suspension clouds, are probably related to occasional westerly storms in winter and spring. Coarse regional dust with modes similar to 25 mu m may derive from near-by sources that keep in longer term suspension. The continuous background dust is differentiated into two robust end members (modes: 5-10 and 2-5 mu m) that may represent different sources, wind directions and/or sediment trapping dynamics from long-range, upper-level westerly and episodic northerly wind transport. According to this study grain-size end members of only fluvial origin contribute small amounts to mean Tibetan lake sedimentation (19 +/- 5%), whereas local to regional aeolian transport and background dust deposition dominate the clastic sedimentation in Tibetan lakes (contributions: 42 +/- 14% and 51 +/- 11%). However, fluvial and alluvial reworking of aeolian material from nearby slopes during summer seems to limit end-member interpretation and should be cross-checked with other proxy information. If not considered as a stand-alone proxy, a high transferability to other regions and sediment archives allows helpful reconstructions of past sedimentation history.
Pollen records from large lakes have been used for quantitative palaeoclimate reconstruction, but the influences that lake size (as a result of species-specific variations in pollen dispersal patterns that smaller pollen grains are more easily transported to lake centre) and taphonomy have on these climatic signals have not previously been systematically investigated. We introduce the concept of pollen source area to pollen-based climate calibration using the north-eastern Tibetan Plateau as our study area. We present a pollen data set collected from large lakes in the arid to semi-arid region of central Asia. The influences that lake size and the inferred pollen source areas have on pollen compositions have been investigated through comparisons with pollen assemblages in neighbouring lakes of various sizes. Modern pollen samples collected from different parts of Lake Donggi Cona (in the north-eastern part of the Tibetan Plateau) reveal variations in pollen assemblages within this large lake, which are interpreted in terms of the species-specific dispersal and depositional patterns for different types of pollen, and in terms of fluvial input components. We have estimated the pollen source area for each lake individually and used this information to infer modern climate data with which to then develop a modern calibration data set, using both the multivariate regression tree (MRT) and weighted-averaging partial least squares (WA-PLS) approaches. Fossil pollen data from Lake Donggi Cona have been used to reconstruct the climate history of the north-eastern part of the Tibetan Plateau since the Last Glacial Maximum (LGM). The meanannual precipitation was quantitatively reconstructed using WA-PLS: extremely dry conditions are found to have dominated the LGM, with annual precipitation of around 100 mm, which is only 32% of present-day precipitation. A gradually increasing trend in moisture conditions during the Late Glacial is terminated by an abrupt reversion to a dry phase that lasts for about 1000 yr and coincides with "Heinrich event 1" in the North Atlantic region. Subsequent periods corresponding to the Bolling/Allerod interstadial, with annual precipitation (P-ann) of about 350 mm, and the Younger Dryas event (about 270 mm P-ann) are followed by moist conditions in the early Holocene, with annual precipitation of up to 400 mm. A drier trend after 9 cal. ka BP is followed by a second wet phase in the middle Holocene, lasting until 4.5 cal. ka BP. Relatively steady conditions with only slight fluctuations then dominate the late Holocene, resulting in the present climatic conditions. The climate changes since the LGM have been primarily driven by deglaciation and fluctuations in the intensity of the Asian summer monsoon that resulted from changes in the Northern Hemisphere summer solar insolation, as well as from changes in the North Atlantic climate through variations in the circulation patterns and intensity of the westerlies.
This study aims to compare impacts of climate change on streamflow in four large representative African river basins: the Niger, the Upper Blue Nile, the Oubangui and the Limpopo. We set up the eco-hydrological model SWIM (Soil and Water Integrated Model) for all four basins individually. The validation of the models for four basins shows results from adequate to very good, depending on the quality and availability of input and calibration data.
For the climate impact assessment, we drive the model with outputs of five bias corrected Earth system models of Coupled Model Intercomparison Project Phase 5 (CMIP5) for the representative concentration pathways (RCPs) 2.6 and 8.5. This climate input is put into the context of climate trends of the whole African continent and compared to a CMIP5 ensemble of 19 models in order to test their representativeness. Subsequently, we compare the trends in mean discharges, seasonality and hydrological extremes in the 21st century. The uncertainty of results for all basins is high. Still, climate change impact is clearly visible for mean discharges but also for extremes in high and low flows. The uncertainty of the projections is the lowest in the Upper Blue Nile, where an increase in streamflow is most likely. In the Niger and the Limpopo basins, the magnitude of trends in both directions is high and has a wide range of uncertainty. In the Oubangui, impacts are the least significant. Our results confirm partly the findings of previous continental impact analyses for Africa. However, contradictory to these studies we find a tendency for increased streamflows in three of the four basins (not for the Oubangui). Guided by these results, we argue for attention to the possible risks of increasing high flows in the face of the dominant water scarcity in Africa. In conclusion, the study shows that impact intercomparisons have added value to the adaptation discussion and may be used for setting up adaptation plans in the context of a holistic approach.
Many mountain belts sustain prolonged snow cover for parts of the year, although enquiries into rates of erosion in these landscapes have focused almost exclusively on the snow-free periods. This raises the question of whether annual snow cover contributes significantly to modulating rates of erosion in high-relief terrain. In this context, the sudden release of snow avalanches is a frequent and potentially relevant process, judging from the physical damage to subalpine forest ecosystems, and the amount of debris contained in avalanche deposits. To quantitatively constrain this visual impression and to expand the sparse literature, we sampled sediment concentrations of n = 28 river-spanning snow-avalanche deposits (snow bridges) in the area around Davos, eastern Swiss Alps, and inferred an orders-of-magnitude variability in specific fine sediment and organic carbon yields (1.8 to 830 t km(-2) yr(-1), and 0.04 to 131 tC km(-2) yr(-1), respectively). A Monte Carlo simulation demonstrates that, with a minimum of free parameters, such variability is inherent to the geometric scaling used for computing specific yields. Moreover, the widely applied method of linearly extrapolating plot scale sample data may be prone to substantial under- or overestimates. A comparison of our inferred yields with previously published work demonstrates the relevance of wet snow avalanches as prominent agents of soil erosion and transporters of biogeochemical constituents to mountain rivers. Given that a number of snow bridges persisted below the insulating debris cover well into the summer months, snow-avalanche deposits also contribute to regulating in-channel sediment and organic debris storage on seasonal timescales. Finally, our results underline the potential shortcomings of neglecting erosional processes in the winter and spring months in mountainous terrain subjected to prominent snow cover.
Making confident statements about the evolution of an ice-sheet shelf system with a numerical model requires the capability to reproduce the migration of the grounding line. Here we show that the shallow-ice approximation/shallow-shelf approximation hybrid-type Parallel Ice Sheet Model (PISM), with its recent improvements, is capable of modeling the grounding line motion in a perturbed ice-sheet shelf system. The model is set up according to the three-dimensional Marine Ice-Sheet Model Intercomparison Project (MISMIP3d), and simulations are carried out across a broad range of spatial resolutions. Using (1) a linear interpolation of the grounding line with locally interpolated basal friction and (2) an improved driving-stress computation across the grounding line, the reversibility of the grounding line (i.e. its retreat after an advance forced by a local perturbation of basal resistance) is captured by the model even at medium and low resolutions (Delta x > 10 km). The transient model response is qualitatively similar to that of higher-order models but reveals a higher initial sensitivity to perturbations on very short timescales. Our findings support the application of PISM to the Antarctic ice sheet from regional up to continental scales and on relatively low spatial resolutions.
Differential exhumation in the Puna Plateau and Eastern Cordillera of NW Argentina is controlled by inherited paleostructures and resulting paleotopography related to the Cretaceous Salta Rift paleomargins. The Ceno zoic deformation front related to the development of the Andean retro-arc orogenic system is generally associated with >4 km of exhumation, which is recorded by Cenozoic apatite fi ssion-track (AFT) and (U-Th-[Sm])/He ages (He ages) in the Eastern Cordillera of NW Argentina. New AFT ages from the top of the Nevado de Cachi document Oligocene (ca. 28 Ma) cooling, which, combined with existing data, indicates exhumation of this range between ca. 28 Ma and ca. 14 Ma. However, some of the highest ranges in the Eastern Cordillera preserve Cretaceous ages indicative of limited Cenozoic exhumation. Samples collected from an similar to 3-km-elevation transect along the northern part of the Sierra de Quilmes paleorift fl ank (Laguna Brava) show AFT ages between ca. 80 and ca. 50 Ma and He ages between ca. 45 and ca. 10 Ma. Another set of samples from an similar to 1-km-elevation transect farther to the southwest (La Quebrada) shows Cretaceous AFT ages between ca. 116 Ma and ca. 76 Ma, and mainly Cretaceous He ages, in agreement with AFT data. Analysis of existing AFT and He ages from the area once occupied by the Salta Rift reveals a pattern characterized by Cretaceous ages along paleorift highs and Cenozoic ages within paleorift hanging-wall basins and later foreland basin depocenters. This pattern is interrupted by the Sierras Pampeanas at similar to 28 degrees S, which record mid-Cenozoic ages. Our data are consistent with a complex inherited pattern of pre-Andean paleostructures, likely associated with paleotopography, which was beveled by the Cenozoic regional foreland basin and reactivated during the late Neogene (ca. <10 Ma), strongly controlling the magnitude of Cenozoic uplift and exhumation and thus cooling age distribution. This, combined with variable lithologic erodibility, resulted in an irregular distribution of themochronological ages.
Effects of climate change are particularly strong in high-mountain regions. Most visibly, glaciers are shrinking at a rapid pace, and as a consequence, glacier lakes are forming or growing. At the same time the stability of mountain slopes is reduced by glacier retreat, permafrost thaw and other factors, resulting in an increasing landslide hazard which can potentially impact lakes and therewith trigger far-reaching and devastating outburst floods. To manage risks from existing or future lakes, strategies need to be developed to plan in time for adequate risk reduction measures at a local level. However, methods to assess risks from future lake outbursts are not available and need to be developed to evaluate both future hazard and future damage potential.
Here a method is presented to estimate future risks related to glacier lake outbursts for a local site in southern Switzerland (Naters, Valais). To generate two hazard scenarios, glacier shrinkage and lake formation modelling was applied, combined with simple flood modelling and field work. Furthermore, a land-use model was developed to quantify and allocate land-use changes based on local-to-regional storylines and three scenarios of land-use driving forces. Results are conceptualized in a matrix of three land-use and two hazard scenarios for the year 2045, and show the distribution of risk in the community of Naters, including high and very high risk areas. The study underlines the importance of combined risk management strategies focusing on land-use planning, on vulnerability reduction, as well as on structural measures (where necessary) to effectively reduce future risks related to lake outburst floods.
Flood hazard projections under climate change are typically derived by applying model chains consisting of the following elements: "emission scenario - global climate model - downscaling, possibly including bias correction hydrological model - flood frequency analysis". To date, this approach yields very uncertain results, due to the difficulties of global and regional climate models to represent precipitation. The implementation of such model chains requires major efforts, and their complexity is high.
We propose for the Mekong River an alternative approach which is based on a shortened model chain: "emission scenario - global climate model - non-stationary flood frequency model". The underlying idea is to use a link between the Western Pacific monsoon and local flood characteristics: the variance of the monsoon drives a non-stationary flood frequency model, yielding a direct estimate of flood probabilities. This approach bypasses the uncertain precipitation, since the monsoon variance is derived from large-scale wind fields which are better represented by climate models. The simplicity of the monsoon-flood link allows deriving large ensembles of flood projections under climate change. We conclude that this is a worthwhile, complementary approach to the typical model chains in catchments where a substantial link between climate and floods is found.
Flood estimation and flood management have traditionally been the domain of hydrologists, water resources engineers and statisticians, and disciplinary approaches abound. Dominant views have been shaped; one example is the catchment perspective: floods are formed and influenced by the interaction of local, catchment-specific characteristics, such as meteorology, topography and geology. These traditional views have been beneficial, but they have a narrow framing. In this paper we contrast traditional views with broader perspectives that are emerging from an improved understanding of the climatic context of floods. We come to the following conclusions: (1) extending the traditional system boundaries (local catchment, recent decades, hydrological/hydraulic processes) opens up exciting possibilities for better understanding and improved tools for flood risk assessment and management. (2) Statistical approaches in flood estimation need to be complemented by the search for the causal mechanisms and dominant processes in the atmosphere, catchment and river system that leave their fingerprints on flood characteristics. (3) Natural climate variability leads to time-varying flood characteristics, and this variation may be partially quantifiable and predictable, with the perspective of dynamic, climate-informed flood risk management. (4) Efforts are needed to fully account for factors that contribute to changes in all three risk components (hazard, exposure, vulnerability) and to better understand the interactions between society and floods. (5) Given the global scale and societal importance, we call for the organization of an international multidisciplinary collaboration and data-sharing initiative to further understand the links between climate and flooding and to advance flood research.
The widely reported paradox of methane oversaturation in oxygenated water challenges the prevailing paradigm that microbial methanogenesis only occurs under anoxic conditions. Using a combination of field sampling, incubation experiments, and modeling, we show that the recurring mid-water methane peak in Lake Stechlin, northeast Germany, was not dependent on methane input from the littoral zone or bottom sediment or on the presence of known micro-anoxic zones. The methane peak repeatedly overlapped with oxygen oversaturation in the seasonal thermocline. Incubation experiments and isotope analysis indicated active methane production, which was likely linked to photosynthesis and/or nitrogen fixation within the oxygenated water, whereas lessening of methane oxidation by light allowed accumulation of methane in the oxygen-rich upper layer. Estimated methane efflux from the surface water was up to 5 mmol m(-2) d(-1). Mid-water methane oversaturation was also observed in nine other lakes that collectively showed a strongly negative gradient of methane concentration within 0-20% dissolved oxygen (DO) in the bottom water, and a positive gradient within >= 20% DO in the upper water column. Further investigation into the responsible organisms and biochemical pathways will help improve our understanding of the global methane cycle.
The paper examines the quality of satellite-abased precipitation estimates for the lower Mahanadi River basin (eastern India). The considered data sets known as 3B42 and 3B42-RT (version 7/7A) are routinely produced by the tropical rainfall measuring mission (TRMM) from passive microwave and infrared recordings. While the 3B42-RT data are disseminated in real time, the gauge-aadjusted 3B42 data set is published with a delay of some months. The quality of the two products was assessed in a two-astep procedure. First, the correspondence between the remotely sensed precipitation rates and rain gauge data was evaluated at the subbasin scale. Second, the quality of the rainfall estimates was assessed by analysing their performance in the context of rainfall-arunoff simulation.
At sub-abasin level (4000 to 16 000 km(2)) the satellite-abased areal precipitation estimates were found to be moderately correlated with the gauge-abased counterparts (R-2 of 0.64-0.74 for 3B42 and 0.59-0.72 for 3B42-RT). Significant discrepancies between TRMM data and ground observations were identified at high-aintensity levels. The rainfall depth derived from rain gauge data is often not reflected by the TRMM estimates (hit rate < 0.6 for ground-abased intensities > 80 mm day(-1)). At the same time, the remotely sensed rainfall rates frequently exceed the gauge-abased equivalents (false alarm ratios of 0.2-0.6). In addition, the real-atime product 3B42-RT was found to suffer from a spatially consistent negative bias.
Since the regionalisation of rain gauge data is potentially associated with a number of errors, the above results are subject to uncertainty. Hence, a validation against independent information, such as stream flow, was essential. In this case study, the outcome of rainfall-arunoff simulation experiments was consistent with the above-mentioned findings. The best fit between observed and simulated stream flow was obtained if rain gauge data were used as model input (Nash-Sutcliffe index of 0.76-0.88 at gauges not affected by reservoir operation). This compares to the values of 0.71-0.78 for the gauge-djusted TRMM 3B42 data and 0.65-0.77 for the 3B42-RT real-atime data. Whether the 3B42-RT data are useful in the context of operational runoff prediction in spite of the identified problems remains a question for further research.
New techniques for high pressure falling sphere viscosimetry in DIA-type large volume presses
(2014)
Magmatism forming the Central Anatolian Volcanic Province of Cappadocia, central Turkey, records the last phase of Neotethyan subduction after similar to 11 Ma. Thirteen large calc-alkaline ignimbrite sheets form marker bands within the volcano-sedimentary succession (the Urgup Formation) and provide a robust chronostratigraphy for paleoecologic evaluation of the interleaved paleosols. This paper evaluates the chronologic record in the context of the radiometric, magnetostratigraphic and lithostratigraphic controls. Previous inconsistencies relating primarily to K/Ar evidence were reason for the initiation of an integrated study which includes Ar-40/Ar-39 dating, palaeomagnetic and stratigraphic evidence. The newly determined Ar-40/Ar-39-ages (Lepetit, 2010) are in agreement with Ar/Ar and U/Pb data meanwhile published by Pauquette and Le Pennec (2012) and Aydar et al. (2012). The Ar-40/Ar-39-ages restrict the end of the Urgup Formation to the late Miocene. The paleosol sequence enclosed by the ignimbrites is thus restricted to the late Miocene, the most intense formation of pedogene calcretes correlating with the Messinian Salinity Crisis.
Modern natural hazards research requires dealing with several uncertainties that arise from limited process knowledge, measurement errors, censored and incomplete observations, and the intrinsic randomness of the governing processes. Nevertheless, deterministic analyses are still widely used in quantitative hazard assessments despite the pitfall of misestimating the hazard and any ensuing risks.
In this paper we show that Bayesian networks offer a flexible framework for capturing and expressing a broad range of uncertainties encountered in natural hazard assessments. Although Bayesian networks are well studied in theory, their application to real-world data is far from straightforward, and requires specific tailoring and adaptation of existing algorithms. We offer suggestions as how to tackle frequently arising problems in this context and mainly concentrate on the handling of continuous variables, incomplete data sets, and the interaction of both. By way of three case studies from earthquake, flood, and landslide research, we demonstrate the method of data-driven Bayesian network learning, and showcase the flexibility, applicability, and benefits of this approach.
Our results offer fresh and partly counterintuitive insights into well-studied multivariate problems of earthquake-induced ground motion prediction, accurate flood damage quantification, and spatially explicit landslide prediction at the regional scale. In particular, we highlight how Bayesian networks help to express information flow and independence assumptions between candidate predictors. Such knowledge is pivotal in providing scientists and decision makers with well-informed strategies for selecting adequate predictor variables for quantitative natural hazard assessments.
This paper employs a complex network approach to determine the topology and evolution of the network of extreme precipitation that governs the organization of extreme rainfall before, during, and after the Indian Summer Monsoon (ISM) season. We construct networks of extreme rainfall events during the ISM (June-September), post-monsoon (October-December), and pre-monsoon (March-May) periods from satellite-derived (Tropical Rainfall Measurement Mission, TRMM) and rain-gauge interpolated (Asian Precipitation Highly Resolved Observational Data Integration Towards the Evaluation of Water Resources, APHRODITE) data sets. The structure of the networks is determined by the level of synchronization of extreme rainfall events between different grid cells throughout the Indian subcontinent. Through the analysis of various complex-network metrics, we describe typical repetitive patterns in North Pakistan (NP), the Eastern Ghats (EG), and the Tibetan Plateau (TP). These patterns appear during the pre-monsoon season, evolve during the ISM, and disappear during the post-monsoon season. These are important meteorological features that need further attention and that may be useful in ISM timing and strength prediction.
Water is an essential input to the majority of human activities. Often, access to sufficient water resources is limited by quality and infrastructure aspects, rather than by resource availability alone, and each activity has different requirements regarding the nature of these aspects. This paper develops an integrated approach to assess the adequacy of water resources for the three major water users: the domestic, agricultural and industrial sectors. Additionally, we include environmental water requirements. We first outline the main determinants of water adequacy for each sector. Subsequently, we present an integrated approach using fuzzy logic, which allows assessing sector-specific as well as overall water adequacy. We implement the approach in two case study settings to exemplify the main features of the approach. Using results from two climate models and two forcing RCPs (representative concentration pathways), as well as population projections, we further assess the impacts of climate change in combination with population growth on the adequacy of water resources. The results provide an important step forward in determining the most relevant factors, impeding adequate access to water, which remains an important challenge in many regions of the world. The methodology allows one to directly identify the factors that are most decisive in determining the adequacy of water in each region, pointing towards the most efficient intervention points to improve conditions. Our findings underline the fact that, in addition to water volumes, water quality is a limitation for all sectors and, especially for the environmental sector, high levels of pollution are a threat to water adequacy.
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.
According to Dooge (1986) intermediate-scale catchments are systems of organized complexity, being too organized and yet too small to be characterized on a statistical/conceptual basis, but too large and too heterogeneous to be characterized in a deterministic manner. A key requirement for building structurally adequate models precisely for this intermediate scale is a better understanding of how different forms of spatial organization affect storage and release of water and energy. Here, we propose that a combination of the concept of hydrological response units (HRUs) and thermodynamics offers several helpful and partly novel perspectives for gaining this improved understanding. Our key idea is to define functional similarity based on similarity of the terrestrial controls of gradients and resistance terms controlling the land surface energy balance, rainfall runoff transformation, and groundwater storage and release. This might imply that functional similarity with respect to these specific forms of water release emerges at different scales, namely the small field scale, the hillslope, and the catchment scale. We thus propose three different types of "functional units" - specialized HRUs, so to speak - which behave similarly with respect to one specific form of water release and with a characteristic extent equal to one of those three scale levels. We furthermore discuss an experimental strategy based on exemplary learning and replicate experiments to identify and delineate these functional units, and as a promising strategy for characterizing the interplay and organization of water and energy fluxes across scales. We believe the thermodynamic perspective to be well suited to unmask equifinality as inherent in the equations governing water, momentum, and energy fluxes: this is because several combinations of gradients and resistance terms yield the same mass or energy flux and the terrestrial controls of gradients and resistance terms are largely independent. We propose that structurally adequate models at this scale should consequently disentangle driving gradients and resistance terms, because this optionally allow sequifinality to be partly reduced by including available observations, e. g., on driving gradients. Most importantly, the thermodynamic perspective yields an energy-centered perspective on rainfall-runoff transformation and evapotranspiration, including fundamental limits for energy fluxes associated with these processes. This might additionally reduce equifinality and opens up opportunities for testing thermodynamic optimality principles within independent predictions of rainfall-runoff or land surface energy exchange. This is pivotal to finding out whether or not spatial organization in catchments is in accordance with a fundamental organizing principle.
Sedimentary proxy records constitute a significant portion of the recorded evidence that allows us to investigate paleoclimatic conditions and variability. However, uncertainties in the dating of proxy archives limit our ability to fix the timing of past events and interpret proxy record intercomparisons. While there are various age-modeling approaches to improve the estimation of the age-depth relations of archives, relatively little focus has been placed on the propagation of the age (and radiocarbon calibration) uncertainties into the final proxy record.
We present a generic Bayesian framework to estimate proxy records along with their associated uncertainty, starting with the radiometric age-depth and proxy-depth measurements, and a radiometric calibration curve if required. We provide analytical expressions for the posterior proxy probability distributions at any given calendar age, from which the expected proxy values and their uncertainty can be estimated. We illustrate our method using two synthetic data sets and then use it to construct the proxy records for groundwater inflow and surface erosion from Lonar lake in central India.
Our analysis reveals interrelations between the uncertainty of the proxy record over time and the variance of proxies along the depth of the archive. For the Lonar lake proxies, we show that, rather than the age uncertainties, it is the proxy variance combined with calibration uncertainty that accounts for most of the final uncertainty. We represent the proxy records as probability distributions on a precise, error-free timescale that makes further time series analyses and intercomparisons of proxies relatively simple and clear. Our approach provides a coherent understanding of age uncertainties within sedimentary proxy records that involve radiometric dating. It can be potentially used within existing age modeling structures to bring forth a reliable and consistent framework for proxy record estimation.
Recurrence-plot-based recurrence networks are an approach used to analyze time series using a complex networks theory. In both approaches - recurrence plots and recurrence networks -, a threshold to identify recurrent states is required. The selection of the threshold is important in order to avoid bias of the recurrence network results. In this paper, we propose a novel method to choose a recurrence threshold adaptively. We show a comparison between the constant threshold and adaptive threshold cases to study period-chaos and even period-period transitions in the dynamics of a prototypical model system. This novel method is then used to identify climate transitions from a lake sediment record.
Climate impacts on human livelihoods: where uncertainty matters in projections of water availability
(2014)
Climate change will have adverse impacts on many different sectors of society, with manifold consequences for human livelihoods and well-being. However, a systematic method to quantify human well-being and livelihoods across sectors is so far unavailable, making it difficult to determine the extent of such impacts. Climate impact analyses are often limited to individual sectors (e.g. food or water) and employ sector-specific target measures, while systematic linkages to general livelihood conditions remain unexplored. Further, recent multi-model assessments have shown that uncertainties in projections of climate impacts deriving from climate and impact models, as well as greenhouse gas scenarios, are substantial, posing an additional challenge in linking climate impacts with livelihood conditions. This article first presents a methodology to consistently measure what is referred to here as AHEAD (Adequate Human livelihood conditions for wEll-being And Development). Based on a trans-disciplinary sample of concepts addressing human well-being and livelihoods, the approach measures the adequacy of conditions of 16 elements. We implement the method at global scale, using results from the Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP) to show how changes in water availability affect the fulfilment of AHEAD at national resolution. In addition, AHEAD allows for the uncertainty of climate and impact model projections to be identified and differentiated. We show how the approach can help to put the substantial inter-model spread into the context of country-specific livelihood conditions by differentiating where the uncertainty about water scarcity is relevant with regard to livelihood conditions - and where it is not. The results indicate that livelihood conditions are compromised by water scarcity in 34 countries. However, more often, AHEAD fulfilment is limited through other elements. The analysis shows that the water-specific uncertainty ranges of the model output are outside relevant thresholds for AHEAD for 65 out of 111 countries, and therefore do not contribute to the overall uncertainty about climate change impacts on livelihoods. In 46 of the countries in the analysis, water-specific uncertainty is relevant to AHEAD. The AHEAD method presented here, together with first results, forms an important step towards making scientific results more applicable for policy decisions.
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.
Deforestation is a prominent anthropogenic cause of erosive overland flow and slope instability, boosting rates of soil erosion and concomitant sediment flux. Conventional methods of gauging or estimating post-logging sediment flux often focus on annual timescales but overlook potentially important process response on shorter intervals immediately following timber harvest. We resolve such dynamics with non-parametric quantile regression forests (QRF) based on high-frequency (3 min) discharge measurements and sediment concentration data sampled every 30-60 min in similar-sized (similar to 0.1 km(2)) forested Chilean catchments that were logged during either the rainy or the dry season. The method of QRF builds on the random forest algorithm, and combines quantile regression with repeated random sub-sampling of both cases and predictors. The algorithm belongs to the family of decision-tree classifiers, which allow quantifying relevant predictors in high-dimensional parameter space. We find that, where no logging occurred, similar to 80% of the total sediment load was transported during extremely variable runoff events during only 5% of the monitoring period. In particular, dry-season logging dampened the relative role of these rare, extreme sediment-transport events by increasing load efficiency during more efficient moderate events. We show that QRFs outperform traditional sediment rating curves (SRCs) in terms of accurately simulating short-term dynamics of sediment flux, and conclude that QRF may reliably support forest management recommendations by providing robust simulations of post-logging response of water and sediment fluxes at high temporal resolution.
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.
Permafrost-affected soils are among the most obvious ecosystems in which current microbial controls on organic matter decomposition are changing as a result of global warming. Warmer conditions in polygonal tundra will lead to a deepening of the seasonal active layer, provoking changes in microbial processes and possibly resulting in exacerbated carbon degradation under increasing anoxic conditions. To identify current microbial assemblages in carbon rich, water saturated permafrost environments, four polygonal tundra sites were investigated on Herschel Island and the Yukon Coast, Western Canadian Arctic. Ion Torrent sequencing of bacterial and archaeal 16S rRNA amplicons revealed the presence of all major microbial soil groups and indicated a local, vertical heterogeneity of the polygonal tundra soil community with increasing depth. Microbial diversity was found to be highest in the surface layers, decreasing towards the permafrost table. Quantitative PCR analysis of functional genes involved in carbon and nitrogen-cycling revealed a high functional potential in the surface layers, decreasing with increasing active layer depth. We observed that soil properties driving microbial diversity and functional potential varied in each study site. These results highlight the small-scale heterogeneity of geomorphologically comparable sites, greatly restricting generalizations about the fate of permafrost-affected environments in a warming Arctic.
Comparing thermal wave function methods for multi-configuration time-dependent Hartree simulations
(2014)
We compare two methods for creating stochastic temperature wave functions that can be used for Multi-Configuration Time-Dependent Hartree (MCTDH) simulations. In the first method, the MCTDH coefficients are chosen randomly, while the other method uses a single Hartree product of random single-particle functions (SPFs). We find that using random SPFs dramatically improves convergence for a model system for surface sticking.
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.
In hydrology, the storage-discharge relationship is a fundamental catchment property. Understanding what controls this relationship is at the core of catchment science. To date, there are no direct methods to measure water storage at catchment scales (10(1)-10(3)km(2)). In this study, we use direct measurements of terrestrial water storage dynamics by means of superconducting gravimetry in a small headwater catchment of the Regen River, Germany, to derive empirical storage-discharge relationships in nested catchments of increasing scale. Our results show that the local storage measurements are strongly related to streamflow dynamics at larger scales (> 100km(2); correlation coefficient=0.78-0.81), but at small scale, no such relationship exists (similar to 1km(2); correlation coefficients=-0.11). The geologic setting in the region can explain both the disconnection between local water storage and headwater runoff, and the connectivity between headwater storage and streams draining larger catchment areas. More research is required to understand what controls the form of the observed storage-discharge relationships at the catchment scale. This study demonstrates that high-precision gravimetry can provide new insights into the complex relationship between state and response of hydrological systems.
Advance in geocomputation
(2014)
In situ UV laser spot Ar-40/Ar-39 analyses of distinct phengite types in eclogite-facies rocks from the Sesia-Lanzo Zone (Western Alps, Italy) were combined with SIMS boron isotope analyses as well as boron (B) and lithium (Li) concentration data to link geochronological information with constraints on fluid-rock interaction. In weakly deformed samples, apparent Ar-40/Ar-39 ages of phengite cores span a range of similar to 20 Ma, but inverse isochrons define two distinct main high-pressure (HP) phengite core crystallization periods of 88-82 and 77-74 Ma, respectively. The younger cores have on average lower B contents (similar to 36 mu g/g) than the older ones (similar to 43-48 mu g/g), suggesting that loss of B and resetting of the Ar isotopic system were related. Phengite cores have variable delta B-11 values (-18 parts per thousand to -10 parts per thousand), indicating the lack of km scale B homogenization during HP crystallization.
Overprinted phengite rims in the weakly deformed samples generally yield younger apparent Ar-40/Ar-39 ages than the respective cores. They also show variable effects of heterogeneous excess 40 Ar incorporation and Ar loss. One acceptable inverse isochron age of 77.1 +/- 1.1 Ma for rims surrounding older cores (82.6 +/- 0.6 Ma) overlaps with the second period of core crystallization. Compared to the phengite cores, all rims have lower B and Li abundances but similar delta B-11 values (-15 parts per thousand to -9 parts per thousand), reflecting internal redistribution of B and Li and internal fluid buffering of the B isotopic composition during rim growth. The combined observation of younger Ar-40/Ar-39 ages and boron loss, yielding comparable values of both parameters only in cores and rims of different samples, is best explained by a selective metasomatic overprint. In low permeability samples, this overprint caused recrystallization of phengite rims, whereas higher permeability in other samples led to complete recrystallization of phengite grains.
Strongly deformed samples from a several km long, blueschist-facies shear zone contain mylonitic phengite that forms a tightly clustered group of relatively young apparent Ar-40/Ar-39 ages (64.7-68.8 Ma), yielding an inverse isochron age of 65.0 +/- 3.0 Ma. Almost complete B and Li removal in mylonitic phengite is due to leaching into a fluid. The B isotopic composition is significantly heavier than in phengites from the weakly deformed samples, indicating an external control by a high-delta B-11 fluid (delta B-11 = + 7 +/- 4 parts per thousand). We interpret this result as reflecting phengite recrystallization related to deformation and associated fluid flow in the shear zone. This event also caused partial resetting of the Ar isotope system and further B loss in more permeable rocks of the adjacent unit. We conclude that geochemical evidence for pervasive or limited fluid flow is crucial for the interpretation of Ar-40/Ar-39 data in partially metasomatized rocks.
In aquatic environments, community dynamics of bacteria, especially actively growing bacteria (AGB), are tightly linked with dissolved organic matter (DOM) quantity and quality. We analyzed the community dynamics of DNA-synthesizing and accordingly AGB by linking an improved bromodeoxyuridine immunocytochemistry approach with fluorescence-activated cell sorting (BrdU-FACS). FACS-sorted cells of even oligotrophic ecosystems in winter were characterized by 16S rRNA gene analysis. In incubation experiments, we examined community shifts of AGB in response to the addition of N-acetyl-glucosamine (NAG), one of the most abundant aminosugars in aquatic systems. Our improved BrdU-FACS analysis revealed that AGB winter communities of oligotrophic Lake Stechlin (northeastern Germany) substantially differ from those of total bacteria and consist of Alpha-, Beta-, Gamma-, Deltaproteobacteria, Actinobacteria, Candidatus OP10 and Chloroflexi. AGB populations with different BrdU-fluorescence intensities and cell sizes represented different phylotypes suggesting that single-cell growth potential varies at the taxon level. NAG incubation experiments demonstrated that a variety of widespread taxa related to Alpha-, Beta-, Gammaproteobacteria, Bacteroidetes, Actinobacteria, Firmicutes, Planctomycetes, Spirochaetes, Verrucomicrobia and Chloroflexi actively grow in the presence of NAG. The BrdU-FACS approach enables detailed phylogenetic studies of AGB and, thus, to identify those phylotypes which are potential key players in aquatic DOM cycling.
Identification and characterization of growing large-scale en-echelon fractures in a salt mine
(2014)
The spatiotemporal seismicity of acoustic emission (AE) events recorded in the Morsleben salt mine is investigated. Almost a year after backfilling of the cavities from 2003, microevents are distributed with distinctive stripe shapes above cavities at different depth levels. The physical forces driving the creation of these stripes are still unknown. This study aims to find the active stripes and track fracture developments over time by combining two different temporal and spatial clustering techniques into a single methodological approach. Anomalous seismicity parameters values like sharp b-value changes for two active stripes are good indicators to explain possible stress accumulation at the stripe tips. We identify the formation of two new seismicity stripes and show that the AE activities in active clusters are migrated mostly unidirectional to eastward and upward. This indicates that the growth of underlying macrofractures is controlled by the gradient of extensional stress. Studying size distribution characteristic in terms of frequency-magnitude distribution and b-value in active phase and phase with constant seismicity rate show that deviations from the Gutenberg-Richter power law can be explained by the inclusion of different activity phases: (1) the inactive period before the formation of macrofractures, which is characterized by a deficit of larger events (higher b-values) and (2) the period of fracture growth characterized by the occurrence of larger events (smaller b-values).
Marine seismology usually relies on temporary deployments of stand alone seismic ocean bottom stations (OBS), which are initialized and synchronized on ship before deployment and re-synchronized and stopped on ship after recovery several months later. In between, the recorder clocks may drift and float at unknown rates. If the clock drifts are large or not linear and cannot be corrected for, seismological applications will be limited to methods not requiring precise common timing. Therefore, for example, array seismological methods, which need very accurate timing between individual stations, would not be applicable for such deployments.
We use an OBS test-array of 12 stations and 75 km aperture, deployed for 10 months in the deep sea (4.5-5.5 km) of the mid-eastern Atlantic. The experiment was designed to analyse the potential of broad-band array seismology at the seafloor. After recovery, we identified some stations which either show unusual large clock drifts and/or static time offsets by having a large difference between the internal clock and the GPS-signal (skew).
We test the approach of ambient noise cross-correlation to synchronize clocks of a deep water OBS array with km-scale interstation distances. We show that small drift rates and static time offsets can be resolved on vertical components with a standard technique. Larger clock drifts (several seconds per day) can only be accurately recovered if time windows of one input trace are shifted according to the expected drift between a station pair before the cross-correlation. We validate that the drifts extracted from the seismometer data are linear to first order. The same is valid for most of the hydrophones. Moreover, we were able to determine the clock drift at a station where no skew could be measured. Furthermore, we find that instable apparent drift rates at some hydrophones, which are uncorrelated to the seismometer drift recorded at the same digitizer, indicate a malfunction of the hydrophone.
Crustal earthquake swarms are an expression of intensive cracking and rock damaging over periods of days, weeks or month in a small source region in the crust. They are caused by longer lasting stress changes in the source region. Often, the localized stressing of the crust is associated with fluid or gas migration, possibly in combination with pre-existing zones of weaknesses. However, verifying and quantifying localized fluid movement at depth remains difficult since the area affected is small and geophysical prospecting methods often cannot reach the required resolution.
We apply a simple and robust method to estimate the velocity ratio between compressional (P) and shear (S) waves (upsilon(P)/upsilon(S)-ratio) in the source region of an earthquake swarm. The upsilon(P)/upsilon(S)-ratio may be unusual small if the swarm is related to gas in a porous or fractured rock. The method uses arrival time difference between P and S waves observed at surface seismic stations, and the associated double differences between pairs of earthquakes. An advantage is that earthquake locations are not required and the method seems lesser dependent on unknown velocity variations in the crust outside the source region. It is, thus, suited for monitoring purposes.
Applications comprise three natural, mid-crustal (8-10 km) earthquake swarms between 1997 and 2008 from the NW-Bohemia swarm region. We resolve a strong temporal decrease of upsilon(P)/upsilon(S) before and during the main activity of the swarm, and a recovery of upsilon(P)/upsilon(S) to background levels at the end of the swarms. The anomalies are interpreted in terms of the Biot-Gassman equations, assuming the presence of oversaturated fluids degassing during the beginning phase of the swarm activity.
The Ceres earthquake of 29 September 1969 is the largest known earthquake in southern Africa. Digitized analog recordings from Worldwide Standardized Seismographic Network stations (Powell and Fries, 1964) are used to retrieve the point source moment tensor and the most likely centroid depth of the event using full waveform modeling. A scalar seismic moment of 2.2-2.4 x 10(18) N center dot m corresponding to a moment magnitude of 6.2-6.3 is found. The analysis confirms the pure strike-slip mechanism previously determined from onset polarities by Green and Bloch (1971). Overall good agreement with the fault orientation previously estimated from local aftershock recordings is found. The centroid depth can be constrained to be less than 15 km. In a second analysis step, we use a higher order moment tensor based inversion scheme for simple extended rupture models to constrain the lateral fault dimensions. We find rupture propagated unilaterally for 4.7 s from east-southwest to west-northwest for about 17 km ( average rupture velocity of about 3: 1 km/s).
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.
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.
Numerical simulation of fluid-flow processes in a 3D high-resolution carbonate reservoir analogue
(2014)
A high-resolution three-dimensional (3D) outcrop model of a Jurassic carbonate ramp was used in order to perform a series of detailed and systematic flow simulations. The aim of this study was to test the impact of small- and large-scale geological features on reservoir performance and oil recovery. The digital outcrop model contains a wide range of sedimentological, diagenetic and structural features, including discontinuity surfaces, shoal bodies, mud mounds, oyster bioherms and fractures. Flow simulations are performed for numerical well testing and secondary oil recovery. Numerical well testing enables synthetic but systematic pressure responses to be generated for different geological features observed in the outcrops. This allows us to assess and rank the relative impact of specific geological features on reservoir performance. The outcome documents that, owing to the realistic representation of matrix heterogeneity, most diagenetic and structural features cannot be linked to a unique pressure signature. Instead, reservoir performance is controlled by subseismic faults and oyster bioherms acting as thief zones. Numerical simulations of secondary recovery processes reveal strong channelling of fluid flow into high-permeability layers as the primary control for oil recovery. However, appropriate reservoir-engineering solutions, such as optimizing well placement and injection fluid, can reduce channelling and increase oil recovery.
Carbonate reservoirs pose significant challenges for reservoir modelling and flow prediction due to heterogeneities in rock properties, limits to seismic resolution and limited constraints on subsurface data. Hence, a systematic and streamlined approach is needed to construct geological models and to quickly evaluate key sensitivities in the flow models. This paper discusses results from a reservoir analogue study of a Middle Jurassic carbonate ramp in the High Atlas Mountains of Morocco that has stratigraphic and structural similarities to selected Middle East reservoirs. For this purpose, high-resolution geological models were constructed from the integration of sedimentological, diagenetic and structural studies in the area. The models are approximately 1200 x 1250 m in size, and only faults (no fractures) with offsets greater than 1 m are included. Novel methods have been applied to test the response of flow simulations to the presence or absence of specific geological features, including proxies for hardgrounds, stylolites, patch reefs, and mollusc banks, as a way to guide the level of detail that is suitable for modelling objectives. Our general conclusion from the study is that the continuity of any geological feature with extreme permeability (high or low) has the most significant impact on flow.
Whilst sophisticated multiphase fluid flow models are routinely employed to understand behaviour of oil and gas reservoirs, high-resolution data describing the three-dimensional (3D) distribution of rock characteristics is rarely available to populate models. We present a new approach to developing a quantitative understanding of the effect of individual controls on the distribution of petrophysical properties and their impact on fluid flow. This involves simulating flow through high-detail permeability architectures generated by forward modelling of the coupled depositional-diagenetic evolution of isolated platforms using CARB3D(+). This workflow is exemplified by an investigation of interactions between subsidence and climate, and their expression in spatial variations in reservoir quality in an isolated carbonate platform of similar size and subsidence history to the Triassic Latemar Platform.
Dissolutional lowering during subaerial exposure controls platform-top graininess via platform top hydrodynamics during the subsequent transgression. Dissolved carbonate is reprecipitated as cements by percolating meteoric waters. However, associated subsurface meteoric dissolution generates significant secondary porosity under a more humid climate. Slower subsidence enhances diagenetic overprinting during repeated exposure events. Single-phase streamline simulations show how early diagenesis develops more permeable fairways within the finer-grained condensed units that can act as thief zones for flow from the grainier but less diagenetically altered cyclic units.
Lacustrine sediments have been widely used to investigate past climatic and environmental changes on millennial to seasonal time scales. Sedimentary archives of lakes in mountainous regions may also record non-climatic events such as earthquakes. We argue herein that a set of 64 annual laminae couplets reconciles a stratigraphically inconsistent accelerator mass spectrometry (AMS) C-14 chronology in a similar to 4-m-long sediment core from Lake Mengda, in the north-eastern Tibetan Plateau. The laminations suggest the lake was formed by a large landslide, triggered by the 1927 Gulang earthquake (M = 8.0). The lake sediment sequence can be separated into three units based on lithologic, sedimentary, and isotopic characteristics. Starting from the bottom of the sequence, these are: (1) unweathered, coarse, sandy valley-floor deposits or landslide debris that pre-date the lake, (2) landslide-induced, fine-grained soil or reworked landslide debris with a high organic content, and (3) lacustrine sediments with low organic content and laminations. These annual laminations provide a high-resolution record of anthropogenic and environmental changes during the twentieth century, recording enhanced sediment input associated with two phases of construction activities. The high mean sedimentation rates of up to 4.8 mm year(-1) underscore the potential for reconstructing such distinct sediment pulses in remote, forested, and seemingly undisturbed mountain catchments.
The Tibetan Plateau is a region that is highly sensitive to recent global warming, but the complexity and heterogeneity of its mountainous landscape can result in variable responses. In addition, the scarcity and brevity of regional instrumental and palaeoecological records still hamper our understanding of past and present patterns of environmental change. To investigate how the remote, high-alpine environments of the Nianbaoyeze Mountains, eastern Tibetan Plateau, are affected by climate change and human activity over the last similar to 600 years, we compared regional tree-ring studies with pollen and diatom remains archived in the dated sediments of Dongerwuka Lake (33.22A degrees N, 101.12A degrees E, 4,307 m a.s.l.). In agreement with previous studies from the eastern Tibetan Plateau, a strong coherence between our two juniper-based tree-ring chronologies from the Nianbaoyeze and the Anemaqin Mountains was observed, with pronounced cyclical variations in summer temperature reconstructions. A positive directional trend to warmer summer temperatures in the most recent decades, was, however, not observed in the tree-ring record. Likewise, our pollen and diatom spectra showed minimal change over the investigated time period. Although modest, the most notable change in the diatom relative abundances was a subtle decrease in the dominant planktonic Cyclotella ocellata and a concurrent increase in small, benthic fragilarioid taxa in the similar to 1820s, suggesting higher ecosystem variability. The pollen record subtly indicates three periods of increased cattle grazing activity (similar to 1400-1480 AD, similar to 1630-1760 AD, after 1850 AD), but shows generally no significant vegetation changes during past similar to 600 years. The minimal changes observed in the tree-ring, diatom and pollen records are consistent with the presence of localised cooling centres that are evident in instrumental and tree-ring data within the southeastern and eastern Tibetan Plateau. Given the minor changes in regional temperature records, our complacent palaeoecological profiles suggest that climatically induced ecological thresholds have not yet been crossed in the Nianbaoyeze Mountains region.
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.
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.
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.