TY - JOUR A1 - Berner, Nadine A1 - Trauth, Martin H. A1 - Holschneider, Matthias T1 - Bayesian inference about Plio-Pleistocene climate transitions in Africa JF - Quaternary science reviews : the international multidisciplinary research and review journal N2 - During the last 5 Ma the Earth's ocean-atmosphere system passed through several major transitions, many of which are discussed as possible triggers for human evolution. A classic in this context is the possible influence of the closure of the Panama Strait, the intensification of Northern Hemisphere Glaciation, a stepwise increase in aridity in Africa, and the first appearance of the genus Homo about 2.5 - 2.7 Ma ago. Apart from the fact that the correlation between these events does not necessarily imply causality, many attempts to establish a relationship between climate and evolution fail due to the challenge of precisely localizing an a priori unknown number of changes potentially underlying complex climate records. The kernel-based Bayesian inference approach applied here allows inferring the location, generic shape, and temporal scale of multiple transitions in established records of Plio-Pleistocene African climate. By defining a transparent probabilistic analysis strategy, we are able to identify conjoint changes occurring across the investigated terrigenous dust records from Ocean Drilling Programme (ODP) sites in the Atlantic Ocean (ODP 659), Arabian (ODP 721/722) and Mediterranean Sea (ODP 967). The study indicates a two-step transition in the African climate proxy records at (2.35-2.10) Ma and (1.70 - 1.50) Ma, that may be associated with the reorganization of the Hadley-Walker Circulation. . KW - Plio-Pleistocene KW - Hadley-Walker Circulation KW - climate transition KW - Bayesian inference KW - time series analysis KW - ODP 659 KW - ODP 721/722 KW - ODP 967 Y1 - 2022 U6 - https://doi.org/10.1016/j.quascirev.2021.107287 SN - 0277-3791 SN - 1873-457X VL - 277 PB - Elsevier CY - Oxford ER - TY - JOUR A1 - Sharma, Shubham A1 - Hainzl, Sebastian A1 - Zöller, Gert A1 - Holschneider, Matthias T1 - Is Coulomb stress the best choice for aftershock forecasting? JF - Journal of geophysical research : Solid earth N2 - The Coulomb failure stress (CFS) criterion is the most commonly used method for predicting spatial distributions of aftershocks following large earthquakes. However, large uncertainties are always associated with the calculation of Coulomb stress change. The uncertainties mainly arise due to nonunique slip inversions and unknown receiver faults; especially for the latter, results are highly dependent on the choice of the assumed receiver mechanism. Based on binary tests (aftershocks yes/no), recent studies suggest that alternative stress quantities, a distance-slip probabilistic model as well as deep neural network (DNN) approaches, all are superior to CFS with predefined receiver mechanism. To challenge this conclusion, which might have large implications, we use 289 slip inversions from SRCMOD database to calculate more realistic CFS values for a layered half-space and variable receiver mechanisms. We also analyze the effect of the magnitude cutoff, grid size variation, and aftershock duration to verify the use of receiver operating characteristic (ROC) analysis for the ranking of stress metrics. The observations suggest that introducing a layered half-space does not improve the stress maps and ROC curves. However, results significantly improve for larger aftershocks and shorter time periods but without changing the ranking. We also go beyond binary testing and apply alternative statistics to test the ability to estimate aftershock numbers, which confirm that simple stress metrics perform better than the classic Coulomb failure stress calculations and are also better than the distance-slip probabilistic model. Y1 - 2020 U6 - https://doi.org/10.1029/2020JB019553 SN - 2169-9313 SN - 2169-9356 VL - 125 IS - 9 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Salamat, Mona A1 - Zare, Mehdi A1 - Holschneider, Matthias A1 - Zöller, Gert T1 - Calculation of Confidence Intervals for the Maximum Magnitude of Earthquakes in Different Seismotectonic Zones of Iran JF - Pure and applied geophysics N2 - The problem of estimating the maximum possible earthquake magnitude m(max) has attracted growing attention in recent years. Due to sparse data, the role of uncertainties becomes crucial. In this work, we determine the uncertainties related to the maximum magnitude in terms of confidence intervals. Using an earthquake catalog of Iran, m(max) is estimated for different predefined levels of confidence in six seismotectonic zones. Assuming the doubly truncated Gutenberg-Richter distribution as a statistical model for earthquake magnitudes, confidence intervals for the maximum possible magnitude of earthquakes are calculated in each zone. While the lower limit of the confidence interval is the magnitude of the maximum observed event, the upper limit is calculated from the catalog and the statistical model. For this aim, we use the original catalog which no declustering methods applied on as well as a declustered version of the catalog. Based on the study by Holschneider et al. (Bull Seismol Soc Am 101(4): 1649-1659, 2011), the confidence interval for m(max) is frequently unbounded, especially if high levels of confidence are required. In this case, no information is gained from the data. Therefore, we elaborate for which settings finite confidence levels are obtained. In this work, Iran is divided into six seismotectonic zones, namely Alborz, Azerbaijan, Zagros, Makran, Kopet Dagh, Central Iran. Although calculations of the confidence interval in Central Iran and Zagros seismotectonic zones are relatively acceptable for meaningful levels of confidence, results in Kopet Dagh, Alborz, Azerbaijan and Makran are not that much promising. The results indicate that estimating mmax from an earthquake catalog for reasonable levels of confidence alone is almost impossible. KW - Maximum magnitude of earthquake KW - Level of confidence KW - Confidence interval Y1 - 2016 U6 - https://doi.org/10.1007/s00024-016-1418-5 SN - 0033-4553 SN - 1420-9136 VL - 174 SP - 763 EP - 777 PB - Springer CY - Basel ER - TY - GEN A1 - Lühr, Hermann A1 - Wicht, Johannes A1 - Gilder, Stuart A. A1 - Holschneider, Matthias T1 - Preface T2 - Magnetic Fields in the Solar System: Planets, Moons and Solar Wind Interactions Y1 - 2018 SN - 978-3-319-64292-5 SN - 978-3-319-64291-8 SN - 0067-0057 VL - 448 SP - V EP - VI PB - Springer CY - Dordrecht ER - TY - JOUR A1 - Prokhorov, Boris E. A1 - Förster, Matthias A1 - Lesur, Vincent A1 - Namgaladze, Alexander A. A1 - Holschneider, Matthias A1 - Stolle, Claudia T1 - Modeling of the ionospheric current system and calculating its JF - Magnetic Fields in the Solar System: Planets, Moons and Solar Wind Interactions N2 - The additional magnetic field produced by the ionospheric current system is a part of the Earth’s magnetic field. This current system is a highly variable part of a global electric circuit. The solar wind and interplanetary magnetic field (IMF) interaction with the Earth’s magnetosphere is the external driver for the global electric circuit in the ionosphere. The energy is transferred via the field-aligned currents (FACs) to the Earth’s ionosphere. The interactions between the neutral and charged particles in the ionosphere lead to the so-called thermospheric neutral wind dynamo which represents the second important driver for the global current system. Both processes are components of the magnetosphere–ionosphere–thermosphere (MIT) system, which depends on solar and geomagnetic conditions, and have significant seasonal and UT variations. The modeling of the global dynamic Earth’s ionospheric current system is the first aim of this investigation. For our study, we use the Potsdam version of the Upper Atmosphere Model (UAM-P). The UAM is a first-principle, time-dependent, and fully self-consistent numerical global model. The model includes the thermosphere, ionosphere, plasmasphere, and inner magnetosphere as well as the electrodynamics of the coupled MIT system for the altitudinal range from 80 (60) km up to the 15 Earth radii. The UAM-P differs from the UAM by a new electric field block. For this study, the lower latitudinal and equatorial electrodynamics of the UAM-P model was improved. The calculation of the ionospheric current system’s contribution to the Earth’s magnetic field is the second aim of this study. We present the method, which allows computing the additional magnetic field inside and outside the current layer as generated by the space current density distribution using the Biot-Savart law. Additionally, we perform a comparison of the additional magnetic field calculation using 2D (equivalent currents) and 3D current distribution. Y1 - 2018 SN - 978-3-319-64292-5 SN - 978-3-319-64291-8 U6 - https://doi.org/10.1007/978-3-319-64292-5_10 SN - 0067-0057 SN - 2214-7985 VL - 448 SP - 263 EP - 292 PB - Springer CY - Dordrecht ER - TY - JOUR A1 - Bärenzung, Julien A1 - Holschneider, Matthias A1 - Wicht, Johannes A1 - Sanchez, Sabrina A1 - Lesur, Vincent T1 - Modeling and predicting the short-term evolution of the geomagnetic field JF - Journal of geophysical research : Solid earth N2 - We propose a reduced dynamical system describing the coupled evolution of fluid flow and magnetic field at the top of the Earth's core between the years 1900 and 2014. The flow evolution is modeled with a first-order autoregressive process, while the magnetic field obeys the classical frozen flux equation. An ensemble Kalman filter algorithm serves to constrain the dynamics with the geomagnetic field and its secular variation given by the COV-OBS.x1 model. Using a large ensemble with 40,000 members provides meaningful statistics including reliable error estimates. The model highlights two distinct flow scales. Slowly varying large-scale elements include the already documented eccentric gyre. Localized short-lived structures include distinctly ageostophic features like the high-latitude polar jet on the Northern Hemisphere. Comparisons with independent observations of the length-of-day variations not only validate the flow estimates but also suggest an acceleration of the geostrophic flows over the last century. Hindcasting tests show that our model outperforms simpler predictions bases (linear extrapolation and stationary flow). The predictability limit, of about 2,000 years for the magnetic dipole component, is mostly determined by the random fast varying dynamics of the flow and much less by the geomagnetic data quality or lack of small-scale information. KW - core flow KW - assimilation KW - prediction KW - length of day Y1 - 2018 U6 - https://doi.org/10.1029/2017JB015115 SN - 2169-9313 SN - 2169-9356 VL - 123 IS - 6 SP - 4539 EP - 4560 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Fiedler, Bernhard A1 - Hainzl, Sebastian A1 - Zöller, Gert A1 - Holschneider, Matthias T1 - Detection of Gutenberg-Richter b-Value Changes in Earthquake Time Series JF - Bulletin of the Seismological Society of America N2 - The Gutenberg-Richter relation for earthquake magnitudes is the most famous empirical law in seismology. It states that the frequency of earthquake magnitudes follows an exponential distribution; this has been found to be a robust feature of seismicity above the completeness magnitude, and it is independent of whether global, regional, or local seismicity is analyzed. However, the exponent b of the distribution varies significantly in space and time, which is important for process understanding and seismic hazard assessment; this is particularly true because of the fact that the Gutenberg-Richter b-value acts as a proxy for the stress state and quantifies the ratio of large-to-small earthquakes. In our work, we focus on the automatic detection of statistically significant temporal changes of the b-value in seismicity data. In our approach, we use Bayes factors for model selection and estimate multiple change-points of the frequency-magnitude distribution in time. The method is first applied to synthetic data, showing its capability to detect change-points as function of the size of the sample and the b-value contrast. Finally, we apply this approach to examples of observational data sets for which b-value changes have previously been stated. Our analysis of foreshock and after-shock sequences related to mainshocks, as well as earthquake swarms, shows that only a portion of the b-value changes is statistically significant. Y1 - 2018 U6 - https://doi.org/10.1785/0120180091 SN - 0037-1106 SN - 1943-3573 VL - 108 IS - 5A SP - 2778 EP - 2787 PB - Seismological Society of America CY - Albany ER - TY - JOUR A1 - Diallo, Mamadou Sanou A1 - Kulesh, Michail A1 - Holschneider, Matthias A1 - Kurennaya, Kristina A1 - Scherbaum, Frank T1 - Instantaneous polarization attributes based on an adaptive approximate covariance method JF - Geophysics N2 - We introduce a method for computing instantaneous-polarization attributes from multicomponent signals. This is an improvement on the standard covariance method (SCM) because it does not depend on the window size used to compute the standard covariance matrix. We overcome the window-size problem by deriving an approximate analytical formula for the cross-energy matrix in which we automatically and adaptively determine the time window. The proposed method uses polarization analysis as applied to multicomponent seismic by waveform separation and filtering. Y1 - 2006 U6 - https://doi.org/10.1190/1.2227522 SN - 0016-8033 SN - 1942-2156 (0nline) VL - 71 IS - 5 SP - V99 EP - V104 PB - SEG CY - Tulsa ER - TY - JOUR A1 - Panet, I. A1 - Chambodut, Aude A1 - Diament, M. A1 - Holschneider, Matthias A1 - Jamet, O. T1 - New insights on intraplate volcanism in French Polynesia from wavelet analysis of GRACE, CHAMP, and sea surface data JF - Journal of geophysical research : Solid earth N2 - [ 1] In this paper, we discuss the origin of superswell volcanism on the basis of representation and analysis of recent gravity and magnetic satellite data with wavelets in spherical geometry. We computed a refined gravity field in the south central Pacific based on the GRACE satellite GGM02S global gravity field and the KMS02 altimetric grid, and a magnetic anomaly field based on CHAMP data. The magnetic anomalies are marked by the magnetic lineation of the seafloor spreading and by a strong anomaly in the Tuamotu region, which we interpret as evidence for crustal thickening. We interpret our gravity field through a continuous wavelet analysis that allows to get a first idea of the internal density distribution. We also compute the continuous wavelet analysis of the bathymetric contribution to discriminate between deep and superficial sources. According to the gravity signature of the different chains as revealed by our analysis, various processes are at the origin of the volcanism in French Polynesia. As evidence, we show a large-scale anomaly over the Society Islands that we interpret as the gravity signature of a deeply anchored mantle plume. The gravity signature of the Cook-Austral chain indicates a complex origin which may involve deep processes. Finally, we discuss the particular location of the Marquesas chain as suggesting that the origin of the volcanism may interfere with secondary convection rolls or may be controlled by lithospheric weakness due to the regional stress field, or else related to the presence of the nearby Tuamotu plateau. Y1 - 2006 U6 - https://doi.org/10.1029/2005JB004141 SN - 0148-0227 VL - 111 PB - Union CY - Washington ER - TY - JOUR A1 - Benmehdi, Sabah A1 - Makarava, Natallia A1 - Benhamidouche, N. A1 - Holschneider, Matthias T1 - Bayesian estimation of the self-similarity exponent of the Nile River fluctuation JF - Nonlinear processes in geophysics N2 - The aim of this paper is to estimate the Hurst parameter of Fractional Gaussian Noise (FGN) using Bayesian inference. We propose an estimation technique that takes into account the full correlation structure of this process. Instead of using the integrated time series and then applying an estimator for its Hurst exponent, we propose to use the noise signal directly. As an application we analyze the time series of the Nile River, where we find a posterior distribution which is compatible with previous findings. In addition, our technique provides natural error bars for the Hurst exponent. Y1 - 2011 U6 - https://doi.org/10.5194/npg-18-441-2011 SN - 1023-5809 VL - 18 IS - 3 SP - 441 EP - 446 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Panet, Isabelle A1 - Kuroishi, Yuki A1 - Holschneider, Matthias T1 - Wavelet modelling of the gravity field by domain decomposition methods: an example over Japan JF - Geophysical journal international N2 - With the advent of satellite gravity, large gravity data sets of unprecedented quality at low and medium resolution become available. For local, high resolution field modelling, they need to be combined with the surface gravity data. Such models are then used for various applications, from the study of the Earth interior to the determination of oceanic currents. Here we show how to realize such a combination in a flexible way using spherical wavelets and applying a domain decomposition approach. This iterative method, based on the Schwarz algorithms, allows to split a large problem into smaller ones, and avoids the calculation of the entire normal system, which may be huge if high resolution is sought over wide areas. A subdomain is defined as the harmonic space spanned by a subset of the wavelet family. Based on the localization properties of the wavelets in space and frequency, we define hierarchical subdomains of wavelets at different scales. On each scale, blocks of subdomains are defined by using a tailored spatial splitting of the area. The data weighting and regularization are iteratively adjusted for the subdomains, which allows to handle heterogeneity in the data quality or the gravity variations. Different levels of approximations of the subdomains normals are also introduced, corresponding to building local averages of the data at different resolution levels. We first provide the theoretical background on domain decomposition methods. Then, we validate the method with synthetic data, considering two kinds of noise: white noise and coloured noise. We then apply the method to data over Japan, where we combine a satellite-based geopotential model, EIGEN-GL04S, and a local gravity model from a combination of land and marine gravity data and an altimetry-derived marine gravity model. A hybrid spherical harmonics/wavelet model of the geoid is obtained at about 15 km resolution and a corrector grid for the surface model is derived. KW - Wavelet transform KW - Satellite geodesy KW - Geopotential theory Y1 - 2011 U6 - https://doi.org/10.1111/j.1365-246X.2010.04840.x SN - 0956-540X VL - 184 IS - 1 SP - 203 EP - 219 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Shebalin, Peter A1 - Narteau, Clement A1 - Holschneider, Matthias A1 - Schorlemmer, Danijel T1 - Short-Term earthquake forecasting using early aftershock statistics JF - Bulletin of the Seismological Society of America N2 - We present an alarm-based earthquake forecast model that uses the early aftershock statistics (EAST). This model is based on the hypothesis that the time delay before the onset of the power-law aftershock decay rate decreases as the level of stress and the seismogenic potential increase. Here, we estimate this time delay from < t(g)>, the time constant of the Omori-Utsu law. To isolate space-time regions with a relative high level of stress, the single local variable of our forecast model is the E-a value, the ratio between the long-term and short-term estimations of < t(g)>. When and where the E-a value exceeds a given threshold (i.e., the c value is abnormally small), an alarm is issued, and an earthquake is expected to occur during the next time step. Retrospective tests show that the EAST model has better predictive power than a stationary reference model based on smoothed extrapolation of past seismicity. The official prospective test for California started on 1 July 2009 in the testing center of the Collaboratory for the Study of Earthquake Predictability (CSEP). During the first nine months, 44 M >= 4 earthquakes occurred in the testing area. For this time period, the EAST model has better predictive power than the reference model at a 1% level of significance. Because the EAST model has also a better predictive power than several time-varying clustering models tested in CSEP at a 1% level of significance, we suggest that our successful prospective results are not due only to the space-time clustering of aftershocks. Y1 - 2011 U6 - https://doi.org/10.1785/0120100119 SN - 0037-1106 VL - 101 IS - 1 SP - 297 EP - 312 PB - Seismological Society of America CY - El Cerrito ER - TY - JOUR A1 - Zoller, Gert A1 - Holschneider, Matthias A1 - Ben-Zion, Yehuda T1 - Quasi-static and quasi-dynamic modeling of earthquake failure at intermediate scales N2 - We present a model for earthquake failure at intermediate scales (space: 100 m-100 km, time: 100 m/nu(shear Y1 - 2004 ER - TY - JOUR A1 - Holschneider, Matthias A1 - Teramo, A. A1 - Bottari, A. A1 - Termini, D. T1 - On the relevance of the spatial distribution of events for seismic hazard evaluation N2 - Seismic hazard evaluation is proposed by a methodological approach that allows the study of the influence of different modelling assumptions relative to the spatial and temporal distribution of earthquakes on the maximum values of expected intensities. In particular, we show that the estimated hazard at a fixed point is very sensitive to the assumed spatial distribution of epicentres and their estimators. As we will see, the usual approach, based on uniformly distributing the epicentres inside each seismogenic zone is likely to be biased towards lower expected intensity values. This will be made more precise later. Recall that the term "bias" means, that the expectation of the estimated quantity ( taken as a random variable on the space of statistics) is different from the expectation of the quantity itself. Instead, our approach, based on an estimator that takes into account the observed clustering of events is essentially unbiased, as shown by a Monte-Carlo simulation, and is configured on a 11011-isotropic macroseismic attenuation model which is independently estimated for each zone Y1 - 2004 SN - 0921-030X ER - TY - JOUR A1 - Narteau, C. A1 - Shebalin, P. A1 - Holschneider, Matthias T1 - Onset of power law aftershock decay rates in southern California N2 - Aftershocks rates seem to follow a power law decay, but the question of the aftershock frequency immediately after an earthquake remains open. We estimate an average aftershock decay rate within one day in southern California by stacking in time different sequences triggered by main shocks ranging in magnitude from 2.5 to 4.5. Then we estimate the time delay before the onset of the power law aftershock decay rate. For the last 20 years, we observe that this time delay suddenly increase after large earthquakes, and slowly decreases at a constant rate during periods of low seismicity. In a band-limited power law model such variations can be explained by different patterns of stress distribution at different stages of the seismic cycle. We conclude that, on regional length scales, the brittle upper crust exhibits a collective behavior reflecting to some extent the proximity of a threshold of fracturing Y1 - 2005 SN - 0094-8276 ER - TY - JOUR A1 - Chambodut, Aude A1 - Panet, I. A1 - Mandea, Mioara A1 - Diament, M. A1 - Holschneider, Matthias A1 - Jamet, O. T1 - Wavelet frames : an alternative to spherical harmonic representation of potential fields N2 - Potential fields are classically represented on the sphere using spherical harmonics. However, this decomposition leads to numerical difficulties when data to be modelled are irregularly distributed or cover a regional zone. To overcome this drawback, we develop a new representation of the magnetic and the gravity fields based on wavelet frames. In this paper, we first describe how to build wavelet frames on the sphere. The chosen frames are based on the Poisson multipole wavelets, which are of special interest for geophysical modelling, since their scaling parameter is linked to the multipole depth (Holschneider et al.). The implementation of wavelet frames results from a discretization of the continuous wavelet transform in space and scale. We also build different frames using two kinds of spherical meshes and various scale sequences. We then validate the mathematical method through simple fits of scalar functions on the sphere, named 'scalar models'. Moreover, we propose magnetic and gravity models, referred to as 'vectorial models', taking into account geophysical constraints. We then discuss the representation of the Earth's magnetic and gravity fields from data regularly or irregularly distributed. Comparisons of the obtained wavelet models with the initial spherical harmonic models point out the advantages of wavelet modelling when the used magnetic or gravity data are sparsely distributed or cover just a very local zone Y1 - 2005 SN - 0956-540X ER - TY - JOUR A1 - Diallo, Mamadou Sanou A1 - Kulesh, Michail A1 - Holschneider, Matthias A1 - Scherbaum, Frank A1 - Adler, Frank T1 - Characterization of polarization attributes of seismic waves using continuous wavelet transforms N2 - Complex-trace analysis is the method of choice for analyzing polarized data. Because particle motion can be represented by instantaneous attributes that show distinct features for waves of different polarization characteristics, it can be used to separate and characterize these waves. Traditional methods of complex-trace analysis only give the instantaneous attributes as a function of time or frequency. However. for transient wave types or seismic events that overlap in time, an estimate of the polarization parameters requires analysis of the time-frequency dependence of these attributes. We propose a method to map instantaneous polarization attributes of seismic signals in the wavelet domain and explicitly relate these attributes with the wavelet-transform coefficients of the analyzed signal. We compare our method with traditional complex-trace analysis using numerical examples. An advantage of our method is its possibility of performing the complete wave-mode separation/ filtering process in the wavelet domain and its ability to provide the frequency dependence of ellipticity, which contains important information on the subsurface structure. Furthermore, using 2-C synthetic and real seismic shot gathers, we show how to use the method to separate different wave types and identify zones of interfering wave modes Y1 - 2006 UR - http://geophysics.geoscienceworld.org/ U6 - https://doi.org/10.1190/1.2194511 SN - 0016-8033 ER - TY - JOUR A1 - Kulesh, Michail A1 - Holschneider, Matthias A1 - Diallo, Mamadou Sanou A1 - Xie, Q. A1 - Scherbaum, Frank T1 - Modeling of wave dispersion using continuous wavelet transforms N2 - In the estimate of dispersion with the help of wavelet analysis considerable emphasis has been put on the extraction of the group velocity using the modulus of the wavelet transform. In this paper we give an asymptotic expression of the full propagator in wavelet space that comprises the phase velocity as well. This operator establishes a relationship between the observed signals at two different stations during wave propagation in a dispersive and attenuating medium. Numerical and experimental examples are presented to show that the method accurately models seismic wave dispersion and attenuation Y1 - 2005 SN - 0033-4553 ER -