Refine
Year of publication
- 2020 (206) (remove)
Document Type
- Article (143)
- Postprint (33)
- Doctoral Thesis (24)
- Monograph/Edited Volume (3)
- Other (2)
- 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)
Understanding the relation between injection-induced seismic moment release and operational parameters is crucial for early identification of possible seismic hazards associated with fluid-injection projects. We conducted laboratory fluid-injection experiments on permeable sandstone samples containing a critically stressed fault at different fluid pressurization rates. The observed fluid-induced fault deformation is dominantly aseismic. Fluid-induced stick-slip and fault creep reveal that total seismic moment release of acoustic emission (AE) events is related to total injected volume, independent of respective fault slip behavior. Seismic moment release rate of AE scales with measured fault slip velocity. For injection-induced fault slip in a homogeneous pressurized region, released moment shows a linear scaling with injected volume for stable slip (steady slip and fault creep), while we find a cubic relation for dynamic slip. Our results highlight that monitoring evolution of seismic moment release with injected volume in some cases may assist in discriminating between stable slip and unstable runaway ruptures.
To find out the future of nowadays reef ecosystem turnover under the environmental stresses such as global warming and ocean acidification, analogue studies from the geologic past are needed. As a critical time of reef ecosystem innovation, the Permian-Triassic transition witnessed the most severe demise of Phanerozoic reef builders, and the establishment of modern style symbiotic relationships within the reef-building organisms. Being the initial stage of this transition, the Middle Permian (Capitanian) mass extinction coursed a reef eclipse in the early Late Permian, which lead to a gap of understanding in the post-extinction Wuchiapingian reef ecosystem, shortly before the radiation of Changhsingian reefs. Here, this thesis presents detailed biostratigraphic, sedimentological, and palaeoecological studies of the Wuchiapingian reef recovery following the Middle Permian (Capitanian) mass extinction, on the only recorded Wuchiapingian reef setting, outcropping in South China at the Tieqiao section.
Conodont biostratigraphic zonations were revised from the Early Permian Artinskian to the Late Permian Wuchiapingian in the Tieqiao section. Twenty main and seven subordinate conodont zones are determined at Tieqiao section including two conodont zone below and above the Tieqiao reef complex. The age of Tieqiao reef was constrained as early to middle Wuchiapingian.
After constraining the reef age, detailed two-dimensional outcrop mapping combined with lithofacies study were carried out on the Wuchiapingian Tieqiao Section to investigate the reef growth pattern stratigraphically as well as the lateral changes of reef geometry on the outcrop scale. Semi-quantitative studies of the reef-building organisms were used to find out their evolution pattern within the reef recovery. Six reef growth cycles were determined within six transgressive-regressive cycles in the Tieqiao section. The reefs developed within the upper part of each regressive phase and were dominated by different biotas. The timing of initial reef recovery after the Middle Permian (Capitanian) mass extinction was updated to the Clarkina leveni conodont zone, which is earlier than previous understanding. Metazoans such as sponges were not the major components of the Wuchiapingian reefs until the 5th and 6th cycles. So, the recovery of metazoan reef ecosystem after the Middle Permian (Capitanian) mass extinction was obviously delayed. In addition, although the importance of metazoan reef builders such as sponges did increase following the recovery process, encrusting organisms such as Archaeolithoporella and Tubiphytes, combined with microbial carbonate precipitation, still played significant roles to the reef building process and reef recovery after the mass extinction.
Based on the results from outcrop mapping and sedimentological studies, quantitative composition analysis of the Tieqiao reef complex were applied on selected thin sections to further investigate the functioning of reef building components and the reef evolution after the Middle Permian (Capitanian) mass extinction. Data sets of skeletal grains and whole rock components were analyzed. The results show eleven biocommunity clusters/eight rock composition clusters dominated by different skeletal grains/rock components. Sponges, Archaeolithoporella and Tubiphytes were the most ecologically important components within the Wuchiapingian Tieqiao reef, while the clotted micrites and syndepositional cements are the additional important rock components for reef cores. The sponges were important within the whole reef recovery. Tubiphytes were broadly distributed in different environments and played a key-role in the initial reef communities. Archaeolithoporella concentrated in the shallower part of reef cycles (i.e., the upper part of reef core) and was functionally significant for the enlargement of reef volume.
In general, the reef recovery after the Middle Permian (Capitanian) mass extinction has some similarities with the reef recovery following the end-Permian mass extinction. It shows a delayed recovery of metazoan reefs and a stepwise recovery pattern that was controlled by both ecological and environmental factors. The importance of encrusting organisms and microbial carbonates are also similar to most of the other post-extinction reef ecosystems. These findings can be instructive to extend our understanding of the reef ecosystem evolution under environmental perturbation or stresses.
A ground motion logic tree for seismic hazard analysis in the stable cratonic region of Europe
(2020)
Regions of low seismicity present a particular challenge for probabilistic seismic hazard analysis when identifying suitable ground motion models (GMMs) and quantifying their epistemic uncertainty. The 2020 European Seismic Hazard Model adopts a scaled backbone approach to characterise this uncertainty for shallow seismicity in Europe, incorporating region-to-region source and attenuation variability based on European strong motion data. This approach, however, may not be suited to stable cratonic region of northeastern Europe (encompassing Finland, Sweden and the Baltic countries), where exploration of various global geophysical datasets reveals that its crustal properties are distinctly different from the rest of Europe, and are instead more closely represented by those of the Central and Eastern United States. Building upon the suite of models developed by the recent NGA East project, we construct a new scaled backbone ground motion model and calibrate its corresponding epistemic uncertainties. The resulting logic tree is shown to provide comparable hazard outcomes to the epistemic uncertainty modelling strategy adopted for the Eastern United States, despite the different approaches taken. Comparison with previous GMM selections for northeastern Europe, however, highlights key differences in short period accelerations resulting from new assumptions regarding the characteristics of the reference rock and its influence on site amplification.
Probabilistic assessment of seismic hazard and risk over a geographical region presents the modeler with challenges in the characterization of the site amplification that are not present in site-specific assessment. Using site-to-site residuals from a ground motion model fit to observations from the Japanese KiK-net database, correlations between measured local amplifications and mappable proxies such as topographic slope and geology are explored. These are used subsequently to develop empirical models describing amplification as a direct function of slope, conditional upon geological period. These correlations also demonstrate the limitations of inferring 30-m shearwave velocity from slope and applying them directly into ground motion models. Instead, they illustrate the feasibility of deriving spectral acceleration amplification factors directly from sets of observed records, which are calibrated to parameters that can be mapped uniformly on a regional scale. The result is a geologically calibrated amplification model that can be incorporated into national and regional seismic hazard and risk assessment, ensuring that the corresponding total aleatory variability reflects the predictive capability of the mapped site proxy.
High-resolution surface velocities and strain for Anatolia from Sentinel-1 InSAR and GNSS data
(2020)
Measurements of present-day surface deformation are essential for the assessment of long-term seismic hazard. The European Space Agency's Sentinel-1 satellites enable global, high-resolution observation of crustal motion from Interferometric Synthetic Aperture Radar (InSAR). We have developed automated InSAR processing systems that exploit the first similar to 5 years of Sentinel-1 data to measure surface motions for the similar to 800,000-km(2) Anatolian region. Our new 3-D velocity and strain rate fields illuminate deformation patterns dominated by westward motion of Anatolia relative to Eurasia, localized strain accumulation along the North and East Anatolian Faults, and rapid vertical signals associated with anthropogenic activities and to a lesser extent extension across the grabens of western Anatolia. We show that automatically processed Sentinel-1 InSAR data can characterize details of the velocity and strain rate fields with high resolution and accuracy over large regions. These results are important for assessing the relationship between strain accumulation and release in earthquakes. <br /> Plain Language Summary Satellite-based measurements of small rates of motion of the Earth's surface made at high spatial resolutions and over large areas are important for many geophysical applications including improving earthquake hazard models. We take advantage of recent advances in geodetic techniques in order to measure surface velocities and tectonic strain accumulation across the Anatolia region, including the highly seismogenic and often deadly North Anatolian Fault. We show that by combining Sentinel-1 Interferometric Synthetic Aperture Radar (InSAR) data with Global Navigation Satellite System (GNSS) measurements we can enhance our view of surface deformation associated with active tectonics, the earthquake cycle, and anthropogenic processes.
Geochemical processes change the microstructure of rocks and thereby affect their physical behaviour at the macro scale. A micro-computer tomography (micro-CT) scan of a typical reservoir sandstone is used to numerically examine the impact of three spatial alteration patterns on pore morphology, permeability and elastic moduli by correlating precipitation with the local flow velocity magnitude. The results demonstrate that the location of mineral growth strongly affects the permeability decrease with variations by up to four orders in magnitude. Precipitation in regions of high flow velocities is characterised by a predominant clogging of pore throats and a drastic permeability reduction, which can be roughly described by the power law relation with an exponent of 20. A continuous alteration of the pore structure by uniform mineral growth reduces the permeability comparable to the power law with an exponent of four or the Kozeny-Carman relation. Preferential precipitation in regions of low flow velocities predominantly affects smaller throats and pores with a minor impact on the flow regime, where the permeability decrease is considerably below that calculated by the power law with an exponent of two. Despite their complete distinctive impact on hydraulics, the spatial precipitation patterns only slightly affect the increase in elastic rock properties with differences by up to 6.3% between the investigated scenarios. Hence, an adequate characterisation of the spatial precipitation pattern is crucial to quantify changes in hydraulic rock properties, whereas the present study shows that its impact on elastic rock parameters is limited. The calculated relations between porosity and permeability, as well as elastic moduli can be applied for upscaling micro-scale findings to reservoir-scale models to improve their predictive capabilities, what is of paramount importance for a sustainable utilisation of the geological subsurface.
Relative pollen productivity (RPP) estimates are fractionate values, often in relation to Poaceae, that allow vegetation cover to be estimated from pollen counts with the help of models. RPP estimates are especially used in the scientific community in Europe and China, with a few studies in North America. Here we present a comprehensive compilation of available northern hemispheric RPP studies and their results arising from 51 publications with 60 sites and 131 taxa. This compilation allows scientists to identify data gaps in need of further RPP analyses but can also aid them in finding an RPP set for their study region. We also present a taxonomically harmonised, unified RPP dataset for the Northern Hemisphere and subsets for North America (including Greenland), Europe (including arctic Russia), and China, which we generated from the available studies. The unified dataset gives the mean RPP for 55 harmonised taxa as well as fall speeds, which are necessary to reconstruct vegetation cover from pollen counts and RPP values. Data are openly available at https://doi.org/10.1594/PANGAEA.922661 (Wieczorek and Herzschuh, 2020).
Different upper tail indicators exist to characterize heavy tail phenomena, but no comparative study has been carried out so far. We evaluate the shape parameter (GEV), obesity index, Gini index and upper tail ratio (UTR) against a novel benchmark of tail heaviness - the surprise factor. Sensitivity analyses to sample size and changes in scale-to-location ratio are carried out in bootstrap experiments. The UTR replicates the surprise factor best but is most uncertain and only comparable between records of similar length. For samples with symmetric Lorenz curves, shape parameter, obesity and Gini indices provide consistent indications. For asymmetric Lorenz curves, however, the first two tend to overestimate, whereas Gini index tends to underestimate tail heaviness. We suggest the use of a combination of shape parameter, obesity and Gini index to characterize tail heaviness. These indicators should be supported with calculation of the Lorenz asymmetry coefficients and interpreted with caution.
In recent years, nature-based solutions are receiving increasing attention in the field of disaster risk reduction and climate change adaptation as inclusive, no regret approaches. Ecosystem-based adaptation (EbA) can mitigate the impacts of climate change, build resilience and tackle environmental degradation thereby supporting the targets set by the 2030 Agenda, the Paris Agreement and the Sendai Framework. Despite these benefits, EbA is still rarely implemented in practice. To better understand the barriers to implementation, this research examines policy-makers' perceptions of EbA, using an extended version of Protection Motivation Theory as an analytical framework. Through semi-structured interviews with policy-makers at regional and provincial level in Central Vietnam, it was found that EbA is generally considered a promising response option, mainly due to its multiple ecosystem-service benefits. The demand for EbA measures was largely driven by the perceived consequences of natural hazards and climate change. Insufficient perceived response efficacy and time-lags in effectiveness for disaster risk reduction were identified as key impediments for implementation. Pilot projects and capacity building on EbA are important means to overcome these perceptual barriers. This paper contributes to bridging the knowledge-gap on political decision-making regarding EbA and can, thereby, promote its mainstreaming into policy plans.
Marked along-strike changes in stratigraphy, mountain belt morphology, basement exhumation, and deformation styles characterize the Andean retroarc; these changes have previously been related to spatiotemporal variations in the subduction angle. We modeled new apatite fission track and apatite (U-Th-Sm)/He data from nine ranges located between 26 degrees S and 28 degrees S. Using new and previously published data, we constructed a Cretaceous to Pliocene paleogeographic model that delineates a four-stage tectonic evolution: extensional tectonics during the Cretaceous (120-75 Ma), the formation of a broken foreland basin between 55 and 30 Ma, reheating due to burial beneath sedimentary rocks (18-13 Ma), and deformation, exhumation, and surface uplift during the Late Miocene and the Pliocene (13-3 Ma). Our model highlights how preexisting upper plate structures control the deformation patterns of broken foreland basins. Because retroarc deformation predates flat-slab subduction, we propose that slab anchoring may have been the precursor of Eocene-Oligocene compression in the Andean retroarc. Our model challenges models which consider broken foreland basins and retroarc deformation in the NW Argentinian Andes to be directly related to Miocene flat subduction.
Marked along-strike changes in stratigraphy, mountain belt morphology, basement exhumation, and deformation styles characterize the Andean retroarc; these changes have previously been related to spatiotemporal variations in the subduction angle. We modeled new apatite fission track and apatite (U-Th-Sm)/He data from nine ranges located between 26 degrees S and 28 degrees S. Using new and previously published data, we constructed a Cretaceous to Pliocene paleogeographic model that delineates a four-stage tectonic evolution: extensional tectonics during the Cretaceous (120-75 Ma), the formation of a broken foreland basin between 55 and 30 Ma, reheating due to burial beneath sedimentary rocks (18-13 Ma), and deformation, exhumation, and surface uplift during the Late Miocene and the Pliocene (13-3 Ma). Our model highlights how preexisting upper plate structures control the deformation patterns of broken foreland basins. Because retroarc deformation predates flat-slab subduction, we propose that slab anchoring may have been the precursor of Eocene-Oligocene compression in the Andean retroarc. Our model challenges models which consider broken foreland basins and retroarc deformation in the NW Argentinian Andes to be directly related to Miocene flat subduction.
Seismological and seismotectonic analysis of the northwestern Argentine Central Andean foreland
(2020)
After a severe M W 5.7 earthquake on October 17, 2015 in El Galpón in the province of Salta NW Argentina, I installed a local seismological network around the estimated epicenter. The network covered an area characterized by inherited Cretaceous normal faults and neotectonic faults with unknown recurrence intervals, some of which may have been reactivated normal faults. The 13 three-component seismic stations recorded data continuously for 15 months.
The 2015 earthquake took place in the Santa Bárbara System of the Andean foreland, at about 17km depth. This region is the easternmost morphostructural region of the central Andes. As a part of the broken foreland, it is bounded to the north by the Subandes fold-and-thrust belt and the Sierras Pampeanas to the south; to the east lies the Chaco-Paraná basin.
A multi-stage morphotectonic evolution with thick-skinned basement uplift and coeval thin-skinned deformation in the intermontane basins is suggested for the study area. The release of stresses associated with the foreland deformation can result in strong earthquakes, as the study area is known for recurrent and historical, destructive earthquakes. The available continuous record reaches back in time, when the strongest event in 1692 (magnitude 7 or intensity IX) destroyed the city of Esteco. Destructive earthquakes and surface deformation are thus a hallmark of this part of the Andean foreland.
With state-of-the-art Python packages (e.g. pyrocko, ObsPy), a semi-automatic approach is followed to analyze the collected continuous data of the seismological network. The resulting 1435 hypocenter locations consist of three different groups: 1.) local crustal earthquakes (nearly half of the events belong to this group), 2.) interplate activity, of regional distance in the slab of the Nazca-plate, and 3.) very deep earthquakes at about 600km depth. My major interest focused on the first event class. Those crustal events are partly aftershock events of the El Galpón earthquake and a second earthquake, in the south of the same fault. Further events can be considered as background seismicity of other faults within the study area. Strikingly, the seismogenic zone encompass the whole crust and propagates brittle deformation down, close to the Moho.
From the collected seismological data, a local seismic velocity model is estimated, using VELEST. After the execution of various stability tests, the robust minimum 1D-velocity model implies guiding values for the composition of the local, subsurface structure of the crust. Afterwards, performing a hypocenter relocation enables the assignment of individual earthquakes to aftershock clusters or extended seismotectonic structures. This allows the mapping of previously unknown seismogenic faults.
Finally, focal mechanisms are modeled for events with acurately located hypocenters, using the newly derived local velocity model. A compressive regime is attested by the majority of focal mechanisms, while the strike direction of the individual seismogenic structures is in agreement with the overall north – south orientation of the Central Andes, its mountain front, and individual mountain ranges in the southern Santa-Bárbara-System.
Millennial-scale palaeoclimate variability has been documented in various terrestrial and marine palaeoclimate proxy records throughout the Northern Hemisphere for the last glacial cycle. Its clear expression and rapid shifts between different states of climate (Greenland Interstadials and Stadials) represents a correlation tool beyond the resolution of e.g. luminescence dating, especially relevant for terrestrial deposits. Usually, comparison of terrestrial proxy datasets and the Greenland ice cores indicates a complex expression of millennial-scale climate variability as recorded in terrestrial geoarchives including loess. Loess is the most widespread terrestrial geoarchive of the Quaternary and especially widespread over Eurasia. However, loess often records a smoothed representation of millennial-scale variability without all fidelity when compared to the Greenland data, this being a relevant limiting feature in integrating loess with other palaeoclimate records. To better understand the loess proxy-response to millennial-scale climate variability, we simulate a proxy signal smoothing by natural processes through application of low-pass filters of delta O-18 data from Greenland, a high-resolution palaeoclimate reference record, alongside speleothem isotope records from the Black Sea-Mediterranean region. We show that low-pass filters represent rather simple models for better constraining the expression of millennial-scale climate variability in low sedimentation environments, and in sediments where proxy-response signals are most likely affected by natural smoothing (by e.g. bioturbation). Interestingly, smoothed datasets from Greenland and the Black Sea-Mediterranean region are most similar in the last similar to 15 ka and between similar to 50-30 ka. Between similar to 30-15 ka, roughly corresponding to the Last Glacial Maximum and the deglaciation, the records show dissimilarities, challenging the construction of robust correlative time-scales in this age range. From our analysis it becomes apparent that patterns of palaeoclimate signals in loess-palaeosol sequences often might be better explained by smoothed Greenland reference data than the original high-resolution Greenland dataset, or other reference data. This opens the possibility to better assess the temporal resolution and palaeoclimate potential of loess-palaeosol sequences in recording supra-regional climate patterns, as well as to securely integrate loess with other chronologically better-resolved palaeoclimate records.
Millennial-scale palaeoclimate variability has been documented in various terrestrial and marine palaeoclimate proxy records throughout the Northern Hemisphere for the last glacial cycle. Its clear expression and rapid shifts between different states of climate (Greenland Interstadials and Stadials) represents a correlation tool beyond the resolution of e.g. luminescence dating, especially relevant for terrestrial deposits. Usually, comparison of terrestrial proxy datasets and the Greenland ice cores indicates a complex expression of millennial-scale climate variability as recorded in terrestrial geoarchives including loess. Loess is the most widespread terrestrial geoarchive of the Quaternary and especially widespread over Eurasia. However, loess often records a smoothed representation of millennial-scale variability without all fidelity when compared to the Greenland data, this being a relevant limiting feature in integrating loess with other palaeoclimate records. To better understand the loess proxy-response to millennial-scale climate variability, we simulate a proxy signal smoothing by natural processes through application of low-pass filters of delta O-18 data from Greenland, a high-resolution palaeoclimate reference record, alongside speleothem isotope records from the Black Sea-Mediterranean region. We show that low-pass filters represent rather simple models for better constraining the expression of millennial-scale climate variability in low sedimentation environments, and in sediments where proxy-response signals are most likely affected by natural smoothing (by e.g. bioturbation). Interestingly, smoothed datasets from Greenland and the Black Sea-Mediterranean region are most similar in the last similar to 15 ka and between similar to 50-30 ka. Between similar to 30-15 ka, roughly corresponding to the Last Glacial Maximum and the deglaciation, the records show dissimilarities, challenging the construction of robust correlative time-scales in this age range. From our analysis it becomes apparent that patterns of palaeoclimate signals in loess-palaeosol sequences often might be better explained by smoothed Greenland reference data than the original high-resolution Greenland dataset, or other reference data. This opens the possibility to better assess the temporal resolution and palaeoclimate potential of loess-palaeosol sequences in recording supra-regional climate patterns, as well as to securely integrate loess with other chronologically better-resolved palaeoclimate records.
The Frasnian-Famennian (F-F) transition of Late Devonian was a critical episode in geological history, recording a major mass extinction event. In this study, we focus on an F-F succession from a deep marine context in Bancheng, southern Guangxi, South China, to investigate coeval changes in pelagic environments of the Paleo-Tethys Ocean. The studied succession is exclusively composed of bedded cherts intercalated with multiple siliceous volcanic ash beds. A SIMS zircon U-Pb Concordia age of 367.8 +/- 2.5 Ma is reported for a tuffaceous layer slightly above the F-F boundary. Geochemical ratios of Al/(Al + Fe + Mn), Ce/Ce*, Y/Ho, and Al, Fe contents in bedded cherts indicate that they are of predominantly biogenic/chemical origin with some terrigenous inputs. Negligible enrichment of redox sensitive elements (Mo, U, V) and low V/Cr ratios (<2) suggest persistently oxic conditions existed in the deep pelagic basin at Bancheng, South China during the F-F transition. These findings call into question the widely held hypothesis that marine anoxia was the primary killing mechanism for the F-F crisis. In contrast, multiple tuffaceous layers throughout the F-F boundary succession indicate frequent volcanic activity, which could have released massive amounts of greenhouse gases into the atmosphere, inducing climate warming. This scenario may have increased continental weathering and riverine fluxes into the ocean, reconciling the increases in Al2O3 content and Al/(Al + Fe + Mn) ratio across the F-F boundary. Documentation of persistently oxic conditions and frequent volcanic activitiy provides new perspectives on the inter-relationship between volcanism, climate, and oceanic redox fluctuation during the F-F biotic crisis.
Increasingly available high-frequency data during storm events, when hydrological dynamics most likely activate nitrate storage-flux exchanges, reveal insights into catchment nitrate dynamics. In this study, we explored impacts of seasonality and landscape gradients on nitrate concentration-discharge (C-Q) hysteresis patterns in the Selke catchment, central Germany, which has heterogeneous combinations of meteorological, hydrogeological and land use conditions. Three nested gauging stations established along the main Selke River captured flow and nitrate export dynamics from the uppermost subcatchment (mixed forest and arable land), middle subcatchment (pure steep forest) and lowermost subcatchment (arable and urban land). We collected continuous high-frequency (15-min) discharge and nitrate concentration data from 2012 to 2017 and analyzed the 223 events detected at all three stations. A dominant hysteresis pattern in the uppermost and middle subcatchments was counter-clockwise and combined with an accretion effect, indicating many proximal and mobilized distal nitrate sources. However, 66% of all events at the catchment outlet experienced a dilution effect, possibly due to mechanisms that vary seasonally. During wetting/wet periods (October-March), it was combined mainly with a counter-clockwise pattern due to the dominance of event runoff volume from the uppermost and middle subcatchments. During drying/dry periods (April-September), however, it was combined mainly with a clockwise pattern due to occasional quick surface flows from lowland near-stream urban areas. In addition, the clockwise hysteresis occurred mainly from May-October during mostly drying/dry periods at all three sites, indicating little distal nitrate transport in response to the low terrestrial hydrological connectivity, especially in the lowermost dry and flat sub-catchment. This comprehensive analysis (i.e., clockwise vs. counter-clockwise, accretion vs. dilution) enables in-depth analysis of nitrate export mechanisms during certain periods under different landscape conditions. Specific combination of C-Q relationships could identify target locations for agricultural management actions that decrease nitrate output. Therefore, we strongly encourage long-term multisite and high-frequency monitoring strategies in heterogeneous nested catchment(s), which can help understand process mechanisms, generate data for physical-based water-quality modeling and provide guidance for water and agricultural management.
The Indus Molasse records orogenic sedimentation associated with uplift and erosion of the southern margin of Asia in the course of ongoing India-Eurasia collision. Detailed field investigation clarifies the nature and extent of the depositional contact between this molasse and the underlying basement units. We report the first dataset on detrital zircon U-Pb ages, Hf isotopes and apatite U-Pb ages for the autochthonous molasse in the Indus Suture Zone. A latest Oligocene depositional age is proposed on the basis of the youngest detrital zircon U-Pb age peak and is consistent with published biostratigraphic data. Multiple provenance indicators suggest exclusively northerly derivation with no input from India in the lowermost parts of the section. The results provide constraints on the uplift and erosion history of the Ladakh Range following the initial India-Asia collision.
In this investigation, we examine the uncertainties using the horizontal-to-vertical spectral ratio (HVSR) technique on earthquake recordings to detect site resonant frequencies at 207 KiK-net sites. Our results show that the scenario dependence of response (pseudospectral acceleration) spectral ratio could bias the estimates of resonant frequencies for sites having multiple significant peaks with comparable amplitudes. Thus, the Fourier amplitude spectrum (FAS) should be preferred in computing HVSR. For more than 80% of the investigated sites, the first peak (in the frequency domain) on the average HVSR curve over multiple sites coincides with the highest peak. However, for sites with multiple peaks, the highest peak frequency (f(p)) is less susceptible to the selection criteria of significant peaks and the extent of smoothing to spectrum than the first peak frequency (f(0)). Meanwhile, in comparison to the surface-to-borehole spectral ratio, f(0) tends to underestimate the predominant frequency (at which the largest amplification occurs) more than f(p). In addition, in terms of characterizing linear site response, f(p) shows a better overall performance than f(0). Based on these findings, we thus recommend that seismic network operators provide f(p) on the average HVSRFAS curve as a priority, ideally together with the average HVSRFAS curve in site characterization.
Evaluation of a novel application of earthquake HVSR in site-specific amplification estimation
(2020)
Ground response analyses (GRA) model the vertical propagations of SH waves through flat-layered media (1DSH) and are widely carried out to evaluate local site effects in practice. Horizontal-to-vertical spectral ratio (HVSR) technique is a cost-effective approach to extract certain site-specific information, e.g., site fundamental frequency (f(0)), but HVSR values cannot be directly used to approximate the levels of S-wave amplifications. Motivated by the work of Kawase et al. (2019), we propose a procedure to correct earthquake HVSR amplitudes for direct amplification estimations. The empirical correction compensates HVSR by generic vertical amplification spectra categorized by the vertical fundamental frequency (f(0v)) via kappa-means clustering. In this investigation, we evaluate the effectiveness of the corrected HVSR in approximating observed linear amplifications in comparison with 1DSH modellings. We select a total of 90 KiK-net (Kiban Kyoshin network) surface-downhole sites which are found to have no velocity contrasts below their boreholes and thus of which surface-to-borehole spectral ratios (SBSRs) can be taken as their empirical transfer functions (ETFs). 1DSH-based theoretical transfer functions (TTFs) are computed in the linear domain considering uncertainties in Vs profiles through randomizations. Five goodness-of-fit metrics are adopted to gauge the closeness between observed (ETF) and predicted (i.e., TTF and corrected HVSR) amplifications in both amplitude and spectral shape over frequencies from f(0) to 25 Hz. We find that the empirical correction to HVSR is highly effective and achieves a "good match" in both spectral shape and amplitude at the majority of the 90 KiK-net sites, as opposed to less than one-third for the 1DSH modelling. In addition, the empirical correction does not require a velocity model, which GRAs require, and thus has great potentials in seismic hazard assessments.