Refine
Year of publication
- 2020 (206) (remove)
Document Type
- Article (142)
- Postprint (33)
- Doctoral Thesis (24)
- Monograph/Edited Volume (3)
- Other (3)
- Review (1)
Is part of the Bibliography
- yes (206)
Keywords
- climate change (7)
- remote sensing (6)
- model (5)
- modelling (5)
- models (5)
- Andes (4)
- Chinese loess (4)
- climate-change (4)
- precipitation (4)
- tectonics (4)
Institute
- Institut für Geowissenschaften (206) (remove)
Seismic receiver arrays have variety of applications in seismology, particularly when the signal enhancement is a prerequisite to detect seismic events, and in situations where installing and maintaining sparse networks are impractical. This thesis has mainly focused on the development of a new approach for seismological source and receiver array design.The proposed approach deals with the array design task as an optimization problem. The criteria and prerequisite constraints in array design task are integrated in objective function definition and evaluation of a optimization process. Three cases are covered in this thesis: (1) a 2-D receiver array geometry optimization, (2) a 3-D source array optimization, and (3) an array application to monitor microseismic data, where the effect of different types of noise are evaluated.
A flexible receiver array design framework implements a customizable scenario modelling and optimization scheme by making use of synthetic seismograms. Using synthetic seismograms to evaluate array performance makes it possible to consider additional constraints, e.g. land ownership, site-specific noise levels or characteristics of the seismic sources under investigation. The use of synthetic array beamforming as an array design criteria is suggested. The framework is customized by designing a 2-D small scale receiver array to monitor earthquake swarm activity in northwest Bohemia/ Vogtland in central Europe. Two sub-functions are defined to verify the accuracy of horizontal slowness estimation; one to suppress aliasing effects due to possible secondary lobes of synthetic array beamforming calculated in horizontal slowness space, and the other to reduce the event's mislocation caused by miscalculation of the horizontal slowness vector. Subsequently, a weighting technique is applied to combine the sub-functions into one single scalar objective function to use in the optimization process.
The idea of optimal array is employed to design a 3-D source array, given a well-located catalog of events. The conditions to make source arrays are formulated in four objective functions and a weighted sum technique is used to combine them in one single scalar function. The criteria are: (1) accurate slowness vector estimation, (2) high waveform coherency, (3) low location error and (4) high energy of coda phases. The method is evaluated by two experiments, (1) a synthetic test using realistic synthetic seismograms, (2) using real seismograms, and for each case optimized SA elements are configured using the data from the Vogtland area.
The location of a possible scatterer in the velocity model, that makes the converted/reflected phases, e.g. sp-phases, is retrieved by a grid search method using the optimized SA. The accuracy of the approach and the obtained results demonstrated that the method is applicable to study the crustal structure and the location of crustal scatterers when the strong converted phases are observed in the data and a well-located catalog is available.
Small aperture arrays are employed in seismology for a variety of applications, ranging from pure event detection to monitor and study of microcosmic activities. The monitoring of microseismicity during temporary human activities is often difficult, as the signal-to-noise ratio is very low and noise is strongly increased during the operation. The combination of small aperture seismic arrays with shallow borehole sensors offers a solution. We tested this monitoring approach at two different sites, (1) accompanying a fracking experiment in sedimentary shale at 4~km depth, and (2) above a gas field under depletion. Arrays recordings are compared with recordings available from shallow borehole sensors and examples of detection and location performance of the array are given. The effect of different types of noise at array and borehole stations are compared and discussed.
Seismological agencies play an important role in seismological research and seismic hazard mitigation by providing source parameters of seismic events (location, magnitude, mechanism), and keeping these data accessible in the long term. The availability of catalogues of seismic source parameters is an important requirement for the evaluation and mitigation of seismic hazards, and the catalogues are particularly valuable to the research community as they provide fundamental long-term data in the geophysical sciences. This work is well motivated by the rising demand for developing more robust and efficient methods for routine source parameter estimation, and ultimately generation of reliable catalogues of seismic source parameters. Specifically, the aim of this work is to develop some methods to determine hypocentre location and temporal evolution of seismic sources based on regional and teleseismic observations more accurately, and investigate the potential of these methods to be integrated in near real-time processing.
To achieve this, a location method that considers several events simultaneously and improves the relative location accuracy among nearby events has been provided. This method tries to reduce the biasing effects of the lateral velocity heterogeneities (or equivalently to compensate for limitations and inaccuracies in the assumed velocity model) by calculating a set of timing corrections for each seismic station. The systematic errors introduced into the locations by the inaccuracies in the assumed velocity structure can be corrected without explicitly solving for a velocity model. Application to sets of multiple earthquakes in complex tectonic environments with strongly heterogeneous structure such as subduction zones and plate boundary region demonstrate that this relocation process significantly improves the hypocentre locations compared to standard locations.
To meet the computational demands of this location process, a new open-source software package has been developed that allows for efficient relocation of large-scale multiple seismic events using arrival time data. Upon that, a flexible location framework is provided which can be tailored to various application cases on local, regional, and global scales. The latest version of the software distribution including source codes, a user guide, an example data set, and a change history, is freely available to the community.
The developed relocation algorithm has been modified slightly and then its performance in a simulated near real-time processing has been evaluated. It has been demonstrated that applying the proposed technique significantly reduces the bias in routine locations and enhance the ability to locate the lower magnitude events using only regional arrival data.
Finally, to return to emphasis on global seismic monitoring, an inversion framework has been developed to determine the seismic source time function through direct waveform fitting of teleseismic recordings. The inversion technique can be systematically applied to moderate- size seismic events and has the potential to be performed in near real-time applications. It is exemplified by application to an abnormal seismic event; the 2017 North Korean nuclear explosion.
The presented work and application case studies in this dissertation represent the first step in an effort to establish a framework for automatic, routine generation of reliable catalogues of seismic event locations and source time functions.
Many institutions struggle to tap into the potential of their large archives of radar reflectivity: these data are often affected by miscalibration, yet the bias is typically unknown and temporally volatile. Still, relative calibration techniques can be used to correct the measurements a posteriori. For that purpose, the usage of spaceborne reflectivity observations from the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) platforms has become increasingly popular: the calibration bias of a ground radar (GR) is estimated from its average reflectivity difference to the spaceborne radar (SR). Recently, Crisologo et al. (2018) introduced a formal procedure to enhance the reliability of such estimates: each match between SR and GR observations is assigned a quality index, and the calibration bias is inferred as a quality-weighted average of the differences between SR and GR. The relevance of quality was exemplified for the Subic S-band radar in the Philippines, which is greatly affected by partial beam blockage.
The present study extends the concept of quality-weighted averaging by accounting for path-integrated attenuation (PIA) in addition to beam blockage. This extension becomes vital for radars that operate at the C or X band. Correspondingly, the study setup includes a C-band radar that substantially overlaps with the S-band radar. Based on the extended quality-weighting approach, we retrieve, for each of the two ground radars, a time series of calibration bias estimates from suitable SR overpasses. As a result of applying these estimates to correct the ground radar observations, the consistency between the ground radars in the region of overlap increased substantially. Furthermore, we investigated if the bias estimates can be interpolated in time, so that ground radar observations can be corrected even in the absence of prompt SR overpasses. We found that a moving average approach was most suitable for that purpose, although limited by the absence of explicit records of radar maintenance operations.
Many institutions struggle to tap into the potential of their large archives of radar reflectivity: these data are often affected by miscalibration, yet the bias is typically unknown and temporally volatile. Still, relative calibration techniques can be used to correct the measurements a posteriori. For that purpose, the usage of spaceborne reflectivity observations from the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) platforms has become increasingly popular: the calibration bias of a ground radar (GR) is estimated from its average reflectivity difference to the spaceborne radar (SR). Recently, Crisologo et al. (2018) introduced a formal procedure to enhance the reliability of such estimates: each match between SR and GR observations is assigned a quality index, and the calibration bias is inferred as a quality-weighted average of the differences between SR and GR. The relevance of quality was exemplified for the Subic S-band radar in the Philippines, which is greatly affected by partial beam blockage.
The present study extends the concept of quality-weighted averaging by accounting for path-integrated attenuation (PIA) in addition to beam blockage. This extension becomes vital for radars that operate at the C or X band. Correspondingly, the study setup includes a C-band radar that substantially overlaps with the S-band radar. Based on the extended quality-weighting approach, we retrieve, for each of the two ground radars, a time series of calibration bias estimates from suitable SR overpasses. As a result of applying these estimates to correct the ground radar observations, the consistency between the ground radars in the region of overlap increased substantially. Furthermore, we investigated if the bias estimates can be interpolated in time, so that ground radar observations can be corrected even in the absence of prompt SR overpasses. We found that a moving average approach was most suitable for that purpose, although limited by the absence of explicit records of radar maintenance operations.
Geomechanical and petrological characterisation of exposed slip zones, Alpine Fault, New Zealand
(2020)
The Alpine Fault is a large, plate-bounding, strike-slip fault extending along the north-western edge of the Southern Alps, South Island, New Zealand. It regularly accommodates large (MW > 8) earthquakes and has a high statistical probability of failure in the near future, i.e., is late in its seismic cycle. This pending earthquake and associated co-seismic landslides are expected to cause severe infrastructural damage that would affect thousands of people, so it presents a substantial geohazard. The interdisciplinary study presented here aims to characterise the fault zone’s 4D (space and time) architecture, because this provides information about its rheological properties that will enable better assessment of the hazard
the fault poses.
The studies undertaken include field investigations of principal slip zone fault gouges exposed
along strike of the fault, and subsequent laboratory analyses of these outcrop and additional borehole samples. These observations have provided new information on (I) characteristic microstructures down to the nanoscale that indicate which deformation mechanisms operated within the rocks, (II) mineralogical information that constrains the fault’s geomechanical behaviour and (III) geochemical compositional information that allows the influence of fluid- related alteration processes on material properties to be unraveled.
Results show that along-strike variations of fault rock properties such as microstructures and mineralogical composition are minor and / or do not substantially influence fault zone architecture. They furthermore provide evidence that the architecture of the fault zone, particularly its fault core, is more complex than previously considered, and also more complex than expected for this sort of mature fault cutting quartzofeldspathic rocks. In particular our results strongly suggest that the fault has more than one principal slip zone, and that these form an anastomosing network extending into the basement below the cover of Quaternary sediments.
The observations detailed in this thesis highlight that two major processes, (I) cataclasis and (II) authigenic mineral formation, are the major controls on the rheology of the Alpine Fault. The velocity-weakening behaviour of its fault gouge is favoured by abundant nanoparticles
promoting powder lubrication and grain rolling rather than frictional sliding. Wall-rock fragmentation is accompanied by co-seismic, fluid-assisted dilatancy that is recorded by calcite cementation. This mineralisation, along with authigenic formation of phyllosilicates, quickly alters the petrophysical fault zone properties after each rupture, restoring fault competency. Dense networks of anastomosing and mutually cross-cutting calcite veins and intensively reworked gouge matrix demonstrate that strain repeatedly localised within the narrow fault gouge. Abundantly undeformed euhedral chlorite crystallites and calcite veins cross-cutting both fault gouge and gravels that overlie basement on the fault’s footwall provide evidence that the processes of authigenic phyllosilicate growth, fluid-assisted dilatancy and associated fault healing are processes active particularly close to the Earth’s surface in this fault zone.
Exposed Alpine Fault rocks are subject to intense weathering as direct consequence of abundant orogenic rainfall associated with the fault’s location at the base of the Southern Alps. Furthermore, fault rock rheology is substantially affected by shallow-depth conditions such as the juxtaposition of competent hanging wall fault rocks on poorly consolidated footwall sediments. This means microstructural, mineralogical and geochemical properties of the exposed fault rocks may differ substantially from those at deeper levels, and thus are not characteristic of the majority of the fault rocks’ history. Examples are (I) frictionally weak smectites found within the fault gouges being artefacts formed at temperature conditions, and imparting petrophysical properties that are not typical for most of fault rocks of the Alpine Fault, (II) grain-scale dissolution resulting from subaerial weathering rather than deformation by pressure-solution processes and (III) fault gouge geometries being more complex than expected for deeper counterparts.
The methodological approaches deployed in analyses of this, and other fault zones, and the major results of this study are finally discussed in order to contextualize slip zone investigations of fault zones and landslides. Like faults, landslides are major geohazards, which highlights the importance of characterising their geomechanical properties. Similarities between faults, especially those exposed to subaerial processes, and landslides, include mineralogical composition and geomechanical behaviour. Together, this ensures failure occurs predominantly by cataclastic processes, although aseismic creep promoted by weak phyllosilicates is not uncommon. Consequently, the multidisciplinary approach commonly used to investigate fault zones may contribute to increase the understanding of landslide faulting processes and the assessment of their hazard potential.
Geomorphology seeks to characterize the forms, rates, and magnitudes of sediment and water transport that sculpt landscapes. This is generally referred to as earth surface processes, which incorporates the influence of biologic (e.g., vegetation), climatic (e.g., rainfall), and tectonic (e.g., mountain uplift) factors in dictating the transport of water and eroded material. In mountains, high relief and steep slopes combine with strong gradients in rainfall and vegetation to create dynamic expressions of earth surface processes. This same rugged topography presents challenges in data collection and process measurement, where traditional techniques involving detailed observations or physical sampling are difficult to apply at the scale of entire catchments. Herein lies the utility of remote sensing. Remote sensing is defined as any measurement that does not disturb the natural environment, typically via acquisition of images in the visible- to radio-wavelength range of the electromagnetic spectrum. Remote sensing is an especially attractive option for measuring earth surface processes, because large areal measurements can be acquired at much lower cost and effort than traditional methods. These measurements cover not only topographic form, but also climatic and environmental metrics, which are all intertwined in the study of earth surface processes. This dissertation uses remote sensing data ranging from handheld camera-based photo surveying to spaceborne satellite observations to measure the expressions, rates, and magnitudes of earth surface processes in high-mountain catchments of the Eastern Central Andes in Northwest Argentina. This work probes the limits and caveats of remote sensing data and techniques applied to geomorphic research questions, and presents important progress at this disciplinary intersection.