Refine
Year of publication
- 2014 (228) (remove)
Document Type
- Article (196)
- Doctoral Thesis (19)
- Preprint (8)
- Review (4)
- Postprint (1)
Keywords
- Earthquake source observations (7)
- Holocene (6)
- Tibetan Plateau (4)
- Erosion (3)
- Himalaya (3)
- Seismicity and tectonics (3)
- Aleatory variability (2)
- Andes (2)
- Baseline shift (2)
- Body waves (2)
Institute
- Institut für Geowissenschaften (228) (remove)
We assessed tropical montane cloud forest (TMCF) sensitivity to natural disturbance by drought, fire, and dieback with a 7300-year-long paleorecord. We analyzed pollen assemblages, charcoal accumulation rates, and higher plant biomarker compounds (average chain length [ACL] of n-alkanes) in sediments from Wai 'anapanapa, a small lake near the upper forest limit and the mean trade wind inversion ('IWI) in Hawai`i. The paleorecord of ACL suggests increased drought frequency and a lower awl elevation from 2555-1323 cal yr B.P. and 606-334 cal yr B.P. Charcoal began to accumulate and a novel fire regime was initiated ca. 880 cal yr B.P., followed by a decreased fire return interval at ca. 550 cal yr B.P. Diebacks occurred at 2931, 2161, 1162, and 306 cal yr B.P., and two of these were independent of drought or fire. Pollen assemblages indicate that on average species composition changed only 2.8% per decade. These dynamics, though slight, were significantly associated with disturbance. The direction of species composition change varied with disturbance type. Drought was associated with significantly more vines and lianas; fire was associated with an increase in the tree fern Sadleria and indicators of open, disturbed landscapes at the expense of epiphytic ferns; whereas stand-scale dieback was associated with an increase in the tree fern Cibotium. Though this cloud forest was dynamic in response to past disturbance, it has recovered, suggesting a resilient TMCF with no evidence of state change in vegetation type (e.g., grassland or shrubland).
In hydrology, the storage-discharge relationship is a fundamental catchment property. Understanding what controls this relationship is at the core of catchment science. To date, there are no direct methods to measure water storage at catchment scales (10(1)-10(3)km(2)). In this study, we use direct measurements of terrestrial water storage dynamics by means of superconducting gravimetry in a small headwater catchment of the Regen River, Germany, to derive empirical storage-discharge relationships in nested catchments of increasing scale. Our results show that the local storage measurements are strongly related to streamflow dynamics at larger scales (> 100km(2); correlation coefficient=0.78-0.81), but at small scale, no such relationship exists (similar to 1km(2); correlation coefficients=-0.11). The geologic setting in the region can explain both the disconnection between local water storage and headwater runoff, and the connectivity between headwater storage and streams draining larger catchment areas. More research is required to understand what controls the form of the observed storage-discharge relationships at the catchment scale. This study demonstrates that high-precision gravimetry can provide new insights into the complex relationship between state and response of hydrological systems.
Crustal earthquake swarms are an expression of intensive cracking and rock damaging over periods of days, weeks or month in a small source region in the crust. They are caused by longer lasting stress changes in the source region. Often, the localized stressing of the crust is associated with fluid or gas migration, possibly in combination with pre-existing zones of weaknesses. However, verifying and quantifying localized fluid movement at depth remains difficult since the area affected is small and geophysical prospecting methods often cannot reach the required resolution.
We apply a simple and robust method to estimate the velocity ratio between compressional (P) and shear (S) waves (upsilon(P)/upsilon(S)-ratio) in the source region of an earthquake swarm. The upsilon(P)/upsilon(S)-ratio may be unusual small if the swarm is related to gas in a porous or fractured rock. The method uses arrival time difference between P and S waves observed at surface seismic stations, and the associated double differences between pairs of earthquakes. An advantage is that earthquake locations are not required and the method seems lesser dependent on unknown velocity variations in the crust outside the source region. It is, thus, suited for monitoring purposes.
Applications comprise three natural, mid-crustal (8-10 km) earthquake swarms between 1997 and 2008 from the NW-Bohemia swarm region. We resolve a strong temporal decrease of upsilon(P)/upsilon(S) before and during the main activity of the swarm, and a recovery of upsilon(P)/upsilon(S) to background levels at the end of the swarms. The anomalies are interpreted in terms of the Biot-Gassman equations, assuming the presence of oversaturated fluids degassing during the beginning phase of the swarm activity.
Flood hazard projections under climate change are typically derived by applying model chains consisting of the following elements: "emission scenario - global climate model - downscaling, possibly including bias correction hydrological model - flood frequency analysis". To date, this approach yields very uncertain results, due to the difficulties of global and regional climate models to represent precipitation. The implementation of such model chains requires major efforts, and their complexity is high.
We propose for the Mekong River an alternative approach which is based on a shortened model chain: "emission scenario - global climate model - non-stationary flood frequency model". The underlying idea is to use a link between the Western Pacific monsoon and local flood characteristics: the variance of the monsoon drives a non-stationary flood frequency model, yielding a direct estimate of flood probabilities. This approach bypasses the uncertain precipitation, since the monsoon variance is derived from large-scale wind fields which are better represented by climate models. The simplicity of the monsoon-flood link allows deriving large ensembles of flood projections under climate change. We conclude that this is a worthwhile, complementary approach to the typical model chains in catchments where a substantial link between climate and floods is found.
The Dansgaard-Oeschger oscillations and Heinrich events described in North Atlantic sediments and Greenland ice are expressed in the climate of the tropics, for example, as documented in Arabian Sea sediments. Given the strength of this teleconnection, we seek to reconstruct its range of environmental impacts. We present geochemical and sedimentological data from core SO130-289KL from the Indus submarine slope spanning the last similar to 80 kyr. Elemental and grain size analyses consistently indicate that interstadials are characterized by an increased contribution of fluvial suspension from the Indus River. In contrast, stadials are characterized by an increased contribution of aeolian dust from the Arabian Peninsula. Decadal-scale shifts at climate transitions, such as onsets of interstadials, were coeval with changes in productivity-related proxies. Heinrich events stand out as especially dry and dusty events, indicating a dramatically weakened Indian summer monsoon, potentially increased winter monsoon circulation, and increased aridity on the Arabian Peninsula. This finding is consistent with other paleoclimate evidence for continental aridity in the northern tropics during these events. Our results strengthen the evidence that circum-North Atlantic temperature variations translate to hydrological shifts in the tropics, with major impacts on regional environmental conditions such as rainfall, river discharge, aeolian dust transport, and ocean margin anoxia.
Grain-size distributions offer powerful proxies of past environmental conditions that are related to sediment sorting processes. However, they are often of multimodal character because sediments can get mixed during deposition. To facilitate the use of grain size as palaeoenvironmental proxy, this study aims to distinguish the main detrital processes that contribute to lacustrine sedimentation across the Tibetan Plateau using grain-size end-member modelling analysis. Between three and five robust grain-size end-member subpopulations were distinguished at different sites from similarly-likely end-member model runs. Their main modes were grouped and linked to common sediment transport and depositional processes that can be associated with contemporary Tibetan climate (precipitation patterns and lake ice phenology, gridded wind and shear stress data from the High Asia Reanalysis) and local catchment configurations. The coarse sands and clays with grain-size modes > 250 mu m and < 2 mu m were probably transported by fluvial processes. Aeolian sands (similar to 200 mu m) and coarse local dust (similar to 60 mu m), transported by saltation and in near-surface suspension clouds, are probably related to occasional westerly storms in winter and spring. Coarse regional dust with modes similar to 25 mu m may derive from near-by sources that keep in longer term suspension. The continuous background dust is differentiated into two robust end members (modes: 5-10 and 2-5 mu m) that may represent different sources, wind directions and/or sediment trapping dynamics from long-range, upper-level westerly and episodic northerly wind transport. According to this study grain-size end members of only fluvial origin contribute small amounts to mean Tibetan lake sedimentation (19 +/- 5%), whereas local to regional aeolian transport and background dust deposition dominate the clastic sedimentation in Tibetan lakes (contributions: 42 +/- 14% and 51 +/- 11%). However, fluvial and alluvial reworking of aeolian material from nearby slopes during summer seems to limit end-member interpretation and should be cross-checked with other proxy information. If not considered as a stand-alone proxy, a high transferability to other regions and sediment archives allows helpful reconstructions of past sedimentation history.
In this study, we suggest a novel approach for the retrieval of regional moment tensors for earthquakes with small to moderate magnitudes. The first modification is the combined inversion of broadband and short-period waveform data. The broadband waveforms are inverted in a frequency range suitable for surface waves, whereas for the short-period data a frequency range suitable for body waves is applied. The second modification is the use of first-motion body-wave polarities to select the most probable solution out of all solutions from inversion. To combine three different criteria for selecting the most probable solution (i.e., residual from inversion, double-couple content of solution, number of nonmatching first-motion body-wave polarities), the L2 norm is applied to the normalized parameters. We chose five earthquakes within the Alborz mountains, Iran, as a case study (3.1 <= M-w <= 4.1). In this area, several factors exacerbate the difficulty of performing inversion for moment tensors, for example, a heterogeneous station network and large azimuthal gaps. We have demonstrated that our approach supplies reliable moment tensors when inversion from broadband data alone fails. In one case, we successfully retrieved a stable solution from short-period waveform data alone. Thus, our approach enables successful determination of seismic moment tensors wherever a sparse network of broadband stations has thus far prevented it.
This article presents comparisons among the five ground-motion models described in other articles within this special issue, in terms of data selection criteria, characteristics of the models and predicted peak ground and response spectral accelerations. Comparisons are also made with predictions from the Next Generation Attenuation (NGA) models to which the models presented here have similarities (e.g. a common master database has been used) but also differences (e.g. some models in this issue are nonparametric). As a result of the differing data selection criteria and derivation techniques the predicted median ground motions show considerable differences (up to a factor of two for certain scenarios), particularly for magnitudes and distances close to or beyond the range of the available observations. The predicted influence of style-of-faulting shows much variation among models whereas site amplification factors are more similar, with peak amplification at around 1s. These differences are greater than those among predictions from the NGA models. The models for aleatory variability (sigma), however, are similar and suggest that ground-motion variability from this region is slightly higher than that predicted by the NGA models, based primarily on data from California and Taiwan.
Understanding the relationships between density and spatio-thermal variations at convergent plate boundaries is important for deciphering the present-day dynamics and evolution of subduction zones. In particular, the interaction between densification due to mineralogical phase transitions and slab pull forces is subject to ongoing investigations. We have developed a two-dimensional subduction zone model that is based on thermodynamic equilibrium assemblage calculations and includes the effects of melting processes on the density distribution in the lithosphere. Our model calculates the "metamorphic density" of rocks as a function of pressure, temperature and chemical composition in a subduction zone down to 250 km. We have used this model to show how the hydration, dehydration, partial melting and fractionation processes of rocks all influence the metamorphic density and greatly depend on the temperature field within the subduction system. These processes are largely neglected by other approaches that reproduce the density distribution within this complex tectonic setting. Our model demonstrates that the initiation of edogitization (i.e., when crustal rocks reach higher densities than the ambient mantle) of the slab is not the only significant process that makes the descending slab denser and generates the slab pull force. Instead, the densification of the lithospheric mantle of the sinking slab starts earlier than eclogitization and contributes significantly to slab pull in the early stages of subduction. Accordingly, the complex metamorphic structure of the slab and the mantle wedge has an important impact on the development of subduction zones. (C) 2014 Elsevier B.V. All rights reserved.
Shallow inland water bodies in Malawi continue to be threatened by various environmental challenges despite their importance to the fisheries industry. Due to the complex interaction between natural and anthropogenic disturbances, disentangling the effect of the two may be a complicated process. The littoral zone of most water bodies is important in environmental reconstructions including pollution and lake level monitoring. This study used a littoral zone, transect-based approach employing multi-proxy palaeolimnological techniques to reconstruct recent environmental change (ca. 100 yrs.) in Lake Malombe in the Malawi Rift, East Africa. The results of the study could inform fisheries management in Lake Malombe, which experienced a catastrophic decline in fish stocks. Results support documentary evidence for the complete desiccation of the lake less than 100 years ago. Subsequently, there is evidence for accelerated eutrophication in the recent past. In light of these results, it is concluded that transect sampling approaches rather than relying on single core measurements, and the need for careful consideration of the types of proxy, are significant considerations in palaeoenvironmental reconstructions.
Recurrence-plot-based recurrence networks are an approach used to analyze time series using a complex networks theory. In both approaches - recurrence plots and recurrence networks -, a threshold to identify recurrent states is required. The selection of the threshold is important in order to avoid bias of the recurrence network results. In this paper, we propose a novel method to choose a recurrence threshold adaptively. We show a comparison between the constant threshold and adaptive threshold cases to study period-chaos and even period-period transitions in the dynamics of a prototypical model system. This novel method is then used to identify climate transitions from a lake sediment record.
In wetlands or riparian areas, water withdrawal by plants with access to groundwater or the capillary fringe often causes diurnal groundwater fluctuations. Various approaches use the characteristics of these fluctuations for estimation of daily groundwater evapotranspiration rates. The objective of this paper was to review the available methods, compare them with measured evapotranspiration and assess their recharge assumptions. For this purpose, we employed data of 85 rain-free days of a weighable groundwater lysimeter situated at a grassland site in the Spreewald wetland in north-east Germany. Measurements of hourly recharge and daily evapotranspiration rates were used to assess the different approaches. Our results showed that a maximum of 50% of the day to day variance of the daily evapotranspiration rates could be explained by the approaches based on groundwater fluctuations. Simple and more complex methods performed similarly. For some of the approaches, there were indications that erroneous assumptions compensated each other (e.g., when overestimated recharge counteracted underestimated storage change). We found that the usage of longer time spans resulted in improved estimates of the daily recharge rates and that the estimates were further enhanced by including two night averages. When derived from fitting estimates of recharge or evapotranspiration with according measurements the specific yield, needed to convert changes in water level to water volumes, differed considerably among the methods (from 0.022 to 0.064). Thus, the specific yield can be seen as correction factor that compensates for inadequate process descriptions.
In March 2010, the project CoCoCo (incipient COntinent-COntinent COllision) recorded a 650 km long amphibian N-S wide-angle seismic profile, extending from the Eratosthenes Seamount (ESM) across Cyprus and southern Turkey to the Anatolian plateau. The aim of the project is to reveal the impact of the transition from subduction to continent-continent collision of the African plate with the Cyprus-Anatolian plate. A visual quality check, frequency analysis and filtering were applied to the seismic data and reveal a good data quality. Subsequent first break picking, finite-differences ray tracing and inversion of the offshore wide-angle data leads to a first-arrival tomographic model. This model reveals (1) P-wave velocities lower than 6.5 km/s in the crust, (2) a variable crustal thickness of about 28 - 37 km and (3) an upper crustal reflection at 5 km depth beneath the ESM. Two land shots on Turkey, also recorded on Cyprus, airgun shots south of Cyprus and geological and previous seismic investigations provide the information to derive a layered velocity model beneath the Anatolian plateau and for the ophiolite complex on Cyprus. The analysis of the reflections provides evidence for a north-dipping plate subducting beneath Cyprus. The main features of this layered velocity model are (1) an upper and lower crust with large lateral changes of the velocity structure and thickness, (2) a Moho depth of about 38 - 45 km beneath the Anatolian plateau, (3) a shallow north-dipping subducting plate below Cyprus with an increasing dip and (4) a typical ophiolite sequence on Cyprus with a total thickness of about 12 km. The offshore-onshore seismic data complete and improve the information about the velocity structure beneath Cyprus and the deeper part of the offshore tomographic model. Thus, the wide-angle seismic data provide detailed insights into the 2-D geometry and velocity structures of the uplifted and overriding Cyprus-Anatolian plate. Subsequent gravity modelling confirms and extends the crustal P-wave velocity model. The deeper part of the subducting plate is constrained by the gravity data and has a dip angle of ~ 28°. Finally, an integrated analysis of the geophysical and geological information allows a comprehensive interpretation of the crustal structure related to the collision process.
Making confident statements about the evolution of an ice-sheet shelf system with a numerical model requires the capability to reproduce the migration of the grounding line. Here we show that the shallow-ice approximation/shallow-shelf approximation hybrid-type Parallel Ice Sheet Model (PISM), with its recent improvements, is capable of modeling the grounding line motion in a perturbed ice-sheet shelf system. The model is set up according to the three-dimensional Marine Ice-Sheet Model Intercomparison Project (MISMIP3d), and simulations are carried out across a broad range of spatial resolutions. Using (1) a linear interpolation of the grounding line with locally interpolated basal friction and (2) an improved driving-stress computation across the grounding line, the reversibility of the grounding line (i.e. its retreat after an advance forced by a local perturbation of basal resistance) is captured by the model even at medium and low resolutions (Delta x > 10 km). The transient model response is qualitatively similar to that of higher-order models but reveals a higher initial sensitivity to perturbations on very short timescales. Our findings support the application of PISM to the Antarctic ice sheet from regional up to continental scales and on relatively low spatial resolutions.
Garnet brought to the surface by late Miocene granitoids at La Galite Archipelago (Central Mediterranean, Tunisia) contains abundant primary melt and fluid inclusions. Microstructural observations and mineral chemistry define the host garnet as a peritectic phase produced by biotite incongruent melting at ~800 degrees C and 0.5GPa, under fluid-present conditions. The trapped melt is leucogranitic with an unexpected metaluminous and almost peralkaline character. Fluid inclusions are one phase at room temperature, and contain a CO2-dominated fluid, with minor H2O, N-2 and CH4. Siderite and an OH-bearing phase were identified by Raman and IR spectroscopy within every analysed inclusion, and are interpreted as products of a post-entrapment carbonation/hydration reaction between the fluid and the host during cooling. The fluid present during anatexis is therefore inferred to have been originally richer in both H2O and CO2. The production of anatectic melt with a metaluminous signature can be explained as the result of partial melting of relatively Al-poor protoliths assisted by CO2-rich fluids.
The Isabena catchment (445 km(2)), Spain, features highly diverse spatial heterogeneity in land use, lithology and rainfall. Consequently, the relative contribution in terms of water and sediment yield varies immensely between its subcatchments, and also temporally. This study presents the synthesis of similar to 2.5 years of monitoring rainfall, discharge and suspended sediment concentration (SSC) in the five main subcatchments of the Isabena and its outlet.
Continuous discharge at the subcatchment outlets, nine tipping bucket rainfall and automatic SSC samplers (complemented by manual samples), were collected from June 2011 until November 2013. The water stage records were converted to discharge using a rating curve derived with Bayesian regression. For reconstructing sediment yields, the data from the intermittent SSC sampling needed to be interpolated. We employed non-parametric multivariate regression (Quantile Regression Forests, QRF) using the discharge and rainfall data plus different aggregation levels of these as ancillary predictors. The subsequent Monte Carlo simulations allowed the determination of monthly sediment yields and their uncertainty.
The Isabena catchment shows high erosion dynamics with great variability in space and time, with stark contrasts even between adjacent subcatchments. The natural conditions make water and sediment monitoring and instrumentation very challenging; the measurement of discharge is particularly prone to considerable uncertainties. The QRF method employed for reconstructing sedigraphs and monthly yields proved well suited for the task.
Permafrost-affected soils are among the most obvious ecosystems in which current microbial controls on organic matter decomposition are changing as a result of global warming. Warmer conditions in polygonal tundra will lead to a deepening of the seasonal active layer, provoking changes in microbial processes and possibly resulting in exacerbated carbon degradation under increasing anoxic conditions. To identify current microbial assemblages in carbon rich, water saturated permafrost environments, four polygonal tundra sites were investigated on Herschel Island and the Yukon Coast, Western Canadian Arctic. Ion Torrent sequencing of bacterial and archaeal 16S rRNA amplicons revealed the presence of all major microbial soil groups and indicated a local, vertical heterogeneity of the polygonal tundra soil community with increasing depth. Microbial diversity was found to be highest in the surface layers, decreasing towards the permafrost table. Quantitative PCR analysis of functional genes involved in carbon and nitrogen-cycling revealed a high functional potential in the surface layers, decreasing with increasing active layer depth. We observed that soil properties driving microbial diversity and functional potential varied in each study site. These results highlight the small-scale heterogeneity of geomorphologically comparable sites, greatly restricting generalizations about the fate of permafrost-affected environments in a warming Arctic.
After more than a decade of multidisciplinary studies of the Central American subduction zone mainly in the framework of two large research programmes, the US MARGINS program and the German Collaborative Research Center SFB 574, we here review and interpret the data pertinent to quantify the cycling of mineral-bound volatiles (H2O, CO2, Cl, S) through this subduction system. For input-flux calculations, we divide the Middle America Trench into four segments differing in convergence rate and slab lithological profiles, use the latest evidence for mantle serpentinization of the Cocos slab approaching the trench, and for the first time explicitly include subduction erosion of forearc basement. Resulting input fluxes are 40-62 (53) Tg/Ma/m H2O, 7.8-11.4 (9.3) Tg/Ma/m CO2, 1.3-1.9 (1.6) Tg/Ma/m Cl, and 1.3-2.1 (1.6) Tg/Ma/m S (bracketed are mean values for entire trench length). Output by cold seeps on the forearc amounts to 0.625-1.25 Tg/Ma/m H2O partly derived from the slab sediments as determined by geochemical analyses of fluids and carbonates. The major volatile output occurs at the Central American volcanic arc that is divided into ten arc segments by dextral strike-slip tectonics. Based on volcanic edifice and widespread tephra volumes as well as calculated parental magma masses needed to form observed evolved compositions, we determine long-term (10(5) years) average magma and K2O fluxes for each of the ten segments as 32-242 (106) Tg/Ma/m magma and 0.28-2.91 (1.38) Tg/Ma/m K2O (bracketed are mean values for entire Central American volcanic arc length). Volatile/K2O concentration ratios derived from melt inclusion analyses and petrologic modelling then allow to calculate volatile fluxes as 1.02-14.3 (6.2) Tg/Ma/m H2O, 0.02-0.45 (0.17) Tg/Ma/m CO2, and 0.07-0.34 (0.22) Tg/Ma/m Cl. The same approach yields long-term sulfur fluxes of 0.12-1.08 (0.54) Tg/Ma/m while present-day open-vent SO2-flux monitoring yields 0.06-2.37 (0.83) Tg/Ma/m S. Input-output comparisons show that the arc water fluxes only account for up to 40 % of the input even if we include an "invisible" plutonic component constrained by crustal growth. With 20-30 % of the H2O input transferred into the deeper mantle as suggested by petrologic modeling, there remains a deficiency of, say, 30-40 % in the water budget. At least some of this water is transferred into two upper-plate regions of low seismic velocity and electrical resistivity whose sizes vary along arc: one region widely envelopes the melt ascent paths from slab top to arc and the other extends obliquely from the slab below the forearc to below the arc. Whether these reservoirs are transient or steady remains unknown.