Refine
Year of publication
- 2011 (141) (remove)
Document Type
- Article (111)
- Doctoral Thesis (24)
- Review (2)
- Habilitation Thesis (1)
- Other (1)
- Postprint (1)
- Preprint (1)
Keywords
- Holocene (6)
- Tibetan Plateau (6)
- Site effects (4)
- Holozän (3)
- Pollen (3)
- Wave propagation (3)
- climate change (3)
- Benzene (2)
- Climate change (2)
- Crustal structure (2)
Institute
- Institut für Geowissenschaften (141) (remove)
P>Computing the magnitude of an earthquake requires correcting for the propagation effects from the source to the receivers. This is often accomplished by performing numerical simulations using a suitable Earth model. In this work, the energy magnitude M(e) is considered and its determination is performed using theoretical spectral amplitude decay functions over teleseismic distances based on the global Earth model AK135Q. Since the high frequency part (above the corner frequency) of the source spectrum has to be considered in computing M(e), the influence of propagation and site effects may not be negligible and they could bias the single station M(e) estimations. Therefore, in this study we assess the inter- and intrastation distributions of errors by considering the M(e) residuals computed for a large data set of earthquakes recorded at teleseismic distances by seismic stations deployed worldwide. To separate the inter- and intrastation contribution of errors, we apply a maximum likelihood approach to the M(e) residuals. We show that the interstation errors (describing a sort of site effect for a station) are within +/- 0.2 magnitude units for most stations and their spatial distribution reflects the expected lateral variation affecting the velocity and attenuation of the Earth's structure in the uppermost layers, not accounted for by the 1-D AK135Q model. The variance of the intrastation error distribution (describing the record-to-record component of variability) is larger than the interstation one (0.240 against 0.159), and the spatial distribution of the errors is not random but shows specific patterns depending on the source-to-station paths. The set of coefficients empirically determined may be used in the future to account for the heterogeneities of the real Earth not considered in the theoretical calculations of the spectral amplitude decay functions used to correct the recorded data for propagation effects.
Swedish long-term soil fertility experiments were used to investigate the effect of texture and fertilization regime on soil electrical conductivity. In one geophysical approach, fields were mapped to characterize the horizontal variability in apparent electrical conductivity down to 1.5 m soil depth using an electromagnetic induction meter (EM38 device). The data obtained were geo-referenced by dGPS. The other approach consisted of measuring the vertical variability in electrical conductivity along transects using a multi-electrode apparatus for electrical resistivity tomography (GeoTom RES/IP device) down to 2 m depth. Geophysical field work was complemented by soil analyses. The results showed that despite 40 years of different fertilization regimes, treatments had no significant effects on the apparent electrical conductivity. Instead, the comparison of sites revealed high and low conductivity soils, with gradual differences explained by soil texture. A significant, linear relationship found between apparent electrical conductivity and soil clay content explained 80% of the variability measured. In terms of soil depth, both low and high electrical conductivity values were measured. Abrupt changes in electrical conductivity within a field revealed the presence of 'deviating areas'. Higher values corresponded well with layers with a high clay content, while local inclusions of coarse-textured materials caused a high variability in conductivity in some fields. The geophysical methods tested provided useful information on the variability in soil texture at the experimental sites. The use of spatial EC variability as a co-variable in statistical analysis could be a complementary tool in the evaluation of experimental results.
The Western Alpine Sesia-Lanzo Zone (SLZ) is a sliver of eclogite-facies continental crust exhumed from mantle depths in the hanging wall of a subducted oceanic slab. Eclogite-facies felsic and basic rocks sampled across the internal SLZ show different degrees of retrograde metamorphic overprint associated with fluid influx. The weakly deformed samples preserve relict eclogite-facies mineral assemblages that show partial fluid-induced compositional re-equilibration along grain boundaries, brittle fractures and other fluid pathways. Multiple fluid influx stages are indicated by replacement of primary omphacite by phengite, albitic plagioclase and epidote as well as partial re-equilibration and/or overgrowths in phengite and sodic amphibole, producing characteristic step-like compositional zoning patterns. The observed textures, together with the map-scale distribution of the samples, suggest open-system, pervasive and reactive fluid flux across large rock volumes above the subducted slab. Thermodynamic modelling indicates a minimum amount of fluid of 0 center dot 1-0 center dot 5 wt % interacting with the wall-rocks. Phase relations and reaction textures indicate mobility of K, Ca, Fe and Mg, whereas Al is relatively immobile in these medium-temperature-high-pressure fluids. Furthermore, the thermodynamic models show that recycling of previously fractionated material, such as in the cores of garnet porphyroblasts, largely controls the compositional re-equilibration of the exhumed rock body.
The amount and composition of subduction zone fluids and the effect of fluid-rock interaction at a slab-mantle interface have been constrained by thermodynamic and trace element modelling of partially overprinted blueschist-facies rocks from the Sesia Zone (Western Alps). Deformation-induced differences in fluid flux led to a partial preservation of pristine mineral cores in weakly deformed samples that were used to quantify Li, B, Stand Pb distribution during mineral growth, -breakdown and modification induced by fluid-rock interaction. Our results show that Li and 13 budgets are fluid-controlled, thus acting as tracers for fluid-rock interaction processes, whereas Stand Pb budgets are mainly controlled by the fluid-induced formation of epidote. Our calculations show that fluid-rock interaction caused significant Li and B depletion in the affected rocks due to leaching effects, which in turn can lead to a drastic enrichment of these elements in the percolating fluid. Depending on available fluid-mineral trace element distribution coefficients modelled fluid rock ratios were up to 0.06 in weakly deformed samples and at least 0.5 to 4 in shear zone mylonites. These amounts lead to time integrated fluid fluxes of up to 1.4-10(2) m(3) m(-2) in the weakly deformed rocks and 1-8-10(3) m(3) m(-2) in the mylonites. Combined thermodynamic and trace element models can be used to quantify metamorphic fluid fluxes and the associated element transfer in complex, reacting rock systems and help to better understand commonly observed fluid-induced trace element trends in rocks and minerals from different geodynamic environments.
Eclogites from the main borehole of the Chinese Continental Scientific Drilling project yield highly precise Lu-Hf garnet-clinopyroxene ages of 216.9 +/- 1.2 Ma (four samples) and 220.5 +/- 2.7 Ma (one sample). The spatial distribution of the rare earth elements in garnet is consistent with the preservation of primary growth zoning, unmodified by diffusion, which supports the interpretation that the Lu-Hf ages date the time of formation of garnet, the major rock forming mineral in the eclogites. The preservation of primary REE-zoning, despite peak metamorphic temperatures around 800-850 degrees C. indicates that the Lu-Hf chronometer is perfectly suitable to date garnet-forming reactions in high grade rocks. The range of Lu-Hf ages for eclogites in the Dabie-Sulu UHP terrane point to episodic rather than continuous growth of garnets and thus punctuated metamorphism during the collision of the North China Block and the Yangtze Block. The U-Pb ages and Hf-isotope systematics of zircon grains from one eclogite sample imply a protracted geologic history of the eclogite precursors that started around 2 Ga and culminated in the UHP metamorphism around 220 Ma.
The pressures required for diamond and coesite formation far exceed conditions reached by even the deepest present-day orogenic crustal roots. Therefore the occurrence of metamorphosed continental crust containing these minerals requires processes other than crustal thickening to have operated in the past. Here we report the first in situ finding of diamond and coesite, characterized by micro-Raman spectroscopy, in high-pressure granulites otherwise indistinguishable from granulites found associated with garnet peridotite throughout the European Variscides. Our discovery confirms the provenance of Europe's first reliable diamond, the "Bohemian diamond," found in A.D. 1870, and also represents the first robust evidence for ultrahigh-pressure conditions in a major Variscan crustal rock type. A process of deep continental subduction is required to explain the metamorphic pressures and the granulite-garnet peridotite association, and thus tectonometamorphic models for these rocks involving a deep orogenic crustal root need to be significantly modified.
Die vorliegende Arbeit basiert auf Forschungen in den Jahren 2007-2009. Sie betrachtet die saisonale Arbeitsmigration aus der polnischen Region Konin, wo die Arbeitsmigration aus ökonomischen Gründen, wie auch in ähnlich strukturierten Gebieten Polens, eine lange Tradition hat, die bis ins 19. Jahrhundert zurückgeht. Sie wird die saisonale Migration ins Ausland mit den ökonomischen, sozialen und räumlichen Auswirkungen aus der Perspektive des Einzelnen und seiner unmittelbaren Umgebung, aber auch der Gesellschaft und Herkunftsgebiet der Migranten betrachtet.
Rainfall, snow-, and glacial melt throughout the Himalaya control river discharge, which is vital for maintaining agriculture, drinking water and hydropower generation. However, the spatiotemporal contribution of these discharge components to Himalayan rivers is not well understood, mainly because of the scarcity of ground-based observations. Consequently, there is also little known about the triggers and sources of peak sediment flux events, which account for extensive hydropower reservoir filling and turbine abrasion. We therefore lack basic information on the distribution of water resources and controls of erosion processes. In this thesis, I employ various methods to assess and quantify general characteristics of and links between precipitation, river discharge, and sediment flux in the Sutlej Valley. First, I analyze daily precipitation data (1998-2007) from 80 weather stations in the western Himalaya, to decipher the distribution of rain- and snowfall. Rainfall magnitude frequency analyses indicate that 40% of the summer rainfall budget is attributed to monsoonal rainstorms, which show higher variability in the orogenic interior than in frontal regions. Combined analysis of rainstorms and sediment flux data of a major Sutlej River tributary indicate that monsoonal rainfall has a first order control on erosion processes in the orogenic interior, despite the dominance of snowfall in this region. Second, I examine the contribution of rainfall, snow and glacial melt to river discharge in the Sutlej Valley (s55,000 km2), based on a distributed hydrological model, which covers the period 2000-2008. To achieve high spatial and daily resolution despite limited ground-based observations the hydrological model is forced by daily remote sensing data, which I adjusted and calibrated with ground station data. The calibration shows that the Tropical Rainfall Measuring Mission (TRMM) 3B42 rainfall product systematically overestimates rainfall in semi-arid and arid regions, increasing with aridity. The model results indicate that snowmelt-derived discharge (74%) is most important during the pre-monsoon season (April to June) whereas rainfall (56%) and glacial melt (17%) dominate the monsoon season (July-September). Therefore, climate change most likely causes a reduction in river discharge during the pre-monsoon season, which especially affects the orogenic interior. Third, I investigate the controls on suspended sediment flux in different parts of the Sutlej catchments, based on daily gauging data from the past decade. In conjunction with meteorological data, earthquake records, and rock strength measurements I find that rainstorms are the most frequent trigger of high-discharge events with peaks in suspended sediment concentrations (SSC) that account for the bulk of the suspended sediment flux. The suspended sediment flux increases downstream, mainly due to increases in runoff. Pronounced erosion along the Himalayan Front occurs throughout the monsoon season, whereas efficient erosion of the orogenic interior is confined to single extreme events. The results of this thesis highlight the importance of snow and glacially derived melt waters in the western Himalaya, where extensive regions receive only limited amounts of monsoonal rainfall. These regions are therefore particularly susceptible to global warming with major implications on the hydrological cycle. However, the sediment discharge data show that infrequent monsoonal rainstorms that pass the orographic barrier of the Higher Himalaya are still the primary trigger of the highest-impact erosion events, despite being subordinate to snow and glacially–derived discharge. These findings may help to predict peak sediment flux events and could underpin the strategic development of preventative measures for hydropower infrastructures.
The East African Plateau provides a spectacular example of geodynamic plateau uplift, active continental rifting, and associated climatic forcing. It is an integral part of the East African Rift System and has an average elevation of approximately 1,000 m. Its location coincides with a negative Bouguer gravity anomaly with a semi-circular shape, closely related to a mantle plume, which influences the Cenozoic crustal development since its impingement in Eocene-Oligocene time. The uplift of the East African Plateau, preceding volcanism, and rifting formed an important orographic barrier and tectonically controlled environment, which is profoundly influenced by climate driven processes. Its location within the equatorial realm supports recently proposed hypotheses, that topographic changes in this region must be considered as the dominant forcing factor influencing atmospheric circulation patterns and rainfall distribution. The uplift of this region has therefore often been associated with fundamental climatic and environmental changes in East Africa and adjacent regions. While the far-reaching influence of the plateau uplift is widely accepted, the timing and the magnitude of the uplift are ambiguous and are still subject to ongoing discussion. This dilemma stems from the lack of datable, geomorphically meaningful reference horizons that could record surface uplift. In order to quantify the amount of plateau uplift and to find evidence for the existence of significant relief along the East African Plateau prior to rifting, I analyzed and modeled one of the longest terrestrial lava flows; the 300-km-long Yatta phonolite flow in Kenya. This lava flow is 13.5 Ma old and originated in the region that now corresponds to the eastern rift shoulders. The phonolitic flow utilized an old riverbed that once drained the eastern flank of the plateau. Due to differential erosion this lava flow now forms a positive relief above the parallel-flowing Athi River, which is mimicking the course of the paleo-river. My approach is a lava-flow modeling, based on an improved composition and temperature dependent method to parameterize the flow of an arbitrary lava in a rectangular-shaped channel. The essential growth pattern is described by a one-dimensional model, in which Newtonian rheological flow advance is governed by the development of viscosity and/or velocity in the internal parts of the lava-flow front. Comparing assessments of different magma compositions reveal that length-dominated, channelized lava flows are characterized by high effusion rates, rapid emplacement under approximately isothermal conditions, and laminar flow. By integrating the Yatta lava flow dimensions and the covered paleo-topography (slope angle) into the model, I was able to determine the pre-rift topography of the East African Plateau. The modeling results yield a pre-rift slope of at least 0.2°, suggesting that the lava flow must have originated at a minimum elevation of 1,400 m. Hence, high topography in the region of the present-day Kenya Rift must have existed by at least 13.5 Ma. This inferred mid-Miocene uplift coincides with the two-step expansion of grasslands, as well as important radiation and speciation events in tropical Africa. Accordingly, the combination of my results regarding the Yatta lava flow emplacement history, its location, and its morphologic character, validates it as a suitable “paleo-tiltmeter” and has thus to be considered as an important topographic and volcanic feature for the topographic evolution in East Africa.
In dieser Arbeit wird das regionale Klimamodell HIRHAM mit einer horizontalen Auflösung von 50 km und 19 vertikalen Schichten erstmals auf den asiatischen Kontinent angewendet, um die indische Monsunzirkulation unter rezenten und paläoklimatischen Bedingungen zu simulieren. Das Integrationsgebiet des Modells erstreckt sich von etwa 0ºN - 50ºN und 42ºE - 110ºE und bedeckt dabei sowohl die hohe Topographie des Himalajas und Tibet Plateaus als auch den nördlichen Indischen Ozean. Das Ziel besteht in der Beschreibung der regionalen Kopplung zwischen der Monsunzirkulation und den orographischen sowie diabatischen Antriebsmechanismen. Eine 44-jährige Modellsimulation von 1958-2001, die am seitlichen und unteren Rand von ECMWF Reanalysen (ERA40) angetrieben wird, bildet die Grundlage für die Validierung der Modellergebnisse mit Beobachtungen auf der Basis von Stations- und Gitterdatensätzen. Der Fokus liegt dabei auf der atmosphärischen Zirkulation, der Temperatur und dem Niederschlag im Sommer- und Wintermonsun, wobei die Qualität des Modells sowohl in Bezug zur langfristigen und dekadischen Klimatologie als auch zur interannuellen Variabilität evaluiert wird. Im Zusammenhang mit einer realistischen Reproduktion der Modelltopographie kann für die Muster der Zirkulation und Temperatur eine gute Übereinstimmung zwischen Modell und Daten nachgewiesen werden. Der simulierte Niederschlag zeigt eine bessere Übereinstimmung mit einem hoch aufgelösten Gitterdatensatz über der Landoberfläche Zentralindiens und in den Hochgebirgsregionen, der den Vorteil des Regionalmodells gegenüber der antreibenden Reanalyse hervorhebt. In verschiedenen Fall- und Sensitivitätsstudien werden die wesentlichen Antriebsfaktoren des indischen Monsuns (Meeresoberflächentemperaturen, Stärke des winterlichen Sibirischen Hochs und Anomalien der Bodenfeuchte) untersucht. Die Ergebnisse machen deutlich, dass die Simulation dieser Mechanismen auch mit einem Regionalmodell sehr schwierig ist, da die Komplexität des Monsunsystems hochgradig nichtlinear ist und die vor allem subgridskalig wirkenden Prozesse im Modell noch nicht ausreichend parametrisiert und verstanden sind. Ein paläoklimatisches Experiment für eine 44-jährige Zeitscheibe im mittleren Holozän (etwa 6000 Jahre vor heute), die am Rand von einer globalen ECHAM5 Simulation angetrieben wird, zeigt markante Veränderungen in der Intensität des Monsuns durch die unterschiedliche solare Einstrahlung, die wiederum Einflüsse auf die SST, die Zirkulation und damit auf die Niederschlagsmuster hat.
Spatial and temporal temperature and moisture patterns across the Tibetan Plateau are very complex. The onset and magnitude of the Holocene climate optimum in the Asian monsoon realm, in particular, is a subject of considerable debate as this time period is often used as an analogue for recent global warming. In the light of contradictory inferences regarding past climate and environmental change on the Tibetan Plateau, I have attempted to explain mismatches in the timing and magnitude of change. Therefore, I analysed the temporal variation of fossil pollen and diatom spectra and the geochemical record from palaeo-ecological records covering different time scales (late Quaternary and the last 200 years) from two core regions in the NE and SE Tibetan Plateau. For interpretation purposes I combined my data with other available palaeo-ecological data to set up corresponding aquatic and terrestrial proxy data sets of two lake pairs and two sets of sites. I focused on the direct comparison of proxies representing lacustrine response to climate signals (e.g., diatoms, ostracods, geochemical record) and proxies representing changes in the terrestrial environment (i.e., terrestrial pollen), in order to asses whether the lake and its catchments respond at similar times and magnitudes to environmental changes. Therefore, I introduced the established numerical technique procrustes rotation as a new approach in palaeoecology to quantitatively compare raw data of any two sedimentary records of interest in order to assess their degree of concordance. Focusing on the late Quaternary, sediment cores from two lakes (Kuhai Lake 35.3°N; 99.2°E; 4150 m asl; and Koucha Lake 34.0°N; 97.2°E; 4540 m asl) on the semi-arid northeastern Tibetan Plateau were analysed to identify post-glacial vegetation and environmental changes, and to investigate the responses of lake ecosystems to such changes. Based on the pollen record, five major vegetation and climate changes could be identified: (1) A shift from alpine desert to alpine steppe indicates a change from cold, dry conditions to warmer and more moist conditions at 14.8 cal. ka BP, (2) alpine steppe with tundra elements points to conditions of higher effective moisture and a stepwise warming climate at 13.6 cal. ka BP, (3) the appearance of high-alpine meadow vegetation indicates a further change towards increased moisture, but with colder temperatures, at 7.0 cal. ka BP, (4) the reoccurrence of alpine steppe with desert elements suggests a return to a significantly colder and drier phase at 6.3 cal. ka BP, and (5) the establishment of alpine steppe-meadow vegetation indicates a change back to relatively moist conditions at 2.2 cal. ka BP. To place the reconstructed climate inferences from the NE Tibetan Plateau into the context of Holocene moisture evolution across the Tibetan Plateau, I applied a five-scale moisture index and average link clustering to all available continuous pollen and non-pollen palaeoclimate records from the Tibetan Plateau, in an attempt to detect coherent regional and temporal patterns of moisture evolution on the Plateau. However, no common temporal or spatial pattern of moisture evolution during the Holocene could be detected, which can be assigned to the complex responses of different proxies to environmental changes in an already very heterogeneous mountain landscape, where minor differences in elevation can result in marked variations in microenvironments. Focusing on the past 200 years, I analysed the sedimentary records (LC6 Lake 29.5°N, 94.3°E, 4132 m asl; and Wuxu Lake 29.9°N, 101.1°E, 3705 m asl) from the southeastern Tibetan Plateau. I found that despite presumed significant temperature increases over that period, pollen and diatom records from the SE Tibetan Plateau reveal only very subtle changes throughout their profiles. The compositional species turnover investigated over the last 200 years appears relatively low in comparison to the species reorganisations during the Holocene. The results indicate that climatically induced ecological thresholds are not yet crossed, but that human activity has an increasing influence, particularly on the terrestrial ecosystem. Forest clearances and reforestation have not caused forest decline in our study area, but a conversion of natural forests to semi-natural secondary forests. The results from the numerical proxy comparison of the two sets of two pairs of Tibetan lakes indicate that the use of different proxies and the work with palaeo-ecological records from different lake types can cause deviant stories of inferred change. Irrespective of the timescale (Holocene or last 200 years) or region (SE or NE Tibetan Plateau) analysed, the agreement in terms of the direction, timing, and magnitude of change between the corresponding terrestrial data sets is generally better than the match between the corresponding lacustrine data sets, suggesting that lacustrine proxies may partly be influenced by in-lake or local catchment processes whereas the terrestrial proxy reflects a more regional climatic signal. The current disaccord on coherent temporal and spatial climate patterns on the Tibetan Plateau can partly be ascribed to the complexity of proxy response and lake systems on the Tibetan Plateau. Therefore, a multi-proxy, multi-site approach is important in order to gain a reliable climate interpretation for the complex mountain landscape of the Tibetan Plateau.
Lithopheric rheology at the Dead Sea Transform Fault constrained by thermo-mechanical modeling
(2011)
We use Global Positioning System (GPS) velocities and kinematic Finite Element models (FE-models) to infer the state of locking between the converging Nazca and South America plates in South-Central Chile (36 degrees S -46 degrees S) and to evaluate its spatial and temporal variability. GPS velocities provide information on earthquake-cycle deformation over the last decade in areas affected by the megathrust events of 1960 (M-w = 9.5) and 2010 (M-w = 8.8). Our data confirm that a change in surface velocity patterns of these two seismotectonic segments can be related to their different stages in the seismic cycle: Accordingly, the northern (2010) segment was in a final stage of interseismic loading whereas the southern (1960) segment is still in a postseismic stage and undergoes a prolonged viscoelastic mantle relaxation. After correcting the signals for mantle relaxation, the residual GPS velocity pattern suggests that the plate interface accumulates slip deficit in a spatially and presumably temporally variable way towards the next great event. Though some similarity exist between locking and 1960 coseismic slip, extrapolating the current, decadal scale slip deficit accumulation towards the similar to 300-yr recurrence times of giant events here does neither yield the slip distribution nor the moment magnitude of the 1960 earthquake. This suggests that either the locking pattern is evolving in time (to reconcile a slip deficit distribution similar to the 1960 earthquake) or that some asperities are not persistent over multiple events. The accumulated moment deficit since 1960 suggests that highly locked patches in the 1960 segment are already capable of producing a M similar to 8 event if triggered to fail by stress transfer from the 2010 event.
Secondary forests are gaining increased importance in tropical landscapes and have recently been reported to act as potential belowground carbon sinks. While economic interest in the management of secondary forests to mitigate carbon emissions is rising, the dynamics of soil carbon stocks under these ecosystems remain poorly understood. Recent studies report conflicting results concerning soil carbon trends as well as multiple confounding factors (e.g. soil type, topography and land-use history) affecting these trends. In this study, organic carbon stocks were measured in the mineral soil up to 20 cm depth of at 24 active pastures, 5-8-year-old, and 12-15-year-old secondary forest sites on former pastures. Additionally, we estimated carbon stocks under a 100-year-old secondary forest and compared them to those of nearby mature forests. Abiotic conditions in the study area were homogenous, enabling us to isolate the effect of land-use change on soil organic carbon stocks. Contrary to our expectations, soil carbon stocks in the top 10 cm did not change with young secondary forest development. Pasture soils stored 24.8 +/- 2.9 Mg ha(-1) carbon (mean +/- standard error) in the top 10 cm, and no accumulation of soil carbon was apparent during the first 15 years of secondary succession. Soil carbon stocks under 100-year-old secondary forests, averaging 43.0 +/- 7.9 Mg ha(-1) (mean +/- standard error), were clearly higher than those recorded at younger sites and approached levels of soil carbon stocks under mature forests. These data indicate that soil carbon stocks in this region of Panama are not affected by the land-use transition from pasture to young secondary regrowth. However, an increase of soil carbon storage might be possible over a longer period of time. Our results support trends observed in other tropical areas and highlight the importance of environmental conditions such as soil properties rather than land-use transitions on soil carbon dynamics. While our understanding of organic carbon dynamics in tropical soils remains limited, these results underscore the challenges of undertaking short-term reforestation projects with the expectation of increasing soil carbon sequestration.
Analyses of radiogenic neodymium (Nd), strontium (Sr), and lead (Pb) isotope compositions of clay-sized detrital sediments allow detailed tracing of source areas of sediment supply and present and past transport of particles by water masses in the eastern Indian Ocean. Isotope signatures in surface sediments range from -21.5 (epsilon Nd), 0.8299 ((87)Sr/(86S)r), and 19.89 ((206)Pb/(204)Pb) off northwest Australia to + 0.7 (epsilon Nd), 0.7069 ((87)Sr/(86)Sr), and 17.44 ((206)Pb/(204)Pb) southwest of Java. The radiogenic isotope signatures primarily reflect petrographic characteristics of the surrounding continental bedrocks but are also influenced by weathering-induced grain size effects of Pb and Sr isotope systems with superimposed features that are caused by current transport of clay-sized particles, as evidenced off Australia where a peculiar isotopic signature characterizes sediments underlying the southward flowing Leeuwin Current and the northward flowing West Australian Current (WAC). Gravity core FR10/95-GC17 off west Australia recorded a major isotopic change from Last Glacial Maximum values of -10 (epsilon Nd), 0.745 ((87S)r/(86)Sr), and 18.8 ((206)Pb/(204)Pb) to Holocene values of -22 (epsilon Nd), 0.8 ((87)Sr/(86)Sr), and 19.3 ((206)Pb/(204)Pb), which documents major climatically driven changes of the WAC and in local riverine particle supply from Australia during the past 20 kyr. In contrast, gravity core FR10/95-GC5 located below the present-day pathway of the Indonesian throughflow (ITF) shows a much smaller isotopic variability, indicating a relatively stable ITF hydrography over most of the past 92 kyr. Only the surface sediments differ significantly in their isotopic composition, indicating substantial changes in erosional sources attributed to a change of the current regime during the past 5 kyr.
Hydrologic modelers often need to know which method of quantitative precipitation estimation (QPE) is best suited for a particular catchment. Traditionally, QPE methods are verified and benchmarked against independent rain gauge observations. However, the lack of spatial representativeness limits the value of such a procedure. Alternatively, one could drive a hydrological model with different QPE products and choose the one which best reproduces observed runoff. Unfortunately, the calibration of conceptual model parameters might conceal actual differences between the QPEs. To avoid such effects, we abandoned the idea of determining optimum parameter sets for all QPE being compared. Instead, we carry out a large number of runoff simulations, confronting each QPE with a common set of random parameters. By evaluating the goodness-of-fit of all simulations, we obtain information on whether the quality of competing QPE methods is significantly different. This knowledge is inferred exactly at the scale of interest-the catchment scale. We use synthetic data to investigate the ability of this procedure to distinguish a truly superior QPE from an inferior one. We find that the procedure is prone to failure in the case of linear systems. However, we show evidence that in realistic (nonlinear) settings, the method can provide useful results even in the presence of moderate errors in model structure and streamflow observations. In a real-world case study on a small mountainous catchment, we demonstrate the ability of the verification procedure to reveal additional insights as compared to a conventional cross validation approach.
Challenges in understanding the hydrologic controls on the mobility of slow-moving landslides
(2011)
Slow-moving landslides are a wide-spread type of active mass movement, can cause severe damages to infrastructure, and may be a precursor of sudden catastrophic slope failures. Pore-water pressure is commonly regarded as the most important among a number of possible factors controlling landslide velocity. We used high-resolution monitoring data to explore the relations of landslide mobility and hydrologic processes at the Heumoser landslide in Austria, which is characterized by continuous slow movement along a shear zone. Movement rates showed a seasonality that was associated with elevated pore-water pressures. Pore pressure monitoring revealed a system of confined and separated aquifers with differing dynamics. Analysis of a simple infinite slope mobility model showed that small variations in parameters, along with measured pore pressure dynamics, provided a perfect match to our observations. Modeling showed a stabilizing effect of snow cover due to the additional load. This finding was supported by a multiple regression model, which further suggested that effective pore pressures at the slip surface were partially differing from the borehole observations and were related to preferential infiltration and subsurface flow in adjacent areas. It appears that in a setting like the Heumoser landslide, hydrologic processes delicately influence slope mobility through their control on pore pressure dynamics and the weight of the landslide body, which challenges observation and modeling. Moreover, it appears that their simplicity, and especially their high sensitivity to parameter variations, limits the conclusions that can be drawn from infinite slope models.
Landscapes in the humid tropics are undergoing a continuous change in land use. Deforestation is still taking its toll on forested areas, but at the same time more and more secondary forests emerge where formerly agricultural lands and pastures are being abandoned. Regarding soil hydrology, the extent to which secondary succession can recover soil hydrological properties disturbed by antecedent deforestation and pasture use is yet poorly understood. We investigated the effect of secondary succession on saturated hydraulic conductivity (Ks) at two soil depths (0-6 and 6-12 cm) using a space-for-time approach in a landscape mosaic in central Panama. The following four land-use classes were studied: pasture (P), secondary forest of 5-8 years of age (SF5), secondary forest of 12-15 years of age (SF12) and secondary forest of more than 100 years of age (SF100), each replicated altogether four times in different micro-catchments across the study region. The hydrological implications of differences in Ks in response to land-use change with land use, especially regarding overland flow generation, were assessed via comparisons with rainfall intensities.
Recovery of Ks could be detected in the 0-6 cm depth after 12 years of secondary succession: P and SF5 held similar Ks values, but differed significantly (alpha = 0.05) from SF12 and SF100 which in turn were indistinguishable. Variability within the land cover classes was large but, due to sufficient replication in the study, Ks recovery could be detected nonetheless. Ks in the 6-12 cm depth did not show any differences between the land cover classes; only Ks of the uppermost soil layer was affected by land-use changes. Overland flow - as inferred from comparisons of Ks with rainfall intensities - is more likely on P and SF5 sites compared to SF12 and 5E100 for the upper sample depth; however, generally low values at the 6-12 cm depth are likely to impede vertical percolation during high rainfall intensities regardless of land use.
We conclude that Ks can recover from pasture use under secondary succession up to pre-pasture levels, but the process may take more than 8 years. In order to gain comprehensive understanding of Ks change with land use and its hydrological implications, more studies with detailed land-use histories and combined measurements of Ks, overland flow, precipitation and throughfall are essential.
Species distribution models are an important tool to predict the impact of global change on species distributional ranges and community assemblages. Although considerable progress has been made in the statistical modeling during the last decade, many approaches still ignore important features of species distributions, such as nonlinearity and interactions between predictors, spatial autocorrelation, and nonstationarity, or at most incorporate only some of these features. Ecologists, however, require a modeling framework that simultaneously addresses all these features flexibly and consistently. Here we describe such an approach that allows the estimation of the global effects of environmental variables in addition to local components dealing with spatiotemporal autocorrelation as well as nonstationary effects. The local components can be used to infer unknown spatiotemporal processes; the global component describes how the species is influenced by the environment and can be used for predictions, allowing the fitting of many well-known regression relationships, ranging from simple linear models to complex decision trees or from additive models to models inspired by machine learning procedures. The reliability of spatiotemporal predictions can be qualitatively predicted by separately evaluating the importance of local and global effects. We demonstrate the potential of the new approach by modeling the breeding distribution of the Red Kite (Milvus milvus), a bird of prey occurring predominantly in Western Europe, based on presence/absence data from two mapping campaigns using grids of 40 km 2 in Bavaria. The global component of the model selected seven environmental variables extracted from the CORINE and WorldClim databases to predict Red Kite breeding. The effect of altitude was found to be nonstationary in space, and in addition, the data were spatially autocorrelated, which suggests that a species distribution model that does not allow for spatially varying effects and spatial autocorrelation would have ignored important processes determining the distribution of Red Kite breeding across Bavaria. Thus, predictions from standard species distribution models that do not allow for real-world complexities may be considerably erroneous. Our analysis of Red Kite breeding exemplifies the potential of the innovative approach for species distribution models. The method is also applicable to modeling count data.
Numerous pollen records across the upper Tibetan Plateau indicate that in the early part of the mid-Holocene, Kobresia-rich high-alpine meadows invaded areas formerly dominated by alpine steppe vegetation rich in Artemisia. We examine climate, land-use, and CO2 concentration changes as potential drivers for this marked vegetation change. The climatic implications of these vegetational shifts are explored by applying a newly developed pollen-based moisture-balance transfer-function to fossil pollen spectra from Koucha Lake on the north-eastern Tibetan Plateau (34.0 degrees N; 97.2 degrees E; 4540 m a.s.l.) and Xuguo Lake on the central Tibetan Plateau (31.97 degrees N; 90.3 degrees E; 4595 m a.s.l.), both located in the meadow-steppe transition zone. Reconstructed moisture-balances were markedly reduced (by similar to 150-180 mm) during the early mid-Holocene compared to the late-Holocene. These findings contradict most other records from the Indian monsoonal realm and also most non-pollen records from the Tibetan Plateau that indicate a rather wet early- and mid-Holocene. The extent and timing of anthropogenic land-use involving grazing by large herbivores on the upper Tibetan Plateau and its possible impacts on high-alpine vegetation are still mostly unknown due to the lack of relevant archaeological evidence. Arguments against a mainly anthropogenic origin of Kobresia high-alpine meadows are the discovery of the widespread expansion of obviously 'natural' Kobresia meadows on the south-eastern Tibetan Plateau during the Lateglacial period indicating the natural origin of this vegetation type and the lack of any concurrence between modern human-driven vegetation shifts and the mid-Holocene compositional changes. Vegetation types are known to respond to atmospheric CO2 concentration changes, at least on glacial-interglacial scales. This assumption is confirmed by our sensitivity study where we model Tibetan vegetation at different CO2 concentrations of 375 (present-day), 260 (early Holocene), and 650 ppm (future scenario) using the BIOME4 global vegetation model. Previous experimental studies confirm that vegetation growing on dry and high sites is particularly sensitive to CO2 changes. Here we propose that the replacement of drought-resistant alpine steppes (that are well adapted to low CO2 concentrations) by mesic Kobresia meadows can, at least, be partly interpreted as a response to the increase of CO2 concentration since 7000 years ago due to fertilization and water-saving effects. Our hypothesis is corroborated by former CO2 fertilization experiments performed on various dry grasslands and by the strong recent expansion of high-alpine meadows documented by remote sensing studies in response to recent CO2 increases.
The Valais units in Savoy (Zone des BrSches de Tarentaise) have been re-mapped in great detail and are subject of combined stratigraphic, structural and petrological investigations summarized in this contribution. The sediments and rare relics of basement, together with Cretaceous age mafic and ultramafic rocks of the Valais palaeogeographical domain, represent the heavily deformed relics of the former distal European margin (External Valais units) and an ocean-continent transition (Internal Valais unit or Versoyen unit) that formed during rifting. This rifting led to the opening of the Valais ocean, a northern branch of the Alpine Tethys. Post-rift sediments referred to as "Valais trilogy" stratigraphically overlie both External and Internal Valais successions above an angular unconformity formed in Barremian to Aptian times, providing robust evidence for the timing of the opening of the Valais ocean. The Valais units in Savoy are part of a second and more external mid-Eocene high-pressure belt in the Alps that sutured the Brian double dagger onnais microcontinent to Europe. Top-N D1-deformation led to the formation of a nappe stack that emplaced the largely eclogite-facies Internal Valais unit (Versoyen) onto blueschist-facies External Valais units. The latter originally consisted of, from internal to external, the Petit St. Bernard unit, the Roc de l'Enfer unit, the MoA >> tiers unit and the Quermoz unit. Ongoing top-N D2-thrusting and folding substantially modified this nappe stack. Post 35 Ma D3 folding led to relatively minor modifications of the nappe stack within the Valais units but was associated with substantial top-WNW thrusting of the Valais units over the Dauphinois units along the Roselend thrust during W-directed indentation of the Adria block contributing to the formation of the arc of the Western Alps.
This study monitors regional changes in the crystallinity of carbonaceous matter (CM) by applying Micro-Raman spectroscopy to a total of 214 metasediment samples (largely so-called Bundnerschiefer) dominantly metamorphosed under blueschist- to amphibolite-facies conditions. They were collected within the northeastern margin of the Lepontine dome and easterly adjacent areas of the Swiss Central Alps. Three-dimensional mapping of isotemperature contours in map and profile views shows that the isotemperature contours associated with the Miocene Barrow-type Lepontine metamorphic event cut across refolded nappe contacts, both along and across strike within the northeastern margin of the Lepontine dome and adjacent areas. Further to the northeast, the isotemperature contours reflect temperatures reached during the Late Eocene subduction-related blueschist-facies event and/or during subsequent near-isothermal decompression; these contours appear folded by younger, large-scale post-nappe-stacking folds. A substantial jump in the recorded maximum temperatures across the tectonic contact between the frontal Adula nappe complex and surrounding metasediments indicates that this contact accommodated differential tectonic movement of the Adula nappe with respect to the enveloping Bundnerschiefer after maximum temperatures were reached within the northern Adula nappe, i.e. after Late Eocene time.
A general mean annual temperature increase accompanied with substantial glacial retreat has been noted on the Tibetan Plateau during the last two centuries but most significantly since the mid 1950s. These climate trends are particularly apparent on the southeastern Tibetan Plateau. However, the Tibetan Plateau (due to its heterogeneous mountain landscape) has very complex and spatially differing temperature and precipitations patterns. As a result, intensive palaeolimnological investigations are necessary to decipher these climatic patterns and to understand ecological responses to recent environmental change. Here we present palaeolimnological results from a (210)Pb/(137)Cs-dated sediment core spanning approximately the last 200 years from a remote high-mountain lake (LC6 Lake, working name) on the southeastern Tibetan Plateau. Sediment profiles of diatoms, organic variables (TOC, C:N) and grain size were investigated. The (210)Pb record suggests a period of rapid sedimentation, which might be linked to major tectonic events in the region ca. 1950. Furthermore, unusually high (210)Pb supply rates over the last 50 years suggest that the lake has possibly been subjected to increasing precipitation rates, sediment focussing and/or increased spring thaw. The majority of diatom taxa encountered in the core are typical of slightly acidic to circumneutral, oligotrophic, electrolyte-poor lakes. Diatom species assemblages were rich, and dominated by Cyclotella sp., Achnanthes sp., Aulacoseira sp. and fragilarioid taxa. Diatom compositional change was minimal over the 200-year period (DCCA = 0.85 SD, p = 0.59); only a slightly more diverse but unstable diatom assemblage was recorded during the past 50 years. The results indicate that large-scale environmental changes recorded in the twentieth century (i.e. increased precipitation and temperatures) are likely having an affect on the LC6 Lake, but so far these impacts are more apparent on the lake geochemistry than on the diatom flora. Local and/or regional peculiarities, such as increasing precipitation and cloud cover, or localized climatic phenomena, such as negative climate feedbacks, might have offset the effects of increasing mean surface temperatures.
This study presents the first Si isotope data of the principle Si pools in soils determined by a UV femtosecond laser ablation system coupled to a multicollector inductively coupled plasma mass spectrometer (MC-ICP-MS). This method reveals accurate and precise Si isotope data on bulk materials, and at high spatial resolution, on the mineral scale. The following Si pools have been investigated: a) the Si source to soils on all major silicate minerals on thin sections from bedrock fragments in the soil profiles; b) bulk soils (particle size <2 mm) after fusion to glass beads with an iridium-strip heater or pressed into powder pellets: c) separated clay fractions as pressed powder pellets and e) separated phytoliths as pressed powder pellets. Multiple analyses of three rock standards, BHVO-2, AGV-1 and RGM-1 as fused glass beads and as pressed powder pellets, reveal delta(30)Si values within the expected range of igneous rocks. The MPI-DING reference glass KL2-G exhibits the same Si isotope composition after remelting by an iridium-strip heater showing that this technique does not alter the isotope composition of the glass.
We used this approach to investigated two immature Cambisols developed on sandstone and paragneiss in the Black Forest (Germany), respectively. Bulk soils show a largely uniform Si isotope signature for different horizons and locations, which is close to those of primary quartz and feldspar with delta(30)Si values around -0.4 parts per thousand. Soil clay formation is associated with limited Si mobility, which preserves initial Si isotope signatures of parental minerals. An exception is the organic horizon of the paragneiss catchment where intense weathering leads to a high mobility of Si and significant negative isotope signatures as low as to -1.00 parts per thousand in bulk soils. Biogenic opal in the form of phytoliths, exhibits negative Si isotope signatures of about -0.4 parts per thousand. These results demonstrate that UV femtosecond laser ablation MC-ICP-MS provides a tool to characterize the Si isotope signature of the principle Si pools left behind after weathering and Si transport have altered soils. These results can now serve as a fingerprint of the residual solids that can be used to explain the isotope composition of dissolved Si in soil solutions and river water, which is mostly enriched in the heavy isotopes.
We present the first winter (December to March) rainfall reconstruction based in a novel proxy, the thickness of annual calcite laminations preserved in Lake La Cruz (central-eastern Spain). A previous calibration analysis between laminae thickness and the instrumental data series (1950 to present) indicated a highly significant correlation with winter rainfall. Therefore this study attempts the winter rainfall reconstruction since the onset of laminations (1579 a.d.) by means of the calibration function previously developed. The verification analysis between inferred annual values and earlier instrumental data (1859-1949) confirms the suitability of this novel proxy and the reliability of the series reconstructed. The reconstructed series show the fluctuating character of winter rainfall in the western Mediterranean area; interdecadal dry periods alternated with wetter periods following, in a board sense, the pattern recorded by documentary sources in other regions of the Iberian Peninsula. At present times regional winter rainfall anomalies are highly correlated with the phase of the North Atlantic Oscillation (NAO). However the time series analysis showed the dominance of nonstationary components at high frequencies of the climate signal over the last four centuries suggesting that the connection between winter rainfall and the NAO has not been stable over time and also other modes of variability, not only NAO, may have conditioned winter rainfall variability.
Recent studies have shown that rhizosphere hydraulic properties may differ from those of the bulk soil. Specifically, mucilage at the root-soil interface may increase the rhizosphere water holding capacity and hydraulic conductivity during drying. The goal of this study was to point out the implications of such altered rhizosphere hydraulic properties for soil-plant water relations. We addressed this problem through modeling based on a steady-rate approach. We calculated the water flow toward a single root assuming that the rhizosphere and bulk soil were two concentric cylinders having different hydraulic properties. Based on our previous experimental results, we assumed that the rhizosphere had higher water holding capacity and unsaturated conductivity than the bulk soil. The results showed that the water potential gradients in the rhizosphere were much smaller than in the bulk soil. The consequence is that the rhizosphere attenuated and delayed the drop in water potential in the vicinity of the root surface when the soil dried. This led to increased water availability to plants, as well as to higher effective conductivity under unsaturated conditions. The reasons were two: (i) thanks to the high unsaturated conductivity of the rhizosphere, the radius of water uptake was extended from the root to the rhizosphere surface; and (ii) thanks to the high soil water capacity of the rhizosphere, the water depletion in the bulk soil was compensated by water depletion in the rhizosphere. We conclude that under the assumed conditions, the rhizosphere works as an optimal hydraulic conductor and as a reservoir of water that can be taken up when water in the bulk soil becomes limiting.
Love wave contribution to the ambient vibration H/V amplitude peak observed with array measurements
(2011)
This study applies array methods to measure the relative proportions of Love and Rayleigh waves in the ambient vibration wavefield. Information on these properties is of special relevance for frequencies around the horizontal-to-vertical (H/V) spectral amplitude ratio peak. The analysis of H/V curves, a popular technique in site characterisation, commonly assumes that the curves represent the frequency-dependent Rayleigh wave ellipticity. For the detailed interpretation of amplitudes or the inversion of the curves, it is therefore necessary to estimate and correct for the contribution of other wave types to the ambient vibration wavefield. I use available ambient vibration array measurements to determine the relative amount of Love and Rayleigh waves on the horizontal components by frequency-dependent analysis of the main propagation and polarisation directions, with a special emphasis on the H/V peak frequency as determined from the same recordings. Tests with synthetic data demonstrate the feasibility of this approach, at least in the presence of dominant source regions. Analysis of the data from 12 measurements at nine European sites, which include shallow as well as deep locations that span a wide range of impedance contrasts at the sediment-bedrock interface, indicates that the relative contribution of Rayleigh waves varies widely with frequency, from close to 0% to more than 70%. While most data sets show relative Rayleigh wave contributions between 40% and 50% around the H/V peak, there are also examples where Love waves clearly dominate the wavefield at the H/V peak, even for a site with a low impedance contrast. Longer-term measurements at one site indicate temporal variations in the relative Rayleigh wave content between day- and nighttime. Results calculated with the method introduced herein generally compare well with results of modified spatial autocorrelation analysis. These two methods might be used in a complimentary fashion, as both rely on different properties of the ambient vibration wavefield. This study illustrates that it is possible to measure the relative Rayleigh wave content of the noise wavefield from array data. Furthermore, the examples presented herein indicate it is important to estimate this property, as the assumption that there are an equal proportion of Love and Rayleigh waves is not always correct.
Salt diapirs are common features of sedimentary basins. If close to the surface, they can bear a significant hazard due to possible dissolution sinkholes, karst formation and collapse dolines or their influence on ground water chemistry. We investigate the potential of ambient vibration techniques to map the 3-D roof morphology of shallow salt diapirs. Horizontal-to-vertical (H/V) spectral peaks are derived at more than 900 positions above a shallow diapir beneath the city area of Hamburg, Germany, and are used to infer the depth of the first strong impedance contrast. In addition, 15 small-scale array measurements are conducted at different positions in order to compute frequency-dependent phase velocities of Rayleigh waves between 0.5 and 25 Hz. The dispersion curves are inverted together with the H/V peak frequency to obtain shear-wave velocity profiles. Additionally, we compare the morphology derived from H/V and array measurements to borehole lithology and a gravity-based 3-D model of the salt diapir. Both methods give consistent results in agreement with major features indicated by the independent data. An important result is that H/V and array measurements are better suited to identify weathered gypsum caprocks or gypsum floaters, while gravity-derived models better sample the interface between sediments and homogeneous salt. We further investigate qualitatively the influence of the 3-D subsurface topography of the salt diapir on the validity of local 1-D inversion results from ambient vibration dispersion curve inversion.
This paper presents the results of a research aimed to quantify suspended sediment transport in three experimental catchments in southern Chile, to compare measured suspended sediment load with estimated erosion using the Universal Soil Loss Equation (USLE) applied in a GIS environment and to validate de Modified Universal Soil Loss Equation (MUSLE) used to estimate suspended sediment loads from forest catchments. The catchments are Los Pinos (94.2 ha), Los Ulmos 1 (12.6 ha) and Los Ulmos 2 (17.7 ha). Soil losses estimated with USLE for the three catchments are higher than those measured in runoff experimental lots under bare soil conditions, which could indicate an overestimation of the LS calculated in GIS and the fact that the USLE model does not compute sediment deposit and storage within the catchment. A statistical significant relation was found between measured and estimated (MUSLE) suspended sediment load, which would indicate that this model could be applied to estimate suspended sediment load from small catchments in southern Chile.
The public promotion of renewable energies is expected to increase the number of biogas plants and stimulate energy crops cultivation (e. g. maize) in Germany. In order to assess the indirect effects of the resulting land-use changes on biodiversity, we developed six land-use scenarios and simulated the responses of six farmland wildlife species with the spatially explicit agent-based model system ALMaSS. The scenarios differed in composition and spatial configuration of arable crops. We implemented scenarios where maize for energy production replaced 15% and 30% of the area covered by other cash crops. Biogas maize farms were either randomly distributed or located within small or large aggregation clusters. The animal species investigated were skylark (Alauda arvensis), grey partridge (Perdix perdix), European brown hare (Lepus europaeus), field vole (Microtus agrestis), a linyphiid spider (Erigone atra) and a carabid beetle (Bembidion lampros). The changes in crop composition had a negative effect on the population sizes of skylark, partridge and hare and a positive effect on the population sizes of spider and beetle and no effect on the population size of vole. An aggregated cultivation of maize amplified these effects for skylark. Species responses to changes in the crop composition were consistent across three differently structured landscapes. Our work suggests that with the compliance to some recommendations, negative effects of biogas-related land-use change on the populations of the six representative farmland species can largely be avoided.
The Upper Cretaceous La Cova limestones (southern Pyrenees, Spain) host a rich and diverse larger foraminiferal fauna, which represents the first diversification of K-strategists after the mass extinction at the Cenomanian-Turonian boundary.
The stratigraphic distribution of the main taxa of larger foraminifera defines two assemblages. The first assemblage is characterised by the first appearance of lacazinids (Pseudolacazina loeblichi) and mean-dropsinids (Eofallotia simplex), by the large agglutinated Montsechiana montsechiensis, and by several species of complex rotalids (Rotorbinella campaniola, Iberorotalia reicheli, Orbitokhatina wondersmitti and Calcarinella schaubi). The second assemblage is defined by the appearance of Lacazina pyrenaica, Palandrosina taxyae and Martiguesia cyclamminiformis.
A late Coniacian-early Santonian age was so far accepted for the La Cova limestones, based on indirect correlation with deep-water fades bearing planktic foraminifers of the Dicarinella concavata zone. Strontium isotope stratigraphy, based on many samples of pristine biotic calcite of rudists and ostreids, indicates that the La Cova limestones span from the early Coniacian to the early-middle Santonian boundary. The first assemblage of larger foraminifera appears very close to the early-middle Coniacian boundary and reaches its full diversity by the middle Coniacian. The originations defining the second assemblage are dated as earliest Santonian: they represent important bioevents to define the Coniacian-Santonian boundary in the shallow-water facies of the South Pyrenean province.
By means of the calibration of strontium isotope stratigraphy to the Geological Time Scale, the larger foraminiferal assemblages of the La Cova limestones can be correlated to the standard biozonal scheme of ammonites, planktonic foraminifers and calcareous nannoplankton. This correlation is a first step toward a larger foraminifera standard biozonation for Upper Cretaceous carbonate platform facies.
This study investigates spatial patterns and temporal dynamics of aquifer-river exchange flow at a reach of the River Leith, UK. Observations of sub-channel vertical hydraulic gradients at the field site indicate the dominance of groundwater up-welling into the river and the absence of groundwater recharge from surface water. However, observed hydraulic heads do not provide information on potential surface water infiltration into the top 0-15 cm of the streambed as these depths are not covered by the existing experimental infrastructure. In order to evaluate whether surface water infiltration is likely to occur outside the 'window of detection', i.e. the shallow streambed, a numerical groundwater model is used to simulate hydrological exchanges between the aquifer and the river. Transient simulations of the successfully validated model (Nash and Sutcliff efficiency of 0.91) suggest that surface water infiltration is marginal and that the possibility of significant volumes of surface water infiltrating into non-monitored shallow streambed sediments can be excluded for the simulation period. Furthermore, the simulation results show that with increasing head differences between river and aquifer towards the end of the simulation period, the impact of streambed topography and hydraulic conductivity on spatial patterns of exchange flow rates decreases. A set of peak flow scenarios with altered groundwater-surface water head gradients is simulated in order to quantify the potential for surface water infiltration during characteristic winter flow conditions following the observation period. The results indicate that, particularly at the beginning of peak flow conditions, head gradients are likely to cause substantial increase in surface water infiltration into the streambed. The study highlights the potential for the improvement of process understanding of hyporheic exchange flow patterns at the stream reach scale by simulating aquifer-river exchange fluxes with a standard numerical groundwater model and a simple but robust model structure and parameterization. Copyright
Magnitude estimation for microseismicity induced during the KTB 2004/2005 injection experiment
(2011)
We determined the magnitudes of 2540 microseismic events measured at one single 3C borehole geophone at the German Deep Drilling Site (known by the German acronym, KTB) during the injection phase 2004/2005. For this task we developed a three-step approach. First, we estimated local magnitudes of 104 larger events with a standard method based on amplitude measurements at near-surface stations. Second, we investigated a series of parameters to characterize the size of these events using the seismograms of the borehole sensor, and we compared them statistically with the local magnitudes. Third, we extrapolated the regression curve to obtain the magnitudes of 2436 events that were only measured at the borehole geophone. This method improved the magnitude of completeness for the KTB data set by more than one order down to M = -2.75. The resulting b-value for all events was 0.78, which is similar to the b-value obtained from taking only the greater events with standard local magnitude estimation from near-surface stations, b = 0.86. The more complete magnitude catalog was required to study the magnitude distribution with time and to characterize the seismotectonic state of the KTB injection site. The event distribution with time was consistent with prediction from theory assuming pore pressure diffusion as the underlying mechanism to trigger the events. The value we obtained for the seismogenic index of -4 suggested that the seismic hazard potential at the KTB site is comparatively low.
Enhancing the resolution and accuracy of surface ground-penetrating radar (GPR) reflection data by inverse filtering to recover a zero-phased band-limited reflectivity image requires a deconvolution technique that takes the mixed-phase character of the embedded wavelet into account. In contrast, standard stochastic deconvolution techniques assume that the wavelet is minimum phase and, hence, often meet with limited success when applied to GPR data. We present a new general-purpose blind deconvolution algorithm for mixed-phase wavelet estimation and deconvolution that (1) uses the parametrization of a mixed-phase wavelet as the convolution of the wavelet's minimum-phase equivalent with a dispersive all-pass filter, (2) includes prior information about the wavelet to be estimated in a Bayesian framework, and (3) relies on the assumption of a sparse reflectivity. Solving the normal equations using the data autocorrelation function provides an inverse filter that optimally removes the minimum-phase equivalent of the wavelet from the data, which leaves traces with a balanced amplitude spectrum but distorted phase. To compensate for the remaining phase errors, we invert in the frequency domain for an all-pass filter thereby taking advantage of the fact that the action of the all-pass filter is exclusively contained in its phase spectrum. A key element of our algorithm and a novelty in blind deconvolution is the inclusion of prior information that allows resolving ambiguities in polarity and timing that cannot be resolved using the sparseness measure alone. We employ a global inversion approach for non-linear optimization to find the all-pass filter phase values for each signal frequency. We tested the robustness and reliability of our algorithm on synthetic data with different wavelets, 1-D reflectivity models of different complexity, varying levels of added noise, and different types of prior information. When applied to realistic synthetic 2-D data and 2-D field data, we obtain images with increased temporal resolution compared to the results of standard processing.
In this study, Nuclear Magnetic Resonance (NMR), a non-destructive measurement technique, has been applied for investigation of iron turn-over processes. In non-invasive laboratory experiments, iron dissolution and precipitation reactions in saturated natural sands were observed spatially and temporally. These processes play an important role in groundwater with varying redox and pH conditions. Redox reactions turning Fe2+ into Fe3+ and Fe3+ into Fe2+ were detected in aqueous solution by the difference in magnetic relaxation times. Furthermore, the spatial distribution of the iron reduction reaction, the consumption and diffusive transfer to and from the reaction sites, was observed in a 1D set-up with natural sands. The achieved spatial resolution was less than one millimetre while repeating measurements every half an hour. It showed the system changing from diffusion-limited to reaction-limited.
In the Aconcagua basin the modernization and intensification of land use have resulted both in an agricultural conversion, and the increase in area cultivated. This process has been very intense as far that have entered new areas beyond the limits of the valley and towards sloping ground. This signals changes in land use over the vegetation coverage models. The analysis of the geomorphology of the coastline in detail, accounts for the occurrence of transformations on the beach and dunes that are spatially less significant, but evolutionarily important, since they show a sediment feeding, which could be correlated with changes in vegetation cover in the basin
While several empirical and theoretical studies have clearly shown the negative effects of climate or landscape changes on population and species survival only few of them addressed combined and correlated consequences of these key environmental drivers. This also includes positive landscape changes such as active habitat management and restoration to buffer the negative effects of deteriorating climatic conditions. In this study, we apply a conceptual spatial modelling approach based on functional types to explore the effects of both positive and negative correlations between changes in habitat and climate conditions on the survival of spatially structured populations. We test the effect of different climate and landscape change scenarios on four different functional types that represent a broad spectrum of species characterised by their landscape level carrying capacity, the local population turnover rates at the patch level (K-strategies vs. r-strategies) and dispersal characterstics. As expected, simulation results show that correlated landscape and climatic changes can accelerate (in case of habitat loss or degradation) or slow down (in case of habitat gain or improvement) regional species extinction. However, the strength of the combined changes depends on local turnover at the patch level, the overall landscape capacity of the species, and its specific dispersal characteristics. Under all scenarios of correlated changes in habitat and climate conditions we found the highest sensitivity for functional types representing species with a low landscape capacity but a high population growth rate and a strong density regulation causing a high turnover at the local patch level.
The relative importance of habitat loss or habitat degradation, in combination with climate deterioration, differed among the functional types. However, an increase in regional capacity revealed a similar response pattern: For all types, habitat improvement led to higher survival times than habitat gain, i.e. the establishment of new habitat patches. This suggests that improving local habitat quality at a regional scale is a more promising conservation strategy under climate change than implementing new habitat patches. This conceptual modelling study provides a general framework to better understand and support the management of populations prone to complex environmental changes.
Mountain ecosystems are commonly regarded as being highly sensitive to global change. Due to the system complexity and multifaceted interacting drivers, however, understanding current responses and predicting future changes in these ecosystems is extremely difficult. We aim to discuss potential effects of global change on mountain ecosystems and give examples of the underlying response mechanisms as they are understood at present. Based on the development of scientific global change research in mountains and its recent structures, we identify future research needs, highlighting the major lack and the importance of integrated studies that implement multi-factor, multi-method, multi-scale, and interdisciplinary research.
Recurrence plots and recurrence quantification analysis have become popular in the last two decades. Recurrence based methods have on the one hand a deep foundation in the theory of dynamical systems and are on the other hand powerful tools for the investigation of a variety of problems. The increasing interest encompasses the growing risk of misuse and uncritical application of these methods. Therefore, we point out potential problems and pitfalls related to different aspects of the application of recurrence plots and recurrence quantification analysis.
Humidity is an important determinant of the mycotoxin production (DON, ZEA) by Fusarium species in the grain ears. From a landscape perspective humidity is not evenly distributed across fields. The topographically-controlled redistribution of water within a single field rather leads to spatially heterogeneous soil water content and air humidity. Therefore we hypothesized that the spatial distribution of mycotoxins is related to these topographically-controlled factors. To test this hypothesis we studied the mycotoxin concentrations at contrasting topographic relief positions, i.e. hilltops and depressions characterized by soils of different soil moisture regimes, on ten winter wheat fields in 2006 and 2007. Maize was the preceding crop and minimum tillage was practiced in the fields. The different topographic positions were associated with moderate differences in DON and ZEA concentrations in 2006, but with significant differences in 2007, with six times higher median ZEA and two times higher median DON detected at depression sites compared to the hilltops. The depression sites correspond to a higher topographic wetness index as well as redoximorphic properties in soil profiles, which empirically supports our hypothesis at least for years showing wetter conditions in sensitive time windows for Fusarium infections.
Large-scale soy agriculture in the southern Brazilian Amazon now rivals deforestation for pasture as the region's predominant form of land use change. Such landscape-level change can have substantial consequences for local and regional hydrology, but these effects remain relatively unstudied in this ecologically and economically important region. We examined how the conversion to soy agriculture influences water balances and stormflows using stream discharge (water yields) and the timing of discharge (stream hydrographs) in small (2.5-13.5 km2) forested and soy headwater watersheds in the Upper Xingu Watershed in the state of Mato Grosso, Brazil. We monitored water yield for 1 year in three forested and four soy watersheds. Mean daily water yields were approximately four times higher in soy than forested watersheds, and soy watersheds showed greater seasonal variability in discharge. The contribution of stormflows to annual streamflow in all streams was low (< 13% of annual streamflow), and the contribution of stormflow to streamflow did not differ between land uses. If the increases in water yield observed in this study are typical, landscape-scale conversion to soy substantially alters water-balance, potentially altering the regional hydrology over large areas of the southern Amazon.
Shell architectures of the larger foraminiferal genera cyclopseudedomia, and Rhapydionina were studied by comparing topotypes of previously described species with new specimens retrieved from Late Cretaceous shallow-water carbonates of Pylos (Peloponnese, Greece), where the three genera are found in association. The megalospheric generation of each genus exhibits a distinctive shell shape in adult specimens (i.e., fan-shaped in Cyclopseudedomia, conical in Rhapydionina, and cylindrical in Cuvillierinella). Although their microspheric adults are similarly thin, flat, and discoidal, they can be identified at the genus level by means of a detailed structural analysis. Cavillierinella shows the septula to be interrupted by a large preseptal space, while Cyclopseudedomia and Rhapydionina exhibit continuous, non-interrupted septula. In addition, Cyclopseudedomia presents only one row of medullar chamberlets, whereas Rhapydionina shows numerous medullar chamberlets distributed in a thick basal layer.
Two new species, Cuvillierinella pylosensis and Rhapydionina fleuryi, are described. The former is a more complex taxon than the type species, C salentina, while the latter corresponds to a more primitive species, R. liburnica. Strontium-isotope stratigraphy indicates an uppermost Campanian-lowermost Maastrichtian age for these new species.
Geophysical datasets sensitive to different physical parameters can be used to improve resolution of Earth's internal structure. Herein, we jointly invert long-period magnetotelluric (MT) data and surface-wave dispersion curves. Our approach is based on a joint inversion using a genetic algorithm for a one-dimensional (1-D) isotropic structure, which we extend to 1-D anisotropic media. We apply our new anisotropic joint inversion to datasets from Central Germany demonstrating the capacity of our joint inversion algorithm to establish a 1-D anisotropic model that fits MT and seismic datasets simultaneously and providing new information regarding the deep structure in Central Germany. The lithosphere/asthenosphere boundary is found at approx. 84 km depth and two main anisotropic layers with coincident most conductive/seismic fast-axis direction are resolved at lower crustal and asthenospheric depths. We also quantify the amount of seismic and electrical anisotropy in the asthenosphere showing an emerging agreement between the two anisotropic coefficients.
End member mixing analysis (EMMA) is a commonly applied method to identify and quantify the dominant runoff producing sources of water. It employs tracers to determine the dimensionality of the hydrologic system. Many EMMA studies have been conducted using two to six tracers, with some of the main tracers being Ca, Na, Cl(-), water isotopes, and alkalinity. Few studies use larger tracer sets including minor trace elements such as Li, Rb, Sr, and Ba. None of the studies has addressed the question of the tracer set size and composition, despite the fact that these determine which and how many end members (EM) will be identified. We examine how tracer set size and composition affects the conceptual model that results from an EMMA. We developed an automatic procedure that conducts EMMA while iteratively changing tracer set size and composition. We used a set of 14 tracers and 9 EMs. The validity of the resulting conceptual models was investigated under the aspects of dimensionality, EM combinations, and contributions to stream water. From the 16,369 possibilities, 23 delivered plausible results. The resulting conceptual models are highly sensitive to the tracer set size and composition. The moderate reproducibility of EM contributions indicates a still missing EM. It also emphasizes that the major elements are not always the most useful tracers and that larger tracer sets have an enhanced capacity to avoid false conclusions about catchment functioning. The presented approach produces results that may not be apparent from the traditional approach and it is a first step to add the idea of statistical significance to the EMMA approach.
Dust has the potential to modify global climate by influencing the radiative balance of the atmosphere and by supplying iron and other essential limiting micronutrients to the ocean(1,2). Indeed, dust supply to the Southern Ocean increases during ice ages, and 'iron fertilization' of the subantarctic zone may have contributed up to 40 parts per million by volume (p. p. m. v.) of the decrease (80-100 p. p. m. v.) in atmospheric carbon dioxide observed during late Pleistocene glacial cycles(3-7). So far, however, the magnitude of Southern Ocean dust deposition in earlier times and its role in the development and evolution of Pleistocene glacial cycles have remained unclear. Here we report a high-resolution record of dust and iron supply to the Southern Ocean over the past four million years, derived from the analysis of marine sediments from ODP Site 1090, located in the Atlantic sector of the subantarctic zone. The close correspondence of our dust and iron deposition records with Antarctic ice core reconstructions of dust flux covering the past 800,000 years (refs 8, 9) indicates that both of these archives record large-scale deposition changes that should apply to most of the Southern Ocean, validating previous interpretations of the ice core data. The extension of the record beyond the interval covered by the Antarctic ice cores reveals that, in contrast to the relatively gradual intensification of glacial cycles over the past three million years, Southern Ocean dust and iron flux rose sharply at the Mid-Pleistocene climatic transition around 1.25 million years ago. This finding complements previous observations over late Pleistocene glacial cycles(5,8,9), providing new evidence of a tight connection between high dust input to the Southern Ocean and the emergence of the deep glaciations that characterize the past one million years of Earth history.
Implementation of cartographic symbols for planetary mapping in geographic information systems
(2011)
The steadily growing international interest in the exploration of planets in our Solar System and many advances in the development of space-sensor technology have led to the launch of a multitude of planetary missions to Mercury, Venus, the Earth's moon, Mars and various Outer-Solar System objects, such as the Jovian and Saturnian satellites. Camera instruments carried along on these missions image surfaces in different wavelength ranges and under different viewing angles, permitting additional data to be derived, such as spectral data or digital terrain models. Such data enable researchers to explore and investigate the development of planetary surfaces by analyzing and interpreting the inventory of surface units and structures. Results of such work are commonly abstracted and represented in thematic, mostly geological and geomorphological, maps. In order to facilitate efficient collaboration among different planetary research disciplines, mapping results need to be prepared, described, managed, archived, and visualized in a uniform way. These tasks have been increasingly carried out by means of computer-based geographic information systems (GIS or Cl systems) which have come to be widely employed in the field of planetary research since the last two decades. In this paper we focus on the simplification of mapping processes, putting specific emphasis on a cartographically correct visualization of planetary mapping data using GIS-based environments. We present and discuss the implementation of a set of standardized cartographic symbols for planetary mapping based on the Digital Cartographic Standard for Geologic Map Symbolization as prepared by the United States Geological Survey (USGS) for the Federal Geographic Data Committee (FGDC). Furthermore, we discuss various options to integrate this symbol catalog into generic GI systems, and more specifically into the Environmental Systems Research Institute's (ESRI) ArcGIS environment, and focus on requirements for symbol definitions in the field of planetary mapping. A symbology of this type can be embedded into any modular GIS environment capable in dealing with external stand-alone as well as database-driven management of symbol sets. Using such a uniform GIS-based symbol catalog will give the research community access to map results already cartographically elaborated, enabling them to create digital maps as a secondary data source in subsequent studies.
The chemical composition of eclogites, found as boulders in a Tertiary conglomerate from the Guajira Peninsula, Colombia suggests that these rocks are mainly metamorphosed basaltic andesites. They are depleted in LILE elements compared to MORB, have a negative Nb-anomaly and flat to enriched REE patterns, suggesting that their protoliths evolved in a subduction related tectonic setting. They show island-arc affinities and are similar to primitive island-arc rocks described in the Caribbean. The geochemical characteristics are comparable to low-grade greenschists from the nearby Etpana Terrane, which are interpreted as part of a Cretaceous intra-oceanic arc. These data support evidence that the eclogites and the Etpana terrane rocks formed from the same volcano-sedimentary sequence. Part of this sequence was accreted onto the margin and another was incorporated into the subduction channel and metamorphosed at eclogite facies conditions. Ar-40-Ar-39 ages of 79.2 +/- 1.1Ma and 82.2 +/- 2.5Ma determined on white micas, separated from two eclogite samples, are interpreted to be related to the cooling of the main metamorphic event. The formation of a common volcano-sedimentary protolith and subsequent metamorphism of these units record the ongoing Late Cretaceous continental subduction of the South American margin within the Caribbean intra-oceanic arc subduction zone. This gave way to an arc-continent collision between the Caribbean and the South American plates, where this sequence was exhumed after the Campanian.
Interannual rainfall variations in equatorial East Africa are tightly linked to the El Nino Southern Oscillation (ENSO), with more rain and flooding during El Nino and droughts in La Nina years, both having severe impacts on human habitation and food security. Here we report evidence from an annually laminated lake sediment record from southeastern Kenya for interannual to centennial-scale changes in ENSO-related rainfall variability during the last three millennia and for reductions in both the mean rate and the variability of rainfall in East Africa during the Last Glacial period. Climate model simulations support forward extrapolation from these lake sediment data that future warming will intensify the interannual variability of East Africa's rainfall.
In this study we investigated the effect of different land use options (wildlife versus livestock) on species richness of plants and reptiles in the protected Kgalagadi Transfrontier Park (KTP) versus adjacent non-protected farmland within the same savanna habitat type (Aoub dune veld) in the southern Kalahari, South Africa. Our results show that both plant and reptile species richness as well as plant cover and reptile abundance was significantly higher in the protected KTP than in the non-protected farmland. The higher proportion of shrub but lower proportions of perennial grass cover, herb cover, and herb species richness in the farmland can be explained by higher stocking rates and the differences in feeding behaviour between native wild ungulates (e.g. Antidorcas marsupialis, Oryx gazella) and livestock (mainly sheep). The reptile's prey availability and microhabitats (perennial grass tussocks and rodent burrows) for thermoregulation and protection against predators were significantly lower in the farmland. To conclude, our results clearly show that long term effects of different land use options (wildlife in protected KTP versus extensive livestock production in the non-protected farmland) even within the same habitat type have led to significant changes in vegetation composition, availability of microhabitat structures and in the reptile community.
In present-day land and marine controlled-source electromagnetic (CSEM) surveys, electromagnetic fields are commonly generated using wires that are hundreds of metres long. Nevertheless, simulations of CSEM data often approximate these sources as point dipoles. Although this is justified for sufficiently large source-receiver distances, many real surveys include frequencies and distances at which the dipole approximation is inaccurate. For 1D layered media, electromagnetic (EM) fields for point dipole sources can be computed using well-known quasi-analytical solutions and fields for sources of finite length can be synthesized by superposing point dipole fields. However, the calculation of numerous point dipole fields is computationally expensive, requiring a large number of numerical integral evaluations. We combine a more efficient representation of finite-length sources in terms of components related to the wire and its end points with very general expressions for EM fields in 1D layered media. We thus obtain a formulation that requires fewer numerical integrations than the superposition of dipole fields, permits source and receiver placement at any depth within the layer stack and can also easily be integrated into 3D modelling algorithms. Complex source geometries, such as wires bent due to surface obstructions, can be simulated by segmenting the wire and computing the responses for each segment separately. We first describe our finite-length wire expressions and then present 1D and 3D examples of EM fields due to finite-length sources for typical land and marine survey geometries and discuss differences to point dipole fields.
Three-dimensional hydrostratigraphic models from ground-penetrating radar and direct-push data
(2011)
Three-dimensional models of hydraulic conductivity and porosity are essential to understand and simulate groundwater flow in heterogeneous geological environments. However, considering the inherent limitations of traditional hydrogeological field methods in terms of resolution, alternative field approaches are needed to establish such 3-D models with sufficient accuracy. In this study, we developed a workflow combining 3-D structural information extracted from ground penetrating radar (GPR) images with 1-D in situ physical-property estimates from direct-push (DP) logging to construct a 3-D hydrostratigraphic model. To illustrate this workflow, we collected an similar to 70 m x 90 m 100 MHz 3-D GPR data set over a shallow sedimentary aquifer system resolving six different GPR facies down to similar to 15 m depth. DP logs of the relative dielectric permittivity, the relative hydraulic conductivity, the cone resistance, the sleeve friction and the pore pressure provided crucial data (1) to establish a GPR velocity model for 3-D depth migration and to check the time-to-depth conversion of the GPR data, and (2) to construct a 3-D hydrostratigraphic model. This model was built by assigning porosity values, which were computed from the DP relative dielectric permittivity logs, and DP relative hydraulic conductivity estimates to the identified GPR facies. We conclude that the integration of 3-D GPR structural images and 1-D DP logs of target physical parameters provides an efficient way for detailed 3-D subsurface characterization as needed, for example, for groundwater flow simulations.
Bayesian networks are a powerful and increasingly popular tool for reasoning under uncertainty, offering intuitive insight into (probabilistic) data-generating processes. They have been successfully applied to many different fields, including bioinformatics. In this paper, Bayesian networks are used to model the joint-probability distribution of selected earthquake, site, and ground-motion parameters. This provides a probabilistic representation of the independencies and dependencies between these variables. In particular, contrary to classical regression, Bayesian networks do not distinguish between target and predictors, treating each variable as random variable. The capability of Bayesian networks to model the ground-motion domain in probabilistic seismic hazard analysis is shown for a generic situation. A Bayesian network is learned based on a subset of the Next Generation Attenuation (NGA) dataset, using 3342 records from 154 earthquakes. Because no prior assumptions about dependencies between particular parameters are made, the learned network displays the most probable model given the data. The learned network shows that the ground-motion parameter (horizontal peak ground acceleration, PGA) is directly connected only to the moment magnitude, Joyner-Boore distance, fault mechanism, source-to-site azimuth, and depth to a shear-wave horizon of 2: 5 km/s (Z2.5). In particular, the effect of V-S30 is mediated by Z2.5. Comparisons of the PGA distributions based on the Bayesian networks with the NGA model of Boore and Atkinson (2008) show a reasonable agreement in ranges of good data coverage.
Complex layered deformation within the Aegean crust and mantle revealed by seismic anisotropy
(2011)
Continental lithosphere can undergo pervasive internal deformation, often distributed over broad zones near plate boundaries. However, because of the paucity of observational constraints on three-dimensional movement at depth, patterns of flow within the lithosphere remain uncertain. Endmember models for lithospheric flow invoke deformation localized on faults or deep shear zones or, alternatively, diffuse, viscous-fluid-like flow. Here we determine seismic Rayleigh-wave anisotropy in the crust and mantle of the Aegean region, an archetypal example of continental deformation. Our data reveal a complex, depth-dependent flow pattern within the extending lithosphere. Beneath the northern Aegean Sea, fast shear wave propagation is in a North-South direction within the mantle lithosphere, parallel to the extensional component of the current strain rate field. In the south-central Aegean, where deformation is weak at present, anisotropic fabric in the lower crust runs parallel to the direction of palaeo-extension in the Miocene. The close match of orientations of regional-scale anisotropic fabric and the directions of extension during the last significant episodes of deformation implies that at least a large part of the extension in the Aegean has been taken up by distributed viscous flow in the lower crust and lithospheric mantle.
The Takab complex is composed of a variety of metamorphic rocks including amphibolites, metapelites, mafic granulites, migmatites and meta-ultramafics, which are intruded by the granitoid. The granitoid magmatic activity occurred in relation to the subduction of the Neo-Tethys oceanic crust beneath the Iranian crust during Tertiary times. The granitoids are mainly granodiorite, quartz monzodiorite, monzonite and quartz diorite. Chemically, the magmatic rocks are characterized by ASI < 1.04, AI < 0.87 and high contents of CaO (up to similar to 14.5 wt %), which are consistent with the I-type magmatic series. Low FeO(t)/(FeO(t)+MgO) values (< 0.75) as well as low Nb, Y and K(2)O contents of the investigated rocks resemble the calc-alkaline series. Low SiO(2), K(2)O/Na(2)O and Al(2)O(3) accompanied by high CaO and FeO contents indicate melting of metabasites as an appropriate source for the intrusions. Negative Ti and Nb anomalies verify a metaluminous crustal origin for the protoliths of the investigated igneous rocks. These are comparable with compositions of the associated mafic migmatites, in the Takab metamorphic complex, which originated from the partial melting of amphibolites. Therefore, crustal melting and a collision-related origin for the Takab calc-alkaline intrusions are proposed here on the basis of mineralogy and geochemical characteristics. The P-T evolution during magmatic crystallization and subsolidus cooling stages is determined by the study of mineral chemistry of the granodiorite and the quartz diorite. Magmatic crystallization pressure and temperature for the quartz-diorite and the granodiorite are estimated to be P similar to 7.8 +/- 2.5 kbar, T similar to 760 +/- 75 degrees C and P similar to 5 +/- 1 kbar, T similar to 700 degrees C, respectively. Subsolidus conditions are consistent with temperatures of similar to 620 degrees C and similar to 600 degrees C, and pressures of similar to 5 kbar and similar to 3.5 kbar for the quartz-diorite and the granodiorite, respectively.
Investigation of a right-laterally offset channel at the Miller's Field paleoseismic site yields a late Holocene slip rate of 26.2 +6.4/-4.3 mm/yr (1 sigma) for the main trace of the San Andreas fault at Parkfield, California. This is the first well-documented geologic slip rate between the Carrizo and creeping sections of the San Andreas fault. This rate is lower than Holocene measurements along the Carrizo Plain and rates implied by far-field geodetic measurements (similar to 35 mm/yr). However, the rate is consistent with historical slip rates, measured to the northwest, along the creeping section of the San Andreas fault (<30 mm/yr). The paleoseismic exposures at the Miller's Field site reveal a pervasive fabric of clay shear bands, oriented clockwise oblique to the San Andreas fault strike and extending into the uppermost stratigraphy. This fabric is consistent with dextral aseismic creep and observations of surface slip from the 28 September 2004 M6 Parkfield earthquake. Together, this slip rate and deformation fabric suggest that the historically observed San Andreas fault slip behavior along the Parkfield section has persisted for at least a millennium, and that significant slip is accommodated by structures in a zone beyond the main San Andreas fault trace.
One of the key challenges in the context of local site effect studies is the determination of frequencies where the shakeability of the ground is enhanced. In this context, the H/V technique has become increasingly popular and peak frequencies of H/V spectral ratio are sometimes interpreted as resonance frequencies of the transmission response. In the present study, assuming that Rayleigh surface wave is dominant in H/V spectral ratio, we analyse theoretically under which conditions this may be justified and when not. We focus on 'layer over half-space' models which, although seemingly simple, capture many aspects of local site effects in real sedimentary structures. Our starting point is the ellipticity of Rayleigh waves. We use the exact formula of the H/V-ratio presented by Malischewsky & Scherbaum (2004) to investigate the main characteristics of peak and trough frequencies. We present a simple formula illustrating if and where H/V-ratio curves have sharp peaks in dependence of model parameters. In addition, we have constructed a map, which demonstrates the relation between the H/V-peak frequency and the peak frequency of the transmission response in the domain of the layer's Poisson ratio and the impedance contrast. Finally, we have derived maps showing the relationship between the H/V-peak and trough frequency and key parameters of the model such as impedance contrast. These maps are seen as diagnostic tools, which can help to guide the interpretation of H/V spectral ratio diagrams in the context of site effect studies.
Severely negative Delta(14)C anomalies from the mid-depth Pacific and the Arabian Sea have been taken as support for the hypothesized deglacial release of a previously isolated, extremely (14)C-deplete deep ocean carbon reservoir. We report box model simulations that cast doubt on both the existence of the hypothesized deep reservoir and its ability to explain the mid-depth Delta(14)C anomalies. First, the degree of ice age isolation needed to substantially reduce the deep Delta(14)C of the deep reservoir causes anoxia and the trapping of alkalinity from CaCO(3) dissolution, the latter increasing atmospheric CO(2). Second, even with a completely (14)C-free deep reservoir, achieving the mid-depth Delta(14)C anomalies of observed duration requires ad hoc stifling of aspects of deep circulation to prevent rapid dissipation of the anomalous (14)C-free carbon to the rest of the ocean and the atmosphere. We suggest that the mid-depth anomalies do not record basin-scale Delta(14)C changes but are instead local phenomena.
In this paper we evaluate different methods to predict soil erosion processes. We derived different layers of predictor variables for the study area in the Northern Chianti, Italy, describing the soil-lithologic complex, land use, and topographic characteristics. For a subcatchment of the Orme River, we mapped erosion processes by interpreting aerial photographs and field observations. These were classified as erosional response units (ERU), i.e. spatial areas of homogeneous erosion processes. The ERU were used as the response variable in the soil erosion modelling process. We applied two models i) bootstrap aggregation (Random Forest: RF), and ii) stochastic gradient boosting (TreeNet: TN) to predict the potential spatial distribution of erosion processes for the entire Orme River catchment. The models are statistically evaluated using training data and a set of performance parameters such as the area under the receiver operating characteristic curve (AUC), Cohen's Kappa, and pseudo R2. Variable importance and response curves provide further insight into controlling factors of erosion. Both models provided good performance in terms of classification and calibration; however, TN outperformed RF. Similar classes such as active and inactive landslides can be discriminated and well interpreted by considering response curves and relative variable importance. The spatial distribution of the predicted erosion susceptibilities generally follows topographic constraints and is similar for both models. Hence, the model-based delineation of ERU on the basis of soil and terrain information is a valuable tool in geomorphology; it provides insights into factors controlling erosion processes and may allow the extrapolation and prediction of erosion processes in unsurveyed areas.
Tectonic and metamorphic data for the Central Alps (Switzerland/Italy) are used to discuss this classic example of a Barrovian metamorphic terrain, notably the evolution of its thermal structure in space and time. Available P-T-t data indicate variable contributions of advective and conductive heat transport during collision and subsequent cooling and exhumation. Some areas experienced a prolonged period of partial melting while other areas, at the same time, show but moderate heating. The Barrow-type metamorphic field gradient observed in the final orogen is the result of two distinct tectonic processes, with their related advective and conductive heat transport processes. The two tectonic processes are (1) accretion of material within a subduction channel related to decompression and emplacement of high-pressure units in the middle crust and (2) wedging and related nappe formation in the continental lower plate. The second process postdates the first one. Wedging and underthrusting of continental lower plate material produces heat input into lower crustal levels, and this process is responsible for predominantly conductive heat transport in the overlying units. The interacting processes lead to different maximum temperatures at different times, producing the final Barrovian metamorphic field gradient. The south experienced rapid cooling, whereas the north shows moderate cooling rates. This discrepancy principally reflects differences in the temperature distribution in the deeper crust prior to cooling. Differences in the local thermal gradient that prevailed before the cooling also determined the relationships between cooling rate and exhumation rate in the different areas. Citation: Berger, A., S. M. Schmid, M. Engi, R. Bousquet, and M. Wiederkehr (2011), Mechanisms of mass and heat transport during Barrovian metamorphism: A discussion based on field evidence from the Central Alps (Switzerland/northern Italy), Tectonics, 30, TC1007, doi:10.1029/2009TC002622.
A common problem in ecology is identifying the relationship between relief and site properties obtainable only by point measurements. The method of Multi-Scale Landscape Analysis (MSLA) identifies such correlations. MSLA combines frequency filtering of the digital elevation model (DEM) with an estimation of the optimum filter coefficients using an optimization procedure. Tested using point data of soil decarbonation from a German young moraine landscape, MSLA provided significant results. Implemented within open source software SAMT. MSLA is comfortable and flexible to use, offering applications for numerous other spatial analysis problems.
Soil moisture at the plot or hill-slope scale is an important link between local vadose zone hydrology and catchment hydrology. However, so far only a few methods are on the way to close this gap between point measurements and remote sensing. One new measurement methodology that could determine integral soil moisture at this scale is the aboveground sensing of cosmic-ray neutrons, more precisely of ground albedo neutrons. The present study performed ground albedo neutron sensing (GANS) at an agricultural field in northern Germany. To test the method it was accompanied by other soil moisture measurements for a summer period with corn crops growing on the field and a later autumn-winter period without crops and a longer period of snow cover. Additionally, meteorological data and aboveground crop biomass were included in the evaluation. Hourly values of ground albedo neutron sensing showed a high statistical variability. Six-hourly values corresponded well with classical soil moisture measurements, after calibration based on one reference dry period and three wet periods of a few days each. Crop biomass seemed to influence the measurements only to minor degree, opposed to snow cover which has a more substantial impact on the measurements. The latter could be quantitatively related to a newly introduced field neutron ratio estimated from neutron counting rates of two energy ranges. Overall, our study outlines a procedure to apply the ground albedo neutron sensing method based on devices now commercially available, without the need for accompanying numerical simulations and suited for longer monitoring periods after initial calibration.
Ostracod shells in surface sediments from Ulungur Lake (Xinjiang, China) belong mainly to Limnocythere inopinata as the dominant species, and Candona neglecta and Darwinula stevensoni as accompanying, less abundant taxa. Shells of an additional nine species were recorded only sporadically. The three most abundant ostracods have wide tolerance ranges in terms of salinity, substrate and water depth. The similarly recorded bivalve Pisidium subtruncatum, and the gastropods Gyraulus chinensis and Radix auricularia belong to the most tolerant representatives of the genera. The bivalve and gastropods, in addition to the ostracod assemblage, reflect the fact that Ulungur Lake has experienced strong lake level and salinity variations due to water withdrawal in the catchment and the counteracting diversion of river waters to the lake in recent decades. The substrate in Ulungur Lake is typically fine-grained, apart from the delta region of the Ulungur River channel, which is marked by relatively coarse-grained detrital sediments barren of ostracod shells. This channel was created 40 years ago to divert water to Ulungur Lake and support its local fisheries and recreational facilities. A reassessed Holocene ostracod record from the lake shows that a significantly higher salinity and lower lake level existed in the early Holocene before 6.0 ka in response to the regional climate. In contrast, a higher lake level and lowest salinity is inferred for the late Holocene period between ca. 3.6 and 1.3 ka before present. Afterwards, the lake level declined and salinity increased in response to regional moisture reduction, although conditions similar to the early Holocene lake status were not re-established. Our surface-sediment-derived data provide a baseline for analysis of future environmental variations due to global climate change and regional water management.
Hydrocarbons can be found in many different habitats and represent an important carbon source for microbes. As fossil fuels, they are also an important economical resource and through natural seepage or accidental release they can be major pollutants. DNA-specific stains and molecular probes bind to hydrocarbons, causing massive background fluorescence, thereby hampering cell enumeration. The cell extraction procedure of Kallmeyer et al. (2008) separates the cells from the sediment matrix. In principle, this technique can also be used to separate cells from oily sediments, but it was not originally optimized for this application. Here we present a modified extraction method in which the hydrocarbons are removed prior to cell extraction. Due to the reduced background fluorescence the microscopic image becomes clearer, making cell identification, and enumeration much easier. Consequently, the resulting cell counts from oily samples treated according to our new protocol are significantly higher than those treated according to Kallmeyer et al. (2008). We tested different amounts of a variety of solvents for their ability to remove hydrocarbons and found that n-hexane and in samples containing more mature oils methanol, delivered the best results. However, as solvents also tend to lyse cells, it was important to find the optimum solvent to sample ratio, at which hydrocarbon extraction is maximized and cell lysis minimized. A volumetric ratio of 1:2-1:5 between a formalin-fixed sediment slurry and solvent delivered highest cell counts. Extraction efficiency was around 30-50% and was checked on both oily samples spiked with known amounts of E. coli cells and oil-free samples amended with fresh and biodegraded oil. The method provided reproducible results on samples containing very different kinds of oils with regard to their degree of biodegradation. For strongly biodegraded oil MeOH turned out to be the most appropriate solvent, whereas for less biodegraded samples n-hexane delivered best results.
Northwest Europe's largest heather-dominated sandy habitats are located in the nature reserve Luneburger Heide, Germany. Yet, even these appear to be losing their ability to support some of their stenotopic species such as the ladybird spider, Eresus kollari Rossi 1846, and are thus becoming increasingly important for the preservation of these species. The habitat requirements of this endangered spider species were investigated in order to obtain data that will help stabilize the last remnants of the species' population in northwest Germany. Several heathland habitats were surveyed by pitfall trapping during the mate-search period of the males. Two statistical methods were applied: logistic regression and boosted regression trees (BRT). Both methods showed that three habitat variables are of prime relevance in predicting the occurrence of E. kollari: a) thickness of the organic layer (a negative effect), b) soil temperature at a depth of 10 cm, and c) Calluna cover in the herb layer (both have positive effect). Our results show that choppering (removing above-ground biomass and most of O-layer) and burning are likely appropriate heathland management measures for the conservation of E. kollari. Such measures improve the species' habitat quality by creating a heterogenic (small-scaled) heathland structure with suitable microhabitats. As Calluna heathlands show a clear senescence of the dominant heather, it is essential that those habitat patches be conserved. Further measures, such as transfer experiments, are recommended.
Climatic forcing on channel profiles in the eastern cordillera of the Coroico Region, Bolivia
(2011)
Orographic precipitation has a large impact on channel morphology and rock uplift via a positive feedback to erosion. We show that in the Eastern Cordillera of Bolivia, channel concavities reach their highest values where annual precipitation increases in the downstream direction, exceeding 3000 mm. The steepest channels are upstream of this zone of high concavity, where precipitation rates are <1000 mm yr(-1). Channels exhibit graded forms both upstream and downstream of this transient reach. We conclude that the prolonged effect of orographic erosion and related tectonic uplift is the preservation of channels with extreme concavities in the Eastern Cordillera.
The aim of this paper is to estimate the Hurst parameter of Fractional Gaussian Noise (FGN) using Bayesian inference. We propose an estimation technique that takes into account the full correlation structure of this process. Instead of using the integrated time series and then applying an estimator for its Hurst exponent, we propose to use the noise signal directly. As an application we analyze the time series of the Nile River, where we find a posterior distribution which is compatible with previous findings. In addition, our technique provides natural error bars for the Hurst exponent.
Despite the importance of rhizosphere properties for water flow from soil to roots, there is limited quantitative information on the distribution of water in the rhizosphere of plants.
Here, we used neutron tomography to quantify and visualize the water content in the rhizosphere of the plant species chickpea (Cicer arietinum), white lupin (Lupinus albus), and maize (Zea mays) 12 d after planting.
We clearly observed increasing soil water contents (h) towards the root surface for all three plant species, as opposed to the usual assumption of decreasing water content. This was true for tap roots and lateral roots of both upper and lower parts of the root system. Furthermore, water gradients around the lower part of the roots were smaller and extended further into bulk soil compared with the upper part, where the gradients in water content were steeper.
Incorporating the hydraulic conductivity and water retention parameters of the rhizosphere into our model, we could simulate the gradual changes of h towards the root surface, in agreement with the observations. The modelling result suggests that roots in their rhizosphere may modify the hydraulic properties of soil in a way that improves uptake under dry conditions.
Substantial investment in climate change research has led to dire predictions of the impacts and risks to biodiversity. The Intergovernmental Panel on Climate Change fourth assessment report(1) cites 28,586 studies demonstrating significant biological changes in terrestrial systems(2). Already high extinction rates, driven primarily by habitat loss, are predicted to increase under climate change(3-6). Yet there is little specific advice or precedent in the literature to guide climate adaptation investment for conserving biodiversity within realistic economic constraints(7). Here we present a systematic ecological and economic analysis of a climate adaptation problem in one of the world's most species-rich and threatened ecosystems: the South African fynbos. We discover a counterintuitive optimal investment strategy that switches twice between options as the available adaptation budget increases. We demonstrate that optimal investment is nonlinearly dependent on available resources, making the choice of how much to invest as important as determining where to invest and what actions to take. Our study emphasizes the importance of a sound analytical framework for prioritizing adaptation investments(4). Integrating ecological predictions in an economic decision framework will help support complex choices between adaptation options under severe uncertainty. Our prioritization method can be applied at any scale to minimize species loss and to evaluate the robustness of decisions to uncertainty about key assumptions.
Overland flow is an important hydrological pathway in many forests of the humid tropics. Its generation is subject to topographic controls at differing spatial scales. Our objective was to identify such controls on the occurrence of overland flow in a lowland tropical rainforest. To this end, we installed 95 overland flow detectors (OFDs) in four nested subcatchments of the Lutzito catchment on Barro Colorado Island, Panama, and monitored the frequency of overland flow occurrence during 18 rainfall events at each OFD location temporal frequency. For each such location, we derived three non-digital terrain attributes and 17 digital ones, of which 15 were based on Digital Elevation Models (DEMs) of three different resolutions. These attributes then served as input into a Random Forest ensemble tree model to elucidate the importance and partial and joint dependencies of topographic controls for overland flow occurrence.
Lutzito features a high median temporal frequency in overland flow occurrence of 0.421 among OFD locations. However, spatial temporal frequencies of overland flow occurrence vary strongly among these locations and the subcatchments of Lutzito catchment. This variability is best explained by (1) microtopography, (2) coarse terrain sloping and (3) various measures of distance-to-channel, with the contribution of all other terrain attributes being small. Microtopographic features such as concentrated flowlines and wash areas produce highest temporal frequencies, whereas the occurrence of overland flow drops sharply for flow distances and terrain sloping beyond certain threshold values.
Our study contributes to understanding both the spatial controls on overland flow generation and the limitations of terrain attributes for the spatially explicit prediction of overland flow frequencies.
The role of bedrock fractures and rock mass strength is often considered a primary influence on the efficiency of surface processes and the morphology of landscapes. Quantifying bedrock characteristics at hillslope scales, however, has proven difficult. Here, we present a new field-based method for quantifying the depth and apparent density of bedrock fractures within the shallow subsurface based on seismic refraction surveys. We examine variations in subsurface fracture patterns in both Fiordland and the Southern Alps of New Zealand to better constrain the influence of bedrock properties in governing rates and patterns of landslides, as well as the morphology of threshold landscapes. We argue that intense tectonic deformation produces uniform bedrock fracturing with depth, whereas geomorphic processes produce strong fracture gradients focused within the shallow subsurface. Additionally, we argue that hillslope strength and stability are functions of both the intact rock strength and the density of bedrock fractures, such that for a given intact rock strength, a threshold fracture-density exists that delineates between stable and unstable rock masses. In the Southern Alps, tectonic forces have pervasively fractured intrinsically weak rock to the verge of instability, such that the entire rock mass is susceptible to failure and landslides can potentially extend to great depths. Conversely, in Fiordland, tectonic fracturing of the strong intact rock has produced fracture densities less than the regional stability threshold. Therefore, bedrock failure in Fiordland generally occurs only after geomorphic fracturing has further reduced the rock mass strength. This dependence on geomorphic fracturing limits the depths of bedrock landslides to within this geomorphically weakened zone.
Rapid population growth and economic development have led to increased anthropogenic pressures on the Tibetan Plateau, causing significant land cover changes with potentially severe ecological consequences. To assess whether or not these pressures are also affecting the remote montane-boreal lakes on the SE Tibetan Plateau, fossil pollen and diatom data from two lakes were synthesized. The interplay of aquatic and terrestrial ecosystem response was explored in respect to climate variability and human activity over the past 200 years. Nonmetric multidimensional scaling and Procrustes rotation analysis were undertaken to determine whether pollen and diatom responses in each lake were similar and synchronous. Detrended canonical correspondence analysis was used to develop quantitative estimates of compositional species turnover. Despite instrumental evidence of significant climatic warming on the southeastern Plateau, the pollen and diatom records indicate very stable species composition throughout their profiles and show only very subtle responses to environmental changes over the past 200 years. The compositional species turnover (0.36-0.94 SD) is relatively low in comparison to the species reorganizations known from the periods during the mid-and early-Holocene (0.64-1.61 SD) on the SE Plateau, and also in comparison to turnover rates of sediment records from climate-sensitive regions in the circum arctic. Our results indicate that climatically induced ecological thresholds are not yet crossed, but that human activity has an increasing influence, particularly on the terrestrial ecosystem in our study area. Synergistic processes of post-Little Ice Age warming, 20th century climate warming and extensive reforestations since the 19th century have initiated a change from natural oak-pine forests to seminatural, likely less resilient pine-oak forests. Further warming and anthropogenic disturbances would possibly exceed the ecological threshold of these ecosystems and lead to severe ecological consequences.
Classical methods to analyze the surface composition of atmosphereless planetary objects from an orbiter are IR and gamma ray spectroscopy and neutron backscatter measurements. The idea to analyze surface properties with an in-situ instrument has been proposed by Johnson et al. (1998). There, it was suggested to analyze Europa's thin atmosphere with an ion and neutral gas spectrometer. Since the atmospheric components are released by sputtering of the moon's surface, they provide a link to surface composition. Here we present an improved, complementary method to analyze rocky or icy dust particles as samples of planetary objects from which they were ejected. Such particles, generated by the ambient meteoroid bombardment that erodes the surface, are naturally present on all atmosphereless moons and planets. The planetary bodies are enshrouded in clouds of ballistic dust particles, which are characteristic samples of their surfaces. In situ mass spectroscopic analysis of these dust particles impacting onto a detector of an orbiting spacecraft reveals their composition. Recent instrumental developments and tests allow the chemical characterization of ice and dust particles encountered at speeds as low as 1 km/s and an accurate reconstruction of their trajectories. Depending on the sampling altitude, a dust trajectory sensor can trace back the origin of each analyzed grain with about 10 km accuracy at the surface. Since the detection rates are of the order of thousand per orbit, a spatially resolved mapping of the surface composition can be achieved. Certain bodies (e.g., Europa) with particularly dense dust clouds, could provide impact statistics that allow for compositional mapping even on single flybys. Dust impact velocities are in general sufficiently high at orbiters about planetary objects with a radius > 1000 km and with only a thin or no atmosphere. In this work we focus on the scientific benefit of a dust spectrometer on a spacecraft orbiting Earth's Moon as well as Jupiter's Galilean satellites. This 'dust spectrometer' approach provides key chemical and isotopic constraints for varying provinces or geological formations on the surfaces, leading to better understanding of the body's geological evolution.
Population-level effects of global warming result from concurrent direct and indirect processes. They are typically described by physiologically structured population models (PSPMs). Therefore, inverse modelling offers a tool to identify parameters of individual physiological processes through population-level data analysis, e. g. the temperature dependence of growth from size-frequency data of a field population. Here, we make use of experiments under laboratory conditions, in mesocosms and field monitoring to determine the temperature dependence of growth and mortality of Gammarus pulex. We found an optimum temperature for growth of approximately 17 degrees C and a related temperature coefficient, Q(10), of 1.5 degrees C(-1), irrespective of whether we classically fitted individual growth curves or applied inverse modelling based on PSPMs to laboratory data. From a comparison of underlying data sets we conclude that applying inverse modelling techniques to population-level data results in meaningful response parameters for physiological processes if additional temperature-driven effects, including within-population interaction, can be excluded or determined independently. If this is not the case, parameter estimates describe a cumulative response, e. g. comprising temperature-dependent resource dynamics. Finally, fluctuating temperatures in natural habitats increased the uncertainty in parameter values. Here, PSPM should be applied for virtual monitoring in order to determine a sampling scheme that comprises important dates to reduce parameter uncertainty.
Aiming at the stimulation of intrinsic microbial activity, pulses of pure oxygen or pressurized air were recurrently injected into groundwater polluted with chlorobenzene. To achieve well-controlled conditions and intensive sampling, a large, vertical underground tank was filled with the local unconfined sandy aquifer material. In the course of two individual gas injections, one using pure oxygen and one using pressurized air, the mass transfer of individual gas species between trapped gas phase and groundwater was studied. Field data on the dissolved gas composition in the groundwater were combined with a kinetic model on gas dissolution and transport in porous media. Phase mass transfer of individual gas components caused a temporary enrichment of nitrogen, and to a lower degree of methane, in trapped gas leading to the formation of excess dissolved nitrogen levels downgradient from the dissolving gas phase. By applying a novel gas sampling method for dissolved gases in groundwater it was shown that dissolved nitrogen can be used as a partitioning tracer to indicate complete gas dissolution in porous media.
The stable hydrogen isotope ratios (delta D) of leaf wax n-alkanes record valuable information on plant and ecosystem water relations. It remains, however, unknown if leaf wax n-alkane delta D values record only environmental variation during the brief period of time of leaf growth or if leaf wax n-alkane delta D values are affected by environmental variability throughout the entire lifespan of a leaf. To resolve these uncertainties, we irrigated Populus trichocarpa trees with a pulse of deuterium-enriched water and used compound-specific stable hydrogen isotope analyses to test if the applied tracer could be recovered from leaf wax n-alkanes of leaves that were at different stages of their development during the tracer application. Our experiment revealed that only leaf wax n-alkanes from leaves that had developed during the time of the tracer application were affected, while leaves that were already fully matured at the time of the tracer application were not. We conclude from our study that under controlled environmental conditions, leaf wax n-alkanes are synthesized only early in the ontogeny of a leaf. Our experiment has implications for the interpretation of leaf wax n-alkane delta D values in an environmental context, as it suggests that these compounds record only a brief period of the environmental variability that a leaf experiences throughout its life.
Questions: Which are the factors that influence forest and shrubland loss and regeneration and their underlying drivers?
Location: Central Chile, a world biodiversity hotspot.
Methods: Using land-cover data from the years 1975, 1985, 1999 and 2008, we fitted classification trees and multiple logistic regression models to account for the relationship between different trajectories of vegetation change and a range of biophysical and socio-economic factors.
Results: The variables that most consistently showed significant effects on vegetation change across all time-intervals were slope and distance to primary roads. We found that forest and shrubland loss on one side and regeneration on the other often displayed opposite patterns in relation to the different explanatory variables. Deforestation was positively related to distance to primary roads and to distance within forest edges and was favoured by a low insolation and a low slope. In turn, forest regeneration was negatively related to the distance to primary roads and positively to the distance to the nearest forest patch, insolation and slope. Shrubland loss was positively influenced by slope and distance to cities and primary roads and negatively influenced by distance to rivers. Conversely, shrubland regeneration was negatively related to slope, distance to cities and distance to primary roads and positively related to distance from existing forest patches and distance to rivers.
Conclusions: This article reveals how biophysical and socioeconomic factors influence vegetation cover change and the underlying social, political and economical drivers. This assessment provides a basis for management decisions, considering the crucial role of perennial vegetation cover for sustaining biodiversity and ecosystem services.
The study of electromagnetic (EM) field sensitivities is useful for assessing the feasibility of controlled-source electromagnetic (CSEM) surveys. Sensitivity calculations are also a principal building block of EM inversion schemes. Sensitivities are formally given by the derivatives of the EM field components with respect to conductivity. For horizontally layered media, these derivatives can be evaluated analytically, offering advantages in computational efficiency and accuracy over numerical evaluation. We present a complete set of explicit analytic expressions for the EM field sensitivities in 1-D VTI-anisotropic media for horizontal and vertical electric and magnetic dipole sources, and also for finite horizontal electric sources. Since our derivations are based on a formulation for EM fields that is quite general in allowing for sources and receivers at any depth, our sensitivity expressions exhibit the same generality. We verify our expressions by comparison to numerical solutions, and finally present application examples that demonstrate the utility and versatility of these expressions for CSEM feasibility studies.
Clearing for large-scale soy production and the displacement of cattle-breeding by soybeans are major features of land-use change in the lowland Amazon that can alter hydrologic properties of soils and the runoff generation over large areas. We measured infiltrability and saturated hydraulic conductivity (Ksat) under natural forest, pasture, and soybeans on Oxisols in a region of rapid soybean expansion in Mato Grosso, Brazil. The forest-pasture conversion reduced infiltrability from 1258 to 100 mm/h and Ksat at all depths. The pasture-soy conversion increased infiltrability from 100 to 469 mm/h (attributed to shallow disking), did not affect Ksat at 12.5 cm, but decreased Ksat at 30 cm from 122 to 80 mm/h, suggesting that soybean cultivation enhances subsoil compaction. Permeability decreased markedly with depth under forest, did not change under pasture, and averaged out at one fourth the forest value under soybeans with a similar pattern of anisotropy. Comparisons of permeability with rainfall intensities indicated that land-use change did not alter the predominantly vertical water movement within the soil. We conclude that this landscape is well buffered against land-use changes regarding near-surface hydrology, even though short-lived ponding and perched water tables may occur locally during high-intensity rainfall on pastures and under soybeans.
Logic trees have become the most popular tool for the quantification of epistemic uncertainties in probabilistic seismic hazard assessment (PSHA). In a logic-tree framework, epistemic uncertainty is expressed in a set of branch weights, by which an expert or an expert group assigns degree-of-belief values to the applicability of the corresponding branch models. Despite the popularity of logic-trees, however, one finds surprisingly few clear commitments to what logic-tree branch weights are assumed to be (even by hazard analysts designing logic trees). In the present paper we argue that it is important for hazard analysts to accept the probabilistic framework from the beginning for assigning logic-tree branch weights. In other words, to accept that logic-tree branch weights are probabilities in the axiomatic sense, independent of one's preference for the philosophical interpretation of probabilities. We demonstrate that interpreting logic-tree branch weights merely as a numerical measure of "model quality," which are then subsequently normalized to sum up to unity, will with increasing number of models inevitably lead to an apparent insensitivity of hazard curves on the logic-tree branch weights, which may even be mistaken for robustness of the results. Finally, we argue that assigning logic-tree branch weights in a sequential fashion may improve their logical consistency.
Previous studies based on fossil pollen data have reported significant changes in vegetation on the alpine Tibetan Plateau during the Holocene. However, since the relative proportions of fossil pollen taxa are largely influenced by individual pollen productivities and the dispersal characteristics, such inferences on vegetation have the potential to be considerably biased. We therefore examined the modern pollen-vegetation relationships for four common pollen species on the Tibetan Plateau, using Extended R-value (ERV) models. Assuming an average radius of 100 m for the sampled lakes, we estimated the relevant source area of pollen (RSAP) to be 2200 m (which represents the distance from the lake). Using Poaceae as the reference taxa (Pollen Productivity Estimate, PPE = 1), ERV Submodel 2 derived relative high PPEs for the steppe and desert taxa: 2.079 +/- 0.432 for Artemisia and 5.379 +/- 1.077 for Chenopodiaceae. Low PPEs were estimated for the Cyperaceae (1.036 +/- 0.012). whose plants are characteristic of the alpine Kobresia meadows. Applying these PPEs to four fossil pollen sequences since the Late Glacial, the plant abundances on the central and north-eastern Tibetan Plateau were quantified using the "Regional Estimates of Vegetation Abundance from Large Sites" (REVEALS) model. The proportions of Artemisia and Chenopodiaceae were greatly reduced compared to their original pollen percentages in the reconstructed vegetation, owing to their high productivities and their dispersal characteristics, while Cyperaceae showed a relative increase in the vegetation reconstruction. The reconstructed vegetation assemblages of the four pollen sequence sites always yielded smaller compositional species turnovers than suggested by the pollen spectra, as revealed by Detrended Canonical Correspondence Analyses (DCCA) of the Holocene sections. The strength of the previously reported vegetation changes may therefore have been overestimated, which indicates the importance of taking into account pollen-vegetation relationships when discussing the potential drivers (such as climate, land use, atmospheric CO(2) concentrations) and implications (such as for land surface-climate feedbacks, carbon storage, and biodiversity) of vegetation change.
The Karaburun Peninsula, which is considered part of the Anatolide-Tauride Block of Turkey, contains clastic and carbonate sequences deposited on the northern margin of Gondwana. The Palaeozoic clastic sequence, which is intruded by the Early Triassic granitoid and tectonically overlies a Mesozoic melange sequence, can be divided into three subunits: a lower clastic subunit consisting of a sandstone-shale alternation, an upper clastic subunit consisting of black chert-bearing shales, sandstone and conglomerate, and a Permo-Carboniferous carbonate subunit. The lower Triassic Karaburun I-type granitoid has a high initial Sr-87/Sr-86 ratio (0.709021-0.709168), and low Nd-143/Nd-144 ratio (0.512004-0.512023) and epsilon Nd (-5.34 to -5.70) isotopic values. Geochronological data indicate a crystallization (intrusion) age of 247.1 +/- 2.0 Ma (Scythian). Geochemically, the acidic magmatism reflects a subduction-related continental-arc basin tectonic setting, which can be linked to the opening of the northern branch of Neo-Tethys as a continental back-arc rifting basin on the northern margin of Gondwana. This can be related to the closure through southward subduction of the Palaeotethys Ocean beneath Gondwana.
A sediment core from a closed basin lake (Lake Kuhai) from the semi-arid northeastern Tibetan Plateau was analysed for its pollen record to infer Lateglacial and post glacial vegetation and climatic change. At Lake Kuhai five major vegetation and climate shifts could be identified: (1) a change from cold and dry to relatively warmer and more moist conditions at 14.8 cal ka BP: (2) a shift to conditions of higher effective moisture and a stepwise warmer climate at 13.6 cal ka BP; (3) a further shift with increased moisture but colder conditions at 7.0 cal ka BP; (4) a return to a significantly colder and drier phase at 6.3 cal ka BP; (5) and a change back to relatively moist conditions at 2.2 cal ka BP. To investigate the response of lake ecosystems to climatic changes, statistical comparisons were made between the lake Kuhai pollen record and a formerly published ostracod and sedimentary record from the same sediment core. Furthermore, the pollen and lacustrine proxies from lake Kuhai were compared to a previously published pollen and lacustrine record from the nearby Lake Koucha. Statistical comparisons were done using non-metric multidimensional scaling and Procrustes rotation. Differences between lacustrine and pollen responses within one site could be identified, suggesting that lacustrine proxies are partly influenced by in-lake or local catchment processes, whereas the terrestrial (pollen) proxy shows a regional climate signal. Furthermore, we found regional differences in proxy response between lake Kuhai and Lake Koucha. Both pollen records reacted in similar ways to major environmental changes, with minor differences in the timing and magnitude of these changes. The lacustrine records were very similar in their timing and magnitude of response to environmental changes; however, the nature of change was at times very distinct. To place the current study in the context of Holocene moisture evolution across the Tibetan Plateau, we applied a five-scale moisture index and average link clustering to all available continuous palaeo-climate records from the Tibetan Plateau to possibly find general patterns of moisture evolution on the Plateau. However, no common regional pattern of moisture evolution during the Holocene could be detected. We assign this to complex responses of different proxies to environmental and atmospheric changes in an already very heterogeneous mountain landscape where minor differences in elevation can cause strong variation in microenvironments.
The Dead Sea Transform (DST) is a major left-lateral strike-slip fault that accommodates the relative motion between the African and Arabian plates, connecting a region of extension in the Red Sea to the Taurus collision zone in Turkey over a length of about 1100 km. The Dead Sea Basin (DSB) is one of the largest basins along the DST. The DSB is a morphotectonic depression along the DST, divided into a northern and a southern sub-basin, separated by the Lisan salt diapir. We report on a receiver function study of the crust within the multidisciplinary geophysical project, DEad Sea Integrated REsearch (DESIRE), to study the crustal structure of the DSB. A temporary seismic network was operated on both sides of the DSB between 2006 October and 2008 April. The aperture of the network is approximately 60 km in the E-W direction crossing the DSB on the Lisan peninsula and about 100 km in the N-S direction. Analysis of receiver functions from the DESIRE temporary network indicates that Moho depths vary between 30 and 38 km beneath the area. These Moho depth estimates are consistent with results of near-vertical incidence and wide-angle controlled-source techniques. Receiver functions reveal an additional discontinuity in the lower crust, but only in the DSB and west of it. This leads to the conclusion that the internal crustal structure east and west of the DSB is different at the present-day. However, if the 107 km left-lateral movement along the DST is taken into account, then the region beneath the DESIRE array where no lower crustal discontinuity is observed would have lain about 18 Ma ago immediately adjacent to the region under the previous DESERT array west of the DST where no lower crustal discontinuity is recognized.
Quantification of total cell abundance is one of the most fundamental parameters in the exploration of subsurface life. Despite all recent advances in molecular techniques, this parameter is usually determined by fluorescence microscopy. In order to obtain reliable and reproducible results, it is important not just to focus on the actual cell enumeration but also to consider the entire chain of processing. Starting with the retrieval of the sample, over subsampling and sample processing to the final step of fluorescence microscopy, there are many potential sources of contamination that have to be assessed and, if possible, avoided. Because some degree of sample contamination will always occur, it is necessary to employ some form of contamination control. Different tracers are available, each one with its specific advantages and drawbacks. In many cases, the problems arise not after the sample has arrived in a well-equipped laboratory with highly trained personnel, but much earlier at the drill site or in a field camp. In this review, I discuss the different aspects of cell enumeration in subsurface sediment, evaluating every step in the long process chain.
Coastal uplift and tsunami effects associated to the 2010 M(w)8.8 Maule earthquake in Central Chile
(2011)
On February 27, 2010 at 03:34:08 AM an M(w)8.8 earthquake, with epicenter located off Cobquecura (73.24 degrees W; 36.29 degrees S), severely hit Central Chile. The tsunami waves that followed this event affected the coastal regions between the cities of Valparaiso and Valdivia, with minor effects as far as Coquimbo. The earthquake occurred along the subduction of the Nazca oceanic plate beneath the South American plate. Coseismic coastal uplift was estimated through observations of bleached lithothamnioids crustose coralline algae, which were exposed after the mainshock between 34.13 degrees S and 38.34 degrees S, suggesting the latitudinal distribution of the earthquake rupture. The measured coastal uplift values varied between 240 +/- 20 cm at sites closer to the trench along the western coast of the Arauco peninsula and 15 +/- 10 cm at sites located farther east. A maximum value of 260 +/- 50 cm was observed at the western coast of Santa Maria Island, which is similar to the reported uplift associated with the 1835 earthquake at Concepcion. Land subsidence values on the order of 0.5 m to 1 m evidenced a change in polarity and position of the coseismic hinge at 110-120 km from the trench. In four sites along the coast we observed a close match between coastal uplift values deduced from bleached lithothamnioids algae and GPS measurements. According to field observations tsunami heights reached ea. 14 m in the coastal area of the Maule Region immediately north of the epicenter, and diminished progressively northwards to 4-2 m near Valparaiso. Along the coast of Cobquecura, tsunami height values were inferior to 2-4 m. More variable tsunami heights of 6-8 m were measured at Dichato-Talcahuano and Tirua-Puerto Saavedra, in the Biobio and Arauco regions, respectively, to the south of the epicenter. According to eyewitnesses, the tsunami reached the coast between 12 to 20 and 30 to 45 minutes in areas located closer and faraway from the earthquake rupture zone, respectively. Destructive tsunami waves arrived also between 2.5 and 4.5 hours after the mainshock, especially along the coast of the Biobio and Arauco regions. The tsunami effects were highly variable along the coast, as a result of geomorphological and bathymetric local conditions, besides potential complexities induced by the main shock.
The Mahneshan Metamorphic Complex (MMC) is one of the Precambrian terrains exposed in the northwest of Iran. The MMC underwent two main phases of deformation (D-1 and D-2) and at least two metamorphic events (M-1 and M-2). Critical metamorphic mineral assemblages in the metapelitic rocks testify to regional metamorphism under amphibolite-facies conditions. The dominant metamorphic mineral assemblage in metapelitic rocks (M-1) is muscovite, biotite I, Garnet I, staurolite, Andalusite I and sillimanite. Peak metamorphism took place at 600-620 degrees C and similar to 7 kbar, corresponding to a depth of ca. 24 km. This was followed by decompression during exhumation of the crustal rocks up to the surface. The decrease of temperature and pressure during exhumation produced retrograde metamorphic assemblages (M-2). Secondary phases such as garnet II biotite It. Andalusite II constrain the temperature and pressure of M, retrograde metamorphism to 520-560 degrees C and 2.5-3.5 kbar, respectively. The geothermal gradient obtained for the peak of metamorphism is 33 degrees C km(-1), which indicates that peak metamorphism was of Barrovian type and occurred under medium-pressure conditions. The MMC followed a 'clockwise' P T path during metamorphism, consistent with thermal relaxation following tectonic thickening. The bulk chemistry of the MMC metapelites shows that their protoliths were deposited at an active continental margin. Together with the presence of palaeo-suture zones and ophiolitic rocks around the high-grade metamorphic rocks of the MMC, these features suggest that the Iranian Precambrian basement formed by an island-arc type cratonization.
Well-preserved primary contact relationships between a Late Proterozoic metasedimentary and the metagranitic core and Palaeozoic cover series of the Menderes Massif have been recognized in the eastern part of the Cine submassif on a regional-scale. Metaconglomerates occur as laterally discontinuous channel-fill bodies close the base of the metaquartzarenite directly above the basement. The pebbles in the metaconglomerates consist mainly of different types of tourmaline-rich leucocratic granitoids, tourmalinite and schist in a sandy matrix. Petrographic features, geochemical compositions and zircon radiometric ages (549.6 +/- 3.7-552.3 +/- 3.1 Ma) of the diagnostic clasts of the metaconglomerates (e.g. leucocratic granitoids and tourmalinites) show excellent agreement with their in situ equivalents (549.0 +/- 5.4 Ma) occurring in the Pan-African basement as stocks and veins.
The correlation between clasts in the metaconglomerates and granitoids of the basement suggests that the primary contact between the basement and cover series is a regional unconformity (supra-Pan-African Unconformity) representing deep erosion of the Pan-African basement followed by the deposition of the cover series. Hence the usage of 'core-cover' terminology in the Menderes Massif is valid. Consequently, these new data preclude the views that the granitic precursors of the leucocratic orthogneisses are Tertiary intrusions.
Wavelet modelling of the gravity field by domain decomposition methods: an example over Japan
(2011)
With the advent of satellite gravity, large gravity data sets of unprecedented quality at low and medium resolution become available. For local, high resolution field modelling, they need to be combined with the surface gravity data. Such models are then used for various applications, from the study of the Earth interior to the determination of oceanic currents. Here we show how to realize such a combination in a flexible way using spherical wavelets and applying a domain decomposition approach. This iterative method, based on the Schwarz algorithms, allows to split a large problem into smaller ones, and avoids the calculation of the entire normal system, which may be huge if high resolution is sought over wide areas. A subdomain is defined as the harmonic space spanned by a subset of the wavelet family. Based on the localization properties of the wavelets in space and frequency, we define hierarchical subdomains of wavelets at different scales. On each scale, blocks of subdomains are defined by using a tailored spatial splitting of the area. The data weighting and regularization are iteratively adjusted for the subdomains, which allows to handle heterogeneity in the data quality or the gravity variations. Different levels of approximations of the subdomains normals are also introduced, corresponding to building local averages of the data at different resolution levels.
We first provide the theoretical background on domain decomposition methods. Then, we validate the method with synthetic data, considering two kinds of noise: white noise and coloured noise. We then apply the method to data over Japan, where we combine a satellite-based geopotential model, EIGEN-GL04S, and a local gravity model from a combination of land and marine gravity data and an altimetry-derived marine gravity model. A hybrid spherical harmonics/wavelet model of the geoid is obtained at about 15 km resolution and a corrector grid for the surface model is derived.
We study changes in effective stress (normal stress minus pore pressure) that occurred in the French Alps during the 2003-2004 Ubaye earthquake swarm. Two complementary data sets are used. First, a set of 974 relocated events allows us to finely characterize the shape of the seismogenic area and the spatial migration of seismicity during the crisis. Relocations are performed by a double-difference algorithm. We compute differences in travel times at stations both from absolute picking times and from cross-correlation delays of multiplets. The resulting catalog reveals a swarm alignment along a single planar structure striking N130 degrees E and dipping 80 degrees W. This relocated activity displays migration properties consistent with a triggering by a diffusive fluid overpressure front. This observation argues in favor of a deep-seated fluid circulation responsible for a significant part of the seismic activity in Ubaye. Second, we analyze time series of earthquake detections at a single seismological station located just above the swarm. This time series forms a dense chronicle of +16,000 events. We use it to estimate the history of effective stress changes during this sequence. For this purpose we model the rate of events by a stochastic epidemic-type aftershock sequence model with a nonstationary background seismic rate lambda(0)(t). This background rate is estimated in discrete time windows. Window lengths are determined optimally according to a new change-point method on the basis of the interevent times distribution. We propose that background events are triggered directly by a transient fluid circulation at depth. Then, using rate-and-state constitutive friction laws, we estimate changes in effective stress for the observed rate of background events. We assume that changes in effective stress occurred under constant shear stressing rate conditions. We finally obtain a maximum change in effective stress close to -8 MPa, which corresponds to a maximum fluid overpressure of about 8 MPa under constant normal stress conditions. This estimate is in good agreement with values obtained from numerical modeling of fluid flow at depth, or with direct measurements reported from fluid injection experiments.