TY - JOUR A1 - Baerenzung, Julien A1 - Holschneider, Matthias A1 - Lesur, Vincent T1 - Bayesian inversion for the filtered flow at the Earth's core-mantle boundary JF - Journal of geophysical research : Solid earth N2 - The inverse problem of determining the flow at the Earth's core-mantle boundary according to an outer core magnetic field and secular variation model has been investigated through a Bayesian formalism. To circumvent the issue arising from the truncated nature of the available fields, we combined two modeling methods. In the first step, we applied a filter on the magnetic field to isolate its large scales by reducing the energy contained in its small scales, we then derived the dynamical equation, referred as filtered frozen flux equation, describing the spatiotemporal evolution of the filtered part of the field. In the second step, we proposed a statistical parametrization of the filtered magnetic field in order to account for both its remaining unresolved scales and its large-scale uncertainties. These two modeling techniques were then included in the Bayesian formulation of the inverse problem. To explore the complex posterior distribution of the velocity field resulting from this development, we numerically implemented an algorithm based on Markov chain Monte Carlo methods. After evaluating our approach on synthetic data and comparing it to previously introduced methods, we applied it to a magnetic field model derived from satellite data for the single epoch 2005.0. We could confirm the existence of specific features already observed in previous studies. In particular, we retrieved the planetary scale eccentric gyre characteristic of flow evaluated under the compressible quasi-geostrophy assumption although this hypothesis was not considered in our study. In addition, through the sampling of the velocity field posterior distribution, we could evaluate the reliability, at any spatial location and at any scale, of the flow we calculated. The flow uncertainties we determined are nevertheless conditioned by the choice of the prior constraints we applied to the velocity field. Y1 - 2014 U6 - https://doi.org/10.1002/2013JB010358 SN - 2169-9313 SN - 2169-9356 VL - 119 IS - 4 SP - 2695 EP - 2720 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Baerenzung, Julien A1 - Holschneider, Matthias A1 - Wicht, Johannes A1 - Lesur, Vincent A1 - Sanchez, Sabrina T1 - The Kalmag model as a candidate for IGRF-13 JF - Earth, planets and space N2 - We present a new model of the geomagnetic field spanning the last 20 years and called Kalmag. Deriving from the assimilation of CHAMP and Swarm vector field measurements, it separates the different contributions to the observable field through parameterized prior covariance matrices. To make the inverse problem numerically feasible, it has been sequentialized in time through the combination of a Kalman filter and a smoothing algorithm. The model provides reliable estimates of past, present and future mean fields and associated uncertainties. The version presented here is an update of our IGRF candidates; the amount of assimilated data has been doubled and the considered time window has been extended from [2000.5, 2019.74] to [2000.5, 2020.33]. KW - Geomagnetic field KW - Secular variation KW - Assimilation KW - Kalman filter KW - Machine learning Y1 - 2020 U6 - https://doi.org/10.1186/s40623-020-01295-y SN - 1880-5981 VL - 72 IS - 1 PB - Springer CY - New York ER - TY - JOUR A1 - Bailey, Iain W. A1 - Ben-Zion, Yehuda A1 - Becker, Thorsten W. A1 - Holschneider, Matthias T1 - Quantifying focal mechanism heterogeneity for fault zones in central and southern California N2 - P>We present a statistical analysis of focal mechanism orientations for nine California fault zones with the goal of quantifying variations of fault zone heterogeneity at seismogenic depths. The focal mechanism data are generated from first motion polarities for earthquakes in the time period 1983-2004, magnitude range 0-5, and depth range 0-15 km. Only mechanisms with good quality solutions are used. We define fault zones using 20 km wide rectangles and use summations of normalized potency tensors to describe the distribution of double-couple orientations for each fault zone. Focal mechanism heterogeneity is quantified using two measures computed from the tensors that relate to the scatter in orientations and rotational asymmetry or skewness of the distribution. We illustrate the use of these quantities by showing relative differences in the focal mechanism heterogeneity characteristics for different fault zones. These differences are shown to relate to properties of the fault zone surface traces such that increased scatter correlates with fault trace complexity and rotational asymmetry correlates with the dominant fault trace azimuth. These correlations indicate a link between the long-term evolution of a fault zone over many earthquake cycles and its seismic behaviour over a 20 yr time period. Analysis of the partitioning of San Jacinto fault zone focal mechanisms into different faulting styles further indicates that heterogeneity is dominantly controlled by structural properties of the fault zone, rather than time or magnitude related properties of the seismicity. Y1 - 2010 UR - http://onlinelibrary.wiley.com/journal/10.1111/(ISSN)1365-246X U6 - https://doi.org/10.1111/j.1365-246X.2010.04745.x SN - 0956-540X 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 - 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 - Bettenbühl, Mario A1 - Rusconi, Marco A1 - Engbert, Ralf A1 - Holschneider, Matthias T1 - Bayesian selection of Markov Models for symbol sequences application to microsaccadic eye movements JF - PLoS one N2 - Complex biological dynamics often generate sequences of discrete events which can be described as a Markov process. The order of the underlying Markovian stochastic process is fundamental for characterizing statistical dependencies within sequences. As an example for this class of biological systems, we investigate the Markov order of sequences of microsaccadic eye movements from human observers. We calculate the integrated likelihood of a given sequence for various orders of the Markov process and use this in a Bayesian framework for statistical inference on the Markov order. Our analysis shows that data from most participants are best explained by a first-order Markov process. This is compatible with recent findings of a statistical coupling of subsequent microsaccade orientations. Our method might prove to be useful for a broad class of biological systems. Y1 - 2012 U6 - https://doi.org/10.1371/journal.pone.0043388 SN - 1932-6203 VL - 7 IS - 9 PB - PLoS CY - San Fransisco ER - TY - JOUR A1 - Boehm, Thorsten A1 - Holschneider, Matthias A1 - Lignieres, Frederic A1 - Petit, Pascal A1 - Rainer, Monica A1 - Paletou, Francois A1 - Wade, Gregg A1 - Alecian, Evelyne A1 - Carfantan, Herve A1 - Blazere, Aurore A1 - Mirouh, Giovanni M. T1 - Discovery of starspots on Vega First spectroscopic detection of surface structures on a normal A-type star JF - Astronomy and astrophysics : an international weekly journal N2 - Context. The theoretically studied impact of rapid rotation on stellar evolution needs to be compared with these results of high-resolution spectroscopy-velocimetry observations. Early-type stars present a perfect laboratory for these studies. The prototype A0 star Vega has been extensively monitored in recent years in spectropolarimetry. A weak surface magnetic field was detected, implying that there might be a (still undetected) structured surface. First indications of the presence of small amplitude stellar radial velocity variations have been reported recently, but the confirmation and in-depth study with the highly stabilized spectrograph SOPHIE/OHP was required. Aims. The goal of this article is to present a thorough analysis of the line profile variations and associated estimators in the early-type standard star Vega (A0) in order to reveal potential activity tracers, exoplanet companions, and stellar oscillations. Methods. Vega was monitored in quasi-continuous high-resolution echelle spectroscopy with the highly stabilized velocimeter SOPHIE/OHP. A total of 2588 high signal-to-noise spectra was obtained during 34.7 h on five nights (2 to 6 of August 2012) in high-resolution mode at R = 75 000 and covering the visible domain from 3895 6270 angstrom. For each reduced spectrum, least square deconvolved equivalent photospheric profiles were calculated with a T-eff = 9500 and log g = 4.0 spectral line mask. Several methods were applied to study the dynamic behaviour of the profile variations (evolution of radial velocity, bisectors, vspan, 2D profiles, amongst others). Results. We present the discovery of a spotted stellar surface on an A-type standard star (Vega) with very faint spot amplitudes Delta F/Fc similar to 5 x 10(-4). A rotational modulation of spectral lines with a period of rotation P = 0.68 d has clearly been exhibited, unambiguously confirming the results of previous spectropolarimetric studies. Most of these brightness inhomogeneities seem to be located in lower equatorial latitudes. Either a very thin convective layer can be responsible for magnetic field generation at small amplitudes, or a new mechanism has to be invoked to explain the existence of activity tracing starspots. At this stage it is difficult to disentangle a rotational from a stellar pulsational origin for the existing higher frequency periodic variations. Conclusions. This first strong evidence that standard A-type stars can show surface structures opens a new field of research and ask about a potential link with the recently discovered weak magnetic field discoveries in this category of stars. KW - starspots KW - stars: early-type KW - stars: rotation KW - stars: oscillations KW - stars: individual: Vega KW - asteroseismology Y1 - 2015 U6 - https://doi.org/10.1051/0004-6361/201425425 SN - 0004-6361 SN - 1432-0746 VL - 577 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Bärenzung, Julien A1 - Holschneider, Matthias A1 - Lesur, Vincent T1 - constraints JF - Journal of geophysical research : Solid earth N2 - Prior information in ill-posed inverse problem is of critical importance because it is conditioning the posterior solution and its associated variability. The problem of determining the flow evolving at the Earth's core-mantle boundary through magnetic field models derived from satellite or observatory data is no exception to the rule. This study aims to estimate what information can be extracted on the velocity field at the core-mantle boundary, when the frozen flux equation is inverted under very weakly informative, but realistic, prior constraints. Instead of imposing a converging spectrum to the flow, we simply assume that its poloidal and toroidal energy spectra are characterized by power laws. The parameters of the spectra, namely, their magnitudes, and slopes are unknown. The connection between the velocity field, its spectra parameters, and the magnetic field model is established through the Bayesian formulation of the problem. Working in two steps, we determined the time-averaged spectra of the flow within the 2001–2009.5 period, as well as the flow itself and its associated uncertainties in 2005.0. According to the spectra we obtained, we can conclude that the large-scale approximation of the velocity field is not an appropriate assumption within the time window we considered. For the flow itself, we show that although it is dominated by its equatorial symmetric component, it is very unlikely to be perfectly symmetric. We also demonstrate that its geostrophic state is questioned in different locations of the outer core. Y1 - 2016 U6 - https://doi.org/10.1002/2015JB012464 SN - 2169-9313 SN - 2169-9356 VL - 121 SP - 1343 EP - 1364 PB - American Geophysical Union CY - Washington 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 - 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 - Cotronei, Mariantonia A1 - Di Salvo, Rosa A1 - Holschneider, Matthias A1 - Puccio, Luigia T1 - Interpolation in reproducing kernel Hilbert spaces based on random subdivision schemes JF - Journal of computational and applied mathematics N2 - In this paper we present a Bayesian framework for interpolating data in a reproducing kernel Hilbert space associated with a random subdivision scheme, where not only approximations of the values of a function at some missing points can be obtained, but also uncertainty estimates for such predicted values. This random scheme generalizes the usual subdivision by taking into account, at each level, some uncertainty given in terms of suitably scaled noise sequences of i.i.d. Gaussian random variables with zero mean and given variance, and generating, in the limit, a Gaussian process whose correlation structure is characterized and used for computing realizations of the conditional posterior distribution. The hierarchical nature of the procedure may be exploited to reduce the computational cost compared to standard techniques in the case where many prediction points need to be considered. KW - Subdivision schemes KW - Interpolation KW - Simulation of Gaussian processes KW - Bayesian inversion Y1 - 2016 U6 - https://doi.org/10.1016/j.cam.2016.08.002 SN - 0377-0427 SN - 1879-1778 VL - 311 SP - 342 EP - 353 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Cotronei, Mariantonia A1 - Holschneider, Matthias T1 - Partial parameterization of orthogonal wavelet matrix filters JF - Journal of computational and applied mathematics N2 - In this paper we propose a procedure which allows the construction of a large family of FIR d x d matrix wavelet filters by exploiting the one-to-one correspondence between QMF systems and orthogonal operators which commute with the shifts by two. A characterization of the class of filters of full rank type that can be obtained with such procedure is given. In particular, we restrict our attention to a special construction based on the representation of SO(2d) in terms of the elements of its Lie algebra. Explicit expressions for the filters in the case d = 2 are given, as a result of a local analysis of the parameterization obtained from perturbing the Haar system. KW - Full rank matrix filters KW - Multichannel wavelets KW - Quadrature mirror filters KW - Vector subdivision schemes Y1 - 2013 U6 - https://doi.org/10.1016/j.cam.2012.11.016 SN - 0377-0427 VL - 243 IS - 4 SP - 113 EP - 125 PB - Elsevier CY - Amsterdam 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 - Diallo, Mamadou Sanou A1 - Kulesh, Michail A1 - Holschneider, Matthias A1 - Scherbaum, Frank T1 - Instantaneous polarization attributes in the time-frequency domain and wavefield separation N2 - We introduce a method of wavefield separation from multicomponent data sets based on the use of the continuous wavelet transform. Our method is a further generalization of the approach proposed by Morozov and Smithson, in that by using the continuous wavelet transform, we can achieve a better separation of wave types by designing the filter in the time-frequency domain. Furthermore, using the instantaneous polarization attributes defined in the wavelet domain, we show how to construct filters tailored to separate different wave types (elliptically or linearly polarized), followed by an inverse wavelet transform to obtain the desired wave type in the time domain. Using synthetic and experimental data, we show how the present method can be used for wavefield separation Y1 - 2005 SN - 0016-8025 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 - 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 - Fiedler, Bernhard A1 - Zöller, Gert A1 - Holschneider, Matthias A1 - Hainzl, Sebastian T1 - Multiple Change-Point Detection in Spatiotemporal Seismicity Data JF - Bulletin of the Seismological Society of America N2 - Earthquake rates are driven by tectonic stress buildup, earthquake-induced stress changes, and transient aseismic processes. Although the origin of the first two sources is known, transient aseismic processes are more difficult to detect. However, the knowledge of the associated changes of the earthquake activity is of great interest, because it might help identify natural aseismic deformation patterns such as slow-slip events, as well as the occurrence of induced seismicity related to human activities. For this goal, we develop a Bayesian approach to identify change-points in seismicity data automatically. Using the Bayes factor, we select a suitable model, estimate possible change-points, and we additionally use a likelihood ratio test to calculate the significance of the change of the intensity. The approach is extended to spatiotemporal data to detect the area in which the changes occur. The method is first applied to synthetic data showing its capability to detect real change-points. Finally, we apply this approach to observational data from Oklahoma and observe statistical significant changes of seismicity in space and time. Y1 - 2018 U6 - https://doi.org/10.1785/0120170236 SN - 0037-1106 SN - 1943-3573 VL - 108 IS - 3A SP - 1147 EP - 1159 PB - Seismological Society of America CY - Albany ER - TY - JOUR A1 - Gaci, Said A1 - Zaourar, Naima A1 - Briqueu, Louis A1 - Holschneider, Matthias T1 - Regularity analysis applied to sonic logs data a case study from KTB borehole site JF - Arabian journal of geosciences N2 - Borehole logs provide in situ information about the fluctuations of petrophysical properties with depth and thus allow the characterization of the crustal heterogeneities. A detailed investigation of these measurements may lead to extract features of the geological media. In this study, we suggest a regularity analysis based on the continuous wavelet transform to examine sonic logs data. The description of the local behavior of the logs at each depth is carried out using the local Hurst exponent estimated by two (02) approaches: the local wavelet approach and the average-local wavelet approach. Firstly, a synthetic log, generated using the random midpoints displacement algorithm, is processed by the regularity analysis. The obtained Hurst curves allowed the discernment of the different layers composing the simulated geological model. Next, this analysis is extended to real sonic logs data recorded at the Kontinentales Tiefbohrprogramm (KTB) pilot borehole (Continental Deep Drilling Program, Germany). The results show a significant correlation between the estimated Hurst exponents and the lithological discontinuities crossed by the well. Hence, the Hurst exponent can be used as a tool to characterize underground heterogeneities. KW - Regularity analysis KW - Wavelet transform KW - Well log KW - Fractal KW - Lithology Y1 - 2011 U6 - https://doi.org/10.1007/s12517-010-0129-y SN - 1866-7511 VL - 4 IS - 1-2 SP - 221 EP - 227 PB - Springer CY - Heidelberg ER - TY - JOUR A1 - Gaci, Said A1 - Zaourar, Naima A1 - Hamoudi, Mehdi A1 - Holschneider, Matthias T1 - Local regularity analysis of strata heterogeneities from sonic logs N2 - Borehole logs provide geological information about the rocks crossed by the wells. Several properties of rocks can be interpreted in terms of lithology, type and quantity of the fluid filling the pores and fractures. Here, the logs are assumed to be nonhomogeneous Brownian motions (nhBms) which are generalized fractional Brownian motions (fBms) indexed by depth-dependent Hurst parameters H(z). Three techniques, the local wavelet approach (LWA), the average-local wavelet approach (ALWA), and Peltier Algorithm (PA), are suggested to estimate the Hurst functions (or the regularity profiles) from the logs. First, two synthetic sonic logs with different parameters, shaped by the successive random additions (SRA) algorithm, are used to demonstrate the potential of the proposed methods. The obtained Hurst functions are close to the theoretical Hurst functions. Besides, the transitions between the modeled layers are marked by Hurst values discontinuities. It is also shown that PA leads to the best Hurst value estimations. Second, we investigate the multifractional property of sonic logs data recorded at two scientific deep boreholes: the pilot hole VB and the ultra deep main hole HB, drilled for the German Continental Deep Drilling Program (KTB). All the regularity profiles independently obtained for the logs provide a clear correlation with lithology, and from each regularity profile, we derive a similar segmentation in terms of lithological units. The lithological discontinuities (strata' bounds and faults contacts) are located at the local extrema of the Hurst functions. Moreover, the regularity profiles are compared with the KTB estimated porosity logs, showing a significant relation between the local extrema of the Hurst functions and the fluid-filled fractures. The Hurst function may then constitute a tool to characterize underground heterogeneities. Y1 - 2010 UR - http://www.copernicus.org/EGU/npg/npg.htm U6 - https://doi.org/10.5194/npg-17-455-2010 SN - 1023-5809 ER - TY - JOUR A1 - Hayn, Michael A1 - Holschneider, Matthias T1 - Directional spherical multipole wavelets N2 - We construct a family of admissible analysis reconstruction pairs of wavelet families on the sphere. The construction is an extension of the isotropic Poisson wavelets. Similar to those, the directional wavelets allow a finite expansion in terms of off-center multipoles. Unlike the isotropic case, the directional wavelets are not a tight frame. However, at small scales, they almost behave like a tight frame. We give an explicit formula for the pseudodifferential operator given by the combination analysis-synthesis with respect to these wavelets. The Euclidean limit is shown to exist and an explicit formula is given. This allows us to quantify the asymptotic angular resolution of the wavelets. Y1 - 2009 UR - http://jmp.aip.org/ U6 - https://doi.org/10.1063/1.3177198 SN - 0022-2488 ER - TY - JOUR A1 - Hayn, Michael A1 - Panet, I. A1 - Diament, M. A1 - Holschneider, Matthias A1 - Mandea, Mioara A1 - Davaille, A. T1 - Wavelet-based directional analysis of the gravity field evidence for large-scale undulations JF - Geophysical journal international N2 - In the eighties, the analysis of satellite altimetry data leads to the major discovery of gravity lineations in the oceans, with wavelengths between 200 and 1400 km. While the existence of the 200 km scale undulations is widely accepted, undulations at scales larger than 400 km are still a matter of debate. In this paper, we revisit the topic of the large-scale geoid undulations over the oceans in the light of the satellite gravity data provided by the GRACE mission, considerably more precise than the altimetry data at wavelengths larger than 400 km. First, we develop a dedicated method of directional Poisson wavelet analysis on the sphere with significance testing, in order to detect and characterize directional structures in geophysical data on the sphere at different spatial scales. This method is particularly well suited for potential field analysis. We validate it on a series of synthetic tests, and then apply it to analyze recent gravity models, as well as a bathymetry data set independent from gravity. Our analysis confirms the existence of gravity undulations at large scale in the oceans, with characteristic scales between 600 and 2000 km. Their direction correlates well with present-day plate motion over the Pacific ocean, where they are particularly clear, and associated with a conjugate direction at 1500 km scale. A major finding is that the 2000 km scale geoid undulations dominate and had never been so clearly observed previously. This is due to the great precision of GRACE data at those wavelengths. Given the large scale of these undulations, they are most likely related to mantle processes. Taking into account observations and models from other geophysical information, as seismological tomography, convection and geochemical models and electrical conductivity in the mantle, we conceive that all these inputs indicate a directional fabric of the mantle flows at depth, reflecting how the history of subduction influences the organization of lower mantle upwellings. KW - Wavelet transform KW - Satellite geodesy KW - Gravity anomalies and Earth structure KW - Pacific Ocean Y1 - 2012 U6 - https://doi.org/10.1111/j.1365-246X.2012.05455.x SN - 0956-540X SN - 1365-246X VL - 189 IS - 3 SP - 1430 EP - 1456 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Holschneider, Matthias A1 - Diallo, Mamadou Sanou A1 - Kulesh, Michail A1 - Ohrnberger, Matthias A1 - Luck, E. A1 - Scherbaum, Frank T1 - Characterization of dispersive surface waves using continuous wavelet transforms N2 - In this paper, we propose a method of surface waves characterization based on the deformation of the wavelet transform of the analysed signal. An estimate of the phase velocity (the group velocity) and the attenuation coefficient is carried out using a model-based approach to determine the propagation operator in the wavelet domain, which depends nonlinearly on a set of unknown parameters. These parameters explicitly define the phase velocity, the group velocity and the attenuation. Under the assumption that the difference between waveforms observed at a couple of stations is solely due to the dispersion characteristics and the intrinsic attenuation of the medium, we then seek to find the set of unknown parameters of this model. Finding the model parameters turns out to be that of an optimization problem, which is solved through the minimization of an appropriately defined cost function. We show that, unlike time-frequency methods that exploit only the square modulus of the transform, we can achieve a complete characterization of surface waves in a dispersive and attenuating medium. Using both synthetic examples and experimental data, we also show that it is in principle possible to separate different modes in both the time domain and the frequency domain Y1 - 2005 SN - 0956-540X ER - TY - JOUR A1 - Holschneider, Matthias A1 - Lesur, Vincent A1 - Mauerberger, Stefan A1 - Baerenzung, Julien T1 - Correlation-based modeling and separation of geomagnetic field components JF - Journal of geophysical research : Solid earth N2 - We introduce a technique for the modeling and separation of geomagnetic field components that is based on an analysis of their correlation structures alone. The inversion is based on a Bayesian formulation, which allows the computation of uncertainties. The technique allows the incorporation of complex measurement geometries like observatory data in a simple way. We show how our technique is linked to other well-known inversion techniques. A case study based on observational data is given. Y1 - 2016 U6 - https://doi.org/10.1002/2015JB012629 SN - 2169-9313 SN - 2169-9356 VL - 121 SP - 3142 EP - 3160 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Holschneider, Matthias A1 - Narteau, C. A1 - Shebalin, P. A1 - Peng, Z. A1 - Schorlemmer, Danijel T1 - Bayesian analysis of the modified Omori law JF - Journal of geophysical research : Solid earth N2 - In order to examine variations in aftershock decay rate, we propose a Bayesian framework to estimate the {K, c, p}-values of the modified Omori law (MOL), lambda(t) = K(c + t)(-p). The Bayesian setting allows not only to produce a point estimator of these three parameters but also to assess their uncertainties and posterior dependencies with respect to the observed aftershock sequences. Using a new parametrization of the MOL, we identify the trade-off between the c and p-value estimates and discuss its dependence on the number of aftershocks. Then, we analyze the influence of the catalog completeness interval [t(start), t(stop)] on the various estimates. To test this Bayesian approach on natural aftershock sequences, we use two independent and non-overlapping aftershock catalogs of the same earthquakes in Japan. Taking into account the posterior uncertainties, we show that both the handpicked (short times) and the instrumental (long times) catalogs predict the same ranges of parameter values. We therefore conclude that the same MOL may be valid over short and long times. Y1 - 2012 U6 - https://doi.org/10.1029/2011JB009054 SN - 2169-9313 SN - 2169-9356 VL - 117 IS - 6089 PB - American Geophysical Union CY - Washington 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 - Holschneider, Matthias A1 - Teschke, Gerd T1 - Existence and computation of optimally localized coherent states JF - Journal of mathematical physics N2 - This paper is concerned with localization properties of coherent states. Instead of classical uncertainty relations we consider "generalized" localization quantities. This is done by introducing measures on the reproducing kernel. In this context we may prove the existence of optimally localized states. Moreover, we provide a numerical scheme for deriving them. Y1 - 2006 U6 - https://doi.org/10.1063/1.2375031 SN - 0022-2488 SN - 1089-7658 VL - 47 IS - 3 SP - 211 EP - 214 PB - Elsevier CY - Melville ER - TY - JOUR A1 - Holschneider, Matthias A1 - Zöller, Gert A1 - Clements, R. A1 - Schorlemmer, Danijel T1 - Can we test for the maximum possible earthquake magnitude? JF - Journal of geophysical research : Solid earth Y1 - 2014 U6 - https://doi.org/10.1002/2013JB010319 SN - 2169-9313 SN - 2169-9356 VL - 119 IS - 3 SP - 2019 EP - 2028 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Holschneider, Matthias A1 - Zöller, Gert A1 - Hainzl, Sebastian T1 - Estimation of the maximum possible magnitude in the framework of a doubly truncated Gutenberg-Richter Model JF - Bulletin of the Seismological Society of America N2 - We discuss to what extent a given earthquake catalog and the assumption of a doubly truncated Gutenberg-Richter distribution for the earthquake magnitudes allow for the calculation of confidence intervals for the maximum possible magnitude M. We show that, without further assumptions such as the existence of an upper bound of M, only very limited information may be obtained. In a frequentist formulation, for each confidence level alpha the confidence interval diverges with finite probability. In a Bayesian formulation, the posterior distribution of the upper magnitude is not normalizable. We conclude that the common approach to derive confidence intervals from the variance of a point estimator fails. Technically, this problem can be overcome by introducing an upper bound (M) over tilde for the maximum magnitude. Then the Bayesian posterior distribution can be normalized, and its variance decreases with the number of observed events. However, because the posterior depends significantly on the choice of the unknown value of (M) over tilde, the resulting confidence intervals are essentially meaningless. The use of an informative prior distribution accounting for pre-knowledge of M is also of little use, because the prior is only modified in the case of the occurrence of an extreme event. Our results suggest that the maximum possible magnitude M should be better replaced by M(T), the maximum expected magnitude in a given time interval T, for which the calculation of exact confidence intervals becomes straightforward. From a physical point of view, numerical models of the earthquake process adjusted to specific fault regions may be a powerful alternative to overcome the shortcomings of purely statistical inference. Y1 - 2011 U6 - https://doi.org/10.1785/0120100289 SN - 0037-1106 VL - 101 IS - 4 SP - 1649 EP - 1659 PB - Seismological Society of America CY - El Cerrito ER - TY - JOUR A1 - Kulesh, Michail A. A1 - Diallo, Mamadou Sanou A1 - Holschneider, Matthias T1 - Wavelet analysis of ellipticity, dispersion, and dissipation properties of Rayleigh waves N2 - This paper is devoted to the digital processing of multicomponent seismograms using wavelet analysis. The goal of this processing is to identify Rayleigh surface elastic waves and determine their properties. A new method for calculating the ellipticity parameters of a wave in the form of a time-frequency spectrum is proposed, which offers wide possibilities for filtering seismic signals in order to suppress or extract the Rayleigh components. A model of dispersion and dissipation of elliptic waves written in terms of wavelet spectra of complex (two-component) signals is also proposed. The model is used to formulate a nonlinear minimization problem that allows for a high-accuracy calculation of the group and phase velocities and the attenuation factor for a propagating elliptic Rayleigh wave. All methods considered in the paper are illustrated with the use of test signals. (c) 2005 Pleiades Publishing, Inc Y1 - 2005 SN - 1063-7710 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 - TY - JOUR A1 - Kurtenbach, E. A1 - Eicker, A. A1 - Mayer-Guerr, T. A1 - Holschneider, Matthias A1 - Hayn, M. A1 - Fuhrmann, M. A1 - Kusche, J. T1 - Improved daily GRACE gravity field solutions using a Kalman smoother JF - Journal of geodynamics N2 - Different GRACE data analysis centers provide temporal variations of the Earth's gravity field as monthly, 10-daily or weekly solutions. These temporal mean fields cannot model the variations occurring during the respective time span. The aim of our approach is to extract as much temporal information as possible out of the given GRACE data. Therefore the temporal resolution shall be increased with the goal to derive daily snapshots. Yet, such an increase in temporal resolution is accompanied by a loss of redundancy and therefore in a reduced accuracy if the daily solutions are calculated individually. The approach presented here therefore introduces spatial and temporal correlations of the expected gravity field signal derived from geophysical models in addition to the daily observations, thus effectively constraining the spatial and temporal evolution of the GRACE solution. The GRACE data processing is then performed within the framework of a Kalman filter and smoother estimation procedure. The approach is at first investigated in a closed-loop simulation scenario and then applied to the original GRACE observations (level-1B data) to calculate daily solutions as part of the gravity field model ITG-Grace2010. Finally, the daily models are compared to vertical GPS station displacements and ocean bottom pressure observations. From these comparisons it can be concluded that particular in higher latitudes the daily solutions contain high-frequent temporal gravity field information and represent an improvement to existing geophysical models. KW - GRACE KW - Daily gravity field KW - Kalman smoother KW - ITG-Grace2010 Y1 - 2012 U6 - https://doi.org/10.1016/j.jog.2012.02.006 SN - 0264-3707 VL - 59 IS - 3 SP - 39 EP - 48 PB - Elsevier CY - Oxford ER - TY - JOUR A1 - Le Mouel, Jean-Louis A1 - Narteau, Clement A1 - Greff-Lefftz, Marianne A1 - Holschneider, Matthias T1 - Dissipation at the core-mantle boundary on a small-scale topography N2 - The parameters of the nutations are now known with a good accuracy, and the theory accounts for most of their values. Dissipative friction at the core-mantle boundary (CMB) and at the inner core boundary is an important ingredient of the theory. Up to now, viscous coupling at a smooth interface and electromagnetic coupling have been considered. In some cases they appear hardly strong enough to account for the observations. We advocate here that the CMB has a small- scale roughness and estimate the dissipation resulting from the interaction of the fluid core motion with this topography. We conclude that it might be significant Y1 - 2006 UR - http://www.agu.org/journals/jb/ U6 - https://doi.org/10.1029/2005jb003846 SN - 0148-0227 ER - TY - JOUR A1 - Lesur, Vincent A1 - Wardinski, Ingo A1 - Bärenzung, Julien A1 - Holschneider, Matthias T1 - On the frequency spectra of the core magnetic field Gauss coefficients JF - Physics of the earth and planetary interiors N2 - From monthly mean observatory data spanning 1957-2014, geomagnetic field secular variation values were calculated by annual differences. Estimates of the spherical harmonic Gauss coefficients of the core field secular variation were then derived by applying a correlation based modelling. Finally, a Fourier transform was applied to the time series of the Gauss coefficients. This process led to reliable temporal spectra of the Gauss coefficients up to spherical harmonic degree 5 or 6, and down to periods as short as 1 or 2 years depending on the coefficient. We observed that a k(-2) slope, where k is the frequency, is an acceptable approximation for these spectra, with possibly an exception for the dipole field. The monthly estimates of the core field secular variation at the observatory sites also show that large and rapid variations of the latter happen. This is an indication that geomagnetic jerks are frequent phenomena and that significant secular variation signals at short time scales - i.e. less than 2 years, could still be extracted from data to reveal an unexplored part of the core dynamics. KW - Geomagnetism KW - Core field KW - Secular variation rate of change KW - Geomagnetic jerks KW - Correlation based modelling Y1 - 2017 U6 - https://doi.org/10.1016/j.pepi.2017.05.017 SN - 0031-9201 SN - 1872-7395 VL - 276 SP - 145 EP - 158 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Makarava, N. A1 - Holschneider, Matthias T1 - Estimation of the Hurst exponent from noisy data: a Bayesian approach JF - EUROPEAN PHYSICAL JOURNAL B N2 - We consider a model based on the fractional Brownian motion under the influence of noise. We implement the Bayesian approach to estimate the Hurst exponent of the model. The robustness of the method to the noise intensity is tested using artificial data from fractional Brownian motion. We show that estimation of the parameters achieved when noise is considered explicitly in the model. Moreover, we identify the corresponding noise-amplitude level that allow to receive the correct estimation of the Hurst exponents in various cases. Y1 - 2012 U6 - https://doi.org/10.1140/epjb/e2012-30221-1 SN - 1434-6028 VL - 85 IS - 8 PB - SPRINGER CY - NEW YORK ER - TY - JOUR A1 - Makarava, Natallia A1 - Benmehdi, Sabah A1 - Holschneider, Matthias T1 - Bayesian estimation of self-similarity exponent JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - In this study we propose a Bayesian approach to the estimation of the Hurst exponent in terms of linear mixed models. Even for unevenly sampled signals and signals with gaps, our method is applicable. We test our method by using artificial fractional Brownian motion of different length and compare it with the detrended fluctuation analysis technique. The estimation of the Hurst exponent of a Rosenblatt process is shown as an example of an H-self-similar process with non-Gaussian dimensional distribution. Additionally, we perform an analysis with real data, the Dow-Jones Industrial Average closing values, and analyze its temporal variation of the Hurst exponent. Y1 - 2011 U6 - https://doi.org/10.1103/PhysRevE.84.021109 SN - 1539-3755 VL - 84 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Makarava, Natallia A1 - Bettenbühl, Mario A1 - Engbert, Ralf A1 - Holschneider, Matthias T1 - Bayesian estimation of the scaling parameter of fixational eye movements JF - epl : a letters journal exploring the frontiers of physics N2 - In this study we re-evaluate the estimation of the self-similarity exponent of fixational eye movements using Bayesian theory. Our analysis is based on a subsampling decomposition, which permits an analysis of the signal up to some scale factor. We demonstrate that our approach can be applied to simulated data from mathematical models of fixational eye movements to distinguish the models' properties reliably. Y1 - 2012 U6 - https://doi.org/10.1209/0295-5075/100/40003 SN - 0295-5075 VL - 100 IS - 4 PB - EDP Sciences CY - Mulhouse ER - TY - JOUR A1 - Makarava, Natallia A1 - Menz, Stephan A1 - Theves, Matthias A1 - Huisinga, Wilhelm A1 - Beta, Carsten A1 - Holschneider, Matthias T1 - Quantifying the degree of persistence in random amoeboid motion based on the Hurst exponent of fractional Brownian motion JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - Amoebae explore their environment in a random way, unless external cues like, e. g., nutrients, bias their motion. Even in the absence of cues, however, experimental cell tracks show some degree of persistence. In this paper, we analyzed individual cell tracks in the framework of a linear mixed effects model, where each track is modeled by a fractional Brownian motion, i.e., a Gaussian process exhibiting a long-term correlation structure superposed on a linear trend. The degree of persistence was quantified by the Hurst exponent of fractional Brownian motion. Our analysis of experimental cell tracks of the amoeba Dictyostelium discoideum showed a persistent movement for the majority of tracks. Employing a sliding window approach, we estimated the variations of the Hurst exponent over time, which allowed us to identify points in time, where the correlation structure was distorted ("outliers"). Coarse graining of track data via down-sampling allowed us to identify the dependence of persistence on the spatial scale. While one would expect the (mode of the) Hurst exponent to be constant on different temporal scales due to the self-similarity property of fractional Brownian motion, we observed a trend towards stronger persistence for the down-sampled cell tracks indicating stronger persistence on larger time scales. Y1 - 2014 U6 - https://doi.org/10.1103/PhysRevE.90.042703 SN - 1539-3755 SN - 1550-2376 VL - 90 IS - 4 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Matuschek, Hannes A1 - Kliegl, Reinhold A1 - Holschneider, Matthias T1 - Smoothing Spline ANOVA Decomposition of Arbitrary Splines: An Application to Eye Movements in Reading JF - PLoS one N2 - The Smoothing Spline ANOVA (SS-ANOVA) requires a specialized construction of basis and penalty terms in order to incorporate prior knowledge about the data to be fitted. Typically, one resorts to the most general approach using tensor product splines. This implies severe constraints on the correlation structure, i.e. the assumption of isotropy of smoothness can not be incorporated in general. This may increase the variance of the spline fit, especially if only a relatively small set of observations are given. In this article, we propose an alternative method that allows to incorporate prior knowledge without the need to construct specialized bases and penalties, allowing the researcher to choose the spline basis and penalty according to the prior knowledge of the observations rather than choosing them according to the analysis to be done. The two approaches are compared with an artificial example and with analyses of fixation durations during reading. Y1 - 2015 U6 - https://doi.org/10.1371/journal.pone.0119165 SN - 1932-6203 VL - 10 IS - 3 PB - PLoS CY - San Fransisco ER - TY - JOUR A1 - Mauerberger, Stefan A1 - Schanner, Maximilian Arthus A1 - Korte, Monika A1 - Holschneider, Matthias T1 - Correlation based snapshot models of the archeomagnetic field JF - Geophysical journal international N2 - For the time stationary global geomagnetic field, a new modelling concept is presented. A Bayesian non-parametric approach provides realistic location dependent uncertainty estimates. Modelling related variabilities are dealt with systematically by making little subjective apriori assumptions. Rather than parametrizing the model by Gauss coefficients, a functional analytic approach is applied. The geomagnetic potential is assumed a Gaussian process to describe a distribution over functions. Apriori correlations are given by an explicit kernel function with non-informative dipole contribution. A refined modelling strategy is proposed that accommodates non-linearities of archeomagnetic observables: First, a rough field estimate is obtained considering only sites that provide full field vector records. Subsequently, this estimate supports the linearization that incorporates the remaining incomplete records. The comparison of results for the archeomagnetic field over the past 1000 yr is in general agreement with previous models while improved model uncertainty estimates are provided. KW - geopotential theory KW - archaeomagnetism KW - magnetic field variations through KW - time KW - palaeomagnetism KW - inverse theory KW - statistical methods Y1 - 2020 U6 - https://doi.org/10.1093/gji/ggaa336 SN - 0956-540X SN - 1365-246X VL - 223 IS - 1 SP - 648 EP - 665 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Moldenhawer, Ted A1 - Moreno, Eduardo A1 - Schindler, Daniel A1 - Flemming, Sven A1 - Holschneider, Matthias A1 - Huisinga, Wilhelm A1 - Alonso, Sergio A1 - Beta, Carsten T1 - Spontaneous transitions between amoeboid and keratocyte-like modes of migration JF - Frontiers in Cell and Developmental Biology N2 - The motility of adherent eukaryotic cells is driven by the dynamics of the actin cytoskeleton. Despite the common force-generating actin machinery, different cell types often show diverse modes of locomotion that differ in their shape dynamics, speed, and persistence of motion. Recently, experiments in Dictyostelium discoideum have revealed that different motility modes can be induced in this model organism, depending on genetic modifications, developmental conditions, and synthetic changes of intracellular signaling. Here, we report experimental evidence that in a mutated D. discoideum cell line with increased Ras activity, switches between two distinct migratory modes, the amoeboid and fan-shaped type of locomotion, can even spontaneously occur within the same cell. We observed and characterized repeated and reversible switchings between the two modes of locomotion, suggesting that they are distinct behavioral traits that coexist within the same cell. We adapted an established phenomenological motility model that combines a reaction-diffusion system for the intracellular dynamics with a dynamic phase field to account for our experimental findings. KW - cell migration KW - amoeboid motility KW - keratocytle-like motility KW - modes of KW - migration KW - D. discoideum KW - actin dynamics Y1 - 2022 U6 - https://doi.org/10.3389/fcell.2022.898351 SN - 2296-634X VL - 10 PB - Frontiers Media CY - Lausanne ER - TY - JOUR A1 - Molkenthin, Christian A1 - Donner, Christian A1 - Reich, Sebastian A1 - Zöller, Gert A1 - Hainzl, Sebastian A1 - Holschneider, Matthias A1 - Opper, Manfred T1 - GP-ETAS: semiparametric Bayesian inference for the spatio-temporal epidemic type aftershock sequence model JF - Statistics and Computing N2 - The spatio-temporal epidemic type aftershock sequence (ETAS) model is widely used to describe the self-exciting nature of earthquake occurrences. While traditional inference methods provide only point estimates of the model parameters, we aim at a fully Bayesian treatment of model inference, allowing naturally to incorporate prior knowledge and uncertainty quantification of the resulting estimates. Therefore, we introduce a highly flexible, non-parametric representation for the spatially varying ETAS background intensity through a Gaussian process (GP) prior. Combined with classical triggering functions this results in a new model formulation, namely the GP-ETAS model. We enable tractable and efficient Gibbs sampling by deriving an augmented form of the GP-ETAS inference problem. This novel sampling approach allows us to assess the posterior model variables conditioned on observed earthquake catalogues, i.e., the spatial background intensity and the parameters of the triggering function. Empirical results on two synthetic data sets indicate that GP-ETAS outperforms standard models and thus demonstrate the predictive power for observed earthquake catalogues including uncertainty quantification for the estimated parameters. Finally, a case study for the l'Aquila region, Italy, with the devastating event on 6 April 2009, is presented. KW - Self-exciting point process KW - Hawkes process KW - Spatio-temporal ETAS model KW - Bayesian inference KW - Sampling KW - Earthquake modeling KW - Gaussian process KW - Data augmentation Y1 - 2022 U6 - https://doi.org/10.1007/s11222-022-10085-3 SN - 0960-3174 SN - 1573-1375 VL - 32 IS - 2 PB - Springer CY - Dordrecht 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 - 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 - 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 - Prokhorov, Boris E. A1 - Foerster, M. A1 - He, M. A1 - Namgaladze, Alexander A. A1 - Holschneider, Matthias T1 - Using MFACE as input in the UAM to specify the MIT dynamics JF - Journal of geophysical research : Space physics N2 - The magnetosphere-ionosphere-thermosphere (MIT) dynamic system significantly depends on the highly variable solar wind conditions, in particular, on changes of the strength and orientation of the interplanetary magnetic field (IMF). The solar wind and IMF interactions with the magnetosphere drive the MIT system via the magnetospheric field-aligned currents (FACs). The global modeling helps us to understand the physical background of this complex system. With the present study, we test the recently developed high-resolution empirical model of field-aligned currents MFACE (a high-resolution Model of Field-Aligned Currents through Empirical orthogonal functions analysis). These FAC distributions were used as input of the time-dependent, fully self-consistent global Upper Atmosphere Model (UAM) for different seasons and various solar wind and IMF conditions. The modeling results for neutral mass density and thermospheric wind are directly compared with the CHAMP satellite measurements. In addition, we perform comparisons with the global empirical models: the thermospheric wind model (HWM07) and the atmosphere density model (Naval Research Laboratory Mass Spectrometer and Incoherent Scatter Extended 2000). The theoretical model shows a good agreement with the satellite observations and an improved behavior compared with the empirical models at high latitudes. Using the MFACE model as input parameter of the UAM model, we obtain a realistic distribution of the upper atmosphere parameters for the Northern and Southern Hemispheres during stable IMF orientation as well as during dynamic situations. This variant of the UAM can therefore be used for modeling the MIT system and space weather predictions. Y1 - 2014 U6 - https://doi.org/10.1002/2014JA019981 SN - 2169-9380 SN - 2169-9402 VL - 119 IS - 8 PB - American Geophysical Union CY - Washington 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 - Rocha, Marcia R. A1 - Vasseur, David A. A1 - Hayn, Michael A1 - Holschneider, Matthias A1 - Gaedke, Ursula T1 - Variability patterns differ between standing stock and process rates JF - Oikos N2 - Standing stocks are typically easier to measure than process rates such as production. Hence, stocks are often used as indicators of ecosystem functions although the latter are generally more strongly related to rates than to stocks. The regulation of stocks and rates and thus their variability over time may differ, as stocks constitute the net result of production and losses. Based on long-term high frequency measurements in a large, deep lake we explore the variability patterns in primary and bacterial production and relate them to those of the corresponding standing stocks, i.e. chlorophyll concentration, phytoplankton and bacterial biomass. We employ different methods (coefficient of variation, spline fitting and spectral analysis) which complement each other for assessing the variability present in the plankton data, at different temporal scales. In phytoplankton, we found that the overall variability of primary production is dominated by fluctuations at low frequencies, such as the annual, whereas in stocks and chlorophyll in particular, higher frequencies contribute substantially to the overall variance. This suggests that using standing stocks instead of rate measures leads to an under- or overestimation of food shortage for consumers during distinct periods of the year. The range of annual variation in bacterial production is 8 times greater than biomass, showing that the variability of bacterial activity (e.g. oxygen consumption, remineralisation) would be underestimated if biomass is used. The P/B ratios were variable and although clear trends are present in both bacteria and phytoplankton, no systematic relationship between stock and rate measures were found for the two groups. Hence, standing stock and process rate measures exhibit different variability patterns and care is needed when interpreting the mechanisms and implications of the variability encountered. Y1 - 2011 U6 - https://doi.org/10.1111/j.1600-0706.2010.18786.x SN - 0030-1299 VL - 120 IS - 1 SP - 17 EP - 25 PB - Wiley-Blackwell CY - Malden ER - TY - JOUR A1 - Ropp, Guillaume A1 - Lesur, Vincent A1 - Bärenzung, Julien A1 - Holschneider, Matthias T1 - Sequential modelling of the Earth’s core magnetic field JF - Earth, Planets and Space N2 - We describe a new, original approach to the modelling of the Earth's magnetic field. The overall objective of this study is to reliably render fast variations of the core field and its secular variation. This method combines a sequential modelling approach, a Kalman filter, and a correlation-based modelling step. Sources that most significantly contribute to the field measured at the surface of the Earth are modelled. Their separation is based on strong prior information on their spatial and temporal behaviours. We obtain a time series of model distributions which display behaviours similar to those of recent models based on more classic approaches, particularly at large temporal and spatial scales. Interesting new features and periodicities are visible in our models at smaller time and spatial scales. An important aspect of our method is to yield reliable error bars for all model parameters. These errors, however, are only as reliable as the description of the different sources and the prior information used are realistic. Finally, we used a slightly different version of our method to produce candidate models for the thirteenth edition of the International Geomagnetic Reference Field. KW - geomagnetic field KW - secular variation KW - Kalman filter KW - IGRF Y1 - 2020 U6 - https://doi.org/10.1186/s40623-020-01230-1 SN - 1880-5981 VL - 72 IS - 1 PB - Springer CY - New York 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 - JOUR A1 - Sanchez, S. A1 - Wicht, J. A1 - Baerenzung, Julien A1 - Holschneider, Matthias T1 - Sequential assimilation of geomagnetic observations BT - perspectives for the reconstruction and prediction of core dynamics JF - Geophysical journal international N2 - High-precision observations of the present-day geomagnetic field by ground-based observatories and satellites provide unprecedented conditions for unveiling the dynamics of the Earth’s core. Combining geomagnetic observations with dynamo simulations in a data assimilation (DA) framework allows the reconstruction of past and present states of the internal core dynamics. The essential information that couples the internal state to the observations is provided by the statistical correlations from a numerical dynamo model in the form of a model covariance matrix. Here we test a sequential DA framework, working through a succession of forecast and analysis steps, that extracts the correlations from an ensemble of dynamo models. The primary correlations couple variables of the same azimuthal wave number, reflecting the predominant axial symmetry of the magnetic field. Synthetic tests show that the scheme becomes unstable when confronted with high-precision geomagnetic observations. Our study has identified spurious secondary correlations as the origin of the problem. Keeping only the primary correlations by localizing the covariance matrix with respect to the azimuthal wave number suffices to stabilize the assimilation. While the first analysis step is fundamental in constraining the large-scale interior state, further assimilation steps refine the smaller and more dynamical scales. This refinement turns out to be critical for long-term geomagnetic predictions. Increasing the assimilation steps from one to 18 roughly doubles the prediction horizon for the dipole from about  tree to six centuries, and from 30 to about  60 yr for smaller observable scales. This improvement is also reflected on the predictability of surface intensity features such as the South Atlantic Anomaly. Intensity prediction errors are decreased roughly by a half when assimilating long observation sequences. KW - Magnetic field variations through time KW - Core dynamics KW - Dynamo: theories and simulations KW - Inverse theory KW - Probabilistic forecasting Y1 - 2019 U6 - https://doi.org/10.1093/gji/ggz090 SN - 0956-540X SN - 1365-246X VL - 217 IS - 2 SP - 1434 EP - 1450 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Schachtschneider, R. A1 - Holschneider, Matthias A1 - Mandea, M. T1 - Error distribution in regional modelling of the geomagnetic field JF - Geophysical journal international N2 - In this study we analyse the error distribution in regional models of the geomagnetic field. Our main focus is to investigate the distribution of errors when combining two regional patches to obtain a global field from regional ones. To simulate errors in overlapping patches we choose two different data region shapes that resemble that scenario. First, we investigate the errors in elliptical regions and secondly we choose a region obtained from two overlapping circular spherical caps. We conduct a Monte-Carlo simulation using synthetic data to obtain the expected mean errors. For the elliptical regions the results are similar to the ones obtained for circular spherical caps: the maximum error at the boundary decreases towards the centre of the region. A new result emerges as errors at the boundary vary with azimuth, being largest in the major axis direction and minimal in the minor axis direction. Inside the region there is an error decay towards a minimum at the centre at a rate similar to the one in circular regions. In the case of two combined circular regions there is also an error decay from the boundary towards the centre. The minimum error occurs at the centre of the combined regions. The maximum error at the boundary occurs on the line containing the two cap centres, the minimum in the perpendicular direction where the two circular cap boundaries meet. The large errors at the boundary are eliminated by combining regional patches. We propose an algorithm for finding the boundary region that is applicable to irregularly shaped model regions. KW - Geopotential theory KW - Magnetic anomalies: modelling and interpretation KW - Satellite magnetics Y1 - 2012 U6 - https://doi.org/10.1111/j.1365-246X.2012.05675.x SN - 0956-540X VL - 191 IS - 3 SP - 1015 EP - 1024 PB - Wiley-Blackwell CY - Hoboken ER - TY - JOUR A1 - Schanner, Maximilian Arthus A1 - Mauerberger, Stefan A1 - Korte, Monika A1 - Holschneider, Matthias T1 - Correlation based time evolution of the archeomagnetic field JF - Journal of geophysical research : JGR ; an international quarterly. B, Solid earth N2 - In a previous study, a new snapshot modeling concept for the archeomagnetic field was introduced (Mauerberger et al., 2020, ). By assuming a Gaussian process for the geomagnetic potential, a correlation-based algorithm was presented, which incorporates a closed-form spatial correlation function. This work extends the suggested modeling strategy to the temporal domain. A space-time correlation kernel is constructed from the tensor product of the closed-form spatial correlation kernel with a squared exponential kernel in time. Dating uncertainties are incorporated into the modeling concept using a noisy input Gaussian process. All but one modeling hyperparameters are marginalized, to reduce their influence on the outcome and to translate their variability to the posterior variance. The resulting distribution incorporates uncertainties related to dating, measurement and modeling process. Results from application to archeomagnetic data show less variation in the dipole than comparable models, but are in general agreement with previous findings. Y1 - 2021 U6 - https://doi.org/10.1029/2020JB021548 SN - 2169-9313 SN - 2169-9356 VL - 126 IS - 7 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Schanner, Maximilian A1 - Korte, Monika A1 - Holschneider, Matthias T1 - ArchKalmag14k: A kalman-filter based global geomagnetic model for the holocene JF - Journal of geophysical research : Solid earth N2 - We propose a global geomagnetic field model for the last 14 thousand years, based on thermoremanent records. We call the model ArchKalmag14k. ArchKalmag14k is constructed by modifying recently proposed algorithms, based on space-time correlations. Due to the amount of data and complexity of the model, the full Bayesian posterior is numerically intractable. To tackle this, we sequentialize the inversion by implementing a Kalman-filter with a fixed time step. Every step consists of a prediction, based on a degree dependent temporal covariance, and a correction via Gaussian process regression. Dating errors are treated via a noisy input formulation. Cross correlations are reintroduced by a smoothing algorithm and model parameters are inferred from the data. Due to the specific statistical nature of the proposed algorithms, the model comes with space and time-dependent uncertainty estimates. The new model ArchKalmag14k shows less variation in the large-scale degrees than comparable models. Local predictions represent the underlying data and agree with comparable models, if the location is sampled well. Uncertainties are bigger for earlier times and in regions of sparse data coverage. We also use ArchKalmag14k to analyze the appearance and evolution of the South Atlantic anomaly together with reverse flux patches at the core-mantle boundary, considering the model uncertainties. While we find good agreement with earlier models for recent times, our model suggests a different evolution of intensity minima prior to 1650 CE. In general, our results suggest that prior to 6000 BCE the data is not sufficient to support global models. Y1 - 2022 U6 - https://doi.org/10.1029/2021JB023166 SN - 2169-9313 SN - 2169-9356 VL - 127 IS - 2 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Schindler, Daniel A1 - Moldenhawer, Ted A1 - Stange, Maike A1 - Lepro, Valentino A1 - Beta, Carsten A1 - Holschneider, Matthias A1 - Huisinga, Wilhelm T1 - Analysis of protrusion dynamics in amoeboid cell motility by means of regularized contour flows JF - PLoS Computational Biology : a new community journal N2 - Amoeboid cell motility is essential for a wide range of biological processes including wound healing, embryonic morphogenesis, and cancer metastasis. It relies on complex dynamical patterns of cell shape changes that pose long-standing challenges to mathematical modeling and raise a need for automated and reproducible approaches to extract quantitative morphological features from image sequences. Here, we introduce a theoretical framework and a computational method for obtaining smooth representations of the spatiotemporal contour dynamics from stacks of segmented microscopy images. Based on a Gaussian process regression we propose a one-parameter family of regularized contour flows that allows us to continuously track reference points (virtual markers) between successive cell contours. We use this approach to define a coordinate system on the moving cell boundary and to represent different local geometric quantities in this frame of reference. In particular, we introduce the local marker dispersion as a measure to identify localized membrane expansions and provide a fully automated way to extract the properties of such expansions, including their area and growth time. The methods are available as an open-source software package called AmoePy, a Python-based toolbox for analyzing amoeboid cell motility (based on time-lapse microscopy data), including a graphical user interface and detailed documentation. Due to the mathematical rigor of our framework, we envision it to be of use for the development of novel cell motility models. We mainly use experimental data of the social amoeba Dictyostelium discoideum to illustrate and validate our approach.
Author summary Amoeboid motion is a crawling-like cell migration that plays an important key role in multiple biological processes such as wound healing and cancer metastasis. This type of cell motility results from expanding and simultaneously contracting parts of the cell membrane. From fluorescence images, we obtain a sequence of points, representing the cell membrane, for each time step. By using regression analysis on these sequences, we derive smooth representations, so-called contours, of the membrane. Since the number of measurements is discrete and often limited, the question is raised of how to link consecutive contours with each other. In this work, we present a novel mathematical framework in which these links are described by regularized flows allowing a certain degree of concentration or stretching of neighboring reference points on the same contour. This stretching rate, the so-called local dispersion, is used to identify expansions and contractions of the cell membrane providing a fully automated way of extracting properties of these cell shape changes. We applied our methods to time-lapse microscopy data of the social amoeba Dictyostelium discoideum. Y1 - 2021 U6 - https://doi.org/10.1371/journal.pcbi.1009268 SN - 1553-734X SN - 1553-7358 VL - 17 IS - 8 PB - PLoS CY - San Fransisco ER - TY - JOUR A1 - Schroeter, M-A A1 - Ritter, M. A1 - Holschneider, Matthias A1 - Sturm, H. T1 - Enhanced DySEM imaging of cantilever motion using artificial structures patterned by focused ion beam techniques JF - Journal of micromechanics and microengineering N2 - We use a dynamic scanning electron microscope (DySEM) to map the spatial distribution of the vibration of a cantilever beam. The DySEM measurements are based on variations of the local secondary electron signal within the imaging electron beam diameter during an oscillation period of the cantilever. For this reason, the surface of a cantilever without topography or material variation does not allow any conclusions about the spatial distribution of vibration due to a lack of dynamic contrast. In order to overcome this limitation, artificial structures were added at defined positions on the cantilever surface using focused ion beam lithography patterning. The DySEM signal of such high-contrast structures is strongly improved, hence information about the surface vibration becomes accessible. Simulations of images of the vibrating cantilever have also been performed. The results of the simulation are in good agreement with the experimental images. KW - FIB patterning KW - structured cantilever KW - AFM KW - modal analysis KW - DySEM Y1 - 2016 U6 - https://doi.org/10.1088/0960-1317/26/3/035010 SN - 0960-1317 SN - 1361-6439 VL - 26 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Schröter, M-A A1 - Holschneider, Matthias A1 - Sturm, H. T1 - Analytical and numerical analysis of imaging mechanism of dynamic scanning electron microscopy JF - Nanotechnology N2 - The direct observation of small oscillating structures with the help of a scanning electron beam is a new approach to study the vibrational dynamics of cantilevers and microelectromechanical systems. In the scanning electron microscope, the conventional signal of secondary electrons (SE, dc part) is separated from the signal response of the SE detector, which is correlated to the respective excitation frequency for vibration by means of a lock-in amplifier. The dynamic response is separated either into images of amplitude and phase shift or into real and imaginary parts. Spatial resolution is limited to the diameter of the electron beam. The sensitivity limit to vibrational motion is estimated to be sub-nanometer for high integration times. Due to complex imaging mechanisms, a theoretical model was developed for the interpretation of the obtained measurements, relating cantilever shapes to interaction processes consisting of incident electron beam, electron-lever interaction, emitted electrons and detector response. Conclusions drawn from this new model are compared with numerical results based on the Euler-Bernoulli equation. Y1 - 2012 U6 - https://doi.org/10.1088/0957-4484/23/43/435501 SN - 0957-4484 VL - 23 IS - 43 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Schröter, M-A A1 - Sturm, H. A1 - Holschneider, Matthias T1 - Phase and amplitude patterns in DySEM mappings of vibrating microstructures JF - Nanotechnology N2 - We use a dynamic scanning electron microscope (DySEM) to analyze the movement of oscillating micromechanical structures. A dynamic secondary electron (SE) signal is recorded and correlated to the oscillatory excitation of scanning force microscope (SFM) cantilever by means of lock-in amplifiers. We show, how the relative phase of the oscillations modulate the resulting real part and phase pictures of the DySEM mapping. This can be used to obtain information about the underlying oscillatory dynamics. We apply the theory to the case of a cantilever in oscillation, driven at different flexural and torsional resonance modes. This is an extension of a recent work (Schroter et al 2012 Nanotechnology 23 435501), where we reported on a general methodology to distinguish nonlinear features caused by the imaging process from those caused by cantilever motion. Y1 - 2013 U6 - https://doi.org/10.1088/0957-4484/24/21/215701 SN - 0957-4484 VL - 24 IS - 21 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Schütz, Nadine A1 - Holschneider, Matthias T1 - Detection of trend changes in time series using Bayesian inference JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - Change points in time series are perceived as isolated singularities where two regular trends of a given signal do not match. The detection of such transitions is of fundamental interest for the understanding of the system's internal dynamics or external forcings. In practice observational noise makes it difficult to detect such change points in time series. In this work we elaborate on a Bayesian algorithm to estimate the location of the singularities and to quantify their credibility. We validate the performance and sensitivity of our inference method by estimating change points of synthetic data sets. As an application we use our algorithm to analyze the annual flow volume of the Nile River at Aswan from 1871 to 1970, where we confirm a well-established significant transition point within the time series. Y1 - 2011 U6 - https://doi.org/10.1103/PhysRevE.84.021120 SN - 1539-3755 SN - 1550-2376 VL - 84 IS - 2 PB - American Physical Society CY - College Park 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 - Shebalin, Peter N. A1 - Narteau, Clement A1 - Zechar, Jeremy Douglas A1 - Holschneider, Matthias T1 - Combining earthquake forecasts using differential probability gains JF - Earth, planets and space N2 - We describe an iterative method to combine seismicity forecasts. With this method, we produce the next generation of a starting forecast by incorporating predictive skill from one or more input forecasts. For a single iteration, we use the differential probability gain of an input forecast relative to the starting forecast. At each point in space and time, the rate in the next-generation forecast is the product of the starting rate and the local differential probability gain. The main advantage of this method is that it can produce high forecast rates using all types of numerical forecast models, even those that are not rate-based. Naturally, a limitation of this method is that the input forecast must have some information not already contained in the starting forecast. We illustrate this method using the Every Earthquake a Precursor According to Scale (EEPAS) and Early Aftershocks Statistics (EAST) models, which are currently being evaluated at the US testing center of the Collaboratory for the Study of Earthquake Predictability. During a testing period from July 2009 to December 2011 (with 19 target earthquakes), the combined model we produce has better predictive performance - in terms of Molchan diagrams and likelihood - than the starting model (EEPAS) and the input model (EAST). Many of the target earthquakes occur in regions where the combined model has high forecast rates. Most importantly, the rates in these regions are substantially higher than if we had simply averaged the models. KW - Probabilistic forecasting KW - Earthquake interaction KW - Forecasting and prediction KW - Statistical seismology Y1 - 2014 U6 - https://doi.org/10.1186/1880-5981-66-37 SN - 1880-5981 VL - 66 PB - Springer CY - Heidelberg ER - TY - JOUR A1 - Shebalin, Peter A1 - Narteau, Clement A1 - Holschneider, Matthias T1 - From alarm-based to rate-based earthquake forecast models JF - Bulletin of the Seismological Society of America N2 - We propose a conversion method from alarm-based to rate-based earthquake forecast models. A differential probability gain g(alarm)(ref) is the absolute value of the local slope of the Molchan trajectory that evaluates the performance of the alarm-based model with respect to the chosen reference model. We consider that this differential probability gain is constant over time. Its value at each point of the testing region depends only on the alarm function value. The rate-based model is the product of the event rate of the reference model at this point multiplied by the corresponding differential probability gain. Thus, we increase or decrease the initial rates of the reference model according to the additional amount of information contained in the alarm-based model. Here, we apply this method to the Early Aftershock STatistics (EAST) model, an alarm-based model in which early aftershocks are used to identify space-time regions with a higher level of stress and, consequently, a higher seismogenic potential. The resulting rate-based model shows similar performance to the original alarm-based model for all ranges of earthquake magnitude in both retrospective and prospective tests. This conversion method offers the opportunity to perform all the standard evaluation tests of the earthquake testing centers on alarm-based models. In addition, we infer that it can also be used to consecutively combine independent forecast models and, with small modifications, seismic hazard maps with short- and medium-term forecasts. Y1 - 2012 U6 - https://doi.org/10.1785/0120110126 SN - 0037-1106 VL - 102 IS - 1 SP - 64 EP - 72 PB - Seismological Society of America CY - Albany 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 - Shin, Seoleun A1 - Zöller, Gert A1 - Holschneider, Matthias A1 - Reich, Sebastian T1 - A multigrid solver for modeling complex interseismic stress fields JF - Computers & geosciences : an international journal devoted to the publication of papers on all aspects of geocomputation and to the distribution of computer programs and test data sets ; an official journal of the International Association for Mathematical Geology N2 - We develop a multigrid, multiple time stepping scheme to reduce computational efforts for calculating complex stress interactions in a strike-slip 2D planar fault for the simulation of seismicity. The key elements of the multilevel solver are separation of length scale, grid-coarsening, and hierarchy. In this study the complex stress interactions are split into two parts: the first with a small contribution is computed on a coarse level, and the rest for strong interactions is on a fine level. This partition leads to a significant reduction of the number of computations. The reduction of complexity is even enhanced by combining the multigrid with multiple time stepping. Computational efficiency is enhanced by a factor of 10 while retaining a reasonable accuracy, compared to the original full matrix-vortex multiplication. The accuracy of solution and computational efficiency depend on a given cut-off radius that splits multiplications into the two parts. The multigrid scheme is constructed in such a way that it conserves stress in the entire half-space. KW - Multigrid KW - Multiple time stepping KW - Strike-slip fault model Y1 - 2011 U6 - https://doi.org/10.1016/j.cageo.2010.11.011 SN - 0098-3004 VL - 37 IS - 8 SP - 1075 EP - 1082 PB - Elsevier CY - Oxford ER - TY - JOUR A1 - Wang, Lifeng A1 - Hainzl, Sebastian A1 - Zöller, Gert A1 - Holschneider, Matthias T1 - Stress- and aftershock-constrained joint inversions for coseismic and postseismic slip applied to the 2004 M6.0 Parkfield earthquake JF - Journal of geophysical research : Solid earth N2 - Both aftershocks and geodetically measured postseismic displacements are important markers of the stress relaxation process following large earthquakes. Postseismic displacements can be related to creep-like relaxation in the vicinity of the coseismic rupture by means of inversion methods. However, the results of slip inversions are typically non-unique and subject to large uncertainties. Therefore, we explore the possibility to improve inversions by mechanical constraints. In particular, we take into account the physical understanding that postseismic deformation is stress-driven, and occurs in the coseismically stressed zone. We do joint inversions for coseismic and postseismic slip in a Bayesian framework in the case of the 2004 M6.0 Parkfield earthquake. We perform a number of inversions with different constraints, and calculate their statistical significance. According to information criteria, the best result is preferably related to a physically reasonable model constrained by the stress-condition (namely postseismic creep is driven by coseismic stress) and the condition that coseismic slip and large aftershocks are disjunct. This model explains 97% of the coseismic displacements and 91% of the postseismic displacements during day 1-5 following the Parkfield event, respectively. It indicates that the major postseismic deformation can be generally explained by a stress relaxation process for the Parkfield case. This result also indicates that the data to constrain the coseismic slip model could be enriched postseismically. For the 2004 Parkfield event, we additionally observe asymmetric relaxation process at the two sides of the fault, which can be explained by material contrast ratio across the fault of similar to 1.15 in seismic velocity. Y1 - 2012 U6 - https://doi.org/10.1029/2011JB009017 SN - 2169-9313 SN - 2169-9356 VL - 117 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Zaourar, N. A1 - Hamoudi, M. A1 - Mandea, M. A1 - Balasis, G. A1 - Holschneider, Matthias T1 - Wavelet-based multiscale analysis of geomagnetic disturbance JF - EARTH PLANETS AND SPACE N2 - The dynamics of external contributions to the geomagnetic field is investigated by applying time-frequency methods to magnetic observatory data. Fractal models and multiscale analysis enable obtaining maximum quantitative information related to the short-term dynamics of the geomagnetic field activity. The stochastic properties of the horizontal component of the transient external field are determined by searching for scaling laws in the power spectra. The spectrum fits a power law with a scaling exponent beta, a typical characteristic of self-affine time-series. Local variations in the power-law exponent are investigated by applying wavelet analysis to the same time-series. These analyses highlight the self-affine properties of geomagnetic perturbations and their persistence. Moreover, they show that the main phases of sudden storm disturbances are uniquely characterized by a scaling exponent varying between 1 and 3, possibly related to the energy contained in the external field. These new findings suggest the existence of a long-range dependence, the scaling exponent being an efficient indicator of geomagnetic activity and singularity detection. These results show that by using magnetogram regularity to reflect the magnetosphere activity, a theoretical analysis of the external geomagnetic field based on local power-law exponents is possible. KW - Geomagnetic field KW - magnetosphere KW - geomagnetic storm KW - multiscale analysis KW - spectral exponent Y1 - 2013 U6 - https://doi.org/10.5047/eps.2013.05.001 SN - 1343-8832 SN - 1880-5981 VL - 65 IS - 12 SP - 1525 EP - 1540 PB - TERRA SCIENTIFIC PUBL CO CY - TOKYO ER - TY - JOUR A1 - Zaourar, Naima A1 - Hamoudi, Mohamed A1 - Holschneider, Matthias A1 - Mandea, Mioara T1 - Fractal dynamics of geomagnetic storms JF - Arabian journal of geosciences N2 - We explore fluctuations of the horizontal component of the Earth's magnetic field to identify scaling behaviour of the temporal variability in geomagnetic data recorded by the Intermagnet observatories during the solar cycle 23 (years 1996 to 2005). In this work, we use the remarkable ability of scaling wavelet exponents to highlight the singularities associated with discontinuities present in the magnetograms obtained at two magnetic observatories for six intense magnetic storms, including the sudden storm commencements of 14 July 2000, 29-31 October and 20-21 November 2003. In the active intervals that occurred during geomagnetic storms, we observe a rapid and unidirectional change in the spectral scaling exponent at the time of storm onset. The corresponding fractal features suggest that the dynamics of the whole time series is similar to that of a fractional Brownian motion. Our findings point to an evident relatively sudden change related to the emergence of persistency of the fractal power exponent fluctuations precedes an intense magnetic storm. These first results could be useful in the framework of extreme events prediction studies. KW - Geomagnetic field KW - Magnetosphere KW - Geomagnetic storm KW - Multiscale analysis KW - Spectral exponent Y1 - 2013 U6 - https://doi.org/10.1007/s12517-011-0487-0 SN - 1866-7511 VL - 6 IS - 6 SP - 1693 EP - 1702 PB - Springer CY - Heidelberg ER - TY - JOUR A1 - Zoller, Gert A1 - Hainzl, Sebastian A1 - Holschneider, Matthias A1 - Ben-Zion, Yehuda T1 - Aftershocks resulting from creeping sections in a heterogeneous fault N2 - We show that realistic aftershock sequences with space-time characteristics compatible with observations are generated by a model consisting of brittle fault segments separated by creeping zones. The dynamics of the brittle regions is governed by static/kinetic friction, 3D elastic stress transfer and small creep deformation. The creeping parts are characterized by high ongoing creep velocities. These regions store stress during earthquake failures and then release it in the interseismic periods. The resulting postseismic deformation leads to aftershock sequences following the modified Omori law. The ratio of creep coefficients in the brittle and creeping sections determines the duration of the postseismic transients and the exponent p of the modified Omori law Y1 - 2005 SN - 0094-8276 ER - TY - JOUR A1 - Zoller, Gert A1 - Holschneider, Matthias A1 - Ben-Zion, Yehuda T1 - The role of heterogeneities as a tuning parameter of earthquake dynamics N2 - We investigate the influence of spatial heterogeneities on various aspects of brittle failure and seismicity in a model of a large strike-slip fault. The model dynamics is governed by realistic boundary conditions consisting of constant velocity motion of regions around the fault, static/kinetic friction laws, creep with depth-dependent coefficients, and 3-D elastic stress transfer. The dynamic rupture is approximated on a continuous time scale using a finite stress propagation velocity ("quasidynamic model''). The model produces a "brittle- ductile'' transition at a depth of about 12.5 km, realistic hypocenter distributions, and other features of seismicity compatible with observations. Previous work suggested that the range of size scales in the distribution of strength-stress heterogeneities acts as a tuning parameter of the dynamics. Here we test this hypothesis by performing a systematic parameter-space study with different forms of heterogeneities. In particular, we analyze spatial heterogeneities that can be tuned by a single parameter in two distributions: ( 1) high stress drop barriers in near- vertical directions and ( 2) spatial heterogeneities with fractal properties and variable fractal dimension. The results indicate that the first form of heterogeneities provides an effective means of tuning the behavior while the second does not. In relatively homogeneous cases, the fault self-organizes to large-scale patches and big events are associated with inward failure of individual patches and sequential failures of different patches. The frequency-size event statistics in such cases are compatible with the characteristic earthquake distribution and large events are quasi-periodic in time. In strongly heterogeneous or near-critical cases, the rupture histories are highly discontinuous and consist of complex migration patterns of slip on the fault. In such cases, the frequency-size and temporal statistics follow approximately power-law relations Y1 - 2005 SN - 0033-4553 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 - Zöller, Gert A1 - Hainzl, Sebastian A1 - Holschneider, Matthias T1 - Recurrence of Large Earthquakes : bayesian inference from catalogs in the presence of magnitude uncertainties N2 - We present a Bayesian method that allows continuous updating the aperiodicity of the recurrence time distribution of large earthquakes based on a catalog with magnitudes above a completeness threshold. The approach uses a recently proposed renewal model for seismicity and allows the inclusion of magnitude uncertainties in a straightforward manner. Errors accounting for grouped magnitudes and random errors are studied and discussed. The results indicate that a stable and realistic value of the aperiodicity can be predicted in an early state of seismicity evolution, even though only a small number of large earthquakes has occurred to date. Furthermore, we demonstrate that magnitude uncertainties can drastically influence the results and can therefore not be neglected. We show how to correct for the bias caused by magnitude errors. For the region of Parkfield we find that the aperiodicity, or the coefficient of variation, is clearly higher than in studies which are solely based on the large earthquakes. Y1 - 2010 UR - http://www.springerlink.com/content/101201 U6 - https://doi.org/10.1007/s00024-010-0078-0 SN - 0033-4553 ER - TY - JOUR A1 - Zöller, Gert A1 - Holschneider, Matthias T1 - Induced seismicity: What is the size of the largest expected earthquake? JF - The bulletin of the Seismological Society of America N2 - The injection of fluids is a well-known origin for the triggering of earthquake sequences. The growing number of projects related to enhanced geothermal systems, fracking, and others has led to the question, which maximum earthquake magnitude can be expected as a consequence of fluid injection? This question is addressed from the perspective of statistical analysis. Using basic empirical laws of earthquake statistics, we estimate the magnitude M-T of the maximum expected earthquake in a predefined future time window T-f. A case study of the fluid injection site at Paradox Valley, Colorado, demonstrates that the magnitude m 4.3 of the largest observed earthquake on 27 May 2000 lies very well within the expectation from past seismicity without adjusting any parameters. Vice versa, for a given maximum tolerable earthquake at an injection site, we can constrain the corresponding amount of injected fluids that must not be exceeded within predefined confidence bounds. Y1 - 2014 U6 - https://doi.org/10.1785/0120140195 SN - 0037-1106 SN - 1943-3573 VL - 104 IS - 6 SP - 3153 EP - 3158 PB - Seismological Society of America CY - Albany ER - TY - JOUR A1 - Zöller, Gert A1 - Holschneider, Matthias T1 - The Maximum Possible and the Maximum Expected Earthquake Magnitude for Production-Induced Earthquakes at the Gas Field in Groningen, The Netherlands JF - Bulletin of the Seismological Society of America N2 - The Groningen gas field serves as a natural laboratory for production-induced earthquakes, because no earthquakes were observed before the beginning of gas production. Increasing gas production rates resulted in growing earthquake activity and eventually in the occurrence of the 2012M(w) 3.6 Huizinge earthquake. At least since this event, a detailed seismic hazard and risk assessment including estimation of the maximum earthquake magnitude is considered to be necessary to decide on the future gas production. In this short note, we first apply state-of-the-art methods of mathematical statistics to derive confidence intervals for the maximum possible earthquake magnitude m(max). Second, we calculate the maximum expected magnitude M-T in the time between 2016 and 2024 for three assumed gas-production scenarios. Using broadly accepted physical assumptions and 90% confidence level, we suggest a value of m(max) 4.4, whereas M-T varies between 3.9 and 4.3, depending on the production scenario. Y1 - 2016 U6 - https://doi.org/10.1785/0120160220 SN - 0037-1106 SN - 1943-3573 VL - 106 SP - 2917 EP - 2921 PB - Seismological Society of America CY - Albany ER - TY - JOUR A1 - Zöller, Gert A1 - Holschneider, Matthias T1 - The Earthquake History in a Fault Zone Tells Us Almost Nothing about m(max) JF - Seismological research letters N2 - In the present study, we summarize and evaluate the endeavors from recent years to estimate the maximum possible earthquake magnitude m(max) from observed data. In particular, we use basic and physically motivated assumptions to identify best cases and worst cases in terms of lowest and highest degree of uncertainty of m(max). In a general framework, we demonstrate that earthquake data and earthquake proxy data recorded in a fault zone provide almost no information about m(max) unless reliable and homogeneous data of a long time interval, including several earthquakes with magnitude close to m(max), are available. Even if detailed earthquake information from some centuries including historic and paleoearthquakes are given, only very few, namely the largest events, will contribute at all to the estimation of m(max), and this results in unacceptably high uncertainties. As a consequence, estimators of m(max) in a fault zone, which are based solely on earthquake-related information from this region, have to be dismissed. Y1 - 2016 U6 - https://doi.org/10.1785/0220150176 SN - 0895-0695 SN - 1938-2057 VL - 87 SP - 132 EP - 137 PB - Seismological Society of America CY - Albany ER - TY - JOUR A1 - Zöller, Gert A1 - Holschneider, Matthias A1 - Hainzl, Sebastian T1 - The Maximum Earthquake Magnitude in a Time Horizon: Theory and Case Studies JF - Bulletin of the Seismological Society of America N2 - We show how the maximum magnitude within a predefined future time horizon may be estimated from an earthquake catalog within the context of Gutenberg-Richter statistics. The aim is to carry out a rigorous uncertainty assessment, and calculate precise confidence intervals based on an imposed level of confidence a. In detail, we present a model for the estimation of the maximum magnitude to occur in a time interval T-f in the future, given a complete earthquake catalog for a time period T in the past and, if available, paleoseismic events. For this goal, we solely assume that earthquakes follow a stationary Poisson process in time with unknown productivity Lambda and obey the Gutenberg-Richter law in magnitude domain with unknown b-value. The random variables. and b are estimated by means of Bayes theorem with noninformative prior distributions. Results based on synthetic catalogs and on retrospective calculations of historic catalogs from the highly active area of Japan and the low-seismicity, but high-risk region lower Rhine embayment (LRE) in Germany indicate that the estimated magnitudes are close to the true values. Finally, we discuss whether the techniques can be extended to meet the safety requirements for critical facilities such as nuclear power plants. For this aim, the maximum magnitude for all times has to be considered. In agreement with earlier work, we find that this parameter is not a useful quantity from the viewpoint of statistical inference. Y1 - 2013 U6 - https://doi.org/10.1785/0120120013 SN - 0037-1106 VL - 103 IS - 2A SP - 860 EP - 875 PB - Seismological Society of America CY - Albany ER - TY - JOUR A1 - Zöller, Gert A1 - Holschneider, Matthias A1 - Hainzl, Sebastian A1 - Zhuang, Jiancang T1 - The largest expected earthquake magnitudes in Japan: The statistical perspective JF - Bulletin of the Seismological Society of America N2 - Earthquake catalogs are probably the most informative data source about spatiotemporal seismicity evolution. The catalog quality in one of the most active seismogenic zones in the world, Japan, is excellent, although changes in quality arising, for example, from an evolving network are clearly present. Here, we seek the best estimate for the largest expected earthquake in a given future time interval from a combination of historic and instrumental earthquake catalogs. We extend the technique introduced by Zoller et al. (2013) to estimate the maximum magnitude in a time window of length T-f for earthquake catalogs with varying level of completeness. In particular, we consider the case in which two types of catalogs are available: a historic catalog and an instrumental catalog. This leads to competing interests with respect to the estimation of the two parameters from the Gutenberg-Richter law, the b-value and the event rate lambda above a given lower-magnitude threshold (the a-value). The b-value is estimated most precisely from the frequently occurring small earthquakes; however, the tendency of small events to cluster in aftershocks, swarms, etc. violates the assumption of a Poisson process that is used for the estimation of lambda. We suggest addressing conflict by estimating b solely from instrumental seismicity and using large magnitude events from historic catalogs for the earthquake rate estimation. Applying the method to Japan, there is a probability of about 20% that the maximum expected magnitude during any future time interval of length T-f = 30 years is m >= 9.0. Studies of different subregions in Japan indicates high probabilities for M 8 earthquakes along the Tohoku arc and relatively low probabilities in the Tokai, Tonankai, and Nankai region. Finally, for scenarios related to long-time horizons and high-confidence levels, the maximum expected magnitude will be around 10. Y1 - 2014 U6 - https://doi.org/10.1785/0120130103 SN - 0037-1106 SN - 1943-3573 VL - 104 IS - 2 SP - 769 EP - 779 PB - Seismological Society of America CY - Albany ER -