Refine
Year of publication
- 2009 (120) (remove)
Document Type
- Article (99)
- Doctoral Thesis (15)
- Other (3)
- Monograph/Edited Volume (2)
- Postprint (1)
Keywords
- Biomarker (2)
- Iran (2)
- Momententensor (2)
- Paläoklima (2)
- Seismotektonik (2)
- biomarker (2)
- moment tensor (2)
- seismotectonics (2)
- stable isotopes (2)
- Africa (1)
Institute
- Institut für Geowissenschaften (120) (remove)
Precision farming overcomes the paradigm of uniform field treatment by site-specific data acquisition and treatment to cope with within-field variability. Precision farming heavily relies on spatially dense information about soil and crop status. While it is often difficult and expensive to obtain precise soil information by traditional soil sampling and laboratory analysis some geophysical methods offer means to obtain subsidiary data in an efficient way. In particular, geoelectrical soil mapping has become widely accepted in precision farming. At present it is the most successful geophysical method providing the spatial distribution of relevant agronomic information that enables us to determine management zones for precision farming. Much work has been done to test the applicability of existing geoelectrical methods and to develop measurement systems applicable in the context of precision farming. Therefore, the aim of this paper was to introduce the basic ideas of precision farming, to discuss current precision farming applied geoelectrical methods and instruments and to give an overview about our corresponding activities during recent years. Different experiments were performed both in the laboratory and in the field to estimate first, electrical conductivity affecting factors, second, relationships between direct push and surface measurements, third, the seasonal stability of electrical conductivity patterns and fourth, the relationship between plant yield and electrical conductivity. From the results of these experiments, we concluded that soil texture is a very dominant factor in electrical conductivity mapping. Soil moisture affects both the level and the dynamic range of electrical conductivity readings. Nevertheless, electrical conductivity measurements can be principally performed independent of season. However, electrical conductivity field mapping does not produce reliable maps of spatial particle size distribution of soils, e.g., necessary to generate input parameters for water and nutrient transport models. The missing step to achieve this aim may be to develop multi-sensor systems that allow adjusting the electrical conductivity measurement from the influence of different soil water contents.
Single station seismic noise measurements were carried out at 192 sites in the western part of Istanbul, Turkey. This extensive survey allowed the fundamental resonance frequency of the sedimentary cover to be mapped, and identify areas prone to site amplification. The results are in good agreement with the geological distribution of sedimentary units, indicating a progressive decrease of the fundamental resonance frequencies from the northeastern part, where the bedrock outcrops, towards the southwestern side, where a thickness of some hundreds meters for the sedimentary cover is estimated. The particular distribution of fundamental resonance frequencies indicates that local amplification of the ground motion might play a significative role in explaining the anomalous damage distribution after the 17 August 1999 Kocaeli Earthquake. Furthermore, 2D array measurements of seismic noise were performed in the metropolitan area with the aim of obtaining a preliminary geophysical characterization of the different sedimentary covers. These measurements allow the estimation of the shear-wave velocity profile for some representative areas and the identification of the presence of strong impedance contrast responsible of seismic ground motion amplification. Comparison of a theoretical site response from an estimated S-wave velocity profile with an empirical one based on earthquake recordings strongly encourages the use of the low cost seismic noise techniques for the study of seismic site effects.
To study the applicability of the passive seismic interferometry technique to near-surface geological studies, seismic noise recordings from a small scale 2-D array of seismic stations were performed in the test site of Nauen (Germany). Rayleigh wave Green's functions were estimated for different frequencies. A tomographic inversion of the traveltimes estimated for each frequency from the Green's functions is then performed, allowing the laterally varying 3-D surfacewave velocity structure below the array to be retrieved at engineering-geotechnical scales. Furthermore, a 2-D S-wave velocity cross-section is obtained by combining 1-D velocity structures derived from the inversion of the dispersion curves extracted at several points along a profile where other geophysical analyses were performed. It is shown that the cross-section from passive seismic interferometry provides a clear image of the local structural heterogeneities that are in excellent agreement with georadar and geoelectrical results. Such findings indicate that the interferometry analysis of seismic noise is potentially of great interest for deriving the shallow 3-D velocity structure in urban areas.
In this work, we analyse continuous measurements of microseisms to assess the reliability of the fundamental resonance frequency estimated by means of the horizontal-to-vertical (H/V) spectral ratio within the 0.1-1 Hz frequency range, using short-period sensors (natural period of 1 s). We apply the H/V technique to recordings of stations installed in two alluvial basins with different sedimentary cover thicknesses-the Lower Rhine Embayment (Germany) and the Gubbio Plain (Central Italy). The spectral ratios are estimated over the time-frequency domain, and we discuss the reliability of the results considering both the variability of the microseism activity and the amplitude of the instrumental noise. We show that microseisms measured by short period sensors allow the retrieval of fundamental resonance frequencies greater than about 0.1-0.2 Hz, with this lower frequency bound depending on the relative amplitude of the microseism signal and the self-noise of the instruments. In particular, we show an example where the considered short-period sensor is connected to instruments characterized by an instrumental noise level which allows detecting only fundamental frequencies greater than about 0.4 Hz. Since the frequency at which the peak of the H/V spectral ratio is biased depends upon the seismic signal-to-instrument noise ratio, the power spectral amplitude of instrumental self- noise should be always considered when interpreting the frequency of the peak as the fundamental resonance frequency of the investigated site.
P>A vertical array of accelerometers was installed in Atakoy (western Istanbul) with the long-term aim of improving our understanding of in situ soil behaviour, to assess the modelling and parametric uncertainties associated with the employed methodologies for strong-motion site-response analysis, and for shallow geological investigations. Geotechnical and geophysical investigations were carried out to define the subsoil structure at the selected site. Data associated with 10 earthquakes (2.7 < M < 4.3) collected during the first months of operation of the array were used to image the upgoing and downgoing waves by deconvolution of waveforms recorded at different depths. Results have shown that the velocity of propagation of the imaged upgoing and downgoing waves in the borehole is consistent with that of S or P waves, depending on the component of ground acceleration analysed but independent of the chosen signal window. In particular, an excellent agreement was found between the observed upgoing and downgoing wave traveltimes and the ones calculated using a model derived by seismic noise analysis of array data. The presence of a smaller pulse on the waveforms obtained by deconvolution of the horizontal components suggests both internal S-wave reflection and S-to-P mode conversion, as well as a not normal incidence of the wavefield. The presence of a pulse propagating with S-wave velocity in the uppermost 25 m in the waveforms obtained by the deconvolution of the vertical components suggests P-to-S mode conversion. These evidences imply that, even when site amplification is mainly related to 1-D effects, the standard practice in engineering seismology of deconvolving the surface recording down to the bedrock using an approximate S-wave transfer function (generally valid for vertical incidence of SH waves) might lead to errors in the estimation of the input ground motion required in engineering calculations. Finally, downgoing waves with significant amplitudes were found down to 70 m and even to 140 m depth. This result provides a warning about the use of shallow borehole recordings as input for the numerical simulation of ground motion and for the derivation of ground motion prediction relationships.
In order to consistently approximate the thermal vertical structure of past upper water columns, Mg/Ca ratios of eight planktonic foraminiferal species with different preferential calcification depths selected from 76 tropical Atlantic and Caribbean sediment-surface samples were calibrated with delta O-18-derived calcification temperatures with an overall range of approximate to 8-28 degrees C. Extending the broad number of species-specific calibrations, which agree well especially with our shallow-dweller calibrations, this study presents new bulk calcite Mg/Ca vs. calcification temperature relationships for shallow-dwelling Globigerinoides ruber pink, thermocline-dwelling Globorotalia menardii, and deep-dwelling Globorotalia truncatulinoides dextral and Globorotalia crassaformis not separately calibrated before. The species-specific temperature sensitivities are relatively similar (approximate to 7- 11% increase in Mg/Ca per 1 degrees C), yet y-axis intercepts vary from 0.23-0.65 for the shallow and thermocline dwellers to 0.83-1.32 for the deep dwellers. Based on these differences, we established a 'warm water' calibration for temperatures > 19 degrees C (Mg/Ca=0.29.exp(0.101.T): r=0.90; shallow and thermocline dwellers) and a 'cold water' calibration for temperatures < 15 degrees C (Mg/Ca=0.84.exp(0.083.T); r=0.85; deep dwellers). These calibrations are offset by approximate to 8 degrees C. This maybe significant for paleotemperature reconstructions, which are afflicted with the problem that similar Mg/Ca offsets are probably characteristic of extinct species used to calculate past temperatures.
This study presents noble gas compositions (He, Ne, Ar, Kr, and Xe) of lavas from several Hawaiian volcanoes. Lavas from the Hawaii Scientific Drilling Project (HSDP) core, surface samples from Mauna Kea, Mauna Loa, Kilauea, Hualalai, Kohala and Haleakala as well as lavas from a deep well on the summit of Kilauea were investigated. Noble gases, especially helium, are used as tracers for mantle reservoirs, based on the assumption that high 3He/4He ratios (>8 RA) represent material from the deep and supposedly less degassed mantle, whereas lower ratios (~ 8 RA) are thought to represent the upper mantle. Shield stage Mauna Kea, Kohala and Kilauea lavas yielded MORB-like to moderately high 3He/4He ratios, while 3He/4He ratios in post-shield stage Haleakala lavas are MORB-like. Few samples show 20Ne/22Ne and 21Ne/22Ne ratios different from the atmospheric values, however, Mauna Kea and Kilauea lavas with excess in mantle Ne agree well with the Loihi-Kilauea line in a neon three-isotope plot, whereas one Kohala sample plots on the MORB correlation line. The values in the 4He/40Ar* (40Ar* denotes radiogenic Ar) versus 4He diagram imply open system fractionation of He from Ar, with a deficiency in 4He. Calculated 4He/40Ar*, 3He/22Nes (22NeS denotes solar Ne) and 4He/21Ne ratios for the sample suite are lower than the respective production and primordial ratios, supporting the observation of a fractionation of He from the heavier noble gases, with a depletion of He with respect to Ne and Ar. The depletion of He is interpreted to be partly due to solubility controlled gas loss during magma ascent. However, the preferential He loss suggests that He is more incompatible than Ne and Ar during magmatic processes. In a binary mixing model, the isotopic He and Ne pattern are best explained by a mixture of a MORB-like end-member with a plume like or primordial end-member with a fractionation in 3He/22Ne, represented by a curve parameter r of 15 (r=(³He/²²Ne)MORB/(³He/²²Ne)PLUME or PRIMORDIAL). Whether the high 3He/4He ratios in Hawaiian lavas are indicative of a primitive component within the Hawaiian plume or are rather a product of the crystal-melt- partitioning behavior during partial melting remains to be resolved.
The seismicity of the Dead Sea fault zone (DSFZ) during the last two millennia is characterized by a number of damaging and partly devastating earthquakes. These events pose a considerable seismic hazard and seismic risk to Syria, Lebanon, Palestine, Jordan, and Israel. The occurrence rates for large earthquakes along the DSFZ show indications to temporal changes in the long-term view. The aim of this thesis is to find out, if the occurrence rates of large earthquakes (Mw ≥ 6) in different parts of the DSFZ are time-dependent and how. The results are applied to probabilistic seismic hazard assessments (PSHA) in the DSFZ and neighboring areas. Therefore, four time-dependent statistical models (distributions), including Weibull, Gamma, Lognormal and Brownian Passage Time (BPT), are applied beside the exponential distribution (Poisson process) as the classical time-independent model. In order to make sure, if the earthquake occurrence rate follows a unimodal or a multimodal form, a nonparametric bootstrap test of multimodality has been done. A modified method of weighted Maximum Likelihood Estimation (MLE) is applied to estimate the parameters of the models. For the multimodal cases, an Expectation Maximization (EM) method is used in addition to the MLE method. The selection of the best model is done by two methods; the Bayesian Information Criterion (BIC) as well as a modified Kolmogorov-Smirnov goodness-of-fit test. Finally, the confidence intervals of the estimated parameters corresponding to the candidate models are calculated, using the bootstrap confidence sets. In this thesis, earthquakes with Mw ≥ 6 along the DSFZ, with a width of about 20 km and inside 29.5° ≤ latitude ≤ 37° are considered as the dataset. The completeness of this dataset is calculated since 300 A.D. The DSFZ has been divided into three sub zones; the southern, the central and the northern sub zone respectively. The central and the northern sub zones have been investigated but not the southern sub zone, because of the lack of sufficient data. The results of the thesis for the central part of the DSFZ show that the earthquake occurrence rate does not significantly pursue a multimodal form. There is also no considerable difference between the time-dependent and time-independent models. Since the time-independent model is easier to interpret, the earthquake occurrence rate in this sub zone has been estimated under the exponential distribution assumption (Poisson process) and will be considered as time-independent with the amount of 9.72 * 10-3 events/year. The northern part of the DSFZ is a special case, where the last earthquake has occurred in 1872 (about 137 years ago). However, the mean recurrence time of Mw ≥ 6 events in this area is about 51 years. Moreover, about 96 percent of the observed earthquake inter-event times (the time between two successive earthquakes) in the dataset regarding to this sub zone are smaller than 137 years. Therefore, it is a zone with an overdue earthquake. The results for this sub zone verify that the earthquake occurrence rate is strongly time-dependent, especially shortly after an earthquake occurrence. A bimodal Weibull-Weibull model has been selected as the best fit for this sub zone. The earthquake occurrence rate, corresponding to the selected model, is a smooth function of time and reveals two clusters within the time after an earthquake occurrence. The first cluster begins right after an earthquake occurrence, lasts about 80 years, and is explicitly time-dependent. The occurrence rate, regarding to this cluster, is considerably lower right after an earthquake occurrence, increases strongly during the following ten years and reaches its maximum about 0.024 events/year, then decreases over the next 70 years to its minimum about 0.0145 events/year. The second cluster begins 80 years after an earthquake occurrence and lasts until the next earthquake occurs. The earthquake occurrence rate, corresponding to this cluster, increases extremely slowly, such as it can be considered as an almost constant rate about 0.015 events/year. The results are applied to calculate the time-dependent PSHA in the northern part of the DSFZ and neighbouring areas.
This thesis presents investigations on sediments from two African lakes which have been recording changes in their surrounding environmental and climate conditions since more than 200,000 years. Focus of this work is the time of the last Glacial and the Holocene (the last ~100,000 years before present [in the following 100 kyr BP]). One important precondition for this kind of research is a good understanding of the present ecosystems in and around the lakes and of the sediment formation under modern climate conditions. Both studies therefore include investigations on the modern environment (including organisms, soils, rocks, lake water and sediments). A 90 m long sediment sequence was investigated from Lake Tswaing (north-eastern South Africa) using geochemical analyses. These investigations document alternating periods of high detrital input and low (especially autochthonous) organic matter content and periods of low detrital input, carbonatic or evaporitic sedimentation and high autochthonous organic matter content. These alternations are interpreted as changes between relatively humid and arid conditions, respectively. Before c. 75 kyr BP, they seem to follow changes in local insolation whereas afterwards they appear to be acyclic and are probably caused by changes in ocean circulation and/or in the mean position of the Inter-Tropical Convergence Zone (ITCZ). Today, these factors have main influence on precipitation in this area where rainfall occurs almost exclusively during austral summer. All modern organisms were analysed for their biomarker and bulk organic and compound-specific stable carbon isotope composition. The same investigations on sediments from the modern lake floor document the mixed input of the investigated individual organisms and reveal additional influences by methanotrophic bacteria. A comparison of modern sediment characteristics with those of sediments covering the time 14 to 2 kyr BP shows changes in the productivity of the lake and the surrounding vegetation which are best explained by changes in hydrology. More humid conditions are indicated for times older than 10 kyr BP and younger than 7.5 kyr BP, whereas arid conditions prevailed in between. These observations agree with the results from sediment composition and indications from other climate archives nearby. The second lake study deals with Lake Challa, a small, deep crater lake on the foot of Mount Kilimanjaro. In this lake form mm-scale laminated sediments which were analyses with micro-XRF scanning for changes in the element composition. By comparing these results with investigations on thin sections, results from ongoing sediment trap studies, meteorological data, and investigations on the surrounding rocks and soils, I develop a model for seasonal variability in the limnology and sedimentation of Lake Challa. The lake appears to be stratified during the warm rain seasons (October – December and March – May) during which detrital material is delivered to the lake and carbonates precipitate. On the lake floor forms a dark lamina with high contents of Fe and Ti and high Ca/Al and low Mn/Fe ratios. Diatoms bloom during the cool and windy season (June – September) when mixing down to c. 60 m depth provides easily bio-available nutrients. Contemporaneously, Fe and Mn-oxides are precipitating which cause high Mn/Fe ratios in the light diatom-rich laminae of the sediments. Trends in the Mn/Fe ratio of the sediments are interpreted to reflect changes in the intensity or duration of seasonal mixing in Lake Challa. This interpretation is supported by parallel changes in the organic matter and biogenic silica content observed in the 22 m long profile recovered from Lake Challa. This covers the time of the last 25 kyr BP. It documents a transition around 16 kyr BP from relatively well-mixed conditions with high detrital input during glacial times to stronger stratified conditions which are probably related to increasing lake levels in Challa and generally more humid conditions in East Africa. Intensified mixing is recorded for the time of the Younger Dryas and the period between 11.4 and 10.7 kyr BP. For these periods, reduced intensity of the SW monsoon and intensified NE monsoon are reported from archives of the Indian-Asian Monsoon region, arguing for the latter as a probable source for wind mixing in Lake Challa. This connection is probably also responsible for contemporaneous events in the Mn/Fe ratios of the Lake Challa sediments and in other records of northern hemisphere monsoon intensity during the Holocene and underlines the close interaction of global low latitude atmospheric circulation.
Submarine landslides can generate local tsunamis posing a hazard to human lives and coastal facilities. Two major related problems are: (i) quantitative estimation of tsunami hazard and (ii) early detection of the most dangerous landslides. This thesis focuses on both those issues by providing numerical modeling of landslide-induced tsunamis and by suggesting and justifying a new method for fast detection of tsunamigenic landslides by means of tiltmeters. Due to the proximity to the Sunda subduction zone, Indonesian coasts are prone to earthquake, but also landslide tsunamis. The aim of the GITEWS-project (German-Indonesian Tsunami Early Warning System) is to provide fast and reliable tsunami warnings, but also to deepen the knowledge about tsunami hazards. New bathymetric data at the Sunda Arc provide the opportunity to evaluate the hazard potential of landslide tsunamis for the adjacent Indonesian islands. I present nine large mass movements in proximity to Sumatra, Java, Sumbawa and Sumba, whereof the largest event displaced 20 km³ of sediments. Using numerical modeling, I compute the generated tsunami of each event, its propagation and runup at the coast. Moreover, I investigate the age of the largest slope failures by relating them to the Great 1977 Sumba earthquake. Continental slopes off northwest Europe are well known for their history of huge underwater landslides. The current geological situation west of Spitsbergen is comparable to the continental margin off Norway after the last glaciation, when the large tsunamigenic Storegga slide took place. The influence of Arctic warming on the stability of the Svalbard glacial margin is discussed. Based on new geophysical data, I present four possible landslide scenarios and compute the generated tsunamis. Waves of 6 m height would be capable of reaching northwest Europe threatening coastal areas. I present a novel technique to detect large submarine landslides using an array of tiltmeters, as a possible tool in future tsunami early warning systems. The dislocation of a large amount of sediment during a landslide produces a permanent elastic response of the earth. I analyze this response with a mathematical model and calculate the theoretical tilt signal. Applications to the hypothetical Spitsbergen event and the historical Storegga slide show tilt signals exceeding 1000 nrad. The amplitude of landslide tsunamis is controlled by the product of slide volume and maximal velocity (slide tsunamigenic potential). I introduce an inversion routine that provides slide location and tsunamigenic potential, based on tiltmeter measurements. The accuracy of the inversion and of the estimated tsunami height near the coast depends on the noise level of tiltmeter measurements, the distance of tiltmeters from the slide, and the slide tsunamigenic potential. Finally, I estimate the applicability scope of this method by employing it to known landslide events worldwide.
The Tibetan Plateau is the largest elevated landmass in the world and profoundly influences atmospheric circulation patterns such as the Asian monsoon system. Therefore this area has been increasingly in focus of palaeoenvironmental studies. This thesis evaluates the applicability of organic biomarkers for palaeolimnological purposes on the Tibetan Plateau with a focus on aquatic macrophyte-derived biomarkers. Submerged aquatic macrophytes have to be considered to significantly influence the sediment organic matter due to their high abundance in many Tibetan lakes. They can show highly 13C-enriched biomass because of their carbon metabolism and it is therefore crucial for the interpretation of δ13C values in sediment cores to understand to which extent aquatic macrophytes contribute to the isotopic signal of the sediments in Tibetan lakes and in which way variations can be explained in a palaeolimnological context. Additionally, the high abundance of macrophytes makes them interesting as potential recorders of lake water δD. Hydrogen isotope analysis of biomarkers is a rapidly evolving field to reconstruct past hydrological conditions and therefore of special relevance on the Tibetan Plateau due to the direct linkage between variations of monsoon intensity and changes in regional precipitation / evaporation balances. A set of surface sediment and aquatic macrophyte samples from the central and eastern Tibetan Plateau was analysed for composition as well as carbon and hydrogen isotopes of n-alkanes. It was shown how variable δ13C values of bulk organic matter and leaf lipids can be in submerged macrophytes even of a single species and how strongly these parameters are affected by them in corresponding sediments. The estimated contribution of the macrophytes by means of a binary isotopic model was calculated to be up to 60% (mean: 40%) to total organic carbon and up to 100% (mean: 66%) to mid-chain n-alkanes. Hydrogen isotopes of n-alkanes turned out to record δD of meteoric water of the summer precipitation. The apparent enrichment factor between water and n-alkanes was in range of previously reported ones (≈-130‰) at the most humid sites, but smaller (average: -86‰) at sites with a negative moisture budget. This indicates an influence of evaporation and evapotranspiration on δD of source water for aquatic and terrestrial plants. The offset between δD of mid- and long-chain n-alkanes was close to zero in most of the samples, suggesting that lake water as well as soil and leaf water are affected to a similar extent by those effects. To apply biomarkers in a palaeolimnological context, the aliphatic biomarker fraction of a sediment core from Lake Koucha (34.0° N; 97.2° E; eastern Tibetan Plateau) was analysed for concentrations, δ13C and δD values of compounds. Before ca. 8 cal ka BP, the lake was dominated by aquatic macrophyte-derived mid-chain n-alkanes, while after 6 cal ka BP high concentrations of a C20 highly branched isoprenoid compound indicate a predominance of phytoplankton. Those two principally different states of the lake were linked by a transition period with high abundances of microbial biomarkers. δ13C values were relatively constant for long-chain n-alkanes, while mid-chain n-alkanes showed variations between -23.5 to -12.6‰. Highest values were observed for the assumed period of maximum macrophyte growth during the late glacial and for the phytoplankton maximum during the middle and late Holocene. Therefore, the enriched values were interpreted to be caused by carbon limitation which in turn was induced by high macrophyte and primary productivity, respectively. Hydrogen isotope signatures of mid-chain n-alkanes have been shown to be able to track a previously deduced episode of reduced moisture availability between ca. 10 and 7 cal ka BP, indicated by a 20‰ shift towards higher δD values. Indications for cooler episodes at 6.0, 3.1 and 1.8 cal ka BP were gained from drops of biomarker concentrations, especially microbial-derived hopanoids, and from coincidental shifts towards lower δ13C values. Those episodes correspond well with cool events reported from other locations on the Tibetan Plateau as well as in the Northern Hemisphere. To conclude, the study of recent sediments and plants improved the understanding of factors affecting the composition and isotopic signatures of aliphatic biomarkers in sediments. Concentrations and isotopic signatures of the biomarkers in Lake Koucha could be interpreted in a palaeolimnological context and contribute to the knowledge about the history of the lake. Aquatic macrophyte-derived mid-chain n-alkanes were especially useful, due to their high abundance in many Tibetan Lakes and their ability to record major changes of lake productivity and palaeo-hydrological conditions. Therefore, they have the potential to contribute to a fuller understanding of past climate variability in this key region for atmospheric circulation systems.
With the increasing availability of observational data from different sources at a global level, joint analysis of these data is becoming especially attractive. For such an analysis – oftentimes with little prior knowledge about local and global interactions between the different observational variables at hand – an exploratory, data-driven analysis of the data may be of particular relevance. In the present work we used generalized additive models (GAM) in an exemplary study of spatio-temporal patterns in the tropospheric NO2-distribution derived from GOME satellite observations (1996 to 2001) at global scale. We focused on identifying correlations between NO2 and local wind fields, a quantity which is of particular interest in the analysis of spatio-temporal interactions. Formulating general functional, parametric relationships between the observed NO2 distribution and local wind fields, however, is difficult – if not impossible. So, rather than following a modelbased analysis testing the data for predefined hypotheses (assuming, for example, sinusoidal seasonal trends), we used a GAM with non-parametric model terms to learn this functional relationship between NO2 and wind directly from the data. The NO2 observations showed to be affected by winddominated processes over large areas. We estimated the extent of areas affected by specific NO2 emission sources, and were able to highlight likely atmospheric transport “pathways”. General temporal trends which were also part of our model – weekly, seasonal and linear changes – showed to be in good agreement with previous studies and alternative ways of analysing the time series. Overall, using a non-parametric model provided favorable means for a rapid inspection of this large spatio-temporal NO2 data set, with less bias than parametric approaches, and allowing to visualize dynamical processes of the NO2 distribution at a global scale.
The Chaco foreland basin was initiated during the late Oligocene as a result of thrusting in the Eastern Cordillera in response to Nazca-South America plate convergence. Foreland basins are the result of the flexural isostatic response of an elastic plate to orogenic and/or thrust sheet loading. We carried out flexural modelling along a W-E profile (21.4 degrees S) to investigate Chaco foreland basin development using new information on ages of foreland basin strata, elastic and sedimentary thicknesses and structural histories. It was possible to reproduce present-day elevation, gravity anomaly, Moho depth, elastic thicknesses, foreland sedimentary thicknesses and the basin geometry. Our model predicted the basin geometry and sedimentary thicknesses for different evolutionary stages. Measured thicknesses and previously proposed depozones were compared with our predictions. Our results shed more light on the Chaco foreland basin evolution and suggest that an apparent decrease in elastic thickness beneath the Eastern Cordillera and the Interandean Zone could have occurred between 14 and 6 Ma.
Crustal structure of the southern Dead Sea basin derived from project DESIRE wide-angle seismic data
(2009)
As part of the DEad Sea Integrated REsearch project (DESIRE) a 235 km long seismic wide-angle reflection/ refraction (WRR) profile was completed in spring 2006 across the Dead Sea Transform (DST) in the region of the southern Dead Sea basin (DSB). The DST with a total of about 107 km multi-stage left-lateral shear since about 18 Ma ago, accommodates the movement between the Arabian and African plates. It connects the spreading centre in the Red Sea with the Taurus collision zone in Turkey over a length of about 1 100 km. With a sedimentary infill of about 10 km in places, the southern DSB is the largest pull-apart basin along the DST and one of the largest pull-apart basins on Earth. The WRR measurements comprised 11 shots recorded by 200 three-component and 400 one-component instruments spaced 300 m to 1.2 km apart along the whole length of the E-W trending profile. Models of the P-wave velocity structure derived from the WRR data show that the sedimentary infill associated with the formation of the southern DSB is about 8.5 km thick beneath the profile. With around an additional 2 km of older sediments, the depth to the seismic basement beneath the southern DSB is about 11 km below sea level beneath the profile. Seismic refraction data from an earlier experiment suggest that the seismic basement continues to deepen to a maximum depth of about 14 km, about 10 km south of the DESIRE profile. In contrast, the interfaces below about 20 km depth, including the top of the lower crust and the Moho, probably show less than 3 km variation in depth beneath the profile as it crosses the southern DSB. Thus the Dead Sea pull-apart basin may be essentially an upper crustal feature with upper crustal extension associated with the left- lateral motion along the DST. The boundary between the upper and lower crust at about 20 km depth might act as a decoupling zone. Below this boundary the two plates move past each other in what is essentially a shearing motion. Thermo-mechanical modelling of the DSB supports such a scenario. As the DESIRE seismic profile crosses the DST about 100 km north of where the DESERT seismic profile crosses the DST, it has been possible to construct a crustal cross-section of the region before the 107 km left-lateral shear on the DST occurred.
Interception measurements and assessment of Gash model performance for a tropical semi-arid region
(2009)
Semi-arid environments usually face water scarcity and conflicts for its use; therefore a complete understanding of the water balance in these regions is desired. To evaluate interception, measurements of precipitation, throughfall and stemflow were carried out in a Brazilian tropical semi-arid experimental watershed with well preserved Caatinga vegetation. Data analysis indicates that interception losses correspond to 13% of total rainfall, representing an important process in the watershed's water balance, where runoff is only 6% of total precipitation. Gash interception model was applied in the region with good results for long term simulation. Nevertheless, the model produced significant but not systematic errors on a daily basis. This was attributed to its incapability of representing the temporal variation of precipitation during the event, which is a major factor affecting interception. Rainfall intensity was shown to be a good parameter to determine an applicability threshold for Gash model in the study area.
The Takab calcareous rocks of northwest Iran crop out in association with a variety of metamorphic rocks including mafic granulites, amphibolites, granitic gneisses, pelitic schists and meta-ultramafic rocks. They can be divided into marbles and calc-silicate rocks on the basis of the dominance of calcite/dolomite and silicate minerals. Dominant peak metamorphic granulite facies assemblage of calc-silicate rocks is Scp + Grt(I) + Cpx + Cal + Qtz +/- Hbl(I). The decrease of temperature and pressure during exhumation produced post-peak metamorphic assemblages. Coronal garnet (Grt II) in the calc-silicate rocks was produced by retrograde reactions consuming plagioclase and clinopyroxene, while peak metamorphic garnet (Grt I) occurs as preserved xenoblastic grains in calcite and/or plagioclase (Pl II). Regional metamorphism took place at 740 degrees C and X-CO2 similar to 0.9. Garnet-clinopyroxene-plagioclase-quartz (GADS) barometry yields a pressure of 8-9 kbar, corresponding to a depth of ca. 24-27 km. This was followed by decompression and hydration during exhumation of the crustal rocks up to the surface. Secondary phases such as garnet (II) hornblende (II), plagioclase (II), zoisite and titanite (II) constrain the temperature and pressure of post-peak metamorphism as similar to 600 degrees C and similar to 6 kbar respectively and a fluid with XCO2 as low as 0.4. Halogens were near-absent during the peak metamorphic stage. The scapolite and hornblende crystallized underpeak metamorphic conditions contain very low fluorine and chlorine, whereas relatively high fluorine (similar to 0.8 wt%) in the titanite (II) and hornblende (II) suggests a possible infiltration of F-rich fluids into the calc-silicate rocks during retrogression. It is interpreted to be related to external fluids released during crystallisation of granitoid magmas and/or leucosome patches in the adjacent migmatites.
Precision farming needs management rules to apply spatially differentiated treatments in agricultural fields. Digital soil mapping (DSM) tools, for example apparent soil electrical conductivity, corrected to 25A degrees C (EC25), and digital elevation models, try to explain the spatial variation in soil type, soil properties (e.g. clay content), site and crop that are determined by landscape characteristics such as terrain, geology and geomorphology. We examined the use of EC25 maps to delineate management zones, and identified the main factors affecting the spatial pattern of EC25 at the regional scale in a study area in eastern Germany. Data of different types were compared: EC25 maps for 11 fields, soil properties measured in the laboratory, terrain attributes, geological maps and the description of 75 soil profiles. We identified the factors that influence EC25 in the presence of spatial autocorrelation and field-specific random effects with spatial linear mixed-effects models. The variation in EC25 could be explained to a large degree (R (2) of up to 61%). Primarily, soil organic matter and CaCO3, and secondarily clay and the presence of gleyic horizons were significantly related to EC25. Terrain attributes, however, had no significant effect on EC25. The geological map unit showed a significant relationship to EC25, and it was possible to determine the most important soil properties affecting EC25 by interpreting the geological maps. Including information on geology in precision agriculture could improve understanding of EC25 maps. The EC25 maps of fields should not be assumed to represent a map of clay content to form a basis for deriving management zones because other factors appeared to have a more important effect on EC25.
Based on Proxy records from western Black Sea cores, we provide a comprehensive Study of climate change during the last glacial maximum and late-glacial period in the Black Sea region. For the first time we present a record of relative changes in precipitation for NW Anatolia based on variations in the terrigenous supply expressed as detrital carbonate concentration. The good correspondence between reconstructed rainfall intensity in NW Anatolia and past western Mediterranean sea Surface temperatures (SSTs) implies that during the glacial period the precipitation variability was controlled, like today, by Mediterranean cyclonic disturbances. Periods of reduced precipitation correlate well with low SSTs in the Mediterranean related to Heinrich events H1 and H2. Stable oxygen isotopes and lithological and mineralogical data point to a significant modification in the dominant freshwater/sediment source concomitant to the meltwater inflow after 16.4 cal ka BP. This change implies intensification of the northern sediment source and, with other records from the Mediterranean region, consistently suggests a reorganization of the atmospheric circulation pattern affecting the hydrology of the European continent. The early deglacial northward retreat of both atmospheric and oceanic polar fronts was responsible for the warming in the Mediterranean region, leading simultaneously to more humid conditions in central and northern Europe.
We analyze the role of megathrust geometry on slip estimation using the 1960 Chile earthquake (M-W = 9.5) as an example. A variable slip distribution for this earthquake has been derived by Barrientos and Ward (1990) applying an elastic dislocation model with a planar fault geometry. Their model shows slip patches at 80-110 km depth, isolated from the seismogenic zone, interpreted as aseismic slip. We invert the same geodetic data set using a finite element model (FEM) with precise geometry derived from geophysical data. Isoparametric FEM is implemented to constrain the slip distribution of curve-shaped elements. Slip resolved by our precise geometry model is limited to the shallow region of the plate interface suggesting that the deep patches of moment were most likely an artifact of the planar geometry. Our study emphasizes the importance of fault geometry on slip estimation of large earthquakes.
Empirical ground-motion models used in seismic hazard analysis are commonly derived by regression of observed ground motions against a chosen set of predictor variables. Commonly, the model building process is based on residual analysis and/or expert knowledge and/or opinion, while the quality of the model is assessed by the goodness-of-fit to the data. Such an approach, however, bears no immediate relation to the predictive power of the model and with increasing complexity of the models is increasingly susceptible to the danger of overfitting. Here, a different, primarily data-driven method for the development of ground-motion models is proposed that makes use of the notion of generalization error to counteract the problem of overfitting. Generalization error directly estimates the average prediction error on data not used for the model generation and, thus, is a good criterion to assess the predictive capabilities of a model. The approach taken here makes only few a priori assumptions. At first, peak ground acceleration and response spectrum values are modeled by flexible, nonphysical functions (polynomials) of the predictor variables. The inclusion of a particular predictor and the order of the polynomials are based on minimizing generalization error. The approach is illustrated for the next generation of ground-motion attenuation dataset. The resulting model is rather complex, comprising 48 parameters, but has considerably lower generalization error than functional forms commonly used in ground-motion models. The model parameters have no physical meaning, but a visual interpretation is possible and can reveal relevant characteristics of the data, for example, the Moho bounce in the distance scaling. In a second step, the regression model is approximated by an equivalent stochastic model, making it physically interpretable. The resulting resolvable stochastic model parameters are comparable to published models for western North America. In general, for large datasets generalization error minimization provides a viable method for the development of empirical ground-motion models.
The Middle Spotted Woodpecker (Dendrocopos medius) is the bird species which Germany has the greatest global responsibility to protect. It is an umbrella species for the entire assemblage of animals associated with mature broadleaved trees, especially oak. Even though well studied in small to medium scale stands, the validity of habitat suitability analysis for this species in larger forests has not previously been proved. Aim of this study was to test suitability of permanent forest inventory plots for modelling its distribution in a 17,000 ha forest landscape and to derive habitat threshold values as a basis for formulating management guidelines. Based on 150 randomly selected 12.5 ha plots we identified mean age and basal area of oaks as the most important habitat factors using a backward selection logistic model. Internal validation showed an AUC of 0.89 and a R-2(N) of 0.58. Determination of thresholds using maximally selected rank statistics found higher probability of occurrence in stands with a mean age >95 years. Above that age the probability increased again in stands with more than 6.4 m(2) basal area oak/ha. Our results show that widely available forest inventory data can serve as a valuable basis for monitoring the Middle Spotted Woodpecker, either within the framework of the Natura 2000 Network, or more generally in integrated forest management with the aim of providing suitable habitats for the entire assemblage of species on old deciduous trees, especially oak.
The Aue-Schwarzenberg Granite Zone (ASGZ), in the western Erzgebirge of Germany, is composed of small, late- Variscan F-poor biotite and two-mica granites. The biotite granites (Aue granite suite, Beierfeld, Bernsbach) are weakly to mildly peraluminous (A/CNK = 1.07-1.14; 70-76 wt% SiO2), display similar Sr-87/Sr-86 initial ratios (0.7065-0.7077; t = 325 Ma), and exhibit a narrow range in epsilon Nd-325 (-2.6 to -3.5). They are closely affiliated compositionally with the biotite granites in the distant, more voluminous Nejdek massif (Czech Republic). The two-mica granites (Schwarzenberg granite suite, Lauter) are Si-rich (74-77 wt% SiO2) and mildly to strongly peraluminous (A/CNK = 1.17- 1.26). The granites from Schwarzenberg Lire distinctly higher in their Sr(i)ratios (0.709-0.713; t = 325 Ma) and possess lower values of epsilon Nd-325 (-4.9 to -5.2) relative to the biotite granites. The Lauter granites have a Nd-isotopic composition between -3.6 and -4.0 (t = 325 Ma). Mean Th-U-total Pb uraninite ages (Ma +/- 2 sigma) obtained for the granites from the Aue Suite (324.3 +/- 3. 1), Beierfeld (323.7 +/- 3.1), Bernsbach (320.7 +/- 2.9), Schwarzenberg (323.3 +/- 2.4), and the Kirchberg granite al Burkersdorf (322.7 +/- 3.5) indicate that magmatism in the ASGZ commenced in the Namurian and took place early within the major episode of granite formation in the Erzgebirge-Vogtland zone (327-318 Ma). Geochemical and mineralogical patterns of variably altered samples imply that the ASGZ granites are unlikely to have significantly contributed to the formation of spatially associated metalliferous ore deposits (Sn, W, Mo, Ph, Zn, Bi, Co, Ni), except for uranium. In particular the Aue granite suite should have served as major Source for U accumulated in the economically important post-granitic deposits of Schneeberg and Schlema-Alberoda.
Considering the increasing number and complexity of ground-motion prediction equations available for seismic hazard assessment, there is a definite need for an efficient, quantitative, and robust method to select and rank these models for a particular region of interest. In a recent article, Scherbaum et al. (2009) have suggested an information- theoretic approach for this purpose that overcomes several shortcomings of earlier attempts at using data-driven ground- motion prediction equation selection procedures. The results of their theoretical study provides evidence that in addition to observed response spectra, macroseismic intensity data might be useful for model selection and ranking. We present here an applicability study for this approach using response spectra and macroseismic intensities from eight Californian earthquakes. A total of 17 ground-motion prediction equations, from different regions, for response spectra, combined with the equation of Atkinson and Kaka (2007) for macroseismic intensities are tested for their relative performance. The resulting data-driven rankings show that the models that best estimate ground motion in California are, as one would expect, Californian and western U. S. models, while some European models also perform fairly well. Moreover, the model performance appears to be strongly dependent on both distance and frequency. The relative information of intensity versus response spectral data is also explored. The strong correlation we obtain between intensity-based rankings and spectral-based ones demonstrates the great potential of macroseismic intensities data for model selection in the context of seismic hazard assessment.
The hydrological cycle is a dynamic system by its nature, but sometimes accelerated through anthropogenic activity. A "hydrological change" (i.e. a water cycle that is significantly changing over a longer period of time) can be very different in character, depending on the specific natural conditions and the underlying spatial and temporal scales. Such changes may affect the availability and quality of water as essential pre-requisites for human development and ecosystem stability. Hydrological extremes, such as floods and droughts, may also be affected, what is also vitally important, because of their profound economic and societal impacts. Anthropogenically induced hydrological change can be attributed to three main external causes: first, the Earth's climate is changing significantly and thus directly affecting the terrestrial hydro-systems via the exchange of energy and heat. The second major issue is the land cover and its management that has been modified fundamentally by conversion of land for agriculture, forestry, and other purposes such as industrialisation and urbanisation. Finally, water resources are being used more than ever for human development, especially for agriculture, industrial activities, and navigation. If the regional terrestrial hydrological cycle is changing and counter-measures are desirable, it is from a scientific perspective mandatory to understand the extent and nature of such changes, and, especially, to identify their possible anthropogenic origin. There are, however, fundamental gaps in our knowledge, in particular about the role of feedbacks between individual processes and compartments of the hydrological cycle or the relevance of the interactions with other sub-systems of our planet, such as the atmosphere or the vegetation. This paper mentions several examples of hydrological change and discusses their identification, interaction processes, and feedback mechanisms, along with modelling issues. The possibilities and limitations of modelling are demonstrated by means of two studies: one from the river-lake system on the Middle-Havel River and one from the catchment of the Wahnbach Reservoir. The applied model systems comprise a series of consecutively coupled individual models (so-called one-way-coupling). Model systems that are able reflect feedback effects (two-way- coupling) are still in the development stage. It became clear that the applied model systems were able to reproduce the observed dynamics of the hydrological cycle and of selected matter fluxes. However, one has to be aware that the simulated time periods and scenarios represent rather moderately transient conditions, what is the justification why the one-way-coupling seems to be applicable. Furthermore, it was shown that the modelling uncertainty is considerably large. Nevertheless, this uncertainty can be distinguished from effects of changed internal systems dynamics or from changed boundary conditions, what is a basis for the usability of such model systems for prognostic purposes.
The Younger Dryas event, which began approximately 12,900 years ago, was a period of rapid cooling in the Northern Hemisphere, driven by large-scale reorganizations of patterns of atmospheric and oceanic circulation(1-3). Environmental changes during this period have been documented by both proxy-based reconstructions(3) and model simulations(4), but there is currently no consensus on the exact mechanisms of onset, stabilization or termination of the Younger Dryas(5-8). Here we present high-resolution records from two sediment cores obtained from Lake Krakenes in western Norway and the Nordic seas. Multiple proxies from Lake Krakenes are indicative of rapid alternations between glacial growth and melting during the later Younger Dryas. Meanwhile, reconstructed sea surface temperature and salinity from the Nordic seas show an alternation between sea-ice cover and the influx of warm, salty North Atlantic waters. We suggest that the influx of warm water enabled the westerly wind systems to drift northward, closer to their present-day positions. The winds thus brought relatively warm maritime air to Northern Europe, resulting in rising temperatures and the melting of glaciers. Subsequent input of this fresh meltwater into the ocean spurred the formation of sea ice, which forced the westerly winds back to the south, cooling Northern Europe. We conclude that rapid alternations between these two states immediately preceded the termination of the Younger Dryas and the permanent transition to an interglacial state.
There are far-reaching conceptual similarities between bi-static surface georadar and post-stack, "zero-offset" seismic reflection data, which is expressed in largely identical processing flows. One important difference is, however, that standard deconvolution algorithms routinely used to enhance the vertical resolution of seismic data are notoriously problematic or even detrimental to the overall signal quality when applied to surface georadar data. We have explored various options for alleviating this problem and have tested them on a geologically well-constrained surface georadar dataset. Standard stochastic and direct deterministic deconvolution approaches proved to be largely unsatisfactory. While least-squares-type deterministic deconvolution showed some promise, the inherent uncertainties involved in estimating the source wavelet introduced some artificial "ringiness". In contrast, we found spectral balancing approaches to be effective, practical and robust means for enhancing the vertical resolution of surface georadar data, particularly, but not exclusively, in the uppermost part of the georadar section, which is notoriously plagued by the interference of the direct air- and groundwaves. For the data considered in this study, it can be argued that band- limited spectral blueing may provide somewhat better results than standard band-limited spectral whitening, particularly in the uppermost part of the section affected by the interference of the air- and groundwaves. Interestingly, this finding is consistent with the fact that the amplitude spectrum resulting from least-squares-type deterministic deconvolution is characterized by a systematic enhancement of higher frequencies at the expense of lower frequencies and hence is blue rather than white. It is also consistent with increasing evidence that spectral "blueness" is a seemingly universal, albeit enigmatic, property of the distribution of reflection coefficients in the Earth. Our results therefore indicate that spectral balancing techniques in general and spectral blueing in particular represent simple, yet effective means of enhancing the vertical resolution of surface georadar data and, in many cases, could turn out to be a preferable alternative to standard deconvolution approaches.
Trace element concentrations in aqueous fluids in equilibrium with haplogranitic melt were determined in situ at elevated P-T conditions using hydrothermal diamond-anvil cells and synchrotron-radiation XRF microanalyses. Time- resolved analyses showed that the Rb and Sr concentrations in the fluids became constant in less than 2000 s at all temperatures (500 to 780 degrees C). Although fluid-melt equilibration was very rapid, the change in the concentration of both elements in the fluid with temperature was fairly small (a slight increase for Rb and a slight decrease for Sr). This permitted partitioning data for Rb and Sr between haplogranitic melt and H2O or NaCl+KCl+HCl aqueous solutions at 750 degrees C and 200 to 700 MPa to be obtained from EMP analyses of the quenched melt and the in situ SR-XRF analyses of the equilibrated fluid. The resulting D-Rb(f/m) and D-Sr(f/m) were 0.01 +/- 0.002 and 0.006 +/- 0.001 for water as starting fluid, and increased to 0.47 +/- 0.08 and 0.23 +/- 0.03 for 3.56 m (NaCl+KCl)+0.04 in HCl at pressures of 224 to 360 MPa. In the experiments with H2O as starting fluid, the partition coefficients increased with pressure, i.e. D- Rb(f/m) from 0.01 +/- 0.002 to 0.22 +/- 0.02 and D-Sr(f/m) from 0.006 0.001 to 0.02 +/- 0.005 with a change in pressure from 360 to 700 MPa. At pressures to 360 MPa, the Rb/Sr ratio in the fluid was found to be independent of the initial salt concentration (Rb/Sr = 1.45 +/- 0.6). This ratio increased to 7.89 +/- 1.95 at 700 MPa in experiments with chloride free fluids, which indicates different changes in the Rb and Sr speciation with pressure.
The low-productivity South Pacific Gyre (SPG) is Earth's largest oceanic province. Its sediment accumulates extraordinarily slowly (0.1-1 m per million years). This sediment contains a living community that is characterized by very low biomass and very low metabolic activity. At every depth in cored SPG sediment, mean cell abundances are 3 to 4 orders of magnitude lower than at the same depths in all previously explored subseafloor communities. The net rate of respiration by the subseafloor sedimentary community at each SPG site is 1 to 3 orders of magnitude lower than the rates at previously explored sites. Because of the low respiration rates and the thinness of the sediment, interstitial waters are oxic throughout the sediment column in most of this region. Consequently, the sedimentary community of the SPG is predominantly aerobic, unlike previously explored subseafloor communities. Generation of H-2 by radiolysis of water is a significant electron-donor source for this community. The per-cell respiration rates of this community are about 2 orders of magnitude higher (in oxidation/reduction equivalents) than in previously explored anaerobic subseafloor communities. Respiration rates and cell concentrations in subseafloor sediment throughout almost half of the world ocean may approach those in SPG sediment.
We present cross-hole P- and S-wave seismic experiments that have been performed along a similar to 100 m long transect for the detailed characterization of a contaminated sedimentary site (Bitterfeld research test site, Germany). We invert the corresponding first break arrival times for the P- and S-wave velocity structure and compare two different strategies to interpret these models in terms of pertinent lithological and geotechnical parameter variations. The first (common) approach is based on directly translating the tomographic velocity models into the parameters of interest (e.g., elastic moduli). The second (zonal) approach first reduces the tomographic parameter information to a limited number of characteristic velocity combinations via k-means cluster analysis. Then, for each zone (cluster) further parameters including uncertainties can be estimated. In the presented case study, Our results indicate that the zonal approach provides an effective means for the integrated interpretation of different co-located data.
The aim of this paper is to characterize the spatio-temporal distribution of Central-Europe seismicity. Specifically, by using a non-parametric statistical approach, the proportional hazard model, leading to an empirical estimation of the hazard function, we provide some constrains on the time behavior of earthquake generation mechanisms. The results indicate that the most conspicuous characteristics of M-w 4.0+ earthquakes is a temporal clustering lasting a couple of years. This suggests that the probability of occurrence increases immediately after a previous event. After a few years, the process becomes almost time independent. Furthermore, we investigate the cluster properties of the seismicity of Central-Europe, by comparing the obtained result with the one of synthetic catalogs generated by the epidemic type aftershock sequences (ETAS) model, which previously have been successfully applied for short term clustering. Our results indicate that the ETAS is not well suited to describe the seismicity as a whole, while it is able to capture the features of the short- term behaviour. Remarkably, similar results have been previously found for Italy using a higher magnitude threshold.
Recognition of an inferred Miocene marine incursion affecting areas from Colombia through Peru and Bolivia and into Argentina is essential to delineate the South American Seaway. In Bolivia, corresponding strata of inferred marine origin have been assigned to the late Miocene Yecua Formation. We carried out high-resolution delta C-13 and delta O-18 isotopic studies on 135 in situ carbonates from 3 outcrops, combined with detailed sedimentologic, paleontologic, and ichnologic analysis. Four less negative delta C-13 excursion levels were recorded that coincide well with beds containing marine body (barnacle) and trace (Ophiomorpha) fossils. These strata are interbedded with red-green beds containing mudcracks, plant roots, gypsum, and trace fossils of the continental Scoyenia ichnofacies. Our data are significant in that they show for the first time four possible short-lived marine incursions in the Bolivian central Andes during the late Miocene. The result is constrained by a new U-Pb date of 7.17 +/- 0.34 Ma at the top of Yecua strata.
Neogene magmatism and its possible causal relationship with hydrocarbon generation in SW Colombia
(2009)
The Cretaceous oil-bearing source and reservoir sedimentary succession in the Putumayo Basin, SW Colombia, was intruded by gabbroic dykes and sills. The petrological and geochemical character of the magmatic rocks shows calc- alkaline tendency, pointing to a subduction-related magmatic event. K/Ar dating of amphibole indicates a Late Miocene to Pliocene age (6.1 +/- A 0.7 Ma) for the igneous episode in the basin. Therefore, we assume the intrusions to be part of the Andean magmatism of the Northern Volcanic Zone (NVZ). The age of the intrusions has significant tectonic and economic implications because it coincides with two regional events: (1) the late Miocene/Pliocene Andean orogenic uplift of most of the sub-Andean regions in Peru, Ecuador and Colombia and (2) a pulse of hydrocarbon generation and expulsion that has reached the gas window. High La/Yb, K/Nb and La/Nb ratios, and the obtained Sr-Nd-Pb isotopic compositions suggest the involvement of subducted sediments and/or the assimilation of oceanic crust of the subducting slab. We discuss the possibility that magma chamber(s) west of the basin, below the Cordillera, did increase the heat flow in the basin causing generation and expulsion of hydrocarbons and CO2.
External climate forcings-such as long-term changes in solar insolation-generate different climate responses in tropical and high latitude regions(1). Documenting the spatial and temporal variability of past climates is therefore critical for understanding how such forcings are translated into regional climate variability. In contrast to the data- richmiddle and high latitudes, high-quality climate-proxy records from equatorial regions are relatively few(2-4), especially from regions experiencing the bimodal seasonal rainfall distribution associated with twice-annual passage of the Intertropical Convergence Zone. Here we present a continuous and well-resolved climate-proxy record of hydrological variability during the past 25,000 years from equatorial East Africa. Our results, based on complementary evidence from seismic-reflection stratigraphy and organic biomarker molecules in the sediment record of Lake Challa near Mount Kilimanjaro, reveal that monsoon rainfall in this region varied at half-precessional (similar to 11,500-year) intervals in phase with orbitally controlled insolation forcing. The southeasterly and northeasterly monsoons that advect moisture from the western Indian Ocean were strengthened in alternation when the inter-hemispheric insolation gradient was at a maximum; dry conditions prevailed when neither monsoon was intensified and modest local March or September insolation weakened the rain season that followed. On sub-millennial timescales, the temporal pattern of hydrological change on the East African Equator bears clear high-northern-latitude signatures, but on the orbital timescale it mainly responded to low-latitude insolation forcing. Predominance of low-latitude climate processes in this monsoon region can be attributed to the low-latitude position of its continental regions of surface air flow convergence, and its relative isolation from the Atlantic Ocean, where prominent meridional overturning circulation more tightly couples low-latitude climate regimes to high-latitude boundary conditions.
The Kurancali ultramafic-mafic cumulate body, an allochthonous ophiolitic sliver in central Anatolia, is characterized by the presence of abundant hydrous phases (phlogopite, pargasite) besides augitic diopside, plagioclase, and accessory amounts of rutile, sphene, apatite, zircon, and calcite. Based on modes of the essential minerals, the olivine-orthopyroxene-free cumulates are grouped as clinopyroxenite, hydrous clinopyroxenite, phlogopitite, hornblendite, layered gabbro, and diorite. Petrographical, mineralogical and geochemical features of the rocks infer crystallization from a hydrous magma having high-K calcalkaline affinity with slightly alkaline character, and point to metasomatised mantle as the magma source. Our evidence implies that the metasomatising component, which modified the composition of the mantle wedge source rock in an intraoceanic subduction zone, was a H2O, alkali and carbonate-rich aluminosilicate fluid and/or melt, probably derived from a subducted slab. We suggest that the metasomatic agents in the subarc mantle led to the generation of a hydrous magma, which produced the Kurancali cumulates in an island-arc basement in a supra-subduction-zone setting during the closure of the Izmir-Ankara-Erzincan branch of the Alpine Neotethys Ocean.
Owen et al. [Owen. R.B.. Potts, R., Behrensmeyer, A.K., Ditchfield, P. 2008. Diatomaceous sediments and environmental change in the Pleistocene Olorgesailie Formation, southern Kenya Rift Valley. Palaeogeography, Palaeochmatology, Palaeoccology, 269. 17-37], Diatomaceous sediments and environmental change in the Pleistocene Olorgesailie Formation. southern Kenya Rift. Palaeogeography Palaeoclimatology Palaeciecology, 269, 17-37) argued that diatom assemblage variations in the Olorgesailie Formation indicate considerable environmental instability with both wetter and drier periods, contradicting the proposed period of lake stability and wet climatic conditions between ca 11 and 0.9 million years ago as proposed by Trauth et al [Trauth, M.H. Maslin, MA.. Demo. A.. Strecker, M R.. 2005. Late Cenozoic moisture history of East Africa. Science 309. 2051-2053., Trauth, M H. Mashn. M.A., Deino, A., Bergner. A G.N.. Diihnforth, M. Strecker. M.R, 2007 High- and low-latitude forcing of Plio-Pleistocene East African climate and hL.man evolution. journal of Human Evolution 53, 475-486] Contrary to the interpretation of our work by O Nen et al. [Owen. R.B, Potts, R. Behrensmeyer, A.K. Ditchfield. P. 2008], we never said that the proposed periods of large lakes were characterized by stable conditions.
Culture studies of denitrifying bacteria predict that denitrification will generate equivalent gradients in the delta N-15 and delta O-18 of deep ocean nitrate. A depth profile of nitrate isotopes from the Hawaii Ocean Time-series Station ALOHA shows less of an increase in delta O-18 than in delta N-15 as one ascends from abyssal waters into the denitrification-impacted mid-depth waters. A box model of the ocean nitrate N and O isotopes indicates that this is the effect of the low latitude nitrate assimilation/regeneration cycle: organic N sinking out of the surface spreads the high-delta N-15 signal of pelagic denitrification into waters well below and beyond the suboxic zone, whereas the nitrate delta O-18 signal of denitrification can only be transmitted by circulation in the interior.
The dual isotopes of deep nitrate as a constraint on the cycle and budget of oceanic fixed nitrogen
(2009)
We compare the output of an 18-box geochemical model of the ocean with measurements to investigate the controls on both the mean values and variation of nitrate delta N-15 and delta O-18 in the ocean interior. The delta O-18 of nitrate is our focus because it has been explored less in previous work. Denitrification raises the delta N-15 and delta O-18 of mean ocean nitrate by equal amounts above their input values for N-2 fixation (for delta N-15) and nitrification (for delta O-18), generating parallel gradients in the delta N-15 and delta O-18 of deep ocean nitrate. Partial nitrate assimilation in the photic zone also causes equivalent increases in the delta N-15 and delta O-18 of the residual nitrate that can be transported into the interior. However, the regeneration and nitrification of sinking N can be said to decouple the N and O isotopes of deep ocean nitrate, especially when the sinking N is produced in a low latitude region, where nitrate consumption is effectively complete. The delta N-15 of the regenerated nitrate is equivalent to that originally consumed, whereas the regeneration replaces nitrate previously elevated in delta O-18 due to denitrification or nitrate assimilation with nitrate having the delta O-18 of nitrification. This lowers the delta O-18 of mean ocean nitrate and weakens nitrate delta O-18 gradients in the interior relative to those in delta N-15. This decoupling is characterized and quantified in the box model, and agreement with data shows its clear importance in the real ocean. At the same time, the model appears to generate overly strong gradients in both delta O-18 and delta N-15 within the ocean interior and a mean ocean nitrate delta O-18 that is higher than measured. This may be due to, in the model, too strong an impact of partial nitrate assimilation in the Southern Ocean on the delta N-15 and delta O-18 of preformed nitrate and/or too little cycling of intermediate-depth nitrate through the low latitude photic zone.
Quantitative information on vegetation and climate history from the late glacial-Holocene on the Tibetan Plateau is extremely rare. Here, we present palynological results of a 4.30-m-long sediment record collected from Koucha Lake in the Bayan Har Mountains, northeastern Tibetan Plateau. Vegetation change has been traced by biomisation, ordination of pollen data, and calculation of pollen ratios. The application of a pollen-climate calibration set from the eastern Tibetan Plateau to Koucha Lake pollen spectra yielded quantitative climate information. The area was covered by alpine desert/steppe, characteristic of a cold and dry climate (with 50% less precipitation than today) between 16,700 and 14,600 cal yr BP. Steppe vegetation, warm (similar to 1 degrees C higher than today) and wet conditions prevailed between 14,600 and 6600 cal yr BR These findings contradict evidence from other monsoon-influenced areas of Asia, where the early Holocene is thought to have been moist. Low effective moisture on the northeastern Tibetan Plateau was likely due to high temperature and evaporation, even though precipitation levels may have been similar to present- day values. The vegetation changed to tundra around 6600 cal yr BP, indicating that wet and cool climate conditions occurred on the northeastern Tibetan Plateau during the second half of the Holocene.
Measurements of benthic foraminiferal cadmium:calcium (Cd/Ca) have indicated that the glacial-interglacial change in deep North Pacific phosphate (PO4) concentration was minimal which has been taken by some, workers as a sign that the biological pump did not store more carbon in the deep glacial ocean. Here we present sedimentary redox- sensitive trace metal records from Ocean Drilling Program (ODP) Site 882 (NW subarctic Pacific, water depth 3244 m) to make inferences about changes in deep North Pacific oxygenation and thus respired carbon storage - over the past 150,000 yr. These observations are complemented with biogenic barium and opal measurements as indicators for past organic carbon export to separate the influences of deep-water oxygen concentration and sedimentary organic carbon respiration on the redox state of the sediment. Our results suggest that the deep subarctic Pacific water mass was deleted in ox en during glacial maxima, though it was not anoxic. We reconcile our results with the existing benthic foraminiferal Cd/Ca by invoking a decrease in the fraction of the deep ocean nutrient inventory that was preformed, rather than remineralized. This change would have corresponded to an increase in the deep Pacific storage of respired carbon, which Would have lowered atmospheric carbon dioxide (CO2) by sequestering CO2 away from the atmosphere and by increasing ocean alkalinity through a transient dissolution event in the deep sea. The magnitude of change in preformed nutrients suggested by the North Pacific data Would have accounted for a majority of the observed decrease in glacial atmospheric PCO2.