Refine
Year of publication
- 2010 (90) (remove)
Document Type
- Article (66)
- Doctoral Thesis (21)
- Other (2)
- Monograph/Edited Volume (1)
Keywords
- Arrayseismologie (2)
- Cosmogenic nuclides (2)
- Erdbeben (2)
- Erdbebenkatalog (2)
- Erdbebenschwarm 2008/09 (2)
- InSAR (2)
- Kosmogene Nuklide (2)
- Standorteffekte (2)
- Vogtland/West Bohemia (2)
- Vogtland/Westböhmen (2)
Institute
- Institut für Geowissenschaften (90) (remove)
Based on technological advances made within the past decades, ground-penetrating radar (GPR) has become a well-established, non-destructive subsurface imaging technique. Catalyzed by recent demands for high-resolution, near-surface imaging (e.g., the detection of unexploded ordnances and subsurface utilities, or hydrological investigations), the quality of today's GPR-based, near-surface images has significantly matured. At the same time, the analysis of oil and gas related reflection seismic data sets has experienced significant advances. Considering the sensitivity of attribute analysis with respect to data positioning in general, and multi-trace attributes in particular, trace positioning accuracy is of major importance for the success of attribute-based analysis flows. Therefore, to study the feasibility of GPR-based attribute analyses, I first developed and evaluated a real-time GPR surveying setup based on a modern tracking total station (TTS). The combination of current GPR systems capability of fusing global positioning system (GPS) and geophysical data in real-time, the ability of modern TTS systems to generate a GPS-like positional output and wireless data transmission using radio modems results in a flexible and robust surveying setup. To elaborate the feasibility of this setup, I studied the major limitations of such an approach: system cross-talk and data delays known as latencies. Experimental studies have shown that when a minimal distance of ~5 m between the GPR and the TTS system is considered, the signal-to-noise ratio of the acquired GPR data using radio communication equals the one without radio communication. To address the limitations imposed by system latencies, inherent to all real-time data fusion approaches, I developed a novel correction (calibration) strategy to assess the gross system latency and to correct for it. This resulted in the centimeter trace accuracy required by high-frequency and/or three-dimensional (3D) GPR surveys. Having introduced this flexible high-precision surveying setup, I successfully demonstrated the application of attribute-based processing to GPR specific problems, which may differ significantly from the geological ones typically addressed by the oil and gas industry using seismic data. In this thesis, I concentrated on archaeological and subsurface utility problems, as they represent typical near-surface geophysical targets. Enhancing 3D archaeological GPR data sets using a dip-steered filtering approach, followed by calculation of coherency and similarity, allowed me to conduct subsurface interpretations far beyond those obtained by classical time-slice analyses. I could show that the incorporation of additional data sets (magnetic and topographic) and attributes derived from these data sets can further improve the interpretation. In a case study, such an approach revealed the complementary nature of the individual data sets and, for example, allowed conclusions about the source location of magnetic anomalies by concurrently analyzing GPR time/depth slices to be made. In addition to archaeological targets, subsurface utility detection and characterization is a steadily growing field of application for GPR. I developed a novel attribute called depolarization. Incorporation of geometrical and physical feature characteristics into the depolarization attribute allowed me to display the observed polarization phenomena efficiently. Geometrical enhancement makes use of an improved symmetry extraction algorithm based on Laplacian high-boosting, followed by a phase-based symmetry calculation using a two-dimensional (2D) log-Gabor filterbank decomposition of the data volume. To extract the physical information from the dual-component data set, I employed a sliding-window principle component analysis. The combination of the geometrically derived feature angle and the physically derived polarization angle allowed me to enhance the polarization characteristics of subsurface features. Ground-truth information obtained by excavations confirmed this interpretation. In the future, inclusion of cross-polarized antennae configurations into the processing scheme may further improve the quality of the depolarization attribute. In addition to polarization phenomena, the time-dependent frequency evolution of GPR signals might hold further information on the subsurface architecture and/or material properties. High-resolution, sparsity promoting decomposition approaches have recently had a significant impact on the image and signal processing community. In this thesis, I introduced a modified tree-based matching pursuit approach. Based on different synthetic examples, I showed that the modified tree-based pursuit approach clearly outperforms other commonly used time-frequency decomposition approaches with respect to both time and frequency resolutions. Apart from the investigation of tuning effects in GPR data, I also demonstrated the potential of high-resolution sparse decompositions for advanced data processing. Frequency modulation of individual atoms themselves allows to efficiently correct frequency attenuation effects and improve resolution based on shifting the average frequency level. GPR-based attribute analysis is still in its infancy. Considering the growing widespread realization of 3D GPR studies there will certainly be an increasing demand towards improved subsurface interpretations in the future. Similar to the assessment of quantitative reservoir properties through the combination of 3D seismic attribute volumes with sparse well-log information, parameter estimation in a combined manner represents another step in emphasizing the potential of attribute-driven GPR data analyses.
In the humid tropics, continuing high deforestation rates are seen alongside an increasing expansion of secondary forests. In order to understand and model the consequences of these dynamic land-use changes for regional water cycles, the response of soil hydraulic properties to forest disturbance and recovery has to be quantified.At a site in the Brazilian Amazonia, we annually monitored soil infiltrability and saturated hydraulic conductivity (K-s) at 12.5, 20 cm, and 50 cm soil depth after manual forest conversion to pasture (year zero to four after pasture establishment), and during secondary succession after pasture abandonment (year zero to seven after pasture abandonment). We evaluated the hydrological consequences of the detected changes by comparing the soil hydraulic properties with site-specific rainfall intensities and hydrometric observations. Within one year after grazing started, infiltrability and K-s at 12.5 and 20 cm depth decreased by up to one order of magnitude to levels which are typical for 20-year-old pasture. In the three subsequent monitoring years, infiltrability and K-s remained stable. Land use did not impact on subsoil permeability. Whereas infiltrability values are large enough to allow all rainwater to infiltrate even after the conversion, the sudden decline of near-surface K-s is of hydrological relevance as perched water tables and overland flow occur more often on pastures than in forests at our study site. After pasture abandonment and during secondary succession, seven years of recovery did not suffice to significantly increase infiltrability and K-s at 12.5 depth although a slight recovery is obvious. At 20 cm soil depth, we detected a positive linear increase within the seven-year time frame but annual means did not differ significantly. Although more than a doubling of infiltrability and K-s is still required to achieve pre-disturbance levels, which will presumably take more than a decade, the observed slight increases of K-s might already decrease the probability of perched water table generation and overland flow development well before complete recovery.
In this investigation we analyze the evolution of forms that constitute the Copiapo River estuarine system during the past 30 years. Through photo interpretation and field work is possible to realize that during the period, the essential forms of the estuary haven't manifest significant changes, on the other hand, estuarine complex such as the estuarine lagoon, sandy beach and dunes have presented modifications. The most significant morphological elements in order to understand these phenomena of changes have been the dynamics of waves and their spatial correlation with the width of the surf zone and dune continuum. The basic shapes present stabilization in their dynamic, which allows us to consider the estuarine system of Copiapo, a paleoestuario.
Soil respiration is the second largest flux in the global carbon cycle, yet the underlying below-ground process, carbon dioxide (CO2) production, is not well understood because it can not be measured in the field. CO2 production has frequently been calculated from the vertical CO2 diffusive flux divergence, known as "soil-CO2 profile method". This relatively simple model requires knowledge of soil CO2 concentration profiles and soil diffusive properties. Application of the method for a tropical lowland forest soil in Panama gave inconsistent results when using diffusion coefficients (D) calculated based on relationships with soil porosity and moisture ("physically modeled" D). Our objective was to investigate whether these inconsistencies were related to (1) the applied interpolation and solution methods and/or (2) uncertainties in the physically modeled profile of D. First, we show that the calculated CO2 production strongly depends on the function used to interpolate between measured CO2 concentrations. Secondly, using an inverse analysis of the soil-CO2 profile method, we deduce which D would be required to explain the observed CO2 concentrations, assuming the model perception is valid. In the top soil, this inversely modeled D closely resembled the physically modeled D. In the deep soil, however, the inversely modeled D increased sharply while the physically modeled D did not. When imposing a constraint during the fit parameter optimization, a solution could be found where this deviation between the physically and inversely modeled D disappeared. A radon (Rn) mass balance model, in which diffusion was calculated based on the physically modeled or constrained inversely modeled D, simulated observed Rn profiles reasonably well. However, the CO2 concentrations which corresponded to the constrained inversely modeled D were too small compared to the measurements. We suggest that, in well-structured soils, a missing description of steady state CO2 exchange fluxes across water-filled pores causes the soil-CO2 profile method to fail. These fluxes are driven by the different diffusivities in inter- vs. intra-aggregate pores which create permanent CO2 gradients if separated by a "diffusive water barrier". These results corroborate other studies which have shown that the theory to treat gas diffusion as homogeneous process, a precondition for use of the soil-CO2 profile method, is inaccurate for pore networks which exhibit spatial separation between CO2 production and diffusion out of the soil.
Investigation of transient soil moisture profiles yields valuable information of near- surface processes. A recently developed reconstruction algorithm based on the telegraph equation allows the inverse estimation of soil moisture profiles along coated, three rod TDR probes. Laboratory experiments were carried out to prove the results of the inversion and to understand the influence of probe rod deformation and solid objects close to the probe in heterogonous media. Differences in rod geometry can lead to serious misinterpretations in the soil moisture profile but have small influence on the average soil moisture along the probe. Solids in the integration volume have almost no effect on average soil moisture but result in locally slightly decreased moisture values. Inverted profiles obtained in a loamy soil with a clay content of about 16% were in good agreement with independent measurements.
Investigation of transient soil moisture profiles yields valuable information of near- surface processes. A recently developed reconstruction algorithm based on the telegraph equation allows the inverse estimation of soil moisture profiles along coated, three rod TDR probes. Laboratory experiments were carried out to prove the results of the inversion and to understand the influence of probe rod deformation and solid objects close to the probe in heterogeneous media. Differences in rod geometry can lead to serious misinterpretations in the soil moisture profile, but have small influence on the average soil moisture along the probe. Solids in the integration volume have almost no effect on average soil moisture, but result in locally slightly decreased moisture values. Inverted profiles obtained in a loamy soil with a clay content of about 16% were in good agreement with independent measurements.
Aim: Fossil pollen spectra from lake sediments on the Tibetan Plateau have been used for qualitative climate reconstruction, but no modern pollen-climate calibration set based on lake sediments is available to infer past climate quantitatively. This study aims to develop such a dataset and apply it to fossil data. Location: The Tibetan Plateau, between 30 and 40 degrees N and 87 and 103 degrees E. Methods: We collected surface sediments from 112 lakes and analysed them palynologically. The lakes span a wide range of mean annual precipitation (P-ann; 31-1022 mm), mean annual temperature (T-ann; -6.5 to 1 degrees C), and mean July temperature (T-July; 2.6-19.7 degrees C). Redundancy analysis showed that the modern pollen spectra are characteristic of their respective vegetation types and local climate. Transfer functions for P-ann, T-ann and T-July were developed with weighted averaging partial least squares. Model performance was assessed by leave-one-out cross-validation. Results: The root mean square errors of prediction (RMSEP) were 104 mm (P-ann), 1.18 degrees C (T-ann) and 1.17 degrees C (T-July). The RMSEPs, when expressed as percentages of the gradient sampled, were 10.6% (P-ann), 15.7% (T-ann) and 11.9% (T-July). These low values indicate the good performance of our models. An application of the models to fossil pollen spectra covering the last c. 50 kyr yielded realistic results for Luanhaizi Lake in the Qilian Mountains on the north-eastern Tibetan Plateau (modern P-ann 480 mm; T-ann-1 degrees C). T-ann and P-ann values similar to present ones were reconstructed for late Marine Isotope Stage 3, with minimum values for the Last Glacial Maximum (c. 300 mm and 2 degrees C below present), and maximum values for the early Holocene (c. 70 mm and 0.5 degrees C greater than present). Main conclusions: The modern pollen-climate calibration set will potentially be useful for quantitative climate reconstructions from lake-sediment pollen spectra from the Tibetan Plateau, an area of considerable climatic and biogeographical importance.
A confocal set-up is presented that improves micro-XRF and XAFS experiment with high-pressure e diamond-anvil cells (DACs) In this experiment a probing volume is defined by the focus of the incoming synchrotron radiation beam and that of a polycapillary X-ray half-lens with a very long working distance, which is placed in front of the fluorescence detector This set-up enhances the quality of the fluorescence and XAFS spectra, and thus the sensitivity for detecting elements at low concentrations. It efficiently suppresses signal from outside the sample chamber, which stems from elastic and inelastic scattering of the incoming beam by the diamond anvils as well as from excitation of fluorescence from the body of the DAC
Situated in an active tectonic region, Santiago de Chile, the country´s capital with more than six million inhabitants, faces tremendous earthquake hazard. Macroseismic data for the 1985 Valparaiso and the 2010 Maule events show large variations in the distribution of damage to buildings within short distances indicating strong influence of local sediments and the shape of the sediment-bedrock interface on ground motion. Therefore, a temporary seismic network was installed in the urban area for recording earthquake activity, and a study was carried out aiming to estimate site amplification derived from earthquake data and ambient noise. The analysis of earthquake data shows significant dependence on the local geological structure with regards to amplitude and duration. Moreover, the analysis of noise spectral ratios shows that they can provide a lower bound in amplitude for site amplification and, since no variability in terms of time and amplitude is observed, that it is possible to map the fundamental resonance frequency of the soil for a 26 km x 12 km area in the northern part of the Santiago de Chile basin. By inverting the noise spectral rations, local shear wave velocity profiles could be derived under the constraint of the thickness of the sedimentary cover which had previously been determined by gravimetric measurements. The resulting 3D model was derived by interpolation between the single shear wave velocity profiles and shows locally good agreement with the few existing velocity profile data, but allows the entire area, as well as deeper parts of the basin, to be represented in greater detail. The wealth of available data allowed further to check if any correlation between the shear wave velocity in the uppermost 30 m (vs30) and the slope of topography, a new technique recently proposed by Wald and Allen (2007), exists on a local scale. While one lithology might provide a greater scatter in the velocity values for the investigated area, almost no correlation between topographic gradient and calculated vs30 exists, whereas a better link is found between vs30 and the local geology. When comparing the vs30 distribution with the MSK intensities for the 1985 Valparaiso event it becomes clear that high intensities are found where the expected vs30 values are low and over a thick sedimentary cover. Although this evidence cannot be generalized for all possible earthquakes, it indicates the influence of site effects modifying the ground motion when earthquakes occur well outside of the Santiago basin. Using the attained knowledge on the basin characteristics, simulations of strong ground motion within the Santiago Metropolitan area were carried out by means of the spectral element technique. The simulation of a regional event, which has also been recorded by a dense network installed in the city of Santiago for recording aftershock activity following the 27 February 2010 Maule earthquake, shows that the model is capable to realistically calculate ground motion in terms of amplitude, duration, and frequency and, moreover, that the surface topography and the shape of the sediment bedrock interface strongly modify ground motion in the Santiago basin. An examination on the dependency of ground motion on the hypocenter location for a hypothetical event occurring along the active San Ramón fault, which is crossing the eastern outskirts of the city, shows that the unfavorable interaction between fault rupture, radiation mechanism, and complex geological conditions in the near-field may give rise to large values of peak ground velocity and therefore considerably increase the level of seismic risk for Santiago de Chile.
The closed-chamber method is the most common approach to determine CH4 fluxes in peatlands. The concentration change in the chamber is monitored over time, and the flux is usually calculated by the slope of a linear regression function. Theoretically, the gas exchange cannot be constant over time but has to decrease, when the concentration gradient between chamber headspace and soil air decreases. In this study, we test whether we can detect this non- linearity in the concentration change during the chamber closure with six air samples. We expect generally a low concentration gradient on dry sites (hummocks) and thus the occurrence of exponential concentration changes in the chamber due to a quick equilibrium of gas concentrations between peat and chamber headspace. On wet (flarks) and sedge- covered sites (lawns), we expect a high gradient and near-linear concentration changes in the chamber. To evaluate these model assumptions, we calculate both linear and exponential regressions for a test data set (n = 597) from a Finnish mire. We use the Akaike Information Criterion with small sample second order bias correction to select the best-fitted model. 13.6%, 19.2% and 9.8% of measurements on hummocks, lawns and flarks, respectively, were best fitted with an exponential regression model. A flux estimation derived from the slope of the exponential function at the beginning of the chamber closure can be significantly higher than using the slope of the linear regression function. Non-linear concentration-overtime curves occurred mostly during periods of changing water table. This could be due to either natural processes or chamber artefacts, e.g. initial pressure fluctuations during chamber deployment. To be able to exclude either natural processes or artefacts as cause of non-linearity, further information, e.g. CH4 concentration profile measurements in the peat, would be needed. If this is not available, the range of uncertainty can be substantial. We suggest to use the range between the slopes of the exponential regression at the beginning and at the end of the closure time as an estimate of the overall uncertainty.