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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 -