Simulating and mapping the risk of surge floods in multiple typhoon scenarios: a case study of Yuhuan County,Zhejiang Province,China

We propose a scenario-based method for simulating and mapping the risk of surge floods for use by local authorities concerned with public safety and urban planning in coastal areas. Focusing on the triad of hazard, vulnerability and adaptation capability, we estimate the comprehensive risk and display its spatial distribution using the raster calculation tool in ArcGIS. The detailed methodology is introduced via a case study of Yuhuan, an island county in Zhejiang Province, China, which is frequently affected by typhoon storm surges. First, we designed 24 typhoon scenarios and modeled the flood process in each scenario using the hydrodynamic module of MIKE 21. Second, flood depth and area were used for hazard assessment; an authorized indicator system of land use categories and a survey of emergency shelters were used for vulnerability and adaptation capability assessment, respectively; and a quantified model was used for assessment of the comprehensive risk. Lastly, we used the GIS raster calculation tool for mapping the risk of storm surges in multiple typhoon scenarios. Our principal findings are as follows: (1) Seawalls are more likely to be overtopped or destroyed by more severe storm surges with increasing typhoon intensity. (2) Most of the residential areas with inadequate emergency shelters are highly vulnerable to flood events. (3) As projected in the risk mapping, if an exceptional typhoon with a central pressure of 915 or 925 hPa made a landfall in Yuhuan, a wide range of areas would be flooded and at high risk. (4) Determining optimal strategies based on identification of risk-inducing factors is the most effective way of promoting safe and sustainable development in coastal cities.     

Sea-surge and wind speed extremes: optimal estimation strategies for planners and engineers

Accurate and precise estimation of return levels is often a key goal of any extreme value analysis. For example, in the UK the British Standards Institution (BSI) incorporate estimates of ‘once-in-50-year wind gust speeds’—or 50-year return levels—into their design codes for new structures; similarly, the Dutch Delta Commission use estimates of the 10,000-year return level for sea-surge to aid the construction of flood defence systems. In this paper, we briefly highlight the shortcomings of standard methods for estimating return levels, including the commonly-adopted block maxima and peaks over thresholds approach, before presenting an estimation framework which we show can substantially increase the precision of return level estimates. Our work allows explicit quantification of seasonal effects, as well as exploiting recent developments in the estimation of the extremal index for handling extremal clustering. From frequentist ideas, we turn to the Bayesian paradigm as a natural approach for building complex hierarchical or spatial models for extremes. Through simulations we show that the return level posterior mean does not have an exceedance probability in line with the intended encounter risk; we also argue that the Bayesian posterior predictive value gives the most satisfactory representation of a return level for use in practice, accounting for uncertainty in parameter estimation and future observations. Thus, where feasible, we propose a Bayesian estimation strategy for optimal return level inference.     

Ground-motion models for average spectral acceleration in a period range: direct and indirect methods

Average spectral acceleration, AvgSA, is defined as the geometric mean of spectral acceleration values over a range of periods and it is a ground motion intensity measure used for structural response prediction. One of its advantages stands on the assumption that its distribution is computable from the available GMPEs for spectral acceleration, GMPE-SA, (called here indirect method) without the need for deriving new specific GMPEs for AvgSA, GMPE-AvgSA, (called here direct method). To what extent this assumption is valid, however, has never been verified. As such, we derived an empirical GMPE-AvgSA based on RESORCE ground motion dataset and we compared its predicted values with those from a GMPE-SA via the indirect approach. As expected, the results show that the indirect approach yields median AvgSA estimates that are identical to those of the direct approach. However, the estimates of AvgSA variance of the two methods are identical only if both the GMPE-SA and their empirical correlation coefficients among different SA ordinates are derived from the same record dataset.     

Sampling uncertainties and source <Emphasis Type="Italic">b</Emphasis> likelihood for the Gutenberg-Richter <Emphasis Type="Italic">b</Emphasis> value from the Aki-Utsu method

The Aki-Utsu method of Gutenberg-Richter (G-R) b value estimation is often misapplied so that estimations not using the G-R histogram are often meaningless because they are not based on adequate samples. We propose a method to estimate the likelihood Pr(b?b _m , N, M ₁, M ₂) that an observed b _m estimate, based on a sample of N magnitudes within an [M ₁????≤?ΔM/2,?M ₂?+?ΔM/2) range, where ΔM?=?0.1 is the usual rounding applied to magnitudes, is due to a “true” source b value, b, and use these likelihoods to estimate source b ranges corresponding to various confidence levels. As an example of application of the method, we estimate the b values before and after the occurrence of a 7.4-magnitude earthquake in the Mexican subduction zone, and find a difference of 0.82 between them with 100% confidence that the b values are different.     

Dimensionality reduction for efficient Bayesian estimation of groundwater flow in strongly heterogeneous aquifers

We focus on the Bayesian estimation of strongly heterogeneous transmissivity fields conditional on data sampled at a set of locations in an aquifer. Log-transmissivity, Y, is modeled as a stochastic Gaussian process, parameterized through a truncated Karhunen–Loève (KL) expansion. We consider Y fields characterized by a short correlation scale as compared to the size of the observed domain. These systems are associated with a KL decomposition which still requires a high number of parameters, thus hampering the efficiency of the Bayesian estimation of the underlying stochastic field. The distinctive aim of this work is to present an efficient approach for the stochastic inverse modeling of fully saturated groundwater flow in these types of strongly heterogeneous domains. The methodology is grounded on the construction of an optimal sparse KL decomposition which is achieved by retaining only a limited set of modes in the expansion. Mode selection is driven by model selection criteria and is conditional on available data of hydraulic heads and (optionally) Y. Bayesian inversion of the optimal sparse KLE is then inferred using Markov Chain Monte Carlo (MCMC) samplers. As a test bed, we illustrate our approach by way of a suite of computational examples where noisy head and Y values are sampled from a given randomly generated system. Our findings suggest that the proposed methodology yields a globally satisfactory inversion of the stochastic head and Y fields. Comparison of reference values against the corresponding MCMC predictive distributions suggests that observed values are well reproduced in a probabilistic sense. In a few cases, reference values at some unsampled locations (typically far from measurements) are not captured by the posterior probability distributions. In these cases, the quality of the estimation could be improved, e.g., by increasing the number of measurements and/or the threshold for the selection of KL modes.     

On the possible correlation between the Gutenberg-Richter parameters of the frequency-magnitude relationship

The study of the Gutenberg-Richter (GR) parameters a and b has been very important to describe and characterize the seismicity over the different seismic provinces around the world. As far as we know, the possible correlation between the GR parameters a and b has not received enough attention. Bayrak et al. reported the a and b values for 27 active seismic regions around the boundaries of the main tectonic plates of the world. From these data, we found that there exists a positive correlation between the a and b parameters (R =?0.85, R² =?0.72). On the other hand, we made around 150 computer runs of a spring-block model proposed by Olami et al. (Phys Rev Lett 68(8):1244–1247, 1992). This model roughly emulates the interaction between two fault planes and it reaches a self-organized critical state. With these simulations, we also found that the a and b parameters are positively correlated. Motivated by these results, we propose an analytical demonstration that indeed a and b are positively correlated. In addition, we discuss on other possible applications of the spring-block model to actual seismicity and to frictional experiments made with sandpapers.     

Effect of zonal E × B plasma drift on electron density in the low-latitude ionospheric <Emphasis Type="Italic">F</Emphasis> region at a solar activity maximum near vernal equinox

The variations in the density of the ionospheric F2 layer maximum (NmF2) under the action of the zonal plasma drift perpendicularly to the magnetic (B) and electric (E) fields in the direction geomagnetic west-geomagnetic east have been studied using the three-dimensional nonstationary theoretical model of electron and ion densities (N _e and N _i ) and temperatures (T _e and T _i ) in the low-latitude and midlatitude ionospheric F region and plasmasphere. The method of numerical calculations of N _e , N _i , T _e , and T _i , including the advantages of the Lagrangian and Eulerian methods, is used in the model. A dipole approximation of the geomagnetic field (B), taking into account the non-coincidence of the geographic and geomagnetic poles and differences between the positions of the Earth’s and geomagnetic dipole centers, is accepted in the calculations. The calculated NmF2 and altitudes of the F2 layer maximum (hmF2) have been compared with these quantities measured at 16 low-latitude ionospheric sounding stations during the geomagnetically quiet period October 11–12, 1958. This comparison made it possible to correct the input model parameters: the NRLMSISE-00 model [O], the meridional component of the neutral wind velocity according to the HWW90 model, and the meridional component of the equatorial plasma drift due to the electric field specified by the empirical model. It has been indicated that the effect of the zonal E × B plasma drift on NmF2 can be neglected under daytime conditions and changes in NmF2 and hmF2 under the action of this drift are insignificant under nighttime conditions north of 25° and south of ?26° geomagnetic latitude. The effect of the zonal E × B plasma drift on NmF2 and hmF2 is most substantial in the nightside ionosphere approximately from ?20° to 20° geomagnetic latitude, and the neglect of this drift results in an up to 2.4-fold underestimation of NmF2. The found dependence of the effect of the zonal E × B plasma drift on NmF2 and hmF2 on geomagnetic latitude is related to the longitudinal asymmetry of B, asymmetry of the neutral wind about the geomagnetic equator, and changes in the meridional E × B plasma drift at a change in geomagnetic longitude.     

Improved random noise attenuation using <Emphasis Type="Italic">f</Emphasis>−<Emphasis Type="Italic">x</Emphasis> empirical mode decomposition and local similarity

Conventional f?x empirical mode decomposition (EMD) is an effective random noise attenuation method for use with seismic profiles mainly containing horizontal events. However, when a seismic event is not horizontal, the use of f?x EMD is harmful to most useful signals. Based on the framework of f?x EMD, this study proposes an improved denoising approach that retrieves lost useful signals by detecting effective signal points in a noise section using local similarity and then designing a weighting operator for retrieving signals. Compared with conventional f?x EMD, f?x predictive filtering, and f?x empirical mode decomposition predictive filtering, the new approach can preserve more useful signals and obtain a relatively cleaner denoised image. Synthetic and field data examples are shown as test performances of the proposed approach, thereby verifying the effectiveness of this method.     

Limits of transversely isotropic elastic parameters for the existence of classical Rayleigh waves

Solutions of P-SV equations of motion in a homogeneous transversely isotropic elastic layer contain a factor exp(±ν _j z), where z is the vertical coordinate and j?=?1, 2. For computing Rayleigh wave dispersion in a multi-layered half space, ν _j is computed at each layer. For a given phase velocity (c), ν _j becomes complex depending on the transversely isotropic parameters. When ν _j is complex, classical Rayleigh waves do not exist and generalised Rayleigh waves propagate along a path inclined to the interface. We use transversely isotropic parameters as α _H , β _V , ξ, ? and η and find their limits beyond which ν _j becomes complex. It is seen that ν _j depends on ? and η, but does not depend on ξ. The complex ν _j occurs when ? is small and η is large. For a given c/β _V , the region of complex ν _j in a ? -η plane increases with the increase of α _H /β _V . Further, for a given α _H /β _V , the complex region of ν _j increases significantly with the decrease of c/β _V . This study is useful to compute dispersion parameters of Rayleigh waves in a layered medium.     

Shear wave velocity-based evaluation and design of stone column improved ground for liquefaction mitigation

The evaluation and design of stone column improvement ground for liquefaction mitigation is a challenging issue for the state of practice. In this paper, a shear wave velocity-based approach is proposed based on the well-defined correlations of liquefaction resistance (CRR)-shear wave velocity (V _s)-void ratio (e) of sandy soils, and the values of parameters in this approach are recommended for preliminary design purpose when site specific values are not available. The detailed procedures of pre- and post-improvement liquefaction evaluations and stone column design are given. According to this approach, the required level of ground improvement will be met once the target V _s of soil is raised high enough (i.e., no less than the critical velocity) to resist the given earthquake loading according to the CRR-V _s relationship, and then this requirement is transferred to the control of target void ratio (i.e., the critical e) according to the V _s-e relationship. As this approach relies on the densification of the surrounding soil instead of the whole improved ground and is conservative by nature, specific considerations of the densification mechanism and effect are given, and the effects of drainage and reinforcement of stone columns are also discussed. A case study of a thermal power plant in Indonesia is introduced, where the effectiveness of stone column improved ground was evaluated by the proposed V _s-based method and compared with the SPT-based evaluation. This improved ground performed well and experienced no liquefaction during subsequent strong earthquakes.     

The spatiotemporal analysis of the minimum magnitude of completeness <Emphasis Type="Italic">Mc</Emphasis> and the Gutenberg–Richter law <Emphasis Type="Italic">b</Emphasis>-value parameter using the earthquake catalog of Greece

Spatiotemporal mapping the minimum magnitude of completeness Mc and b-value of the Gutenberg–Richter law is conducted for the earthquake catalog data of Greece. The data were recorded by the seismic network of the Institute of Geodynamics of the National Observatory of Athens (GINOA) in 1970–2010 and by the Hellenic Unified Seismic Network (HUSN) in 2011–2014. It is shown that with the beginning of the measurements at HUSN, the number of the recorded events more than quintupled. The magnitude of completeness Mc of the earthquake catalog for 1970–2010 varies within 2.7 to 3.5, whereas starting from April 2011 it decreases to 1.5–1.8 in the central part of the region and fluctuates around the average of 2.0 in the study region overall. The magnitude of completeness Mc and b-value for the catalogs of the earthquakes recorded by the old (GINOA) and new (HUSN) seismic networks are compared. It is hypothesized that the magnitude of completeness Mc may affect the b-value estimates. The spatial distribution of the b-value determined from the HUSN catalog data generally agrees with the main geotectonic features of the studied territory. It is shown that the b-value is below 1 in the zones of compression and is larger than or equal to 1 in the zones dominated by extension. The established depth dependence of the b-value is pretty much consistent with the hypothesis of a brittle–ductile transition zone existing in the Earth’s crust. It is assumed that the source depth of a strong earthquake can probably be estimated from the depth distribution of the b-value, which can be used for seismic hazard assessment.     

Peculiarities of Rayleigh wave displacement in a transversely isotropic half space

We consider a transversely isotropic medium with vertical axis of symmetry (VTI). Rayleigh wave displacement components in a homogeneous VTI medium contain exp(±kr_jz), where z is the vertical coordinate, k is the wave number, and j?=?1, 2; r_j may be considered as depth-decay factor. In a VTI medium, r_j can be a real or imaginary as in an isotropic medium, or it can be a complex depending on the elastic parameters of the VTI medium; if complex, r₁ and r₂ are complex conjugates. In a homogeneous VTI half space, Rayleigh wave displacement is significantly modified with a phase shift when r_j changes from real to complex with variation of VTI parameters; at the transition, the displacement becomes zero when r₁?=?r₂. In a liquid layer over a VTI half space, Rayleigh waves are dispersive. Here, also Rayleigh wave displacement significantly modified with a phase shift when r_j changes from real to complex with a decrease of period. At very low period, phase velocity of Rayleigh waves becomes less than P-wave velocity in the liquid layer giving rise to Scholte waves (interface waves). The amplitudes of Scholte waves with a VTI half space are found to be significantly larger than those with an isotropic half space.     

Mechanism of formation of <Emphasis Type="Italic">Q</Emphasis>-disturbances in the <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> region of the equatorial ionosphere

Using model simulations, the morphological picture (revealed earlier) of the disturbances in the F ₂ region of the equatorial ionosphere under quiet geomagnetic conditions (Q-disturbances) is interpreted. It is shown that the observed variations in the velocity of the vertical E × B plasma drift, related to the zonal E _y component of the electric field, are responsible for the formation of Q-disturbances. The plasma recombination at altitudes of the lower part of the F ₂ region and the dependence of the rate of this process on heliogeophysical conditions compose the mechanism of Q-disturbance formation at night. The daytime positive Q-disturbances are caused exclusively by a decrease in the upward E × B drift, and this type of disturbances could be related to the known phenomenon of counter electrojet. Possible causes of formation of the daytime negative Q-disturbances are discussed.     

Global empirical model of critical frequency of the ionospheric <Emphasis Type="Italic">F</Emphasis>2-layer for quiet geomagnetic conditions

Using the foF2 database obtained from satellites and ground-based ionospheric stations, we have constructed a global empirical model of the critical frequency of the ionospheric F2-layer (SDMF2—Satellite and Digisonde Data Model of the F2 layer) for quiet geomagnetic conditions (Kp < 3). The input parameters of this model are the geographical coordinates, UT, day, month, year, and the integral index F10.7 (day, τ = 0.96) of solar activity for a given day. The SDMF2 model was based on the Legendre method for the spatial expansion of foF2 monthly medians to 12 in latitude and 8 in longitude of spherical harmonics. The resulting spatial coefficients have been expanded by the Fourier method in three spherical harmonics with respect to UT. The effect of the saturation of critical frequency of the ionospheric F2-layer at high solar activity was described in the SDMF2 model by foF2 as a logarithmic function of F10.7 (day, τ = 0.96). The difference between the SDMF2 and IRI models is a maximum at low solar activity as well as in the Southern Hemisphere and in the oceans. The testing on the basis of ground-based and satellite data has indicated that the SDMF2 model is more accurate than the IRI model.     

Understanding single-station ground motion variability and uncertainty (sigma): lessons learnt from EUROSEISTEST

Accelerometric data from the well-studied valley EUROSEISTEST are used to investigate ground motion uncertainty and variability. We define a simple local ground motion prediction equation (GMPE) and investigate changes in standard deviation (σ) and its components, the between-event variability (τ) and within-event variability (φ). Improving seismological metadata significantly reduces τ (30–50%), which in turn reduces the total σ. Improving site information reduces the systematic site-to-site variability, φ _S2S (20–30%), in turn reducing φ, and ultimately, σ. Our values of standard deviations are lower than global values from literature, and closer to path-specific than site-specific values. However, our data have insufficient azimuthal coverage for single-path analysis. Certain stations have higher ground-motion variability, possibly due to topography, basin edge or downgoing wave effects. Sensitivity checks show that 3 recordings per event is a sufficient data selection criterion, however, one of the dataset’s advantages is the large number of recordings per station (9–90) that yields good site term estimates. We examine uncertainty components binning our data with magnitude from 0.01 to 2 s; at smaller magnitudes, τ decreases and φ _SS increases, possibly due to κ and source-site trade-offs Finally, we investigate the alternative approach of computing φ _SS using existing GMPEs instead of creating an ad hoc local GMPE. This is important where data are insufficient to create one, or when site-specific PSHA is performed. We show that global GMPEs may still capture φ _SS , provided that: (1) the magnitude scaling errors are accommodated by the event terms; (2) there are no distance scaling errors (use of a regionally applicable model). Site terms (φ _S2S ) computed by different global GMPEs (using different site-proxies) vary significantly, especially for hard-rock sites. This indicates that GMPEs may be poorly constrained where they are sometimes most needed, i.e., for hard rock.     

Validity of paleomagnetic poles and principles of constructing their wander paths: A case study of the East European platform

Methods of selecting paleomagnetic data for the construction of apparent polar wander paths (APWPs) are analyzed. It is shown that the existing criteria of reliability of paleomagnetic data cannot be regarded as evidence for their validity. In other words, no unambiguous dependence exists between the reliability and the closeness of paleomagnetic poles to a hypothetical region crossed by the reliable APWP. A new approach to the construction of paleomagnetic APWPs based on simple principles (principle of space and principle of time) is proposed. Using a numerical implementation of this algorithm, three stable clusters were determined (L _p = 164, F _p = 43; L _p = 144, F _p = 13; and L _p = 170, F _p = ?2); the respective maximum estimates of their ages are 248–251, 345, and 385 Ma. These clusters can be regarded as reliable paleomagnetic poles in the Paleozoic of the East European platform.     

Stochastic models of severe weather watches and warnings: transition probabilities

To alert the public to the possibility of tornado (T), hail (H), or convective wind (C), the National Weather Service (NWS) issues watches (V) and warnings (W). There are severe thunderstorm watches (SV), tornado watches (TV), and particularly dangerous situation watches (PV); and there are severe thunderstorm warnings (SW), and tornado warnings (TW). Two stochastic models are formulated that quantify uncertainty in severe weather alarms for the purpose of making decisions: a one-stage model for deciders who respond to warnings, and a two-stage model for deciders who respond to watches and warnings. The models identify all possible sequences of watches, warnings, and events, and characterize the associated uncertainties in terms of transition probabilities. The modeling approach is demonstrated on data from the NWS Norman, Oklahoma, warning area, years 2000–2007. The major findings are these. (i) Irrespective of its official designation, every warning type {SW, TW} predicts with a significant probability every event type {T, H, C}. (ii) An ordered intersection of SW and TW, defined as reinforced warning (RW), provides additional predictive information and outperforms SW and TW. (iii) A watch rarely leads directly to an event, and most frequently is false. But a watch that precedes a warning does matter. The watch type \(\{SV\), TV, \(PV\}\) is a predictor of the warning type \(\{SW\), RW, \(TW\}\) and of the warning performance: It sharpens the false alarm rate of the warning and the predictive probability of an event, and it increases the average lead time of the warning.     

Regression between earthquake magnitudes having errors with known variances

Recent publications on the regression between earthquake magnitudes assume that both magnitudes are affected by error and that only the ratio of error variances is known. If X and Y represent observed magnitudes, and x and y represent the corresponding theoretical values, the problem is to find the a and b of the best-fit line \(y = a x + b\). This problem has a closed solution only for homoscedastic errors (their variances are all equal for each of the two variables). The published solution was derived using a method that cannot provide a sum of squares of residuals. Therefore, it is not possible to compare the goodness of fit for different pairs of magnitudes. Furthermore, the method does not provide expressions for the x and y. The least-squares method introduced here does not have these drawbacks. The two methods of solution result in the same equations for a and b. General properties of a discussed in the literature but not proved, or proved for particular cases, are derived here. A comparison of different expressions for the variances of a and b is provided. The paper also considers the statistical aspects of the ongoing debate regarding the prediction of y given X. Analysis of actual data from the literature shows that a new approach produces an average improvement of less than 0.1 magnitude units over the standard approach when applied to \(M_{w}\) vs. \(m_{b}\) and \(M_{w}\) vs. \(M_{S}\) regressions. This improvement is minor, within the typical error of \(M_{w}\). Moreover, a test subset of 100 predicted magnitudes shows that the new approach results in magnitudes closer to the theoretically true magnitudes for only 65 % of them. For the remaining 35 %, the standard approach produces closer values. Therefore, the new approach does not always give the most accurate magnitude estimates.     

A new <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> estimation parameter for use in earthquake early warning systems

We propose a method that employs the squared displacement integral (ID2) to estimate earthquake magnitudes in real time for use in earthquake early warning (EEW) systems. Moreover, using τ _c and P _d for comparison, we establish formulas for estimating the moment magnitudes of these three parameters based on the selected aftershocks (4.0 ≤ M _s  ≤ 6.5) of the 2008 Wenchuan earthquake. In this comparison, the proposed ID2 method displays the highest accuracy. Furthermore, we investigate the applicability of the initial parameters to large earthquakes by estimating the magnitude of the Wenchuan M _s 8.0 mainshock using a 3-s time window. Although these three parameters all display problems with saturation, the proposed ID2 parameter is relatively accurate. The evolutionary estimation of ID2 as a function of the time window shows that the estimation equation established with ID2 ^Ref determined from the first 8-s of P wave data can be directly applicable to predicate the magnitudes of 8.0. Therefore, the proposed ID2 parameter provides a robust estimator of earthquake moment magnitudes and can be used for EEW purposes.     

Limitations of <Emphasis Type="Italic">H</Emphasis>-<Emphasis Type="Italic">κ</Emphasis> stacking: ambiguous results caused by crustal layering

Over the past decade, the H-κ stacking technique of Zhu and Kanamori (J Geophys Res 105:2969–2980, 2000) has become a standard tool to determine the crustal thickness H and the bulk crustal v_P/v_S ratio κ from teleseismic receiver functions. It is obvious that unfavorable noise conditions as well as a complex 3D velocity structure can severely hamper the interpretation of receiver-function data. However, we observe that ambiguities can even arise from a simple 1D layered velocity structure which raises a high potential for misinterpretations. To analyze the feasibility and basic limitations of the H-κ stacking method, we conduct a series of tests based on synthetic data. The impact of different given elementary parameters, related either to the velocity structure or to the data processing, is evaluated in a series of eight individual tests. We deliberately exclude complications such as 3D structural variations and/or noise to show that even a simple 1D velocity structure, involving, e.g., an additional inter-crustal discontinuity, can have significant consequences for the interpretation of the results. However, our modeling suggests that more complex crustal structures may lead to even less reliable results. Additionally, our tests illustrate that time shifts of the maxima in the H-κ domain due to the superposition and merging of individual phases can lead to significantly overestimated v_P/v_S ratios. In general, the depth to the Moho (or other discontinuities of interest) is less significantly affected. Our tests indicate the necessity to accurately check delay times derived from the maxima of the H-κ stacks against corresponding phases in the receiver functions. Repeating the stacking with varied weighting factors and filter ranges can help to reduce the ambiguities and to avoid possible misinterpretation.