@article{ZurellBergerCabraletal.2010, author = {Zurell, Damaris and Berger, Uta and Cabral, Juliano Sarmento and Jeltsch, Florian and Meynard, Christine N. and Muenkemueller, Tamara and Nehrbass, Nana and Pagel, J{\"o}rn and Reineking, Bjoern and Schroeder, Boris and Grimm, Volker}, title = {The virtual ecologist approach : simulating data and observers}, issn = {0030-1299}, doi = {10.1111/j.1600-0706.2009.18284.x}, year = {2010}, abstract = {Ecologists carry a well-stocked toolbox with a great variety of sampling methods, statistical analyses and modelling tools, and new methods are constantly appearing. Evaluation and optimisation of these methods is crucial to guide methodological choices. Simulating error-free data or taking high-quality data to qualify methods is common practice. Here, we emphasise the methodology of the 'virtual ecologist' (VE) approach where simulated data and observer models are used to mimic real species and how they are 'virtually' observed. This virtual data is then subjected to statistical analyses and modelling, and the results are evaluated against the 'true' simulated data. The VE approach is an intuitive and powerful evaluation framework that allows a quality assessment of sampling protocols, analyses and modelling tools. It works under controlled conditions as well as under consideration of confounding factors such as animal movement and biased observer behaviour. In this review, we promote the approach as a rigorous research tool, and demonstrate its capabilities and practical relevance. We explore past uses of VE in different ecological research fields, where it mainly has been used to test and improve sampling regimes as well as for testing and comparing models, for example species distribution models. We discuss its benefits as well as potential limitations, and provide some practical considerations for designing VE studies. Finally, research fields are identified for which the approach could be useful in the future. We conclude that VE could foster the integration of theoretical and empirical work and stimulate work that goes far beyond sampling methods, leading to new questions, theories, and better mechanistic understanding of ecological systems.}, language = {en} } @phdthesis{Zubaidah2010, author = {Zubaidah, Teti}, title = {Spatio-temporal characteristics of the geomagnetic field over the Lombok Island, the Lesser Sunda Islands region}, series = {Scientific Technical Report}, volume = {STR10}, journal = {Scientific Technical Report}, number = {07}, publisher = {Deutsches GeoForschungsZentrum GFZ}, address = {Potsdam}, doi = {10.2312/GFZ.b103-10079}, school = {Universit{\"a}t Potsdam}, pages = {xv, 117}, year = {2010}, abstract = {The Lombok Island is part of the Lesser Sunda Islands (LSI) region - Indonesia, situated along the Sunda-Banda Arcs transition. It lies between zones characterized by the highest intensity geomagnetic anomalies of this region, remarkable as one of the eight most important features provided on the 1st edition of World Digital Magnetic Anomaly Map. The seismicity of this region during the last years is high, while the geological and tectonic structures of this region are still not known in detail. Some local magnetic surveys have been conducted previously during 2004-2005. However, due to the lower accuracy of the used equipment and a limited number of stations, the qualities of the previous measurements are questionable for more interpretations. Thus a more detailed study to better characterize the geomagnetic anomaly -spatially and temporally- over this region and to deeply explore the related regional geology, tectonic and seismicity is needed. The intriguing geomagnetic anomalies over this island region vis-{\`a}-vis the socio-cultural situations lead to a study with a special aim to contribute to the assessment of the potential of natural hazards (earthquakes) as well as a new natural resource of energy (geothermal potential). This study is intended to discuss several crucial questions, including: i. The real values and the general pattern of magnetic anomalies over the island, as well as their relation to the regional one. ii. Any temporal changes of regional anomalies over the recent time. iii. The relationships between the anomalies and the geology and tectonic of this region, especially new insights that can be gained from the geomagnetic observations. iv. The relationships between the anomalies and the high seismicity of this region, especially some possible links between their variations to the earthquake occurrence. First, all available geomagnetic data of this region and results of the previous measurements are evaluated. The new geomagnetic surveys carried out in 2006 and 2007/2008 are then presented in detail, followed by the general description of data processing and data quality evaluation. The new results show the general pattern of contiguous negative-positive anomalies, revealing an active arc related subduction region. They agree with earlier results obtained by satellite, aeromagnetic, and marine platforms; and provide a much more detailed picture of the strong anomalies on this island. The temporal characteristics of regional anomalies show a decreasing strength of the dipolar structure, where decreasing of the field intensities is faster than the regional secular variations as defined by the global model (the 10th generation of IGRF). However, some exceptions (increasing of anomalies) have to be noted and further analyzed for several locations. Thereafter, simultaneous magnetic anomalies and gravity models are generated and interpreted in detail. Three profiles are investigated, providing new insights into the tectonics and geological evolution of the Lombok Island. Geological structure of this island can be divided as two main parts with different consecutive ages: an old part (from late Oligocene to late Miocene) in the South and a younger one (from Pliocene to Holocene) in the North. A new subduction in the back arc region (the Flores Thrust zone) is considered mature and active, showing a tendency of progressive subduction during 2005-2008. Geothermal potential in the northern part of this island can be mapped in more detail using these geomagnetic regional survey data. The earlier estimates of reservoir depth can be confirmed further to a depth of about 800 m. Evaluation of temporal changes of the anomalies gives some possible explanations related to the evolution of the back arc region, large stress accumulations over the LSI region, a specific electrical characteristic of the crust of the Lombok Island region, and a structural discontinuity over this island. Based on the results, several possible advanced studies involving geomagnetic data and anomaly investigations over the Lombok Island region can be suggested for the future: i. Monitoring the subduction activity of the back arc region (the Flores Thrust zone) and the accumulated stress over the LSI, that could contribute to middle term hazard assessment with a special attention to the earthquake occurrence in this region. Continuous geomagnetic field measurements from a geomagnetic observatory which can be established in the northern part of the Lombok Island and systematic measurements at several repeat stations can be useful in this regards. ii. Investigating the specific electrical characteristic (high conductivity) of the crust, that is probably related to some aquifer layers or metal mineralization. It needs other complementary geophysical methods, such as magnetotelluric (MT) or preferably DC resistivity measurements. iii. Determining the existence of an active structural fault over the Lombok Island, that could be related to long term hazard assessment over the LSI region. This needs an extension of geomagnetic investigations over the neighbouring islands (the Bali Island in the West and the Sumbawa Island in the East; probably also the Sumba and the Flores islands). This seems possible because the regional magnetic lineations might be used to delineate some structural discontinuities, based on the modelling of contrasts in crustal magnetizations.}, language = {en} } @article{ZimmermannZimmermannLarketal.2010, author = {Zimmermann, Beate and Zimmermann, Alexander and Lark, Richard M. and Elsenbeer, Helmut}, title = {Sampling procedures for throughfall monitoring : a simulation study}, issn = {0043-1397}, doi = {10.1029/2009wr007776}, year = {2010}, abstract = {What is the most appropriate sampling scheme to estimate event-based average throughfall? A satisfactory answer to this seemingly simple question has yet to be found, a failure which we attribute to previous efforts' dependence on empirical studies. Here we try to answer this question by simulating stochastic throughfall fields based on parameters for statistical models of large monitoring data sets. We subsequently sampled these fields with different sampling designs and variable sample supports. We evaluated the performance of a particular sampling scheme with respect to the uncertainty of possible estimated means of throughfall volumes. Even for a relative error limit of 20\%, an impractically large number of small, funnel-type collectors would be required to estimate mean throughfall, particularly for small events. While stratification of the target area is not superior to simple random sampling, cluster random sampling involves the risk of being less efficient. A larger sample support, e.g., the use of trough-type collectors, considerably reduces the necessary sample sizes and eliminates the sensitivity of the mean to outliers. Since the gain in time associated with the manual handling of troughs versus funnels depends on the local precipitation regime, the employment of automatically recording clusters of long troughs emerges as the most promising sampling scheme. Even so, a relative error of less than 5\% appears out of reach for throughfall under heterogeneous canopies. We therefore suspect a considerable uncertainty of input parameters for interception models derived from measured throughfall, in particular, for those requiring data of small throughfall events.}, language = {en} } @article{ZimmermannPapritzElsenbeer2010, author = {Zimmermann, Beate and Papritz, Andreas and Elsenbeer, Helmut}, title = {Asymmetric response to disturbance and recovery : Changes of soil permeability under forest-pasture-forest transitions}, issn = {0016-7061}, doi = {10.1016/j.geoderma.2010.07.013}, year = {2010}, abstract = {In the humid tropics, continuing high deforestation rates are seen alongside an increasing expansion of secondary forests. In order to understand and model the consequences of these dynamic land-use changes for regional water cycles, the response of soil hydraulic properties to forest disturbance and recovery has to be quantified.At a site in the Brazilian Amazonia, we annually monitored soil infiltrability and saturated hydraulic conductivity (K-s) at 12.5, 20 cm, and 50 cm soil depth after manual forest conversion to pasture (year zero to four after pasture establishment), and during secondary succession after pasture abandonment (year zero to seven after pasture abandonment). We evaluated the hydrological consequences of the detected changes by comparing the soil hydraulic properties with site-specific rainfall intensities and hydrometric observations. Within one year after grazing started, infiltrability and K-s at 12.5 and 20 cm depth decreased by up to one order of magnitude to levels which are typical for 20-year-old pasture. In the three subsequent monitoring years, infiltrability and K-s remained stable. Land use did not impact on subsoil permeability. Whereas infiltrability values are large enough to allow all rainwater to infiltrate even after the conversion, the sudden decline of near-surface K-s is of hydrological relevance as perched water tables and overland flow occur more often on pastures than in forests at our study site. After pasture abandonment and during secondary succession, seven years of recovery did not suffice to significantly increase infiltrability and K-s at 12.5 depth although a slight recovery is obvious. At 20 cm soil depth, we detected a positive linear increase within the seven-year time frame but annual means did not differ significantly. Although more than a doubling of infiltrability and K-s is still required to achieve pre-disturbance levels, which will presumably take more than a decade, the observed slight increases of K-s might already decrease the probability of perched water table generation and overland flow development well before complete recovery.}, language = {en} } @article{ZeheGraeffMorgneretal.2010, author = {Zehe, Erwin and Gr{\"a}ff, Thomas and Morgner, Markus and Bauer, Andreas and Bronstert, Axel}, title = {Plot and field scale soil moisture dynamics and subsurface wetness control on runoff generation in a headwater in the Ore Mountains}, issn = {1027-5606}, doi = {10.5194/hess-14-873-2010}, year = {2010}, abstract = {This study presents an application of an innovative sampling strategy to assess soil moisture dynamics in a headwater of the Weisseritz in the German eastern Ore Mountains. A grassland site and a forested site were instrumented with two Spatial TDR clusters (STDR) that consist of 39 and 32 coated TDR probes of 60 cm length. Distributed time series of vertically averaged soil moisture data from both sites/ensembles were analyzed by statistical and geostatistical methods. Spatial variability and the spatial mean at the forested site were larger than at the grassland site. Furthermore, clustering of TDR probes in combination with long-term monitoring allowed identification of average spatial covariance structures at the small field scale for different wetness states. The correlation length of soil water content as well as the sill to nugget ratio at the grassland site increased with increasing average wetness and but, in contrast, were constant at the forested site. As soil properties at both the forested and grassland sites are extremely variable, this suggests that the correlation structure at the forested site is dominated by the pattern of throughfall and interception. We also found a very strong correlation between antecedent soil moisture at the forested site and runoff coefficients of rainfall-runoff events observed at gauge Rehefeld. Antecedent soil moisture at the forest site explains 92\% of the variability in the runoff coefficients. By combining these results with a recession analysis we derived a first conceptual model of the dominant runoff mechanisms operating in this catchment. Finally, we employed a physically based hydrological model to shed light on the controls of soil- and plant morphological parameters on soil average soil moisture at the forested site and the grassland site, respectively. A homogeneous soil setup allowed, after fine tuning of plant morphological parameters, most of the time unbiased predictions of the observed average soil conditions observed at both field sites. We conclude that the proposed sampling strategy of clustering TDR probes is suitable to assess unbiased average soil moisture dynamics in critical functional units, in this case the forested site, which is a much better predictor for event scale runoff formation than pre-event discharge. Long term monitoring of such critical landscape elements could maybe yield valuable information for flood warning in headwaters. We thus think that STDR provides a good intersect of the advantages of permanent sampling and spatially highly resolved soil moisture sampling using mobile rods.}, language = {en} } @book{ZangStephansson2010, author = {Zang, Arno and Stephansson, Ove}, title = {Stress field of the earth's crust}, publisher = {Springer Science+Business Media B.V}, address = {Dordrecht}, isbn = {978-1-402-08443-0}, doi = {10.1007/978-1-4020-8444-7}, pages = {240 S.}, year = {2010}, language = {en} } @article{WulfBookhagenScherler2010, author = {Wulf, Hendrik and Bookhagen, Bodo and Scherler, Dirk}, title = {Seasonal precipitation gradients and their impact on fluvial sediment flux in the Northwest Himalaya}, issn = {0169-555X}, doi = {10.1016/j.geomorph.2009.12.003}, year = {2010}, abstract = {Precipitation in the form of rain and snowfall throughout the Himalaya controls river discharge and erosional processes and, thus, has a first-order control on the fluvial sediment flux. Here, we analyze daily precipitation data (1998-2007) of 80 weather stations from the northwestern Himalaya in order to decipher temporal and spatial moisture gradients. In addition, suspended sediment data allow assessment of the impact of precipitation on the fluvial sediment flux for a 10(3)-km(2) catchment (Baspa). We find that weather stations located at the mountain front receive similar to 80\% of annual precipitation during summer (May-Oct), whereas stations in the orogenic interior, i.e., leeward of the orographic barrier, receive similar to 60\% of annual precipitation during winter (Nov-Apr). In both regions 4-6 rainstorm days account for similar to 40\% of the summer budgets, while rainstorm magnitude-frequency relations, derived from 40-year precipitation time-series, indicate a higher storm variability in the interior than in the frontal region. This high variability in maximum annual rainstorm days in the orogenic interior is reflected by a high variability in extreme suspended sediment events in the Baspa Valley, which strongly affect annual erosion yields. The two most prominent 5-day-long erosional events account for 50\% of the total 5-year suspended sediment flux and coincide with synoptic-scale monsoonal rainstorms. This emphasizes the erosional impact of the Indian Summer Monsoon as the main driving force for erosion processes in the orogenic interior, despite more precipitation falling during the winter season.}, language = {en} } @article{WilkeAppelVinczeetal.2010, author = {Wilke, Max and Appel, Karen and Vincze, Laszlo and Schmidt, Christian and Borchert, Manuela and Pascarelli, Sakura}, title = {A confocal set-up for micro-XRF and XAFS experiments using diamond-anvil cells}, issn = {0909-0495}, doi = {10.1107/S0909049510023654}, year = {2010}, abstract = {A confocal set-up is presented that improves micro-XRF and XAFS experiment with high-pressure e diamond-anvil cells (DACs) In this experiment a probing volume is defined by the focus of the incoming synchrotron radiation beam and that of a polycapillary X-ray half-lens with a very long working distance, which is placed in front of the fluorescence detector This set-up enhances the quality of the fluorescence and XAFS spectra, and thus the sensitivity for detecting elements at low concentrations. It efficiently suppresses signal from outside the sample chamber, which stems from elastic and inelastic scattering of the incoming beam by the diamond anvils as well as from excitation of fluorescence from the body of the DAC}, language = {en} } @article{WilkeO'BrienGerdesetal.2010, author = {Wilke, Franziska Daniela Helena and O'Brien, Patrick J. and Gerdes, Axel and Timmerman, Martin Jan and Sudo, Masafumi and Khan, M. Ahmed}, title = {The multistage exhumation history of the Kaghan Valley UHP series, NW Himalaya, Pakistan from U-Pb and Ar-40/Ar- 39 ages}, issn = {0935-1221}, doi = {10.1127/0935-1221/2010/0022-2051}, year = {2010}, abstract = {Amphibole and mica Ar-40/Ar-39 ages as well as zircon, rutile and titanite U-Pb geochronology of eclogites and associated host rocks from the Higher Himalayan Crystalline Nappes (Indian Plate) in the Upper Kaghan Valley, Pakistan allow distinction of a multistage exhumation history. An Eocene age for peak-pressure metamorphism has been obtained by phengite Ar-40/Ar-39 (47.3 +/- 0.3 Ma) and zircon U-Pb (47.3 +/- 0.4 and 47.4 +/- 0.3 Ma) ages from cover and basement gneisses. A very short-lived metamorphic peak and rapid cooling is documented by an amphibole Ar-40/Ar-39 age of 46.6 +/- 0.5 Ma and a rutile U-Pb age of 44.1 +/- 1.3 Ma from eclogites. Phengite and biotite ages from cover and basement sequences metamorphosed during the Himalayan orogeny are 34.5 +/- 0.2 to 28.1 +/- 0.2 Ma whereas youngest biotites, yielding 23.6 +/- 0.1 and 21.7 +/- 0.2 Ma, probably reflect argon partial resetting. The amphibole age, together with those derived from phengite and zircon demonstrate a rate of initial exhumation of 86-143 mm/a i.e. an extremely rapid transport of the Indian Plate continental crust from ultra-high pressure (UHP) conditions back to crustal levels (47-46 Ma for transport from 140 to 40 km depth). Subsequent exhumation (46-41 Ma, 40-35 km) slowed to about 1 mm/a at the base of the continental crust but increased again later towards slightly higher exhumation rates of ca. 2 mm/a (41-34 Ma, 35- 20 km). This indicates a change from buoyancy-driven exhumation at mantle depths to compression forces related to continent-continent collision and accompanied crustal folding, thrusting and stacking that finally exposed the former deeply-buried rocks.}, language = {en} } @article{WilkeO'BrienAltenbergeretal.2010, author = {Wilke, Franziska Daniela Helena and O'Brien, Patrick J. and Altenberger, Uwe and Konrad-Schmolke, Matthias and Khan, M. Ahmed}, title = {Multi-stage reaction history in different eclogite types from the Pakistan Himalaya and implications for exhumation processes}, issn = {0024-4937}, doi = {10.1016/j.lithos.2009.07.015}, year = {2010}, abstract = {Metabasites were sampled from rock series of the subducted margin of the Indian Plate, the so-called Higher Himalayan Crystalline, in the Upper Kaghan Valley, Pakistan. These vary from corona dolerites, cropping out around Saif- ul-Muluk in the south, to coesite-eclogite close to the suture zone against rocks of the Kohistan arc in the north. Bulk rock major- and trace-element chemistry reveals essentially a single protolith as the source for five different eclogite types, which differ in fabric, modal mineralogy as well as in mineral chemistry. The study of newly-collected samples reveals coesite (confirmed by in situ Raman spectroscopy) in both garnet and omphacite. All eclogites show growth of amphiboles during exhumation. Within some coesite-bearing eclogites the presence of glaucophane cores to barroisite is noted whereas in most samples porphyroblastic sodic-calcic amphiboles are rimmed by more aluminous calcic amphibole (pargasite, tschermakite, and edenite). Eclogite facies rutile is replaced by ilmenite which itself is commonly surrounded by titanite. In addition, some eclogite bodies show leucocratic segregations containing phengite, quartz, zoisite and/or kyanite. The important implication is that the complex exhumation path shows stages of initial cooling during decompression (formation of glaucophane) followed by reheating: a very similar situation to that reported for the coesite-bearing eclogite series of the Tso Morari massif, India, 450 km to the south-east.}, language = {en} } @phdthesis{Wilke2010, author = {Wilke, Franziska Daniela Helena}, title = {Quantifying crystalline exhumation in the Himalaya}, doi = {10.25932/publishup-4119}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-43138}, school = {Universit{\"a}t Potsdam}, pages = {IV, 99}, year = {2010}, abstract = {In 1915, Alfred Wegener published his hypotheses of plate tectonics that revolutionised the world for geologists. Since then, many scientists have studied the evolution of continents and especially the geologic structure of orogens: the most visible consequence of tectonic processes. Although the morphology and landscape evolution of mountain belts can be observed due to surface processes, the driving force and dynamics at lithosphere scale are less well understood despite the fact that rocks from deeper levels of orogenic belts are in places exposed at the surface. In this thesis, such formerly deeply-buried (ultra-) high-pressure rocks, in particular eclogite facies series, have been studied in order to reveal details about the formation and exhumation conditions and rates and thus provide insights into the geodynamics of the most spectacular orogenic belt in the world: the Himalaya. The specific area investigated was the Kaghan Valley in Pakistan (NW Himalaya). Following closure of the Tethyan Ocean by ca. 55-50 Ma, the northward subduction of the leading edge of India beneath the Eurasian Plate and subsequent collision initiated a long-lived process of intracrustal thrusting that continues today. The continental crust of India - granitic basement, Paleozoic and Mesozoic cover series and Permo-Triassic dykes, sills and lavas - has been buried partly to mantle depths. Today, these rocks crop out as eclogites, amphibolites and gneisses within the Higher Himalayan Crystalline between low-grade metamorphosed rocks (600-640°C/ ca. 5 kbar) of the Lesser Himalaya and Tethyan sediments. Beside tectonically driven exhumation mechanisms the channel flow model, that describes a denudation focused ductile extrusion of low viscosity material developed in the middle to lower crust beneath the Tibetan Plateau, has been postulated. To get insights into the lithospheric and crustal processes that have initiated and driven the exhumation of this (ultra-) high-pressure rocks, mineralogical, petrological and isotope-geochemical investigations have been performed. They provide insights into 1) the depths and temperatures to which these rocks were buried, 2) the pressures and temperatures the rocks have experienced during their exhumation, 3) the timing of these processes 4) and the velocity with which these rocks have been brought back to the surface. In detail, through microscopical studies, the identification of key minerals, microprobe analyses, standard geothermobarometry and modelling using an effective bulk rock composition it has been shown that published exhumation paths are incomplete. In particular, the eclogites of the northern Kaghan Valley were buried to depths of 140-100 km (36-30 kbar) at 790-640°C. Subsequently, cooling during decompression (exhumation) towards 40-35 km (17-10 kbar) and 630-580°C has been superseded by a phase of reheating to about 720-650°C at roughly the same depth before final exhumation has taken place. In the southern-most part of the study area, amphibolite facies assemblages with formation conditions similar to the deduced reheating phase indicate a juxtaposition of both areas after the eclogite facies stage and thus a stacking of Indian Plate units. Radiometric dating of zircon, titanite and rutile by U-Pb and amphibole and micas by Ar-Ar reveal peak pressure conditions at 47-48 Ma. With a maximum exhumation rate of 14 cm/a these rocks reached the crust-mantle boundary at 40-35 km within 1 Ma. Subsequent exhumation (46-41 Ma, 40-35 km) decelerated to ca. 1 mm/a at the base of the continental crust but rose again to about 2 mm/a in the period of 41-31 Ma, equivalent to 35-20 km. Apatite fission track (AFT) and (U-Th)/He ages from eclogites, amphibolites, micaschists and gneisses yielded moderate Oligocene to Miocene cooling rates of about 10°C/Ma in the high altitude northern parts of the Kaghan Valley using the mineral-pair method. AFT ages are of 24.5±3.8 to 15.6±2.1 Ma whereas apatite (U-Th)/He analyses yielded ages between 21.0±0.6 and 5.3±0.2 Ma. The southern-most part of the Valley is dominated by younger late Miocene to Pliocene apatite fission track ages of 7.6±2.1 and 4.0±0.5 Ma that support earlier tectonically and petrologically findings of a juxtaposition and stack of Indian Plate units. As this nappe is tectonically lowermost, a later distinct exhumation and uplift driven by thrusting along the Main Boundary Thrust is inferred. A multi-stage exhumation path is evident from petrological, isotope-geochemical and low temperature thermochronology investigations. Buoyancy driven exhumation caused an initial rapid exhumation: exhumation as fast as recent normal plate movements (ca. 10 cm/a). As the exhuming units reached the crust-mantle boundary the process slowed down due to changes in buoyancy. Most likely, this exhumation pause has initiated the reheating event that is petrologically evident (e.g. glaucophane rimmed by hornblende, ilmenite overgrowth of rutile). Late stage processes involved widespread thrusting and folding with accompanied regional greenschist facies metamorphism, whereby contemporaneous thrusting on the Batal Thrust (seen by some authors equivalent to the MCT) and back sliding of the Kohistan Arc along the inverse reactivated Main Mantle Thrust caused final exposure of these rocks. Similar circumstances have been seen at Tso Morari, Ladakh, India, 200 km further east where comparable rock assemblages occur. In conclusion, as exhumation was already done well before the initiation of the monsoonal system, climate dependent effects (erosion) appear negligible in comparison to far-field tectonic effects.}, language = {en} } @article{WichuraBousquetOberhaenslietal.2010, author = {Wichura, Henry and Bousquet, Romain and Oberh{\"a}nsli, Roland and Strecker, Manfred and Trauth, Martin H.}, title = {Evidence for middleUocene uplift of the East African Plateau}, issn = {0091-7613}, doi = {10.1130/G31022.1}, year = {2010}, abstract = {Cenozoic uplift of the East African Plateau has been associated with fundamental climatic and environmental changes in East Africa and adjacent regions. While this influence is widely accepted, the timing and the magnitude of plateau uplift have remained unclear. This uncertainty stems from the lack of datable, geomorphically meaningful reference horizons that could record surface uplift. Here, we document the existence of significant relief along the East African Plateau prior to rifting, as inferred from modeling the emplacement history of one of the longest terrestrial lava flows, the similar to 300-km-long Yatta phonolite flow in Kenya. This 13.5 Ma lava flow originated on the present-day eastern Kenya Rift flank, and utilized a riverbed that once routed runoff from the eastern rim of the plateau. Combining an empirical viscosity model with subsequent cooling and using the Yatta lava flow geometry and underlying paleotopography (slope angle), we found that the prerift slope was at least 0.2 degrees, suggesting that the lava flow originated at a minimum elevation of 1400 m. Hence, high paleotopography in the Kenya Rift region must have existed by at least 13.5 Ma. We infer from this that middle Miocene uplift occurred, which coincides with the two-step expansion of grasslands, as well as important radiation and speciation events in tropical Africa.}, language = {en} } @article{WasiolkaJeltschHenscheletal.2010, author = {Wasiolka, Bernd and Jeltsch, Florian and Henschel, Joh and Blaum, Niels}, title = {Space use of the spotted sand lizard (Pedioplanis l. lineoocellata) under different degradation states}, issn = {0141-6707}, doi = {10.1111/j.1365-2028.2009.01085.x}, year = {2010}, abstract = {Although the effects of grazing-induced savannah degradation on animal diversity are well documented, knowledge of how they affect space use or responding behaviour remains poor. In this study, we analysed space use of the spotted sand lizard (Pedioplanis l. lineoocellata) in degraded versus nondegraded habitats of southern Kalahari savannah habitats. Lizards were radio tracked, daily movement distances recorded and home range sizes calculated. In degraded Kalahari savannah habitats where plant diversity and perennial grass cover are low but shrub cover high, P. lineoocellata moves larger distances (40.88 +/- 6.42 m versus 27.43 +/- 5.08 m) and occupies larger home ranges (646.64 +/- 244.84 m(2) versus 209.15 +/- 109.84 m(2)) than in nondegraded habitats (high plant diversity, high perennial grass cover and low shrub cover). We assume that this increase in daily movement distances and home range sizes is a behavioural plasticity to limited food resources in degraded savannah habitats. Although P. lineoocellata is able to adjust to resource-poor savannah habitats, the increase in the lizard's movement activities is likely to result in a higher predation risk. This is supported by the lower availability of protective vegetation i.e. perennial grass cover. Hence, we conclude that despite behavioural plasticity of P. lineoocellata, overgrazing has a severe negative impact on the space use of P. lineoocellata.}, language = {en} } @article{VasquezAltenbergerRomeretal.2010, author = {Vasquez, M{\´o}nica and Altenberger, Uwe and Romer, Rolf L. and Sudo, Masafumi and Moreno-Murillo, Juan Manuel}, title = {Magmatic evolution of the Andean Eastern Cordillera of Colombia during the Cretaceous : Influence of previous tectonic processes}, issn = {0895-9811}, doi = {10.1016/j.jsames.2009.02.003}, year = {2010}, abstract = {The Eastern Cordillera of the Colombian Andes represents an inverted Cretaceous basin where Cretaceous magmatism is characterized by rare mafic dykes and sills. We use Ar-40/Ar-39, Sr-Nd-Pb isotopes, as well as major and trace elements analyses of Cretaceous intrusions from both flanks of the Eastern Cordillera in combination with structural data to document the complex evolution of the basin. Magmatism, which is diachronous and geochemically diverse, seems to be related to mantle melting beneath the most subsiding segments of each sub-basin during enhanced extensional tectonics. The mafic intrusions display two different compositional series: an alkaline one with OIB-like pattern and a tholeiitic one with MORB-like features. This indicates at least two diverse mantle sources. Trace-element patterns suggest that the intrusions were emplaced in an extensional setting. Ar-40/Ar-39 dating on primary plagioclase and hornblende provides plateau ages between similar to 136 and similar to 74 Ma. The geochemical and temporal diversities show that the emplacement of the magmas was tectonically controlled, each sub-basin reflecting an individual subsidence event.}, language = {en} } @article{TuominenEerikaeinenSchibalskietal.2010, author = {Tuominen, Sakari and Eerik{\"a}inen, Kalle and Schibalski, Anett and Haakana, Markus and Lehtonen, Aleksi}, title = {Mapping Biomass Variables with a Multi-Source Forest Inventory Technique}, issn = {0037-5330}, year = {2010}, abstract = {Map form information on forest biomass is required for estimating bioenergy potentials and monitoring carbon stocks. In Finland, the growing stock of forests is monitored using multi-source forest inventory, where variables are estimated in the form of thematic maps and area statistics by combining information of field measurements, satellite images and other digital map data. In this study, we used the multi-source forest inventory methodology for estimating forest biomass characteristics. The biomass variables were estimated for national forest inventory field plots on the basis of measured tree variables. The plot-level biomass estimates were used as reference data for satellite image interpretation. The estimates produced by satellite image interpretation were tested by cross-validation. The results indicate that the method for producing biomass maps on the basis of biomass models and satellite image interpretation is operationally feasible. Furthermore, the accuracy of the estimates of biomass variables is similar or even higher than that of traditional growing stock volume estimates. The technique presented here can be applied, for example, in estimating biomass resources or in the inventory of greenhouse gases.}, language = {en} } @article{TreydteRiginosJeltsch2010, author = {Treydte, Anna C. and Riginos, Corinna and Jeltsch, Florian}, title = {Enhanced use of beneath-canopy vegetation by grazing ungulates in African savannahs}, issn = {0140-1963}, doi = {10.1016/j.jaridenv.2010.07.003}, year = {2010}, abstract = {The cover of large trees in African savannahs is rapidly declining, mainly due to human land-use practices. Trees improve grass nutrient quality and contribute to species and structural diversity of savannah vegetation. However, the response of herbivores to trees as habitat features is unknown We quantified the habitat use of wild and domestic ungulates in two eastern and southern African savannahs. We assessed grazing intensities and quantified dung depositions beneath and around canopies of different sized trees. Grasses were eaten and dung was deposited twice as frequently beneath large (ca. 5 m in height) and very large trees (7-10 m) than in open grasslands. Small trees (<2.5 m) did not show this trend. Grazing intensity and dung deposition decreased with distance away from trees at both study sites. These results suggest that large trees represent essential habitat features for domestic and wild herbivores. Increased dung depositions beneath large trees may further promote the maintenance of a patchy nutrient distribution in savannahs. Small trees cannot provide the same structural and functional advantages as large trees do. We recommend that land-use practices be promoted which conserve large single-standing trees to benefit the flora and fauna of African savannahs.}, language = {en} } @article{TrauthMaslinDeinoetal.2010, author = {Trauth, Martin H. and Maslin, Mark A. and Deino, Alan L. and Junginger, Annett and Lesoloyia, Moses and Odada, Eric O. and Olago, Daniel O. and Olaka, Lydia A. and Strecker, Manfred and Tiedemann, Ralph}, title = {Human evolution in a variable environment : the amplifier lakes of Eastern Africa}, issn = {0277-3791}, doi = {10.1016/j.quascirev.2010.07.007}, year = {2010}, abstract = {The development of rise Cenozoic East African Rift System (EARS) profoundly re-shaped the landscape and significantly increased the amplitude of short-term environmental response to climate variation. In particular, the development of amplifier lakes in rift basins after three million years ago significantly contributed to this exceptional sensitivity of East Africa to climate change compared to elsewhere on the African continent. Amplifier lakes are characterized by tectonically-formed graben morphologies in combination with an extreme contrast between high precipitation in the elevated parts of the catchment and high evaporation in the lake area. Such amplifier lakes respond rapidly to moderate, precessional-forced climate shifts, and as they do so apply dramatic environmental pressure to the biosphere. Rift basins, when either extremely dry or lake-filled, form important barriers for migration, mixing and competition of different populations of animals and hominins. Amplifier lakes link long-term, high-amplitude tectonic processes and short-term environmental fluctuations. East Africa may have become the place where early humans evolved as a consequence of this strong link between different time scales. (C) 2010 Elsevier Ltd. All rights reserved.}, language = {en} } @article{TomasZitzmannHomannetal.2010, author = {Tomas, Sara and Zitzmann, Max and Homann, Martin and Rumpf, Michael and Amour, Frederic and Benisek, Merle-Friederike and Betzler, Christian and Mutti, Maria}, title = {From ramp to platform : building a 3D model of depositional geometries and facies architectures in transitional carbonates in the Miocene, northern Sardinia}, issn = {0172-9179}, doi = {10.1007/s10347-009-0203-7}, year = {2010}, abstract = {The depositional geometry and facies distribution of an Early Miocene (Burdigalian) carbonate system in the Perfugas Basin (NW Sardinia) comprise a well-exposed example of a transition from a ramp to a steep-flanked platform. The carbonate succession (Sedini Limestone Unit) is composed of two depositional sequences separated by a major erosional unconformity. The lower (sequence 1) records a ramp dominated by heterozoan producers and the upper (sequence 2) is dominated by photozoan producers and displays a gradual steepening of the depositional profile into a steep- flanked platform. This paper shows the process of creating a digital outcrop model including a facies model. This process consists of combining field data sets, including 17 sedimentary logs, and a spatial dataset consisting of differential global positioning system data points measured along key stratigraphic surfaces and sedimentary logs, with the goal of locking traditional field observations into a 3D spatial model. Establishing a precise geometrical framework and visualizing the overall change in the platform geometry and the related vertical and lateral facies variations of the Sedini carbonate platform, allows us to better understand the sedimentary processes leading to the geometrical turn- over of the platform. Furthermore, a detailed facies modeling helps us to gain insight into the detailed depositional dynamics. The final model reproduces faithfully the depositional geometries observed in the outcrops and helps in understanding the relationships between facies and architectural framework at the basin scale. Moreover, it provides the basis to characterize semiquantitatively regional sedimentological features and to make further reservoir and subsurface analogue studies.}, language = {en} } @article{TietjenJeltschZeheetal.2010, author = {Tietjen, Britta and Jeltsch, Florian and Zehe, Erwin and Classen, Nikolaus and Groengroeft, Alexander and Schiffers, Katja and Oldeland, Jens}, title = {Effects of climate change on the coupled dynamics of water and vegetation in drylands}, issn = {1936-0584}, doi = {10.1002/Eco.70}, year = {2010}, abstract = {Drylands worldwide are exposed to a highly variable environment and face a high risk of degradation. The effects of global climate change such as altered precipitation patterns and increased temperature leading to reduced water availability will likely increase this risk. At the same time, an elevated atmospheric CO2 level could mitigate the effects of reduced water availability by increasing the water use efficiency of plants. To prevent degradation of drylands, it is essential to understand the underlying processes that affect water availability and vegetation cover. Since water and vegetation are strongly interdependent in water-limited ecosystems, changes can lead to highly non- linear effects. We assess these effects by developing an ecohydrological model of soil moisture and vegetation cover. The water component of the model simulates the daily dynamics of surface water and water contents in two soil layers. Vegetation is represented by two functional types: shrubs and grasses. These compete for soil water and strongly influence hydrological processes. We apply the model to a Namibian thornbush savanna and evaluate the separate and combined effects of decreased annual precipitation, increased temperature, more variable precipitation and elevated atmospheric CO2 on soil moisture and on vegetation cover. The results show that two main factors control the response of plant types towards climate change, namely a change in water availability and a change in water allocation to a specific plant type. Especially, reduced competitiveness of grasses can lead to a higher risk of shrub encroachment in these systems.}, language = {en} } @article{TannebergerFladePreiksaetal.2010, author = {Tanneberger, Franziska and Flade, Martin and Preiksa, Zydrunas and Schr{\"o}der-Esselbach, Boris}, title = {Habitat selection of the globally threatened Aquatic Warbler Acrocephalus paludicola at the western margin of its breeding range and implications for management}, issn = {0019-1019}, doi = {10.1111/j.1474-919X.2010.01016.x}, year = {2010}, abstract = {The globally threatened Aquatic Warbler Acrocephalus paludicola is an umbrella species for fen mires and is at risk of extinction in its westernmost breeding population due to severe habitat loss. We used boosted regression trees to model Aquatic Warbler habitat selection in order to make recommendations for effective management of the last remnant habitats. Habitat data were collected in the years 2004-2006 in all remaining breeding sites in Pomerania (eastern Germany and western Poland) as well as in recently abandoned sites. Models were validated using data from similar Aquatic Warbler habitats in Lithuania. The probability of occurrence of Aquatic Warblers in late May/early June was positively associated with low isolation from other occupied sites, less eutrophic conditions, a high proportion of area mown early in the preceding year, high availability of vegetation 60-70 cm high, high prey abundance and high habitat heterogeneity. Early summer land management is needed in the more productive sites to prevent habitat deterioration by succession to higher and denser vegetation. As this also poses a serious threat to broods, management that creates a mosaic of early and late used patches is recommended to preserve and restore productive Aquatic Warbler sites. In less productive sites, winter mowing can maintain suitable habitat conditions. Aquatic Warbler-friendly land use supports a variety of other threatened plant and animal species typical of fens and sedge meadows and can meet the economic interests of local land users.}, language = {en} } @phdthesis{Strollo2010, author = {Strollo, Angelo}, title = {Development of techniques for earthquake microzonation studies in different urban environment}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-53807}, school = {Universit{\"a}t Potsdam}, year = {2010}, abstract = {The proliferation of megacities in many developing countries, and their location in areas where they are exposed to a high risk from large earthquakes, coupled with a lack of preparation, demonstrates the requirement for improved capabilities in hazard assessment, as well as the rapid adjustment and development of land-use planning. In particular, within the context of seismic hazard assessment, the evaluation of local site effects and their influence on the spatial distribution of ground shaking generated by an earthquake plays an important role. It follows that the carrying out of earthquake microzonation studies, which aim at identify areas within the urban environment that are expected to respond in a similar way to a seismic event, are essential to the reliable risk assessment of large urban areas. Considering the rate at which many large towns in developing countries that are prone to large earthquakes are growing, their seismic microzonation has become mandatory. Such activities are challenging and techniques suitable for identifying site effects within such contexts are needed. In this dissertation, I develop techniques for investigating large-scale urban environments that aim at being non-invasive, cost-effective and quickly deployable. These peculiarities allow one to investigate large areas over a relative short time frame, with a spatial sampling resolution sufficient to provide reliable microzonation. Although there is a negative trade-off between the completeness of available information and extent of the investigated area, I attempt to mitigate this limitation by combining two, what I term layers, of information: in the first layer, the site effects at a few calibration points are well constrained by analyzing earthquake data or using other geophysical information (e.g., shear-wave velocity profiles); in the second layer, the site effects over a larger areal coverage are estimated by means of single-station noise measurements. The microzonation is performed in terms of problem-dependent quantities, by considering a proxy suitable to link information from the first layer to the second one. In order to define the microzonation approach proposed in this work, different methods for estimating site effects have been combined and tested in Potenza (Italy), where a considerable amount of data was available. In particular, the horizontal-to-vertical spectral ratio computed for seismic noise recorded at different sites has been used as a proxy to combine the two levels of information together and to create a microzonation map in terms of spectral intensity ratio (SIR). In the next step, I applied this two-layer approach to Istanbul (Turkey) and Bishkek (Kyrgyzstan). A similar hybrid approach, i.e., combining earthquake and noise data, has been used for the microzonation of these two different urban environments. For both cities, after having calibrated the fundamental frequencies of resonance estimated from seismic noise with those obtained by analysing earthquakes (first layer), a fundamental frequency map has been computed using the noise measurements carried out within the town (second layer). By applying this new approach, maps of the fundamental frequency of resonance for Istanbul and Bishkek have been published for the first time. In parallel, a microzonation map in terms of SIR has been incorporated into a risk scenario for the Potenza test site by means of a dedicated regression between spectral intensity (SI) and macroseismic intensity (EMS). The scenario study confirms the importance of site effects within the risk chain. In fact, their introduction into the scenario led to an increase of about 50\% in estimates of the number of buildings that would be partially or totally collapsed. Last, but not least, considering that the approach developed and applied in this work is based on measurements of seismic noise, their reliability has been assessed. A theoretical model describing the self-noise curves of different instruments usually adopted in microzonation studies (e.g., those used in Potenza, Istanbul and Bishkek) have been considered and compared with empirical data recorded in Cologne (Germany) and Gubbio (Italy). The results show that, depending on the geological and environmental conditions, the instrumental noise could severely bias the results obtained by recording and analysing ambient noise. Therefore, in this work I also provide some guidelines for measuring seismic noise.}, language = {en} } @article{StreichBeckenRitter2010, author = {Streich, Rita and Becken, Michael and Ritter, Oliver}, title = {Imaging of CO2 storage sites, geothermal reservoirs, and gas shales using controlled-source magnetotellurics : modeling studies}, issn = {0009-2819}, doi = {10.1016/j.chemer.2010.05.004}, year = {2010}, abstract = {To balance the steady decrease of conventional hydrocarbon resources, increased utilization of unconventional and new energy resources, such as shale gas and geothermal energy, is required. Also, the geological sequestration of carbon dioxide is being considered as a technology that may temporarily mitigate the effects of CO2 emission. Sites suitable for shale gas production, geothermal exploration, or CO2 sequestration are commonly characterized by electrical resistivities distinctly different from those of the surrounding rocks. Therefore, electromagnetic methods can be viable tools to help identify target sites suitable for exploration, and to monitor reservoirs during energy production or CO2 injection. Among the wide variety of electromagnetic methods available, controlled-source magnetotelluric (CSMT) may be particularly suitable because of (i) its ability to resolve both electrically resistive and conductive structures, (ii) controlled sources offering noise control and thus facilitating surveys in populated regions, and (iii) the potential of penetration throughout the depth range accessible by drilling. Nevertheless, CSMT has not yet been widely employed because of logistical challenges of field operations and the requirement of complex and highly computer-intensive data processing. With these difficulties gradually being mitigated by recent technological developments, CSMT may now be reconsidered as an exploration tool. Here, we investigate by 1D and 3D numerical simulations the feasibility of detecting gas shales and identifying sites eligible for geothermal exploration or CO2 sequestration from CSMT data. We consider surface-to-surface, borehole-to-surface, and cross-hole configurations of the sources and receivers. Results and conclusions on the detectability of the targets of interest are presented for various exploration and monitoring scenarios, which are roughly representative of the geological setting of the North German Basin.}, language = {en} } @article{StephTiedemannPrangeetal.2010, author = {Steph, Silke and Tiedemann, Ralph and Prange, Matthias and Groeneveld, J{\"u}rgen and Schulz, Michael Karl and Timmermann, Axel and N{\"u}rnberg, Dirk and R{\"u}hlemann, Carsten and Saukel, Cornelia and Haug, Gerald H.}, title = {Early Pliocene increase in thermohaline overturning : a precondition for the development of the modern equatorial Pacific cold tongue}, issn = {0883-8305}, doi = {10.1029/2008pa001645}, year = {2010}, abstract = {Unraveling the processes responsible for Earth's climate transition from an "El Nino-like state" during the warm early Pliocene into a modern-like "La Nina-dominated state" currently challenges the scientific community. Recently, the Pliocene climate switch has been linked to oceanic thermocline shoaling at similar to 3 million years ago along with Earth's final transition into a bipolar icehouse world. Here we present Pliocene proxy data and climate model results, which suggest an earlier timing of the Pliocene climate switch and a different chain of forcing mechanisms. We show that the increase in North Atlantic meridional overturning circulation between 4.8 and 4.0 million years ago, initiated by the progressive closure of the Central American Seaway, triggered overall shoaling of the tropical thermocline. This preconditioned the turnaround from a warm eastern equatorial Pacific to the modern equatorial cold tongue state about 1 million years earlier than previously assumed. Since similar to 3.6-3.5 million years ago, the intensification of Northern Hemisphere glaciation resulted in a strengthening of the trade winds, thereby amplifying upwelling and biogenic productivity at low latitudes.}, language = {en} } @article{SobelSeward2010, author = {Sobel, Edward and Seward, Diane}, title = {Influence of etching conditions on apatite fission-track etch pit diameter}, issn = {0009-2541}, doi = {10.1016/j.chemgeo.2009.12.012}, year = {2010}, abstract = {The temperature range of the apatite fission-track partial annealing zone (PAZ) varies as a function of the kinetic characteristics of the apatite crystal as well as the cooling rate. These kinetic characteristics are controlled largely by the chemical composition of the apatites; fission-track etch pit diameter (D-par) has been shown to be a proxy for quantifying these characteristics. Since some annealing models explicitly use D-par as a kinetic indicator, the question of whether D-par varies according to either operator variability or etching conditions has serious implications. The influence of etching conditions includes both minor variations in temperature and etch time as well as larger variations between the three popular concentrations of nitric acid that are presently in use within the fission- track community: 5 and 5.5 M HNO3 at 21 degrees C for 20 s, and 1.6 M (7 vol.\%) HNO3 at 21 degrees C for 45 or 50 s. We have conducted a systematic study of D-par in two widely used apatite age standards: Durango and Fish Canyon Tuff (FCT). Samples were analyzed by at least 2 operators in 2 laboratories; over 15,000 etch pits were measured. Operator imprecision is small compared to other effects both within and between operators. The measured D-par parameters depend primarily on etching conditions including concentration of etchant, and the time and temperature of etching. D-par size increases linearly with time but non-linearly with temperature. Differences between D-par size of Durango and FCT are significantly smaller for the 1.6 M etch compared to the other two etches, implying that the weak etch has less resolving power for D-par measurements. In addition, when etching conditions are changed, confined track lengths for different apatite compositions and D-par values co-vary in a complex fashion. Therefore, the 5 or 5.5 M etch is recommended for studies that employ D-par as a proxy for chemistry; the 1.6 M etch and other weak etches are not recommended. We propose a linear correction for D-par based on measuring D-par in two widely available AFT age standards, Durango and FCT, and cross-plotting the results with those obtained by R. Donelick. The slope of the resulting curve, which passes through the origin, provides a correction factor. This system yields reasonably good corrected values for 3 test samples using the 5 M and 5.5 M HNO3 etches.}, language = {en} } @article{SobelSeward2010, author = {Sobel, Edward and Seward, Diane}, title = {Erratum to "Influence of etching conditions on apatite fission-track etch pit diameter by E.R. Sobel and D. Seward" [Chem. Geol. 271 (2010) 59-69]}, issn = {0009-2541}, doi = {10.1016/j.chemgeo.2010.08.005}, year = {2010}, language = {en} } @article{SimpsonVanMeirvenneLucketal.2010, author = {Simpson, David and Van Meirvenne, Marc and Luck, Erika and Ruhlmann, J{\"o}rg and Saey, Timothy and Bourgeois, Jean}, title = {Sensitivity of multi-coil frequency domain electromagnetic induction sensors to map soil magnetic susceptibility}, issn = {1351-0754}, doi = {10.1111/j.1365-2389.2010.01261.x}, year = {2010}, abstract = {Magnetic susceptibility is an important indicator of anthropogenic disturbance in the natural soil. This property is often mapped with magnetic gradiometers in archaeological prospection studies. It is also detected with frequency domain electromagnetic induction (FDEM) sensors, which have the advantage that they can simultaneously measure the electrical conductivity. The detection level of FDEM sensors for magnetic structures is very dependent on the coil configuration. Apart from theoretical modelling studies, a thorough investigation with field models has not been conducted until now. Therefore, the goal of this study was to test multiple coil configurations on a test field with naturally enhanced magnetic susceptibility in the topsoil and with different types of structures mimicking real archaeological features. Two FDEM sensors were used with coil separations between 0.5 and 2 m and with three coil orientations. First, a vertical sounding was conducted over the undisturbed soil to test the validity of a theoretical layered model, which can be used to infer the depth sensitivity of the coil configurations. The modelled sounding values corresponded well with the measured data, which means that the theoretical models are applicable to layered soils. Second, magnetic structures were buried in the site and the resulting anomalies measured to a very high resolution. The results showed remarkable differences in amplitude and complexity between the responses of the coil configurations. The 2-m horizontal coplanar and 1.1-m perpendicular coil configurations produced the clearest anomalies and resembled best a gradiometer measurement.}, language = {en} } @article{SimpsonVanMeirvenneLucketal.2010, author = {Simpson, David and Van Meirvenne, Marc and Luck, Erika and Bourgeois, Jean and Ruhlmann, J{\"o}rg}, title = {Prospection of two circular Bronze Age ditches with multi-receiver electrical conductivity sensors (North Belgium)}, issn = {0305-4403}, doi = {10.1016/j.jas.2010.03.017}, year = {2010}, abstract = {Two types of electrical conductivity sensors were evaluated to prospect circular ditches surrounding former Bronze Age burial mounds, complementing aerial photography. The first sensor was based on the electrical resistivity (ER) method, while the second sensor was based on frequency-domain electromagnetic induction (FDEM). Both sensors were designed with multiple receivers, which measure several depth sensitivities simultaneously. First, the sensors were tested on an experimental site where a rectangular structure with limited dimensions was dug in a sandy soil. The structure appeared as a higher conductivity anomaly in the low-conductivity sand. Then, both methods were applied on two Bronze Age sites with different soil properties, which were discovered by aerial photography. The first site, in a sandy soil, gave only very weak anomalies. Soil augering revealed that the ditch filling consisted of the same sandy material as the surrounding, therefore this filling was not able to cause a high-conductivity contrast. Due to its lower sensitivity to noise in the low-conductive range, the ER-sensor produced a more pronounced anomaly than the FDEM-sensor. The second site was located on top of a ridge with a shallow substrate of Tertiary, coastal sediments. The ditch was very clearly visible on the sensor maps as a conductive low. At this location, the soil augering revealed that the ditch was dug through an alternating clay-sand layer and subsequently filled up with silty material from the topsoil. Overall, the shallow receiver separation produced anomalies that were both stronger and that corresponded better to the geometry of the ditches. The other receiver separations provided more information on the natural soil layering, and in the case of the ER-array they could be used to obtain a cross-section of the actual electrical conductivity with 2-D inversion modelling. The results of this study proofed that conductivity sensors can detect Bronze Age ditches, with varying contrast depending on the soil geomorphology. Moreover, the sensor maps combined with soil observations by coring provided insight in the environmental conditions that influence the contrast of the anomalies seen on the aerial photographs and the sensor maps.}, language = {en} } @article{SiliverstovsOetschKemfertetal.2010, author = {Siliverstovs, Boriss and {\"O}tsch, Rainald and Kemfert, Claudia and Jaeger, Carlo and Haas, Armin and Kremers, Hans}, title = {Climate change and modelling of extreme temperatures in Switzerland}, issn = {1436-3240}, doi = {10.1007/s00477-009-0321-3}, year = {2010}, abstract = {This study models maximum temperatures in Switzerland monitored in twelve locations using the generalised extreme value (GEV) distribution. The parameters of the GEV distribution are determined within a Bayesian framework. We find that the parameters of the underlying distribution underwent a substantial change in the beginning of the 1980s. This change is characterised by an increase both in the level and the variability. We assess the likelihood of the heat wave of the summer 2003 using the fitted GEV distribution by accounting for the presence of a structural break. The estimation results do suggest that the heat wave of 2003 is not that statistically improbable if an appropriate methodology is used for dealing with nonstationarity.}, language = {en} } @article{SigmanHainHaug2010, author = {Sigman, Daniel M. and Hain, Mathis P. and Haug, Gerald H.}, title = {The polar ocean and glacial cycles in atmospheric CO2 concentration}, issn = {0028-0836}, doi = {10.1038/Nature09149}, year = {2010}, abstract = {Global climate and the atmospheric partial pressure of carbon dioxide (p(CO2atm)) are correlated over recent glacial cycles, with lower p(CO2atm) during ice ages, but the causes of the p(CO2atm) changes are unknown. The modern Southern Ocean releases deeply sequestered CO2 to the atmosphere. Growing evidence suggests that the Southern Ocean CO2 'leak' was stemmed during ice ages, increasing ocean CO2 storage. Such a change would also have made the global ocean more alkaline, driving additional ocean CO2 uptake. This explanation for lower ice-age p(CO2atm), if correct, has much to teach us about the controls on current ocean processes.}, language = {en} } @phdthesis{Shirzaei2010, author = {Shirzaei, Manoochehr}, title = {Crustal deformation source monitoring using advanced InSAR time series and time dependent inverse modeling}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-50774}, school = {Universit{\"a}t Potsdam}, year = {2010}, abstract = {Crustal deformation can be the result of volcanic and tectonic activity such as fault dislocation and magma intrusion. The crustal deformation may precede and/or succeed the earthquake occurrence and eruption. Mitigating the associated hazard, continuous monitoring of the crustal deformation accordingly has become an important task for geo-observatories and fast response systems. Due to highly non-linear behavior of the crustal deformation fields in time and space, which are not always measurable using conventional geodetic methods (e.g., Leveling), innovative techniques of monitoring and analysis are required. In this thesis I describe novel methods to improve the ability for precise and accurate mapping the spatiotemporal surface deformation field using multi acquisitions of satellite radar data. Furthermore, to better understand the source of such spatiotemporal deformation fields, I present novel static and time dependent model inversion approaches. Almost any interferograms include areas where the signal decorrelates and is distorted by atmospheric delay. In this thesis I detail new analysis methods to reduce the limitations of conventional InSAR, by combining the benefits of advanced InSAR methods such as the permanent scatterer InSAR (PSI) and the small baseline subsets (SBAS) with a wavelet based data filtering scheme. This novel InSAR time series methodology is applied, for instance, to monitor the non-linear deformation processes at Hawaii Island. The radar phase change at Hawaii is found to be due to intrusions, eruptions, earthquakes and flank movement processes and superimposed by significant environmental artifacts (e.g., atmospheric). The deformation field, I obtained using the new InSAR analysis method, is in good agreement with continuous GPS data. This provides an accurate spatiotemporal deformation field at Hawaii, which allows time dependent source modeling. Conventional source modeling methods usually deal with static deformation field, while retrieving the dynamics of the source requires more sophisticated time dependent optimization approaches. This problem I address by combining Monte Carlo based optimization approaches with a Kalman Filter, which provides the model parameters of the deformation source consistent in time. I found there are numerous deformation sources at Hawaii Island which are spatiotemporally interacting, such as volcano inflation is associated to changes in the rifting behavior, and temporally linked to silent earthquakes. I applied these new methods to other tectonic and volcanic terrains, most of which revealing the importance of associated or coupled deformation sources. The findings are 1) the relation between deep and shallow hydrothermal and magmatic sources underneath the Campi Flegrei volcano, 2) gravity-driven deformation at Damavand volcano, 3) fault interaction associated with the 2010 Haiti earthquake, 4) independent block wise flank motion at the Hilina Fault system, Kilauea, and 5) interaction between salt diapir and the 2005 Qeshm earthquake in southern Iran. This thesis, written in cumulative form including 9 manuscripts published or under review in peer reviewed journals, improves the techniques for InSAR time series analysis and source modeling and shows the mutual dependence between adjacent deformation sources. These findings allow more realistic estimation of the hazard associated with complex volcanic and tectonic systems.}, language = {en} } @article{SeglGuanterKaufmannetal.2010, author = {Segl, Karl and Guanter, Luis and Kaufmann, Hermann and Schubert, Josef and Kaiser, Stefan and Sang, Bernhard and Hofer, Stefan}, title = {Simulation of spatial sensor characteristics in the context of the EnMAP Hyperspectral mission}, issn = {0196-2892}, doi = {10.1109/Tgrs.2010.2042455}, year = {2010}, abstract = {The simulation of remote sensing images is a valuable tool for defining future Earth observation systems, optimizing instrument parameters, and developing and validating data-processing algorithms. A scene simulator for optical Earth observation data has been developed within the Environmental Mapping and Analysis Program (EnMAP) hyperspectral mission. It produces EnMAP-like data following a sequential processing approach consisting of five independent modules referred to as reflectance, atmospheric, spatial, spectral, and radiometric modules. From a modeling viewpoint, the spatial module is the most complex. The spatial simulation process considers the satellite-target geometry, which is adapted to the EnMAP orbit and operating characteristics, the instrument spatial response, and the sources of spatial nonuniformity (keystone, telescope distortion and smile, and detector coregistration). The spatial module of the EnMAP scene simulator is presented in this paper. The EnMAP spatial and geometric characteristics will be described, the simulation methodology will be presented in detail, and the capability of the EnMAP simulator will be shown by illustrative examples.}, language = {en} } @phdthesis{Scherler2010, author = {Scherler, Dirk}, title = {Climate variability and glacial dynamics in the Himalaya}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-49871}, school = {Universit{\"a}t Potsdam}, year = {2010}, abstract = {In den Hochgebirgen Asiens bedecken Gletscher eine Fl{\"a}che von ungef{\"a}hr 115,000 km² und ergeben damit, neben Gr{\"o}nland und der Antarktis, eine der gr{\"o}ßten Eisakkumulationen der Erde. Die Sensibilit{\"a}t der Gletscher gegen{\"u}ber Klimaschwankungen macht sie zu wertvollen pal{\"a}oklimatischen Archiven in Hochgebirgen, aber gleichzeitig auch anf{\"a}llig gegen{\"u}ber rezenter und zuk{\"u}nftiger globaler Erw{\"a}rmung. Dies kann vor allem in dicht besiedelten Gebieten S{\"u}d-, Ost- und Zentralasiens zu großen Problem f{\"u}hren, in denen Gletscher- und Schnee-Schmelzw{\"a}sser eine wichtige Ressource f{\"u}r Landwirtschaft und Stromerzeugung darstellen. Eine erfolgreiche Prognose des Gletscherverhaltens in Reaktion auf den Klimawandel und die Minderung der sozio{\"o}konomischen Auswirkungen erfordert fundierte Kenntnisse der klimatischen Steuerungsfaktoren und der Dynamik asiatischer Gletscher. Aufgrund ihrer Abgeschiedenheit und dem erschwerten Zugang gibt es nur wenige glaziologische Gel{\"a}ndestudien, die zudem r{\"a}umlich und zeitlich sehr begrenzt sind. Daher fehlen bisher grundlegende Informationen {\"u}ber die Mehrzahl asiatischer Gletscher. In dieser Arbeit benutze ich verschiedene Methoden, um die Dynamik asiatischer Gletscher auf mehreren Zeitskalen zu untersuchen. Erstens teste ich eine Methode zur pr{\"a}zisen satelliten-gest{\"u}tzten Messung von Gletscheroberfl{\"a}chen-Geschwindigkeiten. Darauf aufbauend habe ich eine umfassende regionale Erhebung der Fliessgeschwindigkeiten und Frontdynamik asiatischer Gletscher f{\"u}r die Jahre 2000 bis 2008 durchgef{\"u}hrt. Der gewonnene Datensatz erlaubt einmalige Einblicke in die topographischen und klimatischen Steuerungsfaktoren der Gletscherfließgeschwindigkeiten in den Gebirgsregionen Hochasiens. Insbesondere dokumentieren die Daten rezent ungleiches Verhalten der Gletscher im Karakorum und im Himalaja, welches ich auf die konkurrierenden klimatischen Einfl{\"u}sse der Westwinddrift im Winter und des Indischen Monsuns im Sommer zur{\"u}ckf{\"u}hre. Zweitens untersuche ich, ob klimatisch bedingte Ost-West Unterschiede im Gletscherverhalten auch auf l{\"a}ngeren Zeitskalen eine Rolle spielen und gegebenenfalls f{\"u}r dokumentierte regional asynchrone Gletschervorst{\"o}ße relevant sind. Dazu habe ich mittels kosmogener Nuklide Oberfl{\"a}chenalter von erratischen Bl{\"o}cken auf Mor{\"a}nen ermittelt und eine glaziale Chronologie f{\"u}r das obere Tons Tal, in den Quellgebieten des Ganges, erstellt. Dieses Gebiet befindet sich in der {\"U}bergangszone von monsunaler zu Westwind beeinflusster Feuchtigkeitszufuhr und ist damit ideal gelegen, um die Auswirkungen dieser beiden atmosph{\"a}rischen Zirkulationssysteme auf Gletschervorst{\"o}ße zu untersuchen. Die ermittelte glaziale Chronologie dokumentiert mehrere Gletscherschwankungen w{\"a}hrend des Endstadiums der letzten Pleistoz{\"a}nen Vereisung und w{\"a}hrend des Holz{\"a}ns. Diese weisen darauf hin, dass Gletscherschwankungen im westlichen Himalaja weitestgehend synchron waren und auf graduelle glaziale-interglaziale Temperaturver{\"a}nderungen, {\"u}berlagert von monsunalen Niederschlagsschwankungen h{\"o}herer Frequenz, zur{\"u}ck zu f{\"u}hren sind. In einem dritten Schritt kombiniere ich Satelliten-Klimadaten mit Eisfluss-Absch{\"a}tzungen und topographischen Analysen, um den Einfluss der Gletscher Hochasiens auf die Reliefentwicklung im Hochgebirge zu untersuchen. Die Ergebnisse dokumentieren ausgepr{\"a}gte meridionale Unterschiede im Grad und im Stil der Vergletscherung und glazialen Erosion in Abh{\"a}ngigkeit von topographischen und klimatischen Faktoren. Gegens{\"a}tzlich zu bisherigen Annahmen deuten die Daten darauf hin, dass das monsunale Klima im zentralen Himalaja die glaziale Erosion schw{\"a}cht und durch den Erhalt einer steilen orographischen Barriere das Tibet Plateau vor lateraler Zerschneidung bewahrt. Die Ergebnisse dieser Arbeit dokumentieren, wie klimatische und topographische Gradienten die Gletscherdynamik in den Hochgebirgen Asiens auf Zeitskalen von 10^0 bis 10^6 Jahren beeinflussen. Die Reaktionszeit der Gletscher auf Klimaver{\"a}nderungen sind eng an Eigenschaften wie Schuttbedeckung und Neigung gekoppelt, welche ihrerseits von den topographischen Verh{\"a}ltnissen bedingt sind. Derartige Einflussfaktoren m{\"u}ssen bei pal{\"a}oklimatischen Rekonstruktion und Vorhersagen {\"u}ber die Entwicklung asiatischer Gletscher ber{\"u}cksichtigt werden. Desweiteren gehen die regionalen topographischen Unterschiede der vergletscherten Gebiete Asiens teilweise auf klimatische Gradienten und den langfristigen Einfluss der Gletscher auf die topographische Entwicklung des Gebirgssystems zur{\"u}ck.}, language = {en} } @article{SchaefliHuss2010, author = {Schaefli, Bettina and Huss, Matthias}, title = {Simulation of high mountainous discharge : how much information do we need?}, issn = {1812-2108}, doi = {10.5194/hessd-7-8661-2010}, year = {2010}, abstract = {The hydrologic cycle of high mountainous catchments is frequently simulated with simple precipitation-discharge models representing the snow accumulation and ablation behavior of a very complex environment with a set of lumped equations accounting for altitudinal temperature and precipitation gradients. In this study, we present a methodology to include sparse snow depths measurements into the calibration process. Based on this methodology, we assess for a case study, the Rhonegletscher catchment (Switzerland), how much observed information we need to reliably calibrate the model, such that it reproduces the dominant system dynamics, discharge, as well as glacier mass balance. Here, we focus on the question whether observed discharge is sufficient as a calibration variable or whether we need annual or even seasonal glacier mass balance data. Introducing seasonally variable accumulation and ablation parameters is sufficient to enable the simple model to reproduce observed seasonal mass balances for the Rhonegletscher. Furthermore, our results suggest that calibrating the hydrological model exclusively on discharge can lead to wrong representations of the intra- annual accumulation and ablation processes and to a strong bias in long term glacier mass balance simulations. Adding only a few annual mass balance observations considerably reduces this bias. Calibrating exclusively on annual balance data can, in turn, lead to wrong seasonal mass balance simulations. Even if these results are case study specific, our conclusions provide valuable new insights into the benefit of different types of observations for calibrating hydrological models in glacier catchments. The presented multi-signal calibration framework and the simple method to calibrate a semi-lumped model on point observations has potential for application in other modeling contexts.}, language = {en} } @phdthesis{Ruch2010, author = {Ruch, Jo{\"e}l}, title = {Volcano deformation analysis in the Lazufre area (central Andes) using geodetic and geological observations}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-47361}, school = {Universit{\"a}t Potsdam}, year = {2010}, abstract = {Large-scale volcanic deformation recently detected by radar interferometry (InSAR) provides new information and thus new scientific challenges for understanding volcano-tectonic activity and magmatic systems. The destabilization of such a system at depth noticeably affects the surrounding environment through magma injection, ground displacement and volcanic eruptions. To determine the spatiotemporal evolution of the Lazufre volcanic area located in the central Andes, we combined short-term ground displacement acquired by InSAR with long-term geological observations. Ground displacement was first detected using InSAR in 1997. By 2008, this displacement affected 1800 km2 of the surface, an area comparable in size to the deformation observed at caldera systems. The original displacement was followed in 2000 by a second, small-scale, neighbouring deformation located on the Lastarria volcano. We performed a detailed analysis of the volcanic structures at Lazufre and found relationships with the volcano deformations observed with InSAR. We infer that these observations are both likely to be the surface expression of a long-lived magmatic system evolving at depth. It is not yet clear whether Lazufre may trigger larger unrest or volcanic eruptions; however, the second deformation detected at Lastarria and the clear increase of the large-scale deformation rate make this an area of particular interest for closer continuous monitoring.}, language = {en} } @phdthesis{Plenkers2010, author = {Plenkers, Katrin}, title = {On the Characteristics of Mining-Induced Seismicity with Magnitudes -57 every year, being more than 15\% of the events worldwide. Earthquakes with magnitude above 9 happen far more infrequently, but with catastrophic effects. The most severe consequences thereby arise from tsunamis triggered by these subduction-related earthquakes, as the Sumatra-Andaman event in 2004 showed. In order to enable efficient tsunami early warning, which includes the estimation of wave heights and arrival times, it is necessary to combine different types of real-time sensor data with numerical models of earthquake sources and tsunami propagation. This thesis was created as a result of the GITEWS project (German Indonesian Tsunami Early Warning System). It is based on five research papers and manuscripts. Main project-related task was the development of a database containing realistic earthquake scenarios for the Sunda Arc. This database provides initial conditions for tsunami propagation modeling used by the simulation system at the early warning center. An accurate discretization of the subduction geometry, consisting of 25x150 subfaults was constructed based on seismic data. Green's functions, representing the deformational response to unit dip- and strike slip at the subfaults, were computed using a layered half-space approach. Different scaling relations for earthquake dimensions and slip distribution were implemented. Another project-related task was the further development of the 'GPS-shield' concept. It consists of a constellation of near field GPS-receivers, which are shown to be very valuable for tsunami early warning. The major part of this thesis is related to the geophysical interpretation of GPS data. Coseismic surface displacements caused by the 2004 Sumatra earthquake are inverted for slip at the fault. The effect of different Earth layer models is tested, favoring continental structure. The possibility of splay faulting is considered and shown to be a secondary order effect in respect to tsunamigenity for this event. Tsunami models based on source inversions are compared to satellite radar altimetry observations. Postseismic GPS time series are used to test a wide parameter range of uni- and biviscous rheological models of the asthenosphere. Steady-state Maxwell rheology is shown to be incompatible with near-field GPS data, unless large afterslip, amounting to more than 10\% of the coseismic moment is assumed. In contrast, transient Burgers rheology is in agreement with data without the need for large aseismic afterslip. Comparison to postseismic geoid observation by the GRACE satellites reveals that even with afterslip, the model implementing Maxwell rheology results in amplitudes being too small, and thus supports a biviscous asthenosphere. A simple approach based on the assumption of quasi-static deformation propagation is introduced and proposed for inversion of coseismic near-field GPS time series. Application of this approach to observations from the 2004 Sumatra event fails to quantitatively reconstruct the rupture propagation, since a priori conditions are not fulfilled in this case. However, synthetic tests reveal the feasibility of such an approach for fast estimation of rupturing properties.}, language = {en} } @phdthesis{Hoff2010, author = {Hoff, Ulrike}, title = {Freshwater diatoms as indicators for Holocene environmental- and climate changes on Kamachatka, Russia}, address = {Potsdam}, pages = {VII, 95, XXXVIII S.}, year = {2010}, language = {en} } @article{HintersbergerThiedeStreckeretal.2010, author = {Hintersberger, Esther and Thiede, Rasmus Christoph and Strecker, Manfred and Hacker, Bradley R.}, title = {East-west extension in the NW Indian Himalaya}, issn = {0016-7606}, doi = {10.1130/B26589.1}, year = {2010}, abstract = {Explaining the presence of normal faults in overall compressive settings is a challenging problem in understanding the tectonics of active mountain belts. The Himalayan-Tibetan orogenic system is an excellent setting to approach this problem because it preserves one of the most dramatic records of long-term, contemporaneous shortening and extension. Over the past decades, several studies have described extensional features, not only in the Tibetan Plateau, but also in the Himalaya. For a long time, the favored model explained the function of the Southern Tibetan detachment system, a major fault zone in the Himalaya, as a decoupling horizon between the regime of crustal shortening forming the Himalayan wedge to the south and the extensional regime of the Tibetan Plateau to the north. However, in recent years, increasing evidence has shown that N-S-trending normal faults in the Central Himalaya crosscut not only the Southern Tibetan detachment system, but also the Main Central thrust. Here, we present new structural data and geologic evidence collected within the NW Indian Himalaya and combine them with previously published seismicity data sets in order to document pervasive E-W extension accommodated along N-S-trending faults extending as far south as the footwall of the Main Central thrust. We conducted a kinematic analysis of fault striations on brittle faults, documented and mapped fault scarps in Quaternary sedimentary deposits using satellite imagery, and made field observations in the Greater Sutlej region (Spiti, Lahul, Kinnaur) and the Garhwal Himalaya. Studies of extensional features within the regionally NW- SE-trending NW Indian Himalaya provide the advantage that arc-parallel and E-W extension can be separated, in contrast to the Central Himalaya. Therefore, our observations of E-W extension in the Indian NW Himalaya are well suited to test the applicability of current tectonic models for the whole Himalaya. We favor the interpretation of E-W extension in the NW Indian Himalaya as a propagation of extension driven by collapse of the Tibetan Plateau.}, language = {en} } @misc{HiemerRoesslerScherbaum2010, author = {Hiemer, Stefan and R{\"o}ßler, Dirk and Scherbaum, Frank}, title = {Catalog of Swarm Earthquakes in Vogtland /West Bohemia in 2008/09}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-51710}, year = {2010}, abstract = {The document contains the catalog of earthquakes in Vogtland /West Bohemia within the period of 2008/10/19 -to- 2009/03/16. The events were recorded by a seismic mini-array operated by the Institute of Earthsciences, University of Postdam.}, language = {en} } @misc{HiemerRoesslerScherbaum2010, author = {Hiemer, Stefan and R{\"o}ßler, Dirk and Scherbaum, Frank}, title = {Catalog of Swarm Earthquakes in Vogtland /West Bohemia in 2008/09}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-53837}, year = {2010}, abstract = {The document contains the catalog of earthquakes in Vogtland /West Bohemia within the period of 2008/10/19 -to- 2009/03/16. The events were recorded by a seismic mini-array operated by the Institute of Earthsciences, University of Postdam.}, language = {en} } @article{HerzschuhMischkeMeyeretal.2010, author = {Herzschuh, Ulrike and Mischke, Steffen and Meyer, Hanno and Plessen, Birgit and Zhang, Chengjun}, title = {Lake nutrient variability inferred from elemental (C, N, S) and isotopic (delta C-13, delta N-15) analyses of aquatic plant macrofossils}, issn = {0277-3791}, doi = {10.1016/j.quascirev.2010.05.011}, year = {2010}, abstract = {This paper aims to highlight the potential of using elemental and stable isotope analyses of aquatic macrophytes in palaeolimnological studies. Potamogeton pectinatus material was collected from modem plants (n=68) and from late glacial and Holocene-aged sediments from Koucha Lake (northeastern Tibetan Plateau; 34.0 degrees N; 97.2 degrees E; 4540 m a.s.l.). It was analyzed for delta C-13(Potamogeton) (modern: -23 to 0 parts per thousand, fossil: -19 to -4 parts per thousand) and delta N-15(Potamogeton) (modern: -11.0 to +13.8 parts per thousand, fossil: -9.5 to +6.7 parts per thousand) in addition to elemental carbon and nitrogen (modem C/N-Potamogeton: 7 to 29; fossil: 13 to 68) and sulfur (fossil: 188-899 mu mol/g dry weight). Fossil data were interpreted in terms of palaeo-nutrient availability and palaeo-productivity based on the modem relationships between various proxies and certain environmental data. Productivity of Potamogeton pectinatus mats at Koucha Lake as indicated by palaeo-epsilon(Potamogeton-TIC) (i.e. the enrichment of delta C-13(Potamogeton) relative to the delta(CTIC)-C-13) was reduced during periods of high conductivity, especially between 10.3 and 7.4 cal kyr BP. Potamogeton pectinatus material from these periods was also characterized by high S-Potamogeton indicating high sulfide concentrations and anoxic conditions within the sediments. However, C/N- Potamogeton ratios and delta N-15(Potamogeton) from the lower core section were found to have been altered by decompositional processes. A pronounced shift in the aquatic productivity of Lake Koucha occurred at similar to 7.4 cal kyr BP when the hydrological conditions shifted towards an open lake system and water depth increased. At this time a strong increase in productivity led to a strong decrease in the water HCO3- concentration as inferred from the application of a epsilon-(Potamogeton-TIC)-InHCO3- transfer function. A comparison of reconstructed productivity changes from Koucha Lake with further environmental proxies suggests that primary productivity changes are probably a function of internal lake dynamics and were only indirectly triggered by climate change.}, language = {en} } @article{HerzschuhMischkeMeyeretal.2010, author = {Herzschuh, Ulrike and Mischke, Steffen and Meyer, Hanno and Plessen, Birgit and Zhang, Chengjun}, title = {Using variations in the stable carbon isotope composition of macrophyte remains to quantify nutrient dynamics in lakes}, issn = {0921-2728}, doi = {10.1007/s10933-009-9365-0}, year = {2010}, abstract = {The apparent isotope enrichment factor epsilon(macrophyte) of submerged plants (epsilon(macrophyte-DIC) = delta C-13(macrophyte) - delta C-13(DIC)) is indicative of dissolved inorganic carbon (DIC) supply in neutral to alkaline waters and is related to variations in aquatic productivity (Papadimitriou et al. in Limnol Oceanogr 50:1084-1095, 2005). This paper aims to evaluate the usage of epsilon(macrophyte) inferred from isotopic analyses of submerged plant fossils in addition to analyses of lake carbonate as a palaeolimnological proxy for former HCO3 (-) concentrations. Stable carbon isotopic analysis of modern Potamogeton pectinatus leaves and its host water DIC from the Tibetan Plateau and Central Yakutia (Russia) yielded values between -23.3 and +0.4aEuro degrees and between +14.0 and +6.5aEuro degrees, respectively. Values of epsilon (Potamogeton-DIC) (range -15.4 to +1.1aEuro degrees) from these lakes are significantly correlated with host water HCO3 (-) concentration (range 78-2,200 mg/l) (r = -0.86; P < 0.001), thus allowing for the development of a transfer function. Palaeo-epsilon (Potamogeton-ostracods) values from Luanhaizi Lake on the NE Tibetan Plateau, as inferred from the stable carbon isotope measurement of fossil Potamogeton pectinatus seeds (range -24 to +2.8aEuro degrees) and ostracods (range -7.8 to +7.5\%) range between -14.8 and 1.6aEuro degrees. Phases of assumed disequilibrium between delta C-13(DIC) and delta C-13(ostracods) known to occur in charophyte swards (as indicated by the deposition of charophyte fossils) were excluded from the analysis of palaeo-epsilon. The application of the epsilon (Potamogeton-DIC)-HCO3 (-) transfer function yielded a median palaeo-HCO3 (-) -concentration of 290 mg/l. Variations in the dissolved organic carbon supply compare well with aquatic plant productivity changes and lake level variability as inferred from a multiproxy study of the same record including analyses of plant macrofossils, ostracods, carbonate and organic content.}, language = {en} } @article{HerzschuhBirksMischkeetal.2010, author = {Herzschuh, Ulrike and Birks, H. John B. and Mischke, Steffen and Zhang, Chengjun and B{\"o}hner, J{\"u}rgen}, title = {A modern pollen-climate calibration set based on lake sediments from the Tibetan Plateau and its application to a Late Quaternary pollen record from the Qilian Mountains}, issn = {0305-0270}, doi = {10.1111/j.1365-2699.2009.02245.x}, year = {2010}, abstract = {Aim: Fossil pollen spectra from lake sediments on the Tibetan Plateau have been used for qualitative climate reconstruction, but no modern pollen-climate calibration set based on lake sediments is available to infer past climate quantitatively. This study aims to develop such a dataset and apply it to fossil data. Location: The Tibetan Plateau, between 30 and 40 degrees N and 87 and 103 degrees E. Methods: We collected surface sediments from 112 lakes and analysed them palynologically. The lakes span a wide range of mean annual precipitation (P-ann; 31-1022 mm), mean annual temperature (T-ann; -6.5 to 1 degrees C), and mean July temperature (T-July; 2.6-19.7 degrees C). Redundancy analysis showed that the modern pollen spectra are characteristic of their respective vegetation types and local climate. Transfer functions for P-ann, T-ann and T-July were developed with weighted averaging partial least squares. Model performance was assessed by leave-one-out cross-validation. Results: The root mean square errors of prediction (RMSEP) were 104 mm (P-ann), 1.18 degrees C (T-ann) and 1.17 degrees C (T-July). The RMSEPs, when expressed as percentages of the gradient sampled, were 10.6\% (P-ann), 15.7\% (T-ann) and 11.9\% (T-July). These low values indicate the good performance of our models. An application of the models to fossil pollen spectra covering the last c. 50 kyr yielded realistic results for Luanhaizi Lake in the Qilian Mountains on the north-eastern Tibetan Plateau (modern P-ann 480 mm; T-ann-1 degrees C). T-ann and P-ann values similar to present ones were reconstructed for late Marine Isotope Stage 3, with minimum values for the Last Glacial Maximum (c. 300 mm and 2 degrees C below present), and maximum values for the early Holocene (c. 70 mm and 0.5 degrees C greater than present). Main conclusions: The modern pollen-climate calibration set will potentially be useful for quantitative climate reconstructions from lake-sediment pollen spectra from the Tibetan Plateau, an area of considerable climatic and biogeographical importance.}, language = {en} } @article{HerzschuhBirksLiuetal.2010, author = {Herzschuh, Ulrike and Birks, H. John B. and Liu, Xingqi and Kubatzki, Claudia and Lohmann, Gerrit}, title = {Retracted: What caused the mid-Holocene forest decline on the eastern Tibet-Qinghai Plateau?}, issn = {1466-822X}, doi = {10.1111/j.1466-8238.2009.00501.x}, year = {2010}, abstract = {Aim: Atmospheric CO2 concentrations depend, in part, on the amount of biomass locked up in terrestrial vegetation. Information on the causes of a broad-scale vegetation transition and associated loss of biomass is thus of critical interest for understanding global palaeoclimatic changes. Pollen records from the north-eastern Tibet-Qinghai Plateau reveal a dramatic and extensive forest decline beginning c. 6000 cal. yr bp. The aim of this study is to elucidate the causes of this regional-scale change from high-biomass forest to low-biomass steppe on the Tibet-Qinghai Plateau during the second half of the Holocene. Location: Our study focuses on the north-eastern Tibet-Qinghai Plateau. Stratigraphical data used are from Qinghai Lake (3200 m a.s.l., 36 degrees 32'-37 degrees 15' N, 99 degrees 36'-100 degrees 47' E). Methods: We apply a modern pollen-precipitation transfer function from the eastern and north-eastern Tibet-Qinghai Plateau to fossil pollen spectra from Qinghai Lake to reconstruct annual precipitation changes during the Holocene. The reconstructions are compared to a stable oxygen-isotope record from the same sediment core and to results from two transient climate model simulations. Results: The pollen-based precipitation reconstruction covering the Holocene parallels moisture changes inferred from the stable oxygen-isotope record. Furthermore, these results are in close agreement with simulated model-based past annual precipitation changes. Main conclusions: In the light of these data and the model results, we conclude that it is not necessary to attribute the broad-scale forest decline to human activity. Climate change as a result of changes in the intensity of the East Asian Summer Monsoon in the mid-Holocene is the most parsimonious explanation for the widespread forest decline on the Tibet-Qinghai Plateau. Moreover, climate feedback from a reduced forest cover accentuates increasingly drier conditions in the area, indicating complex vegetation-climate interactions during this major ecological change.}, language = {en} } @article{HerzschuhBirks2010, author = {Herzschuh, Ulrike and Birks, H. John B.}, title = {Evaluating the indicator value of Tibetan pollen taxa for modern vegetation and climate}, issn = {0034-6667}, doi = {10.1016/j.revpalbo.2010.02.016}, year = {2010}, abstract = {Pollen taxa of known indicator value are of great potential in the qualitative interpretation of pollen diagrams. Here we apply several numerical approaches to a lake-sediment based pollen data-set from the eastern and central Tibetan Plateau (112 samples) to assess the indicator value of Tibetan pollen taxa for modem vegetation types and for modern climate. Results from Multi-Response Permutation Procedures indicate that the differences between groups of pollen spectra originating from the same vegetation type (temperate desert, temperate steppe, alpine desert, alpine steppe, high-alpine meadow, subalpine shrub, and patchy forest) are statistically significant. Indicator Species Analyses identify several indicator taxa for most vegetation types. Multivariate regression tree analysis indicates that about 390 mm of annual precipitation is the most critical threshold for the modern pollen spectra. This roughly separates desert and steppe vegetation from high-alpine meadow, subalpine shrub, and patchy forest vegetation. A strong pollen-climate relationship on the Tibetan Plateau is confirmed by the large number of statistically significant pollen taxa-climate (annual precipitation or/and annual temperature) relationships as evaluated by statistical response- modelling, involving generalised linear models.}, language = {en} } @article{HasslerKreylingBeierkuhnleinetal.2010, author = {Haßler, Sibylle Kathrin and Kreyling, J{\"u}rgen and Beierkuhnlein, Carl and Eisold, Jenny and Samimi, Cyrus and Wagenseil, H. and Jentsch, A.}, title = {Vegetation pattern divergence between dry and wet season in a semiarid savanna : spatio-temporal dynamics of plant diversity in northwest Namibia}, issn = {0140-1963}, doi = {10.1016/j.jaridenv.2010.05.021}, year = {2010}, abstract = {African savannas are primarily used as pastures and are subject to changes in climate and management strategies. For sustainable management of these landscapes ecological knowledge on seasonal and long-term variability in plant community composition and the availability of green biomass is essential. In this study, we assessed the effects of dry and wet season on species richness and beta diversity for three sites along a gradient of increasing vegetation cover and precipitation in northwest Namibia. A hexagonal systematic sampling design was used to record floristic data. The Simple Matching, Soerensen, and multi-plot similarity coefficient and distance decay analyses were applied for examining beta diversity. Analyses were repeated while separating the plots according to the presence of woody vegetation. Species richness nearly doubled from dry to wet season; compositional similarity increased from dry to wet season and with increasing aridity of the study sites: distance decay was more pronounced in the dry season without any link to the precipitation gradient. Woody elements in the landscape, which occur along drainage lines or as tree islands, govern spatial and seasonal plant diversity fluctuations. Monitoring them is important for conservation strategies and for establishing grazing rules that ensure a sustainable use of savanna ecosystems.}, language = {en} } @article{HattermannLevermann2010, author = {Hattermann, Tore and Levermann, Anders}, title = {Response of Southern Ocean circulation to global warming may enhance basal ice shelf melting around Antarctica}, issn = {0930-7575}, doi = {10.1007/s00382-009-0643-3}, year = {2010}, abstract = {We investigate the large-scale oceanic features determining the future ice shelf-ocean interaction by analyzing global warming experiments in a coarse resolution climate model with a comprehensive ocean component. Heat and freshwater fluxes from basal ice shelf melting (ISM) are parameterized following Beckmann and Goosse [Ocean Model 5(2):157-170, 2003]. Melting sensitivities to the oceanic temperature outside of the ice shelf cavities are varied from linear to quadratic (Holland et al. in J Clim 21, 2008). In 1\% per year CO2-increase experiments the total freshwater flux from ISM triples to 0.09 Sv in the linear case and more than quadruples to 0.15 Sv in the quadratic case after 140 years at which 4 x 280 ppm = 1,120 ppm was reached. Due to the long response time of subsurface temperature anomalies, ISM thereafter increases drastically, if CO2 concentrations are kept constant at 1,120 ppm. Varying strength of the Antarctic circumpolar current (ACC) is crucial for ISM increase, because southward advection of heat dominates the warming along the Antarctic coast. On centennial timescales the ACC accelerates due to deep ocean warming north of the current, caused by mixing of heat along isopycnals in the Southern Ocean (SO) outcropping regions. In contrast to previous studies we find an initial weakening of the ACC during the first 150 years of warming. This purely baroclinic effect is due to a freshening in the SO which is consistent with present observations. Comparison with simulations with diagnosed ISM but without its influence on the ocean circulation reveal a number of ISM-related feedbacks, of which a negative ISM-feedback, due to the ISM-related local oceanic cooling, is the dominant one.}, language = {en} } @article{HandySchmidBousquetetal.2010, author = {Handy, Mark R. and Schmid, Stefan M. and Bousquet, Romain and Kissling, Eduard and Bernoulli, Daniel}, title = {Reconciling plate-tectonic reconstructions of Alpine Tethys with the geological-geophysical record of spreading and subduction in the Alps}, issn = {0012-8252}, doi = {10.1016/j.earscirev.2010.06.002}, year = {2010}, abstract = {A new reconstruction of Alpine Tethys combines plate-kinematic modelling with a wealth of geological data and seismic tomography to shed light on its evolution, from sea-floor spreading through subduction to collision in the Alps. Unlike previous models, which relate the fate of Alpine Tethys solely to relative motions of Africa, Iberia and Europe during opening of the Atlantic, our reconstruction additionally invokes independent microplates whose motions are constrained primarily by the geological record. The motions of these microplates (Adria, Iberia, Alcapia, Alkapecia, and Tiszia) relative to both Africa and Europe during Late Cretaceous to Cenozoic time involved the subduction of remnant Tethyan basins during the following three stages that are characterized by contrasting plate motions and driving forces: (1) 131-84 Ma intra-oceanic subduction of the Ligurian part of Alpine Tethys attached to Iberia coincided with Eo-alpine orogenesis in the Alcapia microplate, north of Africa. These events were triggered primarily by foundering of the older (170-131 Ma) Neotethyan subduction slab along the NE margin of the composite African-Adriatic plate; subduction was linked by a sinistral transform system to E-W opening of the Valais part of Alpine Tethys; (2) 84-35 Ma subduction of primarily the Piemont and Valais parts of Alpine Tethys which were then attached to the European plate beneath the overriding African and later Adriatic plates. NW translation of Adria with respect to Africa was accommodated primarily by slow widening of the Ionian Sea; (3) 35 Ma-Recent rollback subduction of the Ligurian part of Alpine Tethys coincided with Western Alpine orogenesis and involved the formation of the Gibraltar and Calabrian arcs. Rapid subduction and arc formation were driven primarily by the pull of the gravitationally unstable, retreating Adriatic and African slabs during slow convergence of Africa and Europe. The upper European-Iberian plate stretched to accommodate this slab retreat in a very mobile fashion, while the continental core of the Adriatic microplate acted as a rigid indenter within the Alpine collisional zone. The subducted lithosphere in this reconstruction can be correlated with slab material imaged by seismic tomography beneath the Alps and Apennines, as well as beneath parts of the Pannonian Basin, the Adriatic Sea, the Ligurian Sea, and the Western Mediterranean. The predicted amount of subducted lithosphere exceeds the estimated volume of slab material residing at depth by some 10-30\%, indicating that parts of slabs may be superposed within the mantle transition zone and/or that some of this subducted lithosphere became seismically transparent.}, language = {en} } @phdthesis{Guse2010, author = {Guse, Bj{\"o}rn Felix}, title = {Improving flood frequency analysis by integration of empirical and probabilistic regional envelope curves}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-49265}, school = {Universit{\"a}t Potsdam}, year = {2010}, abstract = {Flood design necessitates discharge estimates for large recurrence intervals. However, in a flood frequency analysis, the uncertainty of discharge estimates increases with higher recurrence intervals, particularly due to the small number of available flood data. Furthermore, traditional distribution functions increase unlimitedly without consideration of an upper bound discharge. Hence, additional information needs to be considered which is representative for high recurrence intervals. Envelope curves which bound the maximum observed discharges of a region are an adequate regionalisation method to provide additional spatial information for the upper tail of a distribution function. Probabilistic regional envelope curves (PRECs) are an extension of the traditional empirical envelope curve approach, in which a recurrence interval is estimated for a regional envelope curve (REC). The REC is constructed for a homogeneous pooling group of sites. The estimation of this recurrence interval is based on the effective sample years of data considering the intersite dependence among all sites of the pooling group. The core idea of this thesis was an improvement of discharge estimates for high recurrence intervals by integrating empirical and probabilistic regional envelope curves into the flood frequency analysis. Therefore, the method of probabilistic regional envelope curves was investigated in detail. Several pooling groups were derived by modifying candidate sets of catchment descriptors and settings of two different pooling methods. These were used to construct PRECs. A sensitivity analysis shows the variability of discharges and the recurrence intervals for a given site due to the different assumptions. The unit flood of record which governs the intercept of PREC was determined as the most influential aspect. By separating the catchments into nested and unnested pairs, the calculation algorithm for the effective sample years of data was refined. In this way, the estimation of the recurrence intervals was improved, and therefore the use of different parameter sets for nested and unnested pairs of catchments is recommended. In the second part of this thesis, PRECs were introduced into a distribution function. Whereas in the traditional approach only discharge values are used, PRECs provide a discharge and its corresponding recurrence interval. Hence, a novel approach was developed, which allows a combination of the PREC results with the traditional systematic flood series while taking the PREC recurrence interval into consideration. An adequate mixed bounded distribution function was presented, which in addition to the PREC results also uses an upper bound discharge derived by an empirical envelope curve. By doing so, two types of additional information which are representative for the upper tail of a distribution function were included in the flood frequency analysis. The integration of both types of additional information leads to an improved discharge estimation for recurrence intervals between 100 and 1000 years.}, language = {en} } @article{GraeffZeheSchlaegeretal.2010, author = {Gr{\"a}ff, Thomas and Zehe, Erwin and Schl{\"a}ger, Stefan and Morgner, Markus and Bauer, Andreas and Becker, Rolf and Creutzfeldt, Benjamin and Bronstert, Axel}, title = {A quality assessment of spatial TDR soil moisture measurements in homogenous and heterogeneous media with laboratory experiments}, issn = {1812-2108}, doi = {10.5194/hessd-7-269-2010}, year = {2010}, abstract = {Investigation of transient soil moisture profiles yields valuable information of near- surface processes. A recently developed reconstruction algorithm based on the telegraph equation allows the inverse estimation of soil moisture profiles along coated, three rod TDR probes. Laboratory experiments were carried out to prove the results of the inversion and to understand the influence of probe rod deformation and solid objects close to the probe in heterogonous media. Differences in rod geometry can lead to serious misinterpretations in the soil moisture profile but have small influence on the average soil moisture along the probe. Solids in the integration volume have almost no effect on average soil moisture but result in locally slightly decreased moisture values. Inverted profiles obtained in a loamy soil with a clay content of about 16\% were in good agreement with independent measurements.}, language = {en} } @article{GraeffZeheSchlaegeretal.2010, author = {Gr{\"a}ff, Thomas and Zehe, Erwin and Schlaeger, Stefan and Morgner, Markus and Bauer, Andreas and Becker, Rolf and Creutzfeldt, Benjamin and Bronstert, Axel}, title = {A quality assessment of Spatial TDR soil moisture measurements in homogenous and heterogeneous media with laboratory experiments}, issn = {1027-5606}, doi = {10.5194/hess-14-1007-2010}, year = {2010}, abstract = {Investigation of transient soil moisture profiles yields valuable information of near- surface processes. A recently developed reconstruction algorithm based on the telegraph equation allows the inverse estimation of soil moisture profiles along coated, three rod TDR probes. Laboratory experiments were carried out to prove the results of the inversion and to understand the influence of probe rod deformation and solid objects close to the probe in heterogeneous media. Differences in rod geometry can lead to serious misinterpretations in the soil moisture profile, but have small influence on the average soil moisture along the probe. Solids in the integration volume have almost no effect on average soil moisture, but result in locally slightly decreased moisture values. Inverted profiles obtained in a loamy soil with a clay content of about 16\% were in good agreement with independent measurements.}, language = {en} } @article{GerkeKoszinskiKalettkaetal.2010, author = {Gerke, Horst H. and Koszinski, Sylvia and Kalettka, Thomas and Sommer, Michael}, title = {Structures and hydrologic function of soil landscapes with kettle holes using an integrated hydropedological approach}, issn = {0022-1694}, doi = {10.1016/j.jhydrol.2009.12.047}, year = {2010}, abstract = {The hummocky post-glacial soil landscapes with kettle holes as internal drainage systems are characterized by ponds that trap lateral fluxes in topographic depressions. A quantitative description is mostly limited by the unknown complexity of hydraulically relevant soil and sediment structures. This paper is focussing on a structure-based approach to identify relevant field-scale flow and transport processes. Illustrative examples demonstrate extreme variations in water table fluctuation for adjoining kettle holes. Explanations require a pedohydrologic concept of the arable soil landscape. Identification of structures is based on geophysical methods and soil hydraulic measurements. Electrical resistivity imaging yields 0.5 m-scale spatial structures that correspond with soil texture distributions. Electromagnetic induction provides larger-scale field maps that reflect major soil and sediment features. Results of both methods correspond within the limits of the different spatial resolutions. With geophysical exploration methods, colluvial areas with textural differences between upper and deeper soil layers, coarse-textured sediment lenses, and stony colluvial regions around kettle holes are identified as potentially relevant flow structures. The colluvial fringe around the pond seems to be a sensitive area with important lateral exchange fluxes. Tensiometer measurements perpendicular to this boundary indicate hydraulic gradients directed from the pond towards the partially saturated soil. The localized infiltration of trapped water in kettle holes can control large fractions of ground water recharge and may have implications for the fate of agricultural chemicals in post-glacial landscapes. While surface and subsurface hydraulic structures may be inferred using minimal-invasive techniques, better understanding of processes and properties governing lateral exchange fluxes between pond and surrounding soil are required.}, language = {en} } @article{FoersterFoersterOberhaenslietal.2010, author = {F{\"o}rster, Hans-J{\"u}rgen and F{\"o}rster, Andrea and Oberh{\"a}nsli, Roland and Stromeyer, Dietrich}, title = {Lithospheric composition and thermal structure of the Arabian Shield in Jordan}, issn = {0040-1951}, doi = {10.1016/j.tecto.2008.11.014}, year = {2010}, abstract = {In this paper, a unique set of samples from the uppermost crust down to the lithospheric mantle of Jordan is analyzed for composition and petrophysical properties (density. thermal conductivity, radiogenic heat production) These data, covering a vertical section of almost 65 km. are used in conjunction with surface heat flow to generate a detailed and comprehensive lithospheric thermal model that reflects the conditions of the Arabian Shield (AS) prior to the post- Oligocene onset of lithosphere thinning and Voluminous basaltic volcanism. The pre-Miocene model geotherms, based on conductive surface heat flows of 55 and 60 mW m(-2). (a) meet the range of lithosphere-asthenosphere boundary depths of 110-160 km known from seismology, (b) conform to results of thermomechanical models on the on.-in of the Dead Sea basin that started in Miocene time. and (c) are consistent with typical xenolith-derived geotherms for terranes of similar age and lithospheric thickness. Moho temperatures (at depths between 35 and 40 km) of the AS in pre-Miocene times were most likely in the order of 530-650 degrees C, with mantle heat flows averaging between 24 and 29 mW m(-2) Results contradict former views of the late Proterozoic/early Cambrian-stabilized AS being an anomalously cold terrane A "cold" thermal structure inferred from previously measured low surface heat flows (generally <= 45 mW m(-2)) is inconsistent with the thickness, composition, and petrophysical properties of the stable lithosphere of the shield.}, language = {en} } @article{ForbrichKutzbachHormannetal.2010, author = {Forbrich, Inke and Kutzbach, Lars and Hormann, Annabell and Wilmking, Martin}, title = {A comparison of linear and exponential regression for estimating diffusive CH4 fluxes by closed-chambers in peatlands}, issn = {0038-0717}, doi = {10.1016/j.soilbio.2009.12.004}, year = {2010}, abstract = {The closed-chamber method is the most common approach to determine CH4 fluxes in peatlands. The concentration change in the chamber is monitored over time, and the flux is usually calculated by the slope of a linear regression function. Theoretically, the gas exchange cannot be constant over time but has to decrease, when the concentration gradient between chamber headspace and soil air decreases. In this study, we test whether we can detect this non- linearity in the concentration change during the chamber closure with six air samples. We expect generally a low concentration gradient on dry sites (hummocks) and thus the occurrence of exponential concentration changes in the chamber due to a quick equilibrium of gas concentrations between peat and chamber headspace. On wet (flarks) and sedge- covered sites (lawns), we expect a high gradient and near-linear concentration changes in the chamber. To evaluate these model assumptions, we calculate both linear and exponential regressions for a test data set (n = 597) from a Finnish mire. We use the Akaike Information Criterion with small sample second order bias correction to select the best-fitted model. 13.6\%, 19.2\% and 9.8\% of measurements on hummocks, lawns and flarks, respectively, were best fitted with an exponential regression model. A flux estimation derived from the slope of the exponential function at the beginning of the chamber closure can be significantly higher than using the slope of the linear regression function. Non-linear concentration-overtime curves occurred mostly during periods of changing water table. This could be due to either natural processes or chamber artefacts, e.g. initial pressure fluctuations during chamber deployment. To be able to exclude either natural processes or artefacts as cause of non-linearity, further information, e.g. CH4 concentration profile measurements in the peat, would be needed. If this is not available, the range of uncertainty can be substantial. We suggest to use the range between the slopes of the exponential regression at the beginning and at the end of the closure time as an estimate of the overall uncertainty.}, language = {en} } @article{FariasVargasTassaraetal.2010, author = {Far{\"i}as, Marcelo and Vargas, Gabriel and Tassara, Andr{\´e}s and Carretier, S{\´e}bastien and Baize, St{\´e}phane and Melnick, Daniel and Bataille, Klaus}, title = {Land-level changes produced by the M-w 8.8 2010 Chilean earthquake}, issn = {0036-8075}, doi = {10.1126/science.1192094}, year = {2010}, abstract = {We observed vertically displaced coastal and river markers after the 27 February 2010 Chilean earthquake [moment magnitude (Mw) 8.8]. Land-level changes range between 2.5 and -1 meters, evident along an ~500-kilometers- long segment identified here as the maximum length of coseismic rupture. A hinge line located 120 kilometers from the trench separates uplifted areas, to the west, from subsided regions. A simple elastic dislocation model fits these observations well; model parameters give a similar seismic moment to seismological estimates and suggest that most of the plate convergence since the 1835 great earthquake was elastically stored and then released during this event.}, language = {en} } @article{EndrunOhrnbergerSavvaidis2010, author = {Endrun, Brigitte and Ohrnberger, Matthias and Savvaidis, Alexandros}, title = {On the repeatability and consistency of three-component ambient vibration array measurements}, issn = {1570-761X}, doi = {10.1007/s10518-009-9159-9}, year = {2010}, abstract = {Ambient vibration measurements with small, temporary arrays that produce estimates of surface wave dispersion have become increasingly popular as a low-cost, non-invasive tool for site characterisation. An important requirement for these measurements to be meaningful, however, is the temporal consistency and repeatability of the resulting dispersion and spatial autocorrelation curve estimates. Data acquired within several European research projects (NERIES task JRA4, SESAME, and other multinational experiments) offer the chance to investigate the variability of the derived data products. The dataset analysed here consists of repeated array measurements, with several years of time elapsed between them. The measurements were conducted by different groups in different seasons, using different instrumentations and array layouts, at six sites in Greece and Italy. Ambient vibration amplitude spectra and locations of dominant sources vary between the two measurements at each location. Still, analysis indicates that this does not influence the derived dispersion information, which is stable in time and neither influenced by the instrumentation nor the analyst. The frequency range over which the dispersion curves and spatial autocorrelation curves can be reliably estimated depends on the array dimensions (minimum and maximum aperture) used in the specific deployment, though, and may accordingly vary between the repeated experiments. The relative contribution of Rayleigh and Love waves to the wavefield can likewise change between repeated measurements. The observed relative contribution of Rayleigh waves is generally at or below 50\%, with especially low values for the rural sites. Besides, the visibility of higher modes depends on the noise wavefield conditions. The similarity of the dispersion and autocorrelation curves measured at each site indicates that the curves are stable, mainly determined by the sub-surface structure, and can thus be used to derive velocity information with depth. Differences between velocity models for the same site derived from independently determined dispersion and autocorrelation curves-as observed in other studies-are consequently not adequately explained by uncertainties in the measurement part.}, language = {en} } @article{EhretZehe2010, author = {Ehret, Uwe and Zehe, Erwin}, title = {Series distance : an intuitive metric for hydrograph comparison}, issn = {1812-2108}, doi = {10.5194/hessd-7-8387-2010}, year = {2010}, abstract = {Applying metrics for hydrograph comparison is a central task in hydrological modelling, used both in model calibration and the evaluation of simulations or forecasts. Motivated by the shortcomings of standard objective metrics such as the Root Mean Square Error or the Mean Peak Time Error and the advantages of visual inspection as a powerful tool for simultaneous, case-specific and multi-criteria (yet subjective) evaluation, we propose a new objective metric termed Series Distance, which is in close accordance with visual evaluation. The Series Distance is an event-based method and consists of three parts, namely a Threat Score which evaluates overall agreement of event occurrence, and the overall distance of matching observed and simulated events with respect to amplitude and timing. The novelty of the latter two is the way in which matching point pairs on the observed and simulated hydrographs are identified, namely by the same relative position in matching segments (rise or recession) of matching events. Thus, amplitude and timing errors are calculated simultaneously but separately, from point pairs that also match visually, considering complete events rather than only individual points (which is for example the case with metrics related to Peak Time Errors). After presenting the Series Distance theory, we discuss its properties and compare it to those of standard metrics and visual inspection, both at the example of simple, artificial hydrographs and an ensemble of realistic forecasts. The results suggest that the Series Distance compares and evaluates hydrographs in a way comparable to visual inspection, but in an objective, reproducible way.}, language = {en} } @article{DieterElsenbeerTurner2010, author = {Dieter, Daniela and Elsenbeer, Helmut and Turner, Benjamin L.}, title = {Phosphorus fractionation in lowland tropical rainforest soils in central Panama}, issn = {0341-8162}, doi = {10.1016/j.catena.2010.05.010}, year = {2010}, abstract = {Phosphorus availability is commonly assumed to limit productivity in lowland tropical rainforests, yet there is relatively little information on the chemical forms of soil phosphorus in such ecosystems. We used the Hedley sequential fractionation scheme to assess phosphorus chemistry in five soils supporting tropical rainforest on Barro Colorado Island, Republic of Panama. The soils represented a range of orders (Inceptisols, Alfisols, and Oxisols) formed on contrasting geological substrates and topography, but under uniform climate and vegetation. Total phosphorus in surface horizons ranged between 315 and 1114 mg P kg(-1), being lowest on a soil derived from marine sediments and highest on soils derived from andesite. The majority of the phosphorus occurred in recalcitrant forms, although between 14\% and 39\% occurred as organic phosphorus. Readily-available phosphate, as extracted by anion-exchange membranes, occurred in small concentrations (4-13 mg P kg(-1)), although labile phosphorus, defined as phosphate extracted by anion-exchange membrane plus inorganic and organic phosphorus extracted by 0.5 M NaHCO3, accounted for between 4.7\% and 11.4\% of the total soil phosphorus. Our results indicate a strong control of geology and topography on soil phosphorus in tropical rainforests, which may have important implications for understanding the diversity and distribution of plant species in such ecosystems. Further, some of the most common soils on Barro Colorado Island, including those on the 50 ha forest dynamics plot, are rich in phosphorus despite their relatively advanced stage of pedogenesis.}, language = {en} } @phdthesis{DiGiacomo2010, author = {Di Giacomo, Domenico}, title = {Determination of the energy magnitude ME : application to rapid response purposes and insights to regional/local variabilities}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-50768}, school = {Universit{\"a}t Potsdam}, year = {2010}, abstract = {Recent large earthquakes put in evidence the need of improving and developing robust and rapid procedures to properly calculate the magnitude of an earthquake in a short time after its occurrence. The most famous example is the 26 December 2004 Sumatra earthquake, when the limitations of the standard procedures adopted at that time by many agencies failed to provide accurate magnitude estimates of this exceptional event in time to launch early enough warnings and appropriate response. Being related to the radiated seismic energy ES, the energy magnitude ME is a good estimator of the high frequency content radiated by the source which goes into the seismic waves. However, a procedure to rapidly determine ME (that is to say, within 15 minutes after the earthquake occurrence) was required. Here it is presented a procedure able to provide in a rapid way the energy magnitude ME for shallow earthquakes by analyzing teleseismic P‑waves in the distance range 20-98. To account for the energy loss experienced by the seismic waves from the source to the receivers, spectral amplitude decay functions obtained from numerical simulations of Greens functions based on the average global model AK135Q are used. The proposed method has been tested using a large global dataset (~1000 earthquakes) and the obtained rapid ME estimations have been compared to other magnitude scales from different agencies. Special emphasis is given to the comparison with the moment magnitude MW, since the latter is very popular and extensively used in common seismological practice. However, it is shown that MW alone provide only limited information about the seismic source properties, and that disaster management organizations would benefit from a combined use of MW and ME in the prompt evaluation of an earthquake's tsunami and shaking potential. In addition, since the proposed approach for ME is intended to work without knowledge of the fault plane geometry (often available only hours after an earthquake occurrence), the suitability of this method is discussed by grouping the analyzed earthquakes according to their type of mechanism (strike-slip, normal faulting, thrust faulting, etc.). No clear trend is found from the rapid ME estimates with the different fault plane solution groups. This is not the case for the ME routinely determined by the U.S. Geological Survey, which uses specific radiation pattern corrections. Further studies are needed to verify the effect of such corrections on ME estimates. Finally, exploiting the redundancy of the information provided by the analyzed dataset, the components of variance on the single station ME estimates are investigated. The largest component of variance is due to the intra-station (record-to-record) error, although the inter-station (station-to-station) error is not negligible and is of several magnitude units for some stations. Moreover, it is shown that the intra-station component of error is not random but depends on the travel path from a source area to a given station. Consequently, empirical corrections may be used 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 the propagation effects.}, language = {en} } @article{DahmKuehnOhrnbergeretal.2010, author = {Dahm, Torsten and Kuehn, Daniela and Ohrnberger, Matthias and Kroeger, Jens and Wiederhold, Helga and Reuther, Claus-Dieter and Dehghani, Ali and Scherbaum, Frank}, title = {Combining geophysical data sets to study the dynamics of shallow evaporites in urban environments : application to Hamburg, Germany}, issn = {0956-540X}, doi = {10.1111/j.1365-246X.2010.04521.x}, year = {2010}, abstract = {Shallowly situated evaporites in built-up areas are of relevance for urban and cultural development and hydrological regulation. The hazard of sinkholes, subrosion depressions and gypsum karst is often difficult to evaluate and may quickly change with anthropogenic influence. The geophysical exploration of evaporites in metropolitan areas is often not feasible with active industrial techniques. We collect and combine different passive geophysical data as microgravity, ambient vibrations, deformation and hydrological information to study the roof morphology of shallow evaporites beneath Hamburg, Northern Germany. The application of a novel gravity inversion technique leads to a 3-D depth model of the salt diapir under study. We compare the gravity-based depth model to pseudo-depths from H/V measurements and depth estimates from small-scale seismological array data. While the general range and trend of the diapir roof is consistent, a few anomalous regions are identified where H/V pseudo-depths indicate shallower structures not observed in gravity or array data. These are interpreted by shallow residual caprock floaters and zones of increased porosity. The shallow salt structure clearly correlates with a relative subsidence in the order of 2 mm yr(-1). The combined interpretation of roof morphology, yearly subsidence rates, chemical analyses of groundwater and of hydraulic head in aquifers indicates that the salt diapir beneath Hamburg is subject to significant ongoing dissolution that may possibly affect subrosion depressions, sinkhole distribution and land usage. The combined analysis of passive geophysical data may be exemplary for the study of shallow evaporites beneath other urban areas.}, language = {en} } @article{CristianoPetrosinoSaccorottietal.2010, author = {Cristiano, Luigia and Petrosino, Simona and Saccorotti, Gilberto and Ohrnberger, Matthias and Scarpa, Roberto}, title = {Shear-wave velocity structure at Mt. Etna from inversion of Rayleigh-wave dispersion patterns (2 s < T < 20 s)}, issn = {1593-5213}, doi = {10.4401/Ag-4574}, year = {2010}, abstract = {In the present study, we investigated the dispersion characteristics of medium-to-long period Rayleigh waves (2 s < T < 20 s) using both single-station techniques (multiple-filter analysis, and phase-match filter) and multichannel techniques (horizontal slowness [p] and angular frequency [omega] stack, and cross-correlation) to determine the velocity structure for the Mt. Etna volcano. We applied these techniques to a dataset of teleseisms, as regional and local earthquakes recorded by two broad-band seismic arrays installed at Mt. Etna in 2002 and 2005, during two seismic surveys organized by the Istituto Nazionale di Geofisica e Vulcanologia (INGV), sezione di Napoli. The dispersion curves obtained showed phase velocities ranging from 1.5 km/s to 4.0 km/s in the frequency band 0.05 Hz to 0.45 Hz. We inverted the average phase velocity dispersion curves using a non-linear approach, to obtain a set of shear-wave velocity models with maximum resolution depths of 25 km to 30 km. Moreover, the presence of lateral velocity contrasts was checked by dividing the whole array into seven triangular sub-arrays and inverting the dispersion curves relative to each triangle.}, language = {en} } @phdthesis{Creutzfeldt2010, author = {Creutzfeldt, Noah Angelo Benjamin}, title = {The effect of water storages on temporal gravity measurements and the benefits for hydrology}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-48575}, school = {Universit{\"a}t Potsdam}, year = {2010}, abstract = {Temporal gravimeter observations, used in geodesy and geophysics to study variation of the Earth's gravity field, are influenced by local water storage changes (WSC) and - from this perspective - add noise to the gravimeter signal records. At the same time, the part of the gravity signal caused by WSC may provide substantial information for hydrologists. Water storages are the fundamental state variable of hydrological systems, but comprehensive data on total WSC are practically inaccessible and their quantification is associated with a high level of uncertainty at the field scale. This study investigates the relationship between temporal gravity measurements and WSC in order to reduce the hydrological interfering signal from temporal gravity measurements and to explore the value of temporal gravity measurements for hydrology for the superconducting gravimeter (SG) of the Geodetic Observatory Wettzell, Germany. A 4D forward model with a spatially nested discretization domain was developed to simulate and calculate the local hydrological effect on the temporal gravity observations. An intensive measurement system was installed at the Geodetic Observatory Wettzell and WSC were measured in all relevant storage components, namely groundwater, saprolite, soil, top soil and snow storage. The monitoring system comprised also a suction-controlled, weighable, monolith-filled lysimeter, allowing an all time first comparison of a lysimeter and a gravimeter. Lysimeter data were used to estimate WSC at the field scale in combination with complementary observations and a hydrological 1D model. Total local WSC were derived, uncertainties were assessed and the hydrological gravity response was calculated from the WSC. A simple conceptual hydrological model was calibrated and evaluated against records of a superconducting gravimeter, soil moisture and groundwater time series. The model was evaluated by a split sample test and validated against independently estimated WSC from the lysimeter-based approach. A simulation of the hydrological gravity effect showed that WSC of one meter height along the topography caused a gravity response of 52 µGal, whereas, generally in geodesy, on flat terrain, the same water mass variation causes a gravity change of only 42 µGal (Bouguer approximation). The radius of influence of local water storage variations can be limited to 1000 m and 50 \% to 80 \% of the local hydro¬logical gravity signal is generated within a radius of 50 m around the gravimeter. At the Geodetic Observatory Wettzell, WSC in the snow pack, top soil, unsaturated saprolite and fractured aquifer are all important terms of the local water budget. With the exception of snow, all storage components have gravity responses of the same order of magnitude and are therefore relevant for gravity observations. The comparison of the total hydrological gravity response to the gravity residuals obtained from the SG, showed similarities in both short-term and seasonal dynamics. However, the results demonstrated the limitations of estimating total local WSC using hydrological point measurements. The results of the lysimeter-based approach showed that gravity residuals are caused to a larger extent by local WSC than previously estimated. A comparison of the results with other methods used in the past to correct temporal gravity observations for the local hydrological influence showed that the lysimeter measurements improved the independent estimation of WSC significantly and thus provided a better way of estimating the local hydrological gravity effect. In the context of hydrological noise reduction, at sites where temporal gravity observations are used for geophysical studies beyond local hydrology, the installation of a lysimeter in combination with complementary hydrological measurements is recommended. From the hydrological view point, using gravimeter data as a calibration constraint improved the model results in comparison to hydrological point measurements. Thanks to their capacity to integrate over different storage components and a larger area, gravimeters provide generalized information on total WSC at the field scale. Due to their integrative nature, gravity data must be interpreted with great care in hydrological studies. However, gravimeters can serve as a novel measurement instrument for hydrology and the application of gravimeters especially designed to study open research questions in hydrology is recommended.}, language = {en} } @article{ClarkeBurbank2010, author = {Clarke, Brian A. and Burbank, Douglas W.}, title = {Bedrock fracturing, threshold hillslopes, and limits to the magnitude of bedrock landslides}, issn = {0012-821X}, doi = {10.1016/j.epsl.2010.07.011}, year = {2010}, abstract = {Bedrock fracturing and rock strength are widely believed to influence landscape morphology and erosional resistance. Yet, understanding of the quantitative relationship between rock-mass strength and landscape evolution remains limited. Here we present a new application of seismic refraction surveys that uses variations in seismic velocity to interpret differences in bedrock fracture density with depth. We use a comparative study of Fiordland and the western Southern Alps of New Zealand to examine how differences in rock type and bedrock fracturing influence landscape morphology and landslide response to rock uplift. In both regions, slopes appear invariant with differential rock-uplift rates and slope distributions reveal modal hillslope angles of similar to 32 degrees. The majority of landslides initiate on slopes steeper than the modal hillslope angle, however, landslide magnitude-frequency distributions reveal order-of-magnitude differences between the regions, with Fiordland experiencing considerably smaller and less frequent landsliding events. Landslide-driven denudation rates of similar to 9 mm/yr in the western Southern Alps and between similar to 0.1 and 0.3 mm/yr in Fiordland approximate estimates of long-term rock-uplift rates for each region. The invariance of hillslope angles, near-normal slope distributions, predominance of landslide initiation on slopes steeper than modal values, and the apparent balance between rates of uplift and landslide-driven erosion suggest that hillslopes in both regions are at threshold angles. Their similar modal slopes further suggest that both ranges are characterized by equivalent rock-mass strength, despite striking differences in lithology. Additionally, our seismic analysis reveals nearly identical surface p-wave velocities. The unexpected equivalence of both modal slopes and surface velocities between these lithologically distinct ranges is attributed to contrasting degrees of surface fracturing that have differentially affected the intact rock properties, such that they now yield equivalent surface velocities and hillslope-scale strengths. Given that surface fractures help regulate threshold angles by modulating hillslope strength; we propose that shallow seismic velocities may provide a quantitative proxy for rock-mass strength. We define two contrasting fracture and landsliding environments. In Fiordland, dense geomorphic fracturing that is focused within the shallow subsurface appears to limit the depth and magnitude of most bedrock landslides. Conversely, in the western Southern Alps, tectonic forces produce pervasive fracturing with depth that results in larger, and deeper landslides. Our data suggest that bedrock fracturing at the Earth's surface modulates threshold hillslope angles, whereas the depth of bedrock fracturing influences the magnitude and frequency of landslide response to tectonic rock uplift.}, language = {en} } @article{CarrapaHauerSchoenbohmetal.2010, author = {Carrapa, Barbara and Hauer, Joern and Schoenbohm, Lindsay M. and Strecker, Manfred and Schmitt, Axel K. and Villanueva, Arturo and Gomez, Jos{\´e} Sosa}, title = {Dynamics of deformation and sedimentation in the northern Sierras Pampeanas : an integrated study of the Neogene Fiambal{\´a} basin, NW Argentina ; reply}, issn = {0016-7606}, doi = {10.1130/B30134.1}, year = {2010}, language = {en} } @article{CacaceKaiserLewerenzetal.2010, author = {Cacace, Mauro and Kaiser, Bjoern Onno and Lewerenz, Bjoern and Scheck-Wenderoth, Magdalena}, title = {Geothermal energy in sedimentary basins : what we can learn from regional numerical models}, issn = {0009-2819}, doi = {10.1016/j.chemer.2010.05.017}, year = {2010}, abstract = {Understanding the interactions between the different processes that control the geothermal and fluid flow fields in sedimentary basins is crucial for exploitation of geothermal energy. Numerical models provide predictive and feasible information for a correct assessment of geothermal resources especially in areas where data acquisition is demanding. Here, we present results from numerical efforts to characterize the thermal structure and its interaction with the fluid system for the area of the North East German Basin (NEGB). The relative impact of the different (diffusive and advective) processes affecting the hydrothermal setting of the basin are investigated by means of three- dimensional numerical simulations. Lithospheric-scale numerical models are evaluated to understand the specific thermal signature of the relevant factors influencing the present-day conductive geothermal field in the NEGB. Shallow and deep structural controls on the thermal configuration of the basin are addressed and quantified. Interaction between the resulting thermal field and the active fluid system is investigated by means of three-dimensional simulations of coupled fluid flow and heat transport. Factors influencing stability and reliability of modeling predictions are discussed. The main effort is to build a physically consistent model for the basin which integrates the impacts of thermal gradients on the regional fluid regime and their coupling with the main geological units defining the basin.}, language = {en} } @phdthesis{Boeniger2010, author = {B{\"o}niger, Urs}, title = {Attributes and their potential to analyze and interpret 3D GPR data}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-50124}, school = {Universit{\"a}t Potsdam}, year = {2010}, abstract = {Based on technological advances made within the past decades, ground-penetrating radar (GPR) has become a well-established, non-destructive subsurface imaging technique. Catalyzed by recent demands for high-resolution, near-surface imaging (e.g., the detection of unexploded ordnances and subsurface utilities, or hydrological investigations), the quality of today's GPR-based, near-surface images has significantly matured. At the same time, the analysis of oil and gas related reflection seismic data sets has experienced significant advances. Considering the sensitivity of attribute analysis with respect to data positioning in general, and multi-trace attributes in particular, trace positioning accuracy is of major importance for the success of attribute-based analysis flows. Therefore, to study the feasibility of GPR-based attribute analyses, I first developed and evaluated a real-time GPR surveying setup based on a modern tracking total station (TTS). The combination of current GPR systems capability of fusing global positioning system (GPS) and geophysical data in real-time, the ability of modern TTS systems to generate a GPS-like positional output and wireless data transmission using radio modems results in a flexible and robust surveying setup. To elaborate the feasibility of this setup, I studied the major limitations of such an approach: system cross-talk and data delays known as latencies. Experimental studies have shown that when a minimal distance of ~5 m between the GPR and the TTS system is considered, the signal-to-noise ratio of the acquired GPR data using radio communication equals the one without radio communication. To address the limitations imposed by system latencies, inherent to all real-time data fusion approaches, I developed a novel correction (calibration) strategy to assess the gross system latency and to correct for it. This resulted in the centimeter trace accuracy required by high-frequency and/or three-dimensional (3D) GPR surveys. Having introduced this flexible high-precision surveying setup, I successfully demonstrated the application of attribute-based processing to GPR specific problems, which may differ significantly from the geological ones typically addressed by the oil and gas industry using seismic data. In this thesis, I concentrated on archaeological and subsurface utility problems, as they represent typical near-surface geophysical targets. Enhancing 3D archaeological GPR data sets using a dip-steered filtering approach, followed by calculation of coherency and similarity, allowed me to conduct subsurface interpretations far beyond those obtained by classical time-slice analyses. I could show that the incorporation of additional data sets (magnetic and topographic) and attributes derived from these data sets can further improve the interpretation. In a case study, such an approach revealed the complementary nature of the individual data sets and, for example, allowed conclusions about the source location of magnetic anomalies by concurrently analyzing GPR time/depth slices to be made. In addition to archaeological targets, subsurface utility detection and characterization is a steadily growing field of application for GPR. I developed a novel attribute called depolarization. Incorporation of geometrical and physical feature characteristics into the depolarization attribute allowed me to display the observed polarization phenomena efficiently. Geometrical enhancement makes use of an improved symmetry extraction algorithm based on Laplacian high-boosting, followed by a phase-based symmetry calculation using a two-dimensional (2D) log-Gabor filterbank decomposition of the data volume. To extract the physical information from the dual-component data set, I employed a sliding-window principle component analysis. The combination of the geometrically derived feature angle and the physically derived polarization angle allowed me to enhance the polarization characteristics of subsurface features. Ground-truth information obtained by excavations confirmed this interpretation. In the future, inclusion of cross-polarized antennae configurations into the processing scheme may further improve the quality of the depolarization attribute. In addition to polarization phenomena, the time-dependent frequency evolution of GPR signals might hold further information on the subsurface architecture and/or material properties. High-resolution, sparsity promoting decomposition approaches have recently had a significant impact on the image and signal processing community. In this thesis, I introduced a modified tree-based matching pursuit approach. Based on different synthetic examples, I showed that the modified tree-based pursuit approach clearly outperforms other commonly used time-frequency decomposition approaches with respect to both time and frequency resolutions. Apart from the investigation of tuning effects in GPR data, I also demonstrated the potential of high-resolution sparse decompositions for advanced data processing. Frequency modulation of individual atoms themselves allows to efficiently correct frequency attenuation effects and improve resolution based on shifting the average frequency level. GPR-based attribute analysis is still in its infancy. Considering the growing widespread realization of 3D GPR studies there will certainly be an increasing demand towards improved subsurface interpretations in the future. Similar to the assessment of quantitative reservoir properties through the combination of 3D seismic attribute volumes with sparse well-log information, parameter estimation in a combined manner represents another step in emphasizing the potential of attribute-driven GPR data analyses.}, language = {en} } @phdthesis{Bulut2010, author = {Bulut, Fatih}, title = {Seismicity and fault kinematics of the North Anatolian fault zone along the 1999 Izmit rupture and its transition to the Marmara seismic gap}, pages = {VII, 98 S.}, year = {2010}, language = {en} } @article{BrunelleSigmanJaccardetal.2010, author = {Brunelle, Brigitte G. and Sigman, Daniel M. and Jaccard, Samuel Laurent and Keigwin, Lloyd D. and Plessen, Birgit and Schettler, Georg and Cook, Mea S. and Haug, Gerald H.}, title = {Glacial/interglacial changes in nutrient supply and stratification in the western subarctic North Pacific since the penultimate glacial maximum}, issn = {0277-3791}, doi = {10.1016/j.quascirev.2010.03.010}, year = {2010}, abstract = {In piston cores from the open subarctic Pacific and the Okhotsk Sea, diatom-bound delta N-15 (delta N-15(db)), biogenic opal, calcium carbonate, and barium were measured from coretop to the previous glacial maximum (MIS 6). Glacial intervals are generally characterized by high delta N-15(db) (similar to 8 parts per thousand) and low productivity, whereas interglacial intervals have a lower delta N-15(db) (5.7-6.3 parts per thousand) and indicate high biogenic productivity. These data extend the regional swath of evidence for nearly complete surface nutrient utilization during glacial maxima, consistent with stronger upper water column stratification throughout the subarctic region during colder intervals. An early deglacial decline in delta N-15(db) of 2 parts per thousand at similar to 17.5 ka, previously observed in the Bering Sea, is found here in the open subarctic Pacific record and arguably also in the Okhotsk, and a case can be made that a similar decrease in delta N-15(db) occurred in both regions at the previous deglaciation as well. The early deglacial delta N-15(db) decrease, best explained by a decrease in surface nutrient utilization, appears synchronous with southern hemisphere-associated deglacial changes and with the Heinrich 1 event in the North Atlantic. This delta N-15(db) decrease may signal the initial deglacial weakening in subarctic North Pacific stratification and/or a deglacial increase in shallow subsurface nitrate concentration. If the former, it would be the North Pacific analogue to the increase in vertical exchange inferred for the Southern Ocean at the time of Heinrich Event 1. In either case, the lack of any clear change in paleoproductivity proxies during this interval would seem to require an early deglacial decrease in the iron-to-nitrate ratio of subsurface nutrient supply or the predominance of light limitation of phytoplankton growth during the deglaciation prior to Bolling-Allerod warming.}, language = {en} } @article{BreitenbachAdkinsMeyeretal.2010, author = {Breitenbach, Sebastian Franz Martin and Adkins, Jess F. and Meyer, Hanno and Marwan, Norbert and Kumar, Kanikicharla Krishna and Haug, Gerald H.}, title = {Strong influence of water vapor source dynamics on stable isotopes in precipitation observed in Southern Meghalaya, NE India}, issn = {0012-821X}, doi = {10.1016/j.epsl.2010.01.038}, year = {2010}, abstract = {To calibrate delta O-18 time-series from speleothems in the eastern Indian summer monsoon (ISM) region of India, and to understand the moisture regime over the northern Bay of Bengal (BoB) we analyze the delta O-18 and delta D of rainwater, collected in 2007 and 2008 near Cherrapunji, India. delta D values range from + 18.5 parts per thousand to 144.4 parts per thousand, while delta O-18 varies between +0.8 parts per thousand and 18.8 parts per thousand. The Local Meteoric Water Line (LMWL) is found to be indistinguishable from the Global Meteoric Water Line (GMWL). Late ISM (September-October) rainfall exhibits lowest delta O-18 and delta D values, with little relationship to the local precipitation amount. There is a trend to lighter isotope values over the course of the ISM, but it does not correlate with the patterns of temperature and rainfall amount delta O-18 and delta D time-series have to be interpreted with caution in terms of the 'amount effect' in this subtropical region. We find that the temporal trend in delta O-18 reflects increasing transport distance during the ISM, isotopic changes in the northern BoB surface waters during late ISM, and vapor re-equilibration with rain droplets. Using an isotope box model for surface ocean waters, we quantify the potential influence of river runoff on the isotopic composition of the seasonal freshwater plume in the northern BoB. Temporal variations in this source can contribute up to 25\% of the observed changes in stable isotopes of precipitation in NE India. To delineate other moisture sources, we use backward trajectory computations and find a strong correlation between source region and isotopic composition. Palaeoclimatic stable isotope time-series from northeast Indian speleothems likely reflect changes in moisture source and transport pathway, as well as the isotopic composition of the BoB surface water, all of which in turn reflect ISM strength. Stalagmite records from the region can therefore be interpreted as integrated measures of the ISM strength.}, language = {en} } @article{BraeunigerKnappKuehnetal.2010, author = {Braeuniger, Claudia and Knapp, Sonja and Kuehn, Ingolf and Klotz, Stefan}, title = {Testing taxonomic and landscape surrogates for biodiversity in an urban setting}, issn = {0169-2046}, doi = {10.1016/j.landurbplan.2010.07.001}, year = {2010}, abstract = {Cities often have higher species diversity than the surrounding landscape. This diversity is important for both nature conservation and urban planning. The recreation of residents and the protection of species and habitats are simultaneous targets of maintaining urban green spaces. Data about the distribution and richness of species and their habitats have been compiled frequently; however, it is difficult and costly to measure the complete biodiversity of a region, necessitating useful surrogates. We tested species and habitat data in 27 protected areas in a Central German city and asked (1) whether the diversity of selected taxa acts as a surrogate for the diversity of other taxa and total investigated diversity, and (2) whether landscape structure and human impact explain species richness. Landscape structure metrics were based on soil and habitat types; human influence was measured as the degree of hemeroby. We tested and accounted for sample bias prior to analyses. (1) Vascular plant species richness explained total richness and single taxon richness best. (2) The size of a protected area was the most important predictor of species richness. After correcting for the effect of size, shape complexity, isolation, and matrix properties remained significant. Accordingly, the type of data frequently used for urban planning - collected over several years, by various persons for various purposes - is suitable regarding systematic conservation planning for species richness. The surrogate taxa concept applies in urban areas but with restrictions. Additionally, species richness should be examined in the context of both the city and its surrounding countryside.}, language = {en} } @article{BorchertWilkeSchmidtetal.2010, author = {Borchert, Manuela and Wilke, Max and Schmidt, Christian and Cauzid, Jean and Tucoulou, R{\´e}mi}, title = {Partitioning of Ba, La, Yb and Y between haplogranitic melts and aqueous solutions : an experimental study}, issn = {0009-2541}, doi = {10.1016/j.chemgeo.2010.06.009}, year = {2010}, abstract = {Barium, lanthanum, ytterbium, and yttrium partitioning experiments between fluid-saturated haplogranitic melts and aqueous solutions were conducted at 750 to 950 degrees C and 0.2 to 1 GPa to investigate the effects of melt and fluid composition, pressure, and temperature. Partition coefficients were determined using different experimental methods. On one hand quenched experiments were performed, and on the other hand, trace element contents in the aqueous fluid were determined directly using a hydrothermal diamond-anvil cell and synchrotron radiation X-ray fluorescence microanalysis of K-lines. The latter required a high excitation energy of 50 key due to the high energies necessary to excite the K-lines of the studied elements. The data from these two techniques showed good agreement for chloridic solutions, whereas quenching had a significant effect on results of the experiments with only water in the case of Ba. In Cl-free experiments, lanthanum and yttrium, trace element contents were even below detection limit in the quenched fluids, whereas small concentrations were detected in comparable in-situ experiments. This distinct difference is likely due to back reactions between fluid and melt upon cooling. The partitioning data of all elements show no dependence on the temperature and only small dependence on pressure. In contrast, the partitioning is strongly influenced by the composition of the starting fluid and melt. For chloridic fluids, there was a sharp increase in the Ba, La, Y and Yb partition coefficients with the alumina saturation index (ASI). The Ba partition coefficient increased from 0.002 at an ASI of 0.8 to 0.55 at an ASI of 1.07. At higher ASI, it decreased slightly to 0.2 at an ASI of similar to 1.3. Likewise, it was one to two orders of magnitude higher in chloridic fluids compared to those found in H2O experiments. Fluid-melt partition coefficients of La and Y increased from 0.002 at an ASI of similar to 0.8 to similar to 0.1 at an ASI of 1.2. In the same ASI range, the Yb partition coefficient increased to a maximum value of 0.02. Even at high salinities all elements fractionate into the melt. The compositional dependence of the partitioning data imply that both melt composition and fluid composition have a strong influence on trace element behavior and that complexation of Ba. REE and Y tin the fluid is not only controlled by the presence of Cl- in the fluid. Instead, interaction of these elements with major melt components dissolved in the fluid is very likely.}, language = {en} } @phdthesis{Borchert2010, author = {Borchert, Manuela}, title = {Interactions between aqueous fluids and silicate melts : equilibration, partitioning and complexation of trace elements}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-42088}, school = {Universit{\"a}t Potsdam}, year = {2010}, abstract = {The origin and evolution of granites has been widely studied because granitoid rocks constitute a major portion of the Earth's crust. The formation of granitic magma is, besides temperature mainly triggered by the water content of these rocks. The presence of water in magmas plays an important role due to the ability of aqueous fluids to change the chemical composition of the magma. The exsolution of aqueous fluids from melts is closely linked to a fractionation of elements between the two phases. Then, aqueous fluids migrate to shallower parts of the Earth's crust because of it's lower density compared to that of melts and adjacent rocks. This process separates fluids and melts, and furthermore, during the ascent, aqueous fluids can react with the adjacent rocks and alter their chemical signature. This is particularly impor- tant during the formation of magmatic-hydrothermal ore deposits or in the late stages of the evolution of magmatic complexes. For a deeper insight to these processes, it is essential to improve our knowledge on element behavior in such systems. In particular, trace elements are used for these studies and petrogenetic interpretations because, unlike major elements, they are not essential for the stability of the phases involved and often reflect magmatic processes with less ambiguity. However, for the majority of important trace elements, the dependence of the geochemical behavior on temperature, pressure, and in particular on the composition of the system are only incompletely or not at all experimentally studied. Former studies often fo- cus on the determination of fluid-melt partition coefficients (Df/m=cfluid/cmelt) of economically interesting elements, e.g., Mo, Sn, Cu, and there are some partitioning data available for ele- ments that are also commonly used for petrological interpretations. At present, no systematic experimental data on trace element behavior in fluid-melt systems as function of pressure, temperature, and chemical composition are available. Additionally, almost all existing data are based on the analysis of quenched phases. This results in substantial uncertainties, particularly for the quenched aqueous fluid because trace element concentrations may change upon cooling. The objective of this PhD thesis consisted in the study of fluid-melt partition coefficients between aqueous solutions and granitic melts for different trace elements (Rb, Sr, Ba, La, Y, and Yb) as a function of temperature, pressure, salinity of the fluid, composition of the melt, and experimental and analytical approach. The latter included the refinement of an existing method to measure trace element concentrations in fluids equilibrated with silicate melts di- rectly at elevated pressures and temperatures using a hydrothermal diamond-anvil cell and synchrotron radiation X-ray fluorescence microanalysis. The application of this in-situ method enables to avoid the main source of error in data from quench experiments, i.e., trace element concentration in the fluid. A comparison of the in-situ results to data of conventional quench experiments allows a critical evaluation of quench data from this study and literature data. In detail, starting materials consisted of a suite of trace element doped haplogranitic glasses with ASI varying between 0.8 and 1.4 and H2O or a chloridic solution with m NaCl/KCl=1 and different salinities (1.16 to 3.56 m (NaCl+KCl)). Experiments were performed at 750 to 950◦C and 0.2 or 0.5 GPa using conventional quench devices (externally and internally heated pressure vessels) with different quench rates, and at 750◦C and 0.2 to 1.4 GPa with in-situ analysis of the trace element concentration in the fluids. The fluid-melt partitioning data of all studied trace elements show 1. a preference for the melt (Df/m < 1) at all studied conditions, 2. one to two orders of magnitude higher Df/m using chloridic solutions compared to experiments with H2O, 3. a clear dependence on the melt composition for fluid-melt partitioning of Sr, Ba, La, Y, and Yb in experiments using chloridic solutions, 4. quench rate-related differences of fluid-melt partition coefficients of Rb and Sr, and 5. distinctly higher fluid-melt partitioning data obtained from in-situ experiments than from comparable quench runs, particularly in the case of H2O as starting solution. The data point to a preference of all studied trace elements for the melt even at fairly high salinities, which contrasts with other experimental studies, but is supported by data from studies of natural co-genetically trapped fluid and melt inclusions. The in-situ measurements of trace element concentrations in the fluid verify that aqueous fluids will change their composition upon cooling, which is in particular important for Cl free systems. The distinct differences of the in-situ results to quench data of this study as well as to data from the literature signify the im- portance of a careful fluid sampling and analysis. Therefore, the direct measurement of trace element contents in fluids equilibrated with silicate melts at elevated PT conditions represents an important development to obtain more reliable fluid-melt partition coefficients. For further improvement, both the aqueous fluid and the silicate melt need to be analyzed in-situ because partitioning data that are based on the direct measurement of the trace element content in the fluid and analysis of a quenched melt are still not completely free of quench effects. At present, all available data on element complexation in aqueous fluids in equilibrium with silicate melts at high PT are indirectly derived from partitioning data, which involves in these experiments assumptions on the species present in the fluid. However, the activities of chemical components in these partitioning experiments are not well constrained, which is required for the definition of exchange equilibria between melt and fluid species. For example, the melt-dependent variation of partition coefficient observed for Sr imply that this element can not only be complexed by Cl- as suggested previously. The data indicate a more complicated complexation of Sr in the aqueous fluid. To verify this hypothesis, the in-situ setup was also used to determine strontium complexation in fluids equilibrated with silicate melts at desired PT conditions by the application of X-ray absorption near edge structure (XANES) spectroscopy. First results show a strong effect of both fluid and melt composition on the resulting XANES spectra, which indicates different complexation environments for Sr.}, language = {en} } @phdthesis{Borchers2010, author = {Borchers, Andreas}, title = {Glaciomarine sedimentation at the continental margin of Prydz Bay, East Antarctica : implications on palaeoenvironmental changes during the Quaternary}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-52620}, school = {Universit{\"a}t Potsdam}, year = {2010}, abstract = {The Antarctic plays an important role in the global climate system. On the one hand, the Antarctic Ice Sheet is the largest freshwater reservoir on Earth. On the other hand, a major proportion of the global bottom-water formation takes place in Antarctic shelf regions, forcing the global thermohaline circulation. The main goal of this dissertation is to provide new insights into the dynamics and stability of the EAIS during the Quaternary. Additionally, variations in the activity of bottom-water formation and their causes are investigated. The dissertation is a German contribution to the International Polar Year 2007/ 2008 and was funded by the 'Deutsche Forschungsgesellschaft' (DFG) within the scope of priority program 1158 'Antarctic research with comparative studies in Arctic ice regions'. During RV Polarstern expedition ANT-XXIII/9, glaciomarine sediments were recovered from the Prydz Bay-Kerguelen region. Prydz Bay is a key region for the study of East EAIS dynamics, as 16\% of the EAIS are drained through the Lambert Glacier into the bay. Thereby, the glacier transports sediment into Prydz Bay which is then further distributed by calving icebergs or by current transport. The scientific approach of this dissertation is the reconstruction of past glaciomarine environments to infer on the response of the Lambert Glacier-Amery Ice Shelf system to climate shifts during the Quaternary. To characterize the depositional setting, sedimentological methods are used and statistical analyses are applied. Mineralogical and (bio)geochemical methods provide a means to reconstruct sediment provenances and to provide evidence on changes in the primary production in the surface water column. Age-depth models were constructed based on palaeomagnetic and palaeointensity measurements, diatom stratigraphy and radiocarbon dating. Sea-bed surface sediments in the investigation area show distinct variations in terms of their clay minerals and heavy-mineral assemblages. Considerable differences in the mineralogical composition of surface sediments are determined on the continental shelf. Clay minerals as well as heavy minerals provide useful parameters to differentiate between sediments which originated from erosion of crystalline rocks and sediments originating from Permo-Triassic deposits. Consequently, mineralogical parameters can be used to reconstruct the provenance of current-transported and ice-rafted material. The investigated sediment cores cover the time intervals of the last 1.4 Ma (continental slope) and the last 12.8 cal. ka BP (MacRobertson shelf). The sediment deposits were mainly influenced by glacial and oceanographic processes and further by biological activity (continental shelf), meltwater input and possibly gravitational transport. Sediments from the continental slope document two major deglacial events: the first deglaciation is associated with the mid-Pleistocene warming recognized around the Antarctic. In Prydz Bay, the Lambert Glacier-Amery Ice Shelf retreated far to the south and high biogenic productivity commenced or biogenic remains were better preserved due to increased sedimentation rates. Thereafter, stable glacial conditions continued until 400 - 500 ka BP. Calving of icebergs was restricted to the western part of the Lambert Glacier. The deeper bathymetry in this area allows for floating ice shelf even during times of decreased sea-level. Between 400 - 500 ka BP and the last interglacial (marine isotope stage 5) the glacier was more dynamic. During or shortly after the last interglacial the LAIS retreated again due to sea-level rise of 6 - 9 m. Both deglacial events correlate with a reduction in the thickness of ice masses in the Prince Charles Mountains. It indicates that a disintegration of the Amery Ice Shelf possibly led to increased drainage of ice masses from the Prydz Bay hinterland. A new end-member modelling algorithm was successfully applied on sediments from the MacRobertson shelf used to unmix the sand grain size fractions sorted by current activity and ice transport, respectively. Ice retreat on MacRobertson Shelf commenced 12.8 cal. ka BP and ended around 5.5 cal. ka BP. During the Holocene, strong fluctuations of the bottomwater activity were observed, probably related to variations of sea-ice formation in the Cape Darnley polynya. Increased activity of bottom-water flow was reconstructed at transitions from warm to cool conditions, whereas bottom-water activity receded during the mid- Holocene climate optimum. It can be concluded that the Lambert Glacier-Amery Ice Shelf system was relatively stable in terms of climate variations during the Quaternary. In contrast, bottom-water formation due to polynya activity was very sensitive to changes in atmospheric forcing and should gain more attention in future research.}, language = {en} } @article{BommerDouglasScherbaumetal.2010, author = {Bommer, Julian J. and Douglas, John and Scherbaum, Frank and Cotton, Fabrice and Bungum, Hilmar and Faeh, Donat}, title = {On the selection of ground-motion prediction equations for seismic hazard analysis}, issn = {0895-0695}, doi = {10.1785/gssrl.81.5.783}, year = {2010}, language = {en} } @article{BoenigerTronicke2010, author = {Boeniger, Urs and Tronicke, Jens}, title = {Improving the interpretability of 3D GPR data using target-specific attributes : application to tomb detection}, issn = {0305-4403}, doi = {10.1016/j.jas.2009.09.049}, year = {2010}, abstract = {Three-dimensional (3D) ground-penetrating radar (GPR) represents an efficient high-resolution geophysical surveying method allowing to explore archaeological sites in a non-destructive manner. To effectively analyze large 3D GPR data sets, their combination with modern visualization techniques (e.g., 3D isoamplitude displays) has been acknowledged to facilitate interpretation beyond classical time-slice analysis. In this study, we focus on the application of data attributes (namely energy, coherency, and similarity), originally developed for petroleum reservoir related problems addressed by reflection seismology, to emphasize temporal and spatial variations within GPR data cubes. Based on two case studies, we illustrate the potential of such attribute based analyses towards a more comprehensive 3D GPR data interpretation. The main goal of both case studies was to localize and potentially characterize tombs inside medieval chapels situated in the state of Brandenburg, Germany. By comparing the calculated data attributes to the conventionally processed data cubes, we demonstrate the superior interpretability of the coherency and the similarity attribute for target identification and characterization.}, language = {en} } @article{BoenigerTronicke2010, author = {Boeniger, Urs and Tronicke, Jens}, title = {Improving the interpretability of 3D GPR data using target-specific attributes : application to tomb detection (vol 37, pg 360, 2009)}, issn = {0305-4403}, doi = {10.1016/S0305-4403(10)00046-4}, year = {2010}, abstract = {Publisher's not}, language = {en} } @article{BoenigerTronicke2010, author = {Boeniger, Urs and Tronicke, Jens}, title = {Integrated data analysis at an archaeological site : a case study using 3D GPR, magnetic, and high-resolution topographic data}, issn = {0016-8033}, doi = {10.1190/1.3460432}, year = {2010}, abstract = {We have collected magnetic, 3D ground-penetrating-radar (GPR), and topographic data at an archaeological site within the Palace Garden of Paretz, Germany. The survey site covers an area of approximately 35 x 40 m across a hill structure (dips of up to 15 degrees) that is partly covered by trees. The primary goal of this study was to detect and locate the remains of ancient architectural elements, which, from historical records, were expected to be buried in the subsurface at this site. To acquire our geophysical data, we used a recently developed surveying approach that combines the magnetic and GPR instrument with a tracking total station (TTS). Besides efficient data acquisition, this approach provides positional information at an accuracy within the centimeter range. At the Paretz field site, this information was critical for processing and analyzing our geophysical data (in particular, GPR data) and enabled us to generate a high-resolution digital terrain model (DTM) of the surveyed area. Integrated analysis and interpretation based on composite images of the magnetic, 3D GPR, and high-resolution DTM data as well as selected attributes derived from these data sets allowed us to outline the remains of an artificial grotto and temple. Our work illustrates the benefit of using multiple surveying technologies, analyzing and interpreting the resulting data in an integrated fashion. It further demonstrates how modern surveying solutions allow for efficient, accurate data acquisition even in difficult terrain.}, language = {en} } @article{BoenigerTronicke2010, author = {Boeniger, Urs and Tronicke, Jens}, title = {On the potential of kinematic GPR surveying using a self-tracking total station : evaluating system crosstalk and latency}, issn = {0196-2892}, doi = {10.1109/Tgrs.2010.2048332}, year = {2010}, abstract = {In this paper, we present an efficient kinematic ground-penetrating radar (GPR) surveying setup using a self- tracking total station (TTS). This setup combines the ability of modern GPR systems to interface with Global Positioning System (GPS) and the capability of the employed TTS system to immediately make the positioning information available in a standardized GPS data format. Wireless communication between the GPR and the TTS system is established by using gain variable radio modems. Such a kinematic surveying setup faces two major potential limitations. First, possible crosstalk effects between the GPR and the positioning system have to be evaluated. Based on multiple walkaway experiments, we show that, for reasonable field setups, instrumental crosstalk has no significant impact on GPR data quality. Second, we investigate systematic latency (i.e., the time delay between the actual position measurement by TTS and its fusion with the GPR data) and its impact on the positional precision of kinematically acquired 2-D and 3-D GPR data. To quantify latency for our kinematic survey setup, we acquired forward-reverse profile pairs across a well-known subsurface target. Comparing the forward and reverse GPR images using three fidelity measures allows determining the optimum latency value and correcting for it. Accounting for both of these potential limitations allows us to kinematically acquire high- quality and high-precision GPR data using off-the-shelf instrumentation without further hardware modifications. Until now, these issues have not been investigated in detail, and thus, we believe that our findings have significant implications also for other geophysical surveying approaches.}, language = {en} } @article{BloecherCacaceLewerenzetal.2010, author = {Bloecher, Mando Guido and Cacace, Mauro and Lewerenz, Bjoern and Zimmermann, G{\"u}nter}, title = {Three dimensional modelling of fractured and faulted reservoirs : framework and implementation}, issn = {0009-2819}, doi = {10.1016/j.chemer.2010.05.014}, year = {2010}, abstract = {Modelling of coupled physical processes in fractured and faulted media is a major challenge for the geoscience community. Due to the complexity related to the geometry of real fracture networks and fault systems, modelling studies have been mainly restricted either to two dimensional cases or to simplified orthogonal fracture systems consisting of vertical and horizontal fractures. An approach to generate three dimensional meshes for realistic fault geometries is presented. The method enables representation of faults in an arbitrary incline as two dimensional planes within a three dimensional, stratified porous matrix of a generic geometry. Based on a structural geological model, the method creates three dimensional unstructured tetrahedral meshes. These meshes can be used for finite element and finite volume numerical simulations. A simulation of a coupled fluid flow and heat transport problem for a two layered porous medium cut by two crossing faults is presented to test the reliability of the method.}, language = {en} } @article{BlaserKruegerOhrnbergeretal.2010, author = {Blaser, Lilian and Kr{\"u}ger, Frank and Ohrnberger, Matthias and Scherbaum, Frank}, title = {Scaling relations of earthquake source parameter estimates with special focus on subduction environment}, issn = {0037-1106}, doi = {10.1785/0120100111}, year = {2010}, abstract = {Earthquake rupture length and width estimates are in demand in many seismological applications. Earthquake magnitude estimates are often available, whereas the geometrical extensions of the rupture fault mostly are lacking. Therefore, scaling relations are needed to derive length and width from magnitude. Most frequently used are the relationships of Wells and Coppersmith (1994) derived on the basis of a large dataset including all slip types with the exception of thrust faulting events in subduction environments. However, there are many applications dealing with earthquakes in subduction zones because of their high seismic and tsunamigenic potential. There are no well-established scaling relations for moment magnitude and length/width for subduction events. Within this study, we compiled a large database of source parameter estimates of 283 earthquakes. All focal mechanisms are represented, but special focus is set on (large) subduction zone events, in particular. Scaling relations were fitted with linear least-square as well as orthogonal regression and analyzed regarding the difference between continental and subduction zone/oceanic relationships. Additionally, the effect of technical progress in earthquake parameter estimation on scaling relations was tested as well as the influence of different fault mechanisms. For a given moment magnitude we found shorter but wider rupture areas of thrust events compared to Wells and Coppersmith (1994). The thrust event relationships for pure continental and pure subduction zone rupture areas were found to be almost identical. The scaling relations differ significantly for slip types. The exclusion of events prior to 1964 when the worldwide standard seismic network was established resulted in a remarkable effect on strike-slip scaling relations: the data do not show any saturation of rupture width of strike- slip earthquakes. Generally, rupture area seems to scale with mean slip independent of magnitude. The aspect ratio L/W, however, depends on moment and differs for each slip type.}, language = {en} } @article{BauerMoeckNordenetal.2010, author = {Bauer, Klaus and Moeck, Inga and Norden, Ben and Schulze, Alexander and Weber, Michael H. and Wirth, Holger}, title = {Tomographic P wave velocity and vertical velocity gradient structure across the geothermal site Groß Schoenebeck (NE German Basin) : relationship to lithology, salt tectonics, and thermal regime}, issn = {0148-0227}, doi = {10.1029/2009jb006895}, year = {2010}, abstract = {Seismic wide-angle data were collected along a 40-km-long profile centered at the geothermal research well GrSk 3/90 in the Northeast German Basin. Tomographic inversion of travel time data provided a velocity and a vertical velocity gradient model, indicative of Cenozoic to Pre-Permian sediments. Wide-angle reflections are modeled and interpreted as top Zechstein and top Pre-Permian. Changes in velocity gradients are interpreted as the transition from mechanical to chemical compaction at 2-3 km depth, and localized salt structures are imaged, suggesting a previously unknown salt pillow in the southern part of the seismic profile. The Zechstein salt shows decreased velocities in the adjacent salt pillows compared to the salt lows, which is confirmed by sonic log data. This decrease in velocity could be explained by the mobilization of less dense salt, which moved and formed the salt pillows, whereas the denser salt remained in place at the salt lows. We interpret a narrow subvertical low-velocity zone under the salt pillow at GrSk 3/ 90 as a fault in the deep Permian to Pre-Permian. This WNW-ESE trending fault influenced the location of the salt tectonics and led to the formation of a fault-bounded graben in the Rotliegend sandstones with optimal mechanical conditions for geothermal production. Thermal modeling showed that salt pillows are related to chimney effects, a decrease in temperature, and increasing velocity. The assumed variations in salt lithology, density, and strain must thus be even higher to compensate for the temperature effect.}, language = {en} } @article{AlAtikAbrahamsonBommeretal.2010, author = {Al Atik, Linda and Abrahamson, Norman A. and Bommer, Julian J. and Scherbaum, Frank and Cotton, Fabrice and Kuehn, Nicolas}, title = {The variability of ground-motion prediction models and its components}, issn = {0895-0695}, doi = {10.1785/gssrl.81.5.794}, year = {2010}, language = {en} } @article{AdhikariKallmeyer2010, author = {Adhikari, Rishi Ram and Kallmeyer, Jens}, title = {Detection and quantification of microbial activity in the subsurface}, issn = {0009-2819}, doi = {10.1016/j.chemer.2010.05.003}, year = {2010}, abstract = {The subsurface harbors a large fraction of Earth's living biomass, forming complex microbial ecosystems. Without a profound knowledge of the ongoing biologically mediated processes and their reaction to anthropogenic changes it is difficult to assess the long-term stability and feasibility of any type of geotechnical utilization, as these influence subsurface ecosystems. Despite recent advances in many areas of subsurface microbiology, the direct quantification of turnover processes is still in its infancy, mainly due to the extremely low cell abundances. We provide an overview of the currently available techniques for the quantification of microbial turnover processes and discuss their specific strengths and limitations. Most techniques employed so far have focused on specific processes, e.g. sulfate reduction or methanogenesis. Recent studies show that processes that were previously thought to exclude each other can occur simultaneously, albeit at very low rates. Without the identification of the respective processes it is impossible to quantify total microbial activity. Even in cases where all simultaneously occurring processes can be identified, the typically very low rates prevent quantification. In many cases a simple measure of total microbial activity would be a better and more robust measure than assays for several specific processes. Enzyme or molecular assays provide a more general approach as they target key metabolic compounds. Depending on the compound targeted a broader spectrum of microbial processes can be quantified. The two most promising compounds are ATP and hydrogenase, as both are ubiquitous in microbes. Technical constraints limit the applicability of currently available ATP-assays for subsurface samples. A recently developed hydrogenase radiotracer assay has the potential to become a key tool for the quantification of subsurface microbial activity.}, language = {en} }